Article Text

Download PDFPDF

Neighbourhood socioeconomic status modifies the association between anxiety and depression during pregnancy and preterm birth: a Community-based Canadian cohort study
  1. Kamala Adhikari1,
  2. Scott B Patten1,
  3. Tyler Williamson1,
  4. Alka B Patel1,2,
  5. Shahirose Premji3,
  6. Suzanne Tough1,4,
  7. Nicole Letourneau1,4,5,6,
  8. Gerald Giesbrecht1,4,
  9. Amy Metcalfe1,7,8
  1. 1 Department of Community Health Sciences, University of Calgary, Calgary, Alberta, Canada
  2. 2 Applied Research and Evaluation- Primary Health Care, Alberta Health Services, Calgary, Alberta, Canada
  3. 3 School of Nursing, Faculty of Health, York University, Calgary, Alberta, Canada
  4. 4 Department of Pediatrics, University of Calgary, Calgary, Alberta, Canada
  5. 5 Faculty of Nursing, University of Calgary, Calgary, Alberta, Canada
  6. 6 Deprtment of Psychiatry, University of Calgary, Calgary, Alberta, Canada
  7. 7 Department of Obstetrics and Gynecology, University of Calgary, Calgary, Alberta, Canada
  8. 8 Department of Medicine, University of Calgary, Calgary, Alberta, Canada
  1. Correspondence to Dr Kamala Adhikari; kamala.adhikaridahal{at}ucalgary.ca

Abstract

Objective This study examined the association of anxiety alone, depression alone and the presence of both anxiety and depression with preterm birth (PTB) and further examined whether neighbourhood socioeconomic status (SES) modified this association.

Design Cohort study using individual-level data from two community-based prospective pregnancy cohort studies (All Our Families; AOF) and Alberta Pregnancy Outcomes and Nutrition (APrON) and neighbourhood SES data from the 2011 Canadian census.

Setting Calgary, Alberta, Canada.

Participants Overall, 5538 pregnant women who were <27 weeks of gestation and >15 years old were enrolled in the cohort studies between 2008 and 2012. 3341 women participated in the AOF study and 2187 women participated in the APrON study, with 231 women participated in both studies. Women who participated in both studies were only counted once.

Primary and secondary outcome measures PTB was defined as delivery prior to 37 weeks of gestation. Depression was defined as an Edinburgh Postnatal Depression Scale (EPDS) score of ≥13, anxiety was defined as an EPDS-anxiety subscale score of ≥6, and the presence of both anxiety and depression was defined as meeting both anxiety and depression definitions.

Results Overall, 7.3% of women delivered preterm infants. The presence of both anxiety and depression, but neither of these conditions alone, was significantly associated with PTB (OR 1.6, 95% CI 1.1 to 2.3) and had significant interaction with neighbourhood deprivation (p=0.004). The predicted probability of PTB for women with both anxiety and depression was 10.0%, which increased to 15.7% if they lived in the most deprived neighbourhoods and decreased to 1.4% if they lived in the least deprived neighbourhoods.

Conclusions Effects of anxiety and depression on risk of PTB differ depending on where women live. This understanding may guide the identification of women at increased risk for PTB and allocation of resources for early identification and management of anxiety and depression.

  • mental health
  • epidemiology
  • preterm birth
  • neighbourhood socioeconomic status
http://creativecommons.org/licenses/by-nc/4.0/

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, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

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.

Strengths and limitations of this study

  • This study used data from two community-based prospective pregnancy cohort studies that were conducted in a relatively representative sample (compared to a hospital-based or clinic-based sample) of pregnant women.

  • This study performed statistical analysis to examine the relationship between anxiety and depression during pregnancy and preterm birth by analysing the presence of both depressive and anxious symptoms versus isolated depressive or anxious symptoms as risk factors of preterm birth, and further analysing whether neighbourhood socioeconomic status modifies the relationship.

  • This study analysed overall preterm birth as data on spontaneous and iatrogenic preterm birth were not available.

  • The study sample over-represents women from urban areas of Alberta, with high socioeconomic status, thus limiting the generalisability of the findings to urban settings.

  • The use of self-reported anxiety and depression measurement scales may introduce measurement inaccuracy.

Background

