Article Text

Download PDFPDF

Lipid-lowering drugs and risk of new-onset diabetes: a cohort study using Japanese healthcare data linked to clinical data for health screening
  1. Nobuhiro Ooba1,
  2. Soko Setoguchi2,
  3. Tsugumichi Sato3,
  4. Kiyoshi Kubota4
  1. 1 Department of Clinical Pharmacy, Nihon University School of Pharmacy, Funabashi, Chiba, Japan
  2. 2 Institute for Health, Health Care Policy and Aging Research, Rutgers University and Rutgers Robert Wood Johnson Medical School, New Brunswick, New Jersey, USA
  3. 3 Department of Pharmacy, Tokyo University of Science, Yamazaki, Noda, Chiba, Japan
  4. 4 NPO Drug Safety Research Unit Japan, Tokyo, Japan
  1. Correspondence to Dr Kiyoshi Kubota; kubotape-tky{at}umin.net

Abstract

Objective To investigate whether lipid-lowering drugs are associated with new-onset diabetes after adjusting for baseline clinical risk factors for diabetes.

Design A retrospective cohort study.

Setting Japanese employees of large corporations and their dependents using health insurance claims data linked to clinical and laboratory data for annual health screenings.

Participants All persons aged 20 to 74 years with dyslipidaemia between 1 January 2005 and 31 March 2011. We defined the index date as the first date when the person met the criteria for dyslipidaemia. Persons were excluded if they had lipid-lowering drugs, or had a diagnosis, a treatment or a laboratory test result (haemoglobin A1c ≥6.5% or fasting blood glucose ≥126 mg/dL) indicating diabetes during the 6-month period before the index date.

Main outcome measures New-onset diabetes.

Results We identified 68 620 persons with dyslipidaemia. During the mean follow-up period of 1.96 years, 3674 persons started treatment with a lipid-lowering drug: 979 with a low potency statin, 2208 with a high potency statin and 487 with a fibrate. Of 3674 new users of a lipid-lowering drug, 3621 had a period of non-use of any lipid-lowering drugs before starting a lipid-lowering drug. Among statin users, the incidence rate of new-onset diabetes was 124.6 per 1000 person-years compared with 22.6 per 1000 person-years in non-users. After adjusting for confounding factors including clinical data in health screening using Cox proportional hazards models, the HR was 1.91 (95% CI 1.38 to 2.64) for low potency statins and 2.61 (2.11 to 3.23) for high potency statins.

Conclusion The use of statins was associated with a 1.9-fold to 2.6-fold increase in the risk of new-onset diabetes in a Japanese population of working age, despite adjusting for clinical risk factors for diabetes.

  • Cardiac Epidemiology
  • Epidemiology
  • Adverse Events

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

Statistics from Altmetric.com

Request Permissions

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

Strengths and limitations of this study

  • This was not a randomised controlled trial but a retrospective cohort study using health insurance claims data linked to clinical and laboratory data for annual health screenings and subject to the confounding bias even if the bias was adjusted as much as possible.

  • Our study population included those who worked in relatively large corporations and their relatives, aged between 20 and 74 years, which means our findings might not be readily generalised to those working in small corporations or self-employed individuals and old people aged ≥75 years.

  • The diagnosis code for diabetes was not validated in our database, and therefore the outcome in sensitivity analyses was defined using an antidiabetic drug with or without the diagnosis code for diabetes.

  • Although the information on some of potential confounders such as smoking, alcohol consumption and physical exercise was not available, we could adjust for several important risk factors of diabetes such as HbA1c and fasting blood glucose essential in diagnosis of diabetes, as well as body mass index, blood pressure, uric acid, low-density lipoprotein cholesterol cholesterol, high-density lipoprotein cholesterol cholesterol, triglyceride and urinary glucose.

  • Our study suggested that the use of statins was associated with an increased risk of new-onset diabetes, compared with the period of non-use, despite adjusting for risk factors for diabetes (eg, HbA1c and body mass index) obtained in the annual health screening.

  • We also found that the risk of new-onset diabetes using high potency statins (atorvastatin, rosuvastatin and pitavastatin) was higher than that of low potency statins (fluvastatin, pravastatin and simvastatin).

