1,729
Views
13
CrossRef citations to date
0
Altmetric
Review Article

Genetic exchange and reassignment in Porphyromonas gingivalis

, &
Article: 1457373 | Received 24 Aug 2017, Accepted 19 Mar 2018, Published online: 12 Apr 2018

ABSTRACT

Porphyromonas gingivalis is considered a keystone pathogen in adult periodontitis but has also been associated with systemic diseases. It has a myriad of virulence factors that differ between strains. Genetic exchange and intracellular genome rearrangements may be responsible for the variability in the virulence of P. gingivalis. The present review discusses how the exchange of alleles can convert this bacterium from commensalistic to pathogenic and potentially shapes the host-microbe environment from homeostasis to dysbiosis. It is likely that genotypes of P. gingivalis with increased pathogenic adaptations may spread in the human population with features acquired from a common pool of alleles. The exact molecular mechanisms that trigger this exchange are so far unknown but they may be elicited by environmental pressure.

Introduction

Porphyromonas gingivalis is considered a keystone pathogen in periodontitis [Citation1] but it may also be involved in several periodontitis-related systemic diseases [Citation2]. When the prevalence of periodontitis was estimated in Americans >65 years, almost two-thirds (62.3%) had one or more periodontal sites with ≥5 mm of clinical attachment loss and almost 50% had at least one site with probing pocket depth of ≥4 mm [Citation3]. The older adult population is growing rapidly in the Western world and it is estimated that by 2040 the number of adults ≥65 years will be about 50% in the USA [Citation3]. Their oral health status will be very important in order to maintain an adequate quality of life.

P. gingivalis expresses a myriad of virulence factors [Citation4]. Since genetic exchange may cause diversity in specific virulence factors in P. gingivalis, the present review will give an overview on how genetic exchange and rearrangement may occur in this anaerobic bacterium.

Natural competence and conjugation

It is well known that prokaryotes can respond dynamically to changing environments by genetic exchange [Citation5]. P. gingivalis exchanges chromosomal DNA between strains by natural competence and conjugation [Citation6Citation8]. By using antibiotic resistance markers, it was demonstrated that the strains ATCC 33277T and W83 of P. gingivalis were able to exchange chromosomal DNA resulting in offspring, which carried DNA from both donor and recipient [Citation6]. The chimeras thereby created had altered phenotypic features as demonstrated through biofilm formation [Citation9]. Natural competence affected by the competence protein ComF was found to be the predominant form of DNA transfer in P. gingivalis while transfer by a conjugation-like approach occurred at a lower frequency [Citation6]. Subsequently, natural competence mechanisms have been detected in multiple P. gingivalis strains and the DNA uptake is not sensitive to the DNA source or modification status [Citation6]. The major source for horizontal DNA transfer and allelic exchange between strains was found to be extracellular DNA present in P. gingivalis biofilms. Moreover, P. gingivalis may even scavenge DNA from other species as a nutritional source [Citation6].

A high level of reticulation was seen in the phylogenetic networks constituted by 23 strains of P. gingivalis [Citation8]. This indicated extensive horizontal gene transfer between the strains. There were variants in the major virulence factor of the Kgp proteinase, the Kgp C-terminal cleaved adhesion domain and the surface proteins FimA, FimCDE, Mfa l, RagAB, Tpr and PrtT. It was concluded that P. gingivalis uses domain rearrangements and genetic exchange to generate diversity in specific surface virulence factors. When the whole-genome sequences between ATCC 33277 and W83 were determined, extensive rearrangements were found between the two strains, possibly induced by different mobile elements (insertion sequences, miniature inverted-repeat transposable elements, and conjugative transposons) [Citation10]. In ATCC 33277, many of the genes showed higher similarity to the genes of other bacterial species such as Bacteroides fragilis than to genes in W83. This suggested that they had been introduced to ATCC 33277 by horizontal gene transfer [Citation10].

The specific combination of genes in the genome of the dominant P. gingivalis strain may vary with the disease state and the host response. Probably exchange of DNA in biofilms by transformation is important to the ecology and persistence of P. gingivalis in the challenging environment of the periodontal pocket and to evasion of host immune defences and immunization strategies [Citation6,Citation9]. P. gingivalis appears to have both a core genome [Citation7,Citation11] and a flexible, dispensable genome. Sampling of DNA from the flexible genome could imply adaptation for the bacterium to the changing environment [Citation7], which is supported by the genetic diversity of P. gingivalis.

Fimbriae types provide selective advantage in periodontal disease

