Education and income-based inequality in tooth loss among Brazilian adults: does the place you live make a difference?

doi:10.21203/rs.2.18559/v2

Background: Socioeconomic inequalities in tooth loss might be minimized or potentialized by the characteristics of the context where people live. We examined whether there is contextual variation in socioeconomic inequalities in tooth loss across Brazilian municipalities.

Methods: Data from the 2010 National Oral Health Survey of 9,633 adults living in 157 Brazilian municipalities were used. Education and household income were used as the individual socioeconomic indicators. At the municipal level, we used the Municipal Human Development Index (HDI) as our contextual indicator of socioeconomic status (low:<0.699 versus high: >0.70). The Relative (RII) and Slope (SII) Indexes of Inequality, Relative (RCI), and Absolute (ACI) Concentration Indexes were calculated to compare the magnitude of education and income-based inequalities among municipalities with low versus high HDI. Multilevel Poisson regression models with random intercepts and slopes were developed. Results: At the individual level, adults with lower education & income reported more tooth loss. The mean number of missing teeth was 9.52 (95%CI: 7.93-11.13) and 6.95 (95%CI: 6.43-7.49) in municipalities with low and high HDI, respectively. Municipalities with high HDI showed higher relative and absolute education-based inequality. For income-based inequalities, higher SII and RCI was observed in municipalities with lower HDI. A significant cross-level interaction indicated that high-education adults reported fewer missing teeth when they lived in municipalities with high HDI compared to adults with the same education level living in low HDI municipalities. For individuals with the lowest education level, there was no difference in the number of teeth between those from municipalities with high and low HDI.

Conclusions: There was a social gradient in tooth loss by education and income. Living in disadvantaged municipalities cannot overcome the risk associated with low schooling. The protective effect of higher education can be reduced when people live in disadvantaged areas.

Dentistry

Epidemiology

Tooth loss

Socioeconomic Factors

Health Status Disparities

Adult

Multilevel Analysis

Tooth loss is an important oral health problem with consequences for physical and psychosocial health as well as the quality of life [1-3]. It is a good proxy for the cumulative oral health status of individuals, [4] summarizing the impacts of adverse circumstances throughout the life course [5].

Despite a significant decline in the prevalence and incidence of severe tooth loss in the past two decades [6], socioeconomic inequalities in this condition persist across the globe, including Brazil [4, 7-9].

The socioeconomic inequalities in tooth loss have been studied considering compositional (individual) and contextual effects based mainly on the pathways linking income inequality to health [10-14]. Metanalysis research with 11 studies showed that adults with lower levels of income presented greater chance of tooth loss [13]. Considering contextual-level, collective evidence has shown associations between high income-based inequality and worse oral health [15]. Compositional explanations suggest that the effect of income inequality is a result of how individual incomes affect oral health. The proposed mechanisms involve the material resources (eg. the ability to purchase higher quality diets, to afford preventive and regular dental care due to the treatment cost), access to symbolic resources (status and rank within one’s community) and behavioral/cultural explanations, which stress the role of poor health behaviors due to low income (tobacco use, high sugar consumption, infrequent and symptomatic dental visits and poor oral hygiene practices) [14, 16]. Contextual explanations sustain that income inequality has broader psychosocial effects on oral health via stress-induced oral-health-related behaviors [17]. Income inequality also may result from lower social spending on public services and infrastructure, including dental care services and water fluoridation [18].

The income inequality had been a frequent contextual social determinant of oral health in the epidemiologic studies. However, it may not capture all dimensions in which social inequality can occur. Different levels of education can also explain such disparities. Education is a socioeconomic indicator that captures the long-term influences of both early life circumstances on adult health, as well as the influence of adult resources (for example, through employment status) on health [19]. Previous studies have shown a consistent independent effect of lower individual education on tooth loss among adults [20-23] and elders [24]. Education equips individuals with knowledge and skills that are useful for the prevention of disease and attitudes toward healthy behavior. Also, higher educational attainment confers greater prestige and status within the community as well as serving as a credential for employment [14].

Beyond the interest in distinguishing the individual (compositional) source of variation from the contextual on socioeconomic inequalities studies, it is crucial to evaluate whether contextual factors affect every income and education group in the same manner [25]. Hence, it is reasonable to assess whether the variations among areas are similar for different income or education groups or whether the variations differ for different education and income groups. This approach is crucial from an equality policy perspective to target for interventions from specific groups living in specific areas. Additionally, it takes into account that socioeconomic diversity may minimize or potentialize the adverse oral health consequences of those living in an area with concentrated disadvantage [26].

The influence of contextual determinants on oral health after accounting for individual-level factors (compositional effect) had been extensively shown assuming the fixed effects [18, 20-24, 27-33]. In general, these studies had shown higher tooth loss in socially disadvantaged areas independently of individual socioeconomic indicators [22, 24, 32-34]. The evidence also showed attenuation on income inequalities among American adults and Brazilian teenagers by public policies, including dental care services and water fluoridation [18, 35].

A fewer number of studies assessed the interaction between oral health and contextual level factors, and most of them were focused on income inequality [20, 23, 31, 35, 36]. Differential effects of public policies on oral health among individuals with higher education and income had been observed [28, 35]. These results have conducted to discuss if the public policy may both decrease or increase the social gap in oral health, considering the aim, the target population, and the approach of implementation of them (e.g. prioritization of municipalities with better socioeconomic). Can the disadvantage resulting from low education and income be decreased in areas that present more opportunities and conditions for healthy living? Or Do public policies have more effect on groups that are already socially advantaged? The double disadvantage, characterized by concomitantly presenting a worse social situation in the individual and contextual level, has been reported on oral epidemiology[37, 38].

For that, it is crucial to investigate the income and education-based inequality considering the social characteristics of context level. A composite indicator of social development can reflect the contextual differences as a result of public policy. The starting point was to estimate the magnitude of income and education-based inequality according to areas with different social development. The aims of this study were therefore: i) to evaluate income and education-based inequalities in tooth loss comparing Brazilian municipalities with high and low Human Development Index (HDI); ii) to investigate the association between tooth loss and socioeconomic indicators on individual (income and education) and contextual (HDI) levels; iii) to determine if the social environment modifies the relationship between income/education and tooth loss. We hypothesized that tooth loss would be inversely associated with an individual’s household income and education, but that the magnitude of the association would be more pronounced in affluent municipalities (high HDI) than in poor municipalities (low HDI).