Introduction

Hyperlipidaemia is an established risk factor for cardiovascular diseases. In Japan, the estimated number of patients with hyperlipidaemia was approximately 1.9 million at 2011.1 Statins are widely used for primary and secondary prevention of cardiovascular diseases2 and have a good safety profile in general, although rare adverse reactions such as liver dysfunction and rhabdomyolysis can occur.3 According to recent studies, the use of statins is associated with an increased risk of new-onset diabetes.4–6 As with hyperlipidaemia, diabetes is an established risk factor for cardiovascular diseases,7 and therefore the possible increased risk of new-onset diabetes by the drugs used to treat hyperlipidaemia may attenuate the effect of reducing the risk of cardiovascular diseases.

Meta-analyses of randomised trials report that the use of statins was associated with a modestly increased risk of diabetes (overall OR=1.09, 95% CI: 1.02 to 1.17) and the OR for individual statins ranged from 0.98 to 1.18.4 However, randomised trials may not accurately represent actual practice or patients. The results from observational studies are not consistent. A study using data from the women’s health initiative found that high potency statins (HR 1.45) and low potency statins (HR 1.48) were similarly associated with the increased risk of diabetes.5 On the other hand, in a study using the administrative healthcare database of Ontario,6 high potency statins had a slightly higher risk than that of low/moderate potency statins (HR 1.22 vs 0.97/1.11). The association between fibrates and new-onset diabetes indicated inconsistent results in a small number of previous studies.8 9

One potential reason for the inconsistency in assessing the magnitude of diabetes risk, arising from previous observational studies, may be due to the unavailability of information on the important confounders, such as cholesterol level, blood pressure, body mass index (BMI), haemoglobin A1c (HbA1c) and fasting blood glucose (FBG). The inability to control for these factors that potentially introduce confounding by indication or severity may have created an association between the use of a statin with new-onset diabetes variably targeting different study populations in previous studies. It is also unclear whether the risk of new-onset diabetes is specific to statins or common to lipid-lowering drugs (LLDs) including fibrates. In addition, there have been few reports on the association between statins and the incidence of diabetes in the Asian population.10 11 We conducted a cohort study using claims data linked to the health screening data including blood pressure, BMI, HbA1c, FBG and other laboratory test results to assess the association between the use of LLDs and new-onset diabetes, adjusted for each of these factors.

Methods

Data sources

We used three types of data: claims data, data derived from enrolment and data of regular health screening for beneficiaries of dozens of corporate health insurances for workers in private firms and their dependents. Under the Industrial Safety and Health Act, in addition to the Act on Assurance of Medical Care for Elderly People,12 all the insurers are required to provide the regular health screening and health guidance to beneficiaries and maintain and use the health screening data for the health guidance. The three types of data (claims, enrolment and health screenings) are linked by the insurers who provide the health guidance to beneficiaries and evaluate the heath guidance outcome. The data for 215 307 beneficiaries aged between 20 and 74 years in the study period (1 January 2005 to 31 March 2011) were collected and maintained by Japan Medical Data Centre Co.13 Information on healthcare utilisation, such as outpatient visits and hospitalisation, drugs, medical procedures and diagnoses, was available from the claims data. With respect to drug information14 within the claims, more than 20 000 local codes are used to specify the trade name and dosage form for all approved drugs. Codes for medical procedures include order of laboratory tests such as HbA1c, triglyceride (TG), total cholesterol, low-density lipoprotein (LDL) cholesterol and high-density lipoprotein (HDL) cholesterol. The 10th revision of the International Classification of Diseases was used for diagnoses. The enrolment data contained the dates of enrolment and disenrollment. Annual health screening is mandatory for employees irrespective of whether they have already had cardiovascular and other diseases under the Industrial Safety and Health Act, and optional for their dependents.15 The health screening data contained values of BMI, systolic and diastolic blood pressure, total cholesterol, LDL cholesterol, HDL cholesterol, TG, HbA1c, FBG, serum creatinine (SCre), uric acid (UA), urinary glucose and urinary protein.

Study cohort

