Genome wide association scan for chronic periodontitis implicates novel locus
© Feng et al.; licensee BioMed Central Ltd. 2014
Received: 4 March 2014
Accepted: 1 July 2014
Published: 9 July 2014
There is evidence for a genetic contribution to chronic periodontitis. In this study, we conducted a genome wide association study among 866 participants of the University of Pittsburgh Dental Registry and DNA Repository, whose periodontal diagnosis ranged from healthy (N = 767) to severe chronic periodontitis (N = 99).
Genotypingi of over half-million single nucleotide polymorphisms was determined. Analyses were done twice, first in the complete dataset of all ethnicities, and second including only samples defined as self-reported Whites. From the top 100 results, twenty single nucleotide polymorphisms had consistent results in both analyses (borderline p-values ranging from 1E-05 to 1E-6) and were selected to be tested in two independent datasets derived from 1,460 individuals from Porto Alegre, and 359 from Rio de Janeiro, Brazil. Meta-analyses of the Single nucleotide polymorphisms showing a trend for association in the independent dataset were performed.
The rs1477403 marker located on 16q22.3 showed suggestive association in the discovery phase and in the Porto Alegre dataset (p = 0.05). The meta-analysis suggested the less common allele decreases the risk of chronic periodontitis.
Our data offer a clear hypothesis to be independently tested regarding the contribution of the 16q22.3 locus to chronic periodontitis.
KeywordsPeriodontitis Chronic periodontitis Aggressive periodontitis Genetics Medical genetics Genetic polymorphism
Although family studies suggest that environmental factors are the major determinants of variance in chronic periodontitis [1–5], comparisons between reared-together and -apart adult monozygous twins indicate that early family environment has no appreciable influence on periodontal status of adults . Several association studies have been published over the last decade aiming to identify genetic factors contributing to chronic periodontitis ; however, the results are not necessarily the same depending on the population studied .
More recently, a genome wide association study [9, 10] including 1,020 and 4,504 participants self-defined as Whites selected from the Atherosclerosis Risk in Communities (ARIC) longitudinal cohort suggested a few novel loci to be possible contributors to chronic periodontitis although none of them reached formal statistical significance. Additionally, the two lists of associated single nucleotide polymorphisms from the ARIC studied samples [9, 10] did not obviously overlap. Divaris et al.  also included analyses based on bacterial colonization of eight species and the results suggested additional loci that may contribute to individual susceptibility of being colonized by specific bacterial groups.
In this study, we took into consideration the presence of ethnic admixture to investigate the association between genetic variation and chronic periodontitis. A genome-wide association scan for chronic periodontitis was conducted, including analysis adjusted by smoking habits and diabetes status and staged by incrementally adding samples from different ethnicities, to address the role of genes in this disease. Our results offer a clear hypothesis to be independently tested regarding the contribution to the 16q22.3 and 21q22.11 loci to chronic periodontitis.
All patients were participants in the Dental Registry and DNA Repository (DRDR) of the University of Pittsburgh School of Dental Medicine. Starting in September of 2006, all individuals that seek treatment at the University of Pittsburgh School of Dental Medicine have been invited to be part of the registry. They give written informed consent authorizing the extraction of information from their dental records. Also, they provide a saliva sample from which DNA can be extracted. Unstimulated saliva samples were obtained from all participants (individuals were asked to spit) and storedii at room temperature until being processed. No centrifugation was performed in the saliva samples. DNA was extracted according to the manufacturer’s instructions. The University of Pittsburgh Institutional Review Board approves this project and all individuals signed a written informed consent document prior to participation.
Summary of the study populations
Rio de Janeiro
Mean age (Years)
In individuals with
DNA samples were genotyped for 620,901 single nucleotide polymorphisms (SNPs)iii. Details of our power calculations are presented as Additional file 1 (Appendix 2: “Power Calculations of the Discovery Sample” and Additional file 1). The particular SNP array chosen includes SNPs that are representative for individuals of both African and European ancestry , which we considered an important aspect of the design, since the study group was comprised of individuals that are self-reported Whites or Blacks.
Association between periodontitis affection status and each single nucleotide polymorphism across the whole genome was tested using PLINK  and all analyses were adjusted for age, sex, diabetes status (yes or no), and smoking status (smoker or non-smoker), variables that are associated with distinct periodontal disease levels [15–19]. Data on ex-smokers was not consistent in all registry dataset and this variable was not used in the analysis. In the analysis of the complete dataset, we also adjusted for the principal components from an evaluation of population structure as described in the Additional file 1 (Appendix 3: ”Genome Wide Analysis,” Appendix 4: “Adjustment for Ethnicity in the genome Wide Analysis,” Additional file 1). We then repeated these analyses with samples from White individuals only (Additional file 1). To account for multiple testing, a p-value lower than 1E-07 (0.05/473,514) was considered statistically significant.
Follow up samples
Summary of the results of the genome wide association scans and independent analysis for chronic periodontitis
Minor allele frequency
Genome wide scan p-value all samples (99 affected, 767 unaffected)
Genome wide scan p-value whites only (63 affected, 543 unaffected)
Independent Porto Alegre cohort p-value (430 affected, 1,030 unaffected)
Independent Rio de Janeiro cohort p-value (183 affected, 176 unaffected)
For the 20 single nucleotide polymorphisms selected to this independent test, genotyping was carried out using TaqMan chemistry  and end-point analysis.iv All genetic markers were in Hardy-Weinberg equilibrium (data not shown). To determine the association between the disease and any allele or genotype frequency, we used logistic regression adjusted for age, sex, ethnicity, diabetes status, and smoking status using PLINK . Data on ex-smokers was not available in these datasets. The sample from Porto Alegre was also adjusted by body mass index as well since these data were available and this variable has been associated with periodontal diseases. P-values equal or lower than 0.0025 (0.05/20) were considered statistically significant for the follow up study results.
In order to derive a summary statistic for association with the four SNPs that showed a trend for association in either of the follow up studied samples from Brazil, a random-effects meta analysis model was used to estimate the odds ratio for the presence of the associated allele determined by the genome-wide association analysis. Before pooling the data, we estimated Cochran’s Q statistic, which indicates the degree of heterogeneity. There was no significant evidence of heterogeneity overall (Q = 2.7, p = 0.264). A random-effects model was used because it includes components of variance both within and between studies. Moreover, because it generally yields a wider confidence interval than a fixed-effects model, the random-effects mode is more conservative . The complete dataset from Pittsburgh was used. MedCalc version 13 was used (MedCalc Software, Ostend, Belgium).
Genome wide associations study
Follow up studies
We selected 20 single nucleotide polymorphisms that had consistent results in both analyses of the total sample and self-reported Whites only to test in two independent cohorts (Table 2). Inclusion of sex, ethnicity, diabetes status, smoking habits, and body mass index along with age in the model did not substantially change the results and data presented here are based on the simplest model adjusted only by age. The rs1477403 marker located on 16q22.3 was the only one that showed a trend for association in the cohort from Porto Alegre, Brazil [odds ratio = 1.2 (95% confidence interval 1.0-1.47); p = 0.05 for the allele distribution, Table 2]. Three markers in 21q22.11 showed a trend for association (nominal p-values lower than 0.05) with chronic periodontitis in the cohort from Rio de Janeiro (Table 2). These markers are in strong linkage disequilibrium with each other (D’ = 1.0).
In this study, 2,685 DNA samples were analyzed coming from two cohort studies and one case–control dataset. These different study designs explain the variation of periodontal disease frequency in each of the study groups ranging from 11% to 50%.
The first step of our study included a genome-wide analysis. Chronic periodontitis has a prevalence of over 47% in the United States based on NHANES data , and in general lower sample sizes are necessary to study a common disease than a rare disease . However, with the relatively modest number of affected individuals and anticipated statistical power, we implemented two strategies to improve statistical power. We included at least four controls for each case, which is considered the golden standard for the numbers of cases and controls to be collected in a case–control genetics study . The other approaching was to use cases with at least 30% sites of the mouth affected by chronic periodontitis, hence avoiding the inclusion of less severe cases. This approach is thought to maximize the variance of predictor variables (each genetic variant or X), which according to bx ± tn-m-1;α√MSE/nVx(1-R2) where MSE is the mean square error, n is the sample size, Vx is the variance of X, and (1-R2) is the proportion of the variance X not shared by any other variables in the model, will increase power and precision . While no single nucleotide polymorphism exhibited association at genome wide significance, several genomic regions showed suggestive evidence for association. However, only four genetic markers in two loci showed also a trend for association in independent experiments with different population datasets, and only one marker showed association when the samples were pulled. These results are interesting because one experiment was done in a hospital-based cohort which clinical data is obtained from different professionals and greater heterogeneity, and the following experiments were done in population-based cohorts and data were collected with experimental rigor to increase homogeneity. rs1477403 is located at 16q22.3 and in a sequence of an uncharacterized non-coding RNA (LOC100506172). The nucleotide change is not conserved in mouse, chimp, orangutan, or macaque according to the data available at UCSC Genome Browser and is unlikely to have a direct functional role, but this possibility cannot be excluded.
Our approach to select markers to follow up included comparing the top 100 results of the two genome wide scan analyses. We could have prioritized markers based on our initial power calculations. However, a fair assumption for periodontal diseases is that individual gene contributions are small and if we used odds ratio cut-offs lower than 1.5, we would likely have several hundred if not thousand possible markers to follow. Two-stage designs for manipulating ranked SNPs based on p-values have been shown to improve the rankings and to decrease overestimated significance values [25–28].We also performed a met-analysis to help interpret the results of the analyses of the four SNPs in the three population groups. If one population produces a large p-value for a given SNP when two other populations produced small ones for that SNP, it seems there are several possible reasons. One would be that the SNP is truly associated with the trait in the populations conferring the signal, but the SNP is not associated with the trait in the third population. Another possibility is that the SNP is associated with the trait in all populations, but the sample of individuals collected from one of the populations by chance happened to provide low power. A third possibility of course is that the SNP is not associated with the trait, and the two populations that showed a signal were both false positives. If there is, in truth, association in the third population, but the sample happened to display low power, then while the direction of any effect seen in the sample would be expected to be the same, it also seems not unlikely that by chance it might actually be opposite (low p-values mean effect sizes near zero in a given sample, and as such, the “effect” could be in either direction). We hypothesize that the signal for the 16q22.3 SNP is real, despite the individual analysis of one of the Brazilian populations does not indicate association (Figure 2). On the other hand, since the signals for the SNPs in 21q22.11 are not consistent (Figures 2 through 4), we hypothesize the evidence for association with this locus is a false positive. These analyses exemplify the challenge of interpreting results for these kinds of studies.
The first genome wide association study in periodontitis studied the aggressive type of the disease . This study identified an association with a marker in the locus of GLT6D1 and functional experiments suggested that reduced GATA3 binding affinity to the GLT6D1 locus could be a component of the pathophysiology of periodontitis . This locus is not one we are suggesting to be associated with chronic periodontitis. The lack of overlap between our findings and of the others [9, 10] in genome wide scanning for chronic periodontitis and the study of Schäfer et al.  is likely due to the fact that these two conditions have distinct genetic influences. We have previously shown that aggressive periodontitis aggregates in families and its most parsimonious mode of inheritance is a semi-general transmission model that allows the heterozygote transmission to vary . This is very distinct from what we see in chronic periodontitis in which no clear familial aggregation can be detected.
Our study benefits from several strengths including genome-wide single nucleotide polymorphism data and rigorous and thorough assessment of phenotypes. Genotyping and quality control/quality assurance yielded data of exceptional quality. Moreover, as one of the first genome wide association studies for chronic periodontitis reported to date, this study accomplished the principal goal (of the non-hypothesis-based genome wide association study design), of generating interest in genes and genomic regions previously unstudied in the context of oral health. However this study also highlights the challenges of identifying genes involved in common complex disease, namely, that numerous genes, mostly of small effect sizes, are likely to contribute to periodontitis, and that discovery of individual variants may be exceedingly difficult. Our study populations had a mix of individuals of both White and Black heritage and this further complicates any analysis since allele frequency may be disparate between different populations. Even though we carefully took into consideration this factor, we cannot exclude the possibility that the suggestive associations we found are influenced by variation in ethnic background of the samples. While research into the genetics of periodontitis lags behind many other prominent common complex diseases, this study provides a launching pad for future candidate gene and functional studies of periodontal diseases. The public availability of these data via online portals will facilitate the utility of this study in designing future efforts and cross-study collaborations to understand the genetics of periodontal diseases.
Our data offer a clear hypothesis to be independently tested regarding the contribution of the 16q22.3 locus to chronic periodontitis.
iPerformed in a Illumina 610-Quad platform.
iiStored in Oragene DNA Self-Collection kits (DNA Genotek Inc., Ottawa, ON, Canada).
iiiUsing the Illumina Human610-Quadv1_B BeadChip (Illumina Inc., San Diego, CA, USA).
ivPerformed on an Applied Biosystems 7900 HT Sequence Detection System machine (Applied Biosystems Inc., Foster City, CA, USA).
The authors are indebted to all the patients who enthusiastically agreed to be part of this project. Data for this study was provided by the Dental Registry and DNA Repository of the School of Dental Medicine, University of Pittsburgh. Financial support for this work was provided by NIH Grant 5TL1RR024155.
Summary of key findings
A marker in 16q22.3 is associated with chronic periodontitis in several diverse populations and can be of importance to determine the risk for the disease.
- Chung CS, Runck DW, Niswander JD, Bilben SE, Kau MCW: Genetic and epidemiologic studies of oral characteristics in Hawaii’s school children. I. caries and periodontal disease. J Dent Res. 1970, 49: 1374-1385.View ArticlePubMedGoogle Scholar
- Chung CS, Kau MCW, Chung SSC, Schendel SA: A genetic and epidemiologic study of periodontal disease in Hawaii: I. racial and other epidemiologic factors. J Periodontol Res. 1977, 12: 148-159.View ArticleGoogle Scholar
- Chung CS, Kau MCW, Chung SSC, Rao DC: A genetic and epidemiologic study of periodontal disease in Hawaii: II. genetic and environmental influence. Am J Hum Genet. 1977, 29: 76-82.PubMedPubMed CentralGoogle Scholar
- Rao DC, Chung CS, Morton NE: Genetic and environmental determinants of periodontal disease. Am J Med Genet. 1979, 4: 39-45.View ArticlePubMedGoogle Scholar
- Beaty TH, Colyer CR, Chang YC, Liang KY, Graybeal JC, Muhammad NK, Levin LS: Familial aggregation of periodontal indices. J Dent Res. 1993, 72: 544-551.View ArticlePubMedGoogle Scholar
- Michalowicz BS: Genetic and heritable risk factors in periodontal disease. J Periodontol. 1994, 65: 479-488.View ArticlePubMedGoogle Scholar
- Laine ML, Crielaard W, Loos BG: Genetic susceptibility to periodontitis. Periodontol 2000. 2012, 58: 37-68.View ArticlePubMedGoogle Scholar
- Schäfer AS, Richter GM, Nothnagel M, Manke T, Dommisch H, Jacobs G, Arit A, Rosenstiel P, Noack B, Groessner-Schreiber B, Jepsen S, Loos BG, Schreiber S: A genome-wide association study identifies GLT6D1 as a susceptibility locus for periodontitis. Hum Mol Genet. 2010, 19: 553-562.View ArticleGoogle Scholar
- Divaris K, Monda KL, North KE, Olshan AF, Lange EM, Moss K, Barros SP, Beck JD, Offenbacher S: Genome-wide association study of periodontal pathogen colonization. J Dent Res. 2012, 91 (S1): 21S-28S.View ArticlePubMedGoogle Scholar
- Divaris K, Monda KL, North KE, Olshan AF, Reynolds LM, Hsueh WC, Lange EM, Moss K, Barros SP, Weyant RJ, Liu Y, Neuman AB, Beck JD, Offenbacher S: Exploring the genetic basis of chronic periodontitis: a genome-wide association study. Hum Mol Genet. 2013, 22: 2312-2324.View ArticlePubMedPubMed CentralGoogle Scholar
- Susin C, Dalla Vecchia CF, Oppermann RV, Haugejorden O, Albandar JM: Periodontal attachment loss in an urban population of Brazilian adults: effect of demographic, behavioral, and environmental risk indicators. J Periodontol. 2004, 75: 1033-1041.View ArticlePubMedGoogle Scholar
- Susin C, Haas AN, Valle PM, Oppermann RV, Albandar JM: Prevalence and risk indicators for chronic periodontitis in adolescents and young adults in south Brazil. J Clin Periodontol. 2011, 8: 326-333.View ArticleGoogle Scholar
- Tandon A, Patterson N, Reich D: Ancestry informative marker panels for African Americans based on subsets of commercially available SNP arrays. Genet Epidemiol. 2011, 35: 80-83.View ArticlePubMedPubMed CentralGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575.View ArticlePubMedPubMed CentralGoogle Scholar
- Waerhaug J: Prevalence of periodontal disease in Ceylon: association with age, sex, oral hygiene, socio-economic factors, vitamin deficiencies, malnutrition, betel and tobacco comsumption and ethnic group. Final report. Acta Odontol Scand. 1967, 25: 205-231.View ArticlePubMedGoogle Scholar
- Henry JL, Sinkford JC: The economic and social impact of periodontal disease. Pub Health Rep. 1979, 94: 172-181.Google Scholar
- Summers CJ, Oberman A: Association of oral disease with 12 selected variables: I. periodontal disease. J Dent Res. 1968, 47: 457-462.View ArticlePubMedGoogle Scholar
- Ainamo J: The seeming effect of tobacco consumption on the occurrence of periodontal disease and dental caries. Suom Hammaslaak Toim. 1971, 67: 87-94.PubMedGoogle Scholar
- Golomb IM: An evaluation of the relation of diabetes mellitus to periodontal disease. New York State Dent J. 1949, 15: 525-528.Google Scholar
- Ranade K, Chang MS, Ting CT, Pei D, Hsiao CF, Olivier M, Pesich R, Hebert J, Chen YD, Dzau VJ, Curb D, Olshen R, Risch N, Cox DR, Botstein D: High-throughput genotyping with single nucleotide polymorphisms. Genome Res. 2001, 11: 1262-1268.PubMedPubMed CentralGoogle Scholar
- Berlin JA, Laird NM, Sacks HS, Chalmers TC: A comparison of statistical methods for combining event rates from clinical trials. Stat Med. 1989, 8: 141-151.View ArticlePubMedGoogle Scholar
- Eke PI, Dye BA, Wei L, Thornton-Evans GO, Genco RJ, CDC Periodontal Disease Surveillance workgroup: Prevalence of periodontitis in adults in the United States: 2009 and 2010. J Dent Res. 2012, 91: 914-920.View ArticlePubMedGoogle Scholar
- Hong EP, Park JW: Sample size and statistical power calculation in genetic association studies. Genomics Inf. 2012, 10: 117-122.View ArticlePubMed CentralGoogle Scholar
- Mackinnon S: Increasing statistical power in psychological research without increasing sample size. Open Science Collaboration 2013 (http://osc.centerforopenscience.org/2013/11/03/Increasing-statistical-power/)
- Li J: Prioritize and select SNPs for association studies with multi-stage designs. J Comp Biol. 2007, 15: 241-257.View ArticleGoogle Scholar
- Li Q, Yu K, Li Z, Zheng G: Max-rank: a simple and robust genome-wide scan for case–control association studies. Hum Genet. 2008, 123: 617-623.View ArticlePubMedGoogle Scholar
- Stronberg U, Bjork J, Vineis P, Broberg K, Zeggini E: Ranking of genome-wide association scan signals by different measures. Int J Epidemiol. 2009, 38: 1364-1373.View ArticleGoogle Scholar
- Roshan U, Chikkagoudar S, Wei Z, Wang K, Hakopnarson H: Ranking causal variants and associated regions in genome-wide association studies by the support vector machine and random forest. Nucleic Acids Res. 2011, 39: e62-View ArticlePubMedPubMed CentralGoogle Scholar
- Carvalho FM, Tinoco EM, Govil M, Marazita ML, Vieira AR: Aggressive periodontitis is likely influenced by a few small effect genes. J Clin Periodontol. 2009, 36: 468-473.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6831/14/84/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.