Data from the present analysis came from the 2010 National Oral Health Survey (SB Brazil 2010) conducted by the Brazilian Ministry of Health in Brazilian urban areas [39] between February and November 2010. The sample was obtained through the random selection of municipalities and census sectors, via multi-stage cluster sampling with probability of selection proportional to population size. Detailed information on the methods is found in other publications [40, 41]. Data for adults aged between 35 and 44 years were used in this study.

Individual interviews using a structured questionnaire was used to obtain demographic and socioeconomic characteristics. Oral health examinations were conducted in people’s homes by calibrated dentists (kappa>0.65) under natural light following the guidelines of the WHO manual for epidemiological studies [42]. The DMFT (Decayed, Missing, and Filled Teeth) index was used to determine tooth status.

Outcome variable

The outcome variable was the number of teeth (discrete quantitative variable) that was missing for any reason, determined by the sum of codes 4 and 5 of the DMFT index.

Exposures

Education and income were used as measures of socioeconomic position at the individual level. Income was measured as total income received by all family members in the month preceding the survey (in seven categories from “R$250.00 or less” to “R$9500.00 or more”). For our analysis, the monthly household income was converted into multiples of the minimum wage, based on the current value at the time of survey (1 minimum wage = R$ 510.00, USD$303.57) and collapsed into four categories: up to 1, 1 to 2.9, 3 to 4.9, 5 or more times the minimum wage. This grouping was defined to distingue mechanisms in which income exerts its effects. Education was asked as the number of years of formal schooling and classified as less than four years (insufficient education), 4-7 (incomplete elementary education), 8-10 (completed elementary, but incomplete secondary education), and 11 or more (completed secondary, incomplete university education, or college graduate), according to the formal education system in Brazil.

We used the Municipal Human Development Index (HDI) as an indicator of municipal area SES. The Brazilian HDI considers three dimensions: Longevity, Education (access to knowledge, based on average years studied of the population) and Income (living standards and purchasing power of the population according to the Municipal Gross Income per capita) [43]. The HDI was obtained from the 2013 Brazil Atlas of Human Development, which allows a selection based on data extracted from 2010 demographic. The values of HDI of sampled municipalities ranged from 0.481 to 0.847, and the groups are defined as very low (0-0.499), low (0.500 – 0.599), medium (0.600 – 0.699), high (0.700 – 0.799) and very high (0.800 – 1.00). According to this classification, the frequency of sampled municipalities was: low and very low (3.93%), medium (7.88), and high and very high (88.19). In this study municipalities were aggregated into low (<0.699) versus high (>0.70)

Covariates

The covariables at individual level were age (adult: 35-39 and 40-45 years old), sex (female, male), skin color/ethnicity (white, black, yellow and brown/Ameridians), time since the last dental visit (< 12 months, between 1 and 2 years, > 3 years, did not visit). Skin color refers to the classification adopted in the demographic census performed in Brazil: whites, blacks, browns, yellows, and Amerindians. In this study, this variable was dichotomized: white versus blacks, browns, yellows, and Amerindians.

At the contextual level, we included the presence of fluoridated water supply (present or absent). The data regarding fluoridation was obtained on the National Basic Sanitation Survey performed by the Brazilian Institute of Geography and Statistics (IBGE) in 2008 [44]. We also included the estimated coverage of the population by primary care oral health services, which corresponds to the mean monthly number of primary care oral health teams for every 3000 individuals to the total population of the municipality in the analyzed year. Higher oral health services coverage indicated higher potential access to basic dental services. The cut off for this variable was 40% that was the goal to be achieved in the biennium 2010/2011[45]. Data about coverage were obtained from the website of the Department of Information Technology of the Unified Health System (DATASUS).

Statistical analysis

Comparison of income and education-based inequalities between municipalities with high and low HDI

Descriptive analysis was performed to obtain mean tooth loss for each municipality, and the results were shown separately according to the HDI level (high or low). The magnitude of relative and absolute educational and income-based inequalities in the tooth loss was calculated using the Relative Index of Inequality (RII), Slope Index of Inequality (SII) and Relative (RCI) and Absolute Concentration Index (ACI) for municipalities with high and low HDI. RII and SII are summary measures recommended when making comparisons across populations [46]. These indices are regression-based and take the whole socioeconomic distribution into account, rather than only comparing the two most extreme groups. For municipalities with high and low HDI, the population in each education or income category was assigned a modified ridit-score based on the midpoint of the range in the cumulative distribution of the participants in the given category. We used generalized linear models (log-binomial regression), with a logarithmic link function to calculate RIIs (rate ratios) and with an identity link function to calculate SIIs (rate differences) [47]. Both indices were estimated with 95% confidence intervals. The RII can be interpreted as the rate ratio and the SII can be interpreted as the rate difference at the bottom and the top of the educational or income hierarchy. If there is no inequality, RII assumes the value of 1.0. The further the value of RII from 1.0, the higher the level of inequality. RII assumes only positive values, with values larger than one indicating a concentration of the indicator among the advantaged and values smaller than one indicating a concentration of the indicator among the disadvantaged. If there is no inequality, SII takes the value of zero. Greater absolute values indicate higher levels of inequality. Positive values indicate higher coverage in the advantaged subgroups, and negative values indicate higher coverage in the disadvantaged subgroups. RCI and ACI were estimated based on the methodology described by Clarke et al. [48]. The RCI provides a measure of the relative differences among education/ income groups. Considering that the RCI was based on morbidity measure (tooth loss), a negative RCI indicates inequality favoring higher income groups. The ACI quantified the absolute differences in health between education and income groups, and this index is not affected by whether it is measured concerning health or morbidity. As RII/SII are mathematically related to RCI/ACI, they will produce the same rank ordering of health inequality between groups but will differ in scale.

Association between tooth loss and socioeconomic indicators

