ABSTRACT
Aim Our goal was to investigate individual susceptibility to periodontitis by conducting an epigenome-wide association study using peripheral blood.
Materials and Methods For this analysis, we included 1077 African American and 457 European American participants of the Atherosclerosis Risk in Communities (ARIC) study who had completed a dental examination or reported being edentulous at visit 4 and had available data on DNA methylation. DNA methylation levels were compared by periodontal disease severity and edentulism to identify differentially methylated regions (DMRs) and evaluate the CpGs belonging to those DMRs using multinominal logistic regression.
Results We identified a region in gene ZFP57 (6p22.1) that was significantly hypomethylated in severe periodontal disease compared to no/mild periodontal disease in European American participants. A separate region in an unknown gene (located in Chr10: 743,992-744,958) demonstrated significant positive association with edentulism compared to no/mild periodontal disease in African American participants. Four CpGs in a region located within HOXA4 were significantly hypermethylated in severe periodontal disease compared to no/mild periodontal disease in African American participants.
Conclusions Our study highlights epigenetic variations in ZPF57 and HOXA4 that were significantly and reproducibly associated with periodontitis. Future studies should evaluate gene regulatory mechanisms in the candidate regions of these loci.
Scientific Rationale for Study Without altering the DNA sequence, epigenetic effects (e.g., DNA methylation changes) can alter gene activity and influence host response to periodontal infections. Our well-powered study investigates individual susceptibility to periodontitis by conducting a thorough assessment of periodontitis-related DNA methylation levels in blood.
Principal Findings We identified two gene regions, ZPF57 and HOXA4, that are differentially methylated in individuals with compared to those without periodontitis.
Practical implications Studying differential leukocyte DNA methylation patterns may point to candidate regions and underlying gene regulatory mechanisms that play a key role in the progression and/or susceptibility to periodontitis.
Introduction
Over 47% of US adults have some form of periodontitis, and the incidence of periodontitis is rising steeply in adults 50 and older.1-3 Given the public health impact of periodontal disease4,5, studying individuals’ susceptibility to periodontitis could be crucial to building effective disease prevention strategies. In the last two decades, genome-wide association studies have attempted to discover genes that contribute to the pathogenesis of periodontitis; but these studies have found no common causal allele that reached genome-wide significance.6,7 Epigenetic modifications of DNA can affect the genetic blueprint of host responses, being either directly lesion-associated in target tissue or associated with an altered immune response; thus, it is of particular interest to assess epigenetic alterations in whole blood as it may (directly or indirectly) modulate the development and progression of periodontitis.
In the context of periodontal disease, biological mechanisms leading to periodontitis are driven by infections,8 and the host response plays a critical role in individual response. Given that changes in methylation levels can reveal a particular host response, studies of DNA methylation may yield valuable information about disease susceptibility to periodontitis. Two epigenome-wide studies (EWAS) of peripheral blood leukocytes in periodontitis have been published. A cross-sectional study in adult female twins by Kurushima et al.9 identified several CpGs, including one in ZNF804A, associated with gingival bleeding and tooth mobility. A pilot study by Hernandez et al.10, conducted in 8 periodontitis patients and 8 matched controls, reported statistically significant differences in peripheral blood DNA methylation in multiple regions, including the ZNF718, HOXA4, and ZFP57 genes.
To expand on previous work of periodontitis-related DNA methylation changes, we conducted an EWAS to discover and validate aberrant patterns of blood-derived DNA methylation in periodontal disease severity and edentulism among 1534 Black and White older adults who underwent a clinical dental examination as a component of the community-based ARIC study. Our study aimed to identify region- and CpG-specific differential methylation in peripheral blood leukocytes collected from participants with moderate or severe periodontitis compared to those with no or no/mild periodontitis, controlling for relevant health conditions, socioeconomic status, self-reported smoking, and methylation-predicted packyears.
Material and Methods
Study Population
Study participants were enrolled in the Atherosclerosis Risk in Communities (ARIC) Study (RRID: SCR_021769), a prospective cohort study of cardiovascular disease risk that enrolled 15,792 participants between 1987 and 1989 from four US communities (Jackson, MS; Washington County, MD; suburban Minneapolis, MN; and Forsyth County, NC).11,12 Participants underwent a baseline clinical examination (Visit 1) and subsequent follow-up clinical exams (Visits 2–9 completed & Visit 10 planned). Participants provided written informed consent at each visit. The ARIC Study protocol was approved by institutional review boards at the four field centers and coordinating center.
At Visit 4 (1996–1998), participants who had at least one tooth or dental implant and were otherwise eligible and consented underwent a clinical dental examination by trained personnel (n = 6017). Edentulous participants were not eligible for the dental examination. Our analysis focused on a sample of ARIC participants with available data on leukocyte DNA methylation levels collected at Visit 2 (1990-92) or Visit 3 (1993-1995). Among the participants with methylation data, those who attended a dental examination13 (n = 1174) or who self-reported being edentulous (n = 360) during Visit 4 were eligible for our study.
DNA Methylation Assessing and Preprocessing
The procedures for DNA methylation assessment have been previously described.14 Briefly, genomic DNA was extracted from peripheral blood leukocyte samples using the Gentra Puregene Blood Kit (Qiagen) and bisulphite conversion of 1 ug genomic DNA was performed using the EZ-96 DNA Methylation Kit (Deep Well Format) (Zymo Research).14 Details are summarized in the Supplemental Methods File.
Periodontal Measures
Periodontal severity was determined based on standardized dental measurements. Probing depth and gingival recession were measured at six sites on all teeth present. Clinical attachment loss (AL) was calculated based on these two measures.15,16 We acknowledged the complexity of classifying periodontal disease by utilizing two different measures of periodontitis severity. The first definition, “periodontitis (ARIC),” was previously used in ARIC to categorize individuals based on the percentage of sites with AL>= 3 mm into no/mild, moderate, or severe periodontitis.17 The second measure, “periodontitis (CDC-AAP),” was developed for population-based surveillance of periodontitis by the US Centers for Disease Control and Prevention–American Academy of Periodontology (CDC-AAP). It includes AL and pocket depth (PD) measurements and the site location (interproximal) to classify participants into no, mild, moderate, or severe periodontitis.18 We used edentulism as a surrogate for the most severe type of periodontitis. ARIC participants who self-reported as both not having any natural teeth and not having any dental implants were classified as edentulous for this analysis.
Covariate Assessment
Modifiable risk factors associated with periodontitis and edentulism include: tobacco smoking19, nutrition20, psychological stress and depression21, obesity, diabetes22, and metabolic syndrome.23 Data on cigarette smoking status (current, former, never) and cigarette smoking dose was measured at each visit during follow-up. Diabetes at baseline was defined as fasting plasma glucose ≥126□mg/dl, non-fasting plasma glucose ≥200□mg/dl, receiving diabetes treatment, or self-report of a diabetes diagnosis by a physician. A lifetime socioeconomic status (SES) variable, calculated based on 12 factors associated with SES, was used to identify participants with high life-course SES.24 Leukocyte cell composition was measured using a DNA methylation deconvolution procedure (see Supplemental Methods), given its potential to confound analyses of blood-based DNA methylation data. Lastly, we used the R function EstDimRMT to determine the number of principal components with non-zero eigenvalues needed to correct batch effects in the current analysis.25,26
Statistical Analysis
We examined the association between CpG-specific DNA methylation and periodontal disease by conducting an epigenome-wide association analysis (EWAS) using multivariable logistic regression models. Such models were used to estimate odds ratios (OR) of severe vs. no/mild periodontitis (ARIC definition), moderate vs. no/mild periodontitis (ARIC definition), and edentulous vs. no/mild periodontitis (ARIC definition), or severe vs. no periodontitis (CDC definition), mild/moderate vs. no periodontitis (CDC definition), and edentulous vs. no periodontitis (CDC definition) per 1 SD increase in methylation level at single CpG sites. All models were adjusted for age, sex, five surrogate variables for batch effects, smoking category (never, former, current), packyears, BMI, and leukocyte cell composition. In addition, to control for potential residual confounding by smoking, we included a methylation-predicted packyears developed based on a meta-analysis of association results between DNA methylation and cigarette smoking in individuals from 16 cohorts.27 The Benjamini-Hochberg method was used to calculate the false discovery rate (FDR) for multiple comparison adjustments (FDR q-value < 0.05). Lastly, we conducted EWAS separately for African American and European American participants. Since blood samples from African American and European American participants were assayed in separate Illumina runs and at different times, we treated methylation data of African American and European American participants as two separate cohorts.
Using the DMRcate R/Bioconductor package, we conducted differentially methylated region (DMR) analyses to identify DMRs associated with periodontal disease status. A linear model that adjusted for the same covariates as the single CpG analyses was fit using the DMRcate function along with these parameters: lambda = 1,000 and kernel adjustment C = 2.28 Statistically significant DMRs were required to have a minimum of two statistically significant single CpGs and to meet the multiple testing adjustment criteria of FDR (q-value) <0.05 (calculated using the Benjamini-Hochberg method).29
CpGs belonging to the statistically significant DMRs were evaluated using a baseline category logit model using both the CDC-AAP and the ARIC periodontitis definitions. In a prior study, Hernandez et al. conducted a pilot EWAS in 8 periodontitis patients and 8 matched controls and reported statistically significant differences in peripheral blood DNA methylation in multiple regions, including the ZNF718, HOXA4, and ZFP57 genes.10 These three DMRs were the only regions that were found to be significant with multiple consecutive CpGs in the same direction after controlling for cell composition in the model.10 We compiled CpGs from the three DMRs identified by Hernandez et al. and tested them in the ARIC dataset. All models controlled for age, sex, SES, BMI, diabetes status, smoking category (never, former, current), packyears, methylation-predicted packyears smoked, leukocyte cell composition, ARIC field center, and five surrogate variables for batch effects. Analyses were conducted separately for African American and European American participants.
Results
Population Characteristics
Table 1 summarizes the characteristics of ARIC participants included in this study by periodontal disease status defined by the CDC-AAP and by the ARIC periodontitis definitions. The average time between blood drawn and Visit 4 was 6 years. Participants with severe periodontitis were more likely to be male, black, obese, less educated, current smoker, and have lower life-course SES (Table 1). Compared to those of European ancestry, African American participants in this study were more likely to be diabetic, have low SES, are less educated, have fewer teeth, and are more likely to have severe periodontitis (Supplemental Table 1). African American participants were younger, with a median age of 55, compared to 58 years old for European American participants. 28.6% of African American participants had severe periodontal disease, with 45.7% having no/mild periodontal disease, while 11.4% and 60.2% of European American participants had severe or no/mild periodontal disease, respectively.
EWAS & DMRcate Analysis
Using the ARIC periodontitis definition, the single CpG EWAS did not identify differentially methylated CpGs that were statistically significant after multiple comparison corrections (q < 0.05). After adjusting for methylation-predicted packyears, we identified one statistically significant DMR (ENSG00000231601 manual annotation Chr10: 743,992-744,958) in edentulism when compared to no/mild periodontal disease in African American participants and one statistically significant DMR (ZFP57 6p22.1) in severe periodontal disease compared to no/mild periodontal disease in European American participants (Supplemental Table 2). The EWAS and DMRcate results did not vary from those presented above when using the CDC-AAP periodontitis definition for analysis.
Zinc Finger Protein Gene 57 (ZFP57)
The statistically significant DMR located on Chr6 (ZFP57 gene), consisting of 17 CpGs, differed between European American participants with severe periodontal disease and those with no/mild periodontal disease. The 17 CpGs perfectly overlapped with a 21-CpG ZFP57 region which was reported by Hernandez et al. as hypomethylated when comparing periodontitis patients with age-matched and sex-matched periodontally healthy controls living in the Republic of Colombia.10 Overlapping CpGs located in the ZFP57 DMR region were found to be associated with a lower risk of severe periodontal disease compared to no/mild periodontal disease (ARIC definition) in European American participants (Table 2). This association was stronger in European American participants 58 years or older (Table 3). It was also comparable for male and female participants and participants of higher or lower SES (Supplemental Table 3). In addition, there was a statistically significant decreased risk of moderate vs. no/mild periodontal disease among European American participants 58 years or older (Table 3). The ZFP57 CpGs were not significantly associated with risk of moderate or severe periodontal disease, as compared to no/mild periodontal disease in African American participants. However, increasing methylation of three ZFP57 CpGs (cg11383134, cg15570656, and cg13835168) was associated with a significant decreased risk of severe vs. no/mild periodontal disease in African American participants 58 years or older (Table 3). When we used the CDC-AAP definition to categorize periodontal disease status, increasing methylation of the ZFP57 CpGs was found to be associated with a significant decreased risk of mild/moderate vs. no periodontal disease in European American ARIC participants. The pattern of associations for those 58 years and older, male and female participants, participants of higher or lower SES, and for edentulism were all similar to what we found using the ARIC definition (Table 2 and Supplemental Table 3).
HOXA4
A HOXA4 DMR, consisting of 13 CpGs, was found to be hypermethylated in periodontitis patients with age-matched and sex-matched periodontally healthy controls.10 This finding is consistent with Hernandez et al. While the HOXA4 DMR was not statistically significant in the DMRcate analyses using both of our periodontitis definitions, when we tested the CpGs contained in the HOXA4 DMR individually, increasing methylation of cg11015251, cg14359292, cg07317062, and cg08657492 was statistically significantly associated with a higher risk of severe periodontal disease compared to no/mild periodontal disease in African American participants (Table 4). This association was comparable for African American participants 58 years or older and younger than 58-years-old, male and female participants, and participants of higher or lower life-course SES (Supplemental Table 4). In contrast, these HOXA4 CpGs were not significantly associated with risk of moderate or severe periodontal disease, as compared to no/mild periodontal disease, in European American participants (Table 4). In addition, increasing methylation of cg14359292 and cg24169822 was associated with decreased risk of edentulism vs. no periodontal disease among European American participants (Table 4; stronger association in females Supplemental Table 4). Using the CDC-AAP definition, we also found significant positive associations between increasing methylation of selected HOXA4 CpGs and periodontitis that were largely consistent with what we found using the ARIC definition.
ENSG00000231601
Using the ARIC periodontitis definition, our DMRcate analysis identified one statistically significant DMR located on Chr10 (ENSG0000023160 gene), consisting of three CpGs (cg22954052, cg25580864, and cg05495790), that differed between African American participants with edentulism and those with no/mild periodontal disease (Table 5). Among African American participants, increasing methylation of cg22954052 and cg05495790 was significantly associated with increased risks of moderate periodontal disease, severe periodontal disease, and edentulism, when compared to no/mild periodontal disease. These associations were stronger among African American participants under age 58 (Supplemental Table 5). cg22954052 and cg05495790 were associated with an increased risk of edentulism vs. no periodontal disease in African American participants. These CpGs were not associated with PD or edentulism in European American participants. When we used the CDC-AAP definition, increasing methylation of cg22954052 and cg05495790 was also significantly associated with increased risks of mild/moderate vs. no periodontal disease and edentulism vs. no periodontal disease, but not severe PD vs no PD, in African American participants.
Discussion
To our knowledge, this is the first large-scale epigenome-wide DNA methylation study of clinically assessed periodontitis. Using epigenome-wide analyses in 1534 ARIC participants, we identified differential methylation patterns comparing various clinically assessed periodontal disease classifications, suggesting that epigenetic changes may contribute to variation in susceptibility to periodontitis. We replicated findings from a previous study reporting an association between DNA methylation levels measured in whole blood at CpG sites in the Zinc Finger Protein Gene 57 (ZFP57) and periodontitis using two clinically assessed periodontal disease definitions. We also found DNA methylation differences near ENSG00000231601 (Chr10) between African American participants with edentulism and those with no or mild periodontal disease.
Our analysis used two periodontitis classifications, ARIC definition and CDC definition, to better handle the complexity of measuring this disease. Since periodontitis is a chronic disease, patients present a wide range of severity of symptoms and progression. Results varied slightly by periodontitis definition used; findings were stronger with the ARIC definition for each of the regions identified. The ARIC definition for “severe” is based on having extensive periodontal disease (30% of teeth with attachment loss >3 mm). In contrast, the CDC definition for “severe” requires a greater degree of attachment loss (>6 mm) but fewer teeth. Thus, the different results suggest that the genetic regions may be more strongly linked to extensive periodontal disease. Overall, we found comparable results for male and female participants and participants of higher or lower SES, but found greater associations among participants 58 or older, possibly due to a higher prevalence of extensive periodontal disease among older ARIC participants. The differences we observed between African Americans and European Americans may have been due to differences in the two populations. The large majority of African American participants are from the Jackson Field Center. Participants from this region experienced very different exposures during their lifetimes and have limited access to health and dental care (compared to other ARIC participants). These factors, which can influence periodontal disease development and progression, may have impacted findings in a manner similar to differences we observed with different disease classifications. Consequently, it is not possible to determine whether our results are due to genetic effects or correlates of socioeconomic factors.
Prior to the present study, there were only two published studies of genome-wide DNA methylation of peripheral blood leukocytes in periodontitis. Kurushima et al. (2019) conducted an epigenome-wide study using Illumina DNA methylation 450K data and a retrospective dental phenotype collection of self-reported dental mobility and gingival bleeding information obtained from 528 older female twins in the Twins UK cohort (mean age: 58 years old, age range: 19-82 years old; compared to African American ARIC participants mean age: 55 years old, age range: 47-70 years old, European American ARIC participants mean age: 58 years old, age range: 47-68 years old).9 This study detected a hypomethylated region in ZNF804A in individuals who experience gingival bleeding, in conjunction with an increased level of ornithine, a metabolite related to gingival inflammation, in blood metabolomics profile. There were no concordances between our results and the results of Kurushima et al., possibly due to clear differences in study design, including the use of clinically examined vs. self-reported periodontitis data. Hernandez et al. (2021) performed an analysis to identify differentially methylated positions and regions with a small number of clinically diagnosed periodontitis patients and gingivally healthy controls10. Results of Hernandez et al. identified ZNF718, HOXA4, and ZFP57 genes related to systemic immune-related epigenetic patterns in periodontitis10, including a ZFP57 region of 21 CpGs that was found to be hypomethylated when comparing periodontitis patients with age-matched and sex-matched periodontally healthy controls. While Hernandez et al. also carried out clinical examination, its study participants were much younger than ours (age range: 25-55 years old), suggesting a more extreme comparison of phenotype in Hernandez et al. Consistent with Hernandez et al., our study identified a ZFP57 DMR to be periodontitis-related. In addition, the 17 CpGs of this selected DMR perfectly overlapped with the 21-CpG ZFP57 DMR identified by Hernandez and colleagues.
Existing evidence supports the relevance of some of the identified ZFP57 CpGs for periodontitis. ZFP57 is an imprint regulator, and DNA methylation in this exact ZFP57 region has been shown to be environmentally modulated in human blood and related to transient neonatal diabetes, Parkinson’s disease, and post-traumatic stress disorder.30-32 The Zinc Finger Proteins (ZNFs) are also important transcription factors. They are involved in tissue development, and their alterations may promote chronic conditions, including oral cancer and periodontitis.33 Recent studies have shown ZNFs playing a role in the progression of periodontitis.34 For instance, ZNFA20 can inhibit NF-κβ pathway activation via inflammatory cytokines34,35 and thus support the disruption of local inflammatory responses by P. gingivalis, the keystone pathogen of periodontitis.36,37 The replication of the ZNF57 DMR in our study suggests a potential sequence of events in periodontitis development, where hypomethylation of the ZNF57 DMR may upregulate the expression of ZNF57, which is related to antigen presentation and immune response regulation.
Observed DNA methylation levels could precede periodontal disease development or be the consequence of periodontal disease. The main limitation of the present study is that we could not establish whether associations we detected between differential methylation in whole blood and periodontitis were the cause rather than the consequences of periodontitis. While ARIC is a prospective study and the methylation data were from blood samples collected up to six years before ARIC dental examination, a dental examination was not performed concurrently with the blood draw, such that we likely studied the combination of prevalent and incident periodontal disease. Moreover, we could not detect single CpG associations at the stringent epigenome-wide significance level. Another limitation is using the 450K array, which covers only ∼2% of human CpG sites and limits the discovery. Lastly, our use of edentulism as a surrogate for the most severe type of periodontitis is not necessarily accurate, given that the cause of edentulism is multi-faceted, including but not limited to having and being treated for periodontal disease and poor access to dental care. Consequently, our analysis cannot differentiate DNA methylation patterns of those who are edentulous due to having periodontitis vs. those who are edentulous for other reasons.
Despite the limitations, several strengths of this study warrant attention, including the use of clinically examined periodontal phenotypes, large sample size, ethnic diversity, and rigorous analytical methods. We also adjusted for SES, diabetes status, imputed cell composition, self-reporting smoking history, and a methylation-predicted packyears smoked to further reduce concern about confounding. Additionally, we utilized two definitions of periodontal disease to address the complexity of disease classification. Analyses in individuals of different ancestries support the need for future large-scale EWAS of periodontitis to investigate ethnicity-specific and age-specific periodontitis-related CpGs.
Along with previously published EWAS, our results suggest that gene regulation mechanisms may influence the susceptibility to and severity of periodontitis. Our well-powered EWAS and DMRcate analysis in participants drawn from the general population identified several genomic regions for which differential DNA methylation sites were associated with periodontitis. This study supports leukocyte DNA methylation for the evaluation of epigenetic patterns in periodontitis. Our study, in conjunction with a prior pilot study, implicates ZFP57 as a promising candidate for future studies to illuminate the underlying gene regulatory mechanisms linking differential DNA methylation to periodontal disease.
Data Availability
The datasets used in the current study are available from the ARIC study.