*Article* **Prospective Longitudinal Changes in the Periodontal Inflamed Surface Area Following Active Periodontal Treatment for Chronic Periodontitis**

**Yoshiaki Nomura <sup>1</sup> , Toshiya Morozumi 2,\* , Atsushi Saito <sup>3</sup> , Atsutoshi Yoshimura <sup>4</sup> , Erika Kakuta 5, Fumihiko Suzuki <sup>6</sup> , Fusanori Nishimura 7, Hideki Takai 8, Hiroaki Kobayashi 9, Kazuyuki Noguchi 10, Keiso Takahashi 11, Koichi Tabeta 12, Makoto Umeda 13, Masato Minabe 2, Mitsuo Fukuda 14, Naoyuki Sugano 15, Nobuhiro Hanada <sup>1</sup> , Nobuo Yoshinari <sup>16</sup> , Satoshi Sekino 17, Shogo Takashiba 18, Soh Sato 19, Toshiaki Nakamura <sup>10</sup> , Tsutomu Sugaya 20, Yohei Nakayama 8, Yorimasa Ogata <sup>8</sup> , Yukihiro Numabe <sup>17</sup> and Taneaki Nakagawa <sup>21</sup>**


**Abstract:** Periodontal disease is a chronic inflammatory disease of the periodontal tissue. The periodontal inflamed surface area (PISA) is a proposed index for quantifying the inflammatory burden resulting

**Citation:** Nomura, Y.; Morozumi, T.; Saito, A.; Yoshimura, A.; Kakuta, E.; Suzuki, F.; Nishimura, F.; Takai, H.; Kobayashi, H.; Noguchi, K.; et al. Prospective Longitudinal Changes in the Periodontal Inflamed Surface Area Following Active Periodontal Treatment for Chronic Periodontitis. *J. Clin. Med.* **2021**, *10*, 1165. https:// doi.org/10.3390/jcm10061165

Academic Editor: Denis Bourgeois

Received: 22 February 2021 Accepted: 8 March 2021 Published: 10 March 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

from periodontitis lesions. This study aimed to investigate longitudinal changes in the periodontal status as evaluated by the PISA following the active periodontal treatment. To elucidate the prognostic factors of PISA, mixed-effect modeling was performed for clinical parameters, tooth-type, and levels of periodontal pathogens as independent variables. One-hundred-twenty-five patients with chronic periodontitis who completed the active periodontal treatment were followed-up for 24 months, with evaluations conducted at 6-month intervals. Five-times repeated measures of mean PISA values were 130+/−173, 161+/−276, 184+/−320, 175+/−417, and 209+/−469 mm2. Changes in clinical parameters and salivary and subgingival periodontal pathogens were analyzed by mixed-effect modeling. Plaque index, clinical attachment level, and salivary levels of *Porphyromonas gingivalis* were associated with changes in PISA at the patient- and tooth-level. Subgingival levels of *P. gingivalis* and *Prevotella intermedia* were associated with changes in PISA at the sample site. For most patients, changes in PISA were within 10% of baseline during the 24-month follow-up. However, an increase in the number of bleeding sites in a tooth with a deep periodontal pocket increased the PISA value exponentially.

**Keywords:** periodontal inflamed surface area; periodontal pathogen; mixed effect modeling; follow-up study

### **1. Introduction**

Periodontal disease is a chronic inflammatory disease of the periodontal tissue. Its progression has been associated with a regression of the mean clinical parameters of most of the oral sites, including teeth. Most patients are generally stable or exhibit a linear type of progression [1–4]. However, a small fraction exhibit burst-type progressions [1,5,6].

