Article Text

Utility of current obesity thresholds in signaling diabetes risk in the James Bay Cree of Eeyou Istchee
  1. Priya Manjoo1,
  2. David Dannenbaum2,3,
  3. Lawrence Joseph4,
  4. Jill Torrie3,
  5. Kaberi Dasgupta4
  1. 1Department of Medicine, University of Victoria, Victoria, British Columbia, Canada
  2. 2Department of Family Medicine, McGill University, Montreal, Quebec, Canada
  3. 3Public Health Department of the Cree Territory, Cree Board of Health & Social Services of James Bay, Chisasibi, Quebec, Canada
  4. 4Department of Medicine, McGill University, Montreal, Quebec, Canada
  1. Correspondence to Dr Kaberi Dasgupta; kaberi.dasgupta{at}


Objective The anthropometric thresholds signaling type 2 diabetes risk have not been well defined for Aboriginal communities. This study examined current thresholds in terms of ability to capture diabetes risk in the Cree of Eeyou Istchee in northern Quebec, Canada.

Research design and methods The study cohort for this analysis included adult participants from the Nituuchischaayihtitaau Aschii Multi-Community Environment and Health Study with complete data on anthropometric measures, fasting glucose, and insulin. Diabetes risk was defined as Homeostatic Model Assessment of Insulin Resistance (HOMA-IR) value >2. Positive and negative likelihood ratios (PLR, NLR) of existing obesity thresholds were evaluated (WHO; International Diabetes Federation, IDF; Adult Treatment Panel III, ATP III). Receiver operating curves were examined to estimate optimal thresholds. In a sensitivity analysis, diabetes risk was defined as HOMA-IR >2.7.

Results The WHO 30 kg/m2 body mass index (BMI) threshold performed well in women (PLR 5.56, 95% CI 1.95 to 15.9; NLR 0.24, 95% CI 0.19 to 0.31) and men (PLR 7.51, 95% CI 2.94 to 19.2; NLR 0.33, 95% CI 0.27 to 0.41). It was close to the estimated optimal threshold (28.5 kg/m2). The ATP III waist circumference threshold (102 cm) performed well in men (PLR 4.64, 95% CI 2.47 to 8.71; NLR 0.21, 95% CI 0.16 to 0.28) and was close to the estimated optimal threshold (101 cm). With diabetes risk defined at HOMA-IR >2.7, PLR values were slightly lower with narrower 95% CIs and optimal thresholds were slightly higher; PLR values remained above 3. For other current thresholds, estimated optimal values were higher and none had a PLR above 2.

Conclusions A BMI of 30 kg/m2 in women and men, and a 102 cm waist circumference in men, are meaningful obesity thresholds in this Aboriginal population. Other thresholds require a further evaluation.

  • Aboriginal Health
  • Obesity and Type 2
  • Type 2 Diabetes

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See:

Statistics from

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Key messages

  • Being at or above the WHO obesity body mass index threshold of 30 kg/m2 or more signals elevated insulin resistance in Cree women and men.

  • Being at or above the Adult Treatment Panel III abdominal obesity waist circumference of 102 cm or more signals insulin resistance in Cree men.

  • The relevance of other obesity thresholds to insulin resistance in the Cree is not clear.


Worldwide, Aboriginal populations consistently demonstrate high type 2 diabetes risk; a recent review, for example, highlighted a twofold to threefold risk increase in Australia, New Zealand, and the USA.1 Obesity is a key driver.2 The body mass index (BMI, weight in kg/height in m2) is the metric most frequently used to capture obesity, defined by the WHO at 30 kg/m2.3 There may, however, be some ethnocultural differences in what constitutes a meaningful BMI threshold. At a given BMI, some studies indicate that Aboriginal groups, Asians, and South Asians, have a greater burden of cardiometabolic risk factors.4–7 Although not a WHO recommendation, there is some evidence that in Asians, the obesity threshold should be below 27 kg/m2.8 The BMI threshold signaling diabetes risk has not been specifically defined for individual Aboriginal communities or among Aboriginal communities in general.