Worldwide, a total of 15 million births occur preterm (ie, before 37 weeks of gestation), with a global average rate of 11.1%.1 Preterm birth (PTB) is responsible for 35% of neonatal deaths globally.2 Among survivors, it is also a significant risk factor for short and long-term morbidities, such as respiratory distress syndrome, cerebral palsy and learning difficulties.3–5 Despite substantial research and interventions to prevent PTB, the incidence of PTB has not declined and its aetiology remains unclear.1 6 Understanding the risk factors for PTB, such as psychosocial distress and neighbourhood low socioeconomic status (SES), may help identify women at increased risk, and assist in the allocation of resources, ultimately reducing the incidence of PTB.

PTB has been linked to psychosocial distress during pregnancy, specifically anxiety and depression—the most common mental health problems during pregnancy.7–10 However, the association between anxiety and depression during pregnancy and PTB is incompletely understood. Many previous studies on the association between anxiety and depression and PTB were conducted in medical settings (ie, hospital and clinic) with small samples and high rates of attrition.7 9 10 Notably, most of the previous studies analysed anxiety or depression without considering that they may occur in a comorbid state.7–11 Comorbid anxiety and depression is, in fact, common (affecting up to 50% of women with anxiety or depression) and is more likely to involve severe symptoms of anxiety and depression than isolated anxiety or depression.12–14 Thus, comorbid anxiety and depression may pose a higher risk of PTB than isolated anxiety or depression, which may influence the association between anxiety or depression and PTB.

Anxiety and depression are negatively correlated with neighbourhood SES.15 Neighbourhood SES is an area-level measure of SES, which aggregates individual SES (such as income, education and employment status) at a certain geographical level.16 Neighbourhood SES may influence the risk of PTB by exposing women to health benefiting or risk elevating factors, such as access to healthy foods, quality health services, opportunities for leisure activity and social support and exposure to societal stressors, crimes and poor air and water quality.16–19 Low neighbourhood SES may affect an individual’s ability to fulfil daily needs, access resources, make lifestyle choices and cope with different situations.16–19 Thus, the risk of PTB that is associated with anxiety and/or depression during pregnancy may differ by neighbourhood SES. To our knowledge, this has not been examined.

This study examined the association of the presence of anxiety symptoms alone, depression symptoms alone and both anxiety and depression symptoms with PTB. This study further examined whether the presence of anxiety, depression and both anxiety and depression interact with neighbourhood SES to increase the risk of PTB. This may help to determine the subgroups of women who are at increased risk for PTB.

Methods

Data sources

This study combined datasets from two community-based prospective pregnancy cohort studies in Alberta, Canada (n=5528). The All Our Families (AOF) cohort study recruited 3341 pregnant women and the Alberta Pregnancy Outcomes and Nutrition (APrON) cohort study recruited 2187 pregnant women, with 231 women participating in both studies. Women contributed only one pregnancy in the cohort. Women who participated in both studies were only counted once. The description and comparability of these two cohort studies is available elsewhere,20 21 and justifies combining these data sources.22 The AOF study aimed to examine maternal well-being and infant outcomes and the APrON study aimed to investigate the role of prenatal maternal nutrition on maternal mental health and infant outcomes.20 21 Briefly, each cohort study had similar inclusion criteria, sampling design (community-based, non-stratified sampling) and data collection methods.21 Both studies recruited pregnancy cohorts between 2008 and 2012 using community-based recruitment strategies (such as face-to-face recruitment in maternity clinics by research assistants or nurses and recruitment in public places using posters, pamphlets and brochures) and followed them up. The follow-up for mother and child dyad is still ongoing in both studies.21

We obtained two de-identified cohort datasets linked with neighbourhood SES data from SAGE (Secondary Analysis to Generate Evidence), the secure data repository developed by PolicyWise for Children & Families which houses these datasets.

Patient and public involvement

This study used de-identified secondary data. Patient and public were not involved in the design or planning of the study.

Variables

Variables that were deemed similar in the two studies were harmonised and appended into a single new dataset. Women who participated in both studies (n=231) were counted only once. Data on age, ethnicity (white includes all Caucasians and non-white includes all non-Caucasians), maternal SES, parity, body mass index (BMI), smoking status, social support, depression and anxiety were collected at <27 weeks of gestation (in the APrON study) and at <25 weeks of gestation (in the AOF study). BMI was calculated based on the self-reported pre-pregnancy height and weight (ie, immediately before pregnancy). Additionally, depression and anxiety were measured during the third trimester (APrON: 27–42 weeks of gestation; AOF: 34–36 weeks of gestation).20 21

