1 Institut für Hygiene und Mikrobiologie der Bayerischen Julius-Maximilians-Universität, 97080 Würzburg, Germany
2 Institut für Hygiene, 48149 Münster, Germany
3 Poliklinik für Parodontologie der Westfälischen Wilhelms-Universität, 48149 Münster, Germany
Correspondence
Herbert Schmidt
Herbert.Schmidt{at}mailbox.tu-dresden.de
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The nucleotide sequences of the fragments of all Porphyromonas gingivalis genes and their alleles (190 sequences) reported in this paper have been entered into the EMBL/GenBank databases and have received continuous accession numbers from AJ555632 to AJ555821.
Present address: Institut für Medizinische Mikrobiologie und Krankenhaushygiene, Medizinische Hochschule Hannover, Carl-Neuberg-Straße 1, 30625 Hannover, Germany.
Present address: Institut für Medizinische Mikrobiologie und Hygiene, Medizinische Fakultät Carl Gustav Carus, Technische Universität Dresden, Fetscherstraße 74, 01307 Dresden, Germany.
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
By AP-PCR, 73 P. gingivalis strains could be divided into 23 or 45 genotypes, depending on the primers used. Using REA, it has been demonstrated that a single type generally colonized individuals and strains of different volunteers were genetically different (Loos et al., 1993; Menard et al., 1994
). Analysis of 100 P. gingivalis strains from humans and animals revealed 78 different multilocus enzyme electrophoresis (MLEE) types (Loos et al., 1993
). Using ribotyping, 25 isolates from six subjects could be divided into six clonal types, each of which was unique for every subject (Socransky & Martin, 1992
). The consensus of these studies is that there is considerable heterogeneity among P. gingivalis isolates but the intra-individual heterogeneity is low. A couple of studies have been performed to relate clonal types with increased pathogenicity. In none of these studies could a clear association with a given clonal type and the periodontal status be made.
To better understand intraspecies variation, we performed multilocus sequence analysis with 10 chromosomal genes. This method, introduced in 1998, has since been used for phylogenetic analysis of many other bacterial pathogens such as Neisseria meningitidis, Campylobacter jejuni, Helicobacter pylori and Escherichia coli (Maiden et al., 1998; Suerbaum et al., 1998
, 2001
; Suerbaum, 2000
; Falush et al., 2003
). Molecular analysis of housekeeping genes has led to a better understanding of the mode of genetic variation within a bacterial species. Here, we apply this method in a slightly modified form to the periodontal pathogen P. gingivalis.
![]() |
METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
PCR.
A single bacterial colony was suspended in 50 µl of a 0·7 % (w/v) sodium chloride solution and this suspension was subsequently diluted 1 : 100, again in 0·7 % sodium chloride. Five microlitres of this dilution were used as template for PCRs. PCRs were performed in MicroAmp reaction tubes (Applied Biosystems) in a GeneAmp System 9600 thermal cycler (Applied Biosystems) in a total volume of 50 µl, consisting of 5 µl bacterial suspension, 30 pmol each primer, 200 µM each dNTP, 5 µl 10-fold-concentrated Taq DNA polymerase synthesis buffer, 3 µl 25 mM MgCl2 stock solution and 2 units AmpliTaq DNA polymerase (Applied Biosystems). PCR primers and conditions are described in Table 1. Analysis of PCR products was performed by horizontal agarose gel electrophoresis in 0·72 % (w/v) agarose gels in Tris/Borate/EDTA (TBE) buffer. Gels were stained for 10 min in a solution of 1 µg ethidium bromide ml-1, washed for 10 min in deionized water and analysed and photographed under UV light (312 nm).
|
PCR products were sequenced from both strands as described previously (Wittstock et al., 2000), using the same primers as applied for PCRs. Four microlitres of the purified PCR product were used as template for sequencing.
Selection and processing of target genes.
The 10 chromosomal P. gingivalis genes ef-tu (accession no. AB035462), ftsQ (accession no. D84504; Ansai et al., 1995), hagB (accession no. Z35494; Progulske-Fox et al., 1995
), gpdxJ (accession no. X97228), pepO (accession no. AB010440; Awano et al., 1999
), mcmA (accession no. L30136; Jackson et al., 1995
), dnaK (AB015879; Yoshida et al., 1999
), recA (accession no. AF064682; Fletcher et al., 1997
), pga (accession no. X95938) and nah (accession no. X78979) were selected for the investigation described here. Since the genome sequence of P. gingivalis had not been published when this study was performed, we used genes that have been published in the framework of the Porphyromonas gingivalis genome project (http://www.pgingivalis.org/). We primarily selected housekeeping genes and other essential genes. We amplified internal fragments of approximately 500 bp of those genes (Table 1
) and sequenced both strands. From these sequences we were able to use fragments between 310 and 420 bp for analyses.
Sequence analyses.
Sequence electropherograms were analysed with the Lasergene software package (DNAStar). DNA sequences were edited and aligned with Bio Edit, version 4.8.10 (http://jwbrown.mbio.ncsu.edu/BioEdit/bioedit.html) (Hall, 1999) and converted into MEGA and NEXUS files with START [Sequence Type Analysis and Recombinational Tests; keith.jolley{at}ceid.ox.ac.uk (http://outbreak.ceid.ox.ac.uk/software.htm)]. Phylogenetic trees were compiled with MEGA 2.1 (http://www.megasoftware.net/) using the unweighted pair group method with arithmetic averages (UPGMA) (Kumar et al., 1994
). Genetic distances were calculated with CLUSTAL W included in Bio Edit (Thompson et al., 1994
). Split decomposition analysis (Bandelt & Dress, 1992
; Huson, 1998
) was performed with SPLITSTREE 2.0 (http://bibiserv.techfak.uni-bielefeld.de/splits/). Homoplasy tests were performed with HOMOPLASY (Suerbaum et al., 1998
). Numerical analysis of polymorphic sites was done with DNASP (Rozas & Rozas, 1999
) and graphical analysis with HAPPLOT (http://www.shigatox.net/stec/programs/happlot/). A stand-alone UNIX version of LIAN 3.1 to calculate the standardized index of association (
) was kindly provided by Bernhard Haubold, Lion Bioscience AG, Heidelberg, Germany (Haubold et al., 1998
). The program STEPHENS (http://www.shigatox.net/stec/programs/stephens) was used to detect clustered polymorphisms in the sequences analysed.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
The allelic numbers were used to calculate a distance matrix. This matrix was converted into NEXUS format and graphically displayed with SPLITSTREE (Fig. 1). The relationship of P. gingivalis isolates based on multilocus analysis is depicted as a star-like structure with rays of different lengths. This star phylogeny is also consistent with a recombinational population structure and is consistent with the observation that each strain carries its own allelic combination. Isolates Eli47/GDR8 und Eli55/Fr19 are more closely related than the other isolates.
|
|
Linkage equilibrium
Two loci are in linkage equilibrium if genotype frequencies at one locus are independent of genotype frequencies at a second locus, otherwise the two loci are in linkage disequilibrium. Linkage disequilibrium can arise from physical linkage, genetic drift and selection on multilocus genotypes. Linkage equilibrium is a situation that should exist in a population undisturbed by selection, migration, etc., in which all possible combinations of linked genes should be present at equal frequency. Recombination (crossover) shifts a population towards linkage equilibrium; selection shifts a population towards linkage disequilibrium.
Using the program LIAN, we tested the null hypothesis of linkage equilibrium for the P. gingivalis dataset. Linkage equilibrium is characterized by statistical independence of alleles at all loci of multilocus data. LIAN tests for this independent assortment by computing , which is a measure of the degree of linkage in multilocus datasets.
is a function of the rate of recombination and is zero when all alleles are in linkage equilibrium. That means probably that the distribution of alleles occurs independently from each other. The
of the P. gingivalis genes analysed was 0·0898. This low (
) value is indicative of extensive recombination. For comparison, the
value of Campylobacter jejuni was 0·256 (Suerbaum et al., 2001
). Simple IA values, which are not standardized, have been obtained by analysis of MLEE data of Neisseria gonorrhoeae, N. meningitidis and H. influenzae and were 0·004, 1·96 and 5·4, respectively (Maynard Smith et al., 1993
). Using the equation for calculation of
(described in the documentation for LIAN 3.1 and by Haubold et al., 1998
) and the information on the number of loci given in that paper (Maynard Smith et al., 1993
), we could calculate their
values (adjusted to the number of loci investigated) to 0·005, 0·14 and 0·337, respectively. The
value calculated for P. gingivalis is similar to that reported for N. gonorrhoeae and therefore supports our estimation that the genes investigated in P. gingivalis are close to linkage equilibrium.
Analysis of clustered polymorphisms
To get more evidence for the influence of recombination on the observed genetic heterogeneity of P. gingivalis genes, we investigated the alleles for clustered polymorphisms, a hallmark of intragenic recombination. Visual inspection of the distribution of polymorphic sites showed multiple examples where alleles appeared to be mosaics of different clusters of polymorphic sites. The statistical significance of clusters was analysed with the test described by Stephens (1985), implemented in the program STEPHENS by T. Whittam (available at http://www.shigatox.net/stec/programs/stephens/). The program classifies each polymorphic site by how it groups (or partitions) the sequences and determines the probability that the sites supporting a given partition are clustered. For our dataset, significant clustering of polymorphic sites was detected in hagB, pga, ftsQ, dnaK and ef-tu. The results of this analysis were concordant with the clusters detected by visual inspection of graphical outputs of HAPPLOT (Fig. 3
).
|
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The remarkable differences in progression and severity of periodontal disease have been linked to variability in the virulence of strains within the pathogenic species. Studies in animals and humans have provided evidence that certain strains of P. gingivalis are associated with periodontal disease, whereas other strains seem less virulent (Amano et al., 1999; Baker et al., 2000
; Nakagawa et al., 2002
). For example, it has been shown that two FimA genotypes (Type II and Type IV) of P. gingivalis were strongly associated with periodontal disease, whereas others were weakly or not associated with periodontal disease (Amano et al., 1999
, 2000
). Thus, the differentiation of disease-contributing strains by variations of the fimA gene has been considered an alternative approach for risk assessment and treatment planning of P. gingivalis-related periodontal diseases.
Molecular epidemiology attempts bacterial classification on the basis of natural genetic variation. During the last decades, molecular epidemiology has been used as an adjunct and increasingly as a replacement for classical methods.
Three types of population structures are known in bacteria: clonal, panmictic and epidemic. These types are represented by Salmonella enterica, N. gonorrhoeae and N. meningitidis, respectively. In a clonal species, genetic changes would result from the sequential accumulation of mutations after descent from a common ancestor. However, only a few species are truly clonal and horizontal gene transfer has occurred even in S. enterica (Maynard Smith et al., 1993). Panmictic populations may be so variable that identical strains are only found among isolates from direct contacts. Even panmictic populations can contain clonal groupings or geographical specialization. This was not the case in our study, but the number of isolates analysed here is too low to make any conclusions about the occurrence of such clonal groupings (Achtman, 2001
).
The analysis of the population structure of P. gingivalis presented here shows a substantial extent of recombination in P. gingivalis. In a recent study, Frandsen et al. (2001) analysed four housekeeping genes from 57 P. gingivalis strains. They found 41 different sequence types. In that study, the IA was calculated as 0·143 for the 41 sequence types and 0·206 when all 57 isolates were included. When we use the equation for the calculation of
with the information given in the paper of Frandsen et al. (2001)
, which includes the number of loci investigated,
for all 57 isolates is 0·068 and therefore similar to the
value calculated in this paper for P. gingivalis. The authors suggested a non-clonal population structure characterized by recombination (Frandsen et al., 2001
). Our study, based on 10 chromosomal gene fragments from 19 P. gingivalis isolates confirms the importance of recombination in P. gingivalis by multiple population genetic methods.
Recombination has also been described to occur during the chronic colonization of H. pylori in the gastric mucosa (Falush et al., 2001). In this study, the authors could also show that the panmictic population structure of H. pylori may result from frequent recombination during mixed colonization by unrelated strains.
Where does recombination between P. gingivalis strains occur? Periodontal pockets of patients with periodontitis may be a favourable environment for horizontal gene transfer and recombination between P. gingivalis strains. Although it has been shown that P. gingivalis is able to receive foreign DNA by transformation in vitro (Yoshimoto et al., 1993), there are only a few studies available demonstrating indirect evidence for horizontal gene transfer in vivo (Hanley et al., 1999
). Frequent recombination between strains is likely to facilitate the spread of favourable traits. Future experiments on recombination of periodontal pathogens in vivo may broaden our understanding of mechanisms and meaning of genetic variation in this important group of bacteria.
![]() |
ACKNOWLEDGEMENTS |
---|
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Albandar, J. M. & Kingman, A. (1999). Gingival recession, gingival bleeding, and dental calculus in adults 30 years of age and older in the United States, 19881994. J Periodontol 70, 3043.[Medline]
Amano, A., Nakagawa, I., Kataoka, K., Morisaki, I. & Hamada, S. (1999). Distribution of Porphyromonas gingivalis strains with fimA genotypes in periodontitis patients. J Clin Microbiol 37, 14261430.
Amano, A., Kuboniwa, M., Nakagawa, I., Akiyama, S., Morisaki, I. & Hamada, S. (2000). Prevalence of specific genotypes of Porphyromonas gingivalis fimA and periodontal health status. J Dent Res 79, 16641668.[Abstract]
Ansai, T., Yamashita, Y., Awano, S., Shibata, Y., Wachi, M., Nagai, K. & Takehara, T. (1995). A murC gene in Porphyromonas gingivalis 381. Microbiology 141, 20472052.[Abstract]
Awano, S., Ansai, T., Mochizuki, H., Yu, W., Tanzawa, K., Turner, A. J. & Takehara, T. (1999). Sequencing, expression and biochemical characterization of the Porphyromonas gingivalis pepO gene encoding a protein homologous to human endothelin-converting enzyme. FEBS Lett 460, 139144.[CrossRef][Medline]
Baelum, V., Fejerskov, O. & Manji, F. (1988). Periodontal diseases in adult Kenyans. J Clin Periodontol 15, 445452.[Medline]
Baelum, V., Luan, W. M., Chen, X. & Fejerskov, O. (1997a). Predictors of destructive periodontal disease incidence and progression in adult and elderly Chinese. Community Dent Oral Epidemiol 25, 265272.[Medline]
Baelum, V., Luan, W. M., Chen, X. & Fejerskov, O. (1997b). Predictors of tooth loss over 10 years in adult and elderly Chinese. Community Dent Oral Epidemiol 25, 204210.[Medline]
Baker, P. J., Dixon, M., Evans, R. T. & Roopenian, D. C. (2000). Heterogeneity of Porphyromonas gingivalis strains in the induction of alveolar bone loss in mice. Oral Microbiol Immunol 15, 2732.[CrossRef][Medline]
Bandelt, H. J. & Dress, A. W. (1992). Split decomposition: a new and useful approach to phylogenetic analysis of distance data. Mol Phylogenet Evol 1, 242252.[Medline]
Brown, L. J. & Loe, H. (1993). Prevalence, extent, severity and progression of periodontal disease. Periodontol 2000 2, 5771.
Chen, C. & Slots, J. (1994). Clonal analysis of Porphyromonas gingivalis by the arbitrarily primed polymerase chain reaction. Oral Microbiol Immunol 9, 99103.[Medline]
Cutler, C. W., Kalmar, J. R. & Genco, C. A. (1995). Pathogenic strategies of the oral anaerobe, Porphyromonas gingivalis. Trends Microbiol 3, 4551.[CrossRef][Medline]
Falush, D., Kraft, C., Taylor, N. S., Correa, P., Fox, J. G., Achtman, M. & Suerbaum, S. (2001). Recombination and mutation during long-term gastric colonization by Helicobacter pylori: estimates of clock rates, recombination size, and minimal age. Proc Natl Acad Sci U S A 98, 1505615061.
Falush, D., Wirth, T., Linz, B. & 15 other authors (2003). Traces of human migrations in Helicobacter pylori populations. Science 299, 15821585.
Fletcher, H. M., Morgan, R. M. & Macrina, F. L. (1997). Nucleotide sequence of the Porphyromonas gingivalis W83 recA homolog and construction of a recA-deficient mutant. Infect Immun 65, 45924597.[Abstract]
Frandsen, E. V., Poulsen, K., Curtis, M. A. & Kilian, M. (2001). Evidence of recombination in Porphyromonas gingivalis and random distribution of putative virulence markers. Infect Immun 69, 44794485.
Hall, T. A. (1999). Bio Edit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 41, 9598.
Hanley, S. A., Aduse-Opoku, J. & Curtis, M. A. (1999). A 55-kilodalton immunodominant antigen of Porphyromonas gingivalis W50 has arisen via horizontal gene transfer. Infect Immun 67, 11571171.
Haubold, B., Travisano, M., Rainey, P. B. & Hudson, R. R. (1998). Detecting linkage disequilibrium in bacterial populations. Genetics 150, 13411348.
Huson, D. H. (1998). SplitsTree: analyzing and visualizing evolutionary data. Bioinformatics 14, 6873.[Abstract]
Jackson, C. A., Kirszbaum, L., Dashper, S. & Reynolds, E. C. (1995). Cloning, expression and sequence analysis of the genes encoding the heterodimeric methylmalonyl-CoA mutase of Porphyromonas gingivalis W50. Gene 167, 127132.[CrossRef][Medline]
Kato, T., Takahashi, N. & Kuramitsu, H. K. (1992). Sequence analysis and characterization of the Porphyromonas gingivalis prtC gene, which expresses a novel collagenase activity. J Bacteriol 174, 38893895.[Abstract]
Kumar, S., Tamura, K. & Nei, M. (1994). MEGA: Molecular Evolutionary Genetics Analysis software for microcomputers. Comput Appl Biosci 10, 189191.[Abstract]
Kuramitsu, H. K., Yoneda, M. & Madden, T. (1995). Proteases and collagenases of Porphyromonas gingivalis. Adv Dent Res 9, 3740.[Abstract]
Lamont, R. J. & Jenkinson, H. F. (1998). Life below the gum line: pathogenic mechanisms of Porphyromonas gingivalis. Microbiol Mol Biol Rev 62, 12441263.
Loos, B. G., Mayrand, D., Genco, R. J. & Dickinson, D. P. (1990). Genetic heterogeneity of Porphyromonas (Bacteroides) gingivalis by genomic DNA fingerprinting. J Dent Res 69, 14881493.[Abstract]
Loos, B. G., Dyer, D. W., Whittam, T. S. & Selander, R. K. (1993). Genetic structure of populations of Porphyromonas gingivalis associated with periodontitis and other oral infections. Infect Immun 61, 204212.[Abstract]
Maiden, M. C., Bygraves, J. A., Feil, E. & 10 other authors (1998). Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci U S A 95, 31403145.
Maynard Smith, J. & Smith, N. H. (1998). Detecting recombination from gene trees. Mol Biol Evol 15, 590599.[Abstract]
Maynard Smith, J., Smith, N. H., O'Rorke, M. & Spratt, B. G. (1993). How clonal are bacteria? Proc Natl Acad Sci U S A 90, 43844388.[Abstract]
Menard, C., Brousseau, R. & Mouton, C. (1992). Application of polymerase chain reaction with arbitrary primer (AP-PCR) to strain identification of Porphyromonas (Bacteroides) gingivalis. FEMS Microbiol Lett 74, 163168.[Medline]
Menard, C., Gosselin, P., Duhaime, J. F. & Mouton, C. (1994). Polymerase chain reaction using arbitrary primer for the design and construction of a DNA probe specific for Porphyromonas gingivalis. Res Microbiol 145, 595602.[CrossRef][Medline]
Nakagawa, I., Amano, A., Kuboniwa, M., Nakamura, T., Kawabata, S. & Hamada, S. (2002). Functional differences among FimA variants of Porphyromonas gingivalis and their effects on adhesion to and invasion of human epithelial cells. Infect Immun 70, 277285.
Ozmeric, N., Preus, N. R. & Olsen, I. (2000). Genetic diversity of Porphyromonas gingivalis and its possible importance to pathogenicity. Acta Odontol Scand 58, 183187.[CrossRef][Medline]
Pilot, T. & Barmes, D. E. (1987). An update on periodontal conditions in adults, measured by CPITN. Int Dent J 37, 169172.[Medline]
Progulske-Fox, A., Tumwasorn, S., Lepine, G., Whitlock, J., Savett, D., Ferretti, J. J. & Banas, J. A. (1995). The cloning, expression and sequence analysis of a second Porphyromonas gingivalis gene that codes for a protein involved in hemagglutination. Oral Microbiol Immunol 10, 311318.[Medline]
Rozas, J. & Rozas, R. (1999). DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics 15, 174175.
Sheiham, A. (1991). Public health aspects of periodontal diseases in Europe. J Clin Periodontol 18, 362369.[Medline]
Sheiham, A. (1997). Is the chemical prevention of gingivitis necessary to prevent severe periodontitis? Periodontol 2000 15, 1524.
Socransky, S. S. & Martin, L. (1992). Ribotyping of Porphyromonas gingivalis and Campylobacter rectus. J Dent Res 71, 246.
Stephens, J. C. (1985). Statistical methods of DNA sequence analysis: detection of intragenic recombination or gene conversion. Mol Biol Evol 2, 539556.[Abstract]
Suerbaum, S. (2000). Genetic variability within Helicobacter pylori. Int J Med Microbiol 290, 175181.[Medline]
Suerbaum, S. & Achtman, M. (1999). Evolution of Helicobacter pylori: the role of recombination. Trends Microbiol 7, 182184.[CrossRef][Medline]
Suerbaum, S., Smith, J. M., Bapumia, K., Morelli, G., Smith, N. H., Kunstmann, E., Dyrek, I. & Achtman, M. (1998). Free recombination within Helicobacter pylori. Proc Natl Acad Sci U S A 95, 1261912624.
Suerbaum, S., Lohrengel, M., Sonnevend, A., Ruberg, F. & Kist, M. (2001). Allelic diversity and recombination in Campylobacter jejuni. J Bacteriol 183, 25532559.
Thompson, J. D., Higgins, D. G. & Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22, 46734680.[Abstract]
Wittstock, M. (1999). Untersuchungen zur PrtC-Kollagenase des Parodontitis-Erregers Porphyromonas gingivalis. PhD thesis, Institute of Hygiene and Microbiology, University of Würzburg, Germany.
Wittstock, M., Flemmig, T. F., Schmidt, H., Mutters, R. & Karch, H. (1996). Serodiagnosis of Porphyromonas gingivalis infection by immunoblot analysis with recombinant collagenase. J Clin Microbiol 34, 24112413.[Abstract]
Wittstock, M., Schmidt, H., Flemmig, T. F. & Karch, H. (2000). Heterogeneity of the prtC gene of Porphyromonas gingivalis. Oral Microbiol Immunol 15, 3339.[CrossRef][Medline]
Yoshida, A., Nakano, Y., Yamashita, Y., Oho, T., Shibata, Y., Ohishi, M. & Koga, T. (1999). A novel dnaK operon from Porphyromonas gingivalis. FEBS Lett 446, 287291.[CrossRef][Medline]
Yoshimoto, H., Takahashi, Y., Hamada, N. & Umemoto, T. (1993). Genetic transformation of Porphyromonas gingivalis by electroporation. Oral Microbiol Immunol 8, 208212.[Medline]
Received 28 January 2003;
revised 2 May 2003;
accepted 2 June 2003.
HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
J MED MICROBIOL | ALL SGM JOURNALS |