Article Text

Effect of low-protein diet on kidney function in diabetic nephropathy: meta-analysis of randomised controlled trials
  1. Uru Nezu1,
  2. Hiroshi Kamiyama2,
  3. Yoshinobu Kondo3,
  4. Mio Sakuma4,
  5. Takeshi Morimoto4,
  6. Shinichiro Ueda1
  1. 1Department of Clinical Pharmacology & Therapeutics, University of the Ryukyus, Okinawa, Japan
  2. 2Department of Endocrinology & Metabolism, Yokohama City University Graduate School of Medicine, Kanagawa, Japan
  3. 3Department of Endocrinology & Metabolism, Chigasaki Municipal Hospital, Kanagawa, Japan
  4. 4Department of Internal Medicine, Hyogo College of Medicine, Hyogo, Japan
  1. Correspondence to Dr Shinichiro Ueda; blessyou{at}med.u-ryukyu.ac.jp

Abstract

Objective To evaluate the effect of low-protein diet on kidney function in patients with diabetic nephropathy.

Design A systematic review and a meta-analysis of randomised controlled trials.

Data sources MEDLINE, EMBASE, Cochrane Library, ClinicalTrials.gov, International Standard Randomised Controlled Trial Number (ISRCTN) Register and University Hospital Medical Information Network-Clinical Trials Registry (UMIN-CTR) from inception to 10 December 2012. Internet searches were also carried out with general search engines (Google and Google Scholar).

Study selection Randomised controlled trials that compared low-protein diet versus control diet and assessed the effects on kidney function, proteinuria, glycaemic control or nutritional status.

Primary and secondary outcome measures and data synthesis The primary outcome was a change in the glomerular filtration rate (GFR). The secondary outcomes were changes in proteinuria, post-treatment value of glycated haemoglobin A1C (HbA1c) and post-treatment value of serum albumin. The results were summarised as the mean difference for continuous outcomes and pooled by the random effects model. Subgroup analyses and sensitivity analyses were conducted regarding patient characteristics, intervention period, methodological quality and assessment of diet compliance. The assessment of diet compliance was performed based on the actual protein intake ratio (APIR) of the low-protein diet group to the control group.

Results We identified 13 randomised controlled trials enrolling 779 patients. A low-protein diet was associated with a significant improvement in GFR (5.82 ml/min/1.73 m2, 95% CI 2.30 to 9.33, I2=92%; n=624). This effect was consistent across the subgroups of type of diabetes, stages of nephropathy and intervention period. However, GFR was improved only when diet compliance was fair (8.92, 95% CI 2.75 to 15.09, I2=92% for APIR <0.9 and 0.03, 95% CI −1.49 to 1.56, I2=90% for APIR ≥0.9). Proteinuria and serum albumin were not differed between the groups. HbA1c was slightly but significantly decreased in the low-protein diet group (−0.26%, 95% CI −0.35 to −0.18, I2=0%; n=536).

Conclusions Low-protein diet was significantly associated with improvement of diabetic nephropathy. The adverse effects of low-protein diet were not apparent such as worsening of glycaemic control and malnutrition.

This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/3.0/ and http://creativecommons.org/licenses/by-nc/3.0/legalcode

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.

Article summary

Article focus

  • Our research question is whether low-protein diet treatment is beneficial or not in patients with diabetic nephropathy.

  • Our hypothesis was that low-protein diet improves kidney function, but this effect is affected by patients’ actual compliance with diet treatment.

Key messages

  • A low-protein diet significantly improved kidney function in patients with diabetic nephropathy.

  • This effectiveness was observed only when patients’ diet compliance was fair; however, the successful treatment may not need to be as stringent as the current clinical guidelines recommend.

Strengths and limitations of this study

  • Unique and reasonable approach to minimise the bias by interstudy difference in patients’ compliance with diet treatment.

  • Quality of the evidence was not high when assessed by the Grading of Recommendations, Assessment, Development and Evaluation (GRADE) approach.

  • A few missing data were imputed, generating a risk of heterogeneity.

Introduction