Both cohorts used an identical measure of depression, that is, the Edinburgh Postnatal Depression Scale (EPDS). The EPDS is a 10-item self-reported scale with each item ranging from 0 to 3 to assess symptoms of current depression (ie, how women have felt in the past 7 days).23 The EPDS has high internal consistency of 0.87,23 a sensitivity of 78% and specificity of 99% in the obstetric population,24 25 and is the most common scale used to measure antenatal and postnatal depression.26 The recommended standard cut-off score of ≥13 out of 30 points on the EPDS was used to define the presence of clinically significant depression during pregnancy.27 While the EPDS was specifically designed to assess depression, three items (namely items 3, 4 and 5) comprising the anxiety subscale (EDPS-3A) have been suggested as a measure of anxiety by previous studies,28 29 with a sensitivity of 66.7% and specificity of 88.2% in the obstetric population.29 The standard cut-off of ≥6 out of a maximum of 9 is used to define the presence of clinically significant anxiety during pregnancy.29 The cohort studies used different measures of anxiety: the AOF study used the State-Trait Anxiety Inventory and the APrON study used the Symptoms Checklist 90. Thus, the EDPS-3A was chosen as a measure of anxiety to have a consistent measure across studies and to avoid the introduction of misclassification bias related to the use of different tools. Presence of isolated anxiety or depression was defined as meeting the anxiety or depression definition during pregnancy. Presence of both anxiety and depression was defined as meeting both anxiety and depression definitions at the same time point in pregnancy. Births that occurred before 37 weeks of gestation were defined as PTB (both spontaneous and iatrogenic included). PTB was measured at 4 months of postpartum period based on maternal recall of week of gestation at delivery.

Neighbourhood SES data were measured by the Pampalon material deprivation index (derived from the 2011 Statistics Canada census)30 31 which was aggregated at the dissemination area (DA) level. DA is the smallest geographical unit available in the Canadian census, consisting of 400–700 persons.32 The Pampalon material deprivation index is a composite measure of neighbourhood SES that combines the proportion of persons without high school diplomas, the average personal income and the rate of unemployment within the DA. It is used as a deprivation quintile, with quintile one representing the least deprived and quintile five representing the most deprived neighbourhoods.30 Neighbourhood SES was assigned to each cohort based on their postal code of residence at the time of cohort recruitment.

Data analysis

Bivariate analysis was used to identify variables associated with PTB as well as anxiety and/or depression (‘anxiety only’, ‘depression only’ and ‘both anxiety and depression’). The significantly associated (p<0.05) variables were parity, ethnicity and BMI. A multivariable logistic regression model was constructed to examine the association between anxiety and/or depression and PTB. The model also included parity, ethnicity, BMI, maternal age, smoking, social support, maternal education and household total income variables. These variables were selected to adjust for in the model based on our prior knowledge (or conceptual understanding based on literature) that they are associated with both outcome (ie, PTB) and exposure (ie, anxiety and/or depression) but do not reside in the causal pathway of the relationship between anxiety and/or depression and PTB. The underlying hypothetical relationship of the variables have been shown using a direct acyclic diagram (online supplementary figure 1).

Supplemental material

A multilevel multivariable logistic regression model, which assumes the lack of independence of observations and accounts for the variation between groups or areas, was then constructed to examine the effect modification of neighbourhood SES on the association between anxiety and/or depression and PTB. This model included interaction terms in addition to parity, ethnicity, BMI, maternal age, smoking, social support, maternal education and household total income variables. The interaction terms comprised ‘anxiety only’, ‘depression only’ and ‘both anxiety and depression’ combined with each quintile of deprivation indices. Deprivation quintile fourth and fifth were combined as there were few or no cases in some strata.

The presence of significant interactions was identified through the p values associated with beta coefficients of each interaction term. Subsequently, we constructed another model without the interaction terms. A likelihood ratio test was used to compare the goodness of model fit between those two nested models – with and without the interaction terms. Adjusted prediction of PTB (ie, predicted probability of PTB that was evaluated at the average value of covariates, parity, ethnicity, BMI, maternal age, smoking, social support, maternal education and household total income variables, across observations) was estimated using the model with interaction terms. Alpha (α) of <0.05 was used to determine statistical significance. All analyses were performed using STATA/IC V.14.1.

Missing data

