Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Site-level progression of periodontal disease during a follow-up period

  • Yoshiaki Nomura ,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Software, Writing – original draft, Writing – review & editing

    nomura-y@tsurumi-u.ac.jp

    Affiliation Department of Translational Research, Tsurumi University School of Dental Medicine, Yokohama, Japan

  • Toshiya Morozumi,

    Roles Conceptualization, Data curation, Investigation, Supervision

    Affiliation Division of Periodontology, Department of Oral Biological Science, Niigata University Graduate School of Medical and Dental Sciences, Niigata, Japan

  • Taneaki Nakagawa,

    Roles Conceptualization, Funding acquisition

    Affiliation Department of Dentistry and Oral Surgery, School of Medicine, Keio University, Tokyo, Japan

  • Tsutomu Sugaya,

    Roles Investigation

    Affiliation Division of Periodontology and Endodontology, Department of Oral Health Science, Hokkaido University Graduate School of Dental Medicine, Sapporo, Japan

  • Masamitsu Kawanami,

    Roles Investigation

    Affiliation Division of Periodontology and Endodontology, Department of Oral Health Science, Hokkaido University Graduate School of Dental Medicine, Sapporo, Japan

  • Fumihiko Suzuki,

    Roles Investigation

    Affiliation Division of Dental Anesthesiology, Department of Oral Surgery, Ohu University School of Dentistry, Kooriyama, Japan

  • Keiso Takahashi,

    Roles Investigation

    Affiliation Division of Periodontics, Department of Conservative Dentistry, Ohu University School of Dentistry, Kooriyama, Japan

  • Yuzo Abe,

    Roles Investigation

    Affiliation Comprehensive Dental Care, The Nippon Dental University Niigata Hospital, Niigata, Japan

  • Soh Sato,

    Roles Investigation

    Affiliation Department of Periodontology, School of life Dentistry at Niigata, The Nippon Dental University, Niigata, Japan

  • Asako Makino-Oi,

    Roles Investigation

    Affiliation Department of Periodontology, Tokyo Dental College, Tokyo, Japan

  • Atsushi Saito,

    Roles Investigation

    Affiliation Department of Periodontology, Tokyo Dental College, Tokyo, Japan

  • Satomi Takano,

    Roles Investigation

    Affiliation Bunkyo-Dori Dental Clinic, Chiba, Japan

  • Masato Minabe,

    Roles Investigation

    Affiliations Bunkyo-Dori Dental Clinic, Chiba, Japan, Division of Periodontology, Department of Oral Interdisciplinary Medicine, School of Dentistry, Kanagawa Dental University, Yokosuka, Japan

  • Yohei Nakayama,

    Roles Investigation

    Affiliation Department of Periodontology, Nihon University School of Dentistry at Matsudo, Matsudo, Japan

  • Yorimasa Ogata,

    Roles Investigation

    Affiliation Department of Periodontology, Nihon University School of Dentistry at Matsudo, Matsudo, Japan

  • Hiroaki Kobayashi,

    Roles Investigation

    Affiliation Department of Periodontology, Graduate School of Medical and Dental Sciences, Tokyo Medical and Dental University, Tokyo, Japan

  • Yuichi Izumi,

    Roles Investigation

    Affiliation Department of Periodontology, Graduate School of Medical and Dental Sciences, Tokyo Medical and Dental University, Tokyo, Japan

  • Naoyuki Sugano,

    Roles Investigation

    Affiliation Department of Periodontology, Nihon University School of Dentistry, Tokyo, Japan

  • Koichi Ito,

    Roles Investigation

    Affiliation Nihon University School of Dentistry, Tokyo, Japan

  • Satoshi Sekino,

    Roles Investigation

    Affiliation Department of Periodontology, School of Life Dentistry at Tokyo, The Nippon Dental University, Tokyo, Japan

  • Yukihiro Numabe,

    Roles Investigation

    Affiliation Department of Periodontology, School of Life Dentistry at Tokyo, The Nippon Dental University, Tokyo, Japan

  • Chie Fukaya,

    Roles Investigation

    Affiliation Department of Dentistry and Oral Surgery, School of Medicine, Keio University, Tokyo, Japan

  • Nobuo Yoshinari,

    Roles Investigation

    Affiliation Department of Periodontology, School of Dentistry, Matsumoto Dental University, Shiojiri, Japan

  • Mitsuo Fukuda,

    Roles Investigation

    Affiliation Department of Periodontology, School of Dentistry, Aichi Gakuin University, Nagoya, Japan

  • Toshihide Noguchi,

    Roles Investigation

    Affiliation Department of Periodontology, School of Dentistry, Aichi Gakuin University, Nagoya, Japan

  • Tomoo Kono,

    Roles Investigation

    Affiliation Department of Periodontology, Osaka Dental University, Osaka, Japan

  • Makoto Umeda,

    Roles Investigation

    Affiliation Department of Periodontology, Osaka Dental University, Osaka, Japan

  • Osamu Fujise,

    Roles Investigation

    Affiliation Section of Periodontology, Division of Oral Rehabilitation, Faculty of Dental Science, Kyushu University, Fukuoka, Japan

  • Fusanori Nishimura,

    Roles Investigation

    Affiliation Section of Periodontology, Division of Oral Rehabilitation, Faculty of Dental Science, Kyushu University, Fukuoka, Japan

  • Atsutoshi Yoshimura,

    Roles Investigation

    Affiliation Department of Periodontology, Unit of Translational Medicine, Nagasaki University Graduate School of Biomedical Sciences, Nagasaki, Japan

  • Yoshitaka Hara,

    Roles Investigation

    Affiliation Department of Periodontology, Unit of Translational Medicine, Nagasaki University Graduate School of Biomedical Sciences, Nagasaki, Japan

  • Toshiaki Nakamura,

    Roles Investigation

    Affiliation Department of Periodontology, Kagoshima University Graduate School of Medical and Dental Sciences, Kagoshima, Japan

  • Kazuyuki Noguchi,

    Roles Investigation

    Affiliation Department of Periodontology, Kagoshima University Graduate School of Medical and Dental Sciences, Kagoshima, Japan

  • Erika Kakuta,

    Roles Data curation, Investigation

    Affiliation Department of Oral Microbiology, Tsurumi University School of Dental Medicine, Yokohama, Japan

  • Nobuhiro Hanada,

    Roles Formal analysis, Investigation, Writing – review & editing

    Affiliation Department of Translational Research, Tsurumi University School of Dental Medicine, Yokohama, Japan

  • Shogo Takashiba,

    Roles Conceptualization, Investigation, Methodology

    Affiliation Department of Pathophysiology-Periodontal Science, Okayama University Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama, Japan

  • Yasuharu Amitani,

    Roles Data curation, Formal analysis, Methodology, Validation

    Affiliation Department of Mathematics, Tsurumi University School of Dental Medicine, Yokohama, Japan

  •  [ ... ],
  • Hiromasa Yoshie

    Roles Funding acquisition, Methodology, Project administration, Supervision, Writing – review & editing

    Affiliation Division of Periodontology, Department of Oral Biological Science, Niigata University Graduate School of Medical and Dental Sciences, Niigata, Japan

  • [ view all ]
  • [ view less ]

