Abstract
Localized non-inheritable developmental defects of tooth enamel (DDE) are classified as enamel hypoplasia (EH), opacity (OP) and post-eruptive breakdown (PEB) using the Enamel Defects Index. To better understand the etiology of DDE, and in particular possibly modifiable variables, we assessed the linkages amongst exposome variables during the specific time duration of the development of the DDE. In general, the human primary central maxillary incisor teeth develop between 13-14 weeks in utero and 3-4 weeks’ postpartum of a full-term delivery, followed by tooth eruption at about 1 year of age. We utilized existing datasets of mother and child dyad data that encompassed 12 weeks’ gestation through birth and early infancy, and child DDE outcomes from digital images of the erupted primary maxillary central incisor teeth. We applied a Bayesian modeling paradigm to assess the important predictors of EH, OP, and PEB. The results of Gibbs variable selection showed a key set of predictors: mother’s pre-pregnancy body mass index (BMI); maternal serum levels of calcium and phosphorus at gestational week 28; child’s gestational age; and both mother’s and child’s functional vitamin D deficiency (FVDD). In this sample of healthy mothers and children, significant predictors for OP included the child having a gestational period > 36 weeks and FVDD at birth, and for PEB included a mother’s pre-pregnancy BMI < 21.5 and higher serum phosphorus level at week 28.
1. Introduction
Three major types of localized non-inheritable developmental defects of tooth enamel (DDE) are classified in the Enamel Defects Index as enamel hyperplasia (EH), opacities (OP) and post-eruptive breakdown (PEB) (1). Prevalence data for these defects in the primary dentition are scarce, and the results generally reflect global convenience samples with ranges of 4-99%, 2-98% and 6-50% for the three defects, respectively (2–6). One of the consequences that these developmental defects can incur is dental treatment. For example, the enamel surface irregularities of EH can provide niches for cariogenic bacteria leading to dental caries and the subsequent need for dental treatment (7, 8). Dental treatment in these young children with DDE is challenging and expensive (9–13); and 50% of the early childhood dental caries are found in the primary maxillary anterior dentition (14). For OP, the treatment considerations generally revolve around esthetics; and for PEB, the treatment considerations can involve function and/or esthetics (14–16). Prevention efforts for DDE are hampered by the lack of knowledge of potentially modifiable etiological variables during the development of the structurally defective enamel.
To address this need for knowledge of the variables and their interplay during the development of EH, OP and PEB, we focused our study on DDE in the human primary central maxillary incisor teeth. The crowns of these two teeth develop between 13-14 weeks’ in utero and 3-4 weeks’ postpartum of a full-term delivery (17, 18). To encompass the developmental period for these teeth, we assessed pertinent variables from the mother during pregnancy and at delivery; and the child at delivery through 4-6 weeks of infancy. The outcomes were determined from digital images of the erupted primary maxillary incisor teeth. We conducted a secondary analysis of maternal and child data from a randomized controlled clinical trial (RCT) of vitamin D supplementation for healthy mothers during pregnancy and birth, and child follow-up studies that included images of the child’s maxillary central incisor teeth (19–21).
2. Methods
The study sample was mother and child dyads (n=161) who were a self-selected subset from the larger RCT and follow-up studies including dental imaging with detailed description previously published (19–21). From those data resources, we chose variables that had previously been associated with DDE as evidenced by publication (15, 22–30). We also included variables unique to our datasets that purposely included the period of 12 weeks of gestation through 4-6 weeks’ early infancy to encompass the developmental duration for the primary maxillary central incisor teeth. The variables included maternal predictors: mother’s age, pre-pregnancy body mass index (BMI), total count of months of antacid use from weeks 12-36; serum circulating levels of calcium (Ca), phosphorus (P); vitamin D as 25-hydroxyvitamin D (25(OH)D) and 1,25-dihydroxyvitamin D (1,25(OH)2D); and intact parathyroid hormone (iPTH). Child factors included gestational age, cord blood serum Ca, P, iPTH, 25(OH)D, and 1-25(OH)2D levels, and vitamin D binding protein (VDBP) genotype (focusing on 1s, 1f, and 2 genotypes as the three most common VDBP alleles). Early infancy diet was determined as whether or not the child had received formula by 4-6 weeks of age or was exclusively fed breastmilk. Child DDE outcomes were scored for the facial surfaces of the two maxillary central incisor teeth from digital images using the Enamel Defects Index (1). EH was operationalized by less enamel in the frontal (buccal-lingual) plane and PEB as less enamel in the median (mesial-distal) or the transverse (incisal-coronal) plane. Each facial surface was divided into 3 nonoverlapping regions (incisal, middle, and cervical) and the presence of each defect was scored for each region, resulting in child scores as status (binary) and extent (count, with a range of 0-6 per defect type for the 2 teeth per child) (21). Decayed, restored and missing teeth or regions (due to dental caries, trauma or exfoliation) were excluded from this study’s analyses. Intra-examiner agreement for EH was determined as “substantial” at the child level by a comparative rescore of the 6 tooth regions for 15% of the children (κ = 0.779) (31).
2.1 Variable Transformations
We log-transferred the mother’s BMI before model fitting. The functional vitamin D deficiency (FVDD) group was assigned for both the mothers during pregnancy and the children at birth. This criterion was applied based on previous studies that a concurrent maternal serum 25-dihydroxyvitamin D (25(OH)D) < 20ng/mL (deficiency) and intact parathyroid hormone (iPTH) > 65 pg/mL (abnormal; elevated) provided a definitive ratio i.e., ≤ 0.308 to identify FVDD (32). The mothers had monthly study visits from pregnancy week 12 through to delivery. Three time points during pregnancy were chosen for blood chemistry data because of their previously published trajectory of change: week 12 (starting point), week 28 (elbow point), and week 36 (ending point) (21). The child’s VDBP genotype was included as 6 levels (1s, 1f or 2 homozygous or 1s/1f, 1f/2 or 1s/2 heterozygous).
2.2 Bayesian Variable Selection
In our Bayesian paradigm, we used the Gibbs variable selection (GSV) with an indicator (π) to index the probability of including its corresponding predictor (33).
The linear function was: where Xi is the design matrix and βi is the coefficient for the ith predictor. The prior distribution of indicator was constructed in a hierarchical way: A uniform prior for the selection probability was chosen with α = β = 1.
For the binary outcomes, EH status, OP status, and PEB status, we applied the logistic regression with random effects as the following: For the counted outcomes, EH extent, OP extent, and PEB extent, we applied the truncated Poisson regression with random effects as the following: where 7 represents the upper limit of the count, and Bi is the random effect that allows for a random intercept for each subject under a normal prior distribution with mean 0 and variance τ2.
Following the variable selection, the final model was fitted again using only the selected predictors to obtain the corresponding coefficients.
2.3 Variable Imputation and Multinomial Regression
For missing values and outliers in some predictors, we assumed a prior distribution for each to impute values automatically when model fitting. The VDBP type was a multi-level categorical predictor (6 levels) and with some missing values. For the missing values we considered multinomial regression to impute its values with mother’s race as its predictor, which is based on previous studies about the association between race and genotype towards the VDBP types: where a and b are the coefficients in the multinomial regression.
2.4 Fractional Polynomial Selection
After variable selection, we considered non-linear patterns by considering fractional polynomial models. The candidate power set was the commonly used: In particular, the power equals 0 means the logarithm transformation. We did the selection of power by comparing deviance information criterion (DIC), and the model with the smallest DIC value was chosen.
3. Results
3.1 Demographic information
Major characteristics of the 161 mothers who entered this study are presented in Table 1 by their treatment group assignment from the vitamin D supplementation RCT. The final mother-child dyad cohort sizes were 145 for EH, 144 for OP and 143 for PEB because of missingness in the children’s dental measurements, and also each dental defect had a different missing pattern.
3.2 Bayesian variable selection
To assess the importance of predictors in relation to the DDE outcomes, we employed Gibbs variable selection (33). The models with indicators greater than 0.5 were the optimal predictive models (34). We applied a slightly looser criterion (threshold at 0.3 in the first step and 0.4 in the following steps) to choose variables. We did this because the maternal blood chemistry predictors from 3 time points were used in our model and their mutual correlation could impact the indicator’s significance due to their mutual correlation. Variable selection results are shown in Table 3.
The results for status models and extent score models for each defect are similar but not the same, except for the PEB outcome. We found the maternal FVDD was associated with both EH and OP outcomes, while the mother’s BMI was associated with EH and PEB outcomes, and the child’s gestational age week was only selected for the OP outcome.
As seen in Table 3, different time points (maternal weeks 12 and/or 28) for blood chemistries were selected for the models of the 3 outcomes. Moreover, for the OP models, the child’s FVDD and gestational age were selected as well as the maternal FVDD at week 12.
3.3 Effect of predictors
3.3.1 Enamel Hyperplasia Defects
For EH status, the 95% credible interval for all the coefficients included 0, and the same results for EH extent model, but with a narrower interval width (see Appendix A2). The mean values of the coefficients were smaller than 0; namely, these predictors were negatively associated with the EH outcome, even though not dramatically under the 0.05 criteria.
3.3.2 Opacity Defects
For OP status, the coefficient of gestational week had a 95% credible interval above 0, suggesting a positive association with the outcome, i.e., a longer gestation time was associated with a higher risk of having opacity defects in the child’s teeth. Holding others as constant, 1 unit increase in the gestation week increased the odds ratio as high as B0.36 = 1.43. However, the other predictors did not show an association with OP defects measured by binary status.
The OP extent score showed a strong association with gestational week based on the 95% credible interval. Holding others as constant, 1 unit increase in the gestation week increased the mean count (λ) as high as B0.4 = 1.49 time the baseline.
The child’s FVDD was negatively associated with OP extent score; however, the association only marginally approached the 0.05 level (Appendix A3), resulting in some overlaps between groups in Fig 2 (a). Note the ratio of 25(OH)D and iPTH below or equal to 0.308 indicated FVDD, which was then denoted as the FVDD group. We can see the FVDD group was associated with a higher risk of OP defects in child’s teeth. Holding others as constant, the FVDD group increased the average OP score (λ) as high as B1.1 = 3.00 times compared to the non-FVDD group.
3.3.3 Post Eruptive Breakdown Defects
For the PEB status, the coefficient of mother’s BMI was strongly and negatively associated with the outcome under a 95% credible interval. This indicates that a mother’s higher BMI value was associated with a lower risk of having PEB defects in child’s teeth. Holding others as constant, 1 unit increase from 21.5 in BMI decreased the odds ratio as much as e−2.66×(log(22.5)−log(21.5)) = 0.89.
The PEB score also showed a strong association with mother’s BMI. Holding others as constant, 1 unit increase from 21.5 in BMI decreased the average score (λ) as much as e−1.43×(log(22.5)−log(21.5)) = 0.93.
Moreover, the maternal serum P level was strongly and positively associated with both PEB status and extent score, according to the 95% credible interval. Holding others as constant, 1 unit increase in phosphorus increased the odds ratio as much as e1.02 = 2.77 for PEB status and increased the average score (λ) as much as 1.57. (Appendix A4).
Comparing status and extent models for each specific defect, the predictors showed the same direction towards the outcome on average, if they shared the same predictor set. Also, the coefficients for the extent models commonly had a shorter 95% credible interval width.
3.4 Fractional Polynomial Models
For EH status models, the main effects selected were maternal BMI and maternal Ca at week 28. We tested the fractional polynomial (FP) powers on these two predictors. However, the minimum DIC was 202.8 (Appendix A5), which is trivially improved compared with the linear model (DIC = 205.4). Thus, we kept the original linear trend to consider interactions. For EH extent, the main effects selected were maternal BMI and maternal FVDD at week 28. The latter one was a categorical predictor, so we tested fractional polynomial on the maternal BMI only. Like the status model, DIC values indicated the linear trend is a good fit and we kept it due to simplicity. Similarly, the other four models - OP status, OP extent, PEB status and PEB extent - showed the non-linear terms were not necessary. Though non-linear trends were found by mean trajectory plots, e.g., OP status versus mother’s BMI, we kept the linear trend, for simplicity and explanation consideration.
3.5 Models with Interaction terms
We checked the need to add interaction terms after the fractional polynomial selection. For all 6 parallel models we added possible interaction terms with the main effects fixed in the models. The coefficients of interaction terms all had a 95% credible interval including 0 (see Appendix A6). The DIC values for these models were slightly larger or the same compared with the models without interaction effects. Thus, we did not need to build the models with interaction terms.
4. Discussion
Our Bayesian model selection identified maternal and child predictors for the outcome of developmental defects of enamel (EH, OP and PEB) in the primary maxillary central incisor teeth. This unique study contributes knowledge to the etiology of these three major localized defects by having potential predictors that were periodically measured during the specific period of that tooth development. The key predictors for DDE were the mother’s pre-pregnancy BMI, her serum Ca and P levels at 28 weeks of pregnancy, and her FVDD condition at 12 and 28 weeks of pregnancy; and the child’s gestational age in weeks and FVDD at birth. The models for EH and PEB shared the predictor of mother’s pre-pregnancy BMI; whereas the OP model did not include those predictors. The mother’s FVDD at 12 weeks of pregnancy was important for OP and the mother’s FVDD at 28 weeks was important for EH.
New to this study was the predictor FVDD, determined as the ratio of 25(OH)D / iPTH ≤ 0.308 (20) at specific time points, namely at pregnancy weeks 12 and 28, and for the child at birth. This use of a FVDD grouping helped to reduce the number of predictors and also presented a clearer pattern that iPTH dominated the effect of 25(OH)D towards the defects, compared to our previous study that used 25(OH)D and iPTH levels only as independent variables (21). In future studies we will further assess FVDD predictability when used as a continuous variable over the weeks of pregnancy.
A strength of this observational cohort design using secondary analyses of existing data was the inclusion of maternal and child biomarkers of the key components of enamel formation e.g., Ca, P, 25(OH)D, and iPTH with biological plausibility for relation to enamel formation. Most all categories of previously identified high risk factors associated with these dental defects were included. We assessed early infancy feeding habits as exclusively breast milk or formula fed by 4-6 weeks. Our sample included African American, Caucasian, and Hispanic mothers; and the child’s genotype for VDBP. The self-selected sample for participation through to the dental imaging was from a population of healthy mothers that excluded any conditions related to vitamin D regulation (19). And yet in this healthy sample we found DDE rates of 41.4% for EH, 31.2% for OP and 48.3% for PEB. Our defect rates are not comparable to most published data because our mothers and children were very healthy and the outcome measures were restricted to only 2 teeth (PMCI), rather than the full primary dentition. Also, there is the potential confounding in the mothers because of participation in an RCT of vitamin D supplementation. However, the mothers’ characteristics by treatment group were similar and also our key variables were from circulating blood chemistries rather than vitamin D intake.
In consideration of sample size and the impact of outliers, we found that the risk of OP changed slowly when the gestational age was less than 36 weeks, and then the risk for OP increased dramatically beyond 36 weeks. Thus, we are more likely to state the effect of gestational age as positive when the value of gestational week is relatively very high. Similar ideas related to the effect of outliers may also explain the impact of the mother’s BMI and serum chemistry predictors. One explanation is that the range of our predictors were concentrated within a certain interval, making their extreme values to have a greater leverage effect on the results. This may be due to a relatively limited sample size, as well as limited heterogeneity in the healthy cohort.
Our methodology and results provide a roadmap for assessing time appropriate exposure measures during specific tooth development to better understand the etiology of developmental defects of enamel. Though the PMCI are present at birth with the structural abnormalities of EH and OP, these defects are not recognized until after the teeth erupt into the oral cavity. Our ongoing studies continue to examine the interrelationships among the three defects and their exposome to better understand the in utero and early infancy development of these dental defects in the primary maxillary central incisor teeth. Our efforts continue towards improving the understanding of the development of enamel defects and possibly modifiable variables for the future prevention of DDE.
Data Availability
All data produced in the present work are available at the completion of the ongoing research and upon reasonable request to the authors.