
Neuroendocrine tumors (NETs) are slow growing, indolent neoplasms that often have delayed presentation [1]. Distant metastases are common in NETs, developing in 30%–70% patients at the time of presentation, with the liver being the most common site of metastases, accounting for 60%–90% of cases [2]. Presence of neuroendocrine tumors with liver metastases (NETLMs) is an important prognostic indicator of survival regardless of the site of origin. It often results in significant constitutional symptoms [3-5]. The five-year survival of untreated NETLMs is around 30%–40% [5,6]. NETLMs are not infiltrative, but are expansive, pushing the surrounding liver parenchyma. For this reason, survival outcomes after R0 and R1 resections have been shown to be similar to each other. Hence, both R0 and R1 resections are generally considered to be of curative intent [7]. Although curative resection (R0-R1) is ideal, it is only possible in 5%–15% patients as most patients have numerous bilobar metastases that are not amenable to complete resection [8,9]. Several studies have assessed the impact of debulking resections (R2) in patients from whom curative resection is not possible, with some studies also comparing between degrees of debulking defined by the percentage of the tumor that can be resected [9-16]. However, these studies were often based on small sample sizes with inconsistent findings. As such, this systematic review was conducted to evaluate the effect of curative intent surgery (R0-R1) and debulking surgery (R2) on overall survival (OS) in NETLM. This review further examined effects of two hepatic debulking thresholds (≥ 90% and ≥ 70%) in NETLM patients with tumors not amenable to curative resection.
A systematic review of Embase, CINAHL, Medline, Cochrane, and PubMed was undertaken using search a strategy based on a combination of Medical Subject Headings (MeSH) and free-text (neuroendocrine tumor, pancreatic NET, carcinoid, small bowel NET, rectal NET, liver metastasis, neuroendocrine liver metastasis, debulking surgery, cytoreductive surgery, R0, R1, R2 resection, survival) to find studies published before December 2020.
• Participants: Adults (over 18 years) who underwent surgery for resection of NETLM.
• Comparisons: R0-R1 resection vs. debulking surgery (R2 resection) or a comparison between debulking strategies (≥ 90% and ≥ 70%).
• R0 resection: resection with microscopically negative margins.
• R1 resection: resection with microscopically positive margin without any gross residual disease. R0/R1 resections together were considered as resections of curative intent.
• R2 resection: presence of gross residual disease (hepatic or extrahepatic disease) after resection ± ablative therapy.
• Extent of hepatic debulking surgery for NETLM: the extent of debulking was calculated by the surgeon based on imaging studies (pre- and postoperative), intraoperative ultrasonographic assessment, and pathology reports. Patients were classified based on the percentage of gross hepatic disease resected, with common thresholds being ≥ 90% and ≥ 70%.
The primary outcome measure was OS. Secondary outcome measures were other measures of survival and recurrence, including progression-free survival (PFS), disease-free survival (DFS), and disease-specific survival (DSS) if reported.
Three authors (RK, IK, BD) extracted data from the included studies using predefined proformas. The quality of included studies in meta-analysis was assessed using ROBINS-I [17], a tool for assessing the risk of bias. Results are reported using ROBIS tool (Supplementary Table 1, 2). When studies included patients managed by non-surgical treatment, cohorts of patients were excluded from analysis (137 patients; 8%).
All prospective and retrospective studies were considered for the analysis. Abstracts of potentially relevant articles were independently screened by two authors (RK, BD). Full texts of all articles identified as potentially relevant were then reviewed. Reference lists of these studies were also scanned to identify any additional studies not previously identified. When multiple articles from the same group within an overlapping study period were found, only the most recent studies were included to avoid duplication. Any disagreement over the relevance of a study was resolved after discussion. Review articles, editorials, letters/comments, and non-English papers were excluded.
Differences in survival following a curative surgery (R0-R1) and a debulking surgery (R2) were quantified using hazard ratios (HRs) from individual studies. When HRs were not reported, these were estimated from Kaplan-Meier curves using the approach described by Tierney et al. [18]. Numbers at risk were incorporated into this calculation if reported, with constant censoring assumed otherwise. Resulting HRs were then log-transformed and pooled using a random-effects inverse-variance meta-analysis model with Review Manager 5.3 [19]. Survival outcomes following the two hepatic debulking approaches (≥ 90% and ≥ 70%) were also compared in a similar manner if sufficient data were reported. For comparisons when data were inadequate to perform a formal meta-analysis, a descriptive summary of studies was reported instead.
The literature search initially identified a total of 538 studies (Fig. 1), 13 of which met the inclusion criteria. Thus, they were included in the analysis [10,12-16,20-26]. For these 13 included studies, the average patient age was in the range of 51 to 61 years and the majority of primaries were in the small bowel or pancreas. Further details of patient characteristics of included studies are summarized in Table 1 and Supplementary Table 3.
Table 1 . Included studies
Author (year) | Study period | Total (n) | Follow-up time (mon) | R-status | % debulking | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|
Curative (R0-R1) | Debulking (R2) | ≥ 90% | 90%–99% | ≥ 70% | 70%–79% | < 70% | |||||
Graff-Baker et al. (2014) [22] | 2007–2011 | 52 | 38 | 19 (37%) | 33 (63%) | - | 22 (42%) | - | 11 (21%) | - | |
Maxwell et al. (2016) [13] | 1999–2015 | P: 28 SB: 80 | 49 | NR | NR | P: 10 (28%) SB: 32 (32%) | - | P: 18 (50%) SB: 51 (51%) | - | P: 8 (22%) SB: 17 (17%) | |
Woltering et al. (2017) [15] | 2003–2016 | P: 83c) SB: 487c) | Mean: 79 | 99%–100%: P: 36 (43%) SB: 207 (43%) | < 99%: P: 47 (57%) SB: 280 (57%) | P: 15 (18%)d) SB: 113 (23%)d) | P: 32 (39%)d) SB: 48 (10%)d) | P: NRd) SB: 119 (24%)d) | |||
Morgan et al. (2018) [14] | 2006–2016 | 44/42a) | 33 | 24 (55%) | 20 (45%) | - | 12 (27%) | - | 8 (18%) | - | |
Chamberlain et al. (2000) [23] | 1992–1998 | 34b) | 27 | 15 (44%) | 19 (56%) | - | - | - | - | - | |
Ejaz et al. (2018) [24] | 1990–2014 | 612 | 51 | 433 (71%) | 179 (29%) | - | - | - | - | - | |
Elias et al. (2003) [12] | 1985–2000 | 47 | 62 | 37 (79%) | 10 (21%) | - | - | - | - | - | |
Glazer et al. (2010) [10] | 1978–2009 | 140b) | 50 | R0: 117 (84%) | R1/2: 22 (16%) | - | - | - | - | - | |
Scott et al. (2019) [16] | 1999–2017 | 188/184a) | 29 | NR | NR | 54 (31%)e) | - | - | 82 (48%)e) | 36 (21%)e) | |
Nave et al. (2001) [25] | 1983–1996 | 31 | 42 | 10 (32%) | 21 (68%) | - | - | - | - | - | |
Osborne et al. (2006) [26] | 2000–2004 | 61b) | NR | 38 (62%) | 23 (38%) | - | - | - | - | - | |
Que et al. (1995) [20] | 1984–1992 | 74 | 26 | 28 (38%) | 46 (62%) | - | - | - | - | - | |
Wängberg et al. (1996) [21] | NR | 64 | NR | 14 (22%) | 50 (78%) | - | - | - | - | - |
Follow-up times are reported as medians, unless stated otherwise.
NR, not reported; P, pancreatic neuroendocrine tumor; SB, small bowel neuroendocrine tumor.
a)Some patients underwent multiple resections, hence numbers are reported as “no. of resections/patients”. b)The number of patients treated with resection. c)The number of patients for which the volume of tumor resected was recorded. d)Woltering et al. [15] grouped patients as 90%–98% and < 90% for P, and as 90%–98%, 70%–89% and < 70% for SB. e) Scott et al. [16] grouped the patients as > 90%, 70%–90%, and < 70% debulking.
Eleven studies comprising 1,729 patients compared outcomes between resection with curative intent (R0-R1) and debulking (R2) surgery, details of which are reported in Table 2 [10,12,14,15,20-26]. The majority of these studies defined the completeness of a surgery based on the overall R-status, with two studies using liver-specific R-status instead [24,25]. No studies reported a HR with associated 95% confidence interval (CI) from univariable analysis. Hence, these were estimated if possible. Four studies [10,14,22,23] reported no significant differences in OS between curative or debulking surgery on OS. They did not give sufficient data for a HR to be estimated. In three studies [20,24,26], Kaplan–Meier curves comparing OS between R2 vs. R0-R1 resection were reported. Hence, these curves were used to estimate the HR for this comparison. HRs for further three studies were estimated from Kaplan–Meier curves comparing R2 vs. R0 resections [12,21,25]. Woltering et al. [15] reported outcomes for small bowel and pancreatic NETLMs separately. These were treated as two cohorts for analysis [14,15]. In each case, the reference category was 99%–100% debulking (categorized as R0-R1 resection), which was compared to 70%–89% hepatic debulking for small bowel and < 90% debulking for pancreatic NETLMs, respectively.
Table 2 . Overall survival (OS) in curative vs. debulking surgery
Author | Overall | Five-year OS (median) | Univariable analysis of R2 vs. R0-R1 | ||||
---|---|---|---|---|---|---|---|
Extent of debulking (R2) | Curative (R0-R1) | Debulking (R2) | HR (95% CI) | p-value | |||
Chamberlain et al. [23] | N/A | 85% | 63% | - | 0.18 | ||
Ejaz et al. [24] | ≥ 80% (liver-specific) | 85.2% (not reached) | 60.7% (7.3 yr) | 3.63 (2.35–5.62)a) | < 0.001* | ||
Elias et al. [12] | ≥ 97% | R0: 74% R1: 70% | 47% | vs. R0: 2.6 (0.2–28.7)a) vs. R1: 2.5 (0.5–12.5)a) | 0.44b) | ||
Glazer et al. [10] | ≥ 90% | R0: N/A R1: N/A | N/A | - | 0.4 | ||
Morgan et al. [14] | ≥ 90% ≥ 70% | R0: N/A | 90-99%: N/A 70-99%: N/A | - | 0.64 0.45 | ||
Nave et al. [25] | N/A (liver-specific) | R0: 86% | 26% | vs. R0: 3.90 (1.76–8.64)a) | 0.001* | ||
Osborne et al. [26] | ≥ 90% | (Mean: 4.2 yr) | (Mean: 2.7 yr) | 3.10 (0.91–10.52)a) | < 0.01* | ||
Que et al. [20] | N/A | N/A | N/A | 2.73 (0.84–8.93)a) | NS | ||
Wängberg et al. [21] | N/A | R0: 100% | 63% | vs. R0: 3.74 (1.28–10.96)a) | N/A | ||
Woltering et al. [15] (small bowel) | 90%–98% 70%–89% | 95% | 90%–98%: 87% 70%–89%: 89% | vs. 90%–98%: 2.26 (1.29–3.96)a) vs. 70%–89%: 3.27 (2.02–5.29)a) | N/A | ||
Woltering et al. [15] (pancreas) | 90%–98% < 90% | 84% | 90%–98%: 68% < 90%: 56% | vs. 90%–98%: 3.00 (0.57–15.76)a) vs. < 90%: 4.06 (1.11–14.93)a) | N/A | ||
Graff-Baker et al. [22]c) | 90%–99%, 70%–89% | N/Ac) | N/Ac) | - | 0.93c) |
The extent of debulking is based on the overall R-status, except for the stated studies, which used the liver-specific R-status. Survival estimates are reported as rates at five years and/or medians, and are for the combined R0-R1 group and the R2 group, unless stated otherwise. HRs are for debulking (R2) vs. curative (R0-R1), unless stated otherwise.
HR, hazard ratio; CI, confidence interval; N/A, data not reported; NS, non-significant, but p-value was not given.
*Statistically significant at p < 0.05.
a)HRs were not reported, hence were estimated from Kaplan–Meier curves. b)p-value is a comparison of R0 vs. R1 vs. R2. c)Study reported disease-specific survival, rather than overall survival.
As such, a total of eight cohorts from seven studies were included in the meta-analysis of OS by the completeness of surgery (Fig. 2). After pooling these studies, it was found that OS was significantly shorter in debulking (R2) relative to curative intent (R0/R1) resections (p < 0.001), with a pooled HR of 3.49 (95% CI, 2.70–4.51). Effect sizes reported by these studies were similar, with an I2 statistic of 0%. A funnel plot gave no indication of publication bias (Fig. 3). Sensitivity analysis excluding the two studies using a liver-specific rather than overall R-status returned consistent results (HR, 3.28; 95% CI, 2.26–4.77; p < 0.001; I2 = 0%).
Two studies reported multivariable analyses of OS with aim to isolate the independent effect of the extent of resection after accounting for other potentially confounding factors [10,24]. Glazer et al. [10] did not find that the completeness of surgery categorized as R0 vs. R1 vs. R2 resection was a significant predictor of OS in a multivariable model. However, their model used a stepwise approach to variable selection without reporting a p-value or HR. Hence, their finding could not be further interrogated. On the other hand, Ejaz et al. [24] found R2 resection to be an independent predictor of poorer OS after accounting for a range of demographic, tumor-related, and operative factors (HR, 2.92; 95% CI, 1.65–5.17; p < 0.001).
In addition to OS assessment, some studies also reported data for other survival- and recurrence-related outcomes by completeness of resection. Elias et al. [12] found that patients undergoing a curative surgery had significantly longer DFS than those undergoing a debulking surgery (p = 0.003). However, this was not observed by Glazer et al. (p = 0.8) [10], with Graff-Baker et al. [22] reporting no significant differences in PFS (p = 0.38) or DSS (p = 0.93) between curative and debulking groups. Due to small numbers of studies and heterogeneity of reporting, formal meta-analysis of these secondary outcomes was not possible.
Five studies (654 patients) compared outcomes by the degree of hepatic debulking (≥ 90% and ≥ 70%) in the setting of NETLM [13-16,22]. Of these, both Maxwell et al. [13] and Woltering et al. [15] reported outcomes in the small bowel and pancreatic NETLMs separately, giving seven cohorts for analysis (Table 3). However, the majority of these studies did not report sufficient data for HRs or the associated 95% CI to be estimated. As such, formal meta-analysis was not possible for this section of the review. Hence, these studies were instead analyzed using a qualitative approach.
Table 3 . Survival by degree of hepatic debulking
Author | Debulking | Five-year OS (median) | Five-year PFS (median) | ||||||
---|---|---|---|---|---|---|---|---|---|
Group A vs. B | Group A | Group B | p-value | Group A | Group B | p-value | |||
≥ 90% vs. < 90% debulking | |||||||||
Maxwell et al. [13] (small bowel) | ≥ 90% vs. < 90% | Not reached | 9.1 yr | 0.46 | 3.8 yr | 1.6 yr | 0.005* | ||
Maxwell et al. [13] (pancreas) | ≥ 90% vs. < 90% | Not reached | 6.1 yr | 0.14 | 4.4 yr | 1.3 yr | 0.05* | ||
Woltering et al. [15] (small bowel) | 90%–98% vs. 70%–89% | 87% (22.9 yr) | 89% (12.3 yr) | NS | - | - | - | ||
Woltering et al. [15] (pancreas) | 90%–98% vs. < 90% | 68% (6.7 yr) | 56% (6.3 yr) | 0.015* | - | - | - | ||
Morgan et al. [14] | ≥ 90% vs. 70%–89% | NR | NR | 0.29 | NRa) | NRa) | 0.75a) | ||
Graff-Baker et al. [22] | 90%–99% vs. 70%–89% | - | - | - | NRa) | NRa) | 0.74a,b) | ||
Scott et al. [16] | > 90% vs. 70%–90% | Not reached | 11.1 yr | 0.61 | 4.7 yr | 1.7 yr | < 0.01* | ||
≥ 70% vs. < 70% debulking | |||||||||
Maxwell et al. [13] (small bowel) | ≥ 70% vs. < 70% | Not reached | 9.1 yr | 0.18 | 3.2 yr | 1.7 yr | 0.005* | ||
Maxwell et al. [13] (pancreas) | ≥ 70% vs. < 70% | Not reached | 1.7 yr | 0.001* | 3.0 yr | 0.5 yr | < 0.001* | ||
Woltering et al. [15] (small bowel) | 70%–89% vs. < 70% | 89% (12.3 yr) | 64% (7.4 yr) | NR | - | - | - | ||
Scott et al. [16] | ≥ 70% vs. < 70% | 11.2 yr | 3.1 yr | < 0.001* | 1.7 yrc) | 0.9 yr | < 0.001* |
Survival estimates are reported as rates at five years and/or medians, and are for the stated debulking groups.
OS, overall survival; PFS, progression-free survival; NR, data not reported; NS, non-significant, but p-value was not given.
*Statistically significant at p < 0.05.
a)Liver-Specific PFS. b)p-value represents a comparison between three groups: 100% vs. 90%–99% vs. 70%–89%. c)PFS data were reported for the 70%–90% debulking group.
Four studies [13-16] compared OS between debulking strategies. For the threshold of 90% debulking, only Woltering et al. [15] reported a significant difference in OS, with median survival time in the pancreatic NETLM cohort showing a modest reduction from 6.7 to 6.3 years in the 90%–98% debulking group vs. the < 90% debulking group (p = 0.015). Three studies [13,15,16] additionally compared OS after ≥ 70% and < 70% hepatic debulking. Of these, Scott et al. [16] and the pancreatic NETLM cohort of Maxwell et al. [13] both reported significantly shorter OS after < 70% vs. ≥ 70% debulking, with Scott et al. [16] additionally performing a multivariable analysis and finding < 70% debulking to be independently associated with poorer OS. No significant effect was observed in the other two cohorts [13,15]. In addition to assessing debulking using thresholds, Scott et al. [16] also analyzed the degree of debulking as a continuous variable, which was significant (p < 0.01) in a Cox regression model, implying that OS became progressively longer when the proportion of debulking increased.
In addition to the analysis of OS, four studies also assessed the outcome of PFS [13,14,16,22]. Of these, two studies [13,16] reported significantly longer PFS with greater degrees of hepatic debulking, as quantified by both 70% and 90% thresholds. The two other studies [14,22] specifically considered liver-specific PFS. They did not find it to differ significantly between ≥ 90% and < 90% debulking. Scott et al. [16] additionally analyzed the degree of debulking as a continuous variable and found it to be significantly associated with PFS (p < 0.01), with such association persisting on multivariable analysis (p = 0.01).
The optimal management of NETs in the setting of NETLM remains unclear partly due to the heterogeneity of disease behavior and heterogenous reporting of outcomes [27]. Complete surgical resection is the best curative option for NETLM, with reported five- and ten-year OS rates of up to 74% and 51%, respectively, and a median survival three times that of patients with untreated NETLMs. However, curative resection is only possible in 10% to 20% patients. In addition, it is difficult to achieve a curative resection in extensive disease. As such, where curative resection is not feasible, debulking surgery offers an alternative treatment approach as it may reduce the risk of liver failure due to progression of liver disease and provide relief of hormonal symptoms in patients with functional tumors. As such, the present systematic review and meta-analysis compared outcomes between various liver resection strategies for NETLM.
In the first stage of the current analysis, the aim was to assess differences in OS between curative and debulking resections. Pooling of studies identified by the systematic review found that curative surgery (R0/R1) was associated with significantly longer OS than debulking surgery in NETLM, with a pooled HR of 3.49 and consistent effect sizes across studies. The next stage of the analysis aimed to identify effects of different debulking thresholds in those undergoing incomplete resections. Optimal thresholds of debulking suggested in the literature have evolved over time. In 1990, McEntee et al. [28] set the debulking threshold at 90% based on his early experience of 37 patients. Other authors [9,11,20] endorsed this threshold. It became an acceptable oncologic threshold for increasing patient survival. However, more recently, lower thresholds have been proposed, with Chambers et al. [29] reporting a five-year OS of 74% in a cohort with a hepatic debulking threshold of 70%. As such, thresholds of 70% and 90% were used in the current review.
Due to inconsistencies of the extent of debulking used with poor statistical reporting (i.e., absence of HRs and 95% CI) of the identified studies, it was not possible to perform reliable quantitative meta-analysis of these thresholds. Qualitative review of included studies revealed that the majority of them found no significant difference in OS between patients with < 90% debulking and those with ≥ 90% debulking. However, studies comparing < 70% vs. ≥ 70% debulking tended to show that resections below 70% was associated with significantly shorter OS. Analysis of PFS found that the effect of the degree of hepatic debulking was more pronounced for this outcome, with significant differences in PFS consistently being observed for both 90% and 70% thresholds. Scott et al. [16] additionally analyzed the degree of debulking as continuum rather than grouping based on thresholds and found a significant and progressive improvement in both OS and PFS with greater percentage of debulking.
To summarize these findings, the interpretation is that OS becomes progressively shorter as the degree of hepatic debulking decreases. While this would imply that OS will be shorter below the 90% debulking threshold, the magnitude of this difference is insufficient to be clinically (or statistically) relevant. On the other hand, a marked reduction in OS becomes more observable when debulking is below the 70% threshold. With respect to PFS, significant differences are visible even at the 90% debulking threshold, implying that this is insufficient to reduce the risk of recurrence to be in line with that of a curative surgery.
When curative resection of NETLM is not feasible, non-surgical treatment options are an alternative to debulking surgery. Whilst this was not part of the current review, studies have assessed outcomes after trans-arterial therapy, reporting five-year OS of 40% compared to 70% for hepatic resection [26]. The OS was also significantly longer in those who underwent cytoreduction therapy in this study (median, 24 months vs. 43 months). In another propensity score matched study [30], the mean OS was 38 months for the trans-arterial therapy group and 84 months for the surgical group. Yttrium microspheres have been reported to be more promising in long-term disease control of NETLM. A multi-institutional study [31] with 168 patients showed stable disease in 23% and complete response in 3% of patients, with a median OS of 70 months. However, these results have not been reproduced in other studies [32,33].
Another alternative to resection of NETLM is liver transplantation, although this is subject to some debate. Based on European NET guidelines [34] with careful selection of those with young age, stable disease, low Ki67 index, reduced hepatic load, and the absence of extrahepatic disease, studies of liver transplant in NETLM have reported an acceptable five-year OS of over 50% for midgut tumors and up to 50% for pancreatic NETLMs [35,36]. However, the strict selection criteria, the lack of wide acceptance to NETLM as an indication for liver transplant programs, and the limited donor pool remain limiting factors for offering transplantation to this group of patients.
In addition to survival, quality of life is another important outcome to consider when assessing surgical interventions. However, data on quality of life in patients undergoing resection of NETLM are currently sparse. Spolverato et al. [2] reported no difference in the improvement of overall quality of life between surgical and non-surgical groups of patients having an initial treatment for NETLM, although the proportion of patients who reported being dissatisfied with their treatment was significantly lower in the surgical group than in the non-surgical group (5.4% vs. 9.4%; p = 0.001). Patients with a very poor quality of life at the time of the diagnosis were more likely to experience an improvement in quality of life after treatment.
There are several limitations of this review, the majority of which are related to the consistency and quality of reporting of studies identified by the systematic review. The primary limitation was the fact that no study reported HRs with associated 95% CIs for comparisons of interest. As a result, these statistics had to be estimated from Kaplan-Meier curves, which was subject to a margin of error. It might have introduced bias. Another limitation was the fact that studies rarely reported multivariable or adjusted analyses to account for any baseline differences between treatment groups. Many other factors are known to influence survival, including the age of the patient, lymph node metastasis, symptomatic disease, presence of extrahepatic disease, and the site and presence of the primary tumor. The third limitation was the relatively small number of studies identified, particularly for the analysis of debulking thresholds, and small numbers of patients in some of the included studies. Finally, there was some inconsistency of grouping used by studies. For example, some classified curative surgery as a combination of R0 and R1 resections, whilst others reported these as separate groups or reported only R0 resections. This might have resulted in incompatibility of some studies included for comparisons of curative and debulking surgeries as well as comparisons by debulking thresholds. However, when the risk of bias was assessed using ROBIS tool [37], the overall risk was estimated to be low. It is unlikely that randomized trials will be performed to compare outcomes among various treatment modalities. Therefore, future prospective studies on this subject should aim to capture all the above-mentioned patients and tumor-related prognostic factors to allow a uniform and standardized reporting of results. Quality of life and patient-reported outcomes also need to be included in future studies as they have a significant role in selecting long-term treatment options for these patients.
In conclusion, curative intent surgery (R0-R1) is associated with a significantly longer OS than debulking surgery of NETLM. The extent of debulking also appears to influence both OS and PFS, with outcomes being superior with above 70% debulking and a tendency of improved outcomes with above a 90% threshold.
Supplementary data related to this article can be found at https://doi.org/10.14701/ahbps.21-101.
Dr Rebecca Glover, for her initial help with data collection.
None.
No potential conflict of interest relevant to this article was reported.
Conceptualization: BVMD. Data curation: RK, IK, JH. Methodology: JH, BVMD. Visualization: JH, BVMD. Writing - original draft: RK, IK, JH, SR, SP, BVMD. Writing - review & editing: JH, SR, TS, SP, BVMD.
![]() |
![]() |