Skip to main content

Heritable patterns of tooth decay in the permanent dentition: principal components and factor analyses



Dental caries is the result of a complex interplay among environmental, behavioral, and genetic factors, with distinct patterns of decay likely due to specific etiologies. Therefore, global measures of decay, such as the DMFS index, may not be optimal for identifying risk factors that manifest as specific decay patterns, especially if the risk factors such as genetic susceptibility loci have small individual effects. We used two methods to extract patterns of decay from surface-level caries data in order to generate novel phenotypes with which to explore the genetic regulation of caries.


The 128 tooth surfaces of the permanent dentition were scored as carious or not by intra-oral examination for 1,068 participants aged 18 to 75 years from 664 biological families. Principal components analysis (PCA) and factor analysis (FA), two methods of identifying underlying patterns without a priori surface classifications, were applied to our data.


The three strongest caries patterns identified by PCA recaptured variation represented by DMFS index (correlation, r = 0.97), pit and fissure surface caries (r = 0.95), and smooth surface caries (r = 0.89). However, together, these three patterns explained only 37% of the variability in the data, indicating that a priori caries measures are insufficient for fully quantifying caries variation. In comparison, the first pattern identified by FA was strongly correlated with pit and fissure surface caries (r = 0.81), but other identified patterns, including a second pattern representing caries of the maxillary incisors, were not representative of any previously defined caries indices. Some patterns identified by PCA and FA were heritable (h2 = 30-65%, p = 0.043-0.006), whereas other patterns were not, indicating both genetic and non-genetic etiologies of individual decay patterns.


This study demonstrates the use of decay patterns as novel phenotypes to assist in understanding the multifactorial nature of dental caries.

Peer Review reports


Dental caries is a disease affecting most adults and caused by the complex interplay of numerous environmental, behavioral [1, 2], and genetic risk factors [312]. The etiology of dental caries is further complicated by the non-uniform risk across tooth surfaces of the full dentition leading to distinct patterns of dental decay, as previously described [1325]. Patterns of decay have been used to explore caries etiology under the assumption that different risk factors lead to distinct caries patterns. A well-known example is the maxillary anterior pattern of decay (i.e., "baby bottle" caries) in young children due in part to feeding behaviors [20, 24]. Despite prevailing evidence of the importance of caries patterns, the most common indices used for studying the epidemiology of caries are DMFT and DMFS (i.e., counts of the number of decayed, missing, or filled teeth/surfaces), which do not assess specific decay patterns. As global measures of tooth decay, DMFT and DMFS indices may not be optimal for investigating genetic and environmental factors that manifest as specific patterns of caries across the dentition. Separating the global level of caries into components or patterns with distinct etiologies may be critical for identifying risk factors of modest effect sizes, such as specific genetic loci contributing to tooth decay.

Previous descriptions of caries patterns have usually assumed and compared a priori classifications of tooth surfaces [1422, 25], which often differed among studies, leading to inconsistencies that demonstrate the limited utility of a priori surface classifications. A few studies have modeled the patterns of childhood tooth decay without a priori assumptions and have identified distinct patterns reflecting caries of the maxillary incisor surfaces and pit and fissure surfaces, among others [13, 23, 24].

To our knowledge, no assessment of permanent dentition caries patterns in adults without a priori surface classifications has previously been performed. In this study we utilized two related analytic methodologies for identifying the underlying patterns within our dataset: principal components analysis (PCA) and factor analysis (FA). Three specific purposes of this study were (1) to identify the patterns of dental caries in the permanent dentition of adults without a priori assumptions about tooth surface classifications; (2) to determine the relationship between identified patterns of decay and a priori measures of decay such as DMFS index, decay of pit and fissure surfaces, and decay of smooth surfaces; and (3) to assess the heritability of identified patterns of decay.


Recruitment and data collection

The Center for Oral Health Research in Appalachia (COHRA) was created to identify the community-, family-, and individual-level predictors of oral health outcomes in the Appalachian population [26], a vulnerable subpopulation with poorer oral health compared to the greater US population [2729]. COHRA participants were recruited by household as previously described [6, 7, 26], whereby eligible households were required to include at least one biological parent-offspring pair with the child being 1 to 18 years of age. All members of eligible households were invited to participate without regard to their oral health status, demography, or biological or legal relationships. Written informed consent was provided by all adult participants. Assent with parent or guardian written consent was provided on behalf of all child participants. The study was approved by the COHRA research committee and the Institutional Review Boards of the University of Pittsburgh and West Virginia University.

