Immunosuppression and cancer risk in kidney transplant recipients: A retrospective cohort study.

We assessed whether contemporary immunosuppression agents were associated with cancer among kidney transplant recipients (KTR), and if this association varied by age and sex. We studied a retrospective province-wide cohort of primary KTR (1997-2016). Employing multivariable Cox models, we estimated associations of cumulative doses of prednisone, mycophenolate and tacrolimus administered over the past 10 years, lagged by 2 years, with the incidence of primary malignant neoplasms (PMN). We assessed interactions with age and sex. To assess the impact of exposure recency, we used weighted cumulative exposure (WCE) modeling. Among 1064 KTR, 108 (10.2%) developed PMN over median follow-up of 73 months (interquartile range: 32-120). Adjusted hazard ratios (aHRs) and 95% confidence intervals (CIs) of 0.96 (0.64-1.43), 1.34 (0.96-1.86), and 1.06 (0.88-1.29) were estimated for cumulative daily doses of prednisone (5 mg), mycophenolate (1000 mg), and tacrolimus (2 mg) administered continuously over the past 10 years, respectively. PMN risk associated with cumulative tacrolimus exposure was modified by age (interaction p = .035) and was more pronounced in 15-year and 30-year-old KTR (aHRs of 1.57 [1.08-2.28] and 1.31 [1.03-1.66], respectively) in comparison to older KTR. PMN risk increase associated with higher cumulative mycophenolate dose was more pronounced in females (aHR = 1.86 [1.15-3.00]) than in males (aHR = 1.16 [0.74-1.81]; interaction p = .131). WCE analyses suggested increased PMN risk the higher the mycophenolate doses taken 5-10 years ago. A trend toward increased PMN risk with long-term mycophenolate exposure, particularly in females, and more pronounced risk with long-term tacrolimus exposure in younger KTR, identify opportunities for tailored immunosuppression to mitigate cancer risk.

Data-Driven Simulations to Assess the Impact of Study Imperfections in Time-to-Event Analyses.

Quantitative bias analysis (QBA) permits assessment of the expected impact of various imperfections of the available data on the results and conclusions of a particular real-world study. This article extends QBA methodology to multivariable time-to-event analyses with right-censored endpoints, possibly including time-varying exposures or covariates. The proposed approach employs data-driven simulations, which preserve important features of the data at hand while offering flexibility in controlling the parameters and assumptions that may affect the results. First, the steps required to perform data-driven simulations are described, and then two examples of real-world time-to-event analyses illustrate their implementation and the insights they may offer. The first example focuses on the omission of an important time-invariant predictor of the outcome in a prognostic study of cancer mortality, and permits separating the expected impact of confounding bias from non-collapsibility. The second example assesses how imprecise timing of an interval-censored event - ascertained only at sparse times of clinic visits - affects its estimated association with a time-varying drug exposure. The simulation results also provide a basis for comparing the performance of two alternative strategies for imputing the unknown event times in this setting. The R scripts that permit the reproduction of our examples are provided.

Predicted 25-hydroxyvitamin D over the adult life and the risk of ovarian cancer.

The evidence from previous studies of serum 25-hydroxyvitamin D [25(OH)D] and ovarian cancer risk are not conclusive. However, 25(OH)D was generally only measured in late adulthood, which may not capture the etiologically relevant exposure periods. We investigated predicted 25(OH)D over the adult lifetime in relation to ovarian cancer risk in a population-based case-control study conducted from 2011 to 2016 in Montreal, Canada (490 cases, 896 controls). Predicted 25(OH)D was computed using previously validated regression models. Unconditional multivariable logistic regression models were used to estimate adjusted odds ratios (aOR) and 95% confidence intervals (CI) for average predicted 25(OH)D over the adult life and risk. In addition, the relative importance of different periods of past 25(OH)D exposure was explored using a weighted cumulative exposure (WCE) model. For each 20 nmol/L increase in average predicted 25(OH)D over the adult life, the aOR (95% CI) was 0.73 (0.55-0.96). In WCE analyses, the inverse association was strongest for exposures 5 to 20 years and 35 to 55 years prior to diagnosis, with aORs (95% CIs) of 0.82 (0.69-0.94) and 0.79 (0.66-1.02), respectively, for each 20 nmol/L increase in predicted 25(OH)D. These results support an inverse association between 25(OH)D in adulthood and ovarian cancer risk.