Multilevel Poisson regression procedures with unstructured covariances matrix were used to model the two-level structure of individuals (level 1) nested within municipalities (level 2). Multilevel techniques of analyses provide the overall relationship between the individual, compositional factors, and oral health (fixed part), and the variation between areas that cannot be accountable for such factors (random-intercept parameter). Besides, it is possible to assess the variation in certain individual relationships between municipalities (random slope parameters) and the interaction between individual and contextual characteristics (cross-level interactions). A five-step sequential modeling strategy was adopted: i) model 1(empty model): a model without the inclusion of any covariates, in which the variance in tooth loss is inspected between municipalities. A significant random intercept variance indicates the presence of unexplained differences in tooth loss between municipalities. The Wald test evaluated the significance of random intercept, and the Median Rate Ratio measured the heterogeneity among municipalities, according to Austin et al. (2017) [49]. There is no variation between municipalities if the MRR is 1.0, but the higher the MRR, the greater the area-level variation. ii) model 2 (random intercept, fixed effect): considers all the individual-level variables in the fixed part. This model assessed the association between tooth loss and income and education adjusted for covariables at the individual level. The variation between municipalities is allowed for, conditional on the individual, compositional factors. iii) model 3: as model 2, but including the municipalities-level variables. The proportional change in variance (PCV) was calculated according to Merlo et al. (2005) [50], using the following formula: PCV = (variance model 1 – variance model 2)/variance model 1. iv) model 4 (random slope and random intercept model): as model 2, but the model allows both the intercept and slope to be random parameters. Therefore, each municipality has its intercept and slope in which the variability from the overall intercept and slope can be investigated with the addition of individual and municipality variables and their interactions (cross-level interactions). Random slopes for education and income were considered. The comparison of goodness fit between model 2 and model 4 was performed using the LR-test. v) model 5: as model 4, including the municipalities-level variables, in which the cross-level interactions between education/income and HDI were considered. The interaction term represents the change in the slope of education/income on tooth loss across municipalities when HDI changed from low to high. The estimated mean of missing teeth according to individual socioeconomic variables for municipalities with high and low HDI was demonstrated using a graph of the predicted model.

Procedures for complex sample design were used to calculate the proportions of adults according to investigated variables, the age- and sex-adjusted estimates of the outcomes, SII, RII and ACI. Statistical analyses were performed using STATA version 15.0 (StataCorp LP, College Station, Texas, USA).

The Brazilian National Council of Ethics in Research approved the SBBrasil 2010 study, protocol no. 15498, January 7, 2010.

The sample was 9633 adults, residing in 157 municipalities; 146 individuals in 20 municipalities with a sample of fewer than 10 adults were excluded. The characterization of the sample and mean of missing teeth according to investigated variables are shown in table 1.

Table 1: Proportion of Brazilian adults and mean of tooth loss according to individual and municipalities level variables. Brazil. 2010.

	Sample size	Prop. 95%CI*	Tooth loss* Mean (95% CI)
Individual level variables
Sex
Male	3,317	36.74	6.59(6.01, 7.16)
Female	6,316	63.26	7.81(7.17, 8.44)
Age
35-40	5,061	52.39	5.74(5.18, 6.30)
40-45	4,572	47.61	9.13(8.42, 9.84)
Household income (in the minimum wage)
up to 1	1,388	12.96	10.21(9.13, 11.30)
1 to 2.9	4,699	52.90	8.31(7.72, 8.91)
3 to 4.9	1,823	20.39	5.60(4.80, 6.41)
> = 5 mw	1,480	13.75	3.63(2.92, 4.34)
Education (in years of schooling)
0 to 4	840	10.05	12.09(11.77, 13.41)
5 to 8	2,459	29.41	9.34(8.57, 10.11)
9 to 11	3,966	37.84	6.70(6.14, 7.27)
> 12	2,284	22.70	3.72(3.08, 4.35)
Self-reported skin color/ ethnicity
Black + yellow + blown + Ameridians	5,558	50.24	7.93(7.35, 8.52)
White	4,075	49.76	6.77(6.14, 7.42)
Time since last dental visit
Never used	672	7.32	10.13(8.28, 11.98)
< 1 year	4,488	46.69	6.40(5.81, 6.70)
1 to 2 years	2,409	25.14	7.39(6.72, 8.06)
> 3 years	1,882	20.85	8.39(7.56, 9.11)
Municipalities level variables
Human Development Index
Low + Medium	69	23.72	9.52 (7.93, 11.13)
High	88	76.28	6.95 (6.42, 7.49)
Coverage of oral health services
Below of goal	78	60.82	6.96 (6.40,7.52)
Above of goal	79	39.18	8.31(7.20,9.41)
Fluoridation of water supply
No	52	17.44	9.44(8.48,10.42)
Yes	105	82.56	7.06(6.52,7.59)

95% confidence intervals (CI) in brackets. *Estimates considered weighting and complex sampling design.

The mean number of adults per municipalities was 61.36 and ranged from 10 to 496. The mean of HDI was 0.71 (range: 0.481-0.847) and 76.28% presented high HDI. The goal of oral health services coverage was achieved by 39.18 of the municipalities and 82.56 of those had fluoridation of water supply. The mean number of missing teeth among Brazilian adults was 7.36 (SE: 0.26; 95% CI: 6.86-7.87). This value was 9.52 (SE: 0.82; 95% CI: 7.93-11.13) in municipalities with low HDI and 6.95 (SE: 0.27; 95% CI: 6.43-7.49) in those with high HDI (Figure 1).

The inequalities indexes show a higher number of missing teeth in disadvantaged groups (lowest education and income levels). Higher relative and absolute education-based inequality were observed among municipalities with high HDI. For income-based inequalities, higher SII was observed in municipalities with lower HDI (Table 2).

Table 2: Education and income-based inequalities in municipalities with low and high HDI

Education-based inequality			Income-based inequality
	RII	RCI		RII	RCI
Brazil	0.27 (0.22;0.34)	-0.19		0.35 (0.28; 0.43)	-0.14
HDI Low	0.41 (0.27; 0.63)	-0.14	Low	0.42 (0.28; 0.62)	-0.12
HDI High	0.27 (0.21; 0.34)	-0.20	High	0.37 (0.29; 0.46)	-0.13
	SII	ACI		SII	ACI