We identified 94 630 subjects with or without prior cardiovascular diseases who had either total cholesterol ≥220 mg/dL, LDL cholesterol ≥140 mg/dL, HDL cholesterol <40 mg/dL or TG ≥150 mg/dL measured at the annual health screening period between 1 January 2005 and 31 March 2011 (figure 1).16 17 We defined the index date as the first date when one of these criteria was met and the baseline period as the period when the data on claims and health screening were available, after the date of enrolment but before the index date. Subjects were excluded if they had a baseline period of <6 months, had a diagnosis of diabetes, had used an antidiabetic drug or had a test result indicating diabetes (HbA1c≥6.5% or FBG≥126 mg/dL) within the baseline period. Subjects were also excluded if they had an LLD during the 6 months before the index date, if neither HbA1c nor FBG was available in the baseline period or if two or more LLDs prescribed at one time were filled on the index date.

Figure 1

Study flow diagram. FBG, fasting blood glucose; HbA1c, haemoglobin A1c; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; TC, total cholesterol; TG, triglyceride.

Exposures and outcomes data

We categorised person-times of the cohort members into non-use of LLDs and new use of statins or fibrates (bezafibrate and fenofibrate). Statins were further classified into the low (simvastatin, fluvastatin and pravastatin) and high potency statins (atorvastatin, rosuvastatin and pitavastatin),18 as the risk of diabetes may be different between low and high potency statins.6 19

New-onset diabetes defined as a diagnosis code of diabetes or where an oral antidiabetic drug or insulin was used after the index date. In the sensitivity analyses, we defined new-onset diabetes as the use of an antidiabetic drug after the index date only disregarding the diagnosis code. The end of the observation period for the outcome assessment was defined as either the date of disenrollment or the occurrence of new-onset diabetes or discontinuation of an LLD or switching to or addition of another category of LLD or 31 March 2011, whichever came first. Therefore, the period when the patient used a statin and a fibrate at the same time was excluded from the observation. Continuous use of the LLDs was defined using the grace period of up to 30 days before filling the next prescription. We considered that an LLD was used continuously when an LLD was switched to another LLD under the same category of LLDs (fibrates, high-potency statins and low-potency statins). If two or more LLDs were started on the same date (but after the index date), the observation was censored.

Covariates

The following covariates were obtained from the claims data and considered for confounding adjustment: age, sex, the use of other prescribed medications (systemic corticosteroids, thiazide diuretics, beta-adrenergic blocking agents, antipsychotics, antihypertensives and nicotinic acid) and a laboratory test order of HbA1c, TG, LDL cholesterol and HDL cholesterol during the period of 6 months preceding the index date, a prior diagnosis of myocardial infarction, chronic heart failure, cerebrovascular disease, renal, liver and pulmonary diseases, cancer, hypertension and polycystic ovarian syndrome. We selected the covariates that are known risk factors for diabetes (eg, systemic corticosteroid) and other major conditions (eg, cancer).6 7 The Charlson Comorbidity Index20 was also estimated and included as a covariate. Clinical covariates obtained from the health screening data included BMI, blood pressure (diastolic blood pressure and systolic blood pressure), FBG, HbA1c, levels of serum TG, cholesterol (LDL cholesterol and HDL cholesterol), UA and urinary protein and urinary glucose, estimated glomerular filtration rate (e-GFR; classified as e-GFR≥90, 60≤e-GFR<90, 30≤e-GFR <60, e-GFR<30) and presence of metabolic syndrome (MS). We calculated e-GFR using age, sex and measured SCre.21 MS was defined according to the criteria by Sicras et al criteria,22 where BMI (≥28.8 kg/m2) instead of waist circumference (>102 cm in males and >88 cm in females) was used to define MS, as data of waist circumference were not available in our data.

These covariates were measured during the baseline period for all subjects and again during the 6 months prior to the first prescription for lipid-lowering medications for those who initiated an LLD during the follow-up period.

Statistical analysis