Recipient age influences survival after liver transplant: Results of the French national cohort 2007-2017.

BACKGROUND: In recent years, age at liver transplantation (LT) has markedly increased. In the context of organ shortage, we investigated the impact of recipient age on post-transplantation mortality. METHODS: All adult patients who received a first LT between 2007 and 2017 were included in this cross-sectional study. Recipients' characteristics at the time of listing, donor and surgery data, post-operative complications and follow-up of vital status were retrieved from the national transplantation database. The impact of age on 5-year overall mortality post-LT was estimated using a flexible multivariable parametric model which was also used to estimate the association between age and 10-year net survival, accounting for expected age- and sex-related mortality. RESULTS: Among the 7610 patients, 21.4% were aged 60-65 years, and 15.7% over 65. With increasing age, comorbidities increased but severity of liver disease decreased. Older recipient age was associated with decreased observed survival at 5 years after LT (p < .001), with a significant effect particularly during the first 2 years. The linear increase in the risk of death associated with age does not allow any definition of an age's threshold for LT (p = .832). Other covariates associated with an increased risk of 5-year death were dialysis and mechanical ventilation at transplant, transfusion during LT, hepatocellular carcinoma and donor age. Ten-year flexible net survival analysis confirmed these results. CONCLUSION: Although there was a selection process for older recipients, increasing age at LT was associated with an increased risk of death, particularly in the first years after LT.

Childhood body fatness and the risk of epithelial ovarian cancer: A population-based case-control study in Montreal, Canada.

OBJECTIVE: To assess the association between childhood body fatness and epithelial ovarian cancer (EOC), and whether this association differs by type of EOC. METHODS: Using data from a population-based case-control study (497 cases and 902 controls) in Montreal, Canada conducted 2011-2016, we examined the association between childhood body fatness and EOC, overall and separately for invasive vs. borderline EOCs. A figure rating scale was used to measure body fatness at ages 5 and 10. Multivariable logistic regression was used to estimate adjusted odds ratios (aORs) and 95% confidence intervals (95% CI). Quantitative bias analyses were conducted to assess the impact of exposure misclassification and non-participation. RESULTS: The aOR (95% CI) of overall EOC for high vs. low body fatness was 1.07 (0.85-1.34) at age 5 and 1.28 (0.98-1.68) at age 10. The associations were stronger for invasive EOC, specifically the endometrioid histological type. For borderline cancers, the aORs were below the null value with wide confidence intervals. Bias analyses did not reveal a strong influence of non-participation. Non-differential exposure misclassification may have biased aORs towards the null for invasive cancers but did not appear to have an appreciable influence on the aORs for borderline cancers. CONCLUSIONS: Childhood body fatness may be a risk factor for invasive EOC in later adult life. Our study highlights the potential importance of examining early life factors for a comprehensive understanding of EOC development.

Comparison of Different Modeling Approaches for Prescription Opioid Use and Its Association With Adverse Events.

Previous research linking opioid prescribing to adverse drug events has failed to properly account for the time-varying nature of opioid exposure. This study aimed to explore how the risk of opioid-related emergency department visits, readmissions, or deaths (composite outcome) varies with opioid dose and duration, comparing different novel modeling techniques. A prospective cohort of 1,511 hospitalized patients discharged from 2 McGill-affiliated hospitals in Montreal, 2014-2016, was followed from the first postdischarge opioid dispensation until 1 year after discharge. Marginal structural Cox proportional hazards models and their flexible extensions were used to explore the association between time-varying opioid use and the composite outcome. Weighted cumulative exposure models assessed cumulative effects of past use and explored how its impact depends on the recency of exposure. The patient mean age was 69.6 (standard deviation = 14.9) years; 57.7% were male. In marginal structural model analyses, current opioid use was associated with a 71% increase in the hazard of opioid-related adverse events (adjusted hazard ratio = 1.71, 95% confidence interval: 1.21, 2.43). The weighted cumulative exposure results suggested that the risk cumulates over the previous 50 days of opioid consumption. Flexible modeling techniques helped assess how the risk of opioid-related adverse events may be associated with time-varying opioid exposures while accounting for nonlinear relationships and the recency of past use.