Abstract

Periodontal disease is assessed and its progression is determined via observations on a site-by-site basis. Periodontal data are complex and structured in multiple levels; thus, applying a summary statistical approach (i.e., the mean) for site-level evaluations results in loss of information. Previous studies have shown the availability of mixed effects modeling. However, clinically beneficial information on the progression of periodontal disease during the follow-up period is not available.

We conducted a multicenter prospective cohort study. Using mixed effects modeling, we analyzed 18,834 sites distributed on 3,139 teeth in 124 patients, and data were collected 5 times over a 24-month follow-up period. The change in the clinical attachment level (CAL) was used as the outcome variable. The CAL at baseline was an important determinant of the CAL changes, which varied widely according to the tooth surface. The salivary levels of periodontal pathogens, such as Porphyromonas gingivalis and Aggregatibacter actinomycetemcomitans, were affected by CAL progression. “Linear”- and “burst”-type patterns of CAL progression occurred simultaneously within the same patient. More than half of the teeth that presented burst-type progression sites also presented linear-type progression sites, and most of the progressions were of the linear type. Maxillary premolars and anterior teeth tended to show burst-type progression. The parameters identified in this study may guide practitioners in determining the type and extent of treatment needed at the site and patient levels. In addition, these results show that prior hypotheses concerning "burst" and "linear" theories are not valid.

Introduction

Periodontal disease is a widespread chronic inflammatory disease caused by the host immunoresponse and periodontal pathogens. Periodontal therapeutic treatments consist of active therapy and subsequent follow-up examinations. Longitudinal studies have shown that periodontal conditions progress during follow-up. A slight increase in the pocket probing depth (PPD) [16], loss of the clinical attachment level (CAL) [2, 3] and loss of teeth [4, 7, 8] have been reported. Previous studies have also provided useful information on the risk factors for the progression of this disease. At the subject level, smoking, compliance with treatment, diabetes mellitus, age [9], and proportion of deep PPD [10, 11] were associated with CAL loss. Moreover, sites with PPD≥6 mm were at a significantly higher risk of CAL loss [10].

One of the important features of periodontal diseases is the site-specific localization of infectious processes that lead to tissue destruction [12]. Assessments of disease progression are performed via observations and measurements on a site-by-site basis. The progression may vary for different sites on a tooth and be affected by properties of the tooth, and it may also be concomitant to the risk factors at the subject level. Therefore, in follow-up examinations, the clinical interest is in the behavior at each site.

Periodontal data are complex and structured in multiple levels (usually three or four), such as the individual-level, tooth-level, site-level and repeated-level measurements. Summary statistical approaches (i.e., mean or sum scores) or the maximum value of the site-level evaluation have been applied for patient-level evaluations. However, aggregating site-level information causes a loss of information [13].

When several levels are involved in one statistical model, site-level observations are not independent because they are nested within the tooth level and patient level. Traditional statistical models tend to overlook multilevel structures and disturbed independency among observations, which leads to type I errors [14] and potential misinterpretations. Evaluations of periodontal disease progression should consider the level hierarchy using analytical models; thus, mixed effects models have been introduced and employed to properly assess multilevel-structured data.

By applying a mixed effects model for periodontal data, valuable information has been obtained on periodontal disease progression. Age and plaque scores affected the bone height [15, 16]. Dental plaque is considered at the site level, tooth mobility is considered at the tooth level and several factors, such as age and socio-economic status, are considered at the patient level, and they affected changes in the PPD [17]. Tooth position affected the loss of CAL [18, 19]. The "linear" and "burst" theories of periodontal disease progression represent manifestations of the same phenomenon [19]. The prognostic factors at the bone level at baseline were tooth mobility, bone level at a site, and PPD at a site [20]. However, the outcome variables varied, and information regarding CAL loss was limited. Tooth loss is a true endpoint for long-term evaluations; however, in short-term observations, CAL loss is a surrogate major endpoint for follow-up therapy.

As described above, previous studies have applied mixed effects models and generated valuable information for clinical practice. However, important information about periodontal disease remains unknown, such as the sites that are predisposed to progression, the risk factors that affect progression and the estimated values of disease progression. Therefore, the aims of this study were to identify the characteristics of sites with CAL loss over a 24-month follow-up period after active periodontal therapy, and the evaluations included the linear and burst theories of periodontal disease progression.

Methods

Study design

Setting.

This study was performed as a clinical research project for the diagnosis of periodontitis by the Japanese Society of Periodontology. In cooperation with 17 facilities (16 university hospitals and one clinic) in Japan, 163 follow-up patients with chronic periodontitis who had finished active treatment regulated by the Japanese Insurance system were registered between February 2009 and February 2012.

