Mediterranean diet and risk of rheumatoid arthritis: a population-based case-control study

Study population

This study included subjects from the Swedish epidemiological investigation of RA (EIRA), a population-based case-control study, from defined geographical areas of central Sweden. The EIRA study was initiated in 1996 and is still ongoing. The general study design has been described in detail elsewhere [15]. Briefly, the study includes patients with incident RA (cases) (mean time between symptom onset and diagnosis 9.7 months) diagnosed by a rheumatologist according to the either the American College of Rheumatology (ACR) 1987 criteria or the 2010 ACR/European League Against Rheumatism (EULAR) classification criteria for RA. To each patient two controls were randomly selected from the register of the general population, in connection to inclusion of the patient, matched on age, sex and residential area.

A baseline questionnaire was distributed to patients and controls with comprehensive questions on lifestyle factors, educational level and comorbidities. From November 2005 the EIRA questionnaire also included a 124-item food frequency questionnaire (FFQ). The current study therefore included participants from November 2005 to September 2014.Participants with incomplete an FFQ where excluded (n = 21).

Participants provided written informed consent, and ethical approval was obtained from the Regional Ethical Review Board at Karolinska Institutet, Stockholm, Sweden.

Exposure definition

The self-reported FFQ was used to evaluate the participants’ food intake during the last year before inclusion. In the FFQ participants were asked to indicate how often on average they had consumed various foods by using eight predefined frequency categories ranging from “never or less than once per month” to “three or more times per day”. Fat and energy intake was calculated by multiplying the average frequency of consumption of each food by the fat and energy content, respectively, of age-specific portion sizes, according to the Swedish National Food Administration Database [16].

The Mediterranean diet was defined according to the Mediterranean diet score based on the work of Trichopoulou et al. [17]. The score range between 0 and 9 and each person was assigned 1 point if she/he consumed more than the sex-specific median of the population of the beneficial components, such as vegetables, legumes, fruits and nuts, cereal and fish, or if she/he consumed less than the median of red meat and dairy products, which were considered detrimental components. Moreover, 1 point was assigned to men who consumed between 10 g and 50 g of alcohol per day and to women who consumed between 5 g and 25 g per day, and 1 point if the ratio of monounsaturated lipids to saturated lipids was above the median. Hence, a higher score corresponded to higher adherence to the Mediterranean diet. We did not include butter among the dairy products since it was reported in the FFQ only as “ever consumption”. In sensitivity analysis, we considered ever consumption of butter as consumption once per day and added it to the dairy products group.

Modifications to the original Mediterranean score have been proposed throughout the years [18, 19]. In sensitivity analyses, we defined the Mediterranean diet according to the alternate Mediterranean diet score (aMed), which not include dairy products in the score [13]. Moreover, we used a lower upper limit for alcohol intake among men (30 g instead of 50 g per day), according to the modified Mediterranean diet (mMED) score [19].

Covariate definition

All covariate data were collected at baseline. Smoking was categorised as never, current, former or irregular smoker. Body mass index (BMI) was categorised as underweight (< 18.5 kg/m²), normal weight (18.5–24.9 kg/m²), overweight (25–29.9 kg/m²) or obese (≥ 30 kg/m²). Education was categorised as <10 years, 10–12 years, and > 12 years of education. Physical activity was reported as sedentary, moderate, moderate/regular, exercise/regular. Nutritional supplements were defined as intake in the last 5 years of at least one of the following supplements: multivitamins (with or without minerals), vitamin C, vitamin E, vitamin B, vitamin A, vitamin D, calcium, zinc, iron, magnesium, folic acid, β-caroten or omega 3. Additional covariates were considered when analysing women only: parity (quartiles), age at menarche (quartiles), use of oral contraceptive (yes/no) and use of hormone replacement therapy (yes/no). Genotyping of the HLA-DRB1 gene was conducted as previously described [20, 21]. Among HLA-DRB1 genes, DRB1*01, DRB1*04 and DRB1*10 genes were defined as SE alleles. Any genotype containing one or two of these genes was considered as having “any SE allele” versus those not having any of the genes (“no SE alleles”).

Statistical analyses

