Katterinne N. Mendez1†
Anilei Hoare2†
Cristopher Soto2
Isaac Bugueño2
Marcela Olivera2Claudio Meneses3,4
Jose Manuel Pérez-Donoso1
Eduardo Castro-Nallar1*
Denisse Bravo2*- 1Center for Bioinformatics and Integrative Biology, Facultad de Ciencias de la Vida, Universidad Andrés Bello, Santiago, Chile
- 2Laboratorio de Microbiología Oral, Facultad de Odontología, Universidad de Chile, Santiago, Chile
- 3Centro de Biotecnología Vegetal, Facultad de Ciencias de la Vida, Universidad Andrés Bello, Santiago, Chile
- 4FONDAP Center for Genome Regulation, Santiago, Chile
Porphyromonas gingivalis has been extensively associated with both the onset and progression of periodontitis. We previously isolated and characterized two P. gingivalis strains, one from a patient exhibiting severe chronic periodontitis (CP3) and another from a periodontally healthy individual (H3). We previously showed that CP3 and H3 exhibit differences in virulence since H3 showed a lower resistance to cationic peptides compared with CP3, and a lower ability to induce proliferation in gingival epithelial cells. Here, we aimed to determine whether differences in virulence between these two strains are associated with the presence or absence of specific genes encoding virulence factors. We sequenced the whole genomes of both P. gingivalis CP3 and H3 and conducted a comparative analysis regarding P. gingivalis virulence genetic determinants. To do so, we performed a homology search of predicted protein sequences in CP3 and H3 genomes against the most characterized virulence genes for P. gingivalis available in the literature. In addition, we performed a genomic comparison of CP3 and H3 with all the 62 genomes of P. gingivalis found in NCBI's RefSeq database. This approach allowed us to determine the evolutionary relationships of CP3 and H3 with other virulent and avirulent strains; and additionally, to detect variability in presence/absence of virulence genes among P. gingivalis genomes. Our results show genetic variability in the hemagglutinin genes. While CP3 possesses one copy of hagA and two of hagC, H3 has no hagA and only one copy of hagC. Experimentally, this finding is related to lower in vitro hemmaglutination ability of H3 compared to CP3. Moreover, while CP3 encodes a gene for a major fimbrium subunit FimA type 4 (CP3_00160), H3 possess a FimA type 1 (H3_01400). Such genetic differences are in agreement with both lower biofilm formation ability and less intracellular invasion to oral epithelial cells exhibited by H3, compared with the virulent strain CP3. Therefore, here we provide new results on the genome sequences, comparative genomics analyses, and phenotypic analyses of two P. gingivalis strains. The genomics comparison of these two strains with the other 62 genomes included in the analysis provided relevant results regarding genetic determinants and their association with P. gingivalis virulence.
Introduction
Porphyromonas gingivalis is a Gram-negative anaerobic bacterium that significantly contributes to the onset and progression of periodontitis (Hajishengallis, 2015), an oral disease affecting 30–50% worldwide adult population (Petersen and Ogawa, 2012). Besides, P. gingivalis has also been associated with increased risk of developing other prevalent diseases, such as cardiovascular diseases, rheumatoid arthritis, metabolic diseases, cancer, and Alzheimer's disease, among others (Babic et al., 2015; Kumar, 2017).
P. gingivalis has been proposed to be a keystone pathogen in the etiology of periodontitis, since it can promote changes in homeostasis of the commensal microbiota, contributing to the ecological successions occurring in the subgingival area that are associated with inflammation and tissue destruction in in vivo models (Hajishengallis et al., 2011; Hajishengallis, 2015). Nonetheless, P. gingivalis is present in the subgingival biofilm of both diseased subjects and periodontally healthy individuals, although the frequency of detection is higher in the first group (Lau et al., 2004; Kumawat et al., 2016; Kulkarni et al., 2018). From an ecologic point of view, this species is classified as a pathobiont, being naturally part of the subgingival microbiota in periodontally healthy individuals but, under certain environmental conditions, contributing to or exacerbating the disease state (Hajishengallis and Lamont, 2012; Cugini et al., 2013).
Interestingly, P. gingivalis has not only been found in higher loads in the subgingival biofilm in periodontitis subjects compared to healthy individuals, their virulent properties vary among strains isolated from each condition. In previous studies, we reported that clinical isolates of P. gingivalis obtained from healthy (H1, H2, and H3) and periodontitis subjects (CP1, CP2, CP3) show differences in their virulence. For example, strains from periodontitis patients exhibit a higher resistance to antimicrobial peptides (Díaz et al., 2015) and higher ability to decrease apoptosis of host cells and therefore, to persist intracellularly (Soto et al., 2016). Among them, CP3 and H3 respectively, showed the highest and lowest resistance to antimicrobial peptides and ability to both decrease apoptosis and persist intracellularly, respectively. In agreement with these results, CP3 was obtained from a subject showing the most severe clinical signs of periodontitis, while the subject that donated the sample containing H3, showed to be periodontally healthy.
Several virulence factors contribute to differences in the bacterial virulence of P. gingivalis, in particular the Fimbriae and hemagglutinins. Based on the diversity of the fimA gene, fimbriae are classified into six genotypes (I-V and Ib) (Hamada et al., 1994; Nakagawa et al., 2000), being type II fimbriae, followed by type IV, the most frequently found in P. gingivalis strains isolated from periodontitis patients (Amano et al., 2000, 2004; Missailidis et al., 2004). In contrast, healthy adults carry strains coding type I fimA (Amano et al., 2000). This is not surprising since strains expressing type II FimA have been associated with both higher abilities to adhere and invade epithelial cells (Nakagawa et al., 2002; Amano et al., 2004), as well as an increased production of biofilm (Kuboniwa et al., 2009). Particularly, while strains harboring type II and IV FimA produce biofilms with significantly higher biovolumens with clumped and tall colonies, type I FimA strains form biofilms with a dense basal monolayer and disperse microcolonies (Kuboniwa et al., 2009).
Hemagglutinins are surface proteins associated with bacterial adhesion to the host cells and to erythrocyte agglutination. These proteins facilitate the acquisition of nutrients (heme) by the bacterium through erythrocyte binding (Lépine and Progulske-Fox, 1996). P. gingivalis possesses at least eight hemagglutinins encoded by the hag genes (hagA to hagE) (Progulske-Fox et al., 1989, 1995; Nelson et al., 2003). Among them, the hemagglutinins HagA, HagB, and HagC have been described as important virulence factors (Bélanger et al., 2012; Connolly et al., 2017). HagA is involved in adhesion to both gingival epithelial and endothelial cells (Bélanger et al., 2012); as well as contributing to bacterial survival in both epithelial cells in vitro and murine abscesses (Miller et al., 2017). Conversely, HagB and HagC do not contribute to P. gingivalis virulence in vivo (Miller et al., 2017), although HagB has been associated with bacterial adhesion in vitro (Song et al., 2005).
Virulent and avirulent P. gingivalis strains have been compared in previous studies in order to find genetic determinants that explain the differences in virulence. The comparison between the virulent strain W83 of P. gingivalis with the attenuated strain 33277, using microarrays, showed that their chromosomes were very similar, with ~93% of the predicted genes in common. Intriguingly, 7% of the genes showed very low or no signals in ATCC 33277, suggesting that the ORFs analyzed are absent or divergent in 33277 strain (Chen et al., 2004). One of the most important differences found was the presence of an ORF involved in the synthesis of capsular polysaccharide in W83, which was absent in 33277 (Chen et al., 2004). This is an important virulence factor that contributes to inflammation and participate in the interaction of P. gingivalis with other members of the oral biofilm (Laine and van Winkelhoff, 1998; Davey and Duncan, 2006). Moreover, by using comparative genomics, extensive genomic rearrangements were observed in W83 compared with 33277, including inversions, translocations, deletions or replacements (Naito et al., 2008). More directly associated with pathogenic mechanisms, a genomic comparison between an invasive (W83) and a non-invasive (AJW4) strain of P. gingivalis was performed using whole genome microarray, showing a series of hypothetical ORF, which are polymorphic (Dolgilevich et al., 2011). Moreover, some of them were associated to the invasion of endothelial cells in vitro, explaining –at least in part- why the AJW4 strain is non-invasive.
In this work, we performed whole genome sequencing of a hyper virulent strain (CP3) obtained from a Chilean periodontitis patient exhibiting severe clinical signs of inflammation and tissue damage, as well as of an attenuated strain (H3) obtained from a periodontally-healthy individual. In addition, we performed a comparison between CP3 and H3 with all the available genomes of P. gingivalis through whole-genome nucleotide identity and a maximum likelihood phylogenetic inference to resolve the evolutionary relationships among P. gingivalis strains. Using pan-genome analysis and homology searches, we evaluated the presence and absence of genes associated with virulence factors in order to test for differences in gene content contributing to bacterial virulence. We next evaluated whether such differences correlated with different virulence properties in vitro, by performing gingival epithelial cell adhesion/invasion assays, as well as biofilm and hemagglutination experiments in vitro.
Materials and Methods
Strain Information
CP3 strain was isolated from a subgingival plaque pool of a patient with severe periodontitis, while H3 strain was isolated from a periodontally healthy subject at the Facultad de Odontología, Universidad de Chile. Periodontal status was evaluated by a calibrated clinician considering clinical parameters of periodontitis severity, such as clinical attachment level (CAL) and probing depth (PD) (Table S1). While CP3, showed high virulent properties in vitro even when compared to the highly virulent strain W50; H3 showed low virulent properties in vitro (Díaz et al., 2015; Soto et al., 2016)
Bacterial Strains and Culture Conditions
P. gingivalis reference strains ATCC 33277, W50 (ATCC 53978), and clinical isolates CP3 and H3, were grown as described previously (Soto et al., 2016). Briefly, they were grown in an anaerobic atmosphere at 37°C in blood agar (5% defibrinated sheep blood) or in brain-heart infusion broth (BHI; Oxoid), both supplemented with 5 mg/mL hemin-menadione.
Invasion and Adhesion Assays to Gingival Epithelial Cells
P. gingivalis was grown in BHI supplemented with 5 mg/mL hemin-menadione and incubated at 37°C in anaerobiosis to an optical density of 600 nm (OD600) = 0.6–0.8. For invasion assays, gingival epithelial cells (OKF6/TERT2) were infected at a multiplicity of infection (MOI) of 100 for 90 min at 37°C and 5% CO2. After infection, cells were washed with PBS and incubated 2 h with fresh media supplemented with 300 ug/mL gentamicin and 200 μg/mL metronidazole. Then, cells were rinsed with PBS and incubated with 1% saponin, followed by serial dilution of the resulting supernatant, and plating on blood agar plates supplemented with hemin and menadione. After 5–7 days, colony-forming units (CFU) were quantified. Adhesion assays were performed as described above, except that cells were infected at 4°C. Results were compared by applying ANOVA testing and Dunnett's multiple comparison to analyze statistical differences.
DNA Sequencing and Genome Assembly
Total DNA was extracted from the P. gingivalis clinical isolates CP3 and H3 in tubes containing lysis/stabilization buffer. Samples were lysed using bead-beating, and DNA was extracted by a guanidine thiocyanate silica column-based purification method (Almonacid et al., 2017). CP3 and H3 DNA samples were used to prepare paired-end libraries following the Illumina sample preparation guide and sequenced in an Illumina MiSeq at the Centro de Biotecnología Vegetal, Universidad Andrés Bello (Santiago, Chile). For the CP3 sample, two DNA libraries, “L1” and “L2,” were prepared and sequenced (technical replicate). For quality control, reads were screened and filtered for adapter sequences and reads shorter than 50 bp. Reads with low quality bases (<Q20) were trimmed using FastQC v0.11.8 and PRINSEQ v.0.20.4 (Andrews, 2010; Schmieder and Edwards, 2011). Quality-controlled reads of CP3 libraries L1 and L2, and H3 library were used to perform a de novo assembly of both CP3 and H3 genomes, using SPAdes v3.10.0 (Bankevich et al., 2012). The assembled genomes were annotated using Prokka v1.12 (Seemann, 2014). CP3 and H3 genome assemblies were evaluated by the calculation of statistical values (e.g., N50, number of contigs and bases, and GC content) and completeness analysis through the search of bacterial ortholog genes (OrthoDB database), using NGSQCToolkit v2.3.3 and BUSCO v3 (Patel and Jain, 2012; Waterhouse et al., 2017), respectively. Also, quality-controlled reads were aligned against the assembled contigs to calculate CP3 and H3 genomes coverage, using Bowtie v 2.3.4.1 and Samtools v1.7 (Li et al., 2009; Langmead and Salzberg, 2012), respectively.
P. gingivalis Genomic Dataset
In order to perform a comparative genomic analysis of CP3 and H3 with others P. gingivalis strains, a total of 62 P. gingivalis genomes were retrieved from NCBI, RefSeq database (as of February 2018). All downloaded P. gingivalis genomes were reannotated using Prokka as with CP3 and H3 genomes.
Whole-Genome Nucleotide Identity and Phylogenetic Relationships
Average nucleotide identity (ANI) was calculated for all 64-genome P. gingivalis dataset (i.e., 62 genomes from RefSeq plus CP3 and H3 genomes), using the pyani Python3 module (Pritchard et al., 2016) and pheatmap R packages (https://cran.r-project.org/web/packages/pheatmap/index.html) for results' visualization. kSNP3.0 was used for SNP identification in the 64 P. gingivalis genomes and the core polymorphic positions were used for the phylogenetic reconstruction through Maximum Likelihood inference using RAxML v8.2.9. The obtained tree was plotted using FigTree v1.4.4 (Stamatakis, 2014; Gardner et al., 2015) (http://tree.bio.ed.ac.uk/software/figtree/).
Pan-Genome Analysis and Search for Virulence Genes
The pan-genome analysis of the 64-genome P. gingivalis was performed through the Roary v3.7.0 pipeline (Page et al., 2015). For virulence gene identification, a homology search of predicted protein sequences in CP3 and H3 genomes against the most characterized virulence genes identified in P. gingivalis available in the literature was performed, using CRB-Blast v0.6.6 (https://github.com/cboursnell/crb-blast). Also, virulence genes were identified in the 64-genome dataset using pan-genome results to evaluate presence/absence and protein sequence variation of these genes. According to shared and unique genes between CP3 and H3 genomes, a Venn diagram was plotted using the VennDiagram R package (https://cran.r-project.org/package=VennDiagram). The inferred phylogenetic tree was annotated based on the presence/absence of the virulence genes that are part of the accessory genome, using the pheatmap R package and FigTree.
Biofilm Formation Assays
P. gingivalis cultures grown to exponential phase were obtained, diluted, and inoculated to each well of a 96-well flat bottom plate in supplemented BHI medium (and 1% tryptic soy broth, TSB). Then, plates were incubated anaerobically at 37°C for 72 h. After the incubation time, the supernatant was removed, wells were washed twice by immersion in distilled water and the plate was left to air dry for 1 h. Biofilm staining was performed by adding 100 μL of 0.1% safranin and incubating for 15 min, followed by two washes by immersion in distilled water. Finally, the dye retained in the biofilm was eluted incubating with 100 μL of 95% ethanol for 5 min, the elution was transferred to a new 96-well plate and the absorbance was measured at 490 nm. Results were compared by applying ANOVA testing and Dunnett's multiple comparison to analyze statistical differences.
Hemagglutination Assays
One mL of defibrinated horse blood was centrifuged at 2,100 rpm for 5 min. The resulting pellet (red blood cells) was washed 3 times and then diluted in PBS to 2% solution. In parallel, P. gingivalis strains (W50, H3, and CP3) were grown to an exponential phase culture, which was adjusted to DO600 = 2.0 (3 × 108 CFU/mL). Then, two hundred μL of each suspension was added to one well of a 96-well round-bottom plate. After that, each suspension (W50, H3, and CP3) was serially diluted, by taking 100 μL and mixed with 100 μL PBS (1:2 dilution). This step was repeated until 1:64 dilution. Finally, each well was mixed with an equal volume of 2% sheep erythrocytes at 37°C for 3 h. Experiments were repeated three times.
hagC Expression Analysis by Quantitative Reverse Transcription Polymerase Chain Reaction
Total cytoplasmic RNA was isolated from CP3, W50, and H3 strains using the TRIZOL method as described previously (Chomczynski and Mackey, 1995). Reverse transcription was performed using M-MLV reverse transcriptase (Promega, Madison, WI, USA). To quantify the mRNA expression for hagC and 16S rRNA genes, 1,000 ng of RNA were used to synthetize cDNA and amplify by quantitative real-time polymerase chain reaction (qRT-PCR) using appropriate primers (HagC: TTTGCCAAGAATGTGCTGAC and GTCGAGGGCTATGACCTGAG; 16S: AGGCAGCTTGCCATACTGCG and ACTGTTAGCAACTACCGATGT) and the Brilliant II SYBR qPCR master mix (Agilent, Santa Clara, CA, USA) in a Mx3000p qPCR system (Agilent, Santa Clara, CA, USA). The cycling program was as follows: 95°C for 10 min, followed by 45 cycles of 95°C for 10 s, 58°C for 30 s, and 72°C for 20 s. Finally, a melting curve was obtained by incubating at 95°C for 15 s, 60°C for 1 min, and 95°C for 15 s, to detect non-specific product formation and false positive amplification. As an endogenous control, we determined the16S rRNA gene expression levels. Data were statistically analyzed using PRISM software (version 6.0; GraphPad, La Jolla, CA, USA). Comparisons between strains were made using by applying ANOVA testing and Dunnett's multiple comparison to analyze statistical differences.
Data Availability
This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession SGBA00000000 (CP3) and SGAZ00000000 (H3), under BioProject PRJNA521311. The version described in this paper is version SGBA01000000 (CP3) and SGAZ01000000 (H3).
Results
P. gingivalis CP3 and H3 Genomes
A total of 3.43 million sequences were sequenced and an average of 8% of the raw sequence reads were filtered out. Quality-controlled reads of the CP3 L1 and L2 libraries, and H3 library were used to perform the de novo assembly of the 2.25 and 2.31 Mb-length genomes of CP3 and H3, respectively. CP3 and H3 genomes are composed by 118 and 165 contigs, with a N50 value of 37,600 and 39,176 bp, GC content of 48.47 and 48.35%, and predicted protein count of 1,896 and 1,955, respectively. These values are similar with average statistics values for the 62 P. gingivalis available genomes at NCBI's RefSeq database (length 2.33 Mb, protein count 1,896, and GC content 48.4%; as of February 2018). Both assembled genomes exhibit high completeness, having a 95.3% of bacterial orthologs present in the OrthoDB database. The average coverage of CP3 and H3 genomes is 300x and 28x, respectively. Genome features of both strains are displayed in Table S2.
Whole-Genome Nucleotide Identity and Phylogenetic Relationships of P. gingivalis Strains
In order to compare the whole-genome nucleotide identity of CP3 and H3 genomes with other previously sequenced P. gingivalis strains, an ANI analysis using the 64-genome P. gingivalis dataset (i.e., 62 genomes from RefSeq plus CP3 and H3 genomes) was performed (Table S3). The analysis shows high nucleotide identity among them (>98.22%; alignment fraction >80%; Figure 1). The strain WW2842 shows the highest nucleotide identity value with the CP3 genome (99.13%), while the H3 genome shows the highest identity with WW2881 and KCOM 2799 (99.69 and 99.20%, respectively; Table S4). There are no publications reporting virulence studies for WW2842, WW2881 or KCOM 2799 strains. However, WW2842 and WW2881 were isolated from subgingival plaque of patients with advanced periodontal disease in London, United Kingdom (BioProject: PRJNA401301, BioSample: SAMN07602651, BioSample: SAMN07602653). KCOM 2799 also was isolated from the subgingival plaque of a patient with periodontal disease in Gwangju, Korea (BioProject: PRJNA415884, BioSample: SAMN07836928). Moreover, the phylogenetic tree shows that CP3 and WW2842, and H3 and WW2881 and KCOM 2799 are grouped together, while CP3 and H3 (showing a 98.57% of identity), are found in different clades (Figure 2). The phylogeny also shows that CP3 is closer to the virulent strain MP4-504 (To et al., 2016), compared with other virulent strains such as W50, W83, ATCC49417, and A7436, which are more distantly related (Figure 2; Ebersole et al., 1995; Chen et al., 2004). On the other hand, H3 is phylogenetically related to AJW4, a previously reported very low invasive strain (Dolgilevich et al., 2011). However, H3 is found at a greater distance from other two type strains, ATCC 33277 and 381, which are comparatively less virulent than W50, W83 (Ebersole et al., 1995; Chen et al., 2004). These results suggest a heterogeneous nature of P. gingivalis species, not necessarily related to the periodontal status of their reservoirs.
Figure 1. Average nucleotide identity (ANI) in the 64-genome Porphyromonas gingivalis dataset. Hierarchical clustering of the P. gingivalis genomes based on their average nucleotide identity values. The color gradient from blue to red shows the percentage of identity, from lowest to highest, that each pair of genomes shares.
Figure 2. Phylogenetic relationships among Porphyromonas gingivalis strains and presence/absence of virulence genes. (Left) Maximum likelihood tree inferred from a core SNPs alignment (16,645 polymorphic positions) of the 64-genome P. gingivalis dataset. Node labels represent the bootstrap support values. (Right) Heatmap showing the presence (black) or absence (white) of virulence genes in each P. gingivalis genome (only virulence genes belonging to the accessory genome are shown).
From the virulence genes characterized for P. gingivalis in the literature (Table 2), 12 are part of the 64-genome dataset accessory genome and were selected to evaluate their distribution along the phylogeny (Figure 2). Results show that identical patterns of virulence factor presence and absence belong to phylogenetically unrelated and often geographically separated P. gingivalis strains that exhibit similar virulence phenotype. For instance, P. gingivalis 15_9 belongs to a clade with other virulent strains (W83; W50; A7436) and exhibits exactly the same pattern of virulence genes than CP3, though they are not closely related (Figure 2). Likewise, SJD12 strain exhibits the same pattern of virulence genes than H3 while being far apart in the phylogeny. Altogether, these results suggest that similar patterns of presence/absence of virulence genes may arise by convergent evolution and yield similar phenotypes.
P. gingivalis Pan-Genome and Virulence Genes
Since the P. gingivalis strains included in these analyses come from different individuals and many of them are associated with different levels of virulence, a pan-genome analysis of the 64-genome P. gingivalis dataset was performed. The pan-genome is composed of 7,886 genes, of which 1,229 belong to the core genome and 6,657 genes belong to the accessory genome. Specifically, 65.44% of the CP3 genes and 64.08% of the H3 genes belong to the core genome, while 34.56% of the CP3 genes and 35.92% of the H3 genes belong to the accessory genome (Table 1).
Table 1. Pan-genome results of the 64-genome Porphyromonas gingivalis dataset.
Next, we asked whether genes previously associated with virulence in P. gingivalis (Table 2) were either present or absent in both CP3 and H3 genomes. To this end, an exhaustive search of the virulence genes reported for P. gingivalis ATCC 33277 (reference strain) in the literature, and a screening of these genes in CP3 and H3 genomes through a bidirectional blast was conducted (Table 2). Interestingly, most of the virulence genes reported in the genome of the ATCC 33277 strain are also present in both CP3 and H3 genomes, including those encoding virulence factors as important as gingipains and capsule synthesis proteins (Table 2; Figure 2). Also, the pan-genome analysis was used to evaluate the presence/absence of virulence genes in CP3 and H3 genomes. Intriguingly, although both CP3 and H3 possess hemagglutinin genes, hagA and hagC are differentially found in their genomes. While CP3 have six ORFs coding for hemagglutinins (CP3_01418, CP3_01421, CP3_01422, CP3_01423, CP3_01568, CP3_01901), H3 have only four (H3_00216, H3_01383, H3_01678, H3_01679). Three of these genes belong to the CP3 and H3 core genome (>95% identity) and encode non-characterized hemagglutinins (CP3_01422, CP3_01423, CP3_01568; H3_00216, H3_01678, H3_01679). The other three hag copies found in CP3 genome and one of the H3 gene copies belong to their accessory genome (<95% identity). Moreover, there are differences in the type of hemagglutinins found in the accessory genome: while hagA is only present in CP3 (CP3_01901) and no in H3, CP3 possesses two copies of hagC (CP3_01418, CP3_01421), while H3 possesses only one (H3_01383) (Figure 3).
Table 2. Genes involved in Porphyromonas gingivalis virulence.
Figure 3. Common genes between Porphyromonas gingivalis CP3 and H3 and differences in the copy-number of virulence genes. (A) Venn diagram showing the number of common and unique genes of CP3 (blue) and H3 (green) genomes. Values in parentheses are the percentage of genes of the total genes in CP3 and H3 pan genome. (B) Bar plot showing six virulence genes that we detected to have differences in their copy-number between CP3 (blue bars) and H3 (green bars) genomes.
Another interesting virulence-associated gene is fimA, which encodes different major fimbriae types. Interestingly, while type IV fimA is present in CP3 and absent in H3, type I fimA is present in H3 and absent in CP3 (Figures 2, 3). Finally, differences in two other genes were detected between both strains; while epsD gene encoding for a putative glycosyltransferase was only present in CP3, epsJ was found in H3 and not in CP3 (Figures 2, 3).
Hemagglutination Ability of CP3 and H3 Strains
As described above, the CP3 strain possesses one copy of hagA and two of hagC, while H3 possesses only one copy of hagC (Figure 3). In order to evaluate whether there were differences in the hemmaglutination ability between CP3 and H3 that could be associated with these genetic differences, an in vitro hemagglutination assay was performed. Our results show that while CP3 was able to agglutinate red cells until dilution 1:2 (1.5 × 108 CFUs), H3 only agglutinates when 3.0 × 108 CFUs where added. The hemagglutination titer of red blood cells showed by the virulent reference strain W50 was similar compared to CP3 and higher compared to H3 (Figure 4). Therefore, although H3 hemagglutinate red cells, it requires a higher bacterial load compared to CP3 and W50. In order to evaluate if the differences in the copy number of the hagC gene could influence the hemagglutination ability of H3 and CP3, we evaluated its expression by RT-qPCR. Accordingly, the expression of hagC was significantly lower in H3 strain compared to CP3, which suggests that the copy number might explain the observed differences in hemagglutination (Figure S1). In addition, we evaluated the hagC transcript levels of the reference strain W50, determining a higher expression compared to H3.
Figure 4. Hemagglutination activity of Porphyromonas gingivalis strains CP3, H3, and W50. Hemagglutinin activities of P. gingivalis (W50, H3, and CP3) were assessed in cells serially diluted (from left to right) and incubated with sheep erythrocytes (2%) in a round-bottom microtiter plate. Dilution folds were listed on the top. Red arrows indicate the last dilution that showed full agglutination. The figure is representative of three biological replicates.
Adherence and Biofilm Formation Abilities of CP3 and H3 Strains
Since CP3 possesses fimbria type IV (widely associated with epithelial cell invasion) and H3 possesses fimbria type I (found in the periodontal healthy population; Amano et al., 2000), an assay was performed to evaluate if such genetic differences are reflected in their adhesion and invasion ability to oral epithelial cells. Although no differences in attachment level to epithelial cells were observed between H3, CP3 and the W50 reference strain (Figure 5A), CP3 invasion ability was ~6-fold higher than both H3 and the reference strain (Figure 5B).
Figure 5. Cellular adherence (A), invasion (B), and biofilm formation (C) of different strains of Porphyromonas gingivalis exhibiting different fimA genotypes. Adherence (A) or invasion (B) were calculated as: the number of viable CFU in infected cells/the number of CFU of initial bacteria inoculum. Biofilm formation (C) was determined by safranin staining. The dye retained in the biofilm was eluted and measured at 490 nm. Invasion, adherence and biofilm assays were performed in triplicate on at least three independent occasions. Values that differ significantly from controls are indicated (SEM; n = 3; *p ≤ 0.05).
Since there are differences in the fimA genes encoded in CP3 and H3 genomes, the next step was to evaluate their biofilm formation ability. Results show that CP3 (encoding type IV fimA) possesses a significantly higher ability to form biofilms that H3 (encoding type I fimA) (Figure 5C). Interestingly, both strains formed more biofilm than both reference strains ATCC 33277 (type I) and W50 (type IV).
Discussion
P. gingivalis has been largely associated with periodontitis, showing both one of the highest prevalence and abundance among periodontitis-associated species (Griffen et al., 2012; Abusleme et al., 2013; Kistler et al., 2013; Hong et al., 2015; Kirst et al., 2015; Park et al., 2015). Nevertheless, using molecular approaches, P. gingivalis have also been detected in periodontally healthy individuals, although in lower abundance (Abusleme et al., 2013; Díaz et al., 2015; Choi et al., 2018).
We previously studied two P. gingivalis clinical isolates showing differences in their virulence, one from a chronic severe periodontitis patient (CP3), and the other from a periodontally healthy subject (H3) (Díaz et al., 2015; Soto et al., 2016). Interestingly, virulent properties of these strains are in agreement with the clinical signs linked to periodontal damage determined in the donor subjects. Periodontal parameters observed in the subject from whom the CP3 isolate was obtained had an average probing depth (PD) = 4.2 mm and a clinical attachment loss (CAL) = 6.1 mm. Moreover, from the 107 sites analyzed in the oral exam, 47 sites were ≥ 5 mm and 13 sites ≥ 7 mm, indicating a high level of tissue damage (Table S1). In contrast, the periodontally healthy subject (from whom H3 was obtained) had a PD = 2.5 mm and a CAL = 3.6 mm. No sites ≥ 5 mm were observed from 167 sites analyzed, 166 were ≤ 3 mm (Table S1).
In previous works, we showed that CP3 –compared with H3- showed both higher resistance to antimicrobial peptides (Díaz et al., 2015) and ability of decreasing apoptosis of gingival epithelial cells (Soto et al., 2016). A phenotype that was associated with the presence of high molecular O-antigen molecules in CP3 (Díaz et al., 2015). Therefore, in this work we investigated the existence of other genotypic features that could explain differences in virulence that help us to understand the contribution of P. gingivalis to periodontitis.
Our results showed that CP3 and H3 genomes share ~98 to ~99% identity with other available P. gingivalis genomes with an alignment coverage from ~83 to 95%; this translates in a ~64 to 65% of common genes (core genome; >95% identity) and ~34 to 35% of “unique” genes (accessory genome; <95% identity; Figure 1). Although both strains locate in relatively distal branches in the phylogenetic tree, they are not necessarily closer to either more virulent strains (in the case of CP3), or to less virulent ones (in the case of H3). However, relevant genes involved in virulence show different patterns of presence/absence among P. gingivalis strains (Figure 2).
Notably, we were able to correlate genomic differences with virulence features. First, CP3, harboring one copy of the hemagglutinin gene hagA and two copies of hagC, shows higher hemmaglutination ability in vitro compared to H3 (Figures 3, 4), which encodes only one copy of hagC gene (Figure 3). HagA has been widely related erythrocytes adhesion for nutrients acquisition (Lépine and Progulske-Fox, 1996; Lamont and Jenkinson, 2000). Therefore, we speculate that its absence in H3 impairs its hemagglutination ability. Variation of comparable magnitude in the hemagglutination ability of P. gingivalis strains has been previously reported in strains with other genetic backgrounds (Shoji et al., 2013; Zhang et al., 2016; Puth et al., 2019). Interestingly, among P. gingivalis hemagglutinins, hagA, hagD, and hagE show >70% identity, while hagB and hagC correspond to unrelated families, sharing 93% identity (Connolly et al., 2017). Moreover, although HagB and HagC functions have only been partially defined, it has been reported that both of them produce hemagglutination, although they are not the sole responsible for such activity (Lépine and Progulske-Fox, 1996). Interestingly, hemagglutination in P. gingivalis has been associated to other proteins, such as specific domains of gingipains RgpA and Kgp (Li et al., 2011). Also, it was recently reported that pckA contributes to hemagglutination, independently of the gingipains activity (Wu et al., 2018). Thus, in spite of an association between the genomic differences in hagA and hagC, as well as the levels of expression of hagC (Figure S1) with the hemagglutination activity found in this study, the presence of other(s) virulence factor(s) of P. gingivalis that also contribute to this phenotype cannot be discarded.
To our knowledge, no evidence of association between virulence and copy numbers of specific genes has been described in P. gingivalis. However, in another bacterial model, no association between the virulence of Vibrio parahaemolyticus and the number of copies of Vp_PirA- and B-like toxin genes in their genome was found (Pérez-Chaparro et al., 2008). In this context, the question of whether the sole copy of hagC found in H3 is either not functional or not sufficient to produce hemagglutination remains to be elucidated. Nevertheless, our results show that a H3 exhibit a lower expression of hagC compared to CP3, which is associated with a lower hemagglutination capacity of H3 strain, suggesting a functional role for this difference.
Similarly, others have also reported genetic differences that rely not only in presence/absence of specific genes. Several studies have reported extensive genetic heterogeneity among P. gingivalis strains obtained from periodontitis patients, using molecular typing methods such as restriction endonuclease analysis (REA) (Genco and Loos, 1991), restriction fragment length polymorphism (RFLP) (Loos and Dyer, 1992), multilocus enzyme electrophoresis (MEE) (Loos et al., 1993), and arbitrarily primed polymerase chain reaction (AP-PCR) (Ménard et al., 1992). Furthermore, variations in the DNA of P. gingivalis analyzed by heteroduplex assays [that permit to identify polymorphisms in the ribosomal internal spacer region (ISR)], revealed the presence of heteroduplex types; some of them significantly associated with periodontitis and others equally detected in health or disease (Griffen et al., 1999).
Regarding FimA genotypes, type I fimA has been found mainly in the periodontally healthy population (Feng et al., 2014), as well as in the avirulent type-strain ATCC 33277. In turn, type IV fimA is more frequently associated with P. gingivalis obtained from periodontitis patients (Amano et al., 2004) and is also present in the virulent reference strain W83 (Figure 2). Moreover, it has been shown that phenotypes based on the presence/absence of the six variants of nucleotide sequences of fimA gene (types I to VI) induce different cellular responses (Missailidis et al., 2004). Particularly, type I FimA have been associated with biofilms with a dense basal monolayer and disperse microcolonies; whereas type IV FimA, with thicker biofilms (Kuboniwa et al., 2009). This is consistent with our results: H3 lacks gene CP3_00160 coding for the major fimbrium subunit FimA type IV, but exhibits a major fimbrium subunit FimA type I (Figure 3). In contrast, CP3 lacks the gene coding for major fimbrium subunit FimA type I and possesses the gene coding for the major fimbrium subunit FimA type IV. This genotypic difference is associated with the lower biofilm formation ability of H3, as well as with less intracellular invasion to oral epithelial cells compared to the CP3 strain (Figure 5).
Interestingly, among P. gingivalis virulence factors, not only FimA has been related to epithelial cell adhesion/invasion and biofilm formation. HagA promotes adhesion to epithelial cells, and antibodies against its hemagglutinin domain protects against oral infection in a murine model (Frazer et al., 2006; Bélanger et al., 2012). Therefore, it cannot be discarded that the differences in the hemagglutinins encoded in CP3 genome, and no in H3, are also contributing to CP3 adhesion/invasion abilities in our in vitro model. The hagB gene of P. gingivalis, which is 98.6% identical to hagC, encodes HagB that mediates adhesion to host cells (Progulske-Fox et al., 1999; Song et al., 2005). Hence, it is likely that one could “substitute” the other's function. Therefore, as discussed in the context of differences in hemagglutination, it is unlikely that the lack in only one copy of hagC found in H3 has a major effect on its ability of adhere to and invade epithelial cells. Similarly, HagA has been associated with co-aggregation between P. gingivalis and Treponema denticola, which is another periodontitis-associated bacteria (Ito et al., 2010); and HagC mediates multi-species biofilm formation (Connolly et al., 2017). Thus, besides FimA, CP3 enhanced ability to form biofilms could also be attributed to hemagglutinin differences compared to H3.
Unlike previous reports comparing clinical isolates obtained from patients showing different periodontal disease severity, no differences in capsular antigen genotypes were found between CP3 and H3 strains (van Winkelhoff et al., 1993; Laine et al., 1996; Laine and van Winkelhoff, 1998). Such differences have also been found between the reference strains W83 and ATCC 33277; while the gene coding for the synthesis of capsular polysaccharide is present in W83, it is absent in the ATCC 33277 (Chen et al., 2004). Nonetheless, while the glycosyltransferase EpsD was only present in CP3, the glycosyltransferase EpsJ was only found in H3 (Figures 2, 3). Although, no virulence function has been attributed to such proteins in P. gingivalis, they have been related to exopolysaccharide synthesis in other bacterial species. In Bacillus subtilis, the locus epsHIJK is necessary for the synthesis poly-N-acetylglucosamine, which is a key component of the biofilm (Roux et al., 2015). epsJ genes are annotated as putative glycosyltransferases, showing similarity to icaA in Staphylococcus aureus and pgaC in Escherichia coli. EpsD and EpsJ have also been associated to exopolysaccharides synthesis in Lactococcus lactis (van Kranenburg et al., 1999). While EpsD links glucose-1-phosphate from UDP-glucose to a lipid carrier, EpsJ is likely to be involved in releasing the trisaccharide backbone from the lipid carrier. Due to the relevance of exopolysaccharides in biofilm formation and taking in to account the relevance of such structures in P. gingivalis virulence, we cannot discard that these genetic differences between CP3 and H3 could participate in their biofilm formation ability (Figure 4). Altogether, through the analysis of a dataset composed by 64 P. gingivalis genomes, this study provides relevant results regarding genetic determinants and their association with P. gingivalis virulence.
Data Availability
The datasets generated for this study can be found in NCBI PRJNA521311.
Author Contributions
KM performed the comparative genomics analyses, wrote, and reviewed the manuscript. AH analyzed the data, wrote and reviewed the manuscript. CS carried out hemagglutination, adhesion/invasion assays, and RT-qPCR experiments. IB carried out biofilm assays. MO carried out hemagglutination assays. CM performed sequencing of P. gingivalis strains. JP-D contributed to conception and critically reviewed the manuscript. EC-N contributed to the design of the genomic analyses, wrote, and reviewed the manuscript. DB conceived and designed the research, wrote, and reviewed the manuscript. All authors contributed to the final version.
Funding
This work was supported by grants FIOUCh n° 17/020 and CONICYT-FONDAP 15130011 (DB). CONICYT-FONDECYT n° 1151255(JP-D). CONICYT-FONDECYT de iniciación n° 11160905 (EC-N).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
EC-N would like to thank The George Washington University's high-performance computing facility, ColonialOne, for providing data storage, support, and computing power for genomic analyses (colonialone.gwu.edu).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2019.00246/full#supplementary-material
Table S1. Periodontal parameters of the subjects from whom Porphyromonas gingivalis strains CP3 and H3 were obtained. PD, pocket depth; CAL, clinical attachment level.
Table S2. Sequencing, quality control (QC) of reads, and assembly of the Porphyromonas gingivalis strains CP3 and H3 genomes. The first panel “Sequencing and quality control (QC)” shows the number of sequenced reads and the number of reads after QC. The quality-controlled reads were used for de novo assembly of CP3 and H3 genomes. Statistical features of both assembled genomes are shown in the second panel “Genome assembly,” including the number of predicted proteins and the presence of bacterial orthologs (“Proteome completeness”). C, complete orthologs; S, complete and single-copy orthologs; D, complete and duplicated orthologs; F, fragmented orthologs; M, missing orthologs.
Table S3. Publicly available Porphyromonas gingivalis genomes used in this study. Sixty-two P. gingivalis genomes were downloaded from NCBI's RefSeq database in February 2018.
Table S4. Average nucleotide identity and coverage values for the 64-genome Porphyromonas gingivalis dataset. Sheets 1, 2 show the average nucleotide identity and coverage values, respectively. Sheets 3, 4 detail the identity and coverage values, respectively, of all P. gingivalis strain genomes against CP3 and H3 genomes.
References
Abusleme, L., Dupuy, A. K., Dutzan, N., Silva, N., Burleson, J. A., Strausbaugh, L. D., et al. (2013). The subgingival microbiome in health and periodontitis and its relationship with community biomass and inflammation. ISME J. 7, 1016–1025. doi: 10.1038/ismej.2012.174
Aduse-Opoku, J., Slaney, J. M., Hashim, A., Gallagher, A., Gallagher, R. P., Rangarajan, M., et al. (2006). Identification and characterization of the capsular polysaccharide (K-antigen) locus of Porphyromonas gingivalis. Infect. Immun. 74, 449–460. doi: 10.1128/IAI.74.1.449-460.2006
Almonacid, D. E., Kraal, L., Ossandon, F. J., Budovskaya, Y. V., Cardenas, J. P., Bik, E. M., et al. (2017). 16S rRNA gene sequencing and healthy reference ranges for 28 clinically relevant microbial taxa from the human gut microbiome. PLoS ONE 12:e0176555. doi: 10.1371/journal.pone.0176555
Amano, A., Kuboniwa, M., Nakagawa, I., Akiyama, S., Morisaki, I., and Hamada, S. (2000). Prevalence of specific genotypes of Porphyromonas gingivalis fimA and periodontal health status. J. Dent. Res. 79, 1664–1668. doi: 10.1177/00220345000790090501
Amano, A., Nakagawa, I., Okahashi, N., and Hamada, N. (2004). Variations of Porphyromonas gingivalis fimbriae in relation to microbial pathogenesis. J. Periodontal. Res. 39, 136–142. doi: 10.1111/j.1600-0765.2004.00719.x
Andrews, S. (2010). astQC a Quality Control Tool for High Through- put Sequence Data. Available online at: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/
Babic, A., Poole, E. M., Terry, K. L., Cramer, D. W., Teles, R. P., and Tworoger, S. S. (2015). Periodontal bone loss and risk of epithelial ovarian cancer. Cancer Causes Control CCC 26, 941–947. doi: 10.1007/s10552-015-0575-7
Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021
Bélanger, M., Kozarov, E., Song, H., Whitlock, J., and Progulske-Fox, A. (2012). Both the unique and repeat regions of the Porphyromonas gingivalis hemagglutin A are involved in adhesion and invasion of host cells. Anaerobe 18, 128–134. doi: 10.1016/j.anaerobe.2011.10.005
Chaudhuri, S., Pratap, S., Paromov, V., Li, Z., Mantri, C. K., and Xie, H. (2014). Identification of a diguanylate cyclase and its role in Porphyromonas gingivalis virulence. Infect. Immun. 82, 2728–2735. doi: 10.1128/IAI.00084-14
Chen, T., Hosogi, Y., Nishikawa, K., Abbey, K., Fleischmann, R. D., Walling, J., et al. (2004). Comparative whole-genome analysis of virulent and avirulent strains of Porphyromonas gingivalis. J. Bacteriol. 186, 5473–5479. doi: 10.1128/JB.186.16.5473-5479.2004
Choi, H., Kim, E., Kang, J., Kim, H. J., Lee, J. Y., Choi, J., et al. (2018). Real-time PCR quantification of 9 periodontal pathogens in saliva samples from periodontally healthy Korean young adults. J. Periodontal. Implant Sci. 48, 261–271. doi: 10.5051/jpis.2018.48.4.261
Chomczynski, P., and Mackey, K. (1995). Short technical reports. Modification of the TRI reagent procedure for isolation of RNA from polysaccharide- and proteoglycan-rich sources. BioTechniques 19, 942–945.
Ciuraszkiewicz, J., Smiga, M., Mackiewicz, P., Gmiterek, A., Bielecki, M., Olczak, M., et al. (2014). Fur homolog regulates Porphyromonas gingivalis virulence under low-iron/heme conditions through a complex regulatory network. Mol. Oral Microbiol. 29, 333–353. doi: 10.1111/omi.12077
Coats, S. R., To, T. T., Jain, S., Braham, P. H., and Darveau, R. P. (2009). Porphyromonas gingivalis resistance to polymyxin B is determined by the lipid A 4'-phosphatase, PGN_0524. Int. J. Oral Sci. 1, 126–135. doi: 10.4248/IJOS.09062
Connolly, E., Millhouse, E., Doyle, R., Culshaw, S., Ramage, G., and Moran, G. P. (2017). The Porphyromonas gingivalis hemagglutinins HagB and HagC are major mediators of adhesion and biofilm formation. Mol. Oral Microbiol. 32, 35–47. doi: 10.1111/omi.12151
Cugini, C., Klepac-Ceraj, V., Rackaityte, E., Riggs, J. E., and Davey, M. E. (2013). Porphyromonas gingivalis: keeping the pathos out of the biont. J. Oral Microbiol. 5:19804. doi: 10.3402/jom.v5i0.19804
Davey, M. E., and Duncan, M. J. (2006). Enhanced biofilm formation and loss of capsule synthesis: deletion of a putative glycosyltransferase in Porphyromonas gingivalis. J. Bacteriol. 188, 5510–5523. doi: 10.1128/JB.01685-05
Díaz, L., Hoare, A., Soto, C., Bugueño, I., Silva, N., Dutzan, N., et al. (2015). Changes in lipopolysaccharide profile of Porphyromonas gingivalis clinical isolates correlate with changes in colony morphology and polymyxin B resistance. Anaerobe 33, 25–32. doi: 10.1016/j.anaerobe.2015.01.009
Dolgilevich, S., Rafferty, B., Luchinskaya, D., and Kozarov, E. (2011). Genomic comparison of invasive and rare non-invasive strains reveals Porphyromonas gingivalis genetic polymorphisms. J. Oral Microbiol. 3:5764. doi: 10.3402/jom.v3i0.5764
Ebersole, J. L., Kesavalu, L., Schneider, S. L., Machen, R. L., and Holt, S. C. (1995). Comparative virulence of periodontopathogens in a mouse abscess model. Oral Dis. 1, 115–128. doi: 10.1111/j.1601-0825.1995.tb00174.x
Feng, X., Zhang, L., Xu, L., Meng, H., Lu, R., Chen, Z., et al. (2014). Detection of eight periodontal microorganisms and distribution of Porphyromonas gingivalis fimA genotypes in Chinese patients with aggressive periodontitis. J. Periodontol. 85, 150–159. doi: 10.1902/jop.2013.120677
Frazer, L. T., O'Brien-Simpson, N. M., Slakeski, N., Walsh, K. A., Veith, P. D., Chen, C. G., et al. (2006). Vaccination with recombinant adhesins from the RgpA-Kgp proteinase-adhesin complex protects against Porphyromonas gingivalis infection. Vaccine 24, 6542–6554. doi: 10.1016/j.vaccine.2006.06.013
Gardner, S. N., Slezak, T., and Hall, B. G. (2015). kSNP3. 0: SNP detection and phylogenetic analysis of genomes without genome alignment or reference genome. Bioinformatics 31, 2877–2878. doi: 10.1093/bioinformatics/btv271
Genco, R. J., and Loos, B. G. (1991). The use of genomic DNA fingerprinting in studies of the epidemiology of bacteria in periodontitis. J. Clin. Periodontol. 18, 396–405. doi: 10.1111/j.1600-051X.1991.tb02307.x
Griffen, A. L., Beall, C. J., Campbell, J. H., Firestone, N. D., Kumar, P. S., Yang, Z. K., et al. (2012). Distinct and complex bacterial profiles in human periodontitis and health revealed by 16S pyrosequencing. ISME J. 6, 1176–1185. doi: 10.1038/ismej.2011.191
Griffen, A. L., Lyons, S. R., Becker, M. R., Moeschberger, M. L., and Leys, E. J. (1999). Porphyromonas gingivalis strain variability and periodontitis. J. Clin. Microbiol. 37, 4028–4033.
Hajishengallis, G. (2015). Periodontitis: from microbial immune subversion to systemic inflammation. Nat. Rev. Immunol. 15, 30–44. doi: 10.1038/nri3785
Hajishengallis, G., and Lamont, R. J. (2012). Beyond the red complex and into more complexity: the polymicrobial synergy and dysbiosis (PSD) model of periodontal disease etiology. Mol. Oral Microbiol. 27, 409–419. doi: 10.1111/j.2041-1014.2012.00663.x
Hajishengallis, G., Liang, S., Payne, M. A., Hashim, A., Jotwani, R., Eskan, M. A., et al. (2011). Low-abundance biofilm species orchestrates inflammatory periodontal disease through the commensal microbiota and complement. Cell Host Microbe 10, 497–506. doi: 10.1016/j.chom.2011.10.006
Hamada, S., Fujiwara, T., Morishima, S., Takahashi, I., Nakagawa, I., Kimura, S., et al. (1994). Molecular and immunological characterization of the fimbriae of Porphyromonas gingivalis. Microbiol. Immunol. 38, 921–930. doi: 10.1111/j.1348-0421.1994.tb02148.x
Hasegawa, Y., Tribble, G. D., Baker, H. V., Mans, J. J., Handfield, M., and Lamont, R. J. (2008). Role of Porphyromonas gingivalis SerB in gingival epithelial cell cytoskeletal remodeling and cytokine production. Infect. Immun. 76, 2420–2427. doi: 10.1128/IAI.00156-08
Hiratsuka, K., Kiyama-Kishikawa, M., and Abiko, Y. (2010). Hemin-binding protein 35 (HBP35) plays an important role in bacteria-mammalian cells interactions in Porphyromonas gingivalis. Microb. Pathog. 48, 116–123. doi: 10.1016/j.micpath.2010.01.001
Hong, B. Y., Furtado Araujo, M. V., Strausbaugh, L. D., Terzi, E., Ioannidou, E., and Diaz, P. I. (2015). Microbiome profiles in periodontitis in relation to host and disease characteristics. PLoS ONE 10:e0127077. doi: 10.1371/journal.pone.0127077
Ito, R., Ishihara, K., Shoji, M., Nakayama, K., and Okuda, K. (2010). Hemagglutinin/Adhesin domains of Porphyromonas gingivalis play key roles in coaggregation with Treponema denticola. FEMS Immunol. Med. Microbiol. 60, 251–260. doi: 10.1111/j.1574-695X.2010.00737.x
Kirst, M. E., Li, E. C., Alfant, B., Chi, Y. Y., Walker, C., Magnusson, I., et al. (2015). Dysbiosis and alterations in predicted functions of the subgingival microbiome in chronic periodontitis. Appl. Environ. Microbiol. 81, 783–793. doi: 10.1128/AEM.02712-14
Kistler, J. O., Booth, V., Bradshaw, D. J., and Wade, W. G. (2013). Bacterial community development in experimental gingivitis. PLoS ONE 8:e71227. doi: 10.1371/journal.pone.0071227
Kuboniwa, M., Amano, A., Inaba, H., Hashino, E., and Shizukuishi, S. (2009). Homotypic biofilm structure of Porphyromonas gingivalis is affected by FimA type variations. Oral Microbiol. Immunol. 24, 260–263. doi: 10.1111/j.1399-302X.2009.00511.x
Kulkarni, P. G., Gosavi, S., Haricharan, P. B., Malgikar, S., Mudrakola, D. P., Turagam, N., et al. (2018). Molecular detection of Porphyromonas gingivalis in chronic periodontitis patients. J. Contemp. Dental Pract. 19, 992–996. doi: 10.5005/jp-journals-10024-2371
Kumar, P. S. (2017). From focal sepsis to periodontal medicine: a century of exploring the role of the oral microbiome in systemic disease. J. Physiol. 595, 465–476. doi: 10.1113/JP272427
Kumawat, R. M., Ganvir, S. M., Hazarey, V. K., Qureshi, A., and Purohit, H. J. (2016). Detection of Porphyromonas gingivalis and Treponema denticola in chronic and aggressive periodontitis patients: a comparative polymerase chain reaction study. Contemp. Clin. Dentistry 7, 481–486. doi: 10.4103/0976-237X.194097
Laine, M. L., Appelmelk, B. J., and van Winkelhoff, A. J. (1996). Novel polysaccharide capsular serotypes in Porphyromonas gingivalis. J. Periodontal. Res. 31, 278–284. doi: 10.1111/j.1600-0765.1996.tb00494.x
Laine, M. L., and van Winkelhoff, A. J. (1998). Virulence of six capsular serotypes of Porphyromonas gingivalis in a mouse model. Oral Microbiol. Immunol. 13, 322–325. doi: 10.1111/j.1399-302X.1998.tb00714.x
Lamont, R. J., and Jenkinson, H. F. (2000). Subgingival colonization by Porphyromonas gingivalis. Oral Microbiol. Immunol. 15, 341–349. doi: 10.1034/j.1399-302x.2000.150601.x
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9:357. doi: 10.1038/nmeth.1923
Lau, L., Sanz, M., Herrera, D., Morillo, J. M., Martín, C., and Silva, A. (2004). Quantitative real-time polymerase chain reaction versus culture: a comparison between two methods for the detection and quantification of Actinobacillus actinomycetemcomitans, Porphyromonas gingivalis and Tannerella forsythensis in subgingival plaque samples. J. Clin. Periodontol. 31, 1061–1069. doi: 10.1111/j.1600-051X.2004.00616.x
Lépine, G., and Progulske-Fox, A. (1996). Duplication and differential expression of hemagglutinin genes in Porphyromonas gingivalis. Oral Microbiol. Immunol. 11, 65–78. doi: 10.1111/j.1399-302X.1996.tb00339.x
Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079. doi: 10.1093/bioinformatics/btp352
Li, N., and Collyer, C. A. (2011). Gingipains from Porphyromonas gingivalis- Complex domain structures confer diverse functions. Eur. J. Microbiol. Immunol. 1, 41–58. doi: 10.1556/EuJMI.1.2011.1.7
Li, N., Yun, P., Jeffries, C. M., Langley, D., Gamsjaeger, R., Church, W. B., et al. (2011). The modular structure of haemagglutinin/adhesin regions in gingipains of Porphyromonas gingivalis. Mol. Microbiol. 81, 1358–1373. doi: 10.1111/j.1365-2958.2011.07768.x
Loos, B. G., and Dyer, D. W. (1992). Restriction fragment length polymorphism analysis of the fimbrillin locus, fimA, of Porphyromonas gingivalis. J. Dent. Res. 71, 1173–1181. doi: 10.1177/00220345920710050901
Loos, B. G., Dyer, D. W., Whittam, T. S., and Selander, R. K. (1993). Genetic structure of populations of Porphyromonas gingivalis associated with periodontitis and other oral infections. Infect. Immun. 61, 204–212.
Maisetta, G., Brancatisano, F. L., Esin, S., Campa, M., and Batoni, G. (2011). Gingipains produced by Porphyromonas gingivalis ATCC49417 degrade human-beta-defensin 3 and affect peptide's antibacterial activity in vitro. Peptides 32, 1073–1077. doi: 10.1016/j.peptides.2011.02.003
Ménard, C., Brousseau, R., and Mouton, C. (1992). Application of polymerase chain reaction with arbitrary primer (AP-PCR) to strain identification of Porphyromonas (Bacteroides) gingivalis. FEMS Microbiol. Lett. 74, 163–168. doi: 10.1111/j.1574-6968.1992.tb05360.x
Miller, D. P., Hutcherson, J. A., Wang, Y., Nowakowska, Z. M., Potempa, J., Yoder-Himes, D. R., et al. (2017). Genes contributing to Porphyromonas gingivalis fitness in abscess and epithelial cell colonization environments. Front. Cell. Infect. Microbiol. 7:378. doi: 10.3389/fcimb.2017.00378
Missailidis, C. G., Umeda, J. E., Ota-Tsuzuki, C., Anzai, D., and Mayer, M. P. (2004). Distribution of fimA genotypes of Porphyromonas gingivalis in subjects with various periodontal conditions. Oral Microbiol. Immunol. 19, 224–229. doi: 10.1111/j.1399-302X.2004.00140.x
Naito, M., Hirakawa, H., Yamashita, A., Ohara, N., Shoji, M., Yukitake, H., et al. (2008). 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. 15, 215–225. doi: 10.1093/dnares/dsn013
Nakagawa, I., Amano, A., Kimura, R. K., Nakamura, T., Kawabata, S., and Hamada, S. (2000). Distribution and molecular characterization of Porphyromonas gingivalis carrying a new type of fimA gene. J. Clin. Microbiol. 38, 1909–1914.
Nakagawa, I., Amano, A., Kuboniwa, M., Nakamura, T., Kawabata, S., and 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, 277–285. doi: 10.1128/IAI.70.1.277-285.2002
Narimatsu, M., Noiri, Y., Itoh, S., Noguchi, N., Kawahara, T., and Ebisu, S. (2004). Essential role for the gtfA gene encoding a putative glycosyltransferase in the adherence of Porphyromonas gingivalis. Infect. Immun. 72, 2698–2702. doi: 10.1128/IAI.72.5.2698-2702.2004
Nelson, K. E., Fleischmann, R. D., DeBoy, R. T., Paulsen, I. T., Fouts, D. E., Eisen, J. A., et al. (2003). Complete genome sequence of the oral pathogenic Bacterium Porphyromonas gingivalis strain W83. J. Bacteriol. 185, 5591–5601. doi: 10.1128/JB.185.18.5591-5601.2003
Page, A. J., Cummins, C. A., Hunt, M., Wong, V. K., Reuter, S., Holden, M. T., et al. (2015). Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics 31, 3691–3693. doi: 10.1093/bioinformatics/btv421
Park, O. J., Yi, H., Jeon, J. H., Kang, S. S., Koo, K. T., Kum, K. Y., et al. (2015). Pyrosequencing analysis of subgingival microbiota in distinct periodontal conditions. J. Dent. Res. 94, 921–927. doi: 10.1177/0022034515583531
Patel, R. K., and Jain, M. (2012). NGS QC Toolkit: a toolkit for quality control of next generation sequencing data. PLoS ONE 7:e30619. doi: 10.1371/journal.pone.0030619
Pérez-Chaparro, P. J., Gracieux, P., Lafaurie, G. I., Donnio, P. Y., and Bonnaure-Mallet, M. (2008). Genotypic characterization of Porphyromonas gingivalis isolated from subgingival plaque and blood sample in positive bacteremia subjects with periodontitis. J. Clin. Periodontol. 35, 748–753. doi: 10.1111/j.1600-051X.2008.01296.x
Petersen, P. E., and Ogawa, H. (2012). The global burden of periodontal disease: towards integration with chronic disease prevention and control. Periodontol 60, 15–39. doi: 10.1111/j.1600-0757.2011.00425.x
Pritchard, L., Glover, R. H., Humphris, S., Elphinstone, J. G., and Toth, I. K. (2016). Genomics and taxonomy in diagnostics for food security: soft-rotting enterobacterial plant pathogens. Anal. Methods 8, 12–24. doi: 10.1039/C5AY02550H
Progulske-Fox, A., Kozarov, E., Dorn, B., Dunn, W., Burks, J., and Wu, Y. (1999). Porphyromonas gingivalis virulence factors and invasion of cells of the cardiovascular system. J. Periodontal. Res. 34, 393–399. doi: 10.1111/j.1600-0765.1999.tb02272.x
Progulske-Fox, A., Tumwasorn, S., and Holt, S. C. (1989). The expression and function of a Bacteroides gingivalis hemagglutinin gene in Escherichia coli. Oral Microbiol. Immunol. 4, 121–131. doi: 10.1111/j.1399-302X.1989.tb00238.x
Progulske-Fox, A., Tumwasorn, S., Lépine, G., Whitlock, J., Savett, D., Ferretti, J. J., et al. (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, 311–318. doi: 10.1111/j.1399-302X.1995.tb00160.x
Puth, S., Hong, S. H., Na, H. S., Lee, H. H., Lee, Y. S., Kim, S. Y., et al. (2019). A built-in adjuvant-engineered mucosal vaccine against dysbiotic periodontal diseases. Mucosal Immunol. 12, 565–579. doi: 10.1038/s41385-018-0104-6
Roux, D., Cywes-Bentley, C., Zhang, Y. F., Pons, S., Konkol, M., Kearns, D. B., et al. (2015). Identification of poly-N-acetylglucosamine as a major polysaccharide component of the Bacillus subtilis biofilm matrix. J. Biol. Chem. 290, 19261–19272. doi: 10.1074/jbc.M115.648709
Schmieder, R., and Edwards, R. (2011). Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864. doi: 10.1093/bioinformatics/btr026
Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153
Shoji, M., Yukitake, H., Sato, K., Shibata, Y., Naito, M., Aduse-Opoku, J., et al. (2013). Identification of an O-antigen chain length regulator, WzzP, in Porphyromonas gingivalis. MicrobiologyOpen 2, 383–401. doi: 10.1002/mbo3.84
Song, H., Bélanger, M., Whitlock, J., Kozarov, E., and Progulske-Fox, A. (2005). Hemagglutinin B is involved in the adherence of Porphyromonas gingivalis to human coronary artery endothelial cells. Infect. Immun. 73, 7267–7273. doi: 10.1128/IAI.73.11.7267-7273.2005
Soto, C., Bugueño, I., Hoare, A., Gonzalez, S., Venegas, D., Salinas, D., et al. (2016). The Porphyromonas gingivalis O antigen is required for inhibition of apoptosis in gingival epithelial cells following bacterial infection. J. Periodontal. Res. 51, 518–528. doi: 10.1111/jre.12331
Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033
Stathopoulou, P. G., Galicia, J. C., Benakanakere, M. R., Garcia, C. A., Potempa, J., and Kinane, D. F. (2009). Porphyromonas gingivalis induce apoptosis in human gingival epithelial cells through a gingipain-dependent mechanism. BMC Microbiol. 9:107. doi: 10.1186/1471-2180-9-107
Tancharoen, S., Matsuyama, T., Kawahara, K., Tanaka, K., Lee, L. J., Machigashira, M., et al. (2015). Cleavage of host cytokeratin-6 by lysine-specific gingipain induces gingival inflammation in periodontitis patients. PLoS ONE 10:e0117775. doi: 10.1371/journal.pone.0117775
To, T. T., Liu, Q., Watling, M., Bumgarner, R. E., Darveau, R. P., and McLean, J. S. (2016). Draft genome sequence of low-passage clinical isolate Porphyromonas gingivalis MP4-504. Genome Announcements 4:16. doi: 10.1128/genomeA.00256-16
van Kranenburg, R., van Swam, I. I., Marugg, J. D., Kleerebezem, M., and de Vos, W. M. (1999). Exopolysaccharide biosynthesis in Lactococcus lactis NIZO B40: functional analysis of the glycosyltransferase genes involved in synthesis of the polysaccharide backbone. J. Bacteriol. 181, 338–340.
van Winkelhoff, A. J., Appelmelk, B. J., Kippuw, N., and de Graaff, J. (1993). K-antigens in Porphyromonas gingivalis are associated with virulence. Oral Microbiol. Immunol. 8, 259–265. doi: 10.1111/j.1399-302X.1993.tb00571.x
Waterhouse, R. M., Seppey, M., Simão, F. A., Manni, M., Ioannidis, P., Klioutchnikov, G., et al. (2017). BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol. Biol. Evol. 35, 543–548. doi: 10.1093/molbev/msx319
Weinberg, A., Belton, C. M., Park, Y., and Lamont, R. J. (1997). Role of fimbriae in Porphyromonas gingivalis invasion of gingival epithelial cells. Infect. Immun. 65, 313–316.
Wu, L., Zhao, L., Wang, J., Liu, C., Li, Y., and Wu, Y. (2018). pckA-deficient Porphyromonas gingivalis W83 shows reduction in hemagglutination activity and alteration in the distribution of gingipain activity. Eur. J. Oral Sci. 126, 359–366. doi: 10.1111/eos.12565
Yamamoto, R., Noiri, Y., Yamaguchi, M., Asahi, Y., Maezono, H., and Ebisu, S. (2011). Time course of gene expression during Porphyromonas gingivalis strain ATCC 33277 biofilm formation. Appl. Environ. Microbiol. 77, 6733–6736. doi: 10.1128/AEM.00746-11
Yuan, L., Rodrigues, P. H., Bélanger, M., Dunn, W., and Progulske-Fox, A. (2007). The Porphyromonas gingivalis clpB gene is involved in cellular invasion in vitro and virulence in vivo. FEMS Immunol. Med. Microbiol. 51, 388–398. doi: 10.1111/j.1574-695X.2007.00326.x
Yuan, L., Rodrigues, P. H., Bélanger, M., Dunn, W. A., and Progulske-Fox, A. (2008). Porphyromonas gingivalis htrA is involved in cellular invasion and in vivo survival. Microbiology 154, 1161–1169. doi: 10.1099/mic.0.2007/015131-0
Keywords: Porphyromonas gingivalis, periodontitis, comparative genomics, hemagglutination, fimbrium, biofilm
Citation: Mendez KN, Hoare A, Soto C, Bugueño I, Olivera M, Meneses C, Pérez-Donoso JM, Castro-Nallar E and Bravo D (2019) Variability in Genomic and Virulent Properties of Porphyromonas gingivalis Strains Isolated From Healthy and Severe Chronic Periodontitis Individuals. Front. Cell. Infect. Microbiol. 9:246. doi: 10.3389/fcimb.2019.00246
Received: 11 April 2019; Accepted: 24 June 2019;
Published: 10 July 2019.
Edited by:
Ulvi Kahraman Gürsoy, University of Turku, FinlandReviewed by:
J. Christopher Fenno, University of Michigan, United StatesAshu Sharma, University at Buffalo, United States
Nursen Topcuoglu, Istanbul University, Turkey
Copyright © 2019 Mendez, Hoare, Soto, Bugueño, Olivera, Meneses, Pérez-Donoso, Castro-Nallar and Bravo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eduardo Castro-Nallar, ZWR1YXJkby5jYXN0cm8mI3gwMDA0MDt1bmFiLmNs; Denisse Bravo, ZGVuYnJhdm8mI3gwMDA0MDt1Y2hpbGUuY2w=
†These authors have contributed equally to this work