1. Introduction
The controlled fluoridation of drinking water has been implemented in many countries worldwide for the prevention of dental caries [
1,
2]. While studies show that fluoridation provides benefits regarding oral health, ‘anti-fluoridation’ opinions have been voiced in reference to its putative negative health consequences, such as increased risks for autism, certain cancers, or reduced intelligence [
3]. However, a comprehensive review of laboratory-based mechanistic and epidemiological studies suggested that potential health risks related to fluoride exposure may be due to biases of epidemiological studies with low methodological quality [
4].
With respect to blood pressure, a recent systematic review of epidemiological studies indicated a positive association between fluoride exposure and elevated blood pressure among individuals from endemic fluorosis areas in Iran, Turkey and China [
5], but not in other study populations. Again, the quality of the included study was graded as low because objective biomarkers measurements of fluoride status were lacking and potential confounders were not accounted for in the majority of studies [
5]. The only identified study among children was a biomarker study in a well-characterised cross-sectional sample of children aged 10 to 18 years from Mexico City, a non-fluoride endemic area. This study showed that higher plasma fluoride levels were related to higher systolic and diastolic blood pressure levels, and to higher BMI, trunk fat, insulin, and fasting-glucose levels, albeit only among girls [
6].
Here, we used data from the 2013/2014 and 2015/2016 cycles of the National Health and Nutrition Examination Survey (NHANES) to evaluate the association between objectively measured fluoride status, i.e., plasma fluoride levels, and systolic and diastolic blood pressure among 3495 children and adolescents aged 8 to 19 years. Considering the above-mentioned finding of a worse pattern of cardio-metabolic parameters among girls from Mexico reported by Liu et al. [
6], we further investigated whether higher plasma fluoride was associated with BMI, waist circumference, blood lipids (total, HDL, and non-HDL cholesterol), HbA1c, glucose, insulin, and C-reactive protein levels.
2. Materials and Methods
2.1. Study Population
The National Health and Nutrition Examination Survey (NHANES) is an ongoing cross-sectional survey carried out in the USA. Since the 1990s, the study provides detailed health status monitoring in large, nationally representative samples of non-institutionalised individuals aged 1 to 80 years, from survey cycles carried out every two years. Assessments include detailed medical interviews, physical examinations, and collections of urine and blood samples for comprehensive biomarker measurements. An overview of the protocols used for different examinations and interview assessments can be found on the NHANES website [
7]. All NHANES examinations and the public release of the study’s data were approved by the Research Ethics Review Board (ERB) of the US National Center for Health Statistics (NCHS). Children aged 7–17 years provided informed assent, and parents of these participants provided informed consent. Participants aged 18 years and older provided informed consent.
For the present analyses on fluoride status and parameters of cardio-metabolic health, data from the NHANES cycles 2013–2014 and 2015–2016 were used as plasma fluoride was measured among study participants aged 8 to 19 years in these cycles. Out of the 5451 children and adolescents aged 8 to 19 years, who took part in NHANES between 2013 and 2016, fluoride status was measured among 4470, and three blood pressure measurements were obtained for 3495 of these study participants. Thus, the analytical sample for the present analyses consisted of 3495 individuals. Multiple imputation by Fully Conditional Specification was carried out in case of sporadic missing values for covariates (ratio of family income to poverty: n = 59; parental education level: n = 116; plasma cotinine: n = 146; waist circumference and body mass index: n = 25) accounting for the survey design.
The sample sizes were smaller due to missing blood biomarker measurements for analyses on fluoride status and secondary endpoints (blood lipids: n = 3443; glucose and insulin: n = 2173; HbA1c: n = 2212; and C-reactive protein: n = 1823).
2.2. Fluoride and Blood Pressure Measurements
Details on the methods for fluoride measurements can be found on the NHANES website [
8,
9]. In short, fluoride levels were quantified in tap water and plasma samples using an ion-specific electrode. The lower limits of detection for fluoride in water and plasma were 0.10 mg/L and 0.25 µmol/L. Inter-assay coefficients of variation based on repeated measurements of fluoride concentrations in low, medium, and high-quality control samples were below 5% across both NHANES cycles in plasma and water. Details on the measurements of anthropometric parameters and further plasma biomarkers (blood lipids, biomarkers of sugar metabolism, and C-reactive protein) can be found online on the NHANES website [
7]. The NHANES staff took three consecutive blood pressure measurements, after allowing the participants to rest in a seated position for five minutes. The first measurement was not included in the final calculation of average systolic and diastolic blood pressure as it was taken to allow for any external factors influencing the measurement to be removed. Further details on the equipment and protocols used for the blood pressure measurements are available online [
10].
2.3. Statistical Methods
Survey-weighted means (standard errors) and survey-weighted frequencies were obtained for continuous and categorical covariates, respectively, across sex-specific quartiles of plasma fluoride levels in order to describe the study population. The age- and sex-adjusted Spearman’s coefficient was calculated to evaluate the correlation between plasma and water fluoride levels. In addition, this correlation was visualised by a bubble plot.
Associations between blood pressure parameters, secondary endpoints (BMI, waist circumference, biomarkers of lipid and sugar metabolism, and C-reactive protein) and sex-specific quartiles of plasma fluoride levels were evaluated by survey-weighted linear regression models, using log2-transformed endpoint parameters as continuous dependent variables. In addition, plasma fluoride levels were modelled on the log2 scale to test for linear trends in associations with blood pressure parameters and secondary endpoints.
Regression models were first adjusted for age and sex (Model 1). Additional multivariable adjustment was carried out for the second multivariable linear regression model for standing height (cm), waist circumference (cm), BMI (kg/m2), household reference person’s education level (less than 9th grade, 9–11th grade, high school graduate, some college or AA degree and college graduate or above), cotinine levels (ng/mL) as a biomarker of exposure to second-hand smoke, ethnicity (Mexican American, Other Hispanic, Non-Hispanic White, Non-Hispanic Black, Non-Hispanic Asian, and Other Race including multi-racial), ratio of family income to poverty, and fasting duration (minutes). Confounders were selected by literature search.
As the above-mentioned study among Mexican children had shown sex differences in associations between fluoride levels and blood pressure parameters as well as anthropometric indices and blood lipids [
6], we carried out regression analyses stratified by sex in addition to our main analyses. Associations were considered as statistically significant at two-sided
p-values < 0.05. Statistical analyses were carried out using SPSS 26.0.0.0 (IBM, Armonk, NY, USA).
3. Results
3.1. Characteristics of the Study Population
Characteristics of the study population are shown in . Mean ages were 14.0 ± 0.1, 13.4 ± 0.1, 13.4 ± 0.1 and 13.8 ± 0.2 years across quartiles of plasma fluoride. The proportion of Non-Hispanic White participants was slightly greater in the highest fluoride quartile (60.2% vs. 50.1% in the lowest quartile), whereas the proportion of Mexican Americans showed an opposite trend (12.8% vs. 18.4%). The percentage of household reference persons with higher education level (28.6% with at least a college degree) was higher in fluoride quartile 4 compared to quartile 1 (22.7%).
Table 1. Characteristics of the study population across quartiles of plasma fluoride (n = 3495) *.
While most other characteristics showed similar proportions and mean values across quartiles of plasma fluoride, there were tendencies for higher triglyceride (Q4: 102.7 ± 4.3 vs. Q1: 94.2 ± 3.1 mg/dL), water fluoride (Q4: 0.67 ± 0.04 vs. Q1: 0.30 ± 0.03 mg/L) and cotinine levels (Q4: 12.7 ± 2.5 vs. Q1: 4.5 ± 1.0 ng/L) with higher fluoride levels. Fasting duration was slightly shorter in the highest compared to the lowest quartile (Q4: 362.1 ± 13.8 min vs. Q1: 399.2 ± 13.9). There was a statistically significant positive correlation between water fluoride content and plasma fluoride concentrations (Spearman’s r = 0.41, p < 0.01, also visualised in ).
3.2. Associations between Fluoride Status and Cardio-Metabolic Biomarkers
Adjusted mean values of blood pressure across quartiles of plasma fluoride are shown in . There were no statistically significant associations between fluoride levels and systolic or diastolic blood pressure, neither in the age- and sex-adjusted regression model nor in the multivariable-adjusted one.
Table 2. Adjusted mean levels of blood pressure and metabolic biomarkers across sex-specific quartiles of plasma fluoride (n = 3495).
Results from regression analyses on fluoride status and further cardiometabolic endpoint parameters are shown in . With regard to BMI and waist circumference, age- and sex-adjusted regression analyses did not show significant associations. In the multivariable-adjusted model, there were non-significant trends for positive associations with both parameters, with p-values for linear trend of 0.06 and 0.07, respectively. The only parameter for which a significant linear association with fluoride was observed in both regression models was HbA1c, with marginally lower values in the highest quartile (5.2% vs. 5.3% in the lowest quartile). Results from ancillary analyses stratified by sex are presented in . The above-mentioned inverse association between fluoride and HbA1c levels was statistically significant among female participants (p linear trend = 0.007) but not among male participants (p linear trend = 0.11). Moreover, there was a borderline significant inverse association between fluoride and HDL–cholesterol levels among female study participants (p linear trend = 0.04), whereas HDL was not associated with fluoride among male participants.
4. Discussion
In the present study, fluoride exposure objectively measured by plasma fluoride levels was not associated with systolic and diastolic blood pressure in a large study among 8- to 19-year-old children and adolescents from the USA. A statistically significant inverse association was observed between plasma fluoride levels and HbA1c values, although this association was weak in magnitude. Plasma fluoride was not associated with other biomarkers of cardio-metabolic state.
One motivation for the present analyses was a recent study by Liu et al., in which higher plasma fluoride was significantly associated with blood pressure among girls from Mexico City [
6]. The reason for the discrepancy in the results by Liu et al. and those from the present study remain unclear, as average plasma fluoride concentrations in the study from Mexico were at 0.21 µmol/L, i.e., lower than in the present one (0.40 µmol/L). This difference in plasma fluoride is consistent with the fact that community water fluoridation is not common in Mexico, unlike in the USA, where around 73% of the population consume fluoride-enriched water [
6,
11]. Additionally, it should be noted that the significant associations in the study by Liu et al. were only evident among girls but not among boys [
6], while this analysis showed no indication for a sex difference. Thus, further studies are needed to provide more clarity on potential associations between fluoride status and blood pressure among children and adolescents. Among adults, fluoride exposure has been associated with increased blood pressure in fluoride endemic areas [
5,
12], while biomarker studies from other areas are lacking and are needed before conclusions on a potential blood pressure increasing effect of moderate fluoride exposure can be drawn.
Unlike Liu et al., who reported higher glucose and insulin levels among girls with higher plasma fluoride levels [
6], we did not observe significant associations between fluoride status and these biomarkers. By contrast, there was a moderate inverse association between fluoride and HbA1c among female participants of the present study. We further found an inverse association between plasma fluoride and HDL-cholesterol, while no such association was detected by Liu et al. [
6]. The present study and the study by Liu et al. were consistent in that they showed borderline significant tendencies for positive associations between higher fluoride levels, BMI, and waist circumference, although it should be noted that both studies were cross-sectional. Therefore, reverse causality, i.e., higher fluoride intake via food and water among obese children, cannot be ruled out, and prospective studies are needed to investigate a potential obesogenic effect of fluoride exposure.
While this study indicates that plasma fluoride levels are not related to increased cardio-metabolic risks among children and adolescents, we acknowledge that analyses from the same NHANES dataset suggest associations between fluoride and increased levels of biomarkers of kidney and liver function [
13] and sex steroid hormones [
14]. Moreover, higher plasma fluoride levels were associated with sleep disturbances in NHANES [
15], which is in line with findings from a recent cross-sectional study among Canadian adults [
16]. Again, prospective studies are needed to better assess temporal relationships between fluoride exposure and these health outcomes, while the mentioned studies do underline the necessity for an integrated risk-benefit assessment of fluoride exposure taking into account a broad range of health indicators. Interestingly, we found a tendency for higher cotinine levels among study participants with higher plasma fluoride, which is in line with previous studies that suggest that tobacco exposure may increase fluoride levels [
17,
18]. Thus, accounting for smoking as a potential confounder seems important for future analyses of potential adverse health effects of fluoride.
The following methodological strengths and limitations need to be considered when interpreting the present findings. NHANES is representative for the US general population, and importantly, objective biomarkers of fluoride exposure were measured. The sample size facilitated well-powered statistical analyses, and data on important potential confounders such as smoking exposure were available. However, the cross-sectional design of the study needs to be acknowledged. A potential further limitation of this research project is the focus on two NHANES cycles (2013/2014 and 2015/2016), while biomarker studies from other areas of the world are lacking and needed. Finally, the statistically significant associations between fluoride and HbA1c and HDL should be interpreted with caution, given the moderate differences in mean biomarker concentrations across quartiles of fluoride levels and the high number of statistical tests.