Diabetic nephropathy is the leading cause of end-stage renal disease necessitating renal replacement therapy1 ,2 and is also associated with increased risk of cardiovascular mortality.3 It affects up to 40% of diabetic patients4 and the medical cost for treatment piles up to US$16.8 billion in the USA and US$1.2 billion in the UK each year.5

The progression of diabetic nephropathy can be slowed down by optimal glycaemic control 1 ,6 ,7 and that of blood pressure control by renin-angiotensin system blockade.8 ,9 As for the diet therapy, a low-protein diet (LPD) is recommended in clinical guidelines by the American Diabetes Association.10 ,11 This is based on animal studies and several small studies on humans. However, previously conducted randomised controlled trials (RCTs) have not consistently shown the benefits of LPD.

To elucidate this clinical question, several meta-analyses have been published on this topic.12–14 However, some of them included pre–post study and crossover trials which might have obscured the pure effect of intervention. In addition, RCTs of diet treatment have always been facing the difficulty of achieving sufficient compliance. In this view, we conducted a meta-analysis considering the interstudy variance of diet compliance with an additional study involving a considerable number of patients subsequently published.15

Methods

We carried out meta-analysis in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) Statement for systematic reviews and meta-analyses.16

Search strategy

We searched RCTs via PubMed, EMBASE, Cochrane Library, ClinicalTrials.gov, International Standard Randomised Controlled Trial Number (ISRCTN) Register and University Hospital Medical Information Network-Clinical Trials Registry (UMIN-CTR) from inception to 10 December 2012 to identify relevant citations. Internet searches were also carried out with general search engines (such as Google and Google Scholar). Studies written in English evaluating the effect of LPD comparing with control diet among diabetic patients were identified using the search terms ‘protein restriction’ OR ‘low protein diet’ AND ‘diabetes’.

Study selection

We assessed all the identified studies for the criteria of this meta-analysis. Two independent investigators (UN and HK) sorted out the potentially relevant studies first by title and abstract review, and finally judged the eligibility by full-text review. When discrepancies occurred, we discussed in a committee involving four investigators of our research group (UN, MS, TM and SU).

Inclusion criteria are as follows: published in full text, RCTs with a parallel design of LPD among patients with either type 1 diabetes mellitus (T1DM) or type 2 diabetes mellitus (T2DM) and any of the following outcomes are available; glomerular filtration rate (GFR), creatine clearance (CCr), proteinuria, albuminuria, glycated haemoglobin A1C (HbA1c) or serum albumin. RCTs of crossover design were excluded because of the possible carryover effect. As for the studies likely to have multiple reports, we selected the most recent publication after assessing their independency focusing on the patients’ background, intervention details, outcome settings and the results.

Data extraction

We extracted data related to published year, number of patients and their characteristics (age, gender, duration of diabetes mellitus and stages of diabetes nephropathy), details of the prescribed diet, intervention period. Also, we extracted data for patients’ compliance by integrating the data on actual protein intake (g/kg/day, g/day, mg/mg or energy per cent) evaluated in each study, then calculated the LPD to control ratio of actual protein intake (APIR). We utilised these results to assess study quality and subsequent subgroup analyses.

As primary outcome data, we extracted the mean change in GFR (ml/min/1.73 m2) or CCr (ml/min/1.73 m2) from baseline to the end of the diet intervention. As the secondary outcome, we extracted the mean change in proteinuria (g/24 h), albuminuria (mg/24 h, μg/min), urine albumin to creatine ratio (Cre, mg/mmol), post-treatment value of HbA1c (%) and post-treatment value of serum albumin (g/dl). Different digit numbers of proteinuria (mg/24 h) and albuminuria (g/24 h) were converted to fit the above scales.

We unified the notation for nephropathy stage into three ways; normoalbuminuria, microalbuminuria (incipient nephropathy) and macroalbuminuria (overt nephropathy).The value of HbA1c (%) extracted from Japanese articles was converted from Japan Diabetes Society (JDS) to National Glycohemoglobin Standardization Program (NGSP) by an equation announced by JDS as NGSP (%)=1.02 × JDS (%)+0.25%.17

Risk of bias assessment

