Skip to main content

Variability in intensive care unit admission among pregnant and postpartum women in Canada: a nationwide population-based observational study

Abstract

Background

Pregnancy-related critical illness results in approximately 300,000 deaths globally each year. The objective was to describe the variation in ICU admission and the contribution of patient- and hospital-based factors in ICU admission among acute care hospitals for pregnant and postpartum women in Canada.

Methods

A nationwide cohort study between 2004 and 2015, comprising all pregnant or postpartum women admitted to Canadian hospitals. The primary outcome was ICU admission. Secondary outcomes were severe maternal morbidity (a potentially life-threatening condition) and maternal death (during and within 6 weeks after pregnancy). The proportion of total variability in ICU admission rates due to the differences among hospitals was described using the median odds ratio from multi-level logistic regression models, adjusting for individual hospital clusters.

Results

There were 3,157,248 identifiable pregnancies among women admitted to 342 Canadian hospitals. The overall ICU admission rate was 3.2 per 1000 pregnancies. The rate of severe maternal morbidity was 15.8 per 1000 pregnancies, of which 10% of women were admitted to an ICU. The most common severe maternal morbidity events included postpartum hemorrhage (n = 16,364, 0.52%) and sepsis (n = 11,557, 0.37%). Of the 195 maternal deaths (6.2 per 100,000 pregnancies), only 130 (67%) were admitted to ICUs. Patients dying in hospital, without admission to ICU, included those with cardiovascular compromise, hemorrhage, and sepsis. For 2 pregnant women with similar characteristics at different hospitals, the average (median) odds of being admitted to ICU was 1.92 in 1 hospital compared to another. Hospitals admitting the fewest number of pregnant patients had the highest incidence of severe maternal morbidity and mortality. Patient-level factors associated with ICU admission were maternal comorbidity index (OR 1.88 per 1 unit increase, 95%CI 1.86–1.99), urban residence (OR 1.09, 95%CI 1.02–1.16), and residing at the lowest income quintile (OR 1.44, 95%CI 1.34–1.55).

Conclusions

Most women who experience severe maternal morbidity are not admitted to an ICU. There exists a wide hospital-level variability in ICU admission, with patients living in urban locations and patients of lowest income levels most likely to be admitted to ICU. Cardiovascular compromise, hemorrhage, and sepsis represent an opportunity for improved patient care and outcomes.

Introduction

There are approximately 300,000 maternal deaths in the world annually [1], and in a recent population-based study, we have highlighted 18 episodes of severe maternal morbidity (i.e., a potentially life-threatening condition) per 1000 deliveries in Canada, increasing 1.3% annually for the last one decade [2]. However, severe maternal morbidity can often be anticipated according to known maternal risk factors, and outcomes possibly improved with early recognition and timely supportive acute and intensive care [3,4,5,6,7,8].

Variability in both recognition of critical illness and access to intensive care units (ICUs) influences the risk of death in many critical care syndromes [9, 10]. For example, sociodemographic-based variation and health care provider volume of experience have been shown to be important determinants of outcome for patients with acute lung injury, ischemic cardiac disease, and traumatic injury [11,12,13,14]. We therefore aimed to describe variation in care for acutely ill pregnant and postpartum women in Canada—a country with a vast geography and provincially administered health care—and to investigate the relationship between patient factors and hospital-based variability in ICU admission practices and clinical outcomes.

Methods

Study design, data sources, and population

We conducted a nationwide population-based cohort study in Canada spanning the years from 2004 to 2015, using data from the Discharge Abstract Database (DAD) of Canadian Institute of Health Information (CIHI), which includes information from every admission to acute care hospitals across Canada, except the province of Quebec, which maintains a separate database.

Patients in this cohort were hospitalized pregnant and postpartum women within 6 weeks after the end of pregnancy, including hospitalized women with abortion or ectopic pregnancy. Available variables from the DAD included the following: age, DAD-derived variables (socioeconomic status, comorbid conditions as classified by the International Classification of Diseases and Related Health Problems-Tenth Revision, Canada [ICD-10-CA], interventions and procedures occurring in hospital as determined by the Canadian Classification of Health Interventions [CCI]), pregnancy-specific variables, hospital-level variables, and patient outcomes such as ICU admission, vital status at hospital discharge, and events allowing determination of severe maternal morbidity (Additional file 1: Table S1).

The dataset (the DAD of CIHI) has been well-validated in a number of prior studies [15,16,17,18,19,20]. In addition, the institute annually reports data quality documentation, to have transparency of data quality. An episode of pregnancy was created as previously described [2], so that all related admissions over one pregnancy course could be summarized into one pregnancy episode.

Outcomes

The primary outcome was admission to ICU which included all admissions to CIHI-derived intensive (special) care units, including “step-up”/“step-down” ICUs [21]. Secondary outcomes included a modified Canadian definition of severe maternal morbidity and maternal death [22, 23]. In order to improve the specificity of the outcome of severe maternal morbidity for this analysis, the definition was modified to exclude events consisting only of a blood transfusion and diagnoses reflecting only pre-existing conditions such as HIV, chronic congestive failure, chronic hypertension, and those with an admission-related (as opposed to outcome-related) specific diagnosis [e.g., asthma, placenta previa, severe pre-eclampsia, or sickle cell anemia] [2, 8, 22]. Maternal death was defined as death occurring during pregnancy or within 6 weeks after the end of the pregnancy.