The active periodontal treatment aims to reduce inflammatory response by removing pathogenic bacterial deposits. Following the active treatment, supportive periodontal therapy (SPT) is employed to reduce the probability of periodontal disease progression. The long-term successful SPT prevents tooth loss [7]. The reported rate of tooth loss is 10% of teeth [8–15] in 20% of patients [16] over a period of 10 years. Further, tooth loss has been shown to occur in a small fraction of patients during SPT [8,17]. Therefore, monitoring of periodontal conditions during SPT is important for a successful treatment [18]. Several indices have been routinely used for this purpose. The probing pocket depth (PD) and clinical attachment level (CAL) of the periodontal pocket are morphological outcomes of periodontal disease. However, these outcomes do not quantify the proportion of inflamed periodontal tissue. Additionally, summary statistics of these outcomes using the mean or maximum values results in information loss [1]. The mean values of these outcomes obscure the small number of sites that exhibit a substantial progress, and the maximum value represents only one of 168 probing sites. Therefore, site-level evaluation is important in the clinical setting. During the follow-up period, the periodontal tissue should be evaluated in terms of the morphological and inflammatory response status.

The periodontal inflamed surface area (PISA) is a convenient index that quantifies the surface area of the bleeding pocket epithelium in square millimeters and is calculated using conventional clinical parameters of periodontal health, namely, bleeding on probing (BOP) combined with either PD or CAL, and gingival recession [19]. Thus, PISA may represent the inflammation status of the patients as one continuous variable, in which the calculation of mean or maximum value is not necessary. Several studies have utilized the PISA for evaluating the periodontal status. However, these studies have focused on the cross-sectional correlation of the PISA and non-communicable diseases [20–28]. There have only been a few longitudinal studies regarding the changes in PISA [29].

This study aimed to investigate the changes in periodontal status, as evaluated by the PISA during the longitudinal follow-up care after an active treatment, and to identify the factors that affect the PISA.

### **2. Materials and Methods**

### *2.1. Study Design and Ethics Approval*

This study is part of a clinical research project performed by The Japanese Society of Periodontology. A 24-month follow-up study was performed with 163 patients who had completed an active periodontal treatment (such as initial therapy or periodontal surgery) at 17 facilities in Japan. The end point of those active periodontal treatments was periodontal examination. If the lesion was considered to be dormant, such as a probing pocket depth of 4 mm or more without BOP, the disease was considered stable [30]. In addition, those patients were the candidates for this study. The follow-up patients were seen trimonthly, and the supragingival plaque and calculus were removed if detected. The details of changes in the clinical parameters during the follow-ups have been described in our pervious report [1,31,32]. The study was conducted in compliance with the principles outlined in the Helsinki Declaration. A written informed consent was obtained from each study participant, and the protocol was approved by the Institutional Review Board of each participating institution.

### *2.2. Inclusion Criteria*

The inclusion criteria were as follows: The severity of periodontitis on the initial visit was generalized moderate–to-severe chronic periodontitis (generalized periodontitis, stage III or IV, grade B) [33,34], patients who were systemically healthy; aged >30 years; >20 remaining teeth; and not taking antibiotics, anti-inflammatory drugs or immunosuppressive drugs for 3 months before the start of follow-up.

### *2.3. Diagnosis and Evaluation*

Each patient was diagnosed according to the Guidelines of the American Academy of Periodontology [33]. Oral examinations were carried out by one examiner from each institute (T.M., A.Y., F.S., H.T., H.K., K.N., M.M., M.F., N.S., N.Y., S.S., S.S., T.N., T.S., Y.O., Y.N. and T.N.). Each examiner was a periodontist licensed by the Japanese Society of Periodontology. Intra- and inter-examiner calibration sessions were conducted using periodontal disease models (P15FE-500HPRO-S2A1-GSF) at the beginning and middle of the study period. In brief, full-mouth PD and gingival recession were measured twice, and repeatability for the CAL was assessed. The examiner was judged to have made reproducible measurements after reaching a percentage of agreement within ±1 mm between repeated measurements of at least 95% of measurements.

The PD, BOP, and CAL were measured at six sites per tooth (mesiobuccal, buccal, distobuccal, mesiolingual, lingual, and distolingual) with a periodontal probe (CP-12 Color-Coded Probe; Hu-Friedy, Chicago, IL, USA). The PISA was calculated using the Excel sheet program [19]. By inputting the BOP and PD data from six sites per tooth, the PISA values of each tooth were calculated. The sum of the PISA values is automatically calculated. This value is the summary statistics of PISA for each individual.