Briefly, regulated periodontal treatment consists of 1) Periodontal examination at the first visit, 2) Full-mouth supragingival scaling, 3) Periodontal examination as reevaluation, 4) Subgingival scaling and root planing at sites with a probing depth > 4 mm or a probing depth = 4 mm with BOP, 5) Periodontal examination as reevaluation and periodontal surgery if necessary, and 6) Follow-up. In this study, no patients underwent periodontal surgery.

Diagnosis.

Each diagnosis was based on the Guidelines of the American Academy of Periodontology [21]. Oral examinations were carried out by one examiner at each institute (T.M., C.F., T.S., K.T., Y.A., A.M., M.M., Y.N., H.K., N.S., S.S., N.Y., M.F., T.K., O.F., A.Y., and T.N.). All examiners were periodontists licensed by the Japanese Society of Periodontology.

Patient.

All patients were ≥ 30 years of age, systemically healthy, possessed at least 20 teeth and had not taken systemic antibiotics, anti-inflammatory drugs or immunosuppressive drugs in the 3 months before the start of the investigation. Before periodontal treatment and the initiation of this study, the systemic health conditions of the patients were assessed by questionnaires and medical interviews. All patients were systemically healthy.

Sample size.

The optimal sample size required to detect statistically significant differences in P. gingivalis ratios (%) between patients with or without progression during the supportive periodontal therapy period was determined; based on our previous study [22], we calculated that 77 patients were necessary for each group (α = 0.05 and β = 0.80).

Drop out.

Of the initial 163 follow-up patients, 39 (23.9%) failed to complete the program during the study period. Fourteen patients dropped out following the Great East Japan Earthquake (three facilities), nine failed to comply with visits (six facilities) and eight dropped out for personal reasons (four facilities). Others withdrew because of the use of antimicrobial agents for acute periodontal abscesses (seven patients in five facilities) and because of tooth extraction for root fracture (one patient). In total, 124 subjects successfully completed the study protocol. These 124 patients level analysis are described in our previous report [23, 24].

Research data

The following variables were recorded at baseline and after 6, 12, 18 and 24 months. All data are presented in S1 Table. Descriptive analyses of the variables are shown in S2 Table.

Subject-level explanatory variables.

The demographic background of the patients, such as the age, gender and smoking status, were obtained at baseline. The salivary levels of the total bacteria, Porphyromonas gingivalis and Aggregatibacter actinomycetemcomitans were measured using the Invader PLUS Assay [25, 26]. As shown in S1 Table, bacterial levels showed a skewed distribution. Bacterial variables were dichotomized based on previous reports [22, 2729]. To differentiate sites based on whether they had more or less than 3 mm of CAL progression, cut-off values were set to minimize the difference between sensitivity and specificity [2224, 29]. The cut-off value of A. actinomycetemcomitans at baseline was 0.00006% of Total bacteria (Sensitivity: 0.24, Specificity: 0.86), while that of P. gingivalis was 0.0067% (Sensitivity: 0.58, Specificity: 0.67).

Only one patient was a current smoker (20 cigarettes/day, 25 yrs), and three were previous smokers (5 cigarettes/day, 30 yrs; 10 cigarettes/day, 20 yrs; 10 cigarettes/day, 30 yrs); thus, smoking status was not used as an independent variable.

Tooth-level explanatory variable.

Tooth mobility was evaluated using the standard procedure as described in periodontology textbooks [30, 31]. The assessments were classified as follows:

  1. Degree 1: Mobility of the crown of the tooth, 0.2–1 mm in the horizontal direction
  2. Degree 2: Mobility of the crown of the tooth exceeding 1 mm in horizontal direction
  3. Degree 3: Mobility of the crown of tooth in the vertical direction

Because tooth mobility level 3 was only present in one tooth, levels 2 and 3 were combined.

The plaque index (PlI) was recorded at four sites per tooth: mesial, buccal, distal and lingual.

The mean for each tooth was used for the tooth-level variable [31].

Site-level explanatory variable.

The CAL and bleeding on probing (BOP) were recorded at six sites per tooth (mesiobuccal, buccal, distobuccal, mesiolingual, lingual and distolingual). The tooth surfaces were reclassified into eight categories, which included combinations of maxillary or mandibular and buccal (labial), palatal (lingual), approximal or distal. In this study, third molars were included to distinguish whether the distal surface of the 2nd molar is approximal or distal.

The final data consisted of 18,834 sites distributed on 3,139 teeth in 124 patients, and these data were collected 5 times throughout the study period.

Outcome variables.

The CAL (mm) and differences in the CAL between the baseline and after 24 months were used as outcome variables (ΔCAL). Because a change in CAL ≥ 3 mm is a generally accepted threshold for the progression of CAL [27, 28, 32, 33], he changes over 24 months were classified into six categories: a difference of ≤ -3 mm was defined as improved, between -3 mm and -2 mm was defined as slightly improved, between -1 mm and 1 mm was defined as stable, between 1 mm and 2 mm was defined as slightly progressed, and ≥ 3 mm was defined as progressed. In addition, sites with both ≤ -3 mm improvement and ≥ 3 mm progression were defined as fluctuated.

Multilevel modeling

Periodontal data in this study consisted of a three-level hierarchical structure: patient level, tooth level and site level. The salivary levels of periodontal pathogens, such as A. actinomycetemcomitans and P. gingivalis, were treated as patient-level data. Tooth mobility and PlI were treated as tooth-level data. Clinical parameters with or without BOP, CAL at baseline, tooth type, type of tooth surface and the combination of tooth type and surface were treated as site-level data. P-values less than 0.05 were considered statistically significant.

According to the properties of the outcome variables and analyses, file layouts were changed for each analysis. Three data sets were used for the mixed effects modeling. For the first data set (Model 1), the difference in CAL between the baseline and after 24 months (ΔCAL) was used as the outcome variable, baseline data were used as the explanatory variables, and patient and tooth data were used for the grouping variables.