Investigating the Shape and Strength of the Relationship Between Maternal Weight Gain and Gestational Age at Delivery in Twin and Singleton Pregnancies.

Both inadequate and excessive maternal weight gain are correlated with preterm delivery in singleton pregnancies, yet this relationship has not been adequately studied in twins. We investigated the relationship between time-varying maternal weight gain and gestational age at delivery in twin pregnancies and compared it with that in singletons delivered in the same study population. We used serial weight measurements abstracted from charts for twin and singleton pregnancies delivered during 1998-2013 in Pittsburgh, Pennsylvania. Our exposure was time-varying weight gain z score, calculated using gestational age-standardized and prepregnancy body mass index-stratified twin- and singleton-specific charts, and our outcome was gestational age at delivery. Our analyses used a flexible extension of the Cox proportional hazards model that allowed for nonlinear and time-dependent effects. We found a U-shaped relationship between weight gain z score and gestational age at delivery among twin pregnancies (lowest hazard of delivery observed at z score = 1.2), which we attributed to increased hazard of early preterm spontaneous delivery among pregnancies with low weight gain and increased hazard of late preterm delivery without labor among pregnancies with high weight gain. Our findings may be useful for updating provisional guidelines for maternal weight gain in twin pregnancies.

Alcohol intake and the risk of epithelial ovarian cancer.

PURPOSE: To investigate the association between alcohol intake over the lifetime and the risk of overall, borderline, and invasive ovarian cancer. METHODS: In a population-based case-control study of 495 cases and 902 controls, conducted in Montreal, Canada, average alcohol intake over the lifetime and during specific age periods were computed from a detailed assessment of the intake of beer, red wine, white wine and spirits. Multivariable logistic regression was used to estimate odds ratios (ORs) with 95% confidence intervals (CIs) for the association between alcohol intake and ovarian cancer risk. RESULTS: For each one drink/week increment in average alcohol intake over the lifetime, the adjusted OR (95% CI) was 1.06 (1.01-1.10) for ovarian cancer overall, 1.13 (1.06-1.20) for borderline ovarian cancers and 1.02 (0.97-1.08) for invasive ovarian cancers. This pattern of association was similarly observed for alcohol intake in early (15- < 25 years), mid (25- < 40 years) and late adulthood (≥ 40 years), as well as for the intake of specific alcohol beverages over the lifetime. CONCLUSIONS: Our results support the hypothesis that a higher alcohol intake modestly increases the risk of overall ovarian cancer, and more specifically, borderline tumours.

Importance of accounting for timing of time-varying exposures in association studies: Hydrochlorothiazide and non-melanoma skin cancer.

PURPOSE: Hydrochlorothiazide (HCTZ), a widely prescribed antihypertensive drug with photosensitising properties, has been linked with non-melanoma skin cancer (NMSC) risk. However, previous analyses did not fully explore if and how the impact of past HCTZ exposures accumulates with prolonged use and/or depends on time elapsed since exposures. Therefore, we used different models to more comprehensively assess how NMSC risk vary with HCTZ exposure, and explore how the results may depend on modeling strategies. METHODS: We used different parametric models with alternative time-varying exposure metrics, and the flexible weighted cumulative exposure model (WCE) to estimate associations between HCTZ exposures and NMSC risk in a population-based cohort of HCTZ users over 65 years old, in the province of Ontario, Canada. RESULTS: Among 3844 HCTZ users, 273 developed NMSC during up to 8 years of follow-up. In parametric models, based on all exposures, increased duration of past HCTZ use was associated with an increase of NMSC risk but cumulative dose showed no systematic association. Yet, WCE results suggested that only exposures taken 2.5-4 years in the past were associated with the current NMSC hazard. This finding led us to re-define the parametric models, which also confirmed that any HCTZ dose taken outside this time-window were not systematically associated with NMSC incidence. CONCLUSIONS: Our analyses illustrate how flexible modeling may yield new insights into complex temporal relationships between a time-varying drug exposure and risks of adverse events. Duration and recency of antihypertensive agents exposures must be taken into account in evaluating risk and benefits.

Comparison of cohort and nested case-control designs for estimating the effect of time-varying drug exposure on the risk of adverse event in the presence of ties.