Predictors

A hospital-level predictor was created to represent pregnancy-related admission volumes at each hospital. First, hospitals with fewer than 55 deliveries over the study period (approximately corresponding to fewer than 5 deliveries per year) were excluded using a commonly applied principle of limiting analyses and inferences about data containing fewer than 5 observations in a particular time period (in this case, 1 year). Next, hospital quintiles were created according to pregnancy admission volume at each hospital. Other ordinal predictors (e.g., 5 groups according to ICU admission rate and 5 groups according to the annual number of pregnancy admissions at each hospital) were created in order to perform hospital-related sensitivity analyses.

Comorbid conditions necessary for determining the maternal comorbidity index (MCI) were recorded for each patient, as we have previously determined it to be the most valid risk adjustment tool using health administrative data for hospitalized pregnant and postpartum women [24, 25]. The index ranges from 0 to 45, excluding the age-related variable. Age was removed from the index because the age of the index does not include age younger than 35 years, but age and the index were separately entered into the model. Since six conditions of the maternal comorbidity index had potential overlap with severe maternal morbidity outcome conditions, only the pre-existing conditions were included in the maternal comorbidity index, and those occurring after pregnancy onset were considered as outcomes for severe maternal morbidity, to ensure a valid temporal relation [2].

Statistical analyses

We used the mean and standard deviation or median and interquartile range to summarize continuous variables, and counts and proportions to describe categorical variables, as appropriate. Standardized differences were employed to compare those who required ICU admission with those who did not, where 0.2, 0.4, and 0.6 standardized differences were generally equated to small, medium, and large effect sizes, respectively [26]. All variables were screened on the basis of standardized difference and selected for regression models if the standardized difference was greater than 0.1, in addition, to an a priori decision to include the urban or rural location of patient residence and hospital location. To develop random intercept logistic regression models with patient and hospital levels, first, we explored an intercept-only model predicting the outcome of ICU admission with each hospital considered as a cluster (representing differences in ICU admission rates among various hospitals). Then, we added patient-level variables to the intercept-only model and finally hospital-level variables [27, 28].