Using the Cochrane Collaboration’s ‘risk of bias’ tool,18 we assessed the risk of bias of included studies. We assessed seven domains: (1) sequence generation, (2) allocation concealment, (3) blinding of participants and personnel, (4) blinding of outcome assessment, (5) incomplete outcome data, (6) selective outcome reporting and (7) other bias by patients’ diet compliance. Since this study aimed to investigate the clinical effect of dietary intervention which encourages patients’ lifestyle modification, we considered that patients’ diet compliance was the most critical factor to generate risk of bias. Therefore, we categorised studies with APIR over 0.9 as ‘high risk’.

Next, we scored the risk level of each domain from 0 to 2; high risk=2, unclear risk=1 and low risk=0. As an exception, we gave a score of 3 to ‘high risk’ of ‘other bias by diet compliance’ to lay more weight on this domain. Finally, we assessed the included studies’ overall risk of bias by the total score of the seven domains.

Quantitative data synthesis

We summarised results as the mean difference of continuous variables with 95% CIs and combined data by means of a random effects model with inverse variance weighting. GFR and CCr were used interchangeably, since it is commonly used as an estimate of GFR. Since proteinuria and albuminuria were measured in different measurements and scales, we used the standardised mean difference by dividing the mean value by the SD.

If SDs were missing, we obtained them by converting from alternative variance measures such as SEs, CIs and p values. When even such information was not available, we imputed the value using a technique by Follmann et al19 and Abrams et al.20 This technique utilises the correlation coefficient obtained from a study giving detailed information.

We subsequently conducted subgroup analyses for the main outcomes. Prespecified subgroups were based on the patients’ baseline characteristics (body mass index (BMI), type of diabetes and stages of diabetic nephropathy) and study methodology (intervention period, measurement index of proteinuria, overall assessment risk of bias and diet compliance assessed by APIR).

We used I2 statistics to assess statistical heterogeneity among studies. The possible publication bias was assessed by visual asymmetry of a funnel plot. We referred to the Cochrane Handbook V.5.1.018 for methodological guidance. We used Review Manager (RevMan) for Windows Software V.5.1.7 (the Nordic Cochrane Centre, Copenhagen, Denmark).

Quality of evidence

We graded the quality of evidence for the primary outcome using the Grading of Recommendations, Assessment, Development and Evaluation (GRADE) approach21 ,22 using GRADEpro software V.3.6 (for Windows; Jan Brozek, Andrew Oxman, Holger Schünemann, 2008). The quality of the evidence for GFR was low (see online supplementary table A).

Results

Search results

As shown in figure 1, we initially obtained 912 records through electrical database search. Eight hundred and sixty-four were excluded after evaluation of titles and abstracts. After removing 27 duplicates, we selected 21 full-text articles for detailed assessment for eligibility. Among these, we excluded eight studies: one study owing to lack of comparison,23 two studies of non-randomisation trial24 ,25 and three studies of crossover design.26–28 In addition, we excluded two studies which were likely to have multiple reports published by the same author group.29 ,30 We included the recent publication in analyses. Similarly, two studies by Pijls et al were likely to have multiple reports. However, we did not exclude the previous publication, since only that provided the postintervention value of HbA1c, the secondary outcome in our meta-analysis. Finally, we included 13 RCTs reporting the effects of LPD in diabetic patients.15 ,31–42

Figure 1

Flow diagram of the process for study selection.

Characteristics of included studies

The included studies evaluated the effects of LPD in 779 diabetic patients (209 T1DM and 555 T2DM) from Japan,15 Mexico,31 France,32 Italy,33 Australia,34 Denmark,35 Netherlands,36 ,37 ,39 South Africa,38 Italy41 and the USA.40 ,42 Study patients were middle-aged men and women, mostly obese or overweight (table 1). Mean duration of diabetes history was 18 years. T1DM accounted for six studies and T2DM for five studies. Two studies included both T1DM and T2DM patients and provided no separate information. The stage of diabetic nephropathy ranged from normoalbuminuria to macroalbuminuria. Eight studies provided distinct information generated from a single nephropathy stage; however, the remaining five studies reported only the combined results of two neighbour stages. Baseline GFR was 76 ml/min/1.73 m2 and HbA1c was 8.3% in average. An intervention period ranged from 3 to 60 months (18 months in median). All trials clarified random assignment and the methods were generally adequate (see online supplementary figure A and table B). However, allocation concealment was unclear in about half of the studies. With regards to blinding of intervention, only one study31 applied single-blind method. Although the outcome assessment was not blinded to the assessors in any of the studies, the risk of bias is considered to be small since the outcome is objective. Attrition bias was seen in variety. We considered the proportion in the number of incomplete patients and the reasons for dropping out, in order to see whether these were different across the intervention groups. Although selective reporting was not concerned, other biases by insufficient diet compliance were considered in four studies as we describe later.