For the second data set (Model 2), the patterns of CAL changes were analyzed using the 5 CAL measurements as the outcome variables, and the 5 measurements of clinical parameters and salivary levels of periodontal pathogens as explanatory variables. A mixed effects model was constructed using the repeated measures design. For CAL values, a Gamma distribution was chosen as the most suitable option for probability distribution. For this analysis, time was an independent variable. Therefore, all variables were reconstituted by the site and time. Therefore, the total number of samples was 94,170 (n = 41690, 28005 and 24475 for CAL at baseline < 3 mm, 3 mm and > 3 mm, respectively).

For the final data set (Model 3), variables that affect the progression type were analyzed. Before performing mixed effects modeling, a hierarchical cluster analysis of the variables was performed using all 5 measured CAL values. The generated clusters were reclassified into burst- or linear-type progression. Dichotomous variables of these two burst or linear types were used as the outcome variables. A multilevel logistic regression analysis was performed by applying generalized mixed effects modeling. For this analysis, sites that showed burst- or linear type progression were selected from the previous data set, resulting in a sample size of 3440 (Burst type = 564, Linear type = 2875).

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

Models used in this study are described in S1 File.

Ethical approval

The study was conducted in compliance with the principles outlined in the Helsinki Declaration. Informed written consent was obtained from each subject, and the protocol was approved by the Institutional Review Board of each participating institution. The names of the ethics committee members and their reference numbers are listed in S3 Table.

Results

Among the 18,834 sites investigated in this study, 202 sites (1.1%) were improved, 18,272 sites (97.0%) were stable, and 360 (1.9%) sites were progressed during the 24-month follow-up period. At the patient level, 40 (32.3%) patients had only stable sites, 21 (16.9%) had stable and improved sites, 30 (24.2%) had stable and progressed sites, and 33 (26.7%) had stable, improved and progressed sites. At the tooth level, most of the teeth (89.2%) had only stable sites, and only 0.4% of the teeth had both improved and progressed sites. At the site level, the proportions of progressed sites were higher in the maxillary molars than in the other teeth. The descriptive statistics of the differences in CAL between the baseline and after 24 months (ΔCAL) are detailed in S4 Table.

By applying mixed effects modeling, we first analyzed ΔCAL as a dependent variable. When PlI and BOP were included in the fixed effects model, they were not statistically significant and were thus removed from the random effects model. The fixed effects model with all variables at baseline is shown in S5 Table. Table 1 demonstrates the explanatory variables that were significantly associated with ΔCAL in a random intercept model.

thumbnail
Table 1. Multilevel random intercept model for changes in CAL between the baseline and after 24 months (Model 1).

https://doi.org/10.1371/journal.pone.0188670.t001

Both periodontal pathogens were statistically significant. The coefficients of CAL at the baseline were all negative. Almost all coefficients of the tooth surfaces were statistically significant when the lingual surface of the mandibular anterior tooth was used as a reference. The coefficient of the buccal surface of the maxillary molar teeth was the highest. The value of ΔCAL appeared to vary between the tooth surfaces, and the CAL value at baseline distinctively affected the ΔCAL value. The mean predictive values were largely consistent with the observed values (S6 Table).

For the variables of CAL at baseline and tooth type, random slopes were included in the models. Results are shown in S7 Table. The random slope of CAL at baseline was statistically significant. All tooth types were statistically significant when the mandibular anterior tooth was used as a reference. These results indicated that ΔCAL varies according to the CAL at baseline, and the random slope for molars was different from those of premolars and anterior teeth.

The CAL changes are illustrated separately by combinations of tooth types and the CAL at baseline (Fig 1). A baseline CAL of < 3 mm showed gradual progression, while a baseline CAL of > 3 mm showed improvement. Molars with a baseline CAL of 3 mm were progressed, whereas premolars and anterior teeth were stable or improved. Therefore, because the CAL behavior depended on the CAL at the baseline, mixed effects models were constructed individually using the baseline CAL.

thumbnail
Fig 1. Mean values of the CAL changes during the 24-month follow-up period.

CAL changes during the 24-month follow-up period are separately illustrated by the CAL at baseline and by the type of tooth surface. Baseline CAL values are divided into three groups: (A) <3mm; (B) 3 mm; and (C) > 3 mm.

─●─: Maxillary molar, ---■---: Maxillary premolar, ···▲···: Maxillary anterior,

─○─: Maxillary molar, ---□---: Maxillary premolar, ···△···: Maxillary anterior

Baseline CAL values of < 3mm gradually deteriorated, while baseline CAL values of > 3 mm improved. Molars with a baseline CAL of 3 mm progressed, whereas premolars and anterior teeth were stable or improved.

CAL: clinical attachment level.

https://doi.org/10.1371/journal.pone.0188670.g001

The CAL changes were analyzed by the mixed effects model with repeated measures. The tooth surface was used as explanatory variables. For the groups of CAL at baseline < 3 mm and > 3 mm, random intercept was included for the tooth level. For the group of CAL at baseline = 3 mm, random slope was included for the tooth type. The results are listed in Table 2.

thumbnail
Table 2. Multilevel random effects model with repeated measures for (Model 2 (A), (B) and (C)).

https://doi.org/10.1371/journal.pone.0188670.t002

For groups with a baseline CAL > 3 mm, the coefficient of time was negative, which indicated that sites with CAL > 3 mm at baseline were generally improved. For the group with baseline CAL of 3 mm (Model 2 (B)), the coefficients of mandibular molar, maxillary premolar and maxillary molar were statistically significant. This indicated that for these teeth, the slope of their CAL changes across 24 months were different from that of mandibular anterior teeth. In addition, mandibular and maxillary molars had the highest coefficients of all teeth. These results are consistent with Fig 1. The fixed effects models for all explanatory variables used in this study are shown in S8 Table. When tooth type was combined with tooth surface and analyzed similarly to Model 1 (Table 1), the maxillary molar coefficient was higher than those of other tooth surfaces in the random effects models. These results are shown in S9 Table.