Brazil	-8.87 (-10.14; -7.60)	1.42		-7.63 (-8.92; -6.35)	1.06
HDI Low	-8.09 (-11.16; -5.01)	1.31	Low	-8.47 (-11.88; -5.06)	1.16
HDI High	-8.54 (-9.91; -7.17)	1.37	High	-6.91 (-8.28; -5.53)	0.95

*SII: Slope Index of Inequalities and RII: Relative Index of Inequalities; ACI: Absolute Concentration Index; RCI: Relative Concentration Index.

The multilevel crude estimates showed that all individuals variables were significantly associated with tooth loss (Table 3).

Table 3: Crude Count Ratio (95% confidence intervals in brackets) of tooth loss according to multilevel models among 34-44 year-olds in Brazil, 2010.

	Crude Count Ratio (95% CI)
Individual level variables
Sex
Male	1
Female	1.18 (1.16,1.20)
Age group
35-39 years old	1
40-45 years old	1.59 (1.56,1.61)
Household income (in minimum wage)
up to 1	1
1 to 2.9	0.88 (0.86,0.89)
3 to 4.9	0.68 (0.66,0.70)
> = 5 mw	0.46 (0.45,0.47)
Education (in years of study)
0 to 4	1.00
5 to 8	0.89 (0.87,0.92)
9 to 11	0.65 (0.64,0.67)
> 12	0.41 (0.40,0.42)
Self-reported skin color/ ethnicity
White	1
Black + yellow+ brown + Ameridians)	1.16 (1.15,1.18)
Time since last dental visit
Never used	1
< 1 year	0.75 (0.73,0.77)
1 to 2 years	0.86 (0.83,0.89)
> 3 years	1.01 (0.98,1.05)
Municipalities level variables
Municipal Human Development Index
Low + Medium
High	0.71 (0.62,0.81)
Fluoridation of water supply
No	1
Yes	0.77 (0.66,0.89)
Coverage of public oral health services
Below of goal	1
Above of goal	1.18 (1.03,1.36)

Table 4 shows the multilevel models assuming random intercept and fixed effects. The empty model shows that the contextual level variation was significant, suggesting differences in tooth loss among municipalities (LR test: chi2 = 7171.68; Prob >= chibar2 = <0.0001), and the MRR (1.53) also indicated variation across municipalities. The MRR, as well as the between-municipality variance, come down from 20.09 to 14.17 (-28.46%) after the inclusion of individual variables (compositional differences) (Model 2) (Table 4). In model 3, assuming fixed effects, a social gradient in tooth loss for income and education was observed. Municipalities with high HDI presented a lower number of tooth loss (Count Ratio: 0.86, 95% CI: 0.76-0.98). The inclusion of municipalities level variables represented an additional 12.7% decrease in variability of tooth loss among municipalities (Table 4).

Table 4: Count ratios (95% confidence intervals in brackets) of tooth loss in multilevel models with random intercept and fixed effect between individual and municipalities level variables among 35–44-year old in Brazil, 2010.

Parameters	Empty model (Model 1)	Model 2*	Model 3**
Fixed part	Crude Count Ratio (95% CI)	Adjusted Count Ratio (95% CI)	Adjusted Count Ratio (95% CI)
Individual level variables
Constant	8.62 (8.02, 9.26)	8.83 (8.22, 9.48)	10.29 (9.09,11.65)
Education (in years of study)
0 to 4		1	1
5 to 8		0.92 (0.90,0.94)	0.92 (0.90,0.94)
9 to 11		0.72 (0.70,0.74)	0.72 (0.70,0.74)
> 12		0.50 (0.49,0.52)	0.50 (0.49,0.52)
Income (in minimum wages)
Up to 1		1	1
1 to 2.9		0.96 (0.94,0.98)	0.96 (0.94,0.98)
3 to 4.9		0.84 (0.82,0.86)	0.85 (0.82,0.86)
> 5		0.65 (0.63,0.67)	0.65 (0.63,0.68)
Municipalities level variables
Municipal Human Development Index
Low + Medium			1
High			0.86 (0.76,0.98)
Random part
Area level variance (Random intercept)	20.09 (15.88,25.41)	14.17 (11.15,18.02)	12.37 (9.71,15.76)
PCV^&		-29.46%	-12.7%
Median Rate Ratio	MRR = 1.53	MRR = 1.43	MRR = 1.39

Results of Multilevel Poisson Regression Model assuming random intercept and fixed effect. Exponentiated coefficients; ^&PCV: Proportional change in variance. *Model 2: Adjusted for individual-level variables: sex, age group, skin color, and time since the last dental visit. **Model 3: Adjusted for individual and municipalities level variables: sex, age group, skin color, time since the last dental visit, presence of fluoridated water supply, and coverage of public oral health service above the Brazilian goal.