Similarly, Aboriginal group-specific diabetes risk thresholds have not been defined for other anthropometric measures of obesity, including waist circumference (WC) and waist to hip ratio (WHR). These measures are indicators of abdominal, central, or visceral adiposity. Visceral adipose tissue liberates proinflammatory cytokines such as interleukin 6, tumor necrosis factor α, free fatty acids and resistin,9 which promote a state of insulin resistance, a precondition for type 2 diabetes. WC has a closer correlation with visceral adiposity than BMI.10–12 In the large European Interact study, WC and BMI demonstrated independent associations with diabetes incidence.13 WHR has also demonstrated associations with diabetes.14 Ethnocultural differences have been described for relationships of WC and WHR with cardiometabolic risk factors. For example, South Asian individuals have higher glucose levels than European individuals at an equivalent WC.15 In a study comparing anthropometric risk factors in terms of ability to predict cardiometabolic risk factors in Aboriginal Australians, WHR demonstrated the greatest predictive power.16 The International Diabetes Federation (IDF) does specify different visceral adiposity-defining anthropometric thresholds for some ethnocultural groups,17 but none are specific for Aboriginal populations.

The IDF, WHO, and Third Report of the Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults (Adult Treatment Panel III, ATP III), have defined anthropometric thresholds for obesity. The aim of the present analysis was to evaluate these thresholds in a particular Aboriginal group, in terms of ability to detect individuals at risk for type 2 diabetes.

The group studied was the Cree of the eastern James Bay region (Eeyou Istchee) in Canada. In this population, diabetes prevalence rose 10-fold from 2.4% in 1983 to 22.1% in 2011 (non-standardized rates). Between the 1940s and 1970s, geographic displacement was imposed to allow natural resource exploitation (eg, forestry, mining, and hydroelectric development). In these Cree communities, there was a shift in food procurement methods from hunting, fishing, and berry picking, until at least the 1970s, to a high intake of processed foods and beverages. A nomadic way of life was supplanted by permanent year-round settlement within nine communities. Snowshoeing and canoeing activities were replaced by snowmobiles and motorboats. As recently documented by the Truth and Reconciliation commission,18 until 1969, Aboriginal Canadians children were forced to leave their families and attend boarding schools where many were physically, sexually, and psychologically abused; mortality rates were high in what has been described as cultural genocide. All of these factors have likely contributed to the diabetes epidemic faced by Aboriginal communities in Canada.

Insulin resistance, as measured using the Homeostatic Model Assessment of Insulin Resistance (HOMA-IR), signals an increased risk for type 2 diabetes.19 The present analysis addresses a knowledge gap in the relevance of current anthropometric thresholds for signaling diabetes risk in the James Bay Cree, as captured by elevated HOMA-IR.


Data sources

Data were collected as part of the Multi-Community Environment and Health (Nituuchischaayihtitaau Aschii) Study,20–22 which aimed to provide assessment and surveillance of health status of the Cree in a context of hydroelectric dam projects with potential environmental impacts. Procedures were approved by Research Ethics Boards of Laval, McGill, and McMaster Universities, and the Cree Board of Health and Social Services of James Bay; community level consent for procedures was obtained through Band Council Resolutions. Individual participants also provided informed consent. Large scale hydroelectric dam projects and concerns about increasing mercury contamination led to the 1986 and 2001 Mercury Agreements between the Grand Council of the Crees (GCC) the government of Quebec, and Hydro-Québec, the government body responsible for dam construction. Part of the agreement included assessment of health impacts (eg, risk of consumption of contaminated fish, risks of avoiding fish consumption). There was a direct planning and consultation process between the Cree Health Board and the Cree communities for a comprehensive environment, and health study with the Mercury Agreement funds.

Simple random sampling without replacement was performed within each age stratum in order to build a list of participants who were then contacted by local recruiters. The overall participation rate (adults and children) was 48%. Demographic variables, and items relating to lifestyle factors and medical histories, were obtained through interviewer-administered questionnaires in either liyiyiuyimuwin or English. Medical chart reviews were conducted to obtain additional medical history and details of medication use.