Although genetic variations can be seen in many virulence factors of P. gingivalis, most interest has been devoted to the major fimbrial protein FimA [Citation12Citation14]. The gene encoding the protein subunit that is assembled to form the surface fimbriae of P. gingivalis is the fimA gene. Six alleles are found in the fimA gene: types I–V and Ib. Type II and IV fimA have significantly greater adhesive and invasive capabilities and are associated with more cytotoxic strains compared to other fimA type clones [Citation15]. While type Ib, II and IV are most frequently associated with periodontitis [Citation14Citation18], type I and III are less inflammatory and are associated with strains more likely to form localized self-limiting infections [Citation19]. It has been suggested that the fimA type II allele confers a selective advantage in the advanced disease process thus ensuring its distribution throughout the P. gingivalis population [Citation20].

Another type of fimbria expressed by P. gingivalis is Mfa 1. These minor or short fimbriae are mainly composed of polymers of the Mfa 1 protein and are encoded by the mfa 1 gene [Citation13]. Both FimA and Mfa 1 fimbriae enable bacteria to adhere to surfaces, and are important in biofilm formation [Citation21]. A genetic association between the fimA genotype and the presence of both the mfa 1 and the 53k gene has been suggested [Citation22]. The 53k gene encodes the 53K protein (53-kDa major outer membrane protein) of P. gingivalis. This gene was located at a locus that corresponded to the mfa 1 gene [Citation22]. The outer membrane protein has been reported as a major fimbrilin of a novel-type fimbrium [Citation22]. Generally, P. gingivalis expresses either 53K or Mfa 1 fimbriae in addition to FimA fimbriae [Citation22]. However, strains with fimA genotypes III and V express 53K fimbriae almost exclusively, while fimA genotypes I and II tend to express Mfa 1 fimbriae. The downstream gene product Mfa 2 has a role as anchor for the Mfa 1 fimbriae and as regulator of the Mfa 1 filament length [Citation23]. It is likely that the distinct short fimbria-molecules found in different strains as well as long fimbriae both influence initiation and progression of periodontal disease [Citation24].

A survey of the fimbrilin encoding genes among the 35 currently available P. gingivalis genomic sequences revealed that most of the P. gingivalis strains possess both FimA and Mfa 1 fimbrilin encoding genes (). The only exceptions were strain MP4-504 and W83, of which only the FimA fimbrilin encoding gene was found in their genomes. Interestingly, strain W50 has been shown to be very similar to W83. However, the Mfa 1 encoding gene was only detected in W50, even though the sequence is unfinished. The number of genomes encoding type I, type Ib, type II, type III and type IV FimA were 9, 3, 17, 1 and 5, respectively. The number of genomes encoding type I and type II Mfa 1 were 19 and 14, respectively. The only genome encoding the Type III FimA gene was the strain ATCC 49417 with the protein sequence ID SJL32224.1 (). SJL32224.1 is identical to a previously cloned and derived Type III FimA sequence (Genbank Accession: AHG31000.1) from the same strain [Citation7], thus confirming the presence of this type of FimA in this strain. There are only two other P. gingivalis Type III FimA sequences in the NCBI Protein Database – BAA04627 and Q51826, representing a slightly different version of Type III FimA but so far has only been found in strain 6/26 [Citation25].

Table 1. Occurrences of transposase, fimbrilin and CRISPR features in 35 P. gingivalis genomes*

Figure 1. Major and minor fimbrilin proteins encoded in the 35 P. gingivalis genomic sequences. All the fimbrilin proteins were detected in the following manner: The keyword FimA/Mfa1/fimbrilin were used to select from the proteins predicted in the 35 genomes, which were annotated by the NCBI prokaryotic genome annotation pipeline [Citation54]. The keywords were also used to search for P. gingivalis FimA from the NCBI protein databases. These “seed” proteins were then searched against all the proteins in the 35 genomes to identify additional fimbrilin candidates. These potential candidates matched with the seed proteins at ≥ 90% sequence identity and ≥ 90% length coverage. The additional candidates were combined with other annotated fimbrilin proteins and aligned with the “mafftℍ program [Citation55]. The aligned sequences were used to build a phylogenetic tree with the “FastTree” program [Citation56]. For each fimbrilin shown in the tree, the Genbank accession number was printed first, followed with the strain label (in blue), the original functional annotation (by NCBI), and the size of the protein sequences in amino acids (red)

Figure 1. Major and minor fimbrilin proteins encoded in the 35 P. gingivalis genomic sequences. All the fimbrilin proteins were detected in the following manner: The keyword FimA/Mfa1/fimbrilin were used to select from the proteins predicted in the 35 genomes, which were annotated by the NCBI prokaryotic genome annotation pipeline [Citation54]. The keywords were also used to search for P. gingivalis FimA from the NCBI protein databases. These “seed” proteins were then searched against all the proteins in the 35 genomes to identify additional fimbrilin candidates. These potential candidates matched with the seed proteins at ≥ 90% sequence identity and ≥ 90% length coverage. The additional candidates were combined with other annotated fimbrilin proteins and aligned with the “mafftℍ program [Citation55]. The aligned sequences were used to build a phylogenetic tree with the “FastTree” program [Citation56]. For each fimbrilin shown in the tree, the Genbank accession number was printed first, followed with the strain label (in blue), the original functional annotation (by NCBI), and the size of the protein sequences in amino acids (red)