Cohort and nested case-control (NCC) designs are frequently used in pharmacoepidemiology to assess the associations of drug exposure that can vary over time with the risk of an adverse event. Although it is typically expected that estimates from NCC analyses are similar to those from the full cohort analysis, with moderate loss of precision, only few studies have actually compared their respective performance for estimating the effects of time-varying exposures (TVE). We used simulations to compare the properties of the resulting estimators of these designs for both time-invariant exposure and TVE. We varied exposure prevalence, proportion of subjects experiencing the event, hazard ratio, and control-to-case ratio and considered matching on confounders. Using both designs, we also estimated the real-world associations of time-invariant ever use of menopausal hormone therapy (MHT) at baseline and updated, time-varying MHT use with breast cancer incidence. In all simulated scenarios, the cohort-based estimates had small relative bias and greater precision than the NCC design. NCC estimates displayed bias to the null that decreased with a greater number of controls per case. This bias markedly increased with higher proportion of events. Bias was seen with Breslow's and Efron's approximations for handling tied event times but was greatly reduced with the exact method or when NCC analyses were matched on confounders. When analyzing the MHT-breast cancer association, differences between the two designs were consistent with simulated data. Once ties were taken correctly into account, NCC estimates were very similar to those of the full cohort analysis.

Flares after hydroxychloroquine reduction or discontinuation: results from the Systemic Lupus International Collaborating Clinics (SLICC) inception cohort.

OBJECTIVES: To evaluate systemic lupus erythematosus (SLE) flares following hydroxychloroquine (HCQ) reduction or discontinuation versus HCQ maintenance. METHODS: We analysed prospective data from the Systemic Lupus International Collaborating Clinics (SLICC) cohort, enrolled from 33 sites within 15 months of SLE diagnosis and followed annually (1999-2019). We evaluated person-time contributed while on the initial HCQ dose ('maintenance'), comparing this with person-time contributed after a first dose reduction, and after a first HCQ discontinuation. We estimated time to first flare, defined as either subsequent need for therapy augmentation, increase of ≥4 points in the SLE Disease Activity Index-2000, or hospitalisation for SLE. We estimated adjusted HRs (aHRs) with 95% CIs associated with reducing/discontinuing HCQ (vs maintenance). We also conducted separate multivariable hazard regressions in each HCQ subcohort to identify factors associated with flare. RESULTS: We studied 1460 (90% female) patients initiating HCQ. aHRs for first SLE flare were 1.20 (95% CI 1.04 to 1.38) and 1.56 (95% CI 1.31 to 1.86) for the HCQ reduction and discontinuation groups, respectively, versus HCQ maintenance. Patients with low educational level were at particular risk of flaring after HCQ discontinuation (aHR 1.43, 95% CI 1.09 to 1.87). Prednisone use at time-zero was associated with over 1.5-fold increase in flare risk in all HCQ subcohorts. CONCLUSIONS: SLE flare risk was higher after HCQ taper/discontinuation versus HCQ maintenance. Decisions to maintain, reduce or stop HCQ may affect specific subgroups differently, including those on prednisone and/or with low education. Further study of special groups (eg, seniors) may be helpful.

Adapting SIMEX to correct for bias due to interval-censored outcomes in survival analysis with time-varying exposure.

Many clinical and epidemiological applications of survival analysis focus on interval-censored events that can be ascertained only at discrete times of clinic visits. This implies that the values of time-varying covariates are not correctly aligned with the true, unknown event times, inducing a bias in the estimated associations. To address this issue, we adapted the simulation-extrapolation (SIMEX) methodology, based on assessing how the estimates change with the artificially increased time between clinic visits. We propose diagnostics to choose the extrapolating function. In simulations, the SIMEX-corrected estimates reduced considerably the bias to the null and generally yielded a better bias/variance trade-off than conventional estimates. In a real-life pharmacoepidemiological application, the proposed method increased by 27% the excess hazard of the estimated association between a time-varying exposure, representing the 2-year cumulative duration of past use of a hypertensive medication, and the hazard of nonmelanoma skin cancer (interval-censored events). These simulation-based and real-life results suggest that the proposed SIMEX-based correction may help improve the accuracy of estimated associations between time-varying exposures and the hazard of interval-censored events in large cohort studies where the events are recorded only at relatively sparse times of clinic visits/assessments. However, these advantages may be less certain for smaller studies and/or weak associations.