Weight was measured on a beam scale and height to the nearest cm. WC to the nearest 0.5 cm was determined at the end of exhalation using a measuring tape horizontally placed between the last floating rib and the iliac crest. Hip circumference to the nearest 0.5 cm was assessed at the level of the pubic symphysis and the point of greatest posterior extension of the buttocks. Blood pressure was assessed following a 5 min rest period, using mercury sphygmomanometers and appropriately sized cuffs, in accordance with Canadian Hypertension Education Program (CHEP) guidelines. Three measurements were taken, and the final two systolic and diastolic measurements were separately averaged. Following an overnight fast, venous blood was sampled for measurements of glucose (spectrophotometric assay, Vitros 950, Vitros Chemistry, Ortho-Clinical Diagnostics, Rochester, New York, USA), insulin (immunoassay with chemiluminescent detection, Bayer Health Care, Advia Centaur), and total cholesterol, high-density lipoprotein (HDL) and triglycerides (Vitros 950 Chemistry Station, Ortho-Clinical Diagnostics, Raritan, New Jersey, USA). Low-density lipoprotein (LDL) was computed. Some participants additionally underwent an oral glucose tolerance test (OGTT, 75 g oral glucose load).

Cohort for present analysis

We retained adults (≥18 years) with data on all three anthropometric parameters (BMI, WC, WHR), and fasting glucose and insulin measurement, permitting computation of HOMA-IR. We excluded individuals with type 1 diabetes and those without a medical chart review. Medical chart review procedures were implemented only after data collection in two of the nine communities had been completed; these two communities were thus not included in the present analysis.

Statistical analyses