The plaque index (PlI) was recorded at four sites per tooth (mesial, buccal, distal, and lingual). The degree of tooth mobility was scored on a four-point scale (0–3) at the tooth level. All clinical parameters were recorded at baseline and after 6, 12, 18, and 24 months, the third molars were excluded.

### *2.4. Evaluation of Periodontal Pathogens*

Three periodontal pathogens in the saliva and subgingival plaque were measured: *Aggregatibacter actinomycetemcomitans*, *Porphyromonas gingivalis*, and *Prevotella intermedia*. These samples were collected from each patient at every visit. Briefly, whole saliva was collected by asking the patient to chew on a gum base for 5 min [31], and the subgingival plaque from the deepest pockets (except for those at the third molars) was obtained by consecutive insertion of two sterile paper points into the periodontal pocket for 10 s per point [32]. These three periodontal pathogens and the total bacteria in the subgingival plaque and saliva were counted using a modification of the Invader PLUS assay [35,36]. Bacterial ratios (%) for each species were also determined.

Briefly, a quantitative analysis of the total bacterial count and periodontopathic bacterial counts, including *P. gingivalis*, *P. intermedia*, and *A. actinomycetemcomitans*, was performed using a modification of the Invader PLUS assay [31,32,35,36]. The bacterial DNA was extracted from the subgingival plaque samples from the deepest pockets by suspending each plaque sample in 1 mL of phosphate-buffered saline (pH 7.4) and processed using the MagNA Pure LC Total Nucleic Acid Isolation Kit (Roche, Basel, Switzerland) according to the manufacturer's instructions. Similarly, bacterial DNA was extracted from the 100-μL whole saliva samples using the MagNA Pure LC Total Nucleic Acid Isolation Kit. The individual sequences of each bacterial species were obtained from a public database (National Center for Biotechnology Information, Bethesda, MD, USA). Primers for each species were designed based on a region of the 16S rRNA gene. A pair of universal primers and a universal probe were used to determine the total number of bacteria. Primary probes and Invader oligos were designed using the Invader technology creator (HOLOGIC, Madison, WI, USA) and were based on sequences in the amplified regions [31,32,35–37].

The template DNA was added to a 15-μL reaction mixture containing primers for each species [50 μM deoxynucleoside triphosphate (dNTP), 700 nM primary probe, 70 nM Invader oligo, 2.5 U polymerase chain reaction (PCR) enzyme (EagleTaq DNA polymerase, Roche, Basel, Switzerland), and the Invader core reagent kit (Cleavase XI Invader core reagent kit, HOLOGIC, Madison, WI, USA) containing a fluorescence resonance energy transfer (FRET) mix and an enzyme/MgCl2 solution]. The reaction mixture was preheated at 95 ◦C for 20 min, and a two-step polymerase chain reaction (95 ◦C for 1 s and 63 ◦C for 1 min) was performed for 35 cycles using an ABI PRISM 7900 thermocycler (Applied Biosystems, Foster City, CA, USA). The fluorescence values of carboxyfluorescein (FAM; wavelength/bandwidth: Excitation, 485/20 nm; emission, 530/25 nm) were measured at the end of the incubation/extension step at 63 ◦C for each cycle.

The limit of detection for this method was determined for each species with dilutions of bacterial DNA. Standard curves were constructed based on a crossing point determined by the fit-point method.

The proportions of the three pathogens within the total bacterial count were subsequently determined [38,39]. Bacterial ratios (%) and bacterial counts (log10) for each species were also used in various comparison and diagnostic analyses.

### *2.5. Statistical Analysis*

### 2.5.1. Mixed-Effect Modeling

Changes in the PISA were analyzed at the patient-level and tooth-level by mixedeffect modeling [1,40–44]. The data structure of the PISA was five-time repeated measures. The PISA values for each tooth had a hierarchical structure in the patients. Generalized mixed-effect modeling was applied using the following formula:

Model 1: Patient-level PISA

−L1 : *PISA*(*Subject level*)*<sup>i</sup>* =

