- Open Access
- Open Peer Review
Dioxin-like compounds and bone quality in Cree women of Eastern James Bay (Canada): a cross-sectional study
Environmental Health volume 12, Article number: 54 (2013)
Aboriginal populations living in Canada’s northern regions are exposed to a number of persistent organic pollutants through their traditional diet which includes substantial amounts of predator fish species. Exposure to dioxin-like compounds (DLCs) can cause a variety of toxic effects including adverse effects on bone tissue. This descriptive cross-sectional study was conducted to investigate the relationship between plasma concentrations of DLCs and bone quality parameters in Cree women of Eastern James Bay (Canada).
Two hundred and forty-nine Cree women from seven communities in Eastern James Bay (Canada), aged 35 to 74 years old, participated in the study. In order to determine the total DLC concentration in plasma samples of participants, we measured the aryl hydrocarbon receptor-mediated transcriptional activity elicited by plasma sample extracts using a luciferase reporter gene assay. Plasma concentrations of mono-ortho-substituted dioxin-like polychlorinated biphenyls (DL-PCBs) 105, 118 and 156 were measured by gas chromatography–mass spectrometry. Bone quality parameters (speed of sound, m/s; broadband ultrasound attenuation, dB/MHz; stiffness index, %) were assessed by quantitative ultrasound at the right calcaneus with the Achilles InSight system. Several factors known to be associated with osteoporosis were documented by questionnaire. Multiple linear regression models were constructed for the three ultrasound parameters.
DL-PCBs 105 and 118 concentrations, but not the global DLC concentration, were inversely associated with the stiffness index, even after adjusting for several confounding factors. The stiffness index (log) decreased by −0.22% (p=0.0414) and −0.04% (p=0.0483) with an increase of one μg/L in plasma concentrations of DL-PCB 105 and DL-PCB 118, respectively. Other factors, including age, height, smoking status, menopausal status and the percentage of omega-6 polyunsaturated fatty acids (PUFAs) in erythrocyte membranes were negatively associated with one of the ultrasound parameters, while the percentage of omega-3 PUFAs in these membranes and levels of physical activity and education were positively associated with them.
Our results show that an increase in plasma concentrations of DL-PCBs 105 and 118 was negatively associated with stiffness index, a measure of bone quality/strength, in women of this population. In addition to environmental contaminants, future studies should also consider PUFA intake as a factor influencing bone quality.
Osteoporosis is “a systemic skeletal disease characterized by low bone mass and micro-architectural deterioration of bone tissue, with a consequent increase in bone fragility and susceptibility to fractures” . It is a multifactorial polygenic disease in which genetic determinants are modulated by hormonal, environmental, and nutritional factors , is asymptomatic and progress silently with age . Eighty one percent of all fractures in women aged 50 years and older can be attributed to osteoporosis . Osteoporotic fractures, such as those of the hip, spine and wrist, often appear in older people following minor trauma. Hip fractures lead to rehabilitation problems and greatly decrease the quality of life [5, 6].
Very little information is available on bone mineral density (BMD), fracture risk and their determinants among Canada’s First Nations people. The only existing publications focus on Manitoba’s Aboriginal population [7–12]. Retrospective studies using administrative health data have reported that members of First Nations in Manitoba had significantly higher osteoporotic fracture rates than non-Aboriginals [8, 9]. Factors such as low socioeconomic status, diabetes  or the number of comorbidities and alcohol and drug use  were associated with higher fracture rates in these populations. Low calcium and vitamin D intake among Aboriginals, particularly older women , may also play a role.
Moreover, Aboriginal populations living in the northern regions of the globe are exposed to persistent organic pollutants through their traditional diet . These include dioxin-like compounds (DLCs), similar in structure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), which are lipophilic, resistant to biodegradation and therefore bioaccumulate and bioamplify in aquatic food chains, where Aboriginals obtain part of their diet . Yet the findings of in vitro studies  and experimental studies on rats [17–19] show that TCDD, the most toxic compound in the DLC family, has multiple effects on bone tissue, resulting in increased bone fragility. Few human studies have focused on the relationship between exposure to DLCs and bone quality/strength. Furthermore, while populations are exposed to a complex mixture of DLCs, only the plasma concentrations of certain dioxin-like PCBs (DL-PCBs) have been measured in participants in order to study the relationship between their concentrations and BMD or quantitative ultrasound (QUS) parameters [20–22].
The objective of this study was to investigate the relationship between plasma concentrations of DLCs or certain dioxin-like mono-ortho-substituted PCBs (DL-PCBs 105, 118 and 156) and calcaneal QUS parameters in Cree women of Eastern James Bay (Canada). Ultrasound bone measurement is a method used to assess bone strength that provides information that is complementary to BMD . Recent studies demonstrate that QUS at the calcaneus can predict fractures as effectively as DXA in postmenopausal women and men over the age of 65 . Values of QUS parameters are generally lower in osteoporotic patients than in healthy subjects . This technique has several advantages: it is simple, fast, non-invasive, radiation-free and inexpensive. Furthermore, ultrasound bone measurement devices are portable and easy to use in remote, isolated regions, such as those inhabited by northern Aboriginal communities, where measurement of BMD cannot be performed by dual-energy x-ray absorptiometry (DXA).
An environmental and health survey, entitled “Nituuchischaayihtitaau Aschii. Multi-community Environment- and Health Longitudinal Study in Eeyou Istchee”, was conducted from summer 2005 to summer 2009 in seven Cree communities located north of the 49th parallel north, in the Eastern James Bay region of the province of Quebec (Canada). The objective of this study was to investigate the health of these Aboriginal populations and the risk factors affecting them [26, 27].
Participants, stratified by age (0–7 years; 8–14 years; 15–39 years; ≥40 years) and gender, were selected using simple random sampling, without replacement, from the Beneficiary List in each community to create lists of potential participants to be contacted by recruiters. A weighting was attributed to each stratum to ensure that inferences could be drawn about the population. An initial list of subjects to be contacted was drawn up and all these people were contacted and asked to participate in the study. If they refused, a second list of participants was randomly created and recruitment continued in order to ensure the required number of subjects in each age group . Owing to the low response rate (53.5%) in the first community visited, a correction to the predefined samples had to be taken into account in order to maintain the desired statistical power; the sample size in subsequent communities was therefore increased . The total response rate of females (all age categories combined) varied between 34.4% and 70.9% depending on the community (Centre de Recherche du CHU de Québec, unpublished data). With respect to the non-participants, some could not be reached; some were undecided and withdrew, while others refused to participate. Others were excluded (unknown in the community, non-Cree, pregnant, disabled) [28, 29].
All participants (aged 8 years and over) completed a series of questionnaires in Eeyou Ayimuwin or English and administered by the research staff or interviewers (selected in local communities and provided with appropriate training) [30, 31]. They then underwent a clinical examination with biological samples and anthropometric measurements collected by research nurses.
All women aged 35 to 74 years who participated in the clinical examination were eligible for the calcaneal ultrasound measurement. A total of 254 women had QUS parameters measured. The volume of plasma collected from 249 women was sufficient to allow testing for DLCs and DL-PCBs.
Two QUS parameters were measured at the right calcaneus of Cree women using a portable Achilles Insight system (GE Healthcare Lunar, Madison, WI, USA): 1) speed of sound (SOS, m/s) and 2) broadband ultrasound attenuation (BUA, dB/MHz) . The Stiffness Index (SI, %) was automatically calculated by the system from the two parameters using the manufacturer’s formula [SI% = (0.67*BUA) + (0.28*SOS) – 420]. A research nurse inspected and calibrated the system’s membranes daily using the phantom provided by the manufacturer. In vitro accuracy was assessed by taking a number of repeated measurements with the manufacturer’s phantom: the mean coefficients of variation (CVs) recorded for the SOS parameter ranged from 0.05 to 0.12% depending on the community.
Body composition was measured using a bioelectric impedance analyzer (Tanita Corporation of America, IL, USA). Weight, lean body mass and body fat (kg), height, sitting height, waist and hip circumference (cm) were measured by research nurses using standardized techniques.
DLCs and DL-PCBs
Plasma samples were tested for DLCs and DL-PCBs at INSPQ’s Toxicology Laboratory (Quebec City, Canada). The luciferase reporter gene cell assay described by Medehouenou et al.  was used to measure AhR-mediated transcriptional activity elicited by plasma extracts and provide an integrated measure of DLCs. The limit of detection (LOD) was 30 pg TCDD-equivalents (EQ)/L. For quality control purposes, we used human plasma spiked with TCDD to obtain a concentration of 322 pg/L. The mean concentration of this sample tested 22 times was 361 pg TCDD-EQ/L (bias=12.1%) and the CV was 14%.
DL-PCBs (DL-PCB 105: 2,3,3’,4,4’-pentachlorobiphenyl; DL-PCB 118: 2,3’4,4’,5-pentachlorobiphenyl and DL-PCB 156: 2,3,3’,4,4’,5-hexachlorobiphenyl) were measured by gas chromatography–mass spectrometry (GC-MS) as described extensively in the project reports [34, 35]. For samples from the first Cree community investigated, a liquid-liquid extraction was used followed by separation and quantification of PCB congeners by GC-MS using an Agilent 6890 Network GC system (Wilmington, DE, USA), fitted with an Agilent 7683 series automatic liquid sampler and an Agilent 5973 Network mass spectrometer. The recovery rate for the three DL-PCB congeners was over 90%. Coefficients of variation of 5.9%, 4.2% and 11.5% were obtained for DL-PCBs 105, 118 and 156 congeners respectively following repeated measures of standard reference material (SRM) 1589 (National Institute of Standards and Technology, Gaithersburg, MD, USA). The limit of detection for the three DL-PCBs was 0.02 μg/L . For samples from the six other Cree communities investigated, a solid-phase extraction was used followed by GC-MS analysis with the same instrumentation as described above. Coefficients of variation of 7.4%, 7.7% and 9.7% were obtained for the three DL-PCB 105, 118 and 156 congeners respectively, following repeated measures of SRM 1589; the detection limit was 0.01 μg/L .
Concentrations of total lipids in plasma samples were calculated using Phillips et al.’s formula : [Total lipids (g/L) = (0.878*cholesterol mmol/L) + (0.885*triglycerides mmol/L) + 0.623]. Total cholesterol and triglyceride concentrations were determined by standard enzymatic methods. For the total cholesterol analysis, CVs of 1.5% and 1.3% were obtained for reference samples containing 2.87 and 6.67 mmol/L respectively. For the triglycerides analysis, CVs of 2.0% and 1.6% were noted for reference samples with concentrations of 1.07 and 2.32 mmol/L respectively .
Levels of omega-3 and omega-6 PUFAs were measured in erythrocyte membranes by liquid gas chromatography with flame ionisation detection at the Centre de Recherche sur les Maladies Lipidiques, CHU de Québec as previously described . Omega-3 PUFAs refer to the following fatty acids: α-linolenic acid (C18:3n-3) + docosapentaenoic acid (C22:5n-3) + docosahexaenoic acid (C22:6n-3) + eicosapentaenoic acid (C20:5n-3) + (C18:4n-3) + (C20:3n-3) + (C20:4n-3). Omega-6 PUFAs refer to the following fatty acids: linoleic acid (C18:2n-6) + arachidonic acid (C20:4n-6) + (C18:3n-6) + (C20:2n-6) + (C20:3n-6) + (C22:2n-6) + (C22:4n-6) + (C22:5n-6).
Apolipoprotein B (Apo B) was measured by nephelometry using a BN ProSpec station (Dade Behring, Mississauga, ON). Control levels of 0.46 and 1.44 g/L for Apo B showed interassay CVs of 3.2% and 1.8% .
Blood mercury levels (nmol/L) were measured by cold vapour atomic absorption spectrometry (Pharmacia) . Cadmium (nmol/L), selenium (μmol/L) and lead (μmol/L) were measured by inductively coupled plasma mass spectrometry (ICP-MS). The blood samples were diluted in ammonium hydroxide and the metals converted to their elemental form by aspirating the sample into argon plasma before being identified and quantified by mass spectrometry. The samples were analysed using a Perkin Elmer Sciex Elan 6000 ICP-MS instrument). The LODs were 0.04 nmol/L for cadmium, 0.001 μmol/L for lead, 0.49 nmol/L for mercury and 0.09 μmol/L for selenium .
Glucose was measured by a hexokinase enzymatic assay employing the Roche Modular system. Interassay CVs were 1.6% and 1.4% for glucose at control values of 4.74 and 15.66 mmol/L respectively . Alpha-tocopherol (vitamin E) was determined using a Waters high-pressure liquid chromatography system (Lachine, QC) equipped with an autosampler, a reverse phase column (Nucleosil ODS1) and a UV detector . Serum levels of 25-hydroxyvitamin D (25OHD) were measured in the Biochemistry Laboratory at Montreal’s Hôpital St-Luc (CHUM) by a procedure that includes extraction and quantification by competitive radioimmunoassay using the IDS radioimmunoassay kit (Medicorp Inc., Montreal, QC) .
Questionnaires were used to collect information about participants’ socio-demographic characteristics (date and place of birth; level of education: none or primary school/secondary and higher), their lifestyle habits (tobacco use: yes/no; use of calcium and vitamin D supplements in the last 12 months: yes/no; physical activity: active/inactive; alcohol intake in the last 12 months: yes/no; milk consumption in the last month: yes/no), their gynaecological history (menopausal status: menopausal/non-menopausal; parity: yes/no; number of children; use of hormonal contraceptives: yes/no; use of hormone replacement therapy, HRT: yes/no) and their personal fracture history (yes/no). Women were considered menopausal if they had had no menstrual periods for one year before recruitment. Sodium intake (≤ 2300 mg vs > 2300 mg/day) was compiled based on answers to the 24-hour dietary recall and data on the sodium content of foods .
We consulted the medical records of participants to document causes of secondary osteoporosis (CSO: yes/no), namely the use of certain medications in the last year (corticosteroids, heparin, Dilantin, Prednisone, chemotherapy) or the diagnosis of diseases (Cushing’s disease, rheumatoid arthritis, bone cancer, gastrectomy, kidney failure, liver failure, thyroid and parathyroid diseases) that may contribute to bone loss .
Physical activity (PA) was assessed using the short form of the International Physical Activity Questionnaire . A dichotomous variable was created based on the median value of the Total MET (metabolic equivalent)-minutes/week, calculated using the algorithm in the IPAQ document . Values at or below the median (3954 MET-minutes/week) define an inadequate level of PA, or subjects considered inactive, while values above the median define an adequate level of PA, or subjects considered active.
Our database comprised 249 participants. DLC concentrations were below the LOD (30 pg TCDD-EQ/L) for 16 participants. In these cases, we imputed a value between 0 and the LOD selected by simple random sampling with replacement. Of the 44 participants from the first community visited, four had DL-PCB 105 concentrations below the LOD, while one participant had a concentration of DL-PCB 156 below the LOD (LOD = 0.02 μg/L). Plasma concentrations of DL-PCBs 105, 118 and 156 below the LOD (0.01 μg/L) were obtained in 45, 7 and 12 participants respectively among the 205 residents of the other six communities. In these cases, we imputed a value equal to LOD/2, that is, 0.01 μg/L for participants from the first community and 0.005 μg/L for participants from the other communities.
Descriptive statistics (mean, standard error, minimum and maximum for quantitative variables; sample size and % per modality for qualitative variables) of the variables were generated for all participants.
The association between DLCs or each DL-PCB and QUS parameters was examined by simple and multiple linear regressions. Box-Cox procedures were used to resolve problems encountered with the hypotheses of normality and/or homoscedasticity in the multiple regression models for the dependent variables SOS and SI, which were subsequently log-transformed. Pearson’s correlation coefficients were calculated between the QUS parameters and the main exposure variables. The confounding effect of different factors on the relation between DLCs or DL-PCBs and QUS parameters was investigated. Secondary independent variables were tested to determine whether or not they should be included in the initial multiple regression models. Only those with a p value ≤ 0.20 (in simple linear regression analysis) were included. Multicollinearity between the variables of the initial regression models was tested to ensure that variables were not redundant. To avoid multicollinearity caused by the introduction of highly intercorrelated variables in the same linear regression model, we proposed for each dependent variable (SOS, BUA and SI) separate models that comprised only one of the main exposure variables (either DLCs, DL-PCB 105, DL-PCB 118 or DL-PCB 156). All potentially confounding variables were entered in the multiple regression model and those that did not change the value of β coefficient of the main exposure variable by more than 10% were removed one by one, starting with the variable that accounted for the least variation in the dependent variable. In the case of highly intercorrelated anthropometric variables, when one of these variables was not identified as a confounding factor, and therefore not retained in the final model, the initial model was rebuilt with the other variable.The final adjustment of the models was for the total plasma lipid concentration . We also constructed a second set of multivariate models (models II in Additional files 1, 2 and 3) adjusted for the same confounders and covariates. A p value of < 0.05 in a bilateral situation is considered statistically significant. The software used was SAS version 9.2 (SAS Institute Inc., Cary, NC, USA).
The project was approved by the Research Ethics Committees of Université Laval and McGill University, in partnership with McMaster University, and by the Research Committee of the Cree Board of Health and Social Services of James Bay. Participation in the study was voluntary and a consent form was signed by each participant. All information concerning the participants was kept strictly confidential.
With the weighting applied, the sample of 249 women (102 menopausal women and 147 non-menopausal women) could be extrapolated to a population of 2056 women (883 menopausal women and 1173 non-menopausal women).
QUS parameters, plasma DLC and DL-PCB concentrations as well as selected characteristics of the Cree women are shown in Table 1 and Additional file 4. Most participants were non-menopausal, had a secondary level education or higher, had had children, were sedentary, non-smokers and did not drink alcohol. A small proportion of women had used hormonal contraception or HRT, had taken supplements (calcium and vitamin D) in the 12 months prior to the study and had a history of fracture and conditions and/or medications identified as CSO.
Most of the Cree women (96.9%) had a waist circumference ≥ 88 cm (abdominal obesity according to Health Canada) ; the median value was 112 cm. With respect to body mass index (BMI, kg/m2), according to Health Canada criteria , only 3.5% of Cree women had a normal BMI (between 18.5 and 24.9), 19.2% were overweight (BMI between 25.0 and 29.9 kg/m2) and 77.3% were obese (BMI ≥ 30 kg/m2). Among Cree women, 9.0% had a height of less than 153 cm, which has been associated with osteoporosis .
With regard to blood mercury levels, 64.5% of Cree women exceeded the maximum blood mercury level of 16 nmol/L observed in women from the Quebec City area , whereas 13.2% exceeded Canada’s alert level of 99.7 nmol/L .
Only 1.9% of participants exhibited a blood selenium concentration exceeding the maximum level (3.6 μmol/L) recorded in women from the Quebec City area  and 1.8% showed a blood cadmium concentration over the maximum value in the Quebec City inhabitants (55 nmol/L) .
With respect to serum levels of vitamin D (25OHD), 2.3% of Cree women had a level considered as critical (< 27.5 nmol/L) ; 14.6% displayed inadequate levels (< 37.5 nmol/L), while 39.5% showed a level considered minimal (< 50 nmol/L) .
Fasting glucose levels were considered normal (< 6.1 mmol/L) in 56.0% of Cree women, while 44.0% had levels corresponding to hyperglycemia (≥ 6.1 mmol/L). For 26.7% of Cree women, glucose levels could correspond to type 2 diabetes (≥ 7.0 mmol/L) according to the Canadian Diabetes Association .
Pearson’s r correlation coefficients revealed that QUS parameters were negatively and significantly correlated with DLCs, DL-PCBs and age, while the three DL-PCBs were positively and significantly correlated with one another and with age (see Additional file 5).
The final multivariate models of each QUS parameter, constructed for each of the main exposure variables (DLCs or DL-PCBs), are presented in Tables 2, 3 and 4. DL-PCBs 105 and 118 were significantly and negatively associated with SI (log); these associations persisted even after adjustment for several confounding factors. Hence the SI (log) decreased by −0.22% (p=0.0414) and −0.04% (p=0.0483) for every one-unit increase (1 μg/L) in plasma concentrations of DL-PCBs 105 and 118, respectively (Table 4). In QUS parameters models built for each of the main exposure variables and adjusted for the same factors (see Additional files 1, 2 and 3, models II), DL-PCB 105 was negatively and significantly associated with BUA (Additional file 2), whereas DL-PCBs 105 and 118 were negatively and significantly associated with SI (log) (Additional file 3).
Several variables were identified as confounders of associations between exposure variables and QUS parameters. Some confounding factors were significantly associated with SOS (log): age (in all models), level of education (DLC, DL-PCB 105 and DL-PCB 118 models), height (all DL-PCB models), PA (DL-PCB 105 and DL-PCB 118 models), smoking status and menopausal status (DL-PCB 105 model) (Table 2). In general, SOS (log) decreased with increasing age, height, menopausal status and smoking status and SOS (log) increased with increasing level of education and physical activity.
Age (in all models), level of education (DLC, DL-PCB 118 and DL-PCB 156 models) and % of omega-6 PUFAs (DLC model) were significantly associated with BUA (Table 3). In general, BUA decreased with increasing age and % of omega-6 PUFAs, and increased with increasing level of education.
Lastly, age (in all models), level of education (DLC, DL-PCB 118 and DL-PCB 156 models) and % of omega-3 PUFAs (DL-PCB 118 model) were significantly associated with SI (Table 4). In general, SI (log) decreased with increasing age, and increased with an increase in the total % of omega-3 PUFAs and with increasing level of education.
This is the first population study to investigate the relationship between total DLC concentration measured with an AhR responsive reporter gene assay and bone quality parameters . In Cree women aged 35 to 74 years of Eastern James Bay, DL-PCBs 105 and 118 were negatively and significantly associated with SI. However, we did not find any association between total DLC concentrations measured by the reporter gene bioassay and QUS parameters, which suggests that the effect of PCBs is not linked to activation of the AhR signalling pathway. Alternatively, the lack of association could also be due to the lower precision of the bioassay measurement compared to that of GC-MS measurements.
Three studies have been conducted on the relation between DL-PCBs and QUS bone parameters, with quite controversial findings. In a group of Swedish men from the general population (n = 115, age 40–75 years), Glynn et al.  reported a positive association at the limit of significance (β = 0.0044; p = 0.05) between serum concentrations of DL-PCB 167 and whole-body BMD (measured by DXA), after adjustment for several confounding factors. However, no association was observed between DL-PCB 167 and BUA or SOS parameters measured by ultrasound at the left calcaneus . No association was observed between DL-PCBs 105, 118, 156 and QUS parameters or lumbar spine, femoral neck or whole-body BMD .
Exploratory multivariate analyses conducted by Côté et al.  revealed that plasma concentrations of DL-PCB 156 were negatively associated with QUS parameters measured in 153 peri- and postmenopausal Inuit women from Nuuk in Greenland (SOS model: β = −22.68, p = 0.014; BUA model: β = −8.12, p = 0.028; SI model: β = −11.95, p = 0.009) . However, DL-PCBs 105 and 118 were not associated with QUS parameters .
Hodgson et al.  found a negative association between serum DL-PCB 118 concentration and forearm BMD measured using DXA (β = −0.00024, p = 0.002) in the dominant arm of 154 men (60–81 year old) from the Swedish OSCAR cohort study. The odds ratio for low BMD of 1.06 per 10 pg/mL DL-PCB 118 was significant (95% CI: 1.01-1.12) . However, in women (n = 167), DL-PCB 118 was positively associated with BMD (β = 0.00008, p = 0.045). In this study, no association was observed between DL-PCBs 105, 156, 157, 167 and BMD; the authors considered the statistical power to be limited and a large proportion of variance was not explained by regression models .
In additional analyses, we performed logistic regression using medians for SI (91%), DL-PCB 105 (0.03 μg/L) and DL-PCB 118 (0.18 μg/L) as cut-off values between “low” and “high” categories. Statistically-significant odds ratios (OR) of 1.69 for low SI at high DL-PCB 105 exposure (95% CI: 1.34-2.13) and 1.36 for low SI at high DL-PCB 118 exposure (95% CI: 1.06-1.73) were obtained in models adjusted for all variables indicated in Table 4, confirming the associations obtained by multiple linear regression. However, the choice of medians as a cut-off values to define “low” and “high” categories of SI and the main exposure variables is arbitrary and does not convey the same meaning in different populations. Cut-off values for QUS parameters based on pathophysiological considerations are not available in the literature.
The geometric means (GMs) of plasma DL-PCB concentrations among Cree women were higher than those reported in the 2007–2009 Canadian Health Measures Survey for women aged 40 to 79 years and representative of the general Canadian population . However, mean plasma concentrations of DLCs in Cree women estimated by the reporter gene assay were lower than those of women of childbearing age (Inuit, Dene, Caucasian) in other regions in Northern Canada .
Explanatory multivariate analyses revealed that several other factors were significantly associated with QUS parameters. Negative associations were found with age, height, smoking status, menopausal status and omega-6 PUFAs, whereas positive associations were noted with level of education, PA and % of omega-3 PUFAs.
Omega-3 PUFAs measured as a % of total fatty acids in erythrocyte membrane phospholipids in Cree women were positively and significantly associated with SI, while the % of omega-6 PUFAs was negatively associated with the BUA parameter. As observed following additional analyses, the omega-3/omega-6 PUFA ratio was positively and significantly associated with BUA and SI (log). To our knowledge, associations between omega-3 and omega-6 PUFAs and calcaneal QUS parameters have not been reported previously. A high ratio of omega-6 to omega-3 PUFAs has been associated with low BMD , while a higher dietary intake of omega-3 PUFAs had a protective effect against bone loss [51, 52].
Our study has several strengths. Firstly, the representativeness of our population sample was assured by the recruitment strategy and the weighting scheme of the study, which took non-response and refusal to participate rates into account. Hence the findings can be generalized to the entire female Cree population aged 35 to 74 years living in Eastern James Bay communities. Secondly, we took into consideration a large number of potential confounding factors of the association between the main exposure variables (DLCs, DL-PCBs) and QUS parameters. Several of these factors were investigated for the first time (% omega-3 and omega-6 PUFAs, mercury, selenium, lead, glucose, apolipoprotein B). Thirdly, measurement biases in the dependent variables or exposure variables are relatively unlikely. Data collection in the field (ultrasound measurements, anthropometric measurements, biological samples) was performed by research nurses using standardized techniques. Completed questionnaires were reviewed by members of the research team to ensure that all the questions had been completely and properly answered [30, 31]. Individual levels of exposure to contaminants were measured using standardized methods in the same laboratory (INSPQ) which has a quality control system accredited according to ISO/CEI 17025, CAN-P-43 and ILAC G-13 standards. All laboratory analyses (plasma or blood) were carried out, with two exceptions (PUFAs and 25OHD), at the INSPQ laboratory. Measurement accuracy complied with applicable standards.
Our study has also certain limitations. Firstly, the main methodological limitation is its cross-sectional design, with exposure and the dependent variable measured at the same time (a single measurement for each subject), such that the temporal sequence of cause and effect cannot be determined. However, due to the fact that DLCs and DL-PCBs accumulate in the body with age, plasma concentrations in Cree women reflect their lifelong exposure. Secondly, the participation rate in the study varied largely between communities, which may suggest a selection bias. However, this bias is quantitatively unimportant, since it is unlikely that the characteristics of the subjects included in the study are different from those of all eligible persons. The lists of beneficiaries which were used to build the list of contact persons were completed as necessary with information from other local listings; moreover, the investigated communities are small and live in reserves, which facilitated the identification and location of their members. Regarding the non-response and refusal to participate in the study, it is unlikely that these subjects have different levels of plasma DL-PCBs congeners or QUS parameter values, compared with participants; the recruitment of participants was made for a survey of general health and not specifically for the purpose of our study. Thirdly, the number of participants with QUS measurements was relatively small, limiting our capacity to observe an association between DLC concentrations determined by the reported gene assay and QUS parameters, especially considering that this measure is relatively imprecise compared to analytical chemistry based data. Fourthly, while we took a maximum number of potential confounding factors into account given the size of our sample, a residual confounding effect cannot be ruled out. The determinants of bone strength are multigenic and multifactorial and other factors (genetic, nutritional, environmental) not measured in our study could have decreased the residual variance in multiple regression models. Lastly, two variables that are known to influence bone quality could not be thoroughly assessed. The assessment of dietary calcium intake considered only milk consumption. Physical activity was evaluated using the short form of the IPAQ , an instrument that proved difficult to administer in Cree communities, even though its feasibility and validity has been demonstrated in general populations (adults aged 15 to 69 years) of 12 countries . Participants in our study had trouble quantifying their activities, especially the duration and intensity of activities during the week prior to the study . We believe that the proportion of active Cree women was probably overestimated.
We observed negative associations between SI values and plasma concentrations of two dioxin-like PCB congeners in women from this Aboriginal population. Our study does not allow causal inference and residual confounding may still be present. Therefore, we do not believe that dietary recommendations to reduce PCB exposure are warranted at this time. In multivariate models, associations were also noted with dietary and lifestyle factors, suggesting avenues to improve bone quality in this population, such as smoking cessation and increases in omega-3 PUFA intake and physical activity.
aryl hydrocarbon receptor
- Apo B:
bone mineral density
body mass index
dioxin-like polychlorinated biphenyl
broadband ultrasound attenuation
Centre Hospitalier de l’Université de Montréal
- CHU de Québec:
Centre Hospitalier Universitaire de Québec
- 95% CI:
95% confidence interval
causes of secondary osteoporosis
coefficient of variation
decibels per microhertz
dioxin-like polychlorinated biphenyls
dual-energy x-ray absorptiometry
gas chromatography–mass spectrometry
grams per litre
hormone replacement therapy
inductively coupled plasma mass spectrometry
Institut National de Santé Publique de Québec
metres per second
micrograms per litre
micromoles per litre
millimoles per litre
nanomoles per litre
limit of detection
polyunsaturated fatty acid
- pg TCDD-EQ/L:
picograms TCDD-equivalents per litre
speed of sound
standard reference material
WHO: Assessment of fracture risk and its application to screening for postmenopausal osteoporosis. 1994, Geneva: Report of a WHO Study Group, 130-
Gennari L, Merlotti D, De Paola V, Calabrò A, Becherini L, Martini G, Nuti R: Estrogen receptor gene polymorphisms and the genetics of osteoporosis: a HuGE review. Am J Epidemiol. 2005, 161: 307-320. 10.1093/aje/kwi055.
Saag KG, Geusens P: Progress in osteoporosis and fracture prevention: focus on postmenopausal women. Arthritis Res Ther. 2009, 11: 251-10.1186/ar2815.
Bessette L, Ste-Marie LG, Jean S, Davison KS, Beaulieu M, Baranci M, Bessant J, Brown JP: The care gap in diagnosis and treatment of women with a fragility fracture. Osteoporos Int. 2008, 1: 79-86.
Adachi JD, Ioannidis G, Pickard L, Berger C, Prior JC, Joseph L, Hanley DA, Olszynski WP, Murray TM, Anastassiades T, Hopman W, Brown JP, Kirkland S, Joyce C, Papaioannou A, Poliquin S, Tenenhouse A, Papadimitropoulos EA: The association between osteoporotic fractures and health-related quality of life as measured by the Health Utilities Index in the Canadian Multicentre Osteoporosis Study (CaMos). Osteoporos Int. 2003, 14: 895-904. 10.1007/s00198-003-1483-3.
Cummings SR, Melton LJ: Epidemiology and outcomes of osteoporotic fractures. Lancet. 2002, 359: 1761-1767. 10.1016/S0140-6736(02)08657-9.
Leslie WD, Derksen S, Metge C, Lix LM, Salamon EA, Wood SP, Roos LL: Fracture risk among First Nations people: a retrospective matched cohort study. CMAJ. 2004, 171: 869-873.
Leslie WD, Derksen SA, Metge C, Lix LM, Salamon EA, Steiman PW, Roos LL: Demographic risk factors for fracture in First Nations people. Can J Public Health. 2005, 96: 45-50.
Leslie WD, Derksen S, Prior HJ, Lix LM, Metge C, O‘neil J: The interaction of ethnicity and chronic disease as risk factors for osteoporotic fractures: a comparison in Canadian Aboriginals and non-Aboriginals. Osteoporos Int. 2006, 17: 1358-1368. 10.1007/s00198-006-0111-4.
Leslie WD, Metge CJ, Weiler HA, Doupe M, Steiman PW, O’Neil JD: Bone density and bone area in Canadian Aboriginal women: the First Nations Bone Health Study. Osteoporos Int. 2006, 17: 1755-1762. 10.1007/s00198-006-0184-0.
Weiler HA, Leslie WD, Bernstein CN: Parathyroid hormone is predictive of low bone mass in Canadian Aboriginal and White women. Bone. 2008, 42: 498-504. 10.1016/j.bone.2007.11.013.
Leslie WD, Weiler HA, Lix LM, Nyomba BL: Body composition and bone density in Canadian White and Aboriginal women: the First Nations Bone Health Study. Bone. 2008, 42: 990-995. 10.1016/j.bone.2008.01.018.
Weiler HA, Leslie WD, Krahn J, Steiman PW, Metge CJ: Canadian Aboriginal women have a higher prevalence of vitamin D deficiency than non-Aboriginal women despite similar dietary vitamin D intakes. J Nutr. 2007, 137: 461-465.
Van Oostdam J, Donaldson SG, Feeley M, Arnold D, Ayotte P, Bondy G, Chan L, Dewailly E, Furgal CM, Kuhnlein H, Loring E, Muckle G, Myles E, Receveur O, Tracy B, Gill U, Kalhok S: Human health implications of environmental contaminants in Arctic Canada: A review. Sci Total Environ. 2005, 351–352: 165-246.
Lorber M, Patterson D, Huwe J, Kahn H: Evaluation of background exposures of Americans to dioxin-like compounds in the 1990s and the 2000s. Chemosphere. 2009, 77: 640-651. 10.1016/j.chemosphere.2009.08.016.
Carpi D, Korkalainen M, Airoldi L, Fanelli R, Hakansson H, Muhonen V, Tuukkanen J, Viluksela M, Pastorelli R: Dioxin-sensitive proteins in differentiating osteoblasts: effects on bone formation in vitro. Toxicol Sci. 2009, 108: 330-343. 10.1093/toxsci/kfp021.
Lind PM, Wejheden C, Lundberg R, Alvarez-Lloret P, Hermsen SA, Rodriguez-Navarro AB, Larsson S, Rannug A: Short-term exposure to dioxin impairs bone tissue in male rats. Chemosphere. 2009, 75: 680-684. 10.1016/j.chemosphere.2008.12.024.
Herlin M, Kalantari F, Stern N, Sand S, Larsson S, Viluksela M, Tuomisto JT, Tuomisto J, Tuukkanen J, Jämsä T, Lind PM, Håkansson H: Quantitative characterization of changes in bone geometry, mineral density and biomechanical properties in two rat strains with different Ah-receptor structures after long-term exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin. Toxicology. 2010, 273: 1-11. 10.1016/j.tox.2010.04.006.
Finnilä MA, Zioupos P, Herlin M, Miettinen HM, Simanainen U, Håkansson H, Tuukkanen J, Viluksela M, Jämsä T: Effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin exposure on bone material properties. J Biomech. 2010, 43: 1097-1103. 10.1016/j.jbiomech.2009.12.011.
Glynn AW, Michaëlsson K, Lind PM, Wolk A, Aune M, Atuma S, Darnerud PO, Mallmin H: Organochlorines and bone mineral density in Swedish men from the general population. Osteoporos Int. 2000, 11: 1036-1042.
Côté S, Ayotte P, Dodin S, Blanchet C, Mulvad G, Petersen HS, Gingras S, Dewailly E: Plasma organochlorine concentrations and bone ultrasound measurements: a cross-sectional study in peri-and postmenopausal Inuit women from Greenland. Environ Health. 2006, 5: 33-10.1186/1476-069X-5-33.
Hodgson S, Thomas L, Fattore E, Lind PM, Alfven T, Hellström L, Håkansson H, Carubelli G, Fanelli R, Jarup L: Bone mineral density changes in relation to environmental PCB exposure. Environ Health Perspect. 2008, 116: 1162-1166. 10.1289/ehp.11107.
Padilla F, Laugier P: Recent developments in trabecular bone characterization using ultrasound. Curr Osteoporos Rep. 2005, 3: 64-69. 10.1007/s11914-005-0006-x.
Krieg MA, Barkmann R, Gonnelli S, Stewart A, Bauer DC, Del Rio BL, Kaufman JJ, Lorenc R, Miller PD, Olszynski WP, Poiana C, Schott AM, Lewiecki EM, Hans D: Quantitative ultrasound in the management of osteoporosis: the 2007 ISCD Official Positions. J Clin Densitom. 2007, 2008 (11): 163-187.
Danese RD, Licata AA: Ultrasound of the skeleton: review of its clinical applications and pitfalls. Curr Rheumatol Rep. 2001, 3: 245-248. 10.1007/s11926-001-0025-2.
Bonnier-Viger Y, Dewailly E, Egeland GM, Nieboer E, Pereg D: Nituuchischaayihtitaau Aschii. Multi-community Environment- and Health longitudinal study in IIyiyiu Aschii: Mistissini. Technical report: summary of activities, results and recommendations. Edited by: Pereg D, Nieboer E. 2007, Montreal, QC: Cree Board of health and Social Services of James Bay, 389-http://www.creehealth.org/sites/default/files/Environmental%20Health%20Study-%20Technical%20Report%20of%20Mistissini.pdf,
Bonnier-Viger Y, Château-Degat ML, Dewailly E, Egeland GM, Nieboer E: Nituuchischaayihtitaau Aschii multi-community environment and health longitudinal study in Eeyou Istchee: Eastmain and Wemindji Technical Report: summary of 2007 activities, results and recommendations. Public Health Report Series 4 on the Health of the Population. Edited by: Nieboer E, Robinson E, Petrov K. 2011, Chisasibi, QC: Cree Board of health and Social Services of James Bay, 376-http://www.creehealth.org/sites/default/files/Eastmain-Wemindji%20report%20pdf%2015-06-11_0.pdf,
Château-Degat ML, Pereg D, Louchini R: Study Population, Recruitment, Ethics and Confidentiality. Nituuchischaayihtitaau Aschii Multi-community Environment- and Health longitudinal study in IIyiyiu Aschii: Mistissini Technical report: summary of activities, results and recommendations. Edited by: Pereg D, Nieboer E. 2007, Montreal, QC: Cree Board of health and Social Services of James Bay, 45-46.
Côté S, Annasour-Laouan-Sidi E, Château-Degat ML, Nieboer E: Study Population, Recruitment, Ethics and Confidentiality. Nituuchischaayihtitaau Aschii multi-community environment and health longitudinal study in Eeyou Istchee: Eastmain and Wemindji Technical Report: summary of 2007 activities, results and recommendations. Edited by: Nieboer E, Robinson E, Petrov K. 2011, Chisasibi, QC: Public Health Report Series 4 on the Health of the Population. Cree Board of health and Social Services of James Bay, 26-
Côté S: Questionnaires. Nituuchischaayihtitaau Aschii Multi-community Environment- and Health longitudinal study in IIyiyiu Aschii: Mistissini Technical report: summary of activities, results and recommendations. Edited by: Pereg D, Nieboer E. 2007, Montreal, QC: Cree Board of health and Social Services of James Bay, 47-48.
Côté S: Questionnaires. Nituuchischaayihtitaau Aschii multi-community environment and health longitudinal study in Eeyou Istchee: Eastmain and Wemindji Technical Report: summary of 2007 activities, results and recommendations. Edited by: Nieboer E, Robinson E, Petrov K. 2011, Chisasibi, QC: Public Health Report Series 4 on the Health of the Population. Cree Board of health and Social Services of James Bay, 29-30.
Madison, Wisconsin. Achilles InSight - Manuel de l’opérateur.http://www3.gehealthcare.com/en/Products/Categories/Bone_Health/Quantitative_Ultrasound/~/media/Downloads/us/Product/Product-Categories/Bone-Health/Quantitative%20Ultrasound/Achilles/GEHealthcare-Brochure_Achilles-Bone-Ultrasonometers.pdf,
Medehouenou TC, Larochelle C, Dumas P, Dewailly E, Ayotte P: Determinants of AhR-mediated transcriptional activity induced by plasma extracts from Nunavik Inuit adults. Chemosphere. 2010, 80: 75-82. 10.1016/j.chemosphere.2010.04.017.
Pereg D, Ayotte P, Château-Degat ML, Lamarche B, Julien P: Laboratory Analyses. Nituuchischaayihtitaau Aschii Multi-community Environment- and Health longitudinal study in IIyiyiu Aschii: Mistissini Technical report: summary of activities, results and recommendations. Edited by: Pereg D, Nieboer E. 2007, Montreal, QC: Cree Board of health and Social Services of James Bay, 52-60.
Ayotte P, Château-Degat ML, Côté S, Bertrand F, Julien P, Lamarche B: Laboratory Analyses. Nituuchischaayihtitaau Aschii multi-community environment and health longitudinal study in Eeyou Istchee: Eastmain and Wemindji Technical Report: summary of 2007 activities, results and recommendations. Edited by: Nieboer E, Robinson E, Petrov K. 2011, Chisasibi, QC: Public Health Report Series 4 on the Health of the Population. Cree Board of health and Social Services of James Bay, 34-40.
Phillips DL, Pirkle JL, Burse VW, Bernert JT, Henderson LO, Needham LL: Chlorinated hydrocarbon levels in human serum: effects of fasting and feeding. Arch Environ Contam Toxicol. 1989, 18: 495-500. 10.1007/BF01055015.
Votre santé et vous. Sodium.http://www.hc-sc.gc.ca/hl-vs/alt_formats/pdf/iyh-vsv/food-aliment/sodium-fra.pdf,
The Journal of The North American Menopause Society: Management of osteoporosis in postmenopausal women: 2010 Position Statement of the North American Menopause Society. Menopause. 2010, 17: 25-54. 10.1097/gme.0b013e3181c617e6.
Guidelines for data processing and Analysis of the International Physical Activity Questionnaire (IPAQ) – Short and Long Forms. 2005,http://www.ipaq.ki.se/scoring.pdf, November ,
Schisterman EF, Whitcomb BW, Louis GM, Louis TA: Lipid adjustment in the analysis of environmental contaminants and human health risks. Environ Health Perspect. 2005, 113: 853-857. 10.1289/ehp.7640.
Lignes directrices canadiennes pour la classification du poids chez les adultes.http://www.hc-sc.gc.ca/fn-an/alt_formats/hpfb-dgpsa/pdf/nutrition/cg_quick_ref-ldc_rapide_ref-fra.pdf,
Newton JL, Jones DEJ, Wilton K, Pairman J, Parry SW, Francis RM: Calcaneal bone mineral density in older patients who have fallen. Q J Med. 2006, 99: 231-236. 10.1093/qjmed/hcl026.
LeBlanc A, Lapointe S, Beaudet A, Côté I, Dumas P, Labrecque F, Lamy C, Larochelle J, Lepage L, Pelletier F, Weber JF, Levallois P, Gingras S: Étude sur l'établissement de valeurs de référence d'éléments traces et métaux dans le sang, le sérum et l‘urine de la population de la grande région de Québec. 2004, Institut National de Santé Publique du Québec: Direction toxicologie humaine, Direction risques biologiques, environnementaux et occupationnels,http://www.inspq.qc.ca/pdf/publications/289-ValeursReferenceMetaux.pdf,
Santé Québec: A Health Profile of the Inuit. Report of the Santé Québec Health Survey Among the Inuit of Nunavik, 1992. Edited by: Jetté M. 1994, Montréal: Ministère de la Santé et des Services sociaux, Gouvernement du Québec
Mercury. Your health and the environment. A resource tool.http://www.acpo.org.br/MPF_PAPA_LAMPADAS/06_ANEXO_06_health_canada.pdf,
Intaraprasert S, Taneepanichskul S, Theppisai U, Chaturachinda K: Bone density in women receiving Nortplant implants for contraception. J Med Assoc Thai. 1997, 80: 738-741.
Association canadienne du diabète: Lignes directrices de pratique clinique 2008 de l’Association canadienne du diabète pour la prévention et le traitement du diabète au Canada. Can J Diab. 2008, 32: 239-http://www.diabetes.ca/documents/about-diabetes/CPG_FR.pdf,
Paunescu AC: PhD thesis. Les composés dioxin-like mesurés par DR-CALUX et les paramètres osseux évalués par ultrasonographie chez les femmes cries et inuites du Nord-du-Québec et du Groenland. 2012, Université Laval: Département de Médecine Sociale et Préventive
Santé Canada: Rapport sur la biosurveillance humaine des substances chimiques de l’environnement au Canada. Résultats de l’Enquête canadienne sur les mesures de la santé Cycle 1 (2007 à 2009). 2010, 300-http://www.hc-sc.gc.ca/ewh-semt/alt_formats/hecs-sesc/pdf/pubs/contaminants/chms-ecms/report-rapport-fra.pdf,
Weiss LA, Barrett-Connor E, von Mühlen D: Ratio of n-6 to n-3 fatty acids and bone mineral density in older adults: the Rancho Bernardo Study. Am J Clin Nutr. 2005, 81: 934-938.
Kruger MC, Coetzer H, de Winter R, Gerike G, van Papendorp DH: Calcium, gamma-linolenic acid and eicosopentaenoic acid supplementation in senile osteoporosis. Aging Clin Exp Res. 1998, 10: 385-394.
Fernandes G, Bhattacharya A, Rahman M, Zaman K, Banu J: Effects of n-3 fatty acids on autoimmunity and osteoporosis. Front Biosci. 2008, 13: 4015-4020.
Craig CL, Marshall AL, Sjostrom M, Bauman AE, Booth ML, Ainsworth BE, Pratt M, Ekelund U, Yngve A, Sallis JF, Oja P: International Physical Activity Questionnaire: 12-country reliability and validity. Medicine and Science in Sports and Exercise. Med Sci Sports Exerc. 2003, 35: 1381-1395. 10.1249/01.MSS.0000078924.61453.FB.
Egeland G, Denomme D, Johnson-Down L: Physical Activity. Nituuchischaayihtitaau Aschii Multi-community Environment- and Health longitudinal study in IIyiyiu Aschii: Mistissini Technical report: summary of activities, results and recommendations. Edited by: Pereg D, Nieboer E. 2007, Montreal, QC: Cree Board of health and Social Services of James Bay, 64-66.
We would like to thank the Cree women of Eastern James Bay who participated in this health study. We would also like to thank the medical and research staff and the Cree interpreters for their assistance and for data collection. This scientific communication is a report from the Nituuchischaayihtitaau Aschii Multi-Community Environment and Health Longitudinal Study in Iiyiyiu Aschii supported by the Cree people of northern Québec, the Cree First Nations and the Cree Board of Health and Social Services of James Bay through financial contributions from Niskamoon Corporation. ACP is the recipient of a Doctoral Training Award from FRQ-S and a doctoral scholarship from Nasivvik - Centre for Inuit Health and Changing Environments.
The authors declare that they have no competing interests.
ED, PA, SD and EN conceived the study design. ACP participated to data collection, performed the statistical analysis, interpreted the data, and drafted the manuscript. PA interpreted the data and helped draft the manuscript. All authors revised the manuscript critically. All authors approved the final version.
Electronic supplementary material
About this article
Cite this article
Paunescu, A., Dewailly, É., Dodin, S. et al. Dioxin-like compounds and bone quality in Cree women of Eastern James Bay (Canada): a cross-sectional study. Environ Health 12, 54 (2013). https://doi.org/10.1186/1476-069X-12-54
- Calcaneal Ultrasound Parameters
- Dioxin-Like Compounds
- Polychlorinated Biphenyls
- Cree Women