Table 1

Characteristics of included studies

Diet prescription and compliance assessment

Table 2 shows the details of diet prescription and compliance assessment. The prescribed protein level was 0.6–0.8 g/kg/day in LPD and 1.0–1.6 g/kg/day in control. In five studies, patients in the control diet treatment group were instructed to continue their habitual diet instead of setting any numerical goal of protein intake.33 ,35–37 ,39

Table 2

Details of diet prescription and compliance assessment

Diet compliance was assessed in all trials. Ten studies measured 24 h urine urea nitrogen (24 h UUN) and calculated daily protein intake (g/kg/day).15 ,31–33 35–40 They used Maroni's formula, the gold standard of protein intake estimation considering nitrogen loss from a non-urine source.43 ,44 One study by Ben et al reported only the value of 24 UUN (g/day) without using this formula.42 Nine studies conducted additional or alternative assessment, such as spot UUN to Cre (UUN/Cre),34 4 h UUN,41 food questionnaire, food record or recall technique.15 ,31–33 ,37 ,39

APIR ranged from 0.44 to 1.07. When setting the cut-off value of APIR for sufficient compliance as 0.9, it was less than 0.9 in only nine studies (69.2%). Interestingly, all of the T1DM studies showed fair compliance (APIR <0.9).35 ,38–42 However, only two out of five T2DM31 34 study achieved fair diet compliance.31 ,34 Two studies including both T1DM and T2DM patients showed values of 1.07 and 0.69, respectively.32 ,33

Effects of LPD on kidney function

Eleven trials of 624 patients provided the change in kidney function assessed either by GFR or CCr. GFR was significantly increased by 5.82 ml/min/1.73 m2 after LPD (95% CI 2.30 to 9.33 ml/min/1.73 m2; figure 2). We found a significant heterogeneity across the studies (I2=92%, p<0.00001); however, the funnel plot showed no major asymmetricity (see online supplementary figure D).

Figure 2

Effects of low-protein diet on glomerular filtration rate.

Effects of LPD on proteinuria

Twelve studies of 634 patients provided sufficient information regarding change in proteinuria.15 ,31–36 38–42 The standard mean difference showed no significant change in proteinuria after LPD (−0.14, 95% CI −0.74 to 0.46; p=0.65; figure 3). Although we found heterogeneity across the studies (I2=91%, p<0.00001), the funnel plot showed no major asymmetricity (see online supplementary figure E).

Figure 3

Effects of low-protein diet on proteinuria.

Effects of LPD on glycaemic control

Glycaemic control was assessed by the absolute value of HbA1c after diet treatment. In eleven studies which provided sufficient information, HbA1c was slightly but significantly decreased after LPD (−0.26%, 95% CI −0.35 to −0.18; see online supplementary figure B). Although we found heterogeneity across the studies (I2=0%, p<0.00001), the funnel plot showed no major asymmetricity (see online supplementary figure F).

Effects of LPD on nutritional status

Nutritional status was assessed by the absolute value of serum albumin after diet treatment. Only four studies of 179 patients provided sufficient information.32 ,33 ,38 ,39 As a result, serum albumin was not changed after LPD (−0.18, 95% CI −0.53 to 0.17, p=0.32; see online supplementary figure C). Heterogeneity between trials was significant (I2=88%, p<0.00001), and the funnel plot showed asymmetrical appearance (see online supplementary figure G).

Subgroup and sensitivity analyses