Table 5 shows the models whereby the nature of between-municipalities variation was modeled as a function of individual education (random slope). The Likelihood-ratio Test showed that model 4 presented a better fit than model 2 (Likelihood-ratio test: 1204.74, Prob > chi2 = 0.0000), justifying this modeling step. The high-education groups not only have a lower number of tooth loss but also more variance (variance: 0.575 (0.405,0.815) compared to groups with 5-8 years of study (variance: 0.268 (0.271,0.494) and 9-11 years of study (variance: 0.367 (0.272,0.493) (Table 5). The other random slopes were shown as an additional file (Additional file 1). The graph of fitted lines shows how the municipal variation is more pronounced among individuals with higher education (Additional file 2).

Model 5 shows that when the association between education and tooth loss was allowed to vary across municipalities, the main effect of HDI was not significant. The social gradient in tooth loss by individual income and education continue to be observed. Municipalities with fluoridated water showed a lower number of missing teeth (Table 5).

Table 5: Count ratios (95% confidence intervals in brackets) of tooth loss in multilevel models with random intercept and slope (education) between individual and municipalities level variables among 35–44-year old in Brazil, 2010.

	Random intercept and slope (education)
Parameters	Model 3* Adjusted Count ratio (95% CI)	Model 4* Adjusted Count ratio (95% CI)
Fixed part
Individual level variables
Constant	9.18 (8.42,10.00)	9.81 (8.53, 11.28)
Education (in years of study)
0 to 4	1	1
5 to 8	0.82 (0.75,0.91)	0.82 (0.75,0.91)
9 to 11	0.64 (0.57,0.71)	0.64 (0.57,0.71)
> 12	0.44 (0.38,0.50)	0.43 (0.38,0.50)
Income (in minimum wages)
Up to 1	1	1
1 to 2.9	0.94 (0.92,0.96)	0.94 (0.92,0.96)
3 to 4.9	0.83 (0.80,0.85)	0.83 (0.80,0.85)
> 5	0.68 (0.65,0.70)	0.68 (0.65,0.70)
Municipalities level variables
Municipal Human Development Index
Low + Medium
High	-	0.95 (0.85,1.22)
Fluoridation of water supply
No
Yes	-	0.88 (0.78,0.99)
Coverage of public oral health services
Below of goal
Above of goal	-	1.09 (0.98,1.22)
Median Rate Ratio	1.49	1.48

Exponentiated coefficients; *Model 3 and 4: Adjusted for individual level variables: sex, age group, skin color and time since the last dental visit.

The cross-level interaction between education and HDI was significant, indicating that the association between these two variables is not constant across municipalities, with steeper slopes observed in municipalities with high HDI. There is some evidence that the number of missing teeth is lower for the high-education group when they live in municipalities with high HDI compared with those with the same education level living in municipalities with low HDI. For individuals with the lowest education level, there was no difference in the number of teeth comparing those from municipalities with high versus low HDI. The fixed effect of individual and contextual variables did not change with the inclusion of cross-level interactions (Table 6).

When the variation was modeled as a function of individual income, the model was very similar to obtained with fixed effect (Additional file 3). The graph of fitted lines shows that the municipal variation is similar through income levels (Additional file 4). The cross-level interaction between individual income and HDI was not significant (Figure 2).

Table 6: Rate ratios (95% confidence intervals in brackets) of tooth loss in multilevel models with random intercept and slope (education) and cross-level interactions between HDI and education 35-44 years old in Brazil, 2010.

Parameters	.
Fixed part	Adjusted Count Ratio (95% CI)
Individual level variables
Constant	9.18 (8.42,10.00)
Education (in years of study)
0 to 4	1
5-8 years of study	0.87 (0.79,0.95)
9-11 years of study	0.78 (0.67,0.89)
> 12 years of study	0.59 (0.48,0.73)
Income (in minimum wages)
Up to 1
1 to 2.9 minimum wage	0.94 (0.93,0.96)
3 to 4.9 minimum wage	0.83 (0.80,0.85)
> 5 minimum wage	0.67 (0.64,0.69)
Contextual factors (Municipalities level)
Municipal Human Development Index
Low + Medium	1
High	1.04 (0.90-1.19)
Fluoridation of water supply
No
Yes	0.85 (0.75-0.96)
Coverage of public oral health services
Below of goal	1
Above of goal	1.17 (0.96-1.19)
Cross-level interactions
HDI#5-8 years of study	0.99(0.89,1.11)
HDI#9-11 years of study	0.80 (0.66,0.96)
HDI#> 12 years of study	0.67 (0.51,0.88)
Median Rate Ratio (MRR)	1.41

Exponentiated coefficients; *Model was adjusted for individual level variables: sex, age group, skin color and time since the last dental visit.

The tooth loss represents an “end state” in oral health, and the present study reaffirmed the persistent concentration of tooth loss among the most disadvantaged individuals. Over and above individual disparities, our findings also demonstrated that the magnitude of educational and income-based inequalities varied between municipalities with high versus low HDI. Our results showed that educational attainment was “less protective” against tooth loss when people live in disadvantaged areas. Conversely, adults with the lowest levels of education and income had greater tooth loss regardless of whether they lived in a municipality with high or low HDI.

Our study replicated the inverse association between socioeconomic status and tooth loss after adjusting for all covariates at the individual and municipality levels [10-12]. Using different outcomes related to tooth loss (< 9 remaining natural teeth; lack of functional dentition/ inadequate dentition: < 21 natural teeth; count of tooth loss or remaining teeth and edentulousness), previous studies have consistently shown more tooth loss among the more disadvantaged groups. Quasi-experimental evidence that worsened economic circumstances (measured by subjective economic deterioration or housing damage due to disaster damage) are associated with tooth loss was provided by a natural experiment following exposure to the 2011 Great East Japan Earthquake and Tsunami [51]. A meta-analysis of 10 cohort studies and two cross-sectional studies found a significant association between low income and tooth loss [13]. Those living in poverty suffer a greater burden of oral diseases, such as dental caries and periodontitis [52, 53], the main causes of tooth loss in adults. They also have a higher prevalence of systemic conditions, such as diabetes, cardiovascular disease, and obesity [54, 55], that are related to tooth loss. Besides these biological explanations, the conceptual mediators of the effect of SES on oral health include material deprivation, psychosocial distress, and behavioral factors. [14, 17, 19].

In addition to those factors, it has been demonstrated that economic constraints are associated with the type of dental treatment delivered. While subjects in the lower income brackets are more prone to dental extraction, individuals with higher income are more likely to seek periodic routine appointments and conservative dental treatment, resulting in a higher number of retained teeth [56, 57]. In Brazil, the higher tooth loss associated with socioeconomic disadvantage can result from systemic factors related to the organization and delivery of oral health services. The national oral health policy, also known as Smiling Brazil (“Brasil Sorridente”), was implemented only in 2004, ensuring universal coverage [58]. This policy increased the access of adults to oral health services, who historically had significant ongoing unmet needs with only urgent care available in public services [58]. The lack of access to restorative and preventive services in public health services may have contributed to tooth loss, especially among the most disadvantaged individuals. Moreover, previous studies have shown that a change in access to oral health services does not necessarily result in reduced oral health disparities because it takes time for people to change their behaviors, e.g. seeking periodic preventive dental services [59].

Our results also showed differences in tooth loss inequalities according to municipal HDI, with higher absolute and relative education-based inequality found in municipalities with high HDI and higher absolute income-inequality found in municipalities with low HDI. The high education-based inequality among municipalities with high HDI shows the double disadvantage of adults who share the lowest education and living in places with the worst social indicators [38]. These results suggest that public policies may not reach the persons in the same way increasing the gap among social groups.

The association between education and tooth loss was not constant across Brazilian municipalities. The significant cross-interaction term reveals that the magnitude of the association between education and tooth loss was lower in municipalities with high HDI. This result suggests that contextual factors may weaken the potential protective effect of schooling at the individual level. Low HDI municipalities in Brazil are smaller cities with lower wealth, higher rates of urban violence, lower sanitation conditions, and urban infrastructure. These contextual aspects may represent fewer opportunities for healthy living for its residents (less access to health services, less availability of diversified healthy food, fewer leisure options, less culture, and social capital). Our result shows that where one lives has an effect on oral health and that this effect must be evaluated considering the characteristics of the people.

On the other hand, the significant cross-level interactions reinforce the hypothesis that public policies can increase the gap among social groups to the extent that it may not reach people in the same way. The cross-interaction was significant after adjusting by oral health policies (coverage of oral health services) and fluoridation of water supply. DHI is an indicator that incorporating a life expectancy, education, and a modified measure of income in contextual level. However, these indicators may demonstrate the effect of other policies and environmental characteristics that favor access to primary health care in general, sanitation, access to cultural and social opportunities, more social capital, more work and economic security, and more healthy food available. All this contextual factors are potential mediators of the association between HDI and tooth loss and have been pointed out as determinants of oral health [34, 60] in previous studies that had demonstrated the main effect of social-level indicators on tooth loss using multilevel approaches in Brazil, United States, European countries, Australia and Japan [18, 20-24, 27-32].

The analysis performed in this study is advancing, as it may contribute to more targeted policies for the groups with the greatest need, in specific contexts.

The only significant contextual factor we found was the fluoridation of the water supply. Individuals living in municipalities with water fluoridation presented lower tooth loss. The importance of this intervention was previously shown in maintaining functional dentition [21] and tooth retention [23]. Natural experiments applied in oral health context from Brazil showed that adults who accessed fluoridate water <50% of their lifetime presented a higher mean rate ratio of DMFT index compared with those living >75% of their life with residential access to fluoridated water. Longer residential lifetime access to fluoridated water was associated with less dental caries, even in a context of multiple exposures to fluoride [61].

Our models were intentionally parsimonious, adjusting only for key covariates. As with most multilevel studies, our choice of area-unit (municipalities) was made for reasons of sampling and analytic convenience rather than being underpinned by an explicit theory linking area disadvantage and oral health; hence associations among these variables are likely to be underestimated. However, the Brazilian municipalities have to deliver basic education, sanitation, and primary healthcare. Therefore, our chosen contextual factors (HDI, oral health services coverage, and fluoridation of water supply) must fit the area border. We used cross-sectional data where the distinction between current and past exposures cannot be made. Similarly, income can change throughout life resulting in upward and downward social mobility.

There was a social gradient in tooth loss by education and income. Living in disadvantaged municipalities cannot overcome the risk associated with low schooling. The protective effect of higher education can be reduced when people live in disadvantaged areas. This differential effect of education on tooth loss among municipalities with high and low DHI can increase the gap among social groups. The persistent inequality demands for innovative strategies and policies that address the wider social determinants of health, considering the differential effect of public policies among disadvantage and advantage part of society. These policies must protect oral health as a human right through, strengthen intersectoral strategies for poverty reduction, supporting scientific research on the health social determinants, interact communities with public health managers and researchers, support community actions to promote oral health and eliminate barriers of access to oral health care.

HDI - Municipal Human Development Index

RII - Relative Index of Inequality

SII - Slope Index of Inequality

DMFT - Decayed, Missing, and Filled Teeth

IBGE - Brazilian Institute of Geography and Statistics

DATASUS - Department of Information Technology of the Unified Health System

MRR - Median Rate Ratio

PCV - Proportional change in variance

Ethical approval and consent to participate: The Brazilian National Human Research Ethics Committee approved the 2010 NOHS under process number 15,498 on July 1st, 2010. All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Informed consent was obtained from all individual participants included in this study.

Consent for publication: Not applicable

Availability of data and materials: The datasets analyzed during the current study are available in the Ministério da Saúde under request.

Competing Interests: The authors RCF is a member of the editorial board of the BMC Oral Health.

Funding: Raquel Conceição Ferreira received financial support CAPES (PVEX 88881.172365/2018-01) and from FAPEMIG, Brazil (Fundação de Amparo à Pesquisa do Estado de Minas Gerais - PPM-00603-18). The Universidade Federal de Minas Gerais. This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brasil (CAPES) – Finance Code 001. These sponsors had no role in study design; collection, analysis, and interpretation of data; writing of the report; the decision to submit the article for publication.

Authors’contributions: RCF analyzed and interpreted the data, wrote the manuscript. MIBS, LGR, FLC and AEBLM contribute to data analysis and interpretation. IK was the study supervisor, made substantial contributions to: data analysis and interpretation and made critical review of the intellectual content. All authors read and approved the final manuscript.

Acknowledgments: The authors would like to thank “Pró-reitoria de Pesquisa da Universidade Federal de Minas Gerais” for financially supporting the article-processing charge.

Dioguardi M, Gioia GD, Caloro GA, Capocasale G, Zhurakivska K, Troiano G et al. The Association between tooth loss and Alzheimer's Disease: a systematic review with meta-analysis of case control studies. Dent J. 2019;7:E49.
Koka S, Gupta A. Association between missing tooth count and mortality: A systematic review. J Prosthodont Res. 2018;62:134-51.
Haag DG, Peres KG, Balasubramanian M, Brennan DS. Oral conditions and health-related quality of life: a systematic review. J Dent Res. 2017;96:864-74.
Cunha-Cruz J, Hujoel PP, Nadanovsky P. Secular trends in socio-economic disparities in edentulism: USA, 1972-2001. J Dent Res. 2007;86:131-6.
Han DH, Khang YH. Lifecourse socioeconomic position indicators and tooth loss in Korean adults. Community Dent Oral Epidemiol. 2017;45:74-83.
Kassebaum NJ, Bernabe E, Dahiya M, Bhandari B, Murray CJ, Marcenes W. Global burden of severe tooth loss: a systematic review and meta-analysis. J Dent Res. 2014;93:20S-8S.
Bernabe E, Sheiham A. Tooth loss in the United Kingdom--trends in social inequalities: an age-period-and-cohort analysis. PLoS One. 2014;9:e104808.
Elani HW, Harper S, Thomson WM, Espinoza IL, Mejia GC, Ju X et al. Social inequalities in tooth loss: A multinational comparison. Community Dent Oral Epidemiol. 2017;45:266-74.
Peres MA, Barbato PR, Reis SCGB, Freitas CHSdM, Antunes JLF. Tooth loss in Brazil: Analysis of the 2010 Brazilian oral health survey. Revista de Saude Publica. 2013;47:78-89.
Barbato PR, Muller Nagano HC, Zanchet FN, Boing AF, Peres MA. [Tooth loss and associated socioeconomic, demographic, and dental-care factors in Brazilian adults: an analysis of the Brazilian Oral Health Survey, 2002-2003]. Cadernos de saude publica. 2007;23:1803-14.
Jimenez M, Dietrich T, Shih MC, Li Y, Joshipura KJ. Racial/ethnic variations in associations between socioeconomic factors and tooth loss. Community Dent Oral Epidemiol. 2009;37:267-75.
Sabbah W, Tsakos G, Sheiham A, Watt RG. The role of health-related behaviors in the socioeconomic disparities in oral health. Soc Sci Med. 2009;68:298-303.
Seerig LM, Nascimento GG, Peres MA, Horta BL, Demarco FF. Tooth loss in adults and income: Systematic review and meta-analysis. J Dent. 2015;43:1051-9.
Kawachi I, Adler NE, Dow WH. Money, schooling, and health: Mechanisms and causal evidence. Ann N Y Acad Sci. 2010;1186:56-68.
Singh A, Harford J, Schuch HS, Watt RG, Peres MA. Theoretical basis and explanation for the relationship between area-level social inequalities and population oral health outcomes - A scoping review. SSM - population health. 2016;2:451-62.
Singh A, Peres MA, Watt RG. The Relationship between Income and Oral Health: A Critical Review. J Dent Res. 2019:22034519849557.
Bernabe E, Watt RG, Sheiham A, Suominen-Taipale AL, Uutela A, Vehkalahti MM et al. Sense of coherence and oral health in dentate adults: findings from the Finnish Health 2000 survey. Journal of clinical periodontology. 2010;37:981-7.
Bernabe E, Marcenes W. Income inequality and tooth loss in the United States. J Dent Res. 2011;90:724-9.
Galobardes B, Shaw M, Lawlor DA, Lynch JW, Davey Smith G. Indicators of socioeconomic position (part 1). Journal of epidemiology and community health. 2006;60:7-12.
Sanders AE, Turrell G, Slade GD. Affluent neighborhoods reduce excess risk of tooth loss among the poor. J Dent Res. 2008;87:969-73.
Koltermann AP, Giordani JM, Pattussi MP. The association between individual and contextual factors and functional dentition status among adults in Rio Grande do Sul State, Brazil: a multilevel study. Cadernos de saude publica. 2011;27:173-82.
Goulart Mde A, Vettore MV. Is the relative increase in income inequality related to tooth loss in middle-aged adults? J Public Health Dent. 2016;76:65-75.
Barbato PR, Peres MA, Hofelmann DA, Peres KG. Contextual and individual indicators associated with the presence of teeth in adults. Rev Saude Publica. 2015;49:27.
Ito K, Aida J, Yamamoto T, Ohtsuka R, Nakade M, Suzuki K et al. Individual- and community-level social gradients of edentulousness. BMC oral health. 2015;15:34.
Subramania SV, Kawachi I, Kennedy BP. Does the state you live in make a difference? Multilevel analysis of self-rated health in the US. Soc Sci Med. 2001;53:9-19.
Castelli A, Jacobs R, Goddard M, Smith PC. Health, policy and geography: insights from a multi-level modelling approach. Soc Sci Med. 2013;92:61-73.
Turrell G, Sanders AE, Slade GD, Spencer AJ, Marcenes W. The independent contribution of neighborhood disadvantage and individual-level socioeconomic position to self-reported oral health: a multilevel analysis. Community Dent Oral Epidemiol. 2007;35:195-206.
Singh A, Harford J, Antunes JLF, Peres MA. Area-level income inequality and oral health among Australian adults-A population-based multilevel study. PLoS One. 2018;13:e0191438.
Moreira Rda S, Nico LS, Barrozo LV, Pereira JC. Tooth loss in Brazilian middle-aged adults: multilevel effects. Acta odontologica Scandinavica. 2010;68:269-77.
Chalub LL, Martins CC, Ferreira RC, Vargas AM. Functional Dentition in Brazilian Adults: An Investigation of Social Determinants of Health (SDH) Using a Multilevel Approach. PLoS One. 2016;11:e0148859.
Celeste RK, Nadanovsky P, Ponce de Leon A, Fritzell J. The individual and contextual pathways between oral health and income inequality in Brazilian adolescents and adults. Social science & medicine (1982). 2009;69:1468-75.
Koyama S, Aida J, Saito M, Kondo N, Sato Y, Matsuyama Y et al. Community social capital and tooth loss in Japanese older people: a longitudinal cohort study. BMJ open. 2016;6:e010768.
Aida J, Hanibuchi T, Nakade M, Hirai H, Osaka K, Kondo K. The different effects of vertical social capital and horizontal social capital on dental status: a multilevel analysis. Soc Sci Med. 2009;69:512-8.
Barbato PR, Peres KG. Contextual socioeconomic determinants of tooth loss in adults and elderly: a systematic review. Revista brasileira de epidemiologia = Brazilian journal of epidemiology. 2015;18:357-71.
Celeste RK, Nadanovsky P. How much of the income inequality effect can be explained by public policy? Evidence from oral health in Brazil. Health Policy. 2010;97:250-8.
Aguiar VR, Pattussi MP, Celeste RK. The role of municipal public policies in oral health socioeconomic inequalities in Brazil: A multilevel study. Community Dent Oral Epidemiol. 2018;46:245-50.
Vettore MV, Rebelo Vieira JM, J FFG, Martins NMO, Freitas YNL, Lamarca GA et al. Individual- and City-Level Socioeconomic Factors and Tooth Loss among Elderly People: A Cross-Level Multilevel Analysis. International journal of environmental research and public health. 2020;17.
Schulz M, Kunst AE, Brockmann H. High educational attainment moderates the association between dental health-care supply and utilization in Europe. European journal of oral sciences. 2016;124:52-61.
Brasil. SB Brasil 2010: Pesquisa Nacional de Saúde Bucal: resultados principais In:SB Brasil 2010: Pesquisa Nacional de Saúde Bucal: resultados principais 2011. [http://bvsms.saude.gov.br/bvs/publicacoes/pesquisa_nacional_saude_bucal.pdf]
Silva NN, Roncalli, A. G. . Sampling plan, weighting process and design effects of the Brazilian Oral Health Survey. Revista de Saúde Pública. 2013;47:3-11.
Roncalli AG, Silva NN, Nascimento AC, Freitas CHS, Casotti E, Peres KG et al. Aspectos metodológicos do Projeto SBBrasil 2010 de interesse para inquéritos nacionais de saúde. Cadernos de saude publica. 2012;28:S40-S57.
WHO: Oral health surveys: basic methods, 4 ed. Geneva: ORH/EPID; 1997.
WHO: Atlas do Desenvolvimento Humano no Brasil In.: World Health Organization; 2012.
IBGE. Pesquisa Nacional de Saneamento Básico [National Basic Sanitation Survey]. In:Pesquisa Nacional de Saneamento Básico [National Basic Sanitation Survey]. 2008. [https://www.ibge.gov.br/estatisticas/multidominio/9073-pesquisa-nacional-de-saneamento-basico.html]
Brasil. Portaria 3840. In:Portaria 3840. 2010. [http://bvsms.saude.gov.br/bvs/saudelegis/gm/2010/prt3840_07_12_2010.html]
Harper S, Lynch J: Methods for measuring cancer disparities: using data relevant to healthy people 2010 cancer-related objectives. Montreal: Center for Social Epidemiology and Population Health, University of Michigan; 2010.
Spiegelman D, Hertzmark E. Easy SAS calculations for risk or prevalence ratios and differences. American journal of epidemiology. 2005;162:199-200.
Clarke PM, Gerdtham UG, Johannesson M, Bingefors K, Smith L. On the measurement of relative and absolute income-related health inequality. Social science & medicine (1982). 2002;55:1923-8.
Austin PC, Stryhn H, Leckie G, Merlo J. Measures of clustering and heterogeneity in multilevel Poisson regression analyses of rates/count data. Stat Med. 2018;37:572-89.
Merlo J, Yang M, Chaix B, Lynch J, Rastam L. A brief conceptual tutorial on multilevel analysis in social epidemiology: investigating contextual phenomena in different groups of people. J Epidemiol Community Health. 2005;59:729-36.
Matsuyama Y, Aida J, Tsuboya T, Hikichi H, Kondo K, Kawachi I et al. Are Lowered Socioeconomic Circumstances Causally Related to Tooth Loss? A Natural Experiment Involving the 2011 Great East Japan Earthquake. American journal of epidemiology. 2017;186:54-62.
Costa SM, Martins CC, Pinto MQC, Vasconcelos M, Abreu M. Socioeconomic factors and caries in people between 19 and 60 Years of age: an update of a systematic review and meta-analysis of observational studies. International journal of environmental research and public health. 2018;15.
Borrell LN, Crawford ND. Socioeconomic position indicators and periodontitis: examining the evidence. Periodontology 2000. 2012;58:69-83.
Claassen MA, Klein O, Bratanova B, Claes N, Corneille O. A systematic review of psychosocial explanations for the relationship between socioeconomic status and body mass index. Appetite. 2019;132:208-21.
Polzer I, Schwahn C, Volzke H, Mundt T, Biffar R. The association of tooth loss with all-cause and circulatory mortality. Is there a benefit of replaced teeth? A systematic review and meta-analysis. Clinical oral investigations. 2012;16:333-51.
Klock KS. Patients' perceptions of the decision-making process leading to extraction of permanent teeth in Norway. Community Dent Oral Epidemiol. 1995;23:165-9.
Thomson WM, Poulton R, Kruger E, Boyd D. Socio-economic and behavioural risk factors for tooth loss from age 18 to 26 among participants in the Dunedin Multidisciplinary Health and Development Study. Caries research. 2000;34:361-6.
Pucca GA, Jr., Gabriel M, de Araujo ME, de Almeida FC. Ten Years of a National Oral Health Policy in Brazil: Innovation, Boldness, and Numerous Challenges. J Dent Res. 2015;94:1333-7.
Pereira CR, Roncalli AG, Cangussu MC, Noro LR, Patricio AA, Lima KC. [Impact of the Family Health Strategy: an analysis in cities in Northeast Brazil with more than 100,000 inhabitants]. Cadernos de saude publica. 2012;28:449-62.
Singh A, Harford J, Peres MA. Investigating societal determinants of oral health-Opportunities and challenges in multilevel studies. Community Dent Oral Epidemiol. 2018;46:317-27.
Peres MA, Peres KG, Barbato PR, Hofelmann DA. Access to Fluoridated Water and Adult Dental Caries: A Natural Experiment. J Dent Res. 2016;95:868-74.
Valente MIB, Vettore MV. Contextual and individual determinants of periodontal disease: Multilevel analysis based on Andersen's model. Community Dent Oral Epidemiol. 2018;46:161-8.

Additional file 1: Parameters of the random part of the multilevel models with random intercept and slope (education) between individual and municipalities level variables among 35–44-year old in Brazil, 2010.

Additional file 2: Municipalities-specific slopes of education on tooth loss

Additional file 3: Count ratios (95% confidence intervals in brackets) of tooth loss in multilevel models with random intercept and slope (income) between individual and municipalities level variables among 35–44-year old in Brazil, 2010.

Additional file 4: Municipalities-specific slopes of income on tooth loss

Education and income-based inequality in tooth loss among Brazilian adults: does the place you live make a difference?

Status:

Journal Publication

Version 2

Abstract

Figures

Background

Methods

Results

Discussion

Conclusions

List Of Abbreviations

Declarations

References

Additional File Legends

Supplementary Files

Status:

Journal Publication

Version 2