We used the median odds ratio (MOR) to represent the variability of the outcome (e.g., ICU admission) among the hospital clusters of interest [27, 29]. For example, when the MOR = 2, it means that for two pregnant women with the same characteristics, the odds of being admitted to ICU may be twofold higher in one randomly selected hospital compared to another randomly selected hospital from a different cluster when the clusters were ordered by risk. The MOR was computed from the following equation: \( \mathrm{MOR}=\exp \Big(0.675\times \sqrt{2\times \mathrm{hospital}-\mathrm{level}\kern0.34em \mathrm{variance}} \), where the hospital-level variance of the parameter estimate (i.e., β) is derived from the random intercept model. Interval odds ratio-80 (IOR-80) was used to quantify the impact of a variable at the hospital cluster-level on the outcome (e.g., ICU admission) by incorporating between-cluster variability [27, 29]. IOR-80 covers 80% of the distribution of the odds when comparing two subjects with the same value for all other covariates, but in two different hospitals where hospital-level characteristics may differ. We calculated the lower and upper bounds of the IOR-80 from the following equations: \( {\mathrm{IOR}}_{\mathrm{lower}}=\exp \left(\beta -1.2816\times \sqrt{2\times \mathrm{hospital}-\mathrm{level}\kern0.17em \mathrm{variance}}\right) \) and \( \mathrm{IO}{\mathrm{R}}_{\mathrm{upper}}=\exp \left(\beta +1.2816\times \sqrt{2\times \mathrm{hospital}-\mathrm{level}\kern0.17em \mathrm{variance}}\right) \), where β is the regression coefficient for the hospital level. Note that the IOR does not describe the interval around a MOR; instead, the IOR reflects the size and variation of the odds ratios at the hospital level. If the IOR-80 is wide, it implies high variability in the outcome among hospitals compared to the effect of the cluster level variable. If the interval of IOR-80 crosses 1, it infers that a variable at the hospital cluster level is likely not significantly associated with the outcome. Finally, the variance partitioning coefficient (VPC) represents the proportion of total variability in the outcome (e.g., ICU admission rates) due to the differences among hospital clusters. For multi-level logistic regression models, VPC is computed from the following equation: VPC = hospital-level variance/(hospital-level variance + 3.29) [30].

Sensitivity analyses

Two pre-planned sensitivity analyses were conducted. First, we explored the association of hospital volume of pregnancy-related admissions as a hospital-level ordinal variable comprised of five categories, on the outcome of ICU admission. In the second sensitivity analysis, we restricted the primary outcome of ICU admission to only admissions to the highest acuity ICUs, excluding step-up/down units.

The third, post hoc, sensitivity analysis was performed where all analyses were redone using multiple imputation datasets, given the rate of missing data in variables reflecting socioeconomic status (Table 1). Missing variables were assessed based on an algorithm reported elsewhere, to determine characteristics of missingness (i.e., missing completely at random, missing at random, or missing not at random) [31]. We used listwise deletion (entire record excluded from analysis if any single value is missing) for the primary analysis. This method is valid if the data is missing completely at random [32]. As a sensitivity analysis, to investigate the influence of data missing not completely at random, we performed multiple imputation using the fully conditional method with 20 imputed datasets [33, 34]. We imputed missing variables using a fully conditional specification and the discriminant function method to impute the categorical variables: income quintile (14.0% missing), residence rurality (0.7% missing), and hospital rurality (0.6% missing). We created the 20 multiply imputed datasets as this generally provides good efficiency of the parameters (i.e., smaller standard error of the point estimates), in addition to sufficient power to detect differences [35,36,37].

Table 1 General characteristics of pregnant and postpartum women according to whether admitted to ICU or not

All analyses were performed using SAS statistical software, version 9.4 (SAS Institute Inc., Cary, NC) and Excel for Macintosh, version 15.3.9 (Microsoft Corporation, Redmond, WA). All analyses were completed on September 13, 2018. This study was conducted under the data security and privacy policy of CIHI and approved by the institutional review board at Mount Sinai Hospital in Toronto, Canada (10-0305-C) and University of Toronto (#34468).

Results

Incidence and clinical characteristics of pregnancy-related ICU admissions in Canada

During the study period, there were 3,162,303 pregnancies among 2,035,453 mothers, resulting in 10,204 ICU admissions (3.2 cases per 1000 pregnancies) with regional variation (Fig. 1, Table 1). There were 50,012 severe maternal morbidity events (15.8 women with a severe maternal morbidity event per 1000 pregnancies), with 5017 (0.16% overall and 10.0% of those with severe maternal morbidity) pregnant women admitted to ICU (Tables 1 and 2.). A minority of patients experiencing severe maternal morbidity were admitted to ICU (Additional file 2: Table S2, Additional file 3: Table S3). During the study period, there were 195 maternal deaths (6.2 per 100,000 pregnancies (maternal mortality rate 0.01%)); however, only 130 (67%) of these patients were admitted to an ICU prior to death (Table 2). After excluding hospitals with, on average, fewer than 5 deliveries per year, the study cohort consisted of 3,157,248 pregnancies among 342 hospitals (Additional file 4: Table S4, Additional file 5: Table S5).

Fig. 1
figure 1

Incidence of ICU admission by province in Canada. Figure legend: territories cover Yukon, Northwest Territories, and Nunavut

Table 2 Severe maternal morbidity diagnoses: all patients with severe maternal morbidity, those admitted to ICU, and those dying

Patients admitted to an ICU had a higher pre-existing maternal comorbidity index and a resulting greater burden of severe maternal morbidity (Table 1). Among all severe maternal morbidity events, postpartum hemorrhage requiring transfusion and sepsis were the most common (Table 2). Among patients admitted with acute respiratory distress syndrome, cerebral edema or coma, eclampsia, or sepsis, only 76.2%, 57.8%, 21.5%, and 4.7%, respectively, were admitted to ICU. Patients dying in hospital, without admission to ICU, included among others those with some form of cardiovascular compromise (56 patients: obstetric embolism, myocardial insufficiency, or shock), hemorrhage (8 patients), and sepsis (3 patients) (Additional file 6: Table S6).

Hospital volume of pregnancy admission

We created 5 hospital groupings according to the volume of pregnancy admissions during the study period. Hospitals admitting the fewest number of pregnant patients had the highest incidence of severe maternal morbidity (20 cases per 1000 pregnancies) and had the highest maternal mortality (40 cases per 100,000 pregnancies) (Additional file 4: Table S4), despite having similar pre-existing comorbidities. Hospitals admitting the fewest number of pregnant patients were also more likely to be in rural areas and had a higher proportion of patients in lower income quintiles. There were geographic differences in the proportion of pregnancies admitted at low-to-high-volume hospitals, with Ontario having the greatest proportion of pregnant women admitted to high-volume hospitals (Additional file 4: Table S4, Additional file 5: Table S5).

Multi-level logistic regression models

The following were entered into the multi-level logistic regression models: patient-level variables (age, maternal comorbidity index, parity, residence location, transfer status for admission) and hospital-level variables (pregnancy-related admission volume, province and territories, hospital location) (Table 3, Additional file 1: Table S1, Additional file 7: Table S7). Independent patient-level factors associated with ICU admission included maternal comorbidity index (1.88, 95%CI 1.86–1.99), urban residence (1.09, 95%CI 1.02–1.16), and lower incomes. The lower income they had, the higher risk of admitting to ICU (i.e., the lowest income had the highest risk of ICU admission: 1.43 with 95%CI 1.33–1.54). Hospital-level factors associated with ICU admission included specific provinces (Table 3, Additional file 1: Table S1, Additional file 7: Table S7).

Table 3 Estimated odds ratios for multi-level logistic regression model for intensive care unit admission

Variability

Although the overall ICU admission rate was 3.2 cases per 1000 pregnancies, the range in ICU admission variability among all hospitals was wide—from 0 to 500 cases per 1000 pregnancies (Additional file 4: Table S4). There was a minimal variability in ICU admission rate at the provincial level (Additional file 5: Table S5). Using the median odds ratio, for 2 pregnant women with similar characteristics at different hospitals, the average (median) odds of being admitted to ICU was 1.92 in 1 hospital compared to another (Additional file 7: Table S7). The proportion of variability in ICU admission due to hospital-level factors, as denoted by the variance partitioning coefficient, was 33.3%, with substantial variation in the odds of ICU admission across hospitals (Additional file 7: Table S7). Furthermore, ICU admission rate for postpartum readmissions was relatively high (36.5 per 1000 pregnancies) in comparison with other periods of admissions, and ICU admission among abortion-related admissions had the highest variability (median odds ratio, 2.21) (Additional file 8: Table S8).

Sensitivity analyses

All primary and secondary results were similar in sensitivity analyses, employing the number of pregnancy-related hospital admission in continuous variable (Additional file 9: Table S9), considering only the highest acuity ICUs (Additional file 10: Table S10), and using multiple imputation datasets to account for missing data (Additional file 11: Table S11).

Discussion

In this study, we found that only 10% of the 15.8 women per 1000 pregnancies in Canada who experience a severe maternal morbidity event and only two thirds of the 6.2 women who die per 100,000 pregnancies are admitted to ICU. We also found wide hospital-level variability in ICU admission practices. Patients living in urban settings and of lower socioeconomic status were more likely to be admitted to ICU (Table 4).

Table 4 Box summary of key points

Discovering variability in clinical practice for pregnant and postpartum women is important as it presents an opportunity for improved quality of care. For example, identifying large variability in preoperative medical consultation has led to efforts to improve perioperative anesthetic care [38]. Systematic literature search has revealed that ours is the first population-based study to examine the variability of the care of critically ill pregnant and postpartum Canadian women and highlights modifiable factors in ICU admission practices (Table 5, Additional file 12: Table S12). Prior population-based studies of pregnancy-related ICU admissions in addition to ours show variability in ICU admission rates among identified countries, ranging from 3.2 to 39 per 1000 pregnancies (Table 5) [3, 39,40,41,42]. While greater maternal comorbidity was associated with increasing likelihood of ICU admission, a minority of patients with sepsis and eclampsia were admitted to ICU. It may be due to early recovery from a sick condition after the delivery in case delivery was the first step of management and/or different pathway of patient care for pregnant women. In other words, first being admitted to obstetrical general wards may delay an acute care. Although only two thirds of all women who subsequently died were admitted to ICU, some conditions of severe maternal morbidity (e.g., obstetric embolism) and subsequent death could occur very quickly and limit the potential for ICU admission (Additional file 6: Table S6), which were indicated by an additional analysis of maternal death without ICU admission (Additional file 13: Table S13). Maternal death without ICU admission happened to those who were less sick (i.e., less maternal comorbidity index) at the time of hospital admission and less likely to transfer to another institution for acute care, but more likely to develop a life-threatening condition (i.e., severe maternal morbidity) in urban hospitals, compared to maternal death with ICU admission. There were 15 maternal deaths in operating rooms without ICU admission, related to postpartum hemorrhage and obstetric embolism. On the other hand, our results found that other conditions that progress more slowly may inherently have a greater time period for clinicians to consider admission to ICU, and hence a greater potential for intervention and improved outcomes (e.g., sepsis) through quality improvements of implementation of sepsis bundles, for example [43]. These findings highlight a potential opportunity for earlier recognition by clinicians for women at risk and earlier recognition of the onset of critical illness. More practically, in comparison with MCI = 0, increasing maternal comorbidity index has higher odds of admitting to ICU (index = 1, odds ratio (OR) 1.90, 95%CI 1.79–2.02; index = 2 or higher, OR 8.67, 95%CI 8.24–9.12) (Additional file 14: Table S14). Our findings also have relevance to policymakers. Women living in rural and resource-limited locations may experience the greatest risk [2], highlighting additional opportunities for identification of those at risk, with potential solutions of a lower threshold for consultation, referral and increased frequency of follow-up care, outreach of tele-medicine [44, 45], ease of transportation for mothers to specialist care, and/or consideration of a mobile ICU to resource-limited hospitals [46].

Table 5 Research in context. Population-based studies of pregnancy-related ICU admissions

Strengths and weaknesses of the study

Our study has several strengths. First, this is the first study based upon a nationwide cohort to explore hospital variability of pregnancy-related ICU admission using multi-level logistic regression models accounting for hospitals as clusters. Second, we adopted the perspective of an episode of care, to account for the entire pregnancy and postpartum course rather than one admission allowing use to determine outcomes from the entire pregnancy instead of outcomes relating to individual hospital admissions. Third, the maternal comorbidity index was well-validated risk adjustment tool in the DAD of CIHI, to account for known risk factors of maternal death and organ injury [24, 25, 47]. Fourth, maternal age younger than 20 years old and over 35 years old was associated with maternal mortality, and hence properly adjusted in our analyses [8]. Fifth, we employed the best risk adjustment tool for our dataset among existing risk prediction models including ICU admission scores [24, 48] and ensured a valid temporal relationship between the MCI and severe maternal morbidity/ICU admission. Last, we have used existing and well-validated health administrative datasets to describe pregnancies as well as ICU practice across diverse regions and over time [15,16,17,18,19,20].

This study also has limitations. First, certain variables that may be influential for maternal outcomes were unavailable in our dataset, for example, patient characteristics such as body mass index and pregnancy-associated factors such as the use of assisted reproductive technology, among others. However, obesity did not remain in the maternal comorbidity index to account for maternal death and organ injury during the original development [47]. Even among available variables, and considering patient and hospital factors, our predictive models for ICU admission do not explain all variability, highlighting the existence of other potentially important factors. We had a high rate of missing data for certain variables, particularly postal code markers of socioeconomic status. However, there were few differences among known characteristics between groups with missing and non-missing data, and our results did not differ when we used a multiply imputed dataset in our sensitivity analyses. CIHI did not provide structural characteristics for each of all 342 hospitals enrolled in this study such as detailed type of ICUs and existence of neonatal intensive care unit. And hence, groups according to the annual volume of their delivery admission were created, to characterize enrolled hospitals from an obstetric practice point of view. Additionally, ICU utilization and transfer rate in each group of hospitals were also calculated (Additional file 4: Table S4). It is possible that some diagnostic coding error exists, for example, over-coding of ARDS among patients managed in post-anesthesia care unit or labor and delivery room and not the ICU (only 76.2% of patients with ARDS were admitted to ICU). Canadian Coding Standards Evolution Chronicle reports that post-procedural respiratory disorder was used to be coded as ARDS [49]. Over-coding of ARDS may be attributed to this previous coding procedure. Probably, contemporary validation studies of ARDS by utilizing audits of hospital chart is required for future large epidemiological studies, since there was no such specific validation study of ARDS using CIHI DAD after the Berlin definition in 2012. Although the DAD does not contain specific information about the number of ICU beds in each acute care hospital, a previous study estimated the number of ICU beds per 100,000 population, by provinces, in 2009, ranging from 5.5 to 19.3 beds per 100,000 population [50]. Given the provincial ICU beds per 100,000 population in Canada, our study showed that there was a mostly linear relationship between provincial ICU capacity and maternal ICU admission rate (r2 = 0.74, p = 0.002) (Additional file 15: Figure S1), which implies that ICU capacity may influence the admission patterns of pregnancy-related ICU admission as well in Canada. Also, in this analysis, we did not explore neonatal outcomes—this will be an important next step. Finally, although CIHI collects data on each admission at all acute care hospitals across Canada, our findings may not generalize to the province of Quebec which maintains a unique database.

Unanswered questions and future research

We found that only a minority of women who experience severe maternal morbidity are admitted to an intensive care unit in Canada, that there is a wide variability in ICU admission practice across hospitals and provinces, and that the greatest vulnerability exists for women at low socioeconomic status. These represent potentially modifiable factors for patients and should be the focus of quality of care improvement initiatives for clinicians and health systems.

Conclusions

In conclusion of this study, we found that a minority of women who experience a severe maternal morbidity event and only two thirds of those who die are admitted to ICU. There is a wide hospital-level variability in ICU admission practices and that those living in urban locations and those of the lowest socioeconomic status were most likely to be admitted to ICU.

Availability of data and materials

All data and material can be accessed through the manuscript, supplemental information, and direct contact to the corresponding author, if needed.

Abbreviations

CCI:

Canadian Classification of Health Interventions

CI:

Confidence intervals

CIHI:

Canadian Institute of Health Information

DAD:

Discharge Abstract Database

ICD-10CA:

International Classification of Diseases and Related Health Problems Tenth Revision, Canada

ICU:

Intensive care unit

IOR-80:

Interval odds ratio-80

MCI:

Maternal comorbidity index

MOR:

Median odds ratio

OR:

Odds ratio

VPC:

Variance partitioning coefficient

References

  1. Kassebaum NJ, Bertozzi-Villa A, Coggeshall MS, Shackelford KA, Steiner C, Heuton KR, Gonzalez-Medina D, Barber R, Huynh C, Dicker D, Templin T, Wolock TM, Ozgoren AA, Abd-Allah F, Abera SF, Abubakar I, Achoki T, Adelekan A, Ademi Z, Adou AK, Adsuar JC, Agardh EE, Akena D, Alasfoor D, Alemu ZA, Alfonso-Cristancho R, Alhabib S, Ali R, Al Kahbouri MJ, Alla F, et al. Global, regional, and national levels and causes of maternal mortality during 1990-2013: a systematic analysis for the global burden of disease study 2013. Lancet. 2014;384:980–1004.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Aoyama K, Ray JG, Pinto R, Hill A, Scales DC, Lapinsky SE, Fowler RA. Temporal variations in incidence and outcomes of critical illness among pregnant and postpartum women in Canada: a population-based observational study. J Obstet Gynaecol Can. 2019;41:631–40.

    Article  PubMed  Google Scholar 

  3. Wanderer JP, Leffert LR, Mhyre JM, Kuklina EV, Callaghan WM, Bateman BT. Epidemiology of obstetric-related ICU admissions in Maryland: 1999-2008*. Crit Care Med. 2013;41:1844–52.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Pollock W, Rose L, Dennis C. Pregnant and postpartum admissions to the intensive care unit: a systematic review. Intensive Care Med. 2010;36:1465–74.

    Article  PubMed  Google Scholar 

  5. Einav S, Bromiker R, Sela HY. Maternal critical illness. Curr Anesthesiol Rep. 2017;7:55–66.

    Article  Google Scholar 

  6. Souza JP, Gülmezoglu AM, Vogel J, Carroli G, Lumbiganon P, Qureshi Z, Costa MJ, Fawole B, Mugerwa Y, Nafiou I, Neves I, Wolomby-Molondo JJ, Bang HT, Cheang K, Chuyun K, Jayaratne K, Jayathilaka CA, Mazhar SB, Mori R, Mustafa ML, Pathak LR, Perera D, Rathavy T, Recidoro Z, Roy M, Ruyan P, Shrestha N, Taneepanichsku S, Tien NV, Ganchimeg T, et al. Moving beyond essential interventions for reduction of maternal mortality (the WHO Multicountry Survey on Maternal and Newborn Health): a cross-sectional study. Lancet. 2013;381:1747–55.

    Article  PubMed  Google Scholar 

  7. Bateman BT. What’s new in obstetric anesthesia: a focus on maternal morbidity and mortality. Int J Obstet Anesth. 2019;37:68–72.

    Article  PubMed  Google Scholar 

  8. Aoyama K, Pinto R, Ray JG, Hill AD, Scales DC, Lapinsky SE, Hladunewich MA, Seaward GR, Fowler RA. Association of maternal age with severe maternal morbidity and mortality in Canada. JAMA Netw Open. 2019;2:e199875.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Gooch RA, Kahn JM. ICU bed supply, utilization, and health care spending: an example of demand elasticity. JAMA. 2014;311:567–8.

    Article  CAS  PubMed  Google Scholar 

  10. Rubenfeld GD, Rhodes A. How many intensive care beds are enough? Intensive Care Med. 2014;40:451–2.

    Article  PubMed  Google Scholar 

  11. Tu JV, Bowen J, Chiu M, Ko DT, Austin PC, He Y, Hopkins R, Tarride J-E, Blackhouse G, Lazzam C, Cohen EA, Goeree R. Effectiveness and safety of drug-eluting stents in Ontario. N Engl J Med. 2007;357:1393–402.

    Article  CAS  PubMed  Google Scholar 

  12. Nathens AB, Jurkovich GJ, Maier RV, Grossman DC, MacKenzie EJ, Moore M, Rivara FP. Relationship between trauma center volume and outcomes. JAMA. 2001;285:1164–71.

    Article  CAS  PubMed  Google Scholar 

  13. Kahn JM, Goss CH, Heagerty PJ, Kramer AA, O’Brien CR, Rubenfeld GD. Hospital volume and the outcomes of mechanical ventilation. N Engl J Med. 2006;355:41–50.

    Article  CAS  PubMed  Google Scholar 

  14. Yarnell CJ, Fu L, Manuel D, Tanuseputro P, Stukel T, Pinto R, Scales DC, Laupacis A, Fowler RA. Association between immigrant status and end-of-life care in Ontario, Canada. JAMA. 2017;318:1479–88.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Canadian Institute for Health Information. Standards and data submission data quality documentation, Discharge Abstract Database—multi-year information; 2012. p. 1–22.

    Google Scholar 

  16. Canadian Institute for Health Information. Data quality documentation, Discharge Abstract Database - current-year information, 2014–2015; 2015. p. 1–19.

    Google Scholar 

  17. Goel V, Williams J, Anderson G, Blackstien-Hirsch P, Fooks C, Naylor C. Patterns of health care in Ontario; 1996.

    Google Scholar 

  18. Scales DC, Thiruchelvam D, Kiss A, Redelmeier DA. The effect of tracheostomy timing during critical illness on long-term survival. Crit Care Med. 2008;36:2547–57.

    Article  PubMed  Google Scholar 

  19. Juurlink D, Preyra C, Croxford R, Chong A, Austin P, Tu J, Laupacis A. Canadian Institute for Health Information Discharge Abstract Database: a validation study. Toronto Inst Clin Eval Sci. 2006:1–69.

  20. Urquia ML, Stukel TA, Fung K, Glazier RH, Ray JG. Estimating gestational age at birth: a population-based derivation-validation study. Chronic Dis Inj Can. 2011;31:103–8.

    CAS  PubMed  Google Scholar 

  21. Garland A, Yogendran M, Olafson K, Scales DC, McGowan KL, Fransoo R. The accuracy of administrative data for identifying the presence and timing of admission to intensive care units in a Canadian province. Med Care. 2012;50:1–6.

    Article  Google Scholar 

  22. Joseph KS, Liu S, Rouleau J, Kirby RS, Kramer MS, Sauve R, Fraser WD, Young DC, Liston RM. Severe maternal morbidity in Canada, 2003 to 2007: surveillance using routine hospitalization data and ICD-10CA codes. J Obstet Gynaecol Canada. 2010(32):837–46.

    Article  CAS  PubMed  Google Scholar 

  23. Filippi V, Chou D, Barreix M, Say L. A new conceptual framework for maternal morbidity. Int J Gynecol Obstet. 2018;141:4–9.

    Article  Google Scholar 

  24. Aoyama K, D’Souza R, Inada E, Lapinsky SE, Fowler RA. Measurement properties of comorbidity indices in maternal health research: a systematic review. BMC Pregnancy Childbirth. 2017;17:372.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Metcalfe A, Lix LM, Johnson JA, Currie G, Lyon AW, Bernier F, Tough SC. Validation of an obstetric comorbidity index in an external population. BJOG An Int J Obstet Gynaecol. 2015;122:1748–55.

    Article  CAS  Google Scholar 

  26. Yang D, Dalton J. A unified approach to measuring the effect size between two groups using SAS®. SAS Glob Forum. 2012;6:1–6.

  27. Austin PC, Merlo J. Intermediate and advanced topics in multilevel logistic regression analysis. Stat Med. 2017;36:3257–77.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Austin PC, Goel V, Van Walraven C. An introduction to multilevel regression models. Can J Public Heal. 2001;92:150–4.

    Article  CAS  Google Scholar 

  29. Larsen K, Merlo J. Appropriate assessment of neighborhood effects on individual health: integrating random and fixed effects in multilevel logistic regression. Am J Epidemiol. 2005;161:81–8.

    Article  PubMed  Google Scholar 

  30. Merlo J, Chaix B, Ohlsson H, Beckman A, Johnell K, Hjerpe P, Råstam L, Larsen K. A brief conceptual tutorial of multilevel analysis in social epidemiology: using measures of clustering in multilevel logistic regression to investigate contextual phenomena. J Epidemiol Community Health. 2006;60:290–7.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Vesin A, Azoulay E, Ruckly S, Vignoud L, Rusinovà K, Benoit D, Soares M, Azeivedo-Maia P, Abroug F, Benbenishty J, Timsit JF. Reporting and handling missing values in clinical studies in intensive care units. Intensive Care Med. 2013;39:1396–404.

    Article  PubMed  Google Scholar 

  32. Haukoos JS, Newgard CD. Advanced statistics: missing data in clinical research-part 1: an introduction and conceptual framework. Acad Emerg Med. 2007;14:662–8.

    PubMed  Google Scholar 

  33. Newgard CD, Haukoos JS. Advanced statistics: missing data in clinical research-part 2: multiple imputation. Acad Emerg Med. 2007;14:669–78.

    PubMed  Google Scholar 

  34. Sterne JAC, White IR, Carlin JB, Spratt M, Royston P, Kenward MG, Wood AM, Carpenter JR. Multiple imputation for missing data in epidemiological and clinical research: potential and pitfalls. BMJ. 2009;338:b2393.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Lee KJ, Carlin JB. Multiple imputation for missing data: fully conditional specification versus multivariate normal imputation. Am J Epidemiol. 2010;171:624–32.

    Article  PubMed  Google Scholar 

  36. Graham JW, Olchowski AE, Gilreath TD. How many imputations are really needed? Some practical clarifications of multiple imputation theory. Prev Sci. 2007;8:206–13.

    Article  PubMed  Google Scholar 

  37. van Buuren S. Multiple imputation of discrete and continuous data by fully conditional specification. Stat Methods Med Res. 2007;16:219–42.

    Article  PubMed  Google Scholar 

  38. Wijeysundera DN, Austin PC, Beattie WS, Hux JE, Laupacis A. Variation in the practice of preoperative medical consultation for major elective noncardiac surgery. Anesthesiology. 2012;116:25–34.

    Article  PubMed  Google Scholar 

  39. Harrison D, Penny J, Yentis S, Fayek S, Brady A. Case mix, outcome and activity for obstetric admissions to adult, general critical care units: a secondary analysis of the ICNARC Case Mix Programme Database. Crit Care. 2005;9(Suppl 3):S25–37.

    Article  PubMed Central  Google Scholar 

  40. Madan I, Puri I, Jain NJ, Grotegut C, Nelson D, Dandolu V. Characteristics of obstetric intensive care unit admissions in New Jersey. J Matern Neonatal Med. 2009;22:785–90.

    Article  Google Scholar 

  41. Chantry AA, Deneux-Tharaux C, Bonnet MP, Bouvier-Colle MH. Pregnancy-related ICU admissions in France: trends in rate and severity, 2006-2009. Crit Care Med. 2015;43:78–86.

    Article  PubMed  Google Scholar 

  42. van Vliet EOG, Askie LA, Mol BWJ, Oudijk MA. Antiplatelet agents and the prevention of spontaneous preterm birth: a systematic review and meta-analysis. Obstet Gynecol. 2017;129:327–36.

    Article  PubMed  CAS  Google Scholar 

  43. Machado FR, Ferreira EM, Schippers P, de Paula IC, Saes LSV, de Oliveira FI, Tuma P, Nogueira Filho W, Piza F, Guare S, Mangini C, Guth GZ, Azevedo LCP, FGR F, JLG d A, Mansur NS, Salomão R. Implementation of sepsis bundles in public hospitals in Brazil: a prospective study with heterogeneous results. Crit Care. 2017;21(1):268.

  44. Wilcox ME, Adhikari NKJ. The effect of telemedicine in critically ill patients: systematic review and meta-analysis. Crit Care. 2012;16(4):R127.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Kovacevic P, Dragic S, Kovacevic T, Momcicevic D, Festic E, Kashyap R, Niven AS, Dong Y, Gajic O. Impact of weekly case-based tele-education on quality of care in a limited resource medical intensive care unit. Crit Care. 2019;23:1–7.

    Article  Google Scholar 

  46. Wiegersma JS, Droogh JM, Zijlstra JG, Fokkema J, Ligtenberg JJM. Quality of interhospital transport of the critically ill: impact of a mobile intensive care unit with a specialized retrieval team. Crit Care. 2011;15:1–8.

    Article  Google Scholar 

  47. Bateman BT, Mhyre JM, Hernandez-Diaz S, Huybrechts KF, Fischer MA, Creanga AA, Callaghan WM, Gagne JJ. Development of a comorbidity index for use in obstetric patients. Obstet Gynecol. 2013;122:957–65.

    Article  PubMed  Google Scholar 

  48. Aoyama K, D’Souza R, Pinto R, Ray JG, Hill A, Scales DC, Lapinsky SE, Seaward GR, Hladunewich M, Shah PS, Fowler RA. Risk prediction models for maternal mortality: a systematic review and meta-analysis. PLoS One. 2018;13:e0208563.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Canadian Institute for Health Information: Canadian Coding Standards Evolution Chronicle. Ottawa; 2018

  50. Fowler RA, Abdelmalik P, Wood G, Foster D, Gibney N, Bandrauk N, Turgeon AF, Lamontagne F, Kumar A, Zarychanski R, Green R, Bagshaw SM, Stelfox HT, Foster R, Dodek P, Shaw S, Granton J, Lawless B, Hill A, Rose L, Adhikari NK, Scales DC, Cook DJ, Marshall JC, Martin C, Jouvet P, Canadian Critical Care Trials Group, Canadian ICU Capacity Group. Critical care capacity in Canada: results of a national cross-sectional study. Crit Care. 2015;19:133.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

KA received Fellowship funding from the Canadian Institutes of Health Research during 2013 and 2017, to pursue the PhD stream of Clinical Epidemiology & Health Care Research at Institute of Health Policy, Management and Evaluation, University of Toronto.

Funding

KA, JGR, DCS, SEL, GRS, MH, and RAF have received a Canadian Institutes of Health Research Transitional Open Operating Grant for this project (#342397). Funders had no role in the analysis or decision to submit this manuscript for peer review.

Author information

Authors and Affiliations

Authors

Contributions

KA conceived this paper. KA, RP, JGR, AH, DCS, and RAF developed the study protocol. KA, RP, and RAF performed all the analyses. RP, JGR, AH, DCS, SEL, GRS, MH, and RAF provided input on the interpretation. KA and RAF wrote the initial draft of the manuscript. RP, JGR, AH, DCS, SEL, GRS, and MH helped to draft the final version, which was approved by all authors.

Corresponding author

Correspondence to Kazuyoshi Aoyama.

Ethics declarations

Ethics approval and consent to participate

This study was conducted under the data security and privacy policy of CIHI and approved by the institutional review board at Mount Sinai Hospital in Toronto, Canada (10-0305-C) and University of Toronto (#34468).

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests. Full disclosure of interests available to view online as supporting information.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Table S1. List of all variables used in this study.

Additional file 2:

Table S2. The number of severe maternal morbidity events with ICU admission.

Additional file 3:

Table S3. The number of severe maternal morbidity events with maternal death.

Additional file 4:

Table S4. Characteristics of categorized 5 hospital groups according to hospital volume.

Additional file 5:

Table S5. Characteristics of the cohort by province/territory.

Additional file 6:

Table S6. Severe maternal morbidity diagnoses that could abruptly or slowly progress from clinical point of view.

Additional file 7:

Table S7. Estimated regression coefficients and variance components for the multi-level mixed logistic regression models for the outcome of ICU admission [Outcome=ICU admission, main predictors= quintile of hospitals according to the number of pregnancy admission at each hospital].

Additional file 8:

Table S8. Variability of ICU admission stratified by different type of admission.

Additional file 9:

Table S9. Sensitivity analysis: Estimated regression coefficients and variance components for the multi-level mixed logistic regression models with the first modification of main predictors for the outcome of ICU admission [Outcome=ICU admission, main predictors=the number of pregnancy admission at each hospital (continuous variable)].

Additional file 10:

Table S10. Sensitivity analysis: Estimated regression coefficients and variance components for the multi-level mixed logistic regression models with the second modification of main predictors for the outcome of highest acuity ICU admission [Outcome= highest acuity ICU admission, main predictors= quintile of hospitals according to the number of pregnancy admission at each hospital].

Additional file 11:

Table S11. Sensitivity analysis: Estimated regression coefficients and variance components for the multi-level mixed logistic regression models for the outcome of ICU admission within multiple imputation datasets [Outcome=ICU admission, main predictors=quintile of hospitals according to the number of pregnancy admission at each hospital].

Additional file 12:

Table S12. Search strategy in MEDLINE, EMBASE, Cochrane.

Additional file 13:

Table S13. General characteristics of maternal death according to whether admitted to ICU or not.

Additional file 14:

Table S14. Estimated odds ratios for multilevel logistic regression model for Intensive Care Unit admission with categorized Maternal Comorbidity Index.

Additional file 15:

Figure S1. Provincial ICU Beds per 100,000 population and Provincial Intensive Care Unit Admissions per 1000 Live Births.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Aoyama, K., Pinto, R., Ray, J.G. et al. Variability in intensive care unit admission among pregnant and postpartum women in Canada: a nationwide population-based observational study. Crit Care 23, 381 (2019). https://doi.org/10.1186/s13054-019-2660-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13054-019-2660-x

Keywords