Finally, we classified the CAL changes during the 24-month follow-up period into six groups: improved, slightly improved, stable, slightly progressed, progressed and fluctuated, which accounted for 422 (2.2%), 1254 (6.7%), 14866 (78.9%), 1523 (8.1%) 688 (3.7%) and 81 (0.4%) changes, respectively. For the progressed group, a cluster analysis was performed that generated 5 clusters, which accounted for 575 (86.3%), 65 (9.4%), 9 (1.3%), 37 (5.4%) and 2 (0.3%) changes. The changes of these six groups and 5 clusters are portrayed in Fig 2(A) and 2(B). The slope of cluster 1 was moderate, and the slopes of the other clusters were steep. Cluster 1 may correspond to linear-type progressed sites, and the other clusters may correspond to burst-type progressed sites. Cross-tabulations by the number of linear- and burst-type progressed sites at the patient level and the tooth level are shown in S10 Table. Among the 83 patients with at least one progressed site, 43 (51.8%) had only a linear-type progressed site and 40 (51.8%) had both linear- and burst-type progressed sites. One patient had a burst-type progressed site without a linear type. Among the 71 teeth with a burst-type progression site, 32 (45.1%) lacked a linear progression site and 39 (54.9%) had an accompanying linear-type progression site. To investigate the predictors for distinguishing between linear-type and burst-type progression, a multilevel logistic regression analysis with repeated measures was performed. The fixed effects model to distinguish burst-type from linear-type progression using all the independent variables used in this study is shown in S11 Table. Table 3 shows the results of the random intercept model. The buccal surface of maxillary premolars and maxillary anterior teeth were statistically significant variables, and their coefficients were positive. The coefficients of tooth mobility were statistically significant, and their coefficients were negative.

thumbnail
Fig 2. CAL change patterns during the 24-month follow-up period.

  1. (A) Changes of the improved, slightly improved, stable, slightly progressed, progressed and fluctuated categories.
    ···▲···: Improved, ··△···: Slightly improved, ─●─: Stable.
    ---□---: Slightly progressed, ---■---: progressed, ─■─: Fluctuated
    Differences in the CAL changes over 24 months were classified into six categories: ≤ -3 mm, improved; between -3 mm and -2 mm, slightly improved; between -1 mm to 1 mm, stable; between 1 mm and 2 mm, slightly progressed; 3mm, progressed. In addition, cases with both ≤ -3 mm and ≥ 3mm were classified as fluctuated.
  2. (B) CAL progression patterns of the progressed category
    ···▲···: Cluster 1, ···△···: Cluster 2, ─●─: Cluster 3.
    ─□─: Cluster 4, —■—: Cluster 5
    A hierarchical cluster analysis was performed for the progressed type portrayed in Fig. 2(A), and 5 clusters were generated. The slope of cluster 1 was moderate, and the slopes of the other clusters were steep. Cluster 1 may correspond to the linear-type progressed sites, and the other clusters may correspond to the burst-type progressed sites.

https://doi.org/10.1371/journal.pone.0188670.g002

thumbnail
Table 3. Multilevel logistic regression model with repeated measures to distinguish “linear” and “burst” progression during the 24-month follow-up period.

https://doi.org/10.1371/journal.pone.0188670.t003

Discussion

The characteristics of periodontal disease represent heterogeneous phenomena, and the clinical manifestations of the disease caused by interactions between the host and agents may vary according to the site and tooth and may even vary within the same individual and longitudinally over time [33, 34]. Studies of the associations between clinical or etiological parameters and the progression of periodontal disease during follow-up have been conducted at the patient level, tooth level and site level. At the patient level, aggregated summary measurements of site-level information have been utilized and were found to result in a loss of information and subsequent misinterpretation. Furthermore, standard procedures for examining periodontal tissue rely on a site-based examination, and the management of a periodontal lesion is predominantly site specific [35]. Another interest of dental clinicians is the level of disease and its behavior at each site because of the site-specific nature of periodontal disease [3640]. Most sites are inactive and do not show progression. In fact, the number of stable sites in this study was 18,272 (97.0%).

As shown in S3 Table, even using a primitive descriptive analysis and cross-tabulation, a site-level analysis provides beneficial information. At the patient level, 32.3% of the patients had only stable sites, and 26.7% of the patients had both improved and progressed sites during the 24-month follow-up period. At the tooth level, most of the teeth were stable, and an extreme minority (0.4%) had both improved and progressed sites. Progressed and improved teeth were observed in only 6.4% and 3.4% of the patients, respectively. These data indicate that in approximately one-quarter of the patients, progression and improvement proceed simultaneously in the same individual. Therefore, the classification of patients as progressive, stable, or improved may be elusive. At the tooth level, most teeth were stable, and progressed teeth were a small minority. However, in dental practice, identifying these minority downhill teeth or sites is crucial for predicting the disease progression using clinical parameters and efficiently managing the etiological agents.

Multilevel modeling is one of the techniques of mixed effects analyses and has been employed for periodontal data. Previous studies have specified the risk factors that are associated with the progression of periodontal conditions at the patient, tooth and site levels. Novel findings were observed with respect to site-level factors. The coefficients of most tooth surfaces and the CAL at baseline were statistically significant for the ΔCAL value. As shown in Fig 1, for the site with a baseline CAL of 1–2 mm, the CALs were progressed, and the level of progression was highly dependent on the type of tooth, with coefficients varying from 0.06 to 0.186. All tooth type coefficients were statistically significant. Therefore, the sites with a baseline CAL of < 3 mm were generally progressed, and their changes varied among tooth types. For sites with a baseline CAL of 3 mm, the coefficients of the maxillary and mandibular molars were statistically significant and positive. Molars with a baseline CAL of 3 mm tended to be progressed. For the sites with a baseline CAL of > 3 mm, the coefficients of the maxillary molars were statistically significant and positive. Therefore, maxillary molars with a baseline CAL > 3 mm tended to be progressed. These findings were consistent with those of previous reports that have applied multilevel modeling. The coefficients for the tooth type increased from anterior to posterior teeth for the CAL progression [19] and the PPD [17], which was likely because of the inaccessibility of plaque control based on the specific anatomy of multi-rooted teeth. However, the cited studies did not classify their results by maxilla, mandible, and tooth surface. The finding that the CAL of the buccal surface of maxillary molars tended to progress is a new finding. The results from an observational study that was conducted during supportive periodontal therapy (SPT) and did not use a multilevel analysis showed that most of the lost teeth were maxillary second molars, followed by maxillary first molars, and these results reinforce our findings [11].

