Data sources
In this study, we used the data from the Longitudinal Health Insurance Database (LHID) in Taiwan (2000–2015), which is a subset database selected randomly from the National Health Insurance Research Database (NHIRD). The National Health Insurance (NHI) Program in Taiwan, which was launched in 1995, includes approximately 23 million beneficiaries or more than 99% of the entire population in Taiwan [17]. The NHI program covers all necessary medical care (including outpatient and inpatient), dental care, Chinese medicine, and prescription drugs. The LHID contains information on health service utilization for approximately one million beneficiaries who represent approximately 5% of the Taiwanese population. The NHIRD contains patient identification numbers, birthdays, sexes, ICD-9-CM diagnostic codes (up to five each), and outcomes.
Study design and sampled participants
Direction 1, that is, tracing T2DM in patients with periodontitis, comprised patients who were newly diagnosed with periodontitis from January 1, 2000, to December 31, 2015, according to the ICD-9-CM code 523.4. Each enrolled periodontitis patient was required to have made at least three dental visits with the 523.4 code being filed within the previous 1 year. The periodontitis patients were then categorized into two groups. The patients who had NHI order code for sub-gingival curettage/root planning (91006-91008C) or periodontal flap operation (91009B-91010B) were categorized into the treated periodontitis group (Pt group), whereas those without these codes were categorized into the periodontitis without therapy group (Pwt group)[18]. The exclusion criteria for the study were as follows: the patients with periodontitis from January 1, 1998, to December 31, 1999; the patients with T2DM from 1998 to 1999 or before their first visit during which periodontitis was diagnosed; the subjects aged < 40 years; and those who had insufficient medical information or withdrawal from the NHI program throughout the 15-year study period (Fig. 1). Among the 1,936,512 individuals, there were 11,011 periodontitis patients enrolled in the Pt group. After matching for age, sex, and index date, 11,011 periodontitis patients were included in the Pwt group (same exclusion criteria and one-fold propensity score) (Fig. 1). Additionally, another 11,011 dental patients who were not diagnosed with periodontitis were categorized into the non-periodontitis control group.
Direction 2, that is, tracking periodontitis in patients with T2DM who were taking ≥ 2 anti-diabetic medications, consisted of patients who were selected from the medical claim data according to the ICD-9-CM code 250. Patients who were diagnosed with T2DM or periodontitis from 1998 to 1999 and aged < 40 years were excluded. In total, 157,798 patients with T2DM were enrolled in the diabetes group, whereas 157,798 participants were recruited in the non-diabetes group. The outcome of periodontitis was then tracked. Besides, the outcome was subsequently categorized into the treated periodontitis subgroup (Pt) and the periodontitis without therapy subgroup (Pwt).
The covariates included gender, age, insurance premium (< 18,000, 18,000–34,999, and ≥ 35,000 NT$), urbanization level of residence (levels 1–4), and level of care (hospital center, regional, and local hospitals). The urbanization level of residence was based on the population and various other indicators. Briefly, the urbanization level 1 was defined as a population > 1,250,000 and specific designation, while the levels 2, 3, and 4 as populations between 500,000 and 1,249,999; 149,999 and 499,999; and < 149,999, respectively.
The baseline comorbidities included hypertension (ICD-9-CM codes: 401.1, 401.9, 402.10, 402.90, 404.10, 404.90, 405.1, and 405.9), hyperlipidemia (ICD-9 CM code: 272.x), coronary artery disease (CAD; ICD-9 CM code: 410–414), obesity (ICD-9 CM codes: 278.00–278.01), smoking (ICD-9-CM code: 305.1), chronic obstructive pulmonary disease (COPD; ICD-9-CM codes: 490–496), and alcoholism (ICD-9-CM codes: 303 and 305.0); additionally, the revised Charlson Comorbidity Index (CCI_R; CCI removed diabetes mellitus, hypertension, and CAD) was included.
Statistical analysis
Chi-square and Fisher’s exact tests were used to evaluate the differences between categorical variables, whereas the t-test and the one-way analysis of variance with Scheffe’s post hoc test were used for continuous variables. Multivariate Cox proportional hazards regression analysis was used to determine the risk of T2DM and periodontitis (directions 1 and 2, respectively). The results were presented as hazard ratios (HRs) with 95% confidence intervals (CIs). Sensitivity analysis was further used to exclude the diagnosis of dementia in the first 1 or 5 year(s) and to eliminate any potential protopathic bias. The difference in the risk of outcome disease between the study and reference groups was estimated using the Kaplan-Meier method and the log-rank test. A two-tailed p < 0.05 was considered statistically significant.