Ten-year risk of cerebrovascular accidents in incident rheumatoid arthritis: a population-based study of trends over time.

OBJECTIVE: To evaluate secular trends in 10-year risk of incident cerebrovascular accidents (CVA), in incident RA relative to the general population. METHODS: We conducted a retrospective study of a population-based incident cohort with RA onset from 1997 to 2004 in British Columbia, Canada, with matched general population controls (2:1), using administrative health data. RA and general population cohorts were divided according to year of RA onset, defined according to the first RA visit of the case definition. Incident CVA was defined as the first CVA occurring within 10 years from the first RA visit. Secular trend was assessed using delayed-entry Cox models with a two-way interaction term between the year of RA onset and indicator of RA vs general population. Linear, quadratic and spline functions of year of RA onset were compared with assess non-linear effects. The model with the lowest Akaike Information Criterion was selected. RESULTS: Overall, 23 545 RA and 47 090 general population experienced 658 and 1220 incident CVAs, respectively. A spline Cox model with a knot at year of onset 1999 was selected. A significant decline in risk of CVA was observed in individuals with RA onset after 1999 [0.90 (0.86, 0.95); P = 0.0001]. The change in CVA risk over time differed significantly in RA with onset from 1999 onwards compared with the general population (P-value of interaction term = 0.03), but not before 1999 (P = 0.06). CONCLUSION: Our findings suggest that people with RA onset from 1999 onwards, had a significantly greater decline in 10-year risk of CVA compared with the general population.

Spline-based accelerated failure time model.

The accelerated failure time (AFT) model has been suggested as an alternative to the Cox proportional hazards model. However, a parametric AFT model requires the specification of an appropriate distribution for the event time, which is often difficult to identify in real-life studies and may limit applications. A semiparametric AFT model was developed by Komárek et al based on smoothed error distribution that does not require such specification. In this article, we develop a spline-based AFT model that also does not require specification of the parametric family of event time distribution. The baseline hazard function is modeled by regression B-splines, allowing for the estimation of a variety of smooth and flexible shapes. In comprehensive simulations, we validate the performance of our approach and compare with the results from parametric AFT models and the approach of Komárek. Both the proposed spline-based AFT model and the approach of Komárek provided unbiased estimates of covariate effects and survival curves for a variety of scenarios in which the event time followed different distributions, including both simple and complex cases. Spline-based estimates of the baseline hazard showed also a satisfactory numerical stability. As expected, the baseline hazard and survival probabilities estimated by the misspecified parametric AFT models deviated from the truth. We illustrated the application of the proposed model in a study of colon cancer.

Analysis of time-to-event for observational studies: Guidance to the use of intensity models.

This paper provides guidance for researchers with some mathematical background on the conduct of time-to-event analysis in observational studies based on intensity (hazard) models. Discussions of basic concepts like time axis, event definition and censoring are given. Hazard models are introduced, with special emphasis on the Cox proportional hazards regression model. We provide check lists that may be useful both when fitting the model and assessing its goodness of fit and when interpreting the results. Special attention is paid to how to avoid problems with immortal time bias by introducing time-dependent covariates. We discuss prediction based on hazard models and difficulties when attempting to draw proper causal conclusions from such models. Finally, we present a series of examples where the methods and check lists are exemplified. Computational details and implementation using the freely available R software are documented in Supplementary Material. The paper was prepared as part of the STRATOS initiative.

Impact of time-varying cumulative bevacizumab exposures on survival: re-analysis of data from randomized clinical trial in patients with metastatic colo-rectal cancer.