We described the baseline characteristics measured during the 6 months prior to the index date for patients with the period of non-use, as well as the baseline characteristics measured on the date of new use of an LLD for new users of LLDs. We calculated crude incidence rates of new-onset diabetes for users and non-users of LLDs. For missing values in the health screening data, we first conducted a complete case analysis, followed by a multiple imputation using SAS V. 9.4 PROC MI. While approximately 20% of values for FBG, SCre, e-GFR and UA were missing, <8% were missing for BMI, diastolic or systolic blood pressure, HbA1c, LDL cholesterol, HDL cholesterol, TG, urinary protein and urinary glucose and MS. The Cox proportional hazards regression model was used to estimate unadjusted, age-sex-adjusted and multivariate-adjusted effects of the LLDs on new-onset diabetes. In all analyses, the reference used was the risk of new-onset diabetes in person-time of non-use of any LLDs.

We conducted multiple sensitivity analyses. First, we conducted the analyses including the period after switching of the LLD, although the observation was censored when an LLD was discontinued or switched in the primary analysis. In this sensitivity analysis, we used the latest value of the covariates assuming that they were not affected by the preceding LLD. Second, we conducted analyses using the grace period of up to 15 or 60 days. Third, to evaluate the robustness of our findings for potential outcome misclassification, the outcome was defined as the incident use of antidiabetic drugs only, disregarding whether the diagnosis code of diabetes was recorded. Lastly, to assess the impact of time-varying confounding, we used a marginal structural model (MSM) to analyse the data using the covariate measured after the date of new use of an LLD while the same definition as in the primary analysis was used for the outcome and grace period.23 As in the previous study using the database, we used the inverse-probability weight estimated on a monthly basis.24 We considered TG as a possible time-varying confounder because the average TG in the patients with fibrates (around 400 mg/dL) was much higher than that with statins (around 160 mg/dL), and the primary analysis revealed that the risk of new-onset diabetes was associated with TG (in mg/dL) measured on the date of a new use of an LLD (HR 1.0005, 95% CI: 1.0001 to 1.0008). All analyses were performed with SAS V. 9.4 (SAS Institute).

Results

Study population

We identified 94 630 persons with dyslipidaemia from laboratory test results in health screening. Of those, 68 620 (73%) persons met the eligibility criteria and were included in the study cohort (figure 1). Of these with dyslipidaemia, 64 946 (69%) had a period of non-use only, 3621 (4%) had a period of non-use prior to the period of LLD and 53 (0.06%) had a period of LLD only. Thus, 3674 (=3621 + 53) (4%) had a period of new use of an LLD (statin or fibrate) with or without a period of non-use. Among 3674 initiators of an LLD, 979 (27%) started with a low potency statin (pravastatin (n=726), fluvastatin (n=111) and simvastatin (n=142)), 2208 (60%) started with a high potency statin (atorvastatin (n=704), rosuvastatin (n=1016) and pitavastatin (n=488)) and 487 (13%) started with a fibrate (bezafibrate (n=300) and fenofibrate (n=187)). Table 1 shows the baseline characteristics. The mean age of patients with a period of non-use was younger than that of new users of an LLD. The proportion of men in new fibrate users was higher than that in new statin users.

Table 1

Baseline characteristics of patients with period of use/non-use of lipid-lowering drugs

LLDs and new-onset diabetes

The crude incidence rate of new-onset diabetes was 22.6 per 1000 person-years during the period of non-use while that for the use of LLDs was 105.5 for low potency statins, 133.1 for high potency statins and 99.2 per 1000 person-years for fibrates. Compared with non-use, the use of low or high potency statins or fibrates was associated with an increased risk of new-onset diabetes (table 2).

Table 2

Association between lipid-lowering drugs and new-onset diabetes

By the complete case analysis using the Cox proportional hazards regression, the risk of new-onset diabetes for high potency statins (HR 3.42, 95% CI 2.54 to 4.60) and low potency statins (2.04, 1.25 to 3.33) and fibratres (2.35, 1.08 to 5.15) was found to be higher compared with the risk during the period of non-use. The Cox proportional hazards regression for multiple imputation revealed that the risk of new-onset diabetes for high potency statins (adjusted HR 2.61, 95% CI 2.11 to 3.23) was 1.4 to 1.6 times higher than that of low potency statins (1.91, 1.38 to 2.64) and fibrates (1.64, 0.98 to 2.76) (table 2).