Few previous studies have been focused on the diagnostic value of tooth mobility, which is an interval-scale parameter; thus, individual mean values of tooth mobility are not available. Therefore, at the patient level, tooth mobility cannot be used as a dependent and/or independent variable with traditional statistical models. A mixed effects model, which can analyze hierarchical data simultaneously, can use tooth mobility as an independent variable. A previous study showed that tooth mobility is one determinant of PPD [17]. Our results showed that tooth mobility is a statistically significant factor for the ΔCAL and CAL changes. Additionally, tooth mobility was a significant predictor of burst-type progression.

Bleeding on probing (BOP) is a diagnostic criterion of gingival inflammation. Its utility as a clinical indicator for periodontal disease progression has been investigated. At the patient level, patients with a mean BOP of ≤ 20% have a significantly lower risk of CAL progression [41]. Patients with 16% or more BOP sites have a higher chance of losing attachment [42]. On the site level, periodontal pockets with successive incidence of BOP have a chance of losing attachment [42]. However, it has been shown that the positive predictive value of BOP for attachment loss is only 6%, while the negative predictive value is 98% [43]. In the present study, BOP was not a statistically significant predictor of ΔCAL. An absence of BOP may be an indicator of periodontal stability. However, further study is necessary to determine if the presence of BOP is an indicator of periodontal disease progression.

Considerable evidence has been obtained showing that periodontal pathogens are important indicators or predictors of the prognosis of periodontal disease. Several previous studies that employed mixed effects modeling have shown that among hierarchical factors, the contribution of patient-level factors was small compared with that of tooth- or site-level factors. The salivary levels of P. gingivalis and A. actinomycetemcomitans were statistically significant determinants of CAL progression, and their coefficient was not smaller than that of tooth- and site-level factors. However, for the CAL changes, the coefficient of A. actinomycetemcomitans was higher than that of the tooth- and site-level factors. Our previous report showed that the detection rate of A. actinomycetemcomitans was low. Therefore, many patients were not infected; however, carriers of A. actinomycetemcomitans were at a greater risk of CAL progression.

Two competing theories can explain the progression of periodontal disease. One theory is presented by the linear model, in which the overall sites slowly and progressively lose attachment. This theory has been supported by epidemiological studies [44, 45]. Another theory is presented by the burst model, in which multiple random sites show progression within a short period along with remission or repair. Other sites are free of progression throughout the individual’s life [32, 4648].

Several models have been proposed to explain these dynamic conditions of periodontal disease. The “linear” model slightly overestimates changes, and the “burst” model slightly underestimates changes [49, 50]. The lack-of-fit methods could not explain the dynamic conditions [51]. Certain sites improved while others deteriorated in a cyclical manner; thus, the linear and burst theories of periodontal disease progression are a manifestation of the same phenomenon [19]. The cited studies have attempted to resolve the dynamic conditions of all sites via single equations. Therefore, these studies did not obtain conclusive results and do not provide useful information for clinicians. As described in S10 Table, both linear- and burst-type progression occurred simultaneously within the same patients, and more than half of the teeth presented burst-type progression sites accompanied by linear-type progression sites. Thus, appropriate clustering is indispensable. The results showing that the palatal surface of the maxillary premolar and the buccal surface of maxillary anterior teeth tended to exhibit burst-type progression may be useful information for dental clinicians.

This study assessed various predictors of CAL changes. In summary, the important findings were as follows:

  • The baseline CAL was an important determinant of the direction of CAL changes, and the CAL changes varied widely in accordance with the tooth surface;
  • Salivary levels of periodontal pathogens, such as P. gingivalis and A. actinomycetemcomitans, were affected by CAL progression;
  • Linear- and burst-type progression of the CAL occurred simultaneously within the same patient, and more than half of the teeth with a burst-type progression site had an accompanying linear-type progression site; furthermore, most of the progressions were of the linear type, while the palatal surface of the maxillary premolar and buccal surface of maxillary anterior teeth tended to show burst-type progression.

Conclusions

The parameters identified in this study may be useful to practitioners for deciding the type and extent of treatment needed at both the site and patient levels. In addition, prior hypotheses concerning "burst" and "linear" theories are not valid.

Supporting information

S1 Table. All the data analyzed in this study.

https://doi.org/10.1371/journal.pone.0188670.s002

(XLSX)

S2 Table. Descriptive analyses of the variables analyzed in this study.

https://doi.org/10.1371/journal.pone.0188670.s003

(DOCX)

S3 Table. The name of the ethics committees, the committee’s reference number.

https://doi.org/10.1371/journal.pone.0188670.s004

(DOCX)

S4 Table. Descriptive analysis of the CAL change (ΔCAL) by patient, tooth and site level.

https://doi.org/10.1371/journal.pone.0188670.s005

(DOCX)

S5 Table. Fixed effect model for ΔCAL by all the variables.

https://doi.org/10.1371/journal.pone.0188670.s006

(DOCX)

S6 Table. Observed and predictive values of the changes in CAL between the baseline and 24 months.

https://doi.org/10.1371/journal.pone.0188670.s007

(DOCX)

S7 Table. Models for ΔCAL by random slope of tooth type and CAL at baseline.

https://doi.org/10.1371/journal.pone.0188670.s008

(DOCX)

S8 Table. Fixed effect model with repeated measures for the CAL changes during the 24-month follow up periods.

https://doi.org/10.1371/journal.pone.0188670.s009

(DOCX)

S9 Table. Multilevel random effect model with repeated measures for the CAL changes during the 24-month follow up periods.

https://doi.org/10.1371/journal.pone.0188670.s010

(DOCX)