$$\begin{aligned} & \pi\_0 + \pi\_1 (A \g e)\_l + \pi\_2 (S \epsilon \alpha)\_l + \pi\_4 (S \text{alinary levels of } A.a)\_l + \pi\_5 (S \text{alinary levels of } P.g.)\_l \\ & + \pi\_6 (S \text{alinary levels of } P.i)\_l + \pi\_7 (mean \text{CAL})\_l + \pi\_8 (Pli)\_l + \sum\_{w=1}^5 \pi\_9^{(w)} (time)\_j + \varepsilon\_l \end{aligned} \tag{1}$$

$$e\_l \sim N\left(0, \delta\_e^2\right).$$

Fixed effect;

Patient-level: Age, sex, time, salivary levels of *A. actinomycetemcomitans*, *P. gingivalis*, and *P. intermedia*, mean CAL, PlI, time Random effect: Patients, random intercept Covariance Type: AR1 Link functions: Gamma

### Model 2: Tooth-level PISA

−*L*1 : *PISA*(*Tooth level*)*ij* =

$$\begin{cases} \pi\_{0j} + \pi\_{1lj}(A\varrho\varepsilon)\_{lj} + \pi\_{2j}(\text{Sex})\_{lj} + \pi\_{3j}(\text{Saliinary levels of A. a})\_{lj} + \pi\_{4j}(\text{Saliinary levels of P. g})\_{ij} \\ \eta + \pi\_{5j} \text{ (Saliinary levels of P. i)}\_{ij} + \varepsilon\_{ij} \end{cases} \tag{1}$$