Exchange of allele types in fimbriae

Kerr et al. [Citation7] reported that P. gingivalis exchanged fimbriae allele types I and IV into four distinct strain backgrounds via natural competence. A complete exchange of the entire fimA allele was detected but the rate of exchange varied between different strain backgrounds. Also, genetic exchange within other regions of the genetic locus of fimbriae was seen. Exchange of fimbriae allele type caused marked phenotypic changes, e.g. in the quantity of fimbriae produced, membrane blebbing and autoaggregation. When the type I allele was replaced with the type III or IV allele, increased invasion of gingival fibroblasts were seen compared to that achieved with the isogenic parent strain. By using plasmid vectors, Kato et al. [Citation26] substituted fimA of ATCC 33277 with type II fimA, and that of OMZ314 (type II strain) with type I fimA. Replacement of type I fimA with type II increased bacterial adhesion/invasion to epithelial cells, whereas substitution with type I fimA resulted in reduced efficiency. Similarly, Kerr et al. [Citation7] found that substitution of the type IV allele from strain W83 into ATCC 33277 transformed the 33277 strain into a strain phenotypically more similar to the W83 strain. Such strategies may be used by pathogenic strains to modify their interactions with the host and thereby change the homeostatic balance in their favour.

Clonality

The expression of heterogenic virulence properties in P. gingivalis may depend on clonal diversity. In clonal populations of bacteria, genetic diversity primarily arises by point mutations. Point mutations within open reading frames that are limited in size can occur in P. gingivalis at frequencies of 10–9, similar to other bacterial species such as Escherichia coli [Citation20].