The Mediterranean diet score was analysed both as a continuous and as a categorical variable, the latter according to the quartiles of the distribution. Odds ratios (ORs) and their 95% confidence intervals (CIs) were calculated using conditional logistic regression, and multivariable models were adjusted for smoking, BMI, educational level, physical activity, use of dietary supplements, and energy intake. Analyses were stratified by gender, rheumatoid factor (RF) and anti-citrullinated protein antibody (ACPA) status. We additionally stratified by smoking status and presence of at least one copy of the HLA-DRB1 shared epitope allele. In sensitivity analyses, we excluded persons with extreme energy intake (i.e. 3 standard errors from the mean value on the log-transformed scale, n = 49). Moreover, to evaluate the influence of each single item, we analysed data excluding one component of the score at a time.

We estimated the dose–response trend association between the Mediterranean diet score and risk of RA using restricted cubic splines with knots at 10, 50 and 90 centiles of the Mediterranean score (e.g. at 2, 4, and 7) [22].

To analyse the association in a group of patients presumed to have more homogeneous disease aetiology, we restricted the analyses to patients and controls with at least one copy of the HLA-DRB1 shared epitope allele (SE) and who were also current smokers. For this analysis, since patients and controls were no longer matched, we performed logistic regression adjusted for age, gender and residential area, and then additionally adjusted for the confounding variables listed above. Statistical analyses were implemented using SAS (V.9.4) and Stata (V.13.1), and p values ≤0.05 were considered significant.

Patient involvement

There was no patient involvement in the current study.

Study characteristics

Mediterranean diet and risk of RA

Fewer patients with RA than controls had high adherence (score 6–9) to the Mediterranean diet (patients 24.1%, controls 28.2%). This pattern was similar among men and women, but women were more often in the highest category of the score (among patients, women 25.8%, and men 19.4%; among controls, women 29.1% and men 25.9%).

The Mediterranean diet was associated with a decrease in the odds of developing RA (Fig. 1). One unit increase in the score corresponded to an 8% decrease in the odds of having RA (crude OR conditioned on matching factors = 0.92, 95% CI 0.89–0.95). After additional adjustment for smoking, physical activity, education level, BMI, energy intake and use of supplements, the OR was 0.95 (95% CI 0.92–0.99). The restricted cubic splines model showed a dose–response trend with no departure from linearity (i.e. the second coefficient of the spline was not statistically significant).

When stratified by gender, the inverse association between high adherence to the Mediterranean diet and RA was statistically significant among men (OR = 0.49, 95% CI 0.33–0.73), but not among women (OR = 0.94, 95% CI 0.74–1.18). Adding further covariates for women (parity, use of oral contraceptive, use of hormone replacement therapy, age at menarche), did not change the results. A multiplicative interaction between high adherence to the Mediterranean diet and gender was found (p < 0.01).

When stratified by RA subtype, the inverse association between Mediterranean diet score and RA was observed only among seropositive patients. A high score corresponded to lower odds in RF-positive but not RF-negative patients and in ACPA-positive but not ACPA-negative patients. When stratifying by RA subtype and gender, we found that a high Mediterranean diet score was associated with a non-significantly reduced risk of RF-positive RA among women (adjusted OR (ORadj) = 0.81, 95% CI 0.61–1.07), and was similar in ACPA-positive RA (ORadj = 0.85, 95% CI 0.64–1.13), while among men there was a reduced risk in both RF-positive (ORadj = 0.44, 95% CI 0.26–0.74) and ACPA-positive RA (ORadj = 0.46, 95% CI 0.28–0.76).

Similar results, although statistically significant, were found when analysing smokers separately irrespectively of HLA-DRB1 SE allele status (adjOR = 0.62, 95% CI 0.40–0.95, Additional file 1: Table S5). Also, analysing only those with the HLA-DRB1 SE allele similar results were retrieved (adjOR = 0.60, 95% CI 0.44–0.82; Additional file 1: Table S6).

Sensitivity analysis

We excluded one component of the score at a time. The inverse association with RA observed for the highest category of the score ranged between 0.74 (95% CI 0.60–0.91) after exclusion of the fat-ratio group to 0.85 (95% CI 0.70–1.04) after exclusion of alcohol.

Alternative definitions of the diet score had little impact on the observed association. Inclusion of butter did not change the results (OR for a score of 6–9 vs. 0–2 = 0.79, 95% CI 0.64–0.96). When persons with extreme energy intake (i.e. 3 standard errors from the mean value on the log-transformed scale, n = 49) were excluded, the results did not change (OR for a score of 6–9 vs. 0–2 = 0.79, 95% CI 0.65–0.96). We considered a different definition of moderate alcohol consumption for men: we lowered the upper limit of 50 g per day to 30 g per day. The results did not change overall (OR = 0.80, 95% CI 0.66–0.97) nor for men specifically (OR = 0.51, 95% CI 0.35–0.75). Moreover, we applied the aMed score, resulting in an association between the Mediterranean diet and RA that was no longer statistically significant, either overall (OR = 0.83, 95% CI 0.63–1.09) or in men (OR = 0.78, 95% CI 0.43–1.41).