In table 3, the unadjusted and adjusted HRs of new-onset diabetes of the individual LLDs, compared with non-use, are shown. All of the individual statins except for simvastatin were associated with an increased risk of new-onset diabetes. The HR was higher than 1 for the individual fibrates, but the CIs were wide and crossed 1.00.

Table 3

Association between individual lipid-lowering drugs and new-onset diabetes

Figure 2 shows the results of the sensitivity analyses. The results were similar to those of the primary analysis when the periods after a total of 164 occasions of switching of an LLD were included as well as when the grace period was set as either 15 or 60 days. When the outcome was defined by the new use of antidiabetic drugs, the HR was greater than 1 but the CI of the HR was wider than that in the primary analysis.

Figure 2

Sensitivity analysis.*: primary analysis#: for the grace period of 15 days, no case was observed in those with fibrates. Outcome defined by diagnosis and/or drug; outcome defined by a diagnosis code of diabetes and/or new use of an antidiabetic drug; outcome defined by drug; outcome defined by new use of an antidiabetic drug disregarding the diagnosis code.

In the analysis using MSM, the results were similar to those in the primary analysis. The HR with the inverse-probability weight was 2.75 (2.23 to 3.39) for high potency statins, 1.67 (1.19 to 2.36) for low potency statins and 1.53 (0.87 to 2.71) for fibrates.

Discussion

In this retrospective cohort study using data of claims, enrolment and health screening, we found that the use of LLDs (statins and fibrates) was associated with an increased risk of new-onset diabetes, compared with non-use. The adjusted HR of new-onset diabetes between two fibrates (bezafibrate and fenofibrate) was similar, while the HR in the high potency statins (atorvastatin, rosuvastatin and pitavastatin) was higher than that in the low potency statins (fluvastatin, pravastatin and simvastatin). The HR for new-onset diabetes was varied from 1.5 to 3.1 in the individual statins. As the incidence rate of new-onset diabetes was estimated as 22.6 per 1000 person-years during the period of non-use in our study, the absolute increase of the rate potentially due to the use of statins, corresponding to HR of 1.5 to 3.1, may be approximately 10 to 45 per 1000 person-years.

Previous studies on the association between the use of statins and new-onset diabetes have been inconclusive. In some of the observational cohort studies,5 25 and a meta-analysis of randomised placebo-controlled trials,4 the risk of new-onset diabetes was not different between low and high potency statins. On the other hand, in a study using the pharmacy claims by Zaharan et al, the HR (1.3 to 1.4) of new-onset diabetes in high potency statins (atorvastatin and rosuvastatin) was slightly higher than the HR (1.0 to 1.1) in low potency statins (fluvastatin, pravastatin and simvastatin).26 Our findings suggest that the risk of new-onset diabetes in high potency statins may be higher than that of low potency statins, in line with the findings of Zaharan et al.26 These discrepancies may be explained by the different outcome definitions used among studies.

Our study did not find fibrates (fenofibrate and bezafibrate) to be associated with an increased risk of new-onset diabetes compared with that for non-use. However, we cannot rule out a possible increased risk due to small numbers of patients taking fibrates. There have been few studies on the risk of diabetes in patients prescribed with fibrates. In a randomised controlled trial on cardiovascular events and a claims-based study in Taiwan,8 9 fibrates were not associated with an increased risk of new-onset diabetes. However, the observational study by Flory et al indicated that the risk of new-onset diabetes due to bezafibrate use was different to that of other fibrates.27 In our findings, the risk of new-onset diabetes due to fibrates tended to be high. However, the number of patients with new-onset diabetes was small (19) and the HR had a wide CI. Therefore, it is quite possible that the higher risk of fibrates compared with non-use was obtained just by chance in the current study.

Both statins and fibrates are known to affect glucose homeostasis.28 29 The insulin-sensitising action of adiponectin results from a decrease in hepatic gluconeogenesis and an increase in muscle glucose transport.30 Hypoadiponectin may be associated with the incidence of diabetes.31 32 In a meta-analysis comparing fibrates and statins for elevation of circulating adiponectin, a non-significant negative trend was indicated between changes in plasma adiponectin and glucose levels.33 Fibrates34 and statins35 36 lowered the adiponectin level in some studies. However, another meta-analysis suggested that statin and fibrate therapy may increase circulating adiponectin.37 38