In total, 732 households were recruited, which comprised 2,663 individuals from 740 biological kinships of 1 to 20 family members (mean = 4.72 members). Some kinships spanned multiple households, whereas other households contained multiple kinships. Reported familial relationships were validated using panels of ancestry-informative [30] and whole-genome [31] genetic marker data provided by the Center for Inherited Disease Research at Johns Hopkins University and quality checked jointly by study investigators and the Coordinating Center for the NIH Genes and Environment Initiative (GENEVA; [32]).

Dental caries was assessed via visual inspection with a dental explorer during intra-oral dental examinations conducted by dentists or research dental hygienists calibrated with respect to a reference dentist at least once per year. Inter- and intra-examiner concordances of caries assessments were high [7, 26]. Each tooth surface was scored as sound, pre-cavitated, decayed, filled, missing due to decay, or missing due to reasons other than decay, in accordance with the World Health Organization DMFS/dfs scale and in accordance with the NIH/NIDCR-approved protocol for assessing dental caries for research purposes [33]. This method of caries assessment is compatible with that recommended by the PhenX Toolkit (; designed to facilitate combining data across studies), and the National Center for Health Statistics Dental Examiners Procedures Manual (See Section [34]. Third molars were excluded from caries assessment. Edentulous individuals were recruited into the study but were excluded from caries assessment and analysis.

Statistical analysis

The analytic goal of the present study was to explore patterns of dental caries of the permanent dentition in adults. Therefore we excluded children by restricting our study sample to the 1,068 participants aged 18 to 75 years. For each participant, surface-level caries data on 128 surfaces (i.e., 4 surfaces for each incisor and canine, and 5 surfaces for each premolar and molar) were coded as 0 for sound or missing due to reasons other than decay, or coded as 1 for pre-cavitated, decayed, missing due to decay, or filled/restored. Thus, we generated a matrix of 1,068 participants by 128 indicators of surface-level caries affection status. This matrix was used as input for two related methods of extracting patterns within the data: PCA and FA [35].

PCA uses singular value decomposition of the data matrix to extract a set of uncorrelated variables (called principal components scores, PCs) where the first PC (i.e., PC1) explains the greatest possible amount of variability in the data in a single dimension, and the second PC (i.e., PC2) explains the greatest possible amount of remaining variability in the data in a single dimension orthogonal to PC1, and so on. The result is a number of orthogonal PCs equal to the number of original variables (in our data, 128), with successive PCs each explaining less and less of the data variability. Each PC can be defined as a linear combination of the original variables weighted by their loadings. The first several PCs may represent important patterns in the data, essentially assessing underlying signals from a greater number of correlated phenotype measurements. The loadings provide a way of interpreting the PCs in terms of the original variables. In other words, the loadings describe the pattern of carious lesions across the permanent dentition for a given PC, whereas the actual PCs indicate the extent/severity of caries of that decay pattern.

FA is similar to PCA in that it is used to extract latent variables called factor scores (FACs) from an original data matrix. Like PCs, FACs are calculated as linear combinations of the original variables weighted by their loadings, except that the number of FACs used to model the patterns in the data is chosen a priori, and the FACs are not constrained to be orthogonal. In this study, we modeled the caries data matrix using 10 factors. Like PCA, the goal of FA is to generate FACs representing underlying signals in the data matrix that can then be used as phenotypes, in this case, to identify the risk factors for dental caries.

In practice, FA and PCA often perform similarly. However the two methods take opposite perspectives in extracting patterns from a data matrix: PCA assumes that the observed variables provide the basis for the patterns, whereas FA assumes that latent patterns provide the basis for the observed variables. In this way, PCA is often used for dimension reduction, i.e., summarizing the information from a large number of variables with a few variables, whereas FA may better represent underlying "endophenotypes", i.e., unmeasured phenotypes that manifest as the observed variables. For both PCA and FA, the loadings define the patterns of decay and the PCs and FACs describe the severity of disease for their corresponding patterns.

For comparison to the PCs and FACs, we also generated three a priori caries phenotypes: the DMFS index, pit and fissure surface caries (PFS), and smooth surface caries (SMS). These a priori phenotypes are commonly used in the caries literature. DMFS was calculated as the number of pre-cavitated, decayed, missing due to decay, or filled/restored surfaces. PFS and SMS were calculated in the same way as DMFS except that counts were limited to pit and fissure surfaces and smooth surfaces, respectively. Occlusal surfaces of the premolars and molars, buccal surfaces of the maxillary molars, and lingual surfaces of the mandibular molars were considered pit and fissure surfaces. All other tooth surfaces were considered smooth surfaces.

In order to assess the stability of patterns identified by PCA and FA, we performed a sensitivity analysis by repeating PCA and FA on ten random subsets of the data comprised of 80% of the full sample. We compared the PCs and FACs obtained from random subsets to those from the full sample using the Pearson correlation coefficient, r. PCs 1-4 were extremely stable (r = 0.98 to 1.00), PCs 5-9 were stable (r = 0.86 to 0.95), and PC 10 was moderately stable (r = 0.77) across random subsets. FACs 1-6 were stable (r = 0.86 to 0.99), and FACs 7-10 were moderately stable (r = 0.69 to 0.82) across random subsets. Likewise, we assessed the effect of relatives on PCA and FA by repeating these methods in the maximal subset of unrelated individuals. PCs 1-10 and FACs 1-8 from the unrelated sample were highly correlated (r > 0.95) with those from the full sample, whereas FAC9 and FAC10 were moderately correlated (r = 0.57, and 0.81, respectively). Altogether, these results suggest that caries patterns were generally stable and robust to the inclusion of relatives among the sample.

Heritability estimates of PCs and FACs were calculated using the variance components approach. This method models phenotype correlations among all types of relatives as a function of the expected degree of genetic sharing (i.e. that parents and offspring share 50% of their genome, siblings share 50%, half-siblings share 25%, unrelated individuals share 0%, etc.). Details for this method as applied to our study sample have previously been reported [6, 36]. The heritability estimate is interpreted as the proportion of phenotype variance attributable to the cumulative effect of all genes.

All statistical analyses were performed in the R software package (R Foundation for Statistical Computing, Vienna, AU), except heritability estimates which were obtained from genetic modeling performed in SOLAR [37]. Principal components analysis was performed using the prcomp function with default parameters. Factor analysis was performed using the factanal function with the Thomson's regression-based scores option, 10 factors, and other default parameters. Prevalences, correlations, and figures were all generated in R.


Caries prevalences by surface

Surface-level caries data for 1,068 participants (ages 18 to 75 years, mean age of 34.7 years, 63.3% female, 90.0% self-reported white) across 128 tooth surfaces were collected. Tooth surfaces that exhibited evidence of pre-cavitated lesions or decay, were missing due to decay, or had been filled/restored, were considered carious. Tooth surfaces that were sound or missing due to reasons other than decay were considered non-carious. Caries prevalences per surface (i.e. the proportion of the sample exhibiting caries on a particular tooth surface) are shown in Table 1. Surfaces of the anterior maxillary teeth (i.e., incisor and canines) exhibited greater caries prevalences than anterior mandibular teeth; whereas posterior maxillary teeth (i.e., premolars and molars) exhibited lower pravelences rates than posterior mandibular teeth.

Table 1 Caries prevalences per surface across the permanent dentition (N = 1,068)

Principal components analysis

PCA was performed on the surface-level data in order to extract the underlying patterns of caries. PC1 explained 26.3% of the variability in the data, PC2 explained 6.7%, and all other PCs explained < 5% (Figure 1A). Loadings show that except for anterior mandibular surfaces, all other surfaces contribute similarly to PC1 (Figure 1B) representing a near-global pattern/extent of decay. Loadings for PC2 show opposite contributions of smooth surfaces and pit and fissure surfaces (Figure 1C). Loadings for PC3 show opposite contributions of premolar vs. other surfaces and loadings for PC4 show opposite contributions of maxillary vs. mandibular surfaces (see Additional file 1). Loadings for all other PCs show complex patterns of contributions from tooth surfaces that are not easily discernible in the context of PCs 1 to 4, however, general descriptions of the contributing surfaces are summarized in Table 2.

Figure 1

Principal component analysis. (A) proportion of data variance explained by PCs 1-10 (red) and successive PCs (black). (B) Loadings for PC1 ordered by tooth type, from left to right: maxillary incisors (blue), canines (green), premolars (red), molars (black), mandibular incisors (gray), canines (yellow), premolars (magenta), molars (cyan). For each tooth, contributions of surfaces are listed in the following order: buccal, distal, lingual, mesial, and occlusal, if applicable. Teeth ordered from left to right: maxillary 8, 9, 7, 10, 6, 11, 5, 12, 4, 13, 3, 14, 2, 15; mandibular 24, 25, 23, 26, 22, 27, 21, 28, 20, 29, 19, 30, 18, 31. (C) Loadings for PC2, in the same order; smooth surfaces shaded red, pit and fissure surfaces shaded blue. Scatter plots of (D) PC1 vs. DMFS index, (E) PC2 vs. pit and fissure caries, (F) PC3 vs. smooth surface caries.

Table 2 General interpretations of PCA and FA loadings.

PC1 was nearly identical to DMFS index (r = 0.969; p-value < 10-250 [i.e., the minimum p-value reported using the statistics software]; Figure 1D) indicating that the strongest pattern of caries in the data distinguished individuals by global level of decay. PFS, the count of carious pit and fissure surfaces, was very highly correlated with PC2 after subtracting out PC1 (r = 0.947; p-value < 10-250; Figure 1E). SMS, the count of carious smooth surfaces, was highly correlated with PC3 after subtracting out PC2 and PC1 (r = 0.894; p-value < 10-250; Figure 1E). These correlations show that PC1, PC2, and PC3 capture the patterns of dental decay corresponding to a priori phenotypes, DMFS index, excess PFS (for a given DFMS), and excess SMS (for given DMFS and PFS), respectively.

The heritability (h2) of DMFS index and PCs 1-10 were calculated while simultaneously adjusting for the effects of age, age2, and sex (Table 3). DMFS index, PC1, PC5, and PC7 were all strongly heritable (h2 = 37% to 50%; p-values = 0.043 to 0.008) indicating that some patterns of dental decay were due to genetic etiologies. Other PCs were not heritable indicating that some patterns of dental decay were not due to genetics. Covariates age, age2, and sex explained about 10% of variation in PC1 and very little variation for the remaining PCs.

Table 3 Heritability estimates for DMFS index, PCs, and FACs

Factor analysis

FA was also performed on the surface-level data to identify latent patterns of dental decay (Table 2). 10 factors were extracted which cumulatively explained 44.7% of the variability of the data. FAC1 was primarily due to the contributions of molar surfaces, and to a lesser degree, premolar surfaces (see loadings, Figure 2A). FAC1 was moderately correlated with DMFS index (r = 0.593; p-value < 10-250), and strongly correlated with PFS (r = 0.815, p-value < 10-250; Figure 2B). Loadings showed that maxillary incisor surfaces, and to a lesser degree, maxillary canine surfaces, contribute to FAC2 (Figure 2C). FAC2 was moderately correlated with SMS (r = 0.523; p-value < 10-250) and DMFS index (r = 0.453; p-value < 10-250). See Additional file 1 for loadings of all other factors. In general, most FACs showed low correlations with PCs, indicating that the two methods extracted different patterns from the data. Compared with the PCs, which represented contributions from many teeth, a number of FACs primarily represented contributions of individual teeth (e.g., tooth 20 for FAC7, tooth 29 for FAC8, tooth 13 for FAC10).

Figure 2

Factor analysis. (A) Loadings for FAC1 ordered by tooth type, from left to right: maxillary incisors (blue), canines (green), premolars (red), molars (black), mandibular incisors (gray), canines (yellow), premolars (magenta), molars (cyan). For each tooth, contributions of surfaces are listed in the following order: buccal, distal, lingual, mesial, and occlusal, if applicable. Teeth ordered from left to right: maxillary 8, 9, 7, 10, 6, 11, 5, 12, 4, 13, 3, 14, 2, 15; mandibular 24, 25, 23, 26, 22, 27, 21, 28, 20, 29, 19, 30, 18, 31. (B) Scatter plot of FAC1 vs. pit and fissure surface caries. (C) Loadings for FAC2 in the same order.

The heritability estimates of FACs 1-10 are also shown in Table 3. FAC3 and FAC6 were strongly heritable (h2 = 65.3 and 30.2%; p-value = 0.006 and 0.027, respectively), whereas all other FACs were not heritable. These results echo the PCA results, showing that some caries patterns are due to genetic etiologies, whereas others are not. Significance levels for heritability estimates did not meet Bonferroni adjustment (for 20 models, requiring p-values < 0.0025 for family-wise significance); although, correct adjustment for multiple testing is not clear given the prior significant heritability of DMFS, PFS, and SMS indices reported for this sample [6, 36].


We used two related methods of extracting caries patterns in the permanent dentition from surface-level caries data. PCA yielded many moderate-to-weak patterns, possibly indicating a high degree of noise or sporadic (non-patterned) occurrence of dental caries. Moreover, PCs 1-3 closely recaptured the DMFT, PFS, and SMS indices, an observation that suggests these a priori caries phenotypes may reflect the predominant patterns of decay in the permanent dentition, although cumulatively they account for only 37% of the variability. Some PCs were heritable, whereas many were not, which suggests that genetic patterns of decay may be separable from non-genetic patterns. Unlike PCA, FA did not yield factors that clearly recaptured a priori phenotypes, with the exception that FAC1 was correlated with PFS. Maxillary incisors contributed heavily to FAC2, which is consistent with previous studies that used multidimensional scaling [24] and cluster analysis [23] to explore caries patterns in the primary dentition and showed maxillary incisors formed the second cluster (after other smooth surfaces). Ten factors were insufficient to explain the variability of the data, cumulatively accounting for approximately 45%.

Like PCA, FA yielded some factors that were highly heritable indicating that certain caries patterns may be due to genetic etiologies while others may be due to non-genetic etiologies. Because the caries patterns presented in this manuscript are more precisely and agnostically defined than a priori phenotypes, we conservatively conclude that specific patterns represented by FAC3 and FAC6 are heritable, rather than generalizing to broader surface categories such as SMS. Interestingly, the strongest genetic contribution identified was for FAC3, which was 65.3% heritable (compared to 41.8% for D1MFS index) which suggests that FAC3 may be a better phenotype for gene discovery than a priori caries phenotypes. A similar conclusion can be made for PC7 (50.3% heritable). These results are generally consistent with a previous study comparing PCA and FA that showed FA may better capture underlying genetic signals from correlated phenotype measurements (although both methods perform quite similarly) [35]. Non-heritable PCs and FACs, presumably due to effects of non-genetic risk factors, may be preferred phenotypes for future epidemiological studies of environmental risk factors for dental caries.

The severity of caries significantly increased with age (or age2) for most patterns (results not shown). Heritability estimates were calculated while simultaneously modeling age, age2, and sex, although very similar heritability estimates were obtained in unadjusted models for all patterns except PC1 which exhibited decreased heritability when covariates were omitted (results not shown). These results are sensible given that altogether, age, age2 and sex accounted for about 10% of variance in PC1, but very little variance for the other PCs and FACs.

One of the challenges of using agnostic methods such as PCA and FA to identify underlying patterns of dental decay (devoid of a priori surface classifications) is in interpreting the findings. While some patterns, such as PC1 (defined by near-uniform loadings across most tooth surfaces and therefore representing global extent of decay), and FACs 7, 8 and 10 (each defined by contributions of a single pre-molar), are readily interpretable, other PCs and FACs may be difficult to relate back to the original variables. Moreover, there is no clear method of distinguishing biologically relevant patterns attributable to distinct risk factors from sporadic patterns due to noise. Sensitivity analysis showed that patterns represented by PCs 1-9 and FACs 1-6 were stable, whereas PC10 and FACs 7-10 were moderately stable. The overall stability lends credence to the notion that PCs and FACs considered in this study are not due to chance alone.

This study benefits from the large sample of related individuals with detailed surface-level caries assessment, which facilitated caries pattern extraction and heritability estimation. An additional strength of the analysis was using two different but related methods of extracting caries patterns from the data, which, most importantly, did not use a priori pattern definitions.

Despite these strengths, several limitations of this study warrant discussion, including inherent limitations to assigning tooth surfaces as carious or not. First, caries assessment by visual inspection, though suitable for obtaining data on large numbers of individuals and of sufficient quality for research purposes, may under-represent the true level of disease. Moreover, teeth missing due to decay, for which all surfaces count as carious, and approximal lesions which are often treated by two-surface restorations (leading to filled occlusal surfaces despite absence of decay) may cause caries assessment errors. Likewise, the quality of caries assessment may not be uniform across surfaces of the permanent dentition, which may have caused additional "noise" in the caries measurement. Lastly, prophylactic restorations may inflate caries assessment. These limitations are unavoidable for cross-sectional (i.e., single time point) study designs of dental caries. However, appropriate modeling techniques, such as methods of pattern extraction including PCA and FA, may aid in overcoming theses limitations of the caries assessment.


To our knowledge, this study is the first exploration of caries patterns in the permanent dentition in adults without relying on a priori assumptions or surface classifications. Overall, this study demonstrates the utility of methods for extracting caries patterns from surface-level data and reinforces the complexity of dental caries etiology. Because risk factors that manifest as specific decay patterns may otherwise go unobserved with respect to global or other a priori caries phenotypes, the use of patterns as novel phenotypes may assist in understanding the multifactorial nature of dental caries. This study is one of few but much needed efforts to use decay patterns to define new phenotypes for studying dental caries.



Decay missing filled surfaces


Decay missing filled teeth


Principal components analysis


Factor analysis


Principal components score


Factor score


Pit and fissure surface caries


Smooth surface caries.


  1. 1.

    Hunter PB: Risk factors in dental caries. Int Dent J. 1988, 38 (4): 211-217.

    PubMed  Google Scholar 

  2. 2.

    Anderson M: Risk assessment and epidemiology of dental caries: review of the literature. Pediatric Dentistry. 2002, 24 (5): 377-385.

    PubMed  Google Scholar 

  3. 3.

    Shuler CF: Inherited risks for susceptibility to dental caries. J Dent Educ. 2001, 65 (10): 1038-1045.

    PubMed  Google Scholar 

  4. 4.

    Slayton RL, Cooper ME, Marazita ML: Tuftelin, mutans streptococci, and dental caries susceptibility. J Dent Res. 2005, 84 (8): 711-714. 10.1177/154405910508400805.

    Article  PubMed  Google Scholar 

  5. 5.

    Vieira AR, Marazita ML, Goldstein-McHenry T: Genome-wide scan finds suggestive caries loci. J Dent Res. 2008, 87 (5): 435-439. 10.1177/154405910808700506.

    Article  PubMed  Google Scholar 

  6. 6.

    Wang X, Shaffer JR, Weyant RJ, Cuenco KT, Desensi RS, Crout R, McNeil DW, Marazita ML: Genes and their effects on dental caries may differ between primary and permanent dentitions. Caries Res. 2010, 44 (3): 277-284.

    Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Wendell S, Wang X, Brown M, Cooper ME, DeSensi RS, Weyant RJ, Crout R, McNeil DW, Marazita ML: Taste genes associated with dental caries. J Dent Res. 2010, 89 (11): 1198-1202. 10.1177/0022034510381502.

    Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Boraas JC, Messer LB, Till MJ: A genetic contribution to dental caries, occlusion, and morphology as demonstrated by twins reared apart. J Dent Res. 1988, 67 (9): 1150-1155. 10.1177/00220345880670090201.

    Article  PubMed  Google Scholar 

  9. 9.

    Bretz WA, Corby PM, Hart TC, Costa S, Coelho MQ, Weyant RJ, Robinson M, Schork NJ: Dental caries and microbial acid production in twins. Caries Res. 2005, 39 (3): 168-172. 10.1159/000084793.

    Article  PubMed  Google Scholar 

  10. 10.

    Bretz WA, Corby PM, Melo MR, Coelho MQ, Costa SM, Robinson M, Schork NJ, Drewnowski A, Hart TC: Heritability estimates for dental caries and sucrose sweetness preference. Arch Oral Biol. 2006, 51 (12): 1156-1160. 10.1016/j.archoralbio.2006.06.003.

    Article  PubMed  Google Scholar 

  11. 11.

    Bretz WA, Corby PM, Schork NJ, Robinson MT, Coelho M, Costa S, Melo Filho MR, Weyant RJ, Hart TC: Longitudinal analysis of heritability for dental caries traits. J Dent Res. 2005, 84 (11): 1047-1051. 10.1177/154405910508401115.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Werneck RI, Lazaro FP, Cobat A, Grant AV, Xavier MB, Abel L, Alcais A, Trevilatto PC, Mira MT: A major gene effect controls resistance to caries. J Dent Res. 2011, 90 (6): 735-739. 10.1177/0022034510397614.

    Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Batchelor PA, Sheiham A: Grouping of tooth surfaces by susceptibility to caries: a study in 5-16 year-old children. BMC Oral Health. 2004, 4 (1): 2-10.1186/1472-6831-4-2.

    Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Douglass JM, Yi W, Xue ZB, Tinanoff N: Dental caries in preschool Beijing and Connecticut children as described by a new caries analysis system. Community Dent Oral Epidemiol. 1994, 22 (2): 94-99. 10.1111/j.1600-0528.1994.tb01580.x.

    Article  PubMed  Google Scholar 

  15. 15.

    Greenwell AL, Johnsen D, DiSantis TA, Gerstenmaier J, Limbert N: Longitudinal evaluation of caries patterns form the primary to the mixed dentition. Pediatr Dent. 1990, 12 (5): 278-282.

    PubMed  Google Scholar 

  16. 16.

    Johnsen DC: Dental caries patterns in preschool children. Dent Clin North Am. 1984, 28 (1): 3-20.

    PubMed  Google Scholar 

  17. 17.

    Johnsen DC, Bhat M, Kim MT, Hagman FT, Allee LM, Creedon RL, Easley MW: Caries levels and patterns in head start children in fluoridated and non-fluoridated, urban and non-urban sites in Ohio, USA. Community Dent Oral Epidemiol. 1986, 14 (4): 206-210. 10.1111/j.1600-0528.1986.tb01536.x.

    Article  PubMed  Google Scholar 

  18. 18.

    Johnsen DC, Schechner TG, Gerstenmaier JH: Proportional changes in caries patterns from early to late primary dentition. J Public Health Dent. 1987, 47 (1): 5-9. 10.1111/j.1752-7325.1987.tb01951.x.

    Article  PubMed  Google Scholar 

  19. 19.

    Johnsen DC, Schubot D, Bhat M, Jones PK: Caries pattern identification in primary dentition: a comparison of clinician assignment and clinical analysis groupings. Pediatr Dent. 1993, 15 (2): 113-115.

    PubMed  Google Scholar 

  20. 20.

    Johnsen DC, Schultz DW, Schubot DB, Easley MW: Caries patterns in Head Start children in a fluoridated community. J Public Health Dent. 1984, 44 (2): 61-66. 10.1111/j.1752-7325.1984.tb03049.x.

    Article  PubMed  Google Scholar 

  21. 21.

    O'Sullivan DM, Tinanoff N: Maxillary anterior caries associated with increased caries risk in other primary teeth. J Dent Res. 1993, 72 (12): 1577-1580. 10.1177/00220345930720120801.

    Article  PubMed  Google Scholar 

  22. 22.

    O'Sullivan DM, Tinanoff N: The association of early dental caries patterns with caries incidence in preschool children. J Public Health Dent. 1996, 56 (2): 81-83. 10.1111/j.1752-7325.1996.tb02401.x.

    Article  PubMed  Google Scholar 

  23. 23.

    Psoter WJ, Pendrys DG, Morse DE, Zhang HP, Mayne ST: Caries patterns in the primary dentition: cluster analysis of a sample of 5169 Arizona children 5-59 months of age. Int J Oral Sci. 2009, 1 (4): 189-195. 10.4248/IJOS09077.

    Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Psoter WJ, Zhang H, Pendrys DG, Morse DE, Mayne ST: Classification of dental caries patterns in the primary dentition: a multidimensional scaling analysis. Community Dent Oral Epidemiol. 2003, 31 (3): 231-238. 10.1034/j.1600-0528.2003.00044.x.

    Article  PubMed  Google Scholar 

  25. 25.

    Douglass JM, Tinanoff N, Tang JM, Altman DS: Dental caries patterns and oral health behaviors in Arizona infants and toddlers. Community Dent Oral Epidemiol. 2001, 29 (1): 14-22. 10.1034/j.1600-0528.2001.00004.x.

    Article  PubMed  Google Scholar 

  26. 26.

    Polk DE, Weyant RJ, Crout RJ, McNeil DW, Tarter RE, Thomas JG, Marazita ML: Study protocol of the Center for Oral Health Research in Appalachia (COHRA) etiology study. BMC Oral Health. 2008, 8: 18-

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    From the Centers for Disease Control and Prevention. Total tooth loss among persons aged > or = 65 years-selected states, 1995-1997. JAMA. 1999, 281 (14): 1264-1266.

  28. 28.

    Public health and aging: retention of natural teeth among older adults--United States, 2002. MMWR Morb Mortal Wkly Rep. 2003, 52 (50): 1226-1229.

  29. 29.

    Purnell LDCM: Appalachians. Transcultural Health Care: A Culturally Competent Approac. Edited by: Purnell LD PB. 1998, Philadelphia: F. A. Davis

    Google Scholar 

  30. 30.

    Marosy B, Romm J, Hetrick K, Doheny K, Pugh E, Tsai Y: Development of a low cost SNP barcode panel [abstract #2653 F]. Presented at the American Society of Human Genetics 57th Annual Meeting. 2007, San Diego, CA

    Google Scholar 

  31. 31.

    Cornelis MC, Agrawal A, Cole JW, Hansel NN, Barnes KC, Beaty TH, Bennett SN, Bierut LJ, Boerwinkle E, Doheny KF: The Gene, environment association studies consortium (GENEVA): maximizing the knowledge obtained from GWAS by collaboration across studies of multiple conditions. Genet Epidemiol. 2010, 34 (4): 364-372. 10.1002/gepi.20492.

    Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Laurie CC, Doheny KF, Mirel DB, Pugh EW, Bierut LJ, Bhangale T, Boehm F, Caporaso NE, Cornelis MC, Edenberg HJ: Quality control and quality assurance in genotypic data for genome-wide association studies. Genet Epidemiol. 2010, 34 (6): 591-602. 10.1002/gepi.20516.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Drury TF, Horowitz AM, Ismail AI, Maertens MP, Rozier RG, Selwitz RH: Diagnosing and reporting early childhood caries for research purposes. A report of a workshop sponsored by the National Institute of Dental and Craniofacial Research, the Health Resources and Services Administration, and the Health Care Financing Administration. J Public Health Dent. 1999, 59 (3): 192-197. 10.1111/j.1752-7325.1999.tb03268.x.

    Article  PubMed  Google Scholar 

  34. 34.

    National Center for Health Statistics National Health and Nutrition Examination Survey Dental Examiners Procedures Manual. 2001

  35. 35.

    Wang X, Kammerer CM, Anderson S, Lu J, Feingold E: A comparison of principal component analysis and factor analysis strategies for uncovering pleiotropic factors. Genet Epidemiol. 2009, 33 (4): 325-331. 10.1002/gepi.20384.

    Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Shaffer JR, Wang X, DeSensi RS, Wendell S, Weyant RJ, T.Cuenco K, Crout R, McNeil DW, Marazita ML: Genetic susceptibility to dental caries on pit and fissure and smooth surfaces. Caries Res. 2012, 46: 38-46. 10.1159/000335099.

    Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62 (5): 1198-1211. 10.1086/301844.

    Article  PubMed  PubMed Central  Google Scholar 

Pre-publication history

  1. The pre-publication history for this paper can be accessed here:

Download references


The Center for Oral Health Research in Appalachia would like to thank the staff of the following heath care organizations for their efforts toward data collection: University of Pittsburgh, Bradford, Center for Rural Health Practice, Bradford, PA; McKean County Dental Center, Bradford, PA; Cornerstone Care Community Medical and Dental Center, Burgettstown, PA; UPMC Braddock Hospital, Braddock, PA; Camden-on-Gauley Medical Center, WV; Community Health Clinic of Nicholas County, WV; Richwood Area Community Hospital, WV; Summersville Memorial Hospital, WV; and the Webster County Memorial Hospital, WV. We would also like to thank the following organizations for their contributions: the GORGE Connection Rural Health Education Partnership Board, the Webster-Nicholas Rural Health Education Consortium Board, the West Virginia Rural Health Education Partnerships program, the Nicholas and Webster Boards of Education, and the UPMC Braddock Community Advisory Board. We would also like to acknowledge Kai Yu and Paula Trevilatto for their thoughtful comments on this work. Most importantly, we would like to thank the participating families for their contribution toward our goals of understanding and bettering the oral health of rural Appalachian communities.

Support for this study was provided by the National Institute of Dental and Craniofacial Research, including grants R03-DE021425, R01-DE014899 and U01-DE018903, as part of the NIH Genes and Environment Initiative (GENEVA), and by U01-HG004446 for data cleaning by the GENEVA Coordinating Center. Genotyping was provided by the NIDCR through a federal contract from the National Institutes of Health to The Johns Hopkins University, Center for Inherited Disease Research, NIH contract number HHSN268200782096C. Additional support was provided by the University of Pittsburgh School of Dental Medicine, the West Virginia University School of Dentistry and Eberly College of Arts and Sciences. The content presented herein is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of Dental and Craniofacial Research, nor the National Institutes of Health. The funding sources had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information



Corresponding author

Correspondence to John R Shaffer.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

JRS conceived and designed this study; RJW, RC, DWM, and MLM conceived and designed the COHRA initiative; JRS analyzed the data; JRS, EF, WX, KT, DEW, RSD, and MLM managed, cleaned and quality checked the data; JRS, EF, WX, KT, DEW, RSD, DEP, SW, RJW, RC, DWM, and MLM interpreted the results; JRS wrote the manuscript; JRS, EF, XW, KT, DEW, RSD, DEP, SW, RJW, RC, DWM, and MLM read, revised and approved the manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Shaffer, J.R., Feingold, E., Wang, X. et al. Heritable patterns of tooth decay in the permanent dentition: principal components and factor analyses. BMC Oral Health 12, 7 (2012).

Download citation


  • Dental caries genetics
  • Heritability
  • Permanent dentition
  • Pit and fissure surfaces
  • Smooth surfaces
  • Tooth surfaces
  • Principal components analysis
  • Factor analysis
  • Patterns of tooth decay
  • Patterns of dental caries