In non-clonal bacterial populations, however, there is a spectrum of phenotypes from commensal to pathogenic phenotype [Citation27] and a number of allele combinations that can be related to disease [Citation20]. Some allele combinations may be more virulent than others and it is reasonable to think that allele combinations of P. gingivalis differ in health and disease. Non-clonal pathogens are common in bacteria associated with chronic diseases. Based on the P. gingivalis MLST (Multi Locus Sequence Typing) database (https://pubmlst.org/pgingivalis/) the index (IA) for estimating the degree of association between alleles at different loci for representative sequence types, was calculated to be 0.3704 suggesting a weakly clonal population structure [Citation15,Citation28]. By using the standardized index of association (ISA), which takes into account the number of genes investigated, thereby providing a better estimate of the linkage equilibrium than does IA, the value was 0.062. Clonal populations give a score of one, while mono-clonal populations yield zero. This supported a study by Frandsen et al. [Citation29] who found that P. gingivalis has a non-clonal population structure with frequent recombination. These figures also yielded support to the idea of P. gingivalis consisting mainly of mono-clonal strains established through bacterial transfer of genes in vivo, so-called bacterial sex [Citation30]. Consequently, the genetic variation in P. gingivalis strains is considerable suggesting that exchange of genes occurs to improve the chances of P. gingivalis to survive in vivo. This variability can be extensive and comprise entire genes as well as operons [Citation20].

Pathogenicity islands

Many of the features suggesting lateral gene transfer in P. gingivalis are demonstrated through the presence of pathogenicity islands, so called because they contain assemblies of genes for virulence that cause disease [Citation31]. Ranging in size from 10 to 200 kb, pathogenicity islands often carry genes that encode integrases and transposases involved in DNA mobility, and they can be associated with tRNA genes, which are favoured sites for integration of foreign DNA [Citation32]. This was demonstrated by Chen et al. [Citation32] when comparative whole-genome analysis of virulent (W83) and avirulent (ATCC 33277) strains of P. gingivalis was made. In W83, one cluster of genes was associated with two paralogous regions of chromosomes with low G + C content that contained conserved and species-specific hypothetical genes, transposons, insertion sequences and integrases and was located close to tRNA genes, having many of the characteristics of pathogenicity islands acquired by lateral gene transfer.

Rag locus

One pathogenicity island in P. gingivalis is the rag locus [Citation33]. The rag locus of strain W50 encodes RagA (receptor antigen gene A), a predicted TonB-dependent receptor protein, and RagB, a lipoprotein that constitutes an immunodominant outer membrane antigen [Citation33Citation35]. They probably act together on the surface of the organisms to help active transport, mediated through the periplasmic protein TonB or constitute part of a signal transduction system. The rag locus of P. gingivalis may have arisen via horizontal gene transfer from Bacteroides [Citation33,Citation36]. Supporting that horizontal gene transfer had taken place were the low G + C content of the locus, the association with motility elements and a restricted distribution in the species [Citation34]. In P. gingivalis, the exogenous rag locus promoted expression of virulence genes or activation of dormant genes [Citation35].

PCR demonstrated that rag+ P. gingivalis strains were more frequently present in deep periodontal pockets than in shallow pockets of patients with periodontitis [Citation33,Citation35]. While highly virulent P. gingivalis strains always carried the rag locus, avirulent strains did not. In mice, rag-1 isolates more likely caused disease than other rag alleles (rag-2–4) [Citation34]. The rag locus was also seen in a small portion of sites where P. gingivalis was not detected [Citation36], and Frandsen et al. [Citation29] found a random distribution of the ragB locus in 132 P. gingivalis strains.

Reca and pepo housekeeping genes

When MLST was performed of clinical P. gingivalis isolates from 15 subjects with ‘refractory’ periodontitis, not one but several sequence types were detected for most individual periodontal pockets. The allelic variation in sequence types indicated recombination of the recA and pepO housekeeping genes in diseased sites [Citation37]. This could have been the result of recombination between different clones of P. gingivalis within the periodontal pocket.

Conjugative transposons and transposons

Conjugative transposons (Ctns) represent a subgroup of integrative and conjugative elements characterized by their ability to move from one bacteria to another through a process requiring cell-cell contact [Citation38]. They can excise themselves from the genome where they are integrated, transfer themselves to a recipient by conjugation and integrate into the recipient genome [Citation39]. Thus CtnPg1, which was the first complete Ctn found in P. gingivalis, was able to excise itself from chromosomal DNA using an integrase (Int; PGN-0094) encoded in CtnPg1 and produce a circular intermediate [Citation10]. The transfer of CtnPg1 requires a RecA function in the recipient cell. It carried a complete set of genes for conjugative transfer and integration. When CtnPg1 was transferred from ATCC 33277 to W83, chromosomal gene transfer occurred in addition to the CtnPg1 transfer itself.

Many Ctns carry genes for antibiotic resistance, providing e.g. tetracycline resistance and erythromycin resistance [Citation40] and transfer of Ctns often contributes to spread of antibiotic resistance genes. Ctns have been studied most frequently in Bacteroides where many species have large self-transmissible Ctns. ATCC 33277 of P. gingivalis carried four copies of Ctns that were not present in W83 [Citation10]. Several genes were moderately homogenous to the sequences of genes of Ctns from Bacteroides species. However, the sequences of two copies of novel composite transposons detected in the ATCC 33777 genome were identical [Citation10]. Thus the CtnPg1 family differs from the Bacteroides Ctn family and also from the Tn916/Tn1545 family. The CtnPg1 family is distributed among both oral and intestinal bacteria where it can play an important role in horizontal gene transfer, e.g. of antibiotic resistance genes to anaerobic bacteria.

shows the number of complete and partial proteins that are homologous to the transposases, ranging from 8 (strains Ando, F0185, and SDJ5) to 103 (strain A7436). However, the lower numbers are from unfinished draft sequences and are likely under-estimated. In general, P. gingivalis is a species with a high number of transposases, hence the genetic exchange and genomic locus rearrangement in and between strains are expected to be frequent.

Crispr-Cas

The clustered regularly interspaced short palindromic repeats (CRISPR) and their associated proteins (Cas) are adaptive immune systems in prokaryotes present in most Bacteria and Archaea. They provide adaptive immunity against foreign elements such as bacteriophages/viruses, plasmids and transposons. At least 95% of clinical strains from P. gingivalis harbours CRISPR arrays [Citation41]. The presence of transposases and CRISPR-associated sequences in P. gingivalis [Citation11] suggested that some parts of the P. gingivalis genome are mobile so that rearrangement or horizontal gene transfer from one strain to another can take place [Citation20]. In the P. gingivalis W83 genome, all four CRISPR regions were transcribed and one of them was active against dsDNA in vivo [Citation41]. Probably the CRISPR-Cas system has developed as a protective system that verifies and selects DNA entering bacterial cells, helping them to regulate horizontal gene transfer [Citation41].

summarizes the occurrence of the CRISPR-Cas proteins and CRISPR elements detected in the genomes of 35 P. gingivalis strains. Strain 13_1 encodes 19 Cas homologous proteins and thus may contain at least 3 CRISPR systems (a Cas operon contains between 6 to 8 ORFs). On the other hand, JSVI SC001, an environmental strain, contains no detectable Cas proteins, even though 3 CRISPR-like DNA sequence structures were detected. The number of CRISPR systems in the genome may reflect the difficulty of incorporating foreign DNA while the number of direct repeats (DR) detected in the CRISPR elements (last column in ) may indicate the CRISPR activity in the past for the genome and how often the strain had encountered the invasion of foreign DNA.

Proteins and whole genome nucleotide sequences

Protein coding sequences for the virulent W83 strain comprised several genes that may be associated with higher virulence of the strain, including glycosyltransferase, a protein required for capsular polysaccharide biosynthesis, sensor histidine kinase, surface antigen and thiol protease [Citation10].

Phylogenomic comparison based on shared proteins and whole-genome nucleotide sequences consistently showed two groups among 19 P. gingivalis strains examined with closely related members: one consisted of ATCC 33277, 381 and HG66, the other of W83, W50 and A7436 [Citation11]. Among them, ATCC 33277 and W83 were mono-clonal but fell in different groups. Obviously, other factors than just clonality affect the virulence of P. gingivalis strains. Comparative functional genetics showed strain difference in gingipains, attachment, heme, gene mobility, transposases, capsules, CRISPR and phage [Citation11].

shows the homologous DNA segments (syntenies) among the 10 completed P. gingivalis genomic sequences. Strains ATCC 33277 and 381 have been shown to be very close and showed the highest similarity between the syntenies detected (i.e. the coloured rectangles in show nearly identical sizes and genomic placement between these two strains). This is as expected. The order and direction of the syntenies indicated the rearrangement of genomic loci, possibly the result of transposition activities in the past. The blue and red arrowhead lines depicted in exemplify the inversion of two regions of syntenies across all genomes. Interestingly, in the genome of AJW4, an insertion event can be observed in the synteny region indicated by the dash part of the blue arrowhead line. These syntenies represent the core genomic loci that are characteristics for this species while the gaps may have been introduced at different times in different environments.

Figure 2. Syntenic analysis of 9 P. gingivalis genomic sequences. Of the 35 currently available P. gingivalis genomics sequences, eight are considered finished with one final successfully assembled contig. The sequence of JCVI SC001 appears to have a one-contig circular sequence under the Genbank Accession number CM001843, however it is a pseudo-contig generated by ordering the 284 unassembled contigs using the sequence of strain TDC60 as the template. Thus, the syntenies, illustrated as the same color across different genomes, of JCVI SC001 appear to be in similar order of those in TDC60 and may not be the true order in this genome. Syntenies were detected with the “MUAVE” program (version 2.4.0) [Citation57] and were illustrated as rectangles of the same colors across all genomes. Two groups of syntenies were indicated with blue and red arrows underneath to exemplify the different rearrangement and inversion events in different genomes for this species

Figure 2. Syntenic analysis of 9 P. gingivalis genomic sequences. Of the 35 currently available P. gingivalis genomics sequences, eight are considered finished with one final successfully assembled contig. The sequence of JCVI SC001 appears to have a one-contig circular sequence under the Genbank Accession number CM001843, however it is a pseudo-contig generated by ordering the 284 unassembled contigs using the sequence of strain TDC60 as the template. Thus, the syntenies, illustrated as the same color across different genomes, of JCVI SC001 appear to be in similar order of those in TDC60 and may not be the true order in this genome. Syntenies were detected with the “MUAVE” program (version 2.4.0) [Citation57] and were illustrated as rectangles of the same colors across all genomes. Two groups of syntenies were indicated with blue and red arrows underneath to exemplify the different rearrangement and inversion events in different genomes for this species

Insertion segments

Watanabe et al. [Citation42] claimed that the high interspecies diversity in P. gingivalis was caused by frequent insertion sequences (ISs). Frandsen et al. [Citation29] found IS1598 in 69% of 132 P. gingivalis strains examined. Altogether, P. gingivalis has seven different IS elements named ISPg 1–7 [Citation20]. ISPg5 and other IS elements were proposed to contribute to strain diversity in P. gingivalis [Citation36] and they may be used for fingerprinting of strains. More than one ISPg element exists in most strains of P. gingivalis and they can probably be active in vivo [Citation20]. Examples of IS activity related to pathogenicity include transposition [Citation43], activation of antibiotic resistance genes [Citation43Citation45], movement and dissemination of virulence genes [Citation46,Citation47], and their inactivation [Citation48Citation50]. IS transposition can be limited by CRISPR interference as CRISPR limits both IS transposition and the introduction of DNA from other P. gingivalis cells [Citation51]. Thus specific domain rearrangements and genetic exchange to generate diversity in surface virulence factors [Citation8] cannot alone be held responsible for the diversity in virulence, although they probably play an important role.

Plasmid mobilization

CTns are able to mobilize plasmids and other transposons in Bacteroides. Natural plasmids have not been detected in P. (previously Bacteroides) gingivalis. However, plasmid and transposon DNA can be conjugally transferred into P. gingivalis [Citation52]. This does not apply to P. gingivalis strains in general since only some stains could transfer CTns and other accessory mobile elements, while others could not [Citation20]. Thus, plasmid transfer was not possible in W83, possibly because its traAQ locus lacks a traP homolog. This gene is required for plasmid conjugation in Bacteroides [Citation20].

Concluding remarks

As reviewed, P. gingivalis is a species of low clonality that apparently is the result of many accessory genes that contribute to the genomic fluidity, such as the CTns and IS elements. On the other hand, genes that limit the uptake of foreign DNA, e.g. the CRISPR-Cas, are also prevalent in this species. This ‘checks and balances’ governed by the two systems may also be another key to the role of this species in shaping the oral environment in between homeostasis and dysbiosis. It seems possible that genetic transfer and rearrangement makes P. gingivalis adapt to and persist in the host environment and that minor changes in the genetic content can cause significant changes in its phenotype. This may be important for the bacterium’s transfer from commensalism to pathogenicity and provides the strains with unique adaptability in vivo. Even restricted changes in a single gene may result in extensive biological changes important to the survival of P. gingivalis in the host.

Disclosure statement

No potential conflict of interest was reported by the authors.

Additional information

Funding

The European Commission [Grant Number FP7-HEALTH-306029 ‘TRIGGER’] supported this work (IO).

Notes on contributors

Ingar Olsen

Ingar Olsen is professor emeritus and guest researcher at Department of Oral Biology, Faculty of Dentistry, University of Oslo. Senior Research Investigator, Department of Microbiology, Forsyth Institute, Cambridge, MA.Tsute Chen is associate research investigator at The Forsyth Institute, Cambridge, MA, USA and instructor at the Department of Oral Medicine, Infection and Immunity, Harvard School of Dental Medicine, Boston, MA, USA.Gena D. Tribble is associate professor at the Department of Periodontics and Dental Hygiene, School of Dentistry, University of Texas Health Science Center at Houston, Houston, TX, USA

References

  • Hajishengallis G, Darveau RP, Curtis MA. The keystone-pathogen hypothesis. Nat Rev Microbiol. 2012;10:1–9.
  • Olsen I. From the acta prize lecture 2014: the periodontal-systemic connection seen from a microbiological standpoint. Acta Odontol Scand. 2015;73(8):563.
  • Eke PI, Wei L, Borgnakke WS, et al. Periodontitis prevalence in adults ≥ 65 years of age, in the USA. Periodontol. 2000. 2016;72(1):76–95. DOI:10.1111/prd.12145
  • Zenobia C, Hajishengallis G. Porphyromonas gingivalis virulence factors involved in subversion of leukocytes and microbial dysbiosis. Virulence. 2015;6(3):236–243.
  • Bryant J, Chewapreecha C, Bentley SD. Developing insights into the mechanisms of evolution of bacterial pathogens from whole-genome sequences. Future Microbiol. 2012;11:1283–1296. DOI:10.2217/fmb.12.108
  • Tribble GD, Rigney TW, Dao DH, et al. Natural competence is a major mechanism for horizontal DNA transfer in the oral pathogen Porphyromonas gingivalis. MBio. 2012;3(1). pii):e00231–11.
  • Kerr JE, Abramian JR, Dao DH, et al. Genetic exchange of fimbrial alleles exemplifies the adaptive virulence strategy of Porphyromonas gingivalis. PLoS One. 2014;9:e91696.
  • Dashper SG, Mitchell HL, Seers CA, et al. Porphyromonas gingivalis uses specific domain rearrangements and allelic exchange to generate diversity in surface virulence factors. Front Microbiol. 2017;8:48.
  • Tribble GD, Lamont GJ, Progulske-Fox A, et al. Conjugal transfer of chromosomal DNA contributes to genetic variation in the oral pathogen Porphyromonas gingivalis. J Bacteriol. 2007;189(17):6382–6388.
  • Naito M, Hirakawa H, Yamashita A, et al. Determination of the genome sequence of Porphyromonas gingivalis strain ATCC 33277 and genomic comparison with strain W83 revealed extensive genome rearrangements in P. gingivalis. DNA Res. 2008;15(4):215–225.
  • Chen T, Siddiqui H, Olsen I. In silico comparison of 19 Porphyromonas gingivalis strains in genomics, phylogenetics, phylogenomics and functional genomics. Front Cell Infect Microbiol. 2017;7:28.
  • Nagano K, Hasegawa Y, Abiko Y, et al. Porphyromonas gingivalis FimA fimbriae: fimbrial assembly by fimA alone in the fim gene cluster and differential antigenicity among fimA genotypes. PLoS One. 2012;7:e43722.
  • Nagano K, Abiko Y, Yoshida Y, et al. Genetic and antigenic analyses of Porphyromonas gingivalis fimA fimbriae. Mol Oral Microbiol. 2013;28:392–403.
  • Nakagawa I, Amano A, Ohara-Nemoto Y, et al. Identification of a new variant of fimA gene of Porphyromonas gingivalis and its distribution in adults and disabled populations with periodontitis. J Periodontal Res. 2002;37(6):425–432.
  • Enersen M, Olsen I, Kvalheim Ø, et al. fimA genotypes and multilocus sequence types of Porphyromonas gingivalis from patients with periodontitis. J Clin Microbiol. 2008;46:31–42.
  • Missailidis CG, Umeda JE, Ota-Tsuzuki C, et al. Distribution of fimA genotypes of Porphyromonas gingivalis in subjects with various periodontal conditions. Oral Microbiol Immunol. 2004;19:224–229.
  • Zhao L, Wu YF, Meng S, et al. Prevalence of fimA genotypes of Porphyromonas gingivalis and periodontal health status in Chinese adults. J Periodontal Res. 2007;42:511–517.
  • Zheng C, Wu J, Xie H. Differential expression and adherence of Porphyromonas gingivalis FimA genotypes. Mol Oral Microbiol. 2011;26:388–395.
  • Nakano K, Kuboniwa M, Nakagawa I, et al. Comparison of inflammatory changes caused by Porphyromonas gingivalis with distinct fimA genotypes in a mouse abscess model. Oral Microbiol Immunol. 2004;19:205–209.
  • Tribble GD, Kerr JE, Wang BY. Genetic diversity in the oral pathogen Porphyromonas gingivalis: molecular mechanisms and biological consequences. Future Microbiol. 2013;8(5):607–620.
  • Amano A, Nakagawa I, Okahashi N, et al. Variations of Porphyromonas gingivalis fimbriae in relation to microbial pathogenesis. J Periodontal Res. 2004;39:136–142.
  • Nagano K, Hasegawa Y, Yoshida Y, et al. A major fimbrilin variant of Mfa1 fimbriae in Porphyromonas gingivalis. J Dent Res. 2015;94(8):1143–1148.
  • Hasegawa Y, Iwami J, Sato K, et al. Anchoring and length regulation of Porphyromonas gingivalis Mfa1 fimbriae by the downstream gene product Mfa 2. Microbiology. 2009;155(Pt 10):3333–3347.
  • Enersen M, Nakano K, Amano A. Porphyromonas gingivalis fimbriae. J Oral Microbiol. 2013;5. DOI:10.3402/jom.v5i0.20265.
  • Fujiwara T, Morishima S, Takahashi I, et al. Molecular cloning and sequencing of the fimbrilin gene of Porphyromonas gingivalis strains and characterization of recombinant proteins. Biochem Biophys Res Commun. 1993;197(1):241–247.
  • Kato T, Kawai S, Nakano K, et al. Virulence of Porphyromonas gingivalis is altered by substitution of fimbria gene with different genotype. Cell Microbiol. 2007;9(3):753–765.
  • Spratt BG, Maiden MC. Bacterial population genetics, evolution and epidemiology. Philos Trans R Soc Lond B Biol Sci. 1999;354(1384):701–710.
  • Enersen M. Porphyromonas gingivalis: a clonal pathogen?: diversities in housekeeping genes and the major fimbriae gene. J Oral Microbiol. 2011;3. DOI:10.3402/jom.v3i0.8487.
  • Frandsen EV, Poulsen K, Curtis MA, et al. Evidence of recombination in Porphyromonas gingivalis and random distribution of putative virulence markers. Infect Immun. 2001;69(7):4479–4485.
  • Olsen I, Tribble GD, Fiehn N-E, et al. Bacterial sex in dental plaque. J Oral Microbiol. 2013;5. DOI:10.3402/jom.v5i0.20736.
  • Hacker J, Kaper JB. Pathogenicity islands and the evolution of microbes. Annu Rev Microbiol. 2000;54:641–679.
  • Chen T, Hosogi Y, Nishikawa K, et al. Comparative whole-genome analysis of virulent and avirulent strains of Porphyromonas gingivalis. J Bacteriol. 2004;186(16):5473–5479.
  • Curtis MA, Hanley SA, Aduse-Opoku J. The rag locus of Porphyromonas gingivalis: a novel pathogenicity island. J Periodontal Res. 1999;34(7):400–405.
  • Hall LM, Fawell SC, Shi X, et al. Sequence diversity and antigenic variation at the rag locus of Porphyromonas gingivalis. Infect Immun. 2005;73(7):4253–4262.
  • Su Z, Kong F, Wang S, et al. The rag locus of Porphyromonas gingivalis might arise from Bacteroides via horizontal gene transfer. Eur J Clin Microbiol Infect Dis. 2010;29(4):429–437.
  • Califano JV, Kitten T, Lewis JP, et al. Characterization of Porphyromonas gingivalis insertion sequence-like element ISPg5. Infect Immun. 2000;68(9):5247–5253.
  • Enersen M, Olsen I, Caugant DA. Genetic diversity of Porphyromonas gingivalis isolates recovered from single “refractory” periodontitis sites. Appl Environ Microbiol. 2008;74(18):5817–5821.
  • Wozniak RA, Waldor MK. Integrative and conjugative elements: mosaic mobile genetic elements enabling dynamic lateral gene flow. Nat Rev Microbiol. 2010;8(8):552–563.
  • Naito M, Sato K, Shoji M, et al. Characterization of the Porphyromonas gingivalis conjugative transposon CTnPg1: determination of the integration site and the genes essential for conjugal transfer. Microbiology. 2011;157(Pt7):2022–2032.
  • Clewell DB, Flannagan SE. The conjugative transposons of Gram-positive bacteria. In: Clewell DB, editor. Bacterial conjugation. New York: Plenum Press; 1993. p. 369–393.
  • Burmistrz M, Dudek B, Staniec D, et al. Functional analysis of Porphyromonas gingivalis W83 CRISPR-Cas systems. J Bacteriol. 2015;197(16):2631–2641.
  • Watanabe T, Nozawa T, Aikawa C, et al. CRISPR regulation of intraspecies diversification by limiting IS transposition and intercellular recombination. Genome Biol Evol. 2013;5:1099–1114.
  • Lee KY, Hopkins JD, Syvanen M. Direct involvement of IS26 in an antibiotic resistance operon. J Bacteriol. 1990;172(6):3229–3236.
  • Podglajen I, Breuil J, Collatz E. Insertion of a novel DNA sequence, 1S1186, upstream of the silent carbapenemase gene cfiA, promotes expression of carbapenem resistance in clinical isolates of Bacteroides fragilis. Mol Microbiol. 1994;12(1):105–114.
  • Podglajen I, Breuil J, Casin I, et al. Genotypic identification of two groups within the species Bacteroides fragilis by ribotyping and by analysis of PCR-generated fragment patterns and insertion sequence content. J Bacteriol. 1995;177(18):5270–5275.
  • Terai A, Baba K, Shirai H, et al. Evidence for insertion sequence-mediated spread of the thermostable direct hemolysin gene among Vibrio species. J Bacteriol. 1991;173(16):5036–5046.
  • Garcia MI, Labigne A, Le Bouguenec C. Nucleotide sequence of the afimbrial-adhesin-encoding afa-3 gene cluster and its translocation via flanking IS1 insertion sequences. J Bacteriol. 1994;176(24):7601–7613.
  • Ou JT, Baron LS, Rubin FA, et al. Specific insertion and deletion of insertion sequence 1-like DNA element causes the reversible expression of the virulence capsular antigen Vi of Citrobacter freundii in Escherichia coli. Proc Natl Acad Sci U S A. 1988;85(12):4402–4405.
  • Simonet M, Riot B, Fortineau N, et al. Invasin production by Yersinia pestis is abolished by insertion of an IS200-like element within the inv gene. Infect Immun. 1996;64(1):375–379.
  • Lewis JP, Macrina FL. IS195, an insertion sequence-like element associated with protease genes in Porphyromonas gingivalis. Infect Immun. 1998;66(7):3035–3042.
  • Berg DE. Transposon Tn5. In: Berg DE, Howe MM, editors. Mobile DNA. Washington, D.C.:  American Society for Microbiology (ASM); 1989. p. 185–210.
  • Dyer DW, Bilalis G, Michel JH, et al. Conjugal transfer of plasmid and transposon DNA from Escherichia coli into Porphyromonas gingivalis. Biochem Biophys Res Commun. 1992;186(2):1012–1019.
  • Tatusova T, DiCuccio M, Badretdin A, et al. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 2016;44(14):6614–6624.
  • Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–780.
  • Price MN, Dehal PS, Arkin AP. FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5(3):e9490.
  • Darling AE, Mau B, Perna NT. ProgressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One. 2010;5(6): e11147. DOI:10.1371/journal.pone.0011147.
  • Grissa I, Vergnaud G, Pourcel C. The CRISPRdb database and tools to display CRISPRs and to generate dictionaries of spacers and repeats. BMC Bioinformatics. 2007;8:172.