The proportion of missing data for gestational age at delivery (PTB) was 7.5%, for neighbourhood deprivation indices was 7.8%, and for BMI was 6.8%. Other variables had missing data <5%, ranging from 1.3% for depression to 4.4% for household total income. The missing data for these variables occurred due to maternal non-response. Characteristics of groups of women (such as ethnicity, parity, BMI, neighbourhood SES, anxiety and depression) with and without missing data on PTB were compared to assess differences. Multiple imputation was used to address with missing data on the three variables (ie, PTB, BMI and neighbourhood deprivation indices) that had ≥5% missing data .33 Using STATA’s ‘mi Package’, the multiple imputation process was carried out in three steps as recommended by Rubin: imputation, analysis and combination.33 34 The method assumes that the missing data are missing at random and attempts to estimate a missing value within a plausible set of values.33 34 The imputation values (ie, a predictive distribution based on observed data) were estimated using an imputation model (with imputation 50 times).33 34 The imputation model included the variables that were significant with missing data (ie, marital status, duration of stay in Canada, intended pregnancy, alcohol consumption) as well as those that were used in the analysis model (ie, PTB, ethnicity, parity, BMI, maternal age, smoking, social support maternal education, household total income, neighbourhood SES, anxiety and depression, and interaction terms). Sensitivity analysis was done to compare the estimates from the analyses based on multiple imputation and from the analysis restricted to complete case.

Results

Of total 5297 pregnant women, 17.9% of women had anxiety and/or depression: 7.7% of women had both anxiety and depression, followed by 6.0% women had anxiety alone, and 4.2% women had depression alone. Women with both anxiety and depression had a higher rate of PTB (10.6%) compared to those with isolated anxiety (6.5%) or isolated depression (8.2%) or without anxiety and depression (6.9%). A higher proportion of women with a presence of both anxiety and depression (compared to those with anxiety or depression alone) were single, non-white, recent immigrants, had a low household income and were from the most deprived neighbourhoods (p<0.05) (table 1). Mean scores of anxiety (mean=6.6, SD=0.4) and depression (mean=16.2, SD=0.1) were higher among women with both conditions compared to those with anxiety alone (mean=6.1, SD=0.2) or depression alone (mean=14.6, SD=0.1). As shown in table 2, women who delivered preterm infants were more likely to be non-white, obese, primiparous and from the most deprived neighbourhoods. Variables such as maternal ethnicity, parity, BMI, neighbourhood SES, anxiety and depression were significantly associated with the presence of missing data on PTB.

Table 1

Distribution of maternal characteristics across anxiety and depression status during pregnancy

Table 2

Distribution of maternal characteristics across preterm birth status

The presence of both anxiety and depression (adjusted OR (aOR) 1.6, 95% CI 1.1 to 2.3), but neither anxiety alone (aOR 0.8, 95% CI 0.6 to 1.6) nor depression alone (aOR 1.3, 95% CI to 0.8, 2.5), was significantly associated with PTB (table 3). Effect modification was observed between the presence of both anxiety and depression and neighbourhood SES (p=0.004). Compared with women without anxiety and depression, women with both anxiety and depression who lived in quintile three and more deprived neighbourhoods had significantly increased odds of experiencing a preterm delivery (quintile 4 and 5: aOR 2.2, 95% CI 1.3 to, 4.0). In contrast, compared with women without anxiety and depression, women with both anxiety and depression who lived in the least deprived neighbourhood were not at elevated odds of experiencing a preterm delivery (aOR 0.2, 95% CI 0.1 to 1.5) (table 3). The OR estimates from the analyses based on multiple imputation and from the analysis restricted to complete cases were similar, with some confidence intervals being slightly narrower in the multiple imputation analysis.

Table 3

Association between anxiety and depression status during pregnancy and preterm birth*

As shown in table 4, the predicted probability of PTB for women with a presence of both anxiety and depression was 10.0% (95% CI 6.8% to 13.1%). It increased to 15.7% (95% CI 9.5% to 23.2%) if they lived in the most deprived neighbourhoods—an increase of 57.1%—and it decreased to 1.4% (95% CI 0.1% to 4.2%) if they lived in the least deprived neighbourhoods. The predicted probability of PTB for women with depression alone was 9.6% (95% CI 5.2% to 14.1%), which increased to 14.1% (95% CI 2.7% to 25.3%) if they lived in the most deprived neighbourhoods. The predicted probability for women with anxiety alone and women with absence of anxiety and depression remained similar across the neighbourhood deprivation indices.

