Genetic diversity of the obligate intracellular bacterium Chlamydophila pneumoniae by genome-wide analysis of single nucleotide polymorphisms: evidence for highly clonal population structure
- Thomas Rattei†1,
- Stephan Ott†2,
- Michaela Gutacker3,
- Jan Rupp4,
- Matthias Maass5,
- Stefan Schreiber2,
- Werner Solbach4,
- Thierry Wirth6 and
- Jens Gieffers4Email author
© Rattei et al; licensee BioMed Central Ltd. 2007
Received: 12 April 2007
Accepted: 04 October 2007
Published: 04 October 2007
Chlamydophila pneumoniae is an obligate intracellular bacterium that replicates in a biphasic life cycle within eukaryotic host cells. Four published genomes revealed an identity of > 99 %. This remarkable finding raised questions about the existence of distinguishable genotypes in correlation with geographical and anatomical origin.
We studied the genetic diversity of C. pneumoniae by analysing synonymous single nucleotide polymorphisms (sSNPs) that are under reduced selection pressure. We conducted an in silico analysis of the four sequenced genomes, chose 232 representative sSNPs and analysed the loci in 38 C. pneumoniae isolates. We identified 15 different genotypes that were separated in four major clusters. Clusters were not associated with anatomical or geographical origin. However, animal lineages are basal on the C. pneumomiae phylogeny, suggesting a recent transmission to humans through successive bottlenecks some 150,000 years ago. A lack of detectable variation in 17 isolates emphasizes the extraordinary genetic conservation of this species and the high clonality of the population. Moreover, the largest cluster, which encompasses 80% of all analysed strains, is an extremely young clade, that went through an important population expansion some 3,300 years ago.
sSNPs have proven useful as a sensitive marker to gain new insights into genetic diversity, population structure and evolutionary history of C. pneumoniae.
The order Chlamydiales evolved from their free-living ancestors about 500–1000 million years ago (mya) , establishing their intracellular life in lower eukaryotes. Their life cycle consists of infectious elementary bodies and intracellularly replicative reticulate bodies (for recent review see e.g. . The divergence of C. trachomatis and C. pneumoniae 50–200 mya proceeds the appearance of Homo sapiens and, thus, each lineage possessed sufficient biological potential to exploit new hosts  Gene acquisition by horizontal gene transfer (HGT) was not a driving force in the evolution of the Chlamydiales; less than 1 % of the total gene number was estimated to result from HGT . Instead, while adapting to the homeostatic niche within their specific hosts, all Chlamydiaceae species reduced their genome size to little more than 1 Mb. At least 80 % of the genes in any sequenced Chlamydiaceae genome represent orthologs with other Chlamydiaceae species [4–6]. Despite this high level of functional conservation, chlamydiae differ in their host spectrum, tissue tropism and spectrum of host diseases resulting from the infection.
Chlamydophila pneumoniae deserves medical attention as an important respiratory pathogen causing about 10 % of community acquired pneumonias and upper respiratory tract infections like bronchitis, pharyngitis and sinusitis [7, 8]. Virtually every person is believed to be infected at least once during their lifetime. Serological evidence indicates that infections start to occur during childhood and 50 – 60 % of the population has been exposed by 20 years of age . Viable chlamydiae have been isolated from atherosclerotic plaques, implicating a causal role in the development of atherosclerosis . Blood monocytes are believed to be the vector system within the systemic circulation . Moreover, C. pneumoniae has been isolated in the brain and is associated with Alzheimer's disease and Multiple Sclerosis [12, 13].
In contrast to C. trachomatis, no biovars or pathotypes have been described so far for C. pneumoniae, despite of the broad range of disease associations and sites of infection. The genomic sequences of the four isolates AR-39, CWL-029, J138 and TW-183 show a remarkable level of gene conservation and synteny. Overall identity between isolates is > 99% [4–6]. The genomes of the four isolates comprise an equivalent set of coding sequences and differ only by a small number of insertions and deletions and several single nucleotide polymorphisms (SNPs). Analyses of repeated sequences suggested that C. pneumoniae has the highest potential for recombination among fully sequenced Chlamydiaceae and that the majority of recombination hotspots of the C. pneumoniae genome are concentrated in a family of polymorphic proteins . The high degree of conservation is in line with the low rate of HGT as well as a low frequency of genome rearrangements within the species, owing to the unique chlamydial biology and ecological isolation within the intracellular niche. Therefore, point mutations, fixed in the population and accumulating with time, are the major source of genetic variability in C. pneumoniae hence SNPs are among the most sensitive phylogenetic markers to reconstruct the evolutionary history of this species. SNPs occur as two classes of substitutions within genes, referred to as synonymous and non-synonymous SNPs (sSNPs; nsSNPs). nsSNPs result in amino acid replacement and are targets of evolutionary selection. sSNPs make a difference with respect to the codon usage but do not alter the protein sequence and are therefore evolutionary almost neutral. This qualifies sSNPs as a valuable target for population genetic studies. SNPs have been successfully used to reconstruct the phylogeny of the M. tuberculosis complex  and B. anthracis isolates . In the present study, we analysed SNPs of C. pneumoniae to answer questions about the homology of geographically and anatomically different isolates and the existence of distinguishable genotypes. We first conducted an in silico analysis of sSNPs and nsSNPs of the four sequenced C. pneumoniae isolates. To reduce the potential ascertainment bias caused by deriving the SNP events from a limited set of only four completely sequenced genomes, the SNPs selected for sequencing should be as independent as possible. Compared to non-coding and non-synonymous SNPs, synonymous SNPs provide evolutionary almost neutral, abundant and equally distributed sequence markers in the C. pneumoniae genomes. We therefore chose a representative subset of 232 sSNPs and analysed the loci in 38 C. pneumoniae isolates. As chlamydiae are difficult to propagate the available collection is the largest that has even been analysed. By this approach, we were able to gain new insights into the genetic diversity, the population structure and evolutionary history of C. pneumoniae.
SNP-analysis of the reference strains
Origin of the C. pneumoniae isolates
Site of Isolation
frog (Cryptohylax gresshoffi), liver
Central African Republic
peripheral blood monocytes
peripheral blood monocytes
peripheral blood monocytes
Phylogeny of C. pneumoniae
Genotyping of C. pneumoniae
SNP data revealed to be useful to genotype new C. pneumoniae isolates and to allocate them to the described clusters. A subset of sSNPs was defined that discriminates between the genotypes. Additional file 2 gives an algorithm for the identification of different C. pneumoniae isolates. A total of 12 different SNPs (four to five SNPs per isolate) are sufficient to identify these isolates or to relate a new isolate to the characterized ones.
The evolutionary tree of the 38 C. pneumoniae isolates based on 232 sSNPs displays the genetic variability of the species. The topology ot the phylogenetic tree of the four reference isolates (figure 2) is reflected within the tree. There are four major clusters within the tree and 15 different C. pneumoniae genotypes can be defined by the sSNP analysis.
Of all reported studies, this is the highest discriminative analysis of genetic diversity in C. pneumoniae and sSNPs have proven useful to differentiate between C. pneumoniae isolates. Nevertheless, it is most intriguing that even with this sensitive marker we were not able to differentiate 17 out of 38 isolates that originated from different European and American origins and anatomical localisations. This is evidence for a highly clonal population structure and advocates for a recent origin of this bacterial species. Results are in line with the expected widespread homologous recombination in C. pneumoniae.
The genetic clusters do not separate between geographic or anatomical origin of the isolates, even if the Finnish isolates form a distinct subcluster. Isolates that infect the blood or the vasculature cannot be separated from respiratory isolates, thus we hypothesize that they are not genetically different lineages and are rather young. We calculated that the expansion event in cluster IV took place about 3.300 years ago or even more recently assuming a faster clock rate in intracellular bacteriae. Nevertheless, vascular and PBMC isolates have been shown to possess a single copy of the tyr P gene while respiratory isolates have multiple copies. A single tyr P copy only has been proposed to improve the propensity for vascular infection and persistence . We conclude that the multiplication of the tyr P gene might have occurred after the emergence of the genotypes defined by SNPs. Moreover, the identiy of vascular and respirotory isolates in the sSNP analysis does not exclude the possibility of different pathotypes, as sSNPs have no phenotypic correlate. To further study genotype-phenotype correlation, non-synonymous polymorphisms would be more suitable. Given the clonal population structure demonstrated in this study, we assume that polymorphisms defining potential pathotypes can be discovered only by the most sensitive techniques, e.g. the sequence analysis of the whole genome.
We are aware that this sSNP analysis only provides limited information on the species phylogeny and may be biased by the ascertainment of the sequenced SNPs. The SNPs we chose for analysis were identified by four reference isolates. This approach includes an inherent bias, as additional SNPs of other isolates remain unconsidered. The most likely effect of this bias could be a branch collapse and the underestimation of the evolutionary distance. However, we are convinced that the impact of this bias is limited by the fact that the reference isolates come from distant locations, and, thus, represent different genotypic lineages. TW-183 is an ocular strain obtained in Taiwan in 1965, AR-39 was isolated from a patient with acute respiratory symptoms in Seattle, USA in 1983. CWL-029 was isolated from a pneumonia patient in Atlanta, USA in 1988. J-138 was obtained from a boy with pharyngitis in Japan in 1994. This assumption is supported by the fact that each identified cluster, with the exception of the animal isolates, includes one reference isolate and no distant outgroups occurred. Additionally, when we sequenced about 8 % of the genome, we did not identify novel SNPs in human strains that were absent in the reference strains. In contrast, we detected 43 new sSNPs between animal and human isolates. Extrapolating from the 0.1 mega bases sequenced, one can speculate, that about 500 additional sSNPs could be found in the whole genome. Thus, we believe that the evolutionary distance between animal and human strains is highly underestimated in this tree (figure 4).
The major phylogenetic finding is that animal C. pneumoniae isolates from African frog and Australian koala cluster closely together and differ distinctly from human isolates. Highest homologies outside the species are found with C. abortus, a ruminant pathogen that also possess the ability to infect humans. Further analysis showed that the animal isolates are basal and the human isolates are derived. We, thus, hypothesized that C. pneumoniae evolved from an animal pathogen and the C. pneumoniae animal isolates are the ancestors of human isolates. Two previous studies analysing the 16S rRNA gene, the 16S/23S intergenic spacer and domain I of the 23S gene supports this hypothesis: A C. pneumoniae isolate from a horse clustered different from human isolates that were derived [19, 20]. The split between the lineages in our study occured in recent history, about 150,000 years ago. While spread and mixture of human isolates by migration and travel could explain their genetic similarity, it is extremely unlikely, that Australian koala and African frog were able to exchange the intracellular parasite. Thus, it is most intriguing that both isolates differ by only two SNPs despite the long period of evolutionary isolation. We interpret this as evidence for highly clonal population structure as a consequence of a very conserved genome. Additionally, we speculate that the transition into the human or animal population represents an evolutionary bottleneck that contributed to the clonal population structure. A similar evolutionary history was proposed for M tuberculosis based on the conservation of point mutations and deletions [21, 22].
When we analysed the four sequenced genomes in silico, a first remarkable finding was that the abundance of nsSNPs exceeded that of sSNPs (ratio of 1.45:1). In an environment of strong selective pressure one could expect that the number of sSNPs would exceed the number of nsSNPs because the latter are purged from the population by purifying selection. To solve this apparent contradiction we determined the likelihood for a single nucleotide exchange to cause a nsSNP or sSNP for the CWL-029 genome. Absence of any selection would result in a ratio of 3.4:1 between nsSNPs and sSNPs. Thus, the reduced number of observed nsSNP in comparison to the theoretically expected one is evidence for strong selective pressure within the intracellular niche.
We further found that all SNPs are homogeneously distributed around the genome (figure 1). In this respect, C. pneumoniae differs from C. trachomatis, which shows a significant amount of genetic variation between serovars within the "plasticity zone" around the origin of replication [4, 23]. No class of genes (e.g. outer membrane proteins) could be identified that showed more SNPs than the average (data not shown). This finding correlates with the low number of recent HGT events reported for the C. pneumoniae species .
In summary, we have shown that sSNPs are a sensitive tool to reveal the genetic diversity and past demography of C. pneumoniae that remained undetected by other methods. Analysing a good portion of the worldwide available isolates, we were able to divide C. pneumoniae into different genotypes. Nevertheless, a lack of detectable variation in some isolates emphasizes the extraordinary genetic conservation of this species and the high clonality of the population. High selection pressure within the intracellular niche and an evolutionary bottleneck as a consequence of the adaptation to the human or animal host might have contributed to the stability of the genome.
SNPs unique for the indicated reference isolate(s)
Number of SNPs
Isolate(s), identified by SNP
In silico analysis of SNPs
To identify SNPs within the four sequenced isolates, we implemented a Java program which performed the tasks of calculating sequence alignments and extracting the SNPs from the alignments. The genomic sequences and annotations were downloaded from the Genbank database . Starting from the CWL-029 genome and comparing it to the other genomes, syntenic protein coding regions, RNA genes and intergenic regions were identified by Smith-Waterman alignments. The gap-free part of the alignment with the best hit was used for SNP identification, if the fraction of mismatches did not exceed 5%. SNPs were identified from these alignments and classified according to their position on the CWL-029 genome to estimate their specificity. Protein coding SNPs were additionally classified into synonymous and non-synonymous SNPs. The distribution of SNPs around the genome (figure 1) was displayed by GenomeDiagram . The phylogeny of the reference strains was conducted with MEGA 3.1  using the neighbour-joining method with 1000 bootstrap replicates, and distances calculated using the number of different SNPs. The position of polymorphic genes and gene products within metabolic pathways were demonstrated by KEGG pathway assignments to the genes of the CWL-029 genome. They were manually extracted from the pathway maps at the KEGG website . In order to estimate the odds that, for the CWL-029 genome, a single nucleotide mutation would result in a sSNPs or nsSNPs, the number of synonymous and nonsynonymous sites in all coding DNA sequences of the CWL-029 genomes was calculated. Based on the bacterial genetic code we classified each coding nucleotide by the consequence of the three possible mutations. Thus, we identified 842596 sites at which nsSNPs and 245642 sites at which sSNPs could occur.
PCR and sequencing
According to the locations of sSNPs 189 loci evenly distributed around the whole genome were chosen for further analysis. The adjacent up- and downstream region (approx. 500 bp) of each SNP was amplified by PCR (primers provided in additional file 3) according to common protocols. 8 μl of the PCR product were digested with 0.3 U SAP (shrimp alkaline phosphatase) and 1.5 U ExoI (both Amersham Biosciences, Freiburg, Germany) for 15 min. The reaction was stopped by heating for 15 min at 72°C. The sequencing reaction was performed using 1 μl of ABI PRISM™ BigDye™ (Applied Biosystems, Foster City, USA), a 1.6 μM concentration of each sequencing primer (Eurogentec, Seraing, Belgium) and 2 μl of digested PCR product. The reaction conditions were 96°C for 10 s followed by 25 cycles of 95°C for 1 min, and 60°C for 175 s.
Data and phylogenetic analysis
Chromatograms were analysed with Bio Edit  and SNPs were analysed using a Local Blast comparison with the CWL-029 isolate . An additional 43 sSNPs were identified within the sequenced genome fragments (~0.1 mega bases). The SNP data of totally 232 loci were concatenated, resulting in one character string (nucleotide sequence) for each strain. The raw data are provided as additional file 4. Concatenated alignments are given as additional file 5. Phylogenetic analyses were conducted with MEGA 3.1  using the neighbour-joining method with 1000 bootstrap replicates and distances calculated using the number of different SNPs. In order to root the C. pneumoniae phylogeny, we blasted 50 gene fragments that we had already sequenced from the DC-9 and Koala isolates. The highest homologies were obtained for Chlamydophila abortus, an endemic ruminant bacteria that colonizes the placenta. Using this approach we finally obtained a 2 kb alignment of orthologous gene fragments. The dating of different lineages split was calculated based on an E. coli molecular clock rate. Comparisons of homologous protein-coding regions from E. coli and S. enterica indicated an average rate of sequence divergence at synonymous sites of 0.9% per millions years, a value extrapolated from a 120–160 million years divergence time [31, 32].
To determine whether some C. pneumoniae populations underwent recent population expansions, we calculated mismatch distributions and compared these to predicted distributions from models of population expansion. For expanding populations, we converted the parameter tau (τ; calculated from the mismatch distribution) to estimate the time of the expansion (t) using the equation τ = 2μt, where μ is the neutral mutation rate for the locus. The confidence intervals of τ were calculated using a parametric bootstrap approach . Mismatch distributions and τ were calculated in ARLEQUIN 3.0 .
synonymous single nucleotide polymorphism
non-synonymous single nucleotide polymorphism
million years ago
horizontal gene transfer
long branch attraction
polymorphic blood mononuclear cells
J.G. was supported by a grant of the Deutsche Forschungsgemeinschaft (GI 344/2-1) and the University of Lübeck (N23-2003). The excellent technical assistance of K. Roßdeutscher, U. Panknin and S. Pätzmann is gratefully acknowledged. We thank all the colleagues (table 1) who provided us with the isolates.
- Stephens RS: Chlamydial Evolution: A billion years and counting. Chlamydial Infections. Prodeedings of the Tenth International Symposium on Human Chlamydial Infections. Edited by: Schachter, J. 2002, Antalya, TurkeyGoogle Scholar
- Abdelrahman YM, Belland RJ: The chlamydial development cycle. FEMS Microbiol Rev. 2005, 29: 949-59. 10.1016/j.femsre.2005.03.002.PubMedView ArticleGoogle Scholar
- Ortutay C, Gaspari Z, Toth G, Jager E, Vida G, Orosz L, Vellai T: Speciation in Chlamydia: genomewide phylogenetic analyses identified a reliable set of acquired genes. J Mol Evol. 2003, 57: 672-80. 10.1007/s00239-003-2517-3.PubMedView ArticleGoogle Scholar
- Read TD, Brunham RC, Shen C, Gill SR, Heidelberg JF, White O, Hickey EK, Peterson J, Utterback T, Berry : Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39. Nucleic Acids Res. 2000, 28: 1397-406. 10.1093/nar/28.6.1397.PubMed CentralPubMedView ArticleGoogle Scholar
- Shirai M, Hirakawa H, Kimoto M, Tabuchi M, Kishi F, Ouchi K, Shiba T, Ishii K, Hattori M, Kuhara S, Nakazawa T: Comparison of whole genome sequences of Chlamydia pneumoniae J138 from Japan and CWL029 from USA. Nucleic Acids Res. 2000, 28: 2311-4. 10.1093/nar/28.12.2311.PubMed CentralPubMedView ArticleGoogle Scholar
- Kalman S, Mitchell W, Marathe R, Lammel C, Fan J, Hyman RW, Olinger L, Grimwood J, Davis RW, Stephens RS: Comparative genomes of Chlamydia pneumoniae and C. trachomatis. Nat Genet. 1999, 21: 385-9. 10.1038/7716.PubMedView ArticleGoogle Scholar
- Grayston JT: Chlamydia pneumoniae, strain TWAR pneumonia. Annu Rev Med. 1992, 43: 317-23. 10.1146/annurev.me.43.020192.001533.PubMedView ArticleGoogle Scholar
- Kuo CC, Grayston JT: Chlamydia spp. strain TWAR: A newly recognized organism associated with atypical pneumonia and other respiratory infections. Clin Microbial Newsletter. 1988, 10: 137-144. 10.1016/0196-4399(88)90070-0.View ArticleGoogle Scholar
- Saikku P: The epidemiology and significance of Chlamydia pneumoniae. J Infect. 1992, 27-34. 10.1016/0163-4453(92)91913-V. highly variable proteins in the Chlamydophila pneumoniae genome. Nucleic Acids Res 2002, 30:4351–4360., Suppl 1
- Maass M, Bartels C, Kruger S, Krause E, Engel PM, Dalhoff K: Endovascular presence of Chlamydia pneumoniae DNA is a generalized phenomenon in atherosclerotic vascular disease. Atherosclerosis. 1998, S25-30. 10.1016/S0021-9150(98)00117-8. Suppl 1
- Gieffers J, Fullgraf H, Jahn J, Klinger M, Dalhoff K, Katus HA, Solbach W, Maass M: Chlamydia pneumoniae infection in circulating human monocytes is refractory to antibiotic treatment. Circulation. 2001, 103: 351-6.PubMedView ArticleGoogle Scholar
- Balin BJ, Gerard H, Arking EJ, Appelt DM, Branigan PJ, Abrams JT, Whittum-Hudson JA, Hudson AP: Identification and localization of Chlamydia pneumoniae in the Alzheimer's brain. Med Microbiol Immunol (Berl). 1998, 187: 23-42. 10.1007/s004300050071.View ArticleGoogle Scholar
- Gieffers J, Pohl D, Treib J, Dittmann R, Stephan C, Klotz K, Hanefeld F, Solbach W, Haass A, Maass M: Presence of Chlamydia pneumoniae DNA in the cerebral spinal fluid is a common phenomenon in a variety of neurological diseases and not restricted to multiple sclerosis. Ann Neurol. 2001, 49: 585-9. 10.1002/ana.1020.PubMedView ArticleGoogle Scholar
- Rocha EPC, Pradillon O, Bui H, Sayada C, Denamur E: A new family of highly variable proteins in the Chlamydophila pneumoniae genome. Nucleic Acids Res. 2002, 30: 4351-4360. 10.1093/nar/gkf571.PubMed CentralPubMedView ArticleGoogle Scholar
- Gutacker MM, Smoot JC, Migliaccio CA, Ricklefs SM, Hua S, Cousins DV, Graviss EA, Shashkina E, Kreiswirth BN, Musser JM: Genome-wide analysis of synonymous single nucleotide polymorphisms in Mycobacterium tuberculosis complex organisms: resolution of genetic relationships among closely related microbial strains. Genetics. 2002, 162: 1533-43.PubMed CentralPubMedGoogle Scholar
- Pearson T, Busch JD, Ravel J, Read TD, Rhoton SD, U'Ren JM, Simonson TS, Kachur SM, Leadem RR, Cardon ML: Phylogenetic discovery bias in Bacillus anthracis using single-nucleotide polymorphisms from whole-genome sequencing. Proc Natl Acad Sci USA. 2004, 101: 13536-41. 10.1073/pnas.0403844101.PubMed CentralPubMedView ArticleGoogle Scholar
- Harpending HC, Batzer MA, Gurven M, Jorde LB, Rogers AR, Sherry ST: Genetic traces of ancient demography. Proc Natl Acad Sci USA. 1998, 95: 1961-1967. 10.1073/pnas.95.4.1961.PubMed CentralPubMedView ArticleGoogle Scholar
- Gieffers J, Durling L, Ouellette SP, Rupp J, Maass M, Byrne GI, Caldwell HD, Belland RJ: Genotypic differences in the Chlamydia pneumoniae tyrP locus related to vascular tropism and pathogenicity. J Infect Dis. 2003, 188: 1085-93. 10.1086/378692.PubMedView ArticleGoogle Scholar
- Everett KD, Andersen AA: The ribosomal intergenic spacer and domain I of the 23S rRNA gene are phylogenetic markers for Chlamydia spp. Int J Syst Bacteriol. 1997, 47: 461-73.PubMedView ArticleGoogle Scholar
- Pettersson B, Andersson A, Leitner T, Olsvik O, Uhlén M, Storey C, Black CM: Evolutionary relationships among members of the genus Chlamydia based on 16S ribosomal DNA analysis. J Bacteriol. 1997, 179: 4195-205.PubMed CentralPubMedGoogle Scholar
- Sreevatsan S, Pan X, Stockbauer KE, Connell ND, Kreiswirth BN, Whittam TS, Musser JM: Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination. Proc Natl Acad Sci USA. 1997, 94: 9869-74. 10.1073/pnas.94.18.9869.PubMed CentralPubMedView ArticleGoogle Scholar
- Brosch R, Gordon SV, Marmiesse M, Brodin P, Buchrieser C, Eiglmeier K, Garnier T, Gutierrez C, Hewinson G, Kremer : A new evolutionary scenario for the Mycobacterium tuberculosis complex. Proc Natl Acad Sci USA. 2002, 99: 3684-9. 10.1073/pnas.052548299.PubMed CentralPubMedView ArticleGoogle Scholar
- Carlson JH, Hughes S, Hogan D, Cieplak G, Sturdevant DE, McClarty G, Caldwell HD, Belland RJ: Polymorphisms in the Chlamydia trachomatis cytotoxin locus associated with ocular and genital isolates. Infect Immun. 2004, 72: 7063-72. 10.1128/IAI.72.12.7063-7072.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Maass M, Harig U: Evaluation of culture conditions used for isolation of Chlamydia pneumoniae. Am J Clin Pathol. 1995, 103: 141-8.PubMedGoogle Scholar
- NCBI gene bank. [http://www.ncbi.nlm.nih.gov/genomes]
- Pritchard L, White JA, Birch PR, Toth IK: GenomeDiagram: a python package for the visualization of large-scale genomic data. Bioinformatics. 2006, 22: 616-7. 10.1093/bioinformatics/btk021.PubMedView ArticleGoogle Scholar
- Mega Software 3.1. [http://www.megasoftware.net]
- KEGG pathways. [http://www.genome.jp/kegg/]
- Bio Edit. [http://www.mbio.ncsu.edu/BioEdit/bioedit.html]
- Local Blast. [http://www.stdgen.lanl.gov]
- Ochman H, Wilson AC: Evolution in bacteria: evidence for a universal substitution rate in cellular genomes. J Mol Evol. 1987, 26: 74-86. 10.1007/BF02111283.PubMedView ArticleGoogle Scholar
- Wirth T, Falush D, Lan R, Colles F, Mensa P, Wieler LH, Karch H, Reeves PR, Maiden MCJ, Ochman H, Achtman M: Sex and virulence in Escherichia coli : an evolutionary perspective. Mol Mic. 2006, 60: 1136-1151. 10.1111/j.1365-2958.2006.05172.x.View ArticleGoogle Scholar
- Rogers AR: Genetic evidence for a Pleistocene population explosion. Evolution. 1995, 49: 608-615. 10.2307/2410314.View ArticleGoogle Scholar
- Schneider S, Excoffier L: Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA. Genetics. 1999, 3: 1079-1089.Google Scholar
- Excoffier L, Laval G, Schneider S: Evolutionary Bioinformatics Online. 2005Google Scholar