Table 3 shows subgroup analysis according to clinical characteristics and study quality. There were significant differences in change in GFR between the subgroups based on nephropathy stage (p=0.03) and diet compliance (p=0.006). Specifically, GFR was improved in the subgroup of macroalbuminuria and subgroup of fair diet compliance (APIR <0.9). There was no significant difference between subgroups of BMI, type of diabetes, overall risk of bias and diet compliance. As for the change in proteinuria, there were significant differences between the subgroups of BMI (p<0.0001), type of diabetes (p=0.002), nephropathy stage (p=0.001) and measurement index of proteinuria (p<0.00001). There was no significant difference between the subgroups of intervention period, overall risk of bias and diet compliance. The post-treatment value of HbA1c was not differed across the subgroups. As for serum albumin, the sensitivity analysis excluding only one study by Dussol et al, which showed poor diet compliance (APIR 1.02), showed no significant change after LPD (−0.25, 95% CI −0.64 to 0.15; p=0.22, I2=91%, data not shown in table).

Table 3

Subgroup analyses and sensitivity analyses for clinical characteristics and study quality

We conducted sensitivity analyses seeking a pure effect of LPD on diabetic nephropathy. First, we excluded a subgroup of normoalbuminuria patients in the study by Velázquez et al,31 since its separate data were provided. The overall improvement was consistent with regard to GFR (6.31, 95% CI 2.47 to 10.15; p<0.00001, I2=92%) as well as HbA1c (−0.26, 95% CI −0.34 to −0.17; p=0.00001, I2=0%). Subsequently, we excluded an additional three studies that provided combined data of normoalbuminuria and microalbuminuria patients.34 ,36 ,37 In this way, patients without diabetic nephropathy were completely eliminated. As a result, improvement of GFR was still significant (6.32, 95% CI 2.52 to 10.52; p=0.001, I2=93%) and proteinuria improved significantly (−0.62, 95% CI −1.15 to −0.09; p=0.02, I2=84%).

Discussion

A statement of the principal findings

We found a protective effect of dietary intervention by LPD on the course of diabetic nephropathy by improving GFR and proteinuria. In addition, LPD did not worsen either glycaemic control or nutritional status.

Strength and limitations of the study

Although we searched only English publications, the included numbers of studies and patients were larger than those of any of the previous meta-analyses. Another strength is that we took a unique but reasonable approach to minimise the bias by interstudy difference in patients’ compliance to diet treatment. We proposed APIR as the common index that enables the compliance level to be compared across studies. The subgroup analysis based on APIR showed that LPD improved GFR only when intervention was sustainable, which is clinically reasonable.

In addition, APIR was also utilised in the sensitivity analysis of serum albumin, in which three studies with fair diet compliance consistently showed no worsening of nutritional status. Since the number of included studies was limited, the interpretation should be cautious. We need another large scale of RCTs to draw more accurate conclusions in terms of this issue on malnutrition.

This study has some more limitations. First, the quality of the evidence assessed for GFR was not high according to the GRADE approach. Two factors that lowered the grade were the inconsistency of the intervention and the indirectness of the outcome. Although the directions of the intervention were consistent across most studies, there was a small overlap in the CIs, and also the heterogeneity was not negligible. This inconsistency may be partly explained by the difference in the study protocol. GFR was measured in different ways in the RCTs included in this meta-analysis, as is also the case in clinical practice. As for the indirectness of outcome, GFR is a candidate surrogate marker to predict kidney failure or initiation of dialysis. However, we should not oversimplify that increasing GFR always represents a better prognosis of the patients’ kidney function because glomerular hyperfiltration may occur in early stages of renal damage as postulated by Brenner et al.45 In addition, the long-term clinical validity of GFR is not sufficient. We hope that more large-scale prospective studies or another meta-analysis will elucidate the effect of LPD on the change in GFR as well as on more long-term clinical outcomes such as mortality, dialysis or transplantation.

The second limitation was the overlaps of nephropathy stages in the subgroup analysis as shown in table 3. Since several studies enrolled patients in more than two stages of nephropathy without reporting the respective separated data, the subgroup analysis could not make clear subgroups without stage overlaps. Therefore, in order to seek a pure effect of LPD in diabetic nephropathy, we conducted sensitivity analyses by removing a subgroup comparison or the studies that included patients free from diabetic nephropathy. The analyses showed a consistently significant improvement in GFR. Although proteinuria was improved differently from the overall analysis, it is reasonable in two ways. First, it is logically impossible to reduce albuminuria in patients who have not suffered from albuminuria. Second, the relation is clinically compatible between reduction of proteinuria and improvement in kidney function.