Statistical analyses were conducted using R (R V.3.0.2. (Copyright (C) 2013 The R Foundation for Statistical Computing). BMI was calculated by dividing the weight (kg) by the squared height (m2). WHR was calculated by dividing the WC (cm) by the hip circumference (cm).

The β cell function was estimated using HOMA-β=((20×insulin)/(glucose-3.5))%. We computed the HOMA-IR (fasting insulin (µunits/mL)×fasting glucose (mmol/L)/22.5). The constant 22.5 is derived from the product of normal plasma insulin of 5 µunits/mL and normal plasma glucose of 4.5 mmol/L.23 ,24 Thus, for an individual who has a plasma insulin of 5 µunits/mL and glucose of 4.5 mmol/L, HOMA-IR would be 1. While HOMA-IR is a continuous measure of insulin resistance, in the present analyses, ‘high’ insulin resistance was defined as having an HOMA-IR value greater than 2. To confirm that an HOMA-IR of 2 is of clinical importance in the Cree, we examined the odds of dysglycemia (prediabetes or diabetes) at different HOMA-IR thresholds through age-adjusted logistic regression analyses. In a sensitivity analysis, we defined ‘high’ insulin resistance as having an HOMA-IR greater than 2.7. A review of HOMA-IR cut-offs used to signal insulin resistance identified 10 studies, with 5 using a value close to 2, and 3 studies using a value of roughly 2.7.25

Participants were classified as having type 2 diabetes, prediabetes, or neither condition, based on chart review, use of medication, glucose levels and, when available, OGTT results. Diabetes was defined as having a fasting blood glucose of ≥7.0 mmol/L and/or by use of antihyperglycemic medication. Prediabetes was defined with either impaired fasting glucose (6.1 to 6.9 mmol/L) or impaired glucose tolerance (7.8 to 11 mmol/L 2 h after ingesting 75 g of glucose solution orally).26 Descriptive analyses included computation of mean values with SDs or numbers and proportions, as appropriate, stratified by HOMA-IR categories (<2, 2–3, 3–4, >4), separately for women and men. We compared adults retained versus excluded, in terms of the variables available.

We then assessed the test properties of existing obesity and overweight thresholds (men and women separately) in terms of identifying HOMA-IR >2. Thresholds examined included WHO criteria for obesity (BMI ≥30 kg/m2) and overweight (BMI ≥25 kg/m2), and ATP III, IDF (WC ≥80 cm for Asian women and ≥88 cm for European women; ≥90 cm for Asian men and ≥102 for European men), and WHO criteria (WHR >0.85 for women and >0.90 for men) for abdominal obesity. Test properties included sensitivity (test positive/true positive×100: ie, number above the obesity or overweight-defining threshold being examined divided by number with HOMA-IR >2, all multiplied by 100), specificity (test negative/true negative´100: ie, number below the obesity or overweight-defining threshold being examined divided by number with HOMA-IR <2, all multiplied by 100), positive likelihood ratio, PLR (sensitivity/(1−specificity)), and negative likelihood ratio, NLR ((1−sensitivity)/specificity). These were computed with 95% CIs.

In general, shifts from pretest to post-test probability are large with PLRs above 10, moderate between 5 and 10, and small but sometimes clinically important between 2 and 5.27 Similarly, shifts from pretest to post-test probability are large with NLR values below 0.1, moderate between 0.1 and 0.2, and small but sometimes important between 0.2 and 0.5.27 These principles were applied in our interpretation of likelihood ratios.

Receiver operating characteristic (ROC) curve28 ,29 analyses were used to identify the BMI, WC and WHR thresholds associated with the best possible combination of sensitivity and specificity for HOMA IR >2 in the data set (ie, value in the upper left corner of the ROC). PLR and NLR were computed for the optimal thresholds thus defined.


There are nearly 10 000 adults living in Eeyou Istchee. A total of 1405 individuals were consented for data collection procedures, which included medical chart review (figure 1). This included 860 adults. Among these adults, 16 did not have WC data, 53 did not have a chart review, 2 had type 1 diabetes, and in 6, prediabetes/diabetes status was not clear from the data available. Seven hundred and eighty-three adults were retained in the present analyses (783/860, 91%). In adults excluded, mean HOMA-IR was >6 both in women and in men, and characteristics (results not shown) were similar to those in the upper two HOMA-IR categories of those retained (ie, study cohort, table 1).

Table 1

Characteristics, by degree of hepatic insulin resistance, in women and men

The prevalence of diabetes was 12.4% (57) in men and 12.7% (41) in women. The prevalence of prediabetes was 19.4% (63) in men and 24.8% (114) in women. In the study cohort, HOMA-IR distributions in women and men had a right skew, but even among those without diabetes or prediabetes, the mean HOMA-IR was high (5 in women and 4 in men; see online supplementary figure S1A, B). The odds of prediabetes dysglycemia increased twofold above an HOMA-IR threshold of 2, both in age-adjusted and in unadjusted logistic regression analyses (see online supplementary table S1). While mean HOMA-IR was high in women and men, a lower proportion of women had HOMA-IR value of under 2 (4.6% of women, 1; 13.8% of men, table 1). HOMA-IR was higher in older individuals and there was a step-wise increase in all obesity and blood pressure measures, and in medication use, across HOMA-IR categories. While total cholesterol to HDL ratio and LDL-C differed little across groups, HDL was lower and triglycerides higher in the upper two HOMA-IR categories. HOMA-β increased with HOMA-IR.

Body mass index

In women, the WHO BMI obesity threshold of 30 kg/m2 had a sensitivity of 79% (95% CI 75% to 83%) and specificity of 86% (95% CI 64% to 97%) for HOMA-IR >2 (table 2). The PLR was 5.56 (95% CI 1.95 to 15.9) and NLR was 0.24 (95% CI 0.19 to 0.31). In men, this threshold performed similarly with a sensitivity of 70% (95% CI 64% to 75%) and a specificity of 91% (95% CI 78% to 97%), PLR of 7.51 (95% CI 2.94 to 19.2), and NLR of 0.33 (95% 0.27 to 0.41). The WHO BMI overweight threshold of 25 kg/m2 had a higher sensitivity but lower specificity, and lower PLR. The threshold with the maximum sensitivity and specificity (ROC analysis; (see online supplementary figure S2A, B) was 28.5 kg/m2. Defining insulin resistance at HOMA-IR >2.7 instead of at 2, resulted in a higher sensitivity and lower specificity for the obesity threshold of 30 kg/m2 (see online supplementary table S2); PLR values remained above 3 both in women and in men. The optimal threshold rose to 30.8 kg/m2 in women and 29.8 kg/m2 in men.

Table 2

Properties of obesity thresholds for detection of HOMA-IR >2 among Cree

Waist circumference

In women, a WC of 80 cm (ATP III and IDF/Asian) and 88 cm (IDF/European) were both highly sensitive but had low specificities and PLRs, and higher NLRs (table 2); PLR values were also low when the HOMA-IR threshold was raised to 2.7 (see online supplementary table S2). The optimal ROC-defined threshold (see online supplementary figure S2C) was higher at 98 cm (sensitivity of 87% (95% CI 83% to 90%), specificity of 81% (95% CI 58% to 95%), PLR of 4.55 (95% CI 1.88 to 11.0), and NLR of 0.16 (95% CI 0.12 to 0.23)). With an HOMA-IR threshold of 2.7, the optimal WC threshold in women was 105.5 cm (see online supplementary table S2).

In men, a WC of 94 cm (IDF/European) and 90 cm (IDF/Asian men) both had high sensitivity but low specificity (table 2). A WC of 102 cm (ATPIII) had a better balance between sensitivity (82%, 95% CI 77% to 87%) and specificity (82%, 95% CI 68% to 92%), with a PLR of 4.64 (95% CI 2.47 to 8.71) and NLR of 0.21 (95% CI 0.16 to 0.28). With an HOMA-IR threshold of 2.7, a WC of 102 cm or more in men had a PLR of 2.61 and NLR of 0.17. The optimal ROC-defined threshold for WC in men was very similar to the ATPIII threshold at 101 cm (see online supplementary figure S2D). With an HOMA-IR threshold of 2.7, the optimal threshold was higher at 105.5 cm.

Waist to hip ratio

In women (table 2), the WHO 0.85 threshold was sensitive (92%, 95% CI 89% to 95%) and moderately specific (48%, 95% CI 26% to 70%) with a PLR of 1.76 (95% CI 1.17% to 2.65%) and an NLR of 0.16 (95% CI 0.09% to 0.29%) (table 2), similar to the test characteristics of the WHO 0.90 threshold in men. PLR values were also below 2 with an HOMA-IR threshold of 2.7. The optimal ROC-defined threshold for WHR in women was 0.87 (see online supplementary figure S2E) and the optimal value in men was 0.98 (see online supplementary figure S2F). In women, the optimal ROC-defined threshold had a slightly lower sensitivity (87%, 95% CI 83% to 90%) but much higher specificity (81%, 95% CI 58% to 95%) and PLR (4.55, 95% CI 1.88. 11.0); NLR was 0.16 (95% CI 0.12 to 0.23). In men, the optimal ROC-defined threshold also had a lower sensitivity (65%, 95% CI 59% to 70%) and higher specificity (84%, 95% CI 71% to 94%) and PLR (4.17, 95% CI 2.10 to 8.28); NLR was 0.42 (95% CI 0.34 to 0.51). With an HOMA-IR threshold of 2.7, PLR values were also below 2 both in women and in men; the optimal WHR was 0.95 in women and remained at 0.98 in men.


Our analyses indicate that among the Cree, an HOMA-IR >2 is consistent with a twofold greater odds of type 2 diabetes. Compared against this risk indicator, the current WHO 30 kg/m2 BMI threshold demonstrated moderately important PLR values (ie, above 5) and somewhat important NLR values (ie, 0.2 to 0.5) in women (PLR 5.56, 95% CI 1.95 to 15.9; NLR 0.24, 95% CI 0.19 to 0.31) and men (PLR 7.51, 95% CI 2.94 to 19.2; NLR 0.33, 95% CI 0.27 to 0.41). The optimal value was slightly lower (28.5 kg/m2) than the WHO threshold. In men, the ATP III WC threshold (102 cm) had a PLR close to 5 (PLR 4.64, 95% CI 2.47 to 8.71; NLR 0.21, 95% CI 0.16 to 0.28). For WC in women, and WHR in women and men, none of the current thresholds had a PLR value above 5. The findings indicate that for diabetes risk assessment in the Cree, a BMI of above 30 kg/m2 is a risk factor to consider. In men, a WC >102 cm may also be a meaningful threshold in terms of diabetes risk. These thresholds also appeared reasonable when the HOMA-IR cut-off was raised from 2 to 2.7 in sensitivity analyses.

The analyses presented suggest that, for WC in women and WHR in women as well as in men, current thresholds are not associated with a high PLR, and our analyses suggest that optimal thresholds may be higher, although this requires confirmation. In a previous Canadian study,30 Lear and colleagues compared Aboriginal (at least three Aboriginal grandparents) and European-origin Canadians living in Vancouver. They measured anthropometric parameters, blood pressure, and cardiometabolic risk factors, and performed a single slice abdominal CT scan to assess visceral and subcutaneous adipose tissue. At the same BMI and WC, they identified no difference in subcutaneous or visceral adipose tissue areas between Aboriginal and European participants, nor in the prevalence of dyslipidemia, hypertension, impaired fasting glucose or metabolic syndrome. Insulin levels, however, were higher in Aboriginal Canadians, and thus HOMA-IR values, had they been computed, would likely have been higher at a given WC or BMI. It is thus possible that at-risk WC, BMI, and WHR levels, may differ between Aboriginal populations. Similarly, another previous Canadian study7 used factor analysis to ascertain the BMI cut-points equivalent to 30 kg/m2 in South Asian, Chinese, and Aboriginal Canadians (Six Nations, Ontario). Their analysis suggested that a much lower BMI was optimal in these groups (ie, <22 kg/m2) including the Aboriginal Canadians studied. They included a mix of Aboriginal Canadians of Mohawk, Cayuga, Onondaga, Oneida, Seneca and Tuscarora origin, living in the Six Nations community in Ontario; they did not include the Cree. In contrast, we determined the BMI cut-point of 30 kg/m2 to be applicable in terms of risk of elevated HOMA-IR in the James Bay Cree. There is a need for population-specific thresholds to define high risk in groups that may have important cultural, genetic, and epigenetic differences.

Our analyses are limited by their cross-sectional nature; in an ideal study, potential thresholds would be examined in terms of ability to predict the development of prediabetes and diabetes. Given the limitations of the available data, we opted to use an HOMA-IR above 2 as the diabetes risk indicator against which current obesity thresholds were evaluated, and performed a sensitivity analysis using an HOMA-IR of above 2.7. Although the HOMA IR has been extensively tested against other measures of insulin resistance,31 ,32 it is a surrogate marker, and there may be differences between the assays used to derive the HOMA IR equation and the assays of insulin resistance used in our study. However, at the HOMA-IR of >2, there is a clinically important increase in odds of prediabetes or diabetes, as confirmed by our logistic regression analyses (see online supplementary table S1).


Identifying participants at the stage of insulin resistance before the development of diabetes is an important window of opportunity for focused intensive lifestyle change efforts, demonstrated to be effective in large clinical trials.33–36 Anecdotally, clinicians working in Eeyou Istchee (Cree territory) have reported to us that, given the high frequency of a BMI above 30 kg/m2 among their patients, they have started questioning what is ‘normal’ in the community, and considered the possibility that a BMI above this threshold might not have the same importance here as it has in the general population. The analyses presented do suggest that Cree women may be able to tolerate a WC above current thresholds before experiencing an important increase in insulin resistance. This remains to be confirmed in longitudinal studies. A large majority of Cree in Eeyou Istchee have a BMI above 30 kg/m2. Our findings indicate that in this Aboriginal population, a BMI above 30 kg/m2 is associated with a high HOMA-IR. Further, a BMI at this level has a moderately high PLR value for HOMA-IR >2 in women and men. Thus, the WHO obesity threshold demonstrates clinical importance in this regard for this population.


The authors thank the Principal Investigators of the Multi-Community Environment and Health Longitudinal (Nituuchischaayihtitaau Aschii) Study and the Cree Board of Health and Social Services of James Bay for permitting this data analysis.


Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.


  • Contributors PM, DD, and KD formulated the study question and devised the analytic plan. PM and KD performed the literature search and summarized the literature. PM performed the data-analysis in collaboration with LJ and KD. PM and KD co-wrote the manuscript, with review and comment by DD, LJ, and JT. KD is the guarantor and, as such, had full access to all the data in the study, and takes responsibility for the integrity of the data and the accuracy of the data analysis.

  • Funding Data collection for the original study (Nituuchischaayihtitaau Aschii: Multi-Community Environment-and-Health Study in Iiyiyiu Aschii) was funded through an agreement between the Grand Council of the Crees and Hydro Quebec; it was 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.

  • Competing interests PM is the assistant professor of Department of Medicine at the University of British Columbia. DD is the assistant professor of Department of Medicine at McGill University and the medical advisor for chronic diseases with the Cree Board of Health and Social Services of James Bay. LJ is the professor of Department of Epidemiology, Biostatistics, and Occupational Health at McGill University. JT is the assistant director of Public Health Surveillance, Evaluation, Research, and Communication, and Clinical Preventive Practices with the Cree Board of Health and Social Services of James Bay. KD is the associate professor in the Department of Medicine of McGill University and holds a Senior Investigator award from the Fonds de Recherche en santé du Québec.

  • Ethics approval The procedures were approved by Research Ethics Boards of Laval, McGill, and McMaster Universities, and from Board resolution of the Cree Board of Health and Social Services of James Bay; community level consent for the study was obtained through Band Council Resolutions.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement The data are owned by the Cree people of northern Québec and the Cree First Nations, and are held by the Cree Board of Health and Social Services of James Bay. Permission to use these data must be sought from the owners of the data.