Demographic characteristics and the relationship between PLR and clinical features
A total of 238 participants were selected for the final data analysis, after screening according to the inclusion and exclusion criteria. The complete baseline characteristics of these selected participants were presented according to PLR quartiles (Q): Q1: <103; Q2: ≥103 to <137; Q3: ≥137 to <187; and Q4: ≥187. The median age at the time of operation was 60 years, and about 75.21% of the patients were male. In the pooled analysis, there were no significant differences in baseline characteristics between each quartile, with the exceptions of operative approach and tumor size, type, and grade. The higher PLR group tended to undergo open surgery and showed associations with larger tumor size, NEC type, and higher tumor grade. The details of other baseline features, including sex, age, BMI, smoking status, CgA and Syn staining, and TNM tumor stage are shown in Table 1.
Survival analysis
During follow-up, all-cause mortality was noted in 82 patients. Cox proportional hazards regression analyses were performed to assess the associations between PLR and BMI and OS in patients with g-NEN. Three models were constructed to analyze the independent effects of PLR on survival: the crude model, minimally adjusted model (model I), and fully adjusted model (model II). Tumor size and N stage were elected in the fully adjusted model because of a change of >10% from the initial regression coefficient. This study showed that PLR was negatively associated with the survival of patients with g-NEN. The model-based effect size can be explained as a difference in SD of PLR, which is associated with the mortality risk. As shown in Table 2, the PLR increased per SD, and the mortality risk of patients with g-NEN increased by 67%, 63%, and 54% in the crude (HR=1.67, 95% CI 1.32-2.12, P < 0.001), minimally adjusted (HR=1.63, 95% CI 1.28-2.08, P < 0.001), and fully adjusted models (HR=1.54, 95% CI 1.20-1.98, P = 0.001), respectively. Patients with the highest PLR (Q4) had a 2.55-fold increased risk of all-cause mortality compared to those in the lowest PLR group (Q1) in the fully adjusted model. However, the medium PLR groups (Q2 and Q3) were not significantly different from the lowest PLR group (Q1) in the fully adjusted model (all P values > 0.05). Based on this result, the highest PLR group (Q4) and the remaining three quartiles of PLR (Q1-Q3) were further analyzed, and the highest PLR group had a higher mortality risk in the crude, minimally adjusted, and fully adjusted models than the other three groups (all p values < 0.05). The Kaplan-Meier curves illustrate the cumulative mortality risks presented by PLR quartiles (Figure 2A-B). In contrast, BMI was not found to be associated with the survival of patients with g-NEN as a continuous or categorical variable (Table 2).
Stratified subgroup analysis
A stratified analysis can be performed to examine the primary association of interest at different levels of a potential confounding factor. We employed the stratification variables to investigate the effect sizes of PLR on OS. The results showed that BMI had a significant stratified effect at different levels (Table 3). In non-overweight (BMI < 25 kg/m2) patients with g-NEN, the mortality risk increased by 2.02 times (95% CI 1.46-2.78, P < 0.001) per PLR SD increase. In contrast, preoperative PLR was not closely related to the mortality risk in the subgroup of patients with overweight (BMI ≥25 kg/m2). There was a significant interaction between PLR and BMI (P = 0.006). The other clinical features were well balanced, and no meaningful interactions with PLR were found (all p values > 0.05).
PLR and BMI interaction analysis
Table 4 shows the impact of BMI on the predictive efficacy of PLR for the OS of patients with g-NEN. For those with BMI < 25 kg/m2, PLR was negatively associated with OS, and the mortality risk increased by 2.02 (95% CI 1.46-2.78, P < 0.001), 1.96 (95% CI 1.42-2.72, P < 0.001), and 1.95 (95% CI 1.38-2.76, P < 0.001) times per PLR SD increase in the crude, minimally, and fully adjusted models, respectively. In contrast, the mortality risk increased modestly with PLR in the high BMI subgroup (≥ 25 kg/m2), and this increase was not significant. Tests of the interactions between BMI subgroups and continuous PLR on OS were statistically significant (all P values < 0.05) in the crude, minimally, and fully adjusted models. We performed similar analyses to investigate the impact of PLR categories (Q4 vs Q1-3) on the OS of patients with g-NEN, which resulted in findings that were similar to those in the fully adjusted model (P=0.047). However, the interactions between BMI subgroups and PLR categories were not significant in the crude (P = 0.179) and minimally (P = 0.159) adjusted models (Table 4), which indicated that the interactions between PLR and BMI were more obvious when related variables were adjusted. The PLR groups were analyzed according to BMI (<25 or ≥25 kg/m2) in four cohorts as follows: Group 1, PLR Q1-Q3 and BMI < 25 kg/m2; Group 2, Q1-Q3 and BMI ≥ 25 kg/m2; Group 3, PLR Q4 and BMI < 25 kg/m2; and Group 4, PLR Q4 and BMI ≥ 25 kg/m2. Figure 3A shows the cumulative probabilities of deaths according to these four groups (P=0.016). In group 3, patients with g-NEN with higher PLR and lower BMI presented more dramatic mortality risks than those in the remaining groups. We further compared the cumulative probability of death between group 3 and the remaining three groups (group 1, 2 and 4) (Figure 3B), and the cumulative difference in survival rate was more significant (P=0.002).