BACKGROUND: As cancer treatment, biotherapies can be as effective as chemotherapy while reducing the risk of secondary effects, so that they can be taken over longer periods than conventional chemotherapy. Thus, some trials aimed at assessing the benefit of maintaining biotherapies during chemotherapy-free intervals (CFI). For example, the recent PRODIGE9 trial assessed the effect of maintaining bevacizumab during CFI in metastatic colorectal cancer (mCRC) patients. However, its analysis was hindered by a small difference of exposure to the treatment between the randomized groups and by a large proportion of early drop outs, leading to a potentially unbalanced distribution of confounding factors among the trial completers. To address these limitations, we re-analyzed the PRODIGE9 data to assess the effects of different exposure metrics on all-cause mortality of patients with mCRC using methods originally developed for observational studies. METHODS: To account for the actual patterns of drug use by individual patients and for possible cumulative effects, we used five alternative time-varying exposure metrics: (i) cumulative dose, (ii) quantiles of the cumulative dose, (iii) standardized cumulative dose, (iv) Theoretical Blood Concentration (TBC), and (v) Weighted Cumulative Exposure (WCE). The last two metrics account for the timing of drug use. Treatment effects were estimated using adjusted Hazard Ratio from multivariable Cox proportional hazards models. RESULTS: After excluding 112 patients who died during the induction period, we analyzed data on 382 patients, among whom 320 (83.8%) died. All time-varying exposures improved substantially the model's fit to data, relative to using only the time-invariant randomization group. All exposures indicated a protective effect for higher cumulative bevacizumab doses. The best-fitting WCE and TBC models accounted for both the cumulative effects and the different impact of doses taken at different times. CONCLUSIONS: All time-varying analyses, regardless of the exposure metric used, consistently suggested protective effects of higher cumulative bevacizumab doses. However, the results may partly reflect the presence of a confusion bias. Complementing the main ITT analysis of maintenance trials with an analysis of potential cumulative effects of treatment actually taken can provide new insights, but the results must be interpreted with caution because they do not benefit from the randomization. TRIAL REGISTRATION: clinicaltrials.gov, NCT00952029 . Registered 8 August 2009.

Pelvic floor morphometrical and functional changes immediately after pelvic floor muscle training and at 1-year follow-up, in older incontinent women.

AIMS: To compare the effects of group-based and individual physiotherapy for stress or mixed urinary incontinence (UI) on pelvic floor morphometry, pelvic floor muscle (PFM) function, and related self-efficacy, immediately after treatment and at 1 year. METHODS: This is a planned secondary analysis of the group rehabilitation or individual physiotherapy study, an assessor-blinded, randomized, noninferiority trial. Eligible participants included 362 community-dwelling older women with symptoms of stress/mixed UI. After learning how to contract PFMs, participants completed 12 weeks of PFM training, either individually (one-on-one) or as part of a group (eight women). Pelvic floor transperineal ultrasound volumes (morphometry), PFM intravaginal dynamometric data (function), and self-efficacy in performing PFM exercises were acquired at baseline, posttreatment, and at 1 year. RESULTS: Groups were comparable at all time points. Immediately posttreatment, both groups demonstrated significant changes in pelvic floor morphometry during coughs, and in PFM function during contractions and coughs. Participants also reported improved self-efficacy in performing PFM exercises. Results were sustained at 1 year. When participants coughed, pelvic floor structures were better supported (reflected by less caudal movement of the puborectalis sling and a smaller opening of the levator hiatus) in a pattern consistent with the "knack" strategy. Furthermore, both interventions resulted in stronger, faster, more coordinated, and more endurant PFMs. CONCLUSION: In older women with stress or mixed UI, both individual and group-based PFM training resulted in comparable improvements in overall PFM function, pelvic floor morphometry during coughs, and related self-efficacy in performing PFM exercises, which were sustained at 1 year.

Nurse staffing practices and adverse events in acute care hospitals: The research protocol of a multisite patient-level longitudinal study.