The third limitation was that a part of the missing information for outcome was imputed, generating a risk of heterogeneity. We estimated unknown SDs by using the correlation coefficient obtained from included studies in this meta-analysis. However, the correlation coefficient was drawn from a study or a subgroup of macroalbuminuria.31 ,32 Therefore, this estimation might not have been proper substitutes for studies or subgroups of patients with normoalbuminuria or microalbuminuria. Sensitivity analysis excluding the study that enrolled patients with normoalbuminuria still includes microalbuminuric patients, and thus the interpretation should be cautious.

Comparison with other studies

There are three meta-analyses on this issue pooling the data from RCTs. The meta-analysis by Pedrini et al14 reported the beneficial effects of LPD; however, they combined RCTs and non-randomised crossover trials. In addition, they used a composite outcome of GFR or albuminuria. The meta-analyses by Pan et al12 and Robertson et al13 did not show significant effects on kidney function. The different result is explainable by the difference in the pooled study number and population size. The meta-analysis by Robertson et al especially have pooled the data from only seven RCTs, since they focused on a study including T1DM patients. Consistently, in our analysis, GFR in T1DM patients was improved but not statistically significantly.

Pan et al’s meta-analysis included two reports by Pijls et al.37 However, the patient’s background in these two reports was almost identical as shown in table 1. What was different was the number of patients and the intervention period, which was large and longer in a recent publication. Our reviewers discussed in the committee and concluded that the previous publication might be the interim analysis of a longer project. Therefore, although both these studies are listed in our meta-analysis, their results are not used simultaneously in the same outcome analysis. We extracted data on GFR and albuminuria from the recent publication. Data on HbA1c was extracted from the previous publication since it was not reported in the recent one. We believe this strategy will not interfere with excluding the duplicate publication bias as warned in section 10.2.2.1. of the Cochrane Handbook for Systematic Reviews of Interventions V.5.1.0.18

Another difference between our study and Pan et al’s meta-analysis is that we added two newly conducted RCTs by Koya et al15 and Velázquez et al.31 However, the level of 24 h UUN was 1 in the study by Koya et al, which was no less than that in the control group (APIR=1). In addition, the intervention period of the study by Velázquez et al was only 4 months, which might be insufficient to detect the change in GFR as discussed by Zeller et al.40 We need more large-scale RCTs of sufficient length and sufficiently compliant for more conclusive evidence regarding the effect of LPD on GFR.

Messages for clinicians

This meta-analysis showed that LPD improved the kidney function of patients with diabetic nephropathy only when their diet compliance was fair. This finding lets clinicians reaffirm the importance of long-term sustainability of dietary intervention. We think we are quite aware of the importance, but we have also been experiencing the challenge presented by this task. However, this study result showed that the protein restriction need not be as stringent as we have been thinking. The cut-off value of APIR for fair compliance was 0.9 in this study, which was modest compared with the current clinical guideline.11 This finding casts a new light on our management of diet treatment.

Conclusion and future research

A diet intervention by LPD has modest but significant effects on the course of kidney prognosis in patients with diabetic nephropathy, especially when the intervention is sustainable regarding patients’ compliance. This result of meta-analysis questions whether LPD prevents or delays more important clinical outcomes such as kidney failure, initiation of dialysis and death. Further meta-analyses that focus on these outcomes are needed.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

    Files in this Data Supplement:

Footnotes

  • Contributors UN had the idea for the study and designed the method of this meta-analysis including the inclusion and exclusion criteria, conducted data collection and extraction, wrote the first draft of the report and did the statistical analysis with guidance from MS and TM. UN and HK searched the articles and assessed their eligibility. When discrepancies occurred, the eligibility of the articles was discussed by the committee consisting of UN, MS, TM and SU. TM did the major revision and also made comments. All other authors commented on the draft and approved the final version of the manuscript. UN and SU 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.

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

  • Data sharing statement No additional data are available.