S10 Table. Cross-tabulation of the number of “linear”- and “burst”-type progressed sites in 83 patients and 395 teeth with at least one progression site.

https://doi.org/10.1371/journal.pone.0188670.s011

(DOCX)

S11 Table. Fixed effect model with repeated measures for the CAL changes during the 24-month follow up periods.

https://doi.org/10.1371/journal.pone.0188670.s012

(DOCX)

References

  1. 1. Axtelius B, So derfeldt B, Edwardsson S, Attström R. Therapy resistant periodontitis (II). Compliance and general and dental health experiences. J Clin Periodontol. 1997; 24: 646–653. pmid:9378837
  2. 2. Buchmann R, Muller RF, Heineke A, Lange DE. Actinobacillus actinomycetemcomitans in destructive periodontal disease. Three-year follow-up results. J Periodontol. 2000; 71, 444–453. pmid:10776933
  3. 3. De Sanctis M, Zucchelli G. Interleukin-1 gene polymorphisms and long-term stability following guided tissue regeneration therapy. J Periodontol. 2000; 71, 606–613. pmid:10807125
  4. 4. König J, Plagmann HC, Langenfeld N, Kocher T. Retrospective comparison of clinical variables between compliant and non-compliant patients. J Clin Periodontol. 2001; 28, 227–232. pmid:11284535
  5. 5. Lang NP, Tonetti MS. Periodontal risk assessment (PRA) for patients in supportive periodontal therapy (SPT). Oral Health Prev Dent. 2003; 1, 7–16. pmid:15643744
  6. 6. Serino G, Rosling B, Ramberg P, Hellstro MK, Socransky SS, Lindhe J. The long-term effect of systemic tetracycline used as an adjunct to non-surgical treatment of advanced periodontitis. J Clin Periodontol. 2001; 28, 446–452. pmid:11350508
  7. 7. Cattabriga M, Pedrazzoli V, Cattabriga A, Pannutti E, Trapani M, Verrocchi GC. Tetracycline fiber used alone or with scaling and root planing in periodontal maintenance patients: clinical results. Quintessence Int. 1996; 27, 395–400. pmid:8941833
  8. 8. Cattabriga M, Rotundo R, Muzzi L, Nieri M, Verrocchi G, Cario F et al.Retrospective evaluation of the influence of the interleukin-1 genotype on radiographic bone levels in treated periodontal patients over 10 years. J Periodontol. 2001; 72, 767–773. pmid:11453239
  9. 9. Eickholz P, Kaltschmitt J, Berbig J, Reitmeir P, Pretzl B. Tooth loss after active periodontal therapy. 1: patient-related factors for risk, prognosis, and quality of outcome. J Clin Periodontol. 2008; 35, 165–174. pmid:18199150
  10. 10. Armitage G. Manual periodontal probing in supportive periodontal treatment. Periodontology 2000. 1996; 12, 33–39. pmid:9567991
  11. 11. König J, Plagmann HC, Rühling A, Kocher T. Tooth loss and pocket probing depths in compliant periodontally treated patients: a retrospective analysis. J Clin Periodontol. 2002; 29, 1092–1100. pmid:12492910
  12. 12. Socransky SS, Haffajee AD. The nature of periodontal diseases. Ann Periodontol. 1997 2, 3–10. pmid:9151538
  13. 13. Cho Y, Kim HY. Analysis of periodontal data using mixed effects models. J. Periodontal Implant Sci. 2015; 45, 2–7. pmid:25722920
  14. 14. Sterne JA, Curtis MA, Gillett IR, Griffiths GS, Maiden MF, Wilton JM, Johnson NW. Statistical models for data from periodontal research. J Clin Periodontol. 1990; 17, 129–137. pmid:2180989
  15. 15. Albandar JM, Goldstein H. Multi-level statistical models in studies of periodontal diseases. J Periodontol. 1992; 63, 690–695. pmid:1507050
  16. 16. Axtelius B, Söderfeldt B, Attström R. A multilevel analysis of factors affecting pocket probing depth in patients responding differently to periodontal treatment. J Clin Periodontol. 1999; 26, 67–76. pmid:10048639
  17. 17. Tu YK, Gilthorpe MS, Griffiths GS, Maddick IH, Eaton KA, Johnson NW. The application of multilevel modeling in the analysis of longitudinal periodontal data—part I: absolute levels of disease. J Periodontol. 2004;75, 127–136. pmid:15025224
  18. 18. Tu YK, Gilthorpe MS, Griffiths GS, Maddick IH, Eaton KA, Johnson NW. The application of multilevel modeling in the analysis of longitudinal periodontal data—part II: changes in disease levels over time. J Periodontol. 2004; 75, 137–145. pmid:15025225
  19. 19. Gilthorpe MS, Zamzuri AT, Griffiths GS, Maddick IH, Eaton KA, Johnson NW. Unification of the "burst" and "linear" theories of periodontal disease progression: a multilevel manifestation of the same phenomenon. J Dent Res. 2003;82, 200–205. pmid:12598549
  20. 20. Nieri M, Muzzi L, Cattabriga M, Rotundo R, Cairo F, Pini Prato GP. The prognostic value of several periodontal factors measured as radiographic bone level variation: a 10-year retrospective multilevel analysis of treated and maintained periodontal patients. J Periodontol. 2002; 73, 1485–1493. pmid:12546099
  21. 21. Krebs KA, Clem DS 3rd. American Academy of Periodontology, Guidelines for the management of patients with periodontal diseases. J Periodontol 2006;77:1607–1611.
  22. 22. Nomura Y, Shimada Y, Hanada N, Numabe Y, Kamoi K, Sato T et al. Salivary biomarkers for predicting the progression of chronic periodontitis. Arch Oral Biol 2012;57:413–420. pmid:22030151
  23. 23. Morozumi T, Nakagawa T, Nomura Y, Sugaya T, Kawanami M, Suzuki F et al. Salivary pathogen and serum antibody to assess the progression of chronic periodontitis: a 24-mo prospective multicenter cohort study. J Periodontal Res. 2016; 51, 768–778. pmid:26791469
  24. 24. Kakuta E, Nomura Y, Morozumi T, Nakagawa T, Nakamura T, Noguchi K, et al. Assessing the progression of chronic periodontitis using subgingival pathogen levels: a 24-month prospective multicenter cohort study. BMC Oral Health. 2017; 17: 46. pmid:28093069
  25. 25. Tadokoro K, Yamaguchi T, Kawamura K, Shimizu H, Egashira T, Minabe M, et al. Rapid quantification of periodontitis-related bacteria using a novel modification of Invader PLUS technologies. Microbiol Res 2010; 165, 43–49. pmid:18718748
  26. 26. Tada A, Takeuchi H, Shimizu H, Tadokoro K, Tanaka K, Kawamura K, et al. Quantification of periodontopathic bacteria in saliva using the invader assay. Jpn J Infect Dis 2012; 65, 415–423. pmid:22996215. 22996215
  27. 27. Beck JD. Methods of assessing risk for periodontitis and developing multifactorial models. J Periodontol. 1994;65:468–78. pmid:8046563
  28. 28. Ogawa H, Yoshihara A, Hirotomi T, Ando Y, Miyazaki H. Risk factors for periodontal disease progression among elderly people. J Clin Periodontol 2002; 29: 592–597. pmid:12354083
  29. 29. Kudo C, Naruishi K, Maeda H, Abiko Y, Hino T, Iwata M, et al. Assessment of the plasma/serum IgG test to screen for periodontitis. J Dent Res 2012; 91: 1190–1195. pmid:23018816
  30. 30. Muhlemann HR. Periodontometry, a method for measuring tooth mobility. Oral Surg Oral Med Oral Pathol. 1951; 10, 1220–1233.
  31. 31. Lang Niklaus P., Lindhe Jan. Clinical Periodontology and Implant Dentistry. 6th ed. Wiley-Blackwell: 2015.
  32. 32. Socransky SS, Haffajee AD, Goodson JM, Lindhe J. New concepts of destructive periodontal disease. J Clin Periodontol. 1984; 11, 21–32. pmid:6582072
  33. 33. Tonetti MS, Claffey N; European Workshop in Periodontology group C. Advances in the progression of periodontitis and proposal of definitions of a periodontitis case and disease progression for use in risk factor research. Group C consensus report of the 5th European Workshop in Periodontology. J Clin Periodontol. 2005; 32 Suppl 6:210–213. pmid:16128839
  34. 34. Persson GR. Site-based versus subject-based periodontal diagnosis. Periodontol 2000. 2005; 39,145–163. pmid:16135069
  35. 35. Mancl LA, Leroux BG, DeRouen TA. Between-subject and within-subject statistical information in dental research. J Dent Res. 2000; 79,1778–1781. pmid:11077994
  36. 36. Lindhe J, Haffajee AD, Socransky SS. Progression of periodontal disease in adult subjects in the absence of periodontal therapy. J Clin Periodontol. 1983;10,433–442. pmid:6577035
  37. 37. Haffajee AD, Socransky SS, Ebersole JL. Survival analysis of periodontal sites before and after periodontal therapy. J Clin Periodontol. 1985;12:553–567. pmid:3860518
  38. 38. Haffajee AD, Socransky SS. Attachment level changes in destructive periodontal diseases. J Clin Periodontol. 1986; 13, 461–475. pmid:3522651
  39. 39. Osborn J. The choice of computational unit in the statistical analysis of unbalanced clinical trials. J Clin Periodontol. 1987; 14, 519–523. pmid:3316295
  40. 40. Claffey N, Egelberg J. Clinical indicators of probing attachment loss following initial periodontal treatment in advanced periodontitis patients. J Clin Periodontol. 1995; 22, 690–696. pmid:7593699
  41. 41. Joss A, Adler R, Lang NP. Bleeding on probing. A parameter for monitoring periodontal conditions in clinical practice. J Clin Periodontol. 1994; 21, 402–408. pmid:8089242
  42. 42. Lang NP, Joss A, Orsanic T, Gusberti FA, Siegrist BE. Bleeding on probing. A predictor for the progression of periodontal disease? J Clin Periodontol. 1986; 13, 590–596. pmid:3489010
  43. 43. Lang NP, Adler R, Joss A, Nyman S. Absence of bleeding on probing. An indicator of periodontal stability. J Clin Periodontol. 1990; 17, 714–721. pmid:2262585
  44. 44. Axtelius B, Söderfeldt B, Attström R. A multilevel analysis of factors affecting pocket probing depth in patients responding differently to periodontal treatment. J Clin Periodontol. 1999; 26, 67–76. pmid:10048639
  45. 45. Löe H, Anerud A, Boysen H. The natural history of periodontal disease in man: prevalence, severity, and extent of gingival recession. J Periodontol 1992; 63, 489–495. pmid:1625148
  46. 46. Löe H, Anerud A, Boysen H, Smith M. The natural history of periodontal disease in man. The rate of periodontal destruction before 40 years of age. J Periodontol 1978; 49:607–620. pmid:282430
  47. 47. Cohen ME. Bursts of periodontal destruction and remission, percolation phase shifts, and chaos. J Periodontal Res. 1993; 28, 429–436. pmid:8254460
  48. 48. Cohen ME, Ralls SA. Distributions of periodontal attachment levels. Mathematical models and implications. J Periodontol. 1988; 59, 254–258. pmid:3164376
  49. 49. Ralls SA, Cohen ME. Problems in identifying "bursts" of periodontal attachment loss. J Periodontol. 1986; 57, 746–752. pmid:3467060
  50. 50. Allen KM, Hausmann E. Model fit and measurement outcome in attachment measurements: a simulation study. J Periodontal Res. 1995; 30, 15–22. pmid:7722843
  51. 51. Hujoel PP, Leroux BG. Evaluating the burst hypothesis at a site-specific level using the lack-of-fit test. J Periodontol. 1998; 69, 357–362. pmid:9579622