$$-L2: \pi\_{0j} = \beta\_{0j} + \beta\_{1j}(Mean\ CAL) + \beta\_{2j}(Pli) + \sum\_{m=1}^{3} \beta\_{3j}^{(m)}(troch\ mobility)\_{j} + \sum\_{m=1}^{6} \beta\_{4j}^{(m)}(Tool\ types)\_{j} + \sum\_{m=1}^{5} \beta\_{05}^{(m)}(tû\#)\_{j} + r\_{j} \tag{2}$$

$$e\_{lj} \sim N\left(0, \delta\_c^2\right), r\_{oj} \sim N\left(0, \delta\_r^2\right).$$

Fixed effect;

Patient-level: Age, sex, days, salivary levels of *A. actinomycetemcomitans*, *P. gingivalis*, and *P. intermedia* Tooth-level: Mean CAL, PlI, tooth mobility, tooth type, time Random effect: Patients, random intercept Covariance Type: AR1 Link functions: Gamma

### Model 3: PISA of the sampling teeth for subgingival periodontal pathogens

−L1 : *PISA*(*Sampling teeth*)*<sup>i</sup>* = *π*<sup>0</sup> + *π*1(% *of A*.*a in subgingival taotal bacteria*)*<sup>i</sup>* +*π*2(% *of A*.*a in subgingival taotal bacteria*)*<sup>i</sup>* + *π*3(% *of P*. *g in subgingival taotal bacteria* )*<sup>i</sup>* <sup>+</sup>*π*4(% *of P*. *i in subgingival taotal bacteria*)*<sup>i</sup>* <sup>+</sup> <sup>5</sup> ∑ *m*=1 *<sup>π</sup>*(*m*) <sup>5</sup> (*time*)*<sup>j</sup>* + *ε<sup>i</sup>* (1)

$$e\_l \sim N\left(0, \delta\_e^2\right).$$

Fixed effect;

Patient-level: Proportion of *A. actinomycetemcomitans*, *P. gingivalis*, and *P. intermedia* in the total subgingival bacteria, time Random effect: Patients Covariance Type: AR1

Link functions: Gamma

### 2.5.2. Cluster Analysis

A cluster analysis was carried out using K-means. The data used for the cluster analysis were sequential data obtained from five-time repeated measures of PISA at the patient level. The parameter settings for the cluster analysis were: Distance measure: Log likelihood, and the criteria of clusters were Bayesian Information Criterion (BIC) [−2 ln (L) + k ln (L)]. The generated cluster was used for the analysis.

All analyses were performed using SPSS Statistics version 24.0 (IBM, Tokyo, Japan).

### **3. Results**

### *3.1. Patient Characteristics*

Among the 163 patients enrolled in the study, 38 patients were disqualified for lack of data due to a variety of causes such as earthquake, failure to comply with visit, use of antimicrobial agents for acute periodontal abscesses, and tooth extraction for root fracture. Thus, data from a total of 125 patients, 3107 teeth, and their five-time repeated measure PISA values were obtained. A total of 15,535 data were analyzed from the teeth. The study population comprised 50 men and 75 women and the mean age was 59.2 ± 8.7 years. Clinical examinations and measurement of periodontal pathogens were carried out at every visit. Five-time repeated measures of these parameters were analyzed.

### *3.2. Analysis of the Change in PISA at the Patient-Level*

Factors correlated with the PISA were identified using generalized mixed-effect modeling. The results are shown in Table 1. The salivary levels of *P*. *gingivalis*, mean CAL, and PlI were identified as significant factors affecting the PISA. Time was not statistically significant.

**Table 1.** Patient-level analysis of factors affecting the periodontal inflamed surface area (PISA).


Fitness index, BIC: 1317; AICC, 1326; PlI: Plaque index; CAL: Clinical attachment level; *A. a*: *Aggregatibacter actinomycetemcomitans*; *P. g*: *Porphyromonas gingivalis*; *P. i: Prevotella intermedia*; CI: Confidence interval; AICC = −2 *L* + 2 *k* (*k* − 1)/(n – *k* − 1);BIC = −2 *L* + *k* ln (n); *L*: Likelihood; *k*: Number of parameters; n: Sample size.

### *3.3. Analysis of the Change in PISA at the Tooth-Level*

The sum of the PISA of each tooth was calculated as a dependent variable. To determine the effect of tooth type, a multilevel mixed-effect model was used for the analysis. As shown in Table 2, the salivary levels of *P. gingivalis* were identified as a significant factor at the patient-level. The CAL, the PlI, and tooth mobility >2 were identified as significant factors at the tooth-level. Regarding the tooth type, the maxillary molar and premolar tended to increase the PISA during the follow-up period. Time was not statistically significant.

### **Table 2.** Tooth-level analysis of factors affecting the PISA.


Fitness index, AICC: 11,732.53; BIC: 11,745.40; PlI: Plaque index; CAL: Clinical attachment level; *A. a*: *Aggregatibacter actinomycetemcomitans*; *P. g: Porphyromonas gingivalis*; *P. i*: *Prevotella intermedia;* AICC = −2 *L* + 2 *k* (*k* − 1)/ (n – *k* − 1); BIC = −2 *L* + *k* ln (n); *L*: Likelihood; *k*: Number of parameters; n: Sample size.

### *3.4. Classification of the Changes in the PISA*

The change in the PISA was dependent on the tooth-type (Figure S1). The PISA of maxillary molars increased throughout the follow-up period, while the PISA of other tooth-types were relatively stable.

A cluster analysis was carried out to characterize the changes in the PISA. The results of the cluster analysis and changes in the PISA in each cluster are shown in Figure 1. Most of the patients were classified as Cluster 1, and the PISA of this cluster was stable. A minority of patients were classified to other clusters, which exhibited drastic changes in the PISA during the follow-up period.

**Figure 1.** Changes in the PISA according to the clusters. For 117 (93.6%) patients, the PISA was stable during the 24-month follow-up period. However, a minority of patients exhibited a drastic increase or fluctuation in the PISA.

A minority of the teeth may have had an effect on the overall PISA. The teeth with PISA >500 were identified and the changes in the PISA of these teeth are illustrated in Figure 2. A total of seven teeth (maxillary 2nd molars or mandibular 2nd molars) had PISA >500 at least once during the follow-up period. The drastic increase in the overall PISA was primarily dependent on the increase in the number of BOP sites.

**Figure 2.** Changes in the PISA at the tooth-level in teeth with PISA > 500 during the 24-month follow-up period. Red lines indicate the maxillary 2nd molars and blue lines indicate mandibular 2nd molars. Seven teeth exceed 500 PISA. The numbers indicate the mean probing depth of each tooth and number of bleeding sites on probing. In these cases, drastic changes in the PISA were predominantly dependent on the increase in the number of bleeding sites on probing.

### *3.5. Association between Periodontal Pathogens and the Change in the PISA*

Mixed-effect modeling was employed to evaluate the local association between periodontal pathogens and the tooth-level PISA. The results are shown in Table 3. The PISA of the sample teeth was included as a dependent variable. The periodontal pathogens, *P. gingivalis* and *P. intermedia*, were observed to be significantly associated with the PISA of the sample teeth. Time was not a significant factor for PISA. The list of the sample teeth is shown in Table S1.


**Table 3.** Association between subgingival periodontal pathogens and the changes in the PISA of the sample sites.

Fitness index: For PISA, AICC: 734, BIC: 741; CAL: Clinical attachment level, *A. a*: *Aggregatibacter actinomycetemcomitans*, *P. g: Porphyromonas gingivalis*, *P. i*: *Prevotella intermedia*, CI: Confidence interval. AICC = − 2 *L* + 2 *k* (*k* − 1)/ (n – *k* − 1), BIC = −2 *L* + *k* ln (n). *L*: Likelihood; *k*: Number of parameters; n: Sample size.

### **4. Discussion**

This study aimed to investigate the changes in the periodontal status, as evaluated by the PISA during a longitudinal follow-up, and to identify the factors that associate with the PISA.

The mixed-effect models used in this study did not include BOP and PD as independent variables as the PISA was calculated using these variables. The PlI and CAL were observed to be significantly associated with the PISA at both the patient-level and tooth-level. Thus, the results indicated that the PISA was associated with the oral hygiene status and the morphology of periodontal tissue. A previous report has indicated that loss of CAL is associated with poor oral hygiene [45], change in the PD is associated with BOP [46], the PlI is associated with tooth loss [47], and changes in the PD are associated with tooth-type and tooth mobility [48]. These findings suggest that changes in the clinical parameters after an active periodontal treatment are interactive, and that deteriorations in these clinical parameters ultimately lead to tooth loss.

As shown in Tables 1–3, time was not significantly associated with the PISA. This finding suggests that the PISA was a stable clinical parameter during the follow-up period. In contrast, it has been previously reported that time is significantly associated with changes in the CAL, with a positive coefficient in the mixed-effect modeling [1]. During the followup period in our study, the CAL deteriorated gradually. In the site-level analysis, most of the CAL measurements were stable. A small fraction of sites with progressive loss of CAL affected the overall CAL. As shown in Figure 1, most of the patients were classified as Cluster 1 and the PISA of this cluster was stable. However, a drastic increase in the PISA in one tooth had a large effect on the PISA at the patient-level. In cases of deep periodontal pockets, the increase in the number of bleeding sites increased the PISA value exponentially. In this sense, the PISA is more sensitive than BOP%. As shown in Figure 2, the drastic increase in the PISA was derived from the increase in bleeding sites rather than deep periodontal pockets. Thus, PISA is a sensitive clinical parameter for identifying the progression of inflamed periodontal tissue.

The salivary levels of *P. gingivalis* were significantly associated with the PISA at the patient-level and tooth-level. Additionally, the subgingival levels of *P. gingivalis* and *P. intermedia* were significantly associated with the PISA of the periodontal pocket of the sample sites. Previous studies have shown that bacteria levels in the subgingival plaque and in whole saliva are also significantly associated with the PISA [49–54]. Further, levels of *P. gingivalis* in the subgingival plaque from the deepest pockets have shown to be associated with the progression of periodontitis [32], and levels of *P. gingivalis* and *P. intermedia* in the subgingival plaque have shown to be associated with the progression of periodontal inflammation [55]. *P. intermedia* is also correlated with BOP and deeper pockets [56,57]. One study reports that reinfection of the treated sites by *P. gingivalis* and/or *P. intermedia* diminishes the effect of therapy during the follow-up [58], and that an increase in periodontal pathogens increases the PISA. Several studies have also indicated that the adjunctive use of antimicrobial agents during the follow-up enhances the effect of mechanical debridement [59]. Therefore, when the PISA increases during the follow-up period, monitoring periodontal pathogen levels and the adjunctive use of antimicrobial agents may be useful for successful long-term treatment.

Our study has some limitations. First, there was no data regarding the initial mean PISA values of the patients before the periodontal treatment. Thus, we were unable to observe any change in the PISA over the entire treatment period of the patient. Second, all patients were systemically healthy. We enrolled patients with periodontitis and without systemic disease as the effect of systemic inflammation can serve as a bias for PISA. Therefore, the association between PISA and the etiology of periodontal disease due to the acquired risk factors was not elucidated. Third, intra- and inter-examiner calibration sessions were conducted using a periodontal disease model rather than in patients due to the large number of examiners required in a multicenter study [60].

### **5. Conclusions**

For most of the patients, changes in the PISA during the follow-up period after an active periodontal treatment were within 10% of the baseline. However, an increase in the number of bleeding sites in a tooth with a deep periodontal pocket was associated with an exponential increase in the PISA.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/2077-038 3/10/6/1165/s1. Figure S1: Changes in PISA by tooth-type during the follow-up period; Table S1: List of sample teeth for the periodontal pathogen.

**Author Contributions:** Conceptualization, Y.N. (Yoshiaki Nomura) and T.M.; methodology, Y.N. (Yoshiaki Nomura); software, Y.N. (Yoshiaki Nomura); validation, E.K. and N.H.; formal analysis, Y.N. (Yoshiaki Nomura); investigation, T.M., A.S., A.Y., F.S., F.N., H.T., H.K., K.N., K.T., M.U., M.M., M.F., N.S., N.Y., S.S. (Satoshi Sekino), S.S. (Soh Sato), T.N. (Toshiaki Nakamura), T.S., Y.N. (Yohei Nakayama), Y.O., Y.N. (Yukihiro Numabe) and T.N. (Taneaki Nakagawa); data curation, S.T.; writing—original draft preparation, Y.N. (Yoshiaki Nomura); writing—review and editing, Y.N. (Yoshiaki Nomura), T.M. and K.T.; visualization, Y.N. (Yoshiaki Nomura); supervision, T.M.; project administration, T.M. and T.N. (Taneaki Nakagawa); funding acquisition, T.N. (Taneaki Nakagawa). All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by a clinical research project grant from the Japanese Society of Periodontology for the diagnosis of periodontitis.

**Institutional Review Board Statement:** The names of the ethics committee members and their reference numbers are as follows: The regional ethical committee of the Faculty of Dentistry, Niigata University (20-R17-08-06); Keio University School of Medicine, Ethics Committee (20080096); Ethical Committee of Kagoshima University Medical and Dental Hospital (20-58); The Ethics Committee, Nagasaki University Graduate School of Biomedical Sciences (0846-2); Ethical Committee of Kyushu University Faculty of Dental Science (20-11); The Ethics Committee of Osaka Dental University (80712); Ethics Committee of Aichi Gakuin University, School of Dentistry (158); The Ethics Committee of Matsumoto Dental University (0090); The Institutional Review Board of Nippon Dental University (2-1-22); Ethical Committee of Nihon University School of Dentistry (EP08D016); Dental Research Ethics Committee of Tokyo Medical and Dental University (660); Ethics Committee in Nihon University School of Dentistry at Matsudo (EC 08-014); Ethics Committee of Tokyo Dental College (208); The Ethical Review Committee of The Nippon Dental University School of Life Dentistry at Niigata (151); Ohu University Research Ethics Committee (52); Institutional Review Board for Clinical Research of Hokkaido University Hospital (008-0113).

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

**Acknowledgments:** The authors thank Hiromasa Yoshie, Toshihide Noguchi, Masamitsu Kawanami, Koichi Ito, Yuichi Izumi, Yoshitaka Hara, Osamu Fujise, Yuzo Abe, Tomoo Kono, Asako Makino-Oi, and Chie Fukaya, for their advice and helpful comments.

**Conflicts of Interest:** The authors declare no conflict of interest.

### **References**