Table 4

Predicted marginal prevalence of preterm birth*

Discussion

Main findings

This study examined the association of anxiety alone, depression alone and the presence of both anxiety and depression during pregnancy with PTB, using data from two community-based pregnancy cohort studies in Alberta, Canada. The study found an association between the presence of both anxiety and depression and PTB which significantly differed according to neighbourhood SES. Women with both anxiety and depression were more likely to deliver preterm infants if they lived in a relatively more deprived neighbourhood compared to if they lived in a less deprived neighbourhood. For women with both anxiety and depression, the absolute predicted probability of delivering preterm infants was 16% if these women lived in the most deprived neighbourhood, and it was 1% if they lived in the least deprived neighbourhood. Overall, the findings suggest the importance of neighbourhoods on maternal health (in general) and more specifically PTB.

Interpretation

Although few previous studies assessed the association between the presence of both anxiety and depression during pregnancy and PTB, our finding is consistent with their findings that the presence of both anxiety and depression increases the likelihood of PTB.14 35 36 This may be related to the additive effects of prenatal depression and anxiety and the effects of severity of anxiety and depressive symptoms. Previous studies conducted in the general population and in pregnant women found a higher score of anxiety or depression symptoms among those with both anxiety and depression than those with isolated anxiety or depression.36 37 It is also reported in previous studies that individuals with both anxiety and depression have longer depressive episodes, worse psychosocial impairment, poorer response to medication, compromised quality of life and increased suicidality than those with isolated anxiety or depression.12 35 37 Thus, the presence of both anxiety and depression during pregnancy may lead to an increased risk of poor birth outcomes, including PTB, relative to depression or anxiety alone.

Our study did not find an association between anxiety alone or depression alone and PTB, which is consistent with a previous pregnancy cohort study that analysed isolated anxiety or depression separately from the presence of both or comorbid anxiety and depression.36 However, the finding is inconsistent with several previous studies that analysed anxiety or depression intermixing with the presence of both conditions.8 10 It is possible that the association described in the literature requires high levels of anxiety or depression, which is more likely present in the presence of both anxiety and depression symptoms or disorders. Thus, the associations found in previous studies may have been confounded by the presence of both anxiety and depression symptoms or comorbid anxiety and depression disorders. The increased risk of PTB associated with the presence of both anxiety and depression (but not with isolated anxiety or depression) may, in part, explain the inconsistencies across previous findings on the association between prenatal anxiety or depression and PTB. Similarly, previous studies did not analyse the association stratified by neighbourhood SES, meaning that these studies averaged the association across neighbourhood SES, which may also explain the inconsistencies across previous studies findings.

A strong association between the presence of both anxiety and depression and PTB among women living in a relatively more deprived neighbourhood may reflect that, besides individual level risk factors, PTB is related to neighbourhood factors.16–18 For example, women living in deprived neighbourhoods often have less access to healthy foods, quality health services and opportunities for leisure activity, and have more exposure to societal stressors and crimes.16–19 Anxious and depressed women living in less advantaged areas may interpret the deprivation associated stressors more acutely and have less support or are less able to manage or cope with their stressors, making them severely emotionally distressed compared to those living in more advantaged areas.8 11 38 39 Consequently, the elevated risk of delivering preterm is more likely to occur in this group of women. However, it is important to note that, the relationship between mental illness and impoverishment is difficult to interpret as causal, given the bi-directional relationship between them. Furthermore, in our study, the group of women with both anxiety and depression (who often have severe symptoms of anxiety or depression) in the least deprived neighbourhoods had an exceptionally low rate of PTB. The observed association between the presence of both anxiety and depression and PTB among women living in a relatively more deprived neighbourhood seems to depend on this result. Thus, the replication of this finding seems important.

Strengths and limitations