Our study has several limitations. First, this was not a randomised controlled trial but an observational cohort study and therefore we cannot exclude the possibility that the results are biased due to the residual confoundings although the potential confounding effects were adjusted in the analysis as much as possible. In particular, the information on some of known risk factors including smoking, alcohol consumption and physical exercise was not available and we could not adjust for them. Second, our study population included only those who are working in relatively large corporations and their dependents, aged 20 to 74 years. In addition, as the health screening was not compulsory for the employee’s relatives, the proportion of the relatives who had the health screening data was small (approximately 30%) and those who had the health screening may be systematically different from those who did not. For example, those who had the health screening may be more health orientated. Therefore, our findings might not be generalisable to those working in small corporations and self-employed, of older ages, and those who did not have the health screening. Third, the diagnosis code for diabetes was not validated in our database. However, the results of our sensitivity analyses, where the outcome was defined using an antidiabetic drug only, were similar to those in the primary analysis where the outcome was defined by a diagnosis code of diabetes or the use of an antidiabetic drug. In general, the outcome defined by the initiation of a drug may be better than the outcome defined by the diagnosis code.39 40 Fourth, the mean follow-up period (approximately 4 to 5 months) of LLDs in this study was relatively short, and we could not determine the risks associated with the longer use from the current study. A comparison with the previous studies with a longer follow-up period may also be difficult. However, in the population-based study by Carter et al,6 where the risk of new-onset diabetes was found to be increased among patients treated with statins, the median follow-up period was also short (<1 year). There were case reports where the patient had new-onset diabetes in 1 to 2 months after the administration of a statin.41 Fifth, the relatively high frequency of physician visits during the period of drug treatment as compared with that during the period of non-use might have introduced a detection bias. However, some of our findings such as the difference in the magnitude of the risk of new-onset diabetes between high and low potency statins could not be accounted for by the detection bias. Sixth, we could not examine the dose–response relationships, as the number of patients who had new-onset diabetes was not large enough to examine the details, such as any risk between the use of individual statins and dose–response relationships.

Our study has several strengths. We could adjust for HbA1c and FBG, essential in the diagnosis of diabetes, as well as for BMI, blood pressure, UA, LDL cholesterol, HDL cholesterol, TG, proteinuria and urinary glucose measured on the index date and on the date of new use of an LLD.7 These variables have been unavailable in most previous studies using claims data. Statins were associated with an increased risk of new-onset diabetes, even after being adjusted for all available confounders, including covariates, although there may be residual confounders as information was not available on some potential confounders such as smoking, alcohol consumption and physical exercise.

In conclusion, we found that the use of LLDs (low and high potency statins and potentially fibrates as well) in a Japanese workers and their dependents aged 20 to 74 was associated with new-onset diabetes, compared with non-use. The risk of new-onset diabetes in high potency statins was higher than that in low potency statins. Further studies are needed to understand the association between new-onset diabetes and fibrates. More studies are needed to evaluate the association between new-onset diabetes and the use of LLDs in the elderly, and in association with individual LLDs.

References

  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 20.
  21. 21.
  22. 22.
  23. 23.
  24. 24.
  25. 25.
  26. 26.
  27. 27.
  28. 28.
  29. 29.
  30. 30.
  31. 31.
  32. 32.
  33. 33.
  34. 34.
  35. 35.
  36. 36.
  37. 37.
  38. 38.
  39. 39.
  40. 40.
  41. 41.

Footnotes

  • Contributors All authors contributed to the conception and design of the study. NO, TS and KK contributed to the acquisition of the data. NO, SS and KK contributed to the analysis and the interpretation of the data. NO, SS and KK drafted the manuscript. NO and KK are the guarantors.

  • Funding This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.

  • Competing interests None declared.

  • Ethics approval This study was approved by the Ethical Committee of Nihon University School of Pharmacy (No. 14011).

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

  • Data sharing statement No additional data are available.