AIMS: We describe an innovative research protocol to: (a) examine patient-level longitudinal associations between nurse staffing practices and the risk of adverse events in acute care hospitals and; (b) determine possible thresholds for safe nurse staffing. DESIGN: A dynamic cohort of adult medical, surgical and intensive care unit patients admitted to 16 hospitals in Quebec (Canada) between January 2015-December 2019. METHODS: Patients in the cohort will be followed from admission until 30-day postdischarge to assess exposure to selected nurse staffing practices in relation to the subsequent occurrence of adverse events. Five staffing practices will be measured for each shift of an hospitalization episode, using electronic payroll data, with the following time-varying indicators: (a) nursing worked hours per patient; (b) skill mix; (c) overtime use; (d) education mix and; and (e) experience. Four high-impact adverse events, presumably associated with nurse staffing practices, will be measured from electronic health record data retrieved at the participating sites: (a) failure-to-rescue; (b) in-hospital falls; (c) hospital-acquired pneumonia and; and (d) venous thromboembolism. To examine the associations between the selected nurse staffing exposures and the risk of each adverse event, separate multivariable Cox proportional hazards frailty regression models will be fitted, while adjusting for patient, nursing unit and hospital characteristics, and for clustering. To assess for possible staffing thresholds, flexible non-linear spline functions will be fitted. Funding for the study began in October 2019 and research ethics/institutional approval was granted in February 2020. DISCUSSION: To our knowledge, this study is the first multisite patient-level longitudinal investigation of the associations between common nurse staffing practices and the risk of adverse events. It is hoped that our results will assist hospital managers in making the most effective use of the scarce nursing resources and in identifying staffing practices that minimize the occurrence of adverse events.

Lifetime recreational moderate-to-vigorous physical activity and ovarian cancer risk: A case-control study.

Results of epidemiologic studies of physical activity and ovarian cancer risk are inconsistent. Few have attempted to measure physical activity over the lifetime or in specific age windows, which may better capture etiologically relevant exposures. We examined participation in moderate-to-vigorous recreational physical activity (MVPA) in relation to ovarian cancer risk. In a population-based case-control study conducted in Montreal, Canada from 2011 to 2016 (485 cases and 887 controls), information was collected on lifetime participation in various recreational physical activities, which was used to estimate MVPA for each participant. MVPA was represented as average energy expenditure over the lifetime and in specific age-periods in units of metabolic equivalents (METs)-hours per week. Odds ratios (OR) and 95% confidence intervals (CI) for the relation between average MVPA and ovarian cancer risk were estimated using multivariable logistic regression models. Confounding was assessed using directed acyclic graphs combined with a change-in-estimate approach. The adjusted OR (95% CI) for each 28.5 MET-hr/week increment of lifetime recreational MVPA was 1.11 (0.99-1.24) for ovarian cancer overall. ORs for individual age-periods were weaker. When examined by menopausal status, the OR (95% CI) for lifetime MVPA was 1.21 (1.00-1.45) for those diagnosed before menopause and 1.04 (0.89-1.21) for those diagnosed postmenopausally. The suggestive positive associations were stronger for invasive ovarian cancers and more specifically for high-grade serous carcinomas. These results do not support a reduced ovarian cancer risk associated with MVPA.

Associations of 4 Nurse Staffing Practices With Hospital Mortality.

BACKGROUND: Cross-sectional studies of hospital-level administrative data have suggested that 4 nurse staffing practices-using adequate staffing levels, higher proportions of registered nurses (RNs) (skill mix), and more educated and experienced RNs-are each associated with reduced hospital mortality. To increase the validity of this evidence, patient-level longitudinal studies assessing the simultaneous associations of these staffing practices with mortality are required. METHODS: A dynamic cohort of 146,349 adult medical, surgical, and intensive care patients admitted to a Canadian University Health Center was followed for 7 years (2010-2017). We used a multivariable Cox proportional hazards model to estimate the associations between patients' time-varying cumulative exposure to measures of RN understaffing, skill mix, education, and experience, each relative to nursing unit and shift means, and the hazard of in-hospital mortality, while adjusting for patient and nursing unit characteristics, and modeling the current nursing unit of hospitalization as a random effect. RESULTS: Overall, 4854 in-hospital deaths occurred during 3,478,603 patient-shifts of follow-up (13.95 deaths/10,000 patient-shifts). In multivariable analyses, every 5% increase in the cumulative proportion of understaffed shifts was associated with a 1.0% increase in mortality (hazard ratio: 1.010; 95% confidence interval: 1.002-1.017; P=0.009). Moreover, every 5% increase in the cumulative proportion of worked hours by baccalaureate-prepared RNs was associated with a 2.0% reduction of mortality (hazard ratio: 0.980; 95% confidence interval: 0.965-0.995, P=0.008). RN experience and skill mix were not significantly associated with mortality. CONCLUSION: Reducing the frequency of understaffed shifts and increasing the proportion of baccalaureate-prepared RNs are associated with reduced hospital mortality.