To our knowledge, few studies have directly examined the presence of both depressive and anxious symptoms versus isolated depressive or anxious symptoms as risk factors of PTB, and no studies have examined neighbourhood SES as a modifier to the relationship between anxiety and/or depression and PTB. This study is important given its focus on the most common psychological condition (ie, comorbid anxiety and depression) and the importance of identification of specific groups of women who may benefit the most from the preventive interventions. This study used two community-based prospective pregnancy cohort studies. This provided an opportunity to describe PTB across the several strata of anxiety, depression and both anxiety and depression and neighbourhood SES in a relatively representative sample (compared, for example, to a hospital-based or clinic-based sample) of pregnant women. However, even using the two cohorts, some strata had few cases of preterm infants, which may have led to the observed imprecise and/or insignificant estimates (specifically in a group with depression alone). As these cohorts over-represent women with high SES,21 40 41 it limits the generalisability of the findings to other demographic groups. While the use of prospective measurement of depression and anxiety reduces the chance of misclassifications due to recall bias, the use of self-reported anxiety and depression measurement scales may have introduced measurement inaccuracy. Specifically, the EPDS-3A scale has not been validated in a pregnant population and it tends to provide high false-positive results based on its validation on during the postpartum period.28 29 Furthermore, the EPDS-3A is a subscale of the EPDS. The standard cut-off point for the EPDS excluding the items of the EPDS-3A has not been established. While the use of a single scale may overestimate the presence of anxiety and/or depression, being able to identify combined anxiety and depression group using a single scale is advantageous as it facilitates for intervention design. While we examined the association between anxiety and/or depression and PTB analysing the influence of several potential confounders, other confounders such as antidepressant use, other psychiatric conditions and medical risk factors that may influence the associations were not considered since they were not available in the study’s data sources. Similarly, we were not able to separate out spontaneous and iatrogenic PTB in the model – the association might be stronger for spontaneous PTB. Overall, replication of this study addressing these limitations may further the understanding on risk factors and preventive strategies of PTB.

We defined neighbourhoods using the smallest area (ie, DA) where people living in the smallest area are more likely to be similar for the outcomes, used multilevel analysis that accounts for area-level variation, and adjusted for individual level variables, an appropriate analytical approach for multilevel data. However, it is difficult to interpret the influence of neighbourhood SES using area-based variables, where women living in the same area share the same value for the variable. Individuals who live in the same area may also experience different contextual influences from many other areal units, and the timing and duration in which individuals experienced these contextual influences is also uncertain.

Conclusions

Our study found that the presence of both prenatal anxiety and depression increases the likelihood of PTB and the effect of this combination is stronger for women living in low SES neighbourhoods compared to women living in high SES neighbourhoods. The finding may help to inform development of intervention strategies (such as timely screening and management of anxiety and depression) that focus on the most deprived neighbourhood. Furthermore, future research that examines the influence of severity of anxiety and depression on risk of PTB may further the understanding on risk factors and preventive strategies of PTB.

Acknowledgments

KA was awarded the Vanier Canada Graduate Scholarship and the Alberta Innovates Graduate Studentship Award. Amy Metcalfe is supported by a Canadian Institutes of Health Research New Investigator Award. We acknowledge the All Our Families and the Alberta Pregnancy Outcomes and Nutrition cohort study teams for providing permission to use their data. We acknowledge SAGE (Secondary Analysis to Generate Evidence), the secure data repository developed by PolicyWise for Children and Families, which houses these datasets, for providing access to these datasets.

References

Footnotes

  • Contributors KA involved in the conception and design of the study. KA is also responsible for conducting the analysis, interpreting the data and drafting the manuscript. AM provided overall supervision to Kamala in conducting this study and contributed to conception and study design, interpretation of data, provided intellectual content and revisions to manuscript. SP, TW, ABP, SP, ST, NL and GG were involved in the conception and design of the study and provided interpretation and intellectual content to subsequent drafts of the manuscript. All authors read and approved the final draft.

  • Funding KA was awarded the Vanier Canada Graduate Scholarship (Award code: 201611CGV-3 82 013–2 67 341) by the Canadian Institutes of Health Research and the Alberta Innovates Graduate Studentship Award (Award code: 201611CGV-3 82 013–2 67 341) by the Alberta Innovates to conduct this study.

  • Disclaimer The funding agencies have no role in study design, data access, analysis and interpretation; manuscript writing and in the decision to submit the article for publication.

  • Competing interests None declared.

  • Patient consent for publication Not required.

  • Ethics approval Ethics approval for this study was obtained from the Conjoint Health Research Ethics Board at the University of Calgary (REB16-2548_REN1). This study used secondary data and all the data were anonymised; therefore, did not require informed consent.

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

  • Data availability statement Data may be obtained from a third party and are not publicly available. Additional data such as statistical codes, supplementary tables and technical appendix are available upon request (by emailing Kamala Adhikari: kamala.adhikaridahal@ucalgary.ca)