Main findings

In this large population-based case-control study, we observed that the Mediterranean diet score was inversely associated with risk of RA. However, this association was observed only among men and among individuals with seropositive RA, and could not readily be explained by other risk factors. A tendency toward lower risk of RA was found among smokers who carried the HLA-DRB1 shared epitope allele and who had a high Mediterranean diet score.

Previous research

Previous studies of diet and risk of RA have been focused foremost on individual components, and the main dietary compounds that have been linked to a potentially decreased risk of RA include moderate alcohol intake and [2, 3] and consumptions of fish [4, 5], both of which are components of the Mediterranean diet. Currently, only one prospective cohort study from the NHS (pooled data from NHS I and II) has specifically examined the association between the Mediterranean diet and risk of RA [13]. This study included 174,638 female nurses and identified 913 incident cases of RA over 20 years of follow up [13]. Similar to our result among women, that study did not find any overall association between the Mediterranean diet and risk of RA. A nested case-control study in the Swedish Västerbotten intervention programme (a prospective population-based study [14]) investigated the effect of diet, looking at food groups, macronutrients and dietary pattern, including the Mediterranean diet. This study comprised 386 incident RA cases and 1886 matched controls and did not find an association between the Mediterranean diet and RA. Similar to our current study, the majority of the participants were women (70%), but no subgroup analysis by sex was performed.

Sundström et al. [14] used the Mediterranean-diet score of Trichopoulou et al. [17] to define the Mediterranean diet, as we did in the current study, while Hu et al. [13] used a modified version of the score (aMed score) [18]. The main differences between the original Mediterranean diet score [17] and the modified version [18] are the exclusion of dairy products from the aMed score and the different definition of moderate alcohol intake. In the original score moderate alcohol consumption was defined as 10–50 g per day for men and 5–25 g per day for women, whereas in the aMed it is defined as 10–15 g per day for men and 5–15 g per day for women. Indeed, when we applied the aMed score in our study, the association was no longer significant. Similarly, when we excluded one item of the score at a time, the association was no longer statistically significant when alcohol was excluded from the analyses, thus implying an influence of alcohol consumption in the prevention of RA, as shown elsewhere [3, 23].

Strengths and limitations

Although the EIRA study is one of the largest population-based case-control studies including incident cases of RA, with extensive information, making it possible to account for many potential confounders, there are limitations of this study. First, the exposure was self-reported via a single-time-administered FFQ assessing dietary intake 1 year prior to RA diagnosis (at inclusion for controls). Although the FFQ used in this and other studies has proven to have high validity [24–26], any case-control study based on self-report may be prone to recall bias. Reverse causation may also be present as patients with early symptoms of RA might already have changed their diet during the year before diagnosis. However, the FFQ was administered close to onset of the first symptom of RA (mean duration 9.7 months). Moreover, a recent Swedish study suggested that that patients newly diagnosed with RA did not make any major dietary modifications after their diagnosis [27]. Another limitation is that the Mediterranean diet score (as well as other ranking scores) does not reflect the absolute Mediterranean diet, but the distribution of the score 0–9 in the studied population. This means that two persons with exactly the same nutritional intake but from two different study populations might end up with different scores. Moreover, even within the same study populations, the scores for women and men are calculated differently. Despite the extensive details of variables and potential confounders included in the EIRA study, there still might be residual confounding.

Mechanisms

Although the Mediterranean diet has been shown to be effective in reducing overall mortality, cardiovascular disease and cancer [9–11], the biological mechanism is not fully understood. Potential mechanisms include lowering of blood pressure, lipids and inflammatory markers, via bioactive compounds that are rich in the Mediterranean diet, including antioxidants, monounsaturated and omega-3 fatty acids, phytosterols and fibre [28, 29].

In the current study, we found an inverse association between the Mediterranean diet and risk of RA but only among men, and only in seropositive RA. A difference in the effect of Mediterranean diet by gender has been shown before for cardiometabolic variables [30, 31], suggesting that indeed men could be more likely to be beneficially influenced by the Mediterranean diet. Moreover, differences in the aetiology of seropositive and seronegative RA have been shown before, in association with both genetic [32] and environmental risk factors [23, 33].