Pioglitazone use and risk of bladder cancer: population based cohort studyBMJ 2016; 352 doi: https://doi.org/10.1136/bmj.i1541 (Published 30 March 2016) Cite this as: BMJ 2016;352:i1541
- Marco Tuccori, postdoctoral fellow1 2,
- Kristian B Filion, assistant professor of medicine1 2 3,
- Hui Yin, statistician1,
- Oriana H Yu, endocrinologist1 4,
- Robert W Platt, professor of biostatistics1 2 5 6,
- Laurent Azoulay, associate professor of oncology1 7
- 1Centre for Clinical Epidemiology, Lady Davis Institute, Jewish General Hospital, Montreal, Quebec, Canada
- 2Department of Epidemiology, Biostatistics, and Occupational Health, McGill University, Montreal
- 3Division of Clinical Epidemiology, McGill University
- 4Division of Endocrinology, Jewish General Hospital, Montreal
- 5Department of Pediatrics, McGill University, Montreal
- 6Research Institute of the McGill University Health Centre, Montreal
- 7Department of Oncology, McGill University
- Correspondence to: L Azoulay
- Accepted 29 February 2016
Objective To determine whether pioglitazone compared with other antidiabetic drugs is associated with an increased risk of bladder cancer in people with type 2 diabetes.
Design Population based cohort study.
Setting General practices contributing data to the United Kingdom Clinical Practice Research Datalink.
Participants A cohort of 145 806 patients newly treated with antidiabetic drugs between 1 January 2000 and 31 July 2013, with follow-up until 31 July 2014.
Main outcome measures The use of pioglitazone was treated as a time varying variable, with use lagged by one year for latency purposes. Cox proportional hazards models were used to estimate adjusted hazard ratios with 95% confidence intervals of incident bladder cancer associated with pioglitazone overall and by both cumulative duration of use and cumulative dose. Similar analyses were conducted for rosiglitazone, a thiazolidinedione not previously associated with an increased risk of bladder cancer.
Results The cohort generated 689 616 person years of follow-up, during which 622 patients were newly diagnosed as having bladder cancer (crude incidence 90.2 per 100 000 person years). Compared with other antidiabetic drugs, pioglitazone was associated with an increased risk of bladder cancer (121.0 v 88.9 per 100 000 person years; hazard ratio 1.63, 95% confidence interval 1.22 to 2.19). Conversely, rosiglitazone was not associated with an increased risk of bladder cancer (86.2 v 88.9 per 100 000 person years; 1.10, 0.83 to 1.47). Duration-response and dose-response relations were observed for pioglitazone but not for rosiglitazone.
Conclusion The results of this large population based study indicate that pioglitazone is associated with an increased risk of bladder cancer. The absence of an association with rosiglitazone suggests that the increased risk is drug specific and not a class effect.
Pioglitazone, an antidiabetic drug belonging to the thiazolidinedione class, has been shown to improve glycaemic levels in people with type 2 diabetes.1 However, in 2005 the PROactive randomised controlled trial unexpectedly showed an imbalance in the number of cases of bladder cancer with pioglitazone compared with placebo.2 In contrast, this imbalance was never observed in randomised controlled trials of rosiglitazone, the other approved drug belonging to the thiazolidinedione class.1 3
The findings of the PROactive trial were subsequently corroborated in some,4 5 6 7 8 9 10 but not all, observational studies.11 12 13 14 15 16 17 18 19 Indeed, in the five year interim analysis of a large observational study using the Kaiser Permanente Northern California database,4 the use of pioglitazone for 24 months or more was associated with an increased risk of bladder cancer (hazard ratio 1.4, 95% confidence interval 1.03 to 2.0). However, in the final analysis of the Kaiser Permanente Northern California study, which used the same cohort4 with follow-up extended to 10 years, the use of pioglitazone was no longer significantly associated with an increased risk of bladder cancer in a duration-response fashion.20 These null findings are also consistent with those of another large multicohort study.19 The apparent heterogeneity in this literature may be due to methodological limitations, such as the inclusion of prevalent users,5 6 10 11 12 13 14 18 time lag bias,15 immortal time bias,10 14 18 and no consideration of disease latency.8 10 12 17 18
Given these discrepant findings, the methodological shortcoming of previous studies examining this association, and the apparent loss of an association in studies with longer follow-up,20 additional studies are needed to investigate further the association between pioglitazone and bladder cancer. In a large, population based study we assessed the association between the use of pioglitazone and bladder cancer in people with type 2 diabetes.
This study was conducted using the United Kingdom Clinical Practice Research Datalink (CPRD). This database contains the complete primary care medical record of more than 13 million people.21 The Read code classification is used to record medical diagnoses and procedures, and a coded drug dictionary based on the UK Prescription Pricing Authority Dictionary is used to record prescriptions. The CPRD collects information on anthropometric variables such as body mass index and lifestyle variables such as smoking. Data collected in the CPRD have been previously validated and shown to be of high quality.22 23 Furthermore, cancer diagnoses have been found to be highly consistent with those recorded in the UK national cancer data repository.24
We assembled a base cohort composed of all people newly treated for type 2 diabetes, defined as receiving a first ever prescription for a non-insulin antidiabetic drug (metformin, sulfonylureas, prandial glucose regulators, thiazolidinediones, acarbose, dipeptidyl peptidase-4 (DPP-4) inhibitors, glucagon-like peptide (GLP-1) agonists, sodium-glucose cotransporter-2 (SGLT2) inhibitors) between 1 January 1988 and 31 July 2013. Patients were required to be at least 40 years of age and to have at least one year of CPRD medical history before that first prescription. We excluded patients prescribed insulin any time before their first non-insulin antidiabetic prescription (as these may represent those with an advanced form of type 2 diabetes), and patients with a diagnosis of gestational diabetes and polycystic ovary syndrome (as these are other indications for metformin).
Using the base cohort, we identified all patients who initiated a new antidiabetic drug class on or after 1 January 2000 (the year pioglitazone and rosiglitazone entered the UK market) until 31 July 2013. These patients included those newly treated with an antidiabetic drug class, as well as those who switched to or added-on an antidiabetic drug class not previously used in their treatment. Cohort entry was defined by the date of this new prescription. We excluded all patients with a diagnosis of bladder cancer (including malignant, in situ, and benign lesions) at any time before cohort entry, as well as those with less than one year of follow-up after cohort entry. The latter was necessary for latency considerations, as short term drug use are unlikely to cause incident bladder cancer.
All patients were followed from the year after cohort entry until a first ever diagnosis of bladder cancer (malignant and in situ), or censored on death from any cause, end of registration with the general practice, or end of the study period (31 July 2014), whichever occurred first.
Use of thiazolidinediones
In the models we entered the use of thiazolidinediones as a time varying variable and classified it according to one of the four mutually exclusive categories: pioglitazone use, rosiglitazone use, pioglitazone and rosiglitazone use (mainly switchers), and no thiazolidinedione use. Patients were considered unexposed to thiazolidinediones until the time of the first thiazolidinedione prescription and thereafter considered exposed, after accounting for a one year lag period. This lag period was necessary to take into account a latency time window and to minimise possible detection bias around the time of treatment initiation. This was considered the primary exposure definition.
In secondary analyses, we determined whether there was a duration-response and dose-response relation between pioglitazone and incidence of bladder cancer. The duration-response relation was assessed in terms of cumulative duration of use, which was defined, in a time dependent fashion, as the total number of years of use, calculated by summing the durations of all prescriptions received between cohort entry and the time of the event. This variable was then classified using the same categories used in the interim analysis of the Kaiser Permanente Northern California study4: <1 year, 1 to 2 years, and >2 years of use. We also assessed cumulative duration on a continuous scale using a restricted cubic spline model with five knots.25 Dose-response was assessed in terms of cumulative dose, which was calculated in a time dependent fashion as the sum of all doses received up until the date of the event. This variable was also categorised using the same cut-offs used in previous studies4 5: ≤10 500 mg, 10 501-28 000 mg, and >28 000 mg. We assessed the linear trend for cumulative duration of use and dose by considering these variables as continuous in the models.
For comparison purposes we also assessed whether there was a duration-response and dose-response relation with rosiglitazone, in terms of cumulative duration of use (categorically (<1 year, 1 to 2 years, and >2 years) and continuously using restricted cubic spline modelling) and cumulative dose (categorised on the basis of the distribution of use in thirds in the cohort).
All models were adjusted for several variables measured at cohort entry: age, sex, year of cohort entry, body mass index (<30 kg/m2, ≥30 kg/m2, unknown; last measure before cohort entry), smoking status (ever, never, unknown), alcohol related disorders (based on diagnoses for alcohol misuse, alcoholic cirrhosis of the liver, alcoholic hepatitis and failure, and other related disorders), haemoglobin A1c (≤7.4%, >7.4%, unknown; last recorded value before cohort entry), duration of treated diabetes (defined as the time between the first ever non-insulin prescription and cohort entry), previous bladder conditions (cystitis and bladder stones), history of cancer (other than non-melanoma skin cancer), presence of at least one urine protein test in the year before cohort entry, and Charlson comorbidity score26 (adapted so as not to include previous cancer, to avoid duplicate adjustment).
Descriptive statistics were used to summarise the characteristics of pioglitazone users, rosiglitazone users, and non-thiazolidinedione users at cohort entry. We calculated crude incidence rates of bladder cancer, with 95% confidence intervals based on the Poisson distribution overall and for each exposure category.
Time dependent Cox proportional hazards models, with duration of follow-up as the underlying time axis, were used to estimate adjusted hazard ratios and 95% confidence intervals for bladder cancer associated with the use of pioglitazone compared with no thiazolidinedione use. We also conducted two secondary analyses to assess whether there were duration-response and dose-response relations with pioglitazone and risk of bladder cancer in terms of cumulative duration of use and cumulative dose (as defined previously). Identical analyses were done for rosiglitazone. For all models we used the model proposed by Fine and Gray to account for competing risks due to death from any cause.27 We examined the Schoenfeld residuals for the time fixed covariates and found no important departures from the proportional hazards assumption.
We conducted nine sensitivity analyses to assess the robustness of our findings. Firstly, given uncertainties related to the latency time window, we repeated the primary analysis with lag periods of zero and two years. Secondly, we repeated the primary analysis after considering a stricter definition for drug use based on receiving at least four prescriptions within a 12 month window. Thirdly, we repeated the analysis after excluding patients with a history of bladder conditions at any time before cohort entry and censoring on a new diagnosis during follow-up. Fourthly, we repeated the primary analysis after additionally censoring on a new diagnosis of benign bladder lesions, in situ bladder cancer, liver failure, and heart failure (the last two were also additional exclusion criteria). Indeed, thiazolidinediones are contraindicated or not recommended for the two last conditions, the presence of which may lead to thiazolidinedione discontinuation or may influence treatment decisions. Fifthly, in 2011 several regulatory actions were issued because of the potential association between pioglitazone and bladder cancer.28 We performed a sensitivity analysis censoring follow-up to 31 December 2010, as it is possible that patients starting or continuing pioglitazone after that date may have been more carefully screened for bladder cancer. Sixthly, we repeated the primary analysis using multiple imputation for variables with missing values (that is, body mass index, smoking, and haemoglobin A1c).29 30 Seventhly, we additionally adjusted the models for the time dependent use of other antidiabetic drugs (metformin, sulfonylureas, incretin based drugs (GLP-1 analogues or DPP-4 inhibitors), insulin, and other oral hypoglycaemic drugs) during follow-up, lagged by one year for latency considerations. Eighthly, to account for potential time dependent confounding during the 14.5 year study period, we repeated the primary analysis using a marginal structural Cox proportional hazards model with inverse probability of treatment and censoring weighting (see the supplementary file for a detailed description of this method). Finally, we assessed the strength of an unmeasured confounder needed to move the estimated hazard ratio to the null using the “rule out” method proposed by Schneeweiss.31
Head to head comparison of pioglitazone with rosiglitazone
To assess further whether an association between pioglitazone and bladder cancer is a drug specific compared with a class effect, we conducted two additional analyses that directly compared pioglitazone with rosiglitazone. In the first approach, we contrasted the use of pioglitazone with the use of rosiglitazone by repeating our primary analysis with the latter as the reference category. In the second approach, we used the study cohort to assemble a subcohort of patients starting pioglitazone or rosiglitazone between 1 January 2000 and 31 July 2013, with follow-up until 31 July 2014. As with the primary analysis, all patients were required to have at least one year of follow-up after their first prescription for a thiazolidinedione. Consequently, cohort entry was set as the year after the first thiazolidinedione prescription during the study period. All patients were followed until a first ever diagnosis of bladder cancer, or censored on death from any cause, switching to another thiazolidinedione, end of registration with the general practice, or end of the study period, whichever occurred first. The model was adjusted for high dimensional propensity scores,32 which included the prespecified variables listed previously along with another 500 empirically defined variables measured at the time of the first thiazolidinedione prescription. All analyses were conducted with SAS version 9.4 (SAS Institute, Cary, NC).
No patients were involved in setting the research question or the outcome measure, nor were they involved in developing plans for design or implementation of the study. No patients were asked to advise on interpretation or writing up of results. There are no plans to disseminate the results of the research to study participants or the relevant patient community.
A total of 145 806 patients met the study inclusion criteria (see supplementary figure 1). Overall, the cohort was followed for a mean of 4.7 (SD 3.4) years, generating 689 616 person years of follow-up. Overall, 622 patients received a diagnosis of bladder cancer during follow-up, yielding a crude incidence rate of 90.2 (95% confidence interval 83.2 to 97.6) per 100 000 person years. Among patients with an event, the median time between cohort entry and an incident diagnosis of bladder cancer was 4.4 (interquartile range 2.5-6.5) years.
Table 1⇓ presents the characteristics of the cohort overall and stratified by pioglitazone users versus non-thiazolidinedione users at baseline. Compared with non-thiazolidinedione users, pioglitazone users were less likely to be obese but more likely to have increased haemoglobin A1c levels, to have undergone urine protein testing before cohort entry, had a longer duration of treated diabetes, and were more likely to have previous bladder conditions. Pioglitazone users were also more likely to have used sulfonylureas and less likely to have received metformin compared with non-users of thiazolidinediones. The baseline characteristics of rosiglitazone users are similar and shown in supplementary table 1.
Table 2⇓ shows the results of the primary and secondary analyses for pioglitazone. Compared with no thiazolidinedione use, the use of pioglitazone was associated with an increased risk of incident bladder cancer (121.0 v 88.9 per 100 000 person years; adjusted hazard ratio 1.63, 95% confidence interval 1.22 to 2.19). A duration-response relation was observed (P<0.01 for trend) with use of pioglitazone for more than two years associated with an increased risk of bladder cancer (adjusted hazard ratio 1.78, 95% confidence interval 1.21 to 2.64). In the restricted cubic spline analysis, the risk of bladder cancer was increased after 1.8 years of pioglitazone use, and continued to increase with longer durations of use, although this did not achieve statistical significance owing to a relatively small number of events among patients with longer duration of use (fig 1⇓). A dose-response relation was also present (P=0.01 for trend), with cumulative doses less than 10 500 mg (adjusted hazard ratio 1.63, 95% confidence interval 1.02 to 2.60) and more than 28 000 mg (1.70, 1.04 to 2.78) being associated with an increased risk of bladder cancer.
Overall, the use of rosiglitazone was not associated with an increased risk of incident bladder cancer (86.2 v 88.9 per 100 000 person years, adjusted hazard ratio 1.10, 95% confidence interval 0.83 to 1.47; see supplementary table 2). Similarly, there was no evidence of a duration-response relation in terms of cumulative duration of use when it was classified as a categorical variable (P=0.7 for trend; see supplementary table 2) or when it was considered as a continuous variable (see supplementary figure 2). Finally, there was no evidence of a dose-response relation in terms of cumulative dose (P=0.7 for trend; see supplementary table 2).
The results of the sensitivity analyses are summarised in figure 2⇓ and presented in supplementary tables 3 and 4. In all sensitivity analyses, the use of pioglitazone was consistently associated with an increased risk of bladder cancer, with adjusted hazard ratios ranging between 1.46 and 1.76. In contrast, the use of rosiglitazone was not associated with an increased risk of bladder cancer, with adjusted hazard ratios ranging between 1.01 and 1.16 and all estimates accompanied by 95% confidence intervals that included unity. Supplementary Figure 3 shows the exposure-confounder and confounder-disease associations (right of the curve) necessary to reduce the observed hazard ratio of 1.63 down to the null.
Head to head comparison of pioglitazone with rosiglitazone
Table 3⇓ presents the results of the comparisons between pioglitazone and rosiglitazone. In the first of these analyses (main model), pioglitazone use compared with rosiglitazone use was associated with an increased risk of bladder cancer (adjusted hazard ratio 1.48, 95% confidence interval 1.01 to 2.16). Similar findings were observed in the second analysis conducted within the thiazolidinedione subcohort (hazard ratio adjusted for high dimensional propensity scores 1.46, 95% confidence interval 0.94 to 2.27; see supplementary figure 4 and supplementary table 5 for cohort description).
In this large population based cohort study with up to 14.5 years of follow-up, pioglitazone was associated with an overall 63% increased risk of incident bladder cancer. There was also evidence of a duration-response and dose-response relations. In contrast, rosiglitazone was not associated with an increased risk of bladder cancer either overall or by cumulative duration of use and dose. Our findings remained consistent in several sensitivity analyses.
Strengths and weaknesses of this study
This study has several strengths. Firstly, we assembled a large population based cohort of patients newly treated with antidiabetic drugs and followed for up to 14.5 years, thus enabling the identification of a substantial number of patients with bladder cancer. Secondly, the inclusion of new users eliminated biases related to prevalent users.33 Thirdly, we considered a lag period to account for a minimum latency between use of thiazolidinediones and the development of bladder cancer. Fourthly, we defined exposure in a time dependent fashion, thereby eliminating immortal time bias.34 Fifthly, all analyses took into account competing risks due to deaths from any cause, an important consideration given the cardiovascular risk reported for thiazolidinediones in previous studies.1 Finally, the results remained consistent in several sensitivity analyses, thus confirming the robustness of our findings.
This study has some limitations. Firstly, residual confounding from unmeasured variables (eg, diet, physical activity, occupational exposure, pelvic radiation, family history of cancer, and race/ethnicity) is possible. However, the rule out method31 shows that a hypothetical unmeasured confounder would need to be strongly associated with both the exposure (odds ratio >3.7) and the outcome (relative risk >5.0) to move the point estimate down to the null. As the aforementioned variables are modestly associated with the outcome and it is unclear if they are associated with the exposure, we do not believe that residual confounding is a likely explanation for the observed association. Secondly, misclassification of drug use is possible, since the CPRD records prescriptions written by general practitioners and not those written by specialists. However, although some specialists may have been responsible for patients starting thiazolidinediones, general practitioners are likely to have been those prescribing repeat prescriptions for these drugs. Thus, misclassification is likely to minimal and, if present, would lead to an underestimation of the association. Finally, although cancers of the urinary tract have been shown to be well recorded in the CPRD,24 misclassification is possible. However, we expect this potential misclassification to be non-differential between patients using the different antidiabetic drugs included in the study.
Comparison with other studies
Several observational studies have investigated the association between pioglitazone and bladder cancer.4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 Overall, these studies have generated conflicting findings, with seven reporting statistically significant increased risks4 5 6 7 8 9 10 (ranging from 20% to 225%) and nine reporting null associations.11 12 13 14 15 16 17 18 19 The discrepancy between these studies is likely due to certain methodological shortcomings. Indeed, in three studies,10 14 18 the definition for drug use may have introduced immortal time bias, a bias resulting from the misclassification of unexposed person time as exposed person time, which may have led to a spurious underestimate of the association. In another study,15 time lag bias was introduced by comparing pioglitazone with insulin, the latter being a drug typically used at a more advanced stage of the disease, where the risk of cancer, including bladder cancer, may be higher. Prevalent users of antidiabetic drugs were included in 11 studies,4 5 6 10 11 12 13 14 18 19 20 which can be problematic in this context given the relatively rapid onset of bladder cancer after starting pioglitazone. Finally, in five studies,8 10 12 17 18 a minimum time between starting pioglitazone and the diagnosis of bladder cancer was not considered in the analyses, an important consideration given the latency of bladder cancer.
Our findings are consistent with those of the recently published Insulin Resistance Intervention after Stroke (IRIS) trial, which randomised 3895 people without diabetes to either pioglitazone or placebo.35 After a median follow-up of 4.8 years, pioglitazone was associated with a decreased risk of a composite endpoint of stroke or myocardial infarction (hazard ratio 0.76, 95% confidence interval 0.62 to 0.93). However, there was an increased number of bladder cancer events in the pioglitazone group compared with the placebo group (12 (0.6%) v 8 (0.4%), respectively).35 This imbalance was observed despite efforts to exclude patients with a history of or at high risk of bladder cancer (that is, defined by the presence of macroscopic haematuria, use of cyclophosphamide, or previous radiation to the pelvis).36 Thus, although this imbalance did not reach statistical significance,35 it mirrors the imbalance observed in the PROactive trial2 and is consistent with the effect size reported in our study.
Biological plausibility and implications
The biological plausibility of a rapid development of bladder cancer after starting pioglitazone has been debated, since many events observed in the PROactive trial occurred within one year of starting treatment.37 It is possible that these were prevalent cases and not attributable to pioglitazone,37 or promoted by pioglitazone in patients susceptible to developing bladder cancer.38 In our study, the use of a one year lag period ensured that all bladder cancer events had to occur at least one year after starting treatment. However, in sensitivity analyses, removing the lag period attenuated the hazard ratio (1.49, 95% confidence interval 1.13 to 1.97), whereas applying a two year lag period increased the hazard ratio (1.73, 1.26 to 2.39). Moreover, when assessed in a restricted cubic spline model, the risk tended to increase with longer durations of use. Taken together, our findings do not rule out a tumour promoting effect but also suggest that the risk may increase with longer use.
An important finding of our study is the absence of an association between rosiglitazone and bladder cancer. It is important to note that both pioglitazone and rosiglitazone entered the UK market the same year (2000) and both were intended for the same target population.39 Given their similarities, it is unlikely that confounding by indication or detection bias can explain the association observed with pioglitazone. In the head to head comparison, pioglitazone was associated with close to a 50% increased risk of bladder cancer compared with rosiglitazone. Of note, although the biological mechanism for pioglitazone induced bladder cancer is not clear, this imbalance in the risk of bladder cancer between these two thiazolidinediones could likely be explained by pharmacological differences. Indeed, unlike rosiglitazone, which is selective for the peroxisome proliferator activated receptor (PPAR) γ, pioglitazone has a dual PPARα/γ activity.40 41 This is particularly important, as PPARα/γ activation in rat models has been shown to increase the expression of carcinogenic biomarkers in the bladder, which has not been observed with the selective activation of PPARγ.42 43 44 Although differences in PPAR activity are possible explanations for the observed association, additional studies are needed to better understand the biological mechanism behind the possible pioglitazone specific effect on the bladder.
The use of pioglitazone is associated with an increased risk of bladder cancer, which varies in a duration dependent and dose dependent fashion. In contrast, rosiglitazone was not associated with an increased risk of bladder cancer in any analysis, suggesting the risk is drug specific and not a class effect.
What is already known on this topic
The association between the use of pioglitazone and bladder cancer is controversial, with studies reporting contradictory findings
Additional observational studies with longer follow-up are needed to assess whether this drug is associated with an increased risk of bladder cancer
What this study adds
In this large population based study, the use of pioglitazone was associated with an overall 63% increased risk of bladder cancer, with the risk increasing with increasing duration of use and dose
In contrast, the use of rosiglitazone was not associated with an increased risk, with no evidence of a duration-response or dose-response relation
These findings suggest that the association observed with pioglitazone is likely to be a drug specific and not a class effect
KBF holds a Canadian Institutes of Health Research new investigator award, and RWP is supported by a Chercheur-National Award of the Fonds de Recherche du Quebec-Santé (FRQS; Quebec Foundation for Health Research).
Contributors: MT, KBF, OHY, RWP, and LA conceived and designed the study. LA acquired the data. MT, HY, and LA carried out the statistical analysis. All authors analysed and interpreted the data. MT drafted the manuscript. All authors critically revised the manuscript for important intellectual content. LA obtained funding. LA supervised the study and is the guarantor.
Funding: This study is funded by the Canadian Institutes of Health Research.
Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf and declare: this study was funded by the Canadian Institutes of Health Research; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work..
Ethical approval: The study protocol was approved by the independent scientific advisory committee of the CPRD (protocol 11_099A) and by the research ethics board of the Jewish General Hospital, Montreal, Canada.
Data sharing: No additional data available.
Transparency: The lead author (LA) affirms that this manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned (and, if relevant, registered) have been explained.
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 3.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/3.0/.