Skip to main content

ORIGINAL RESEARCH article

Front. Nutr., 22 February 2024
Sec. Clinical Nutrition
Volume 11 - 2024 | https://doi.org/10.3389/fnut.2024.1335052

Prognostic value of bioelectrical impedance analysis in head and neck cancer patients undergoing radiotherapy: a VALOR® study

Inmaculada Prior-Sánchez1* Aura Dulcinea Herrera-Martínez2,3 María Teresa Zarco-Martín4,5 Rocío Fernández-Jiménez6,7,8* Montserrat Gonzalo-Marín9,10 Araceli Muñoz-Garach5,11,12 Francisco Javier Vilchez-López13,14 Manuel Cayón-Blanco14,15 Ana Villarrubia-Pozo16 Concepción Muñoz-Jiménez2,3 Felisa Pilar Zarco-Rodríguez17 Juana María Rabat-Restrepo18 Luis Miguel Luengo-Pérez19,20 Hatim Boughanem6,21,22 María José Martínez-Ramírez1 Jose Manuel García-Almeida6,7,8,21
  • 1Endocrinology and Nutrition Department, Jaen University Hospital, Jaen, Spain
  • 2Maimonides Institute for Biomedical Research of Cordoba (IMIBIC), Cordoba, Spain
  • 3Endocrinology and Nutrition Service, Reina Sofia University Hospital, Cordoba, Spain
  • 4Department of Endocrinology and Nutrition, San Cecilio University Hospital, Granada, Spain
  • 5Granada Biosanitary Research Institute (Ibs. Granada), Granada, Spain
  • 6Malaga Biomedical Research Institute and BIONAND Platform, Endocrinology and Nutrition Department, Hospital Virgen de la Victoria de Malaga, Malaga, Spain
  • 7Department of Endocrinology and Nutrition, Quironsalud Malaga Hospital, Malaga, Spain
  • 8Department of Medicine and Dermatology, Malaga University, Malaga, Spain
  • 9Endocrinology and Nutrition Department, Malaga Regional University Hospital, Malaga, Spain
  • 10Malaga Biomedical Research Institute and BIONAND Platform, Malaga, Spain
  • 11Department of Endocrinology and Nutrition, Virgen de las Nieves University Hospital, Granada, Spain
  • 12Network Biomedical Research Center Physiopathology of Obesity and Nutrition (CiberOBN), Carlos III Health Institute, Madrid, Spain
  • 13Endocrinology and Nutrition Department, Hospital Universitario Puerta del Mar, Cadiz, Spain
  • 14Biomedical Research and Innovation Institute of Cadiz, Cadiz, Spain
  • 15Endocrinology and Nutrition Department, Hospital Universitario de Jerez de la Frontera, Cadiz, Spain
  • 16Department of Endocrinology and Nutrition, Seville Institute of Biomedicine (IBIS), Virgen del Rocio University Hospital, Seville, Spain
  • 17Department of Endocrinology and Nutrition, Valme University Hospital, Seville, Spain
  • 18Department of Endocrinology and Nutrition, Macarena University Hospital, Seville, Spain
  • 19Department of Endocrinology and Nutrition, Badajoz University Hospital, Seville, Spain
  • 20Department of Biomedical Sciences, Universidad de Extremadura, Badajoz, Spain
  • 21Spanish Biomedical Research Center in Physiopathology of Obesity and Nutrition (CIBERObn), Carlos III Health Institute, Madrid, Spain
  • 22Unidad de Gestión Clinica Medicina Interna, Lipids and Atherosclerosis Unit, Maimonides Institute for Biomedical Research in Córdoba (IMIBIC), Reina Sofia University Hospital, University of Córdoba, Cordoba, Spain

Introduction: Bioelectrical impedance analysis (BIA) serves as a method to estimate body composition. Parameters such as phase angle (PA), standardized phase angle (SPA), body mass cell (BCM), BCM index (BCMI), and fat-free mass (FFM) might significantly impact the prognosis of head and neck cancer (HNC) patients. The present study aimed to investigate whether bioelectrical parameters can be used to predict survival in the HNC population and establish the optimal cutoff points for predictive accuracy.

Methods: A multicenter observational study was performed across 12 tertiary hospitals in Andalusia (a region from the south of Spain). A total of 494 patients diagnosed with HNC between 2020 and 2022 at different stages were included in this study, with a minimum follow-up period of 12 months. The BIA assessment was carried out during the first 2 weeks of radical radiotherapy treatment with chemotherapy or other systemic treatments. A multivariate logistic regression analysis of overall survival, complications, hospital admission, and palliative care and its relationship with BIA nutritional assessment was performed.

Results: Significant prognostic factors identified in the multivariable analysis encompassed phase angle (PA), standardized phase angle (SPA), body cell mass (BCM), and BCM index (BCMI). Lower PA and BCM values were significantly associated with adverse clinical outcomes. A BCM threshold above 17 kg/m2 was the most significant predictor for predicting survival within the overall HNC population. The PA values of <5.1° in male and <4.8° in female patients showed the best predictive potential for mortality. Increased PA (as a continuous variable) demonstrated a significantly reduced risk for mortality (OR, 0.64; 95% CI, 0.43–0.94; p < 0.05) and a decreased likelihood of hospital admission (OR, 0.75; 95% CI, 0.52–1.07; p < 0.05). Higher BCM correlated with a lower risk of mortality (OR, 0.88; 95% CI, 0.80–0.96; p < 0.01) and a diminished probability of hospital admission (OR, 0.91; 95% CI, 0.83–0.99; p < 0.05).

Conclusion: BIA is a crucial tool in the nutritional assessment of HNC patients. BCM and PA are the main bioelectrical parameters used to predict clinical outcomes in this population. Future studies are needed to validate BIA variables in a large cohort to ensure whether early intensification of nutritional treatment would improve survival.

Background

Patients with head and neck cancer (HNC) represent a clinical, nutritional, and metabolically heterogeneous group of tumors. This population have a high prevalence of involuntary weight loss and muscle mass (MM) and a disease-related malnutrition (DRM) rate that can reach 30–55% (1). DRM is a negative prognostic factor that leads to altered body composition, decreased physical and mental function, and impaired clinical outcomes from the disease, such as greater morbidity, poorer tolerance to cancer treatments, increased health costs, and lower quality of life and life expectancy (26). The importance of quantifying MM and detecting DRM has initiated a new consensus for the diagnosis of malnutrition and the search for techniques that allow early detection of changes in body composition (7). Recent scientific studies revealed the need to use other tools for the diagnosis and evaluation of malnutrition in cancer patients. These tools provide a more precise diagnosis and allow categorizing and quantifying the degree of DRM and MM loss. In this context, bioelectrical impedance analysis (BIA) is a simple and an inexpensive technique, which is available in most health centers, that provides information about body composition. By passing an electric current through the body, resistance (R) and reactance (Xc) are measured. BIA data provide information regarding different body compartments, including fat mass (FM), fat-free mass (FFM), and appendicular skeletal muscle mass (ASMM), and hydration levels, encompassing total body water (TBW), extracellular water (ECW), and intracellular water (ICW). All these parameters can help with the nutritional diagnosis. BIA not only estimates body composition and MM but also proves useful for objectifying longitudinal changes in body composition parameters (8, 9). Although BIA equipment does not measure MM directly, it can be employed for the diagnosis of malnutrition, sarcopenia, and cachexia (7, 10, 11). In 2010 (10), and after a new update in 2019 (11), the European Working Group on Sarcopenia in Older People (EWGSOP) established the diagnostic criteria for sarcopenia employed in clinical practice using state-of-the-art techniques, including the BIA. However, these definitions and cutoff points have been developed based on the research conducted in older adults and have not been extensively investigated in cancer populations. The group recognizes that more studies are necessary to validate the measurements in specific subgroups, such as those observed in patients with HNC. The new vector systems also calculate the values of phase angle (PA), standardized phase angle (SPA), and body cell mass (BCM), which are closely related to all metabolically active tissues of the body, including lean muscle mass (8, 12). PA could be considered both a prognostic and a nutritional factor in critical patients (13, 14), in various tumors (1520), and in patients with HNC (2, 2123). However, different PA cutoff points have been proposed. On the other hand, there is a lack of data concerning the predictive role of BCM and other bioelectrical parameters determined by BIA in the HNC population.

This study aimed to explore the prognostic value of BIA parameters in a cohort of HNC patients treated with RT. The secondary objective was to set up the predictive cutoff points of nutritional bioelectrical parameters and their correlation to clinical outcomes in this population.

Materials and methods

Study design

This longitudinal observational study of routine clinical practice was conducted by the “VALOR: VALoracion morfofuncional en el paciente OncoRadioterapico” group at the Departments of Endocrinology and Nutrition in 12 hospitals in Andalusia, Spain [“Virgen de la Victoria” University Hospital and Malaga Regional Hospital (Málaga, Spain), Jaen University Hospital (Jaen, Spain), “San Cecilio” Hospital and “Virgen de las Nieves” University Hospital (Granada, Spain), “Puerta del Mar” University Hospital and Jerez de la Frontera University Hospital (Cádiz, Spain), “Virgen del Rocio” University Hospital, “Virgen Macarena” and “Virgen de Valme” University Hospital (Seville, Spain), “Reina Sofía” University Hospital (Córdoba, Spain), and Badajoz University Hospital (Badajoz, Spain)] between 2020 and 2022. Before beginning the study, the approval of the local ethics committee was obtained (reference code: 2381-M1-22).

Patient selection

A total of 494 patients diagnosed with HNC at different stages were included in the study. The diagnosis was confirmed through medical records and pathological examinations, and the biopsy samples were classified by pathologists according to the histological features following the guidelines outlined in the “World Health Organization Classification of Tumors of the Digestive System” (2016) (24). All patients were included during the first 2 weeks of radical radiotherapy treatment with chemotherapy and other systemic treatments. All patients had at least a 1-year follow-up. This follow-up entailed a clinical visit at 3, 6, and 12 months during the first year. Patients were excluded from the study if they had received more than 2 weeks of radiotherapy at the time of inclusion and if they declined to undergo nutritional measurements by BIA due to reasons related to ethnicity, extensive skin lesions, extravasation of fluids through the route and local hematomas, amputation, or having a life expectancy of less than 3 months. All patients provided written informed consent.

Demographic and anthropometric variables

Demographic and anthropometric variables including age (years), sex (male/female), BMI (kg/m2), weight (kg), and weight loss (%) were collected.

Clinicopathological variables

TNM stage (14), chemotherapy (yes/no), complications (yes/no): dermatitis (yes/no), dysphagia (yes/no), mucositis (yes/no), asthenia (yes/no), unplanned admissions (yes/no; days), exitus (yes/no), and persistence/free disease were included. The Eastern Cooperative Oncology Group (ECOG) scale was used to evaluate the functional status and quality of life (0: asymptomatic, normal activity, 1: symptomatic, may wander, 2: bedridden <50% of the day, with minimal assistance, 3: bedridden >50% of the day, with notable assistance, 4: bedridden all day, severely limited, and 5: deceased).

Biochemical variables

Bimolecular markers such as albumin (g/ dL), prealbumin (mg/dL), proteins (g/dL), total cholesterol (mg/dL), creatinine (mg/dL), glomerular filtration rate (ml/min), urea (mg/dL), glucose (mg/dL), C-reactive protein (CRP, mg/dL), TSH (μUI/mL), and HbA1c (%) were measured.

Body composition measurements

Body weight was determined using the BIA device’s weight scale measured to the nearest 0.1 kg, with patients standing in the center of the platform without shoes and wearing only underwear. Height was measured using a seca stadiometer (Hamburg, Germany). Abdominal circumference was measured while patients lay on their backs at the level of the belly button. Body composition analysis was performed using a 50-kHz phase-sensitive impedance analyzer (BIA 101 Whole Body Bioimpedance Vector Analyzer, AKERN, Florence, Italy) that delivers 800 μA using tetrapolar electrodes positioned on the right hand and right foot. The values of phase angle (PA,°), standardized phase angle (SPA), body cell mass (BCM, kg), body cell mass index (BCMI, kg/m2), fat mass (FM, kg), fat mass index (FMI, kg/m2), fat-free mass index (FFMI, kg/m2), appendicular skeletal muscle mass (ASMM, kg), skeletal muscle index (SMI, kg/m2), total body water (TBW, kg), extracellular water (ECW, kg), intracellular water (ICW, kg), hydration levels (%), reactance (Xc, Ω/m), and resistance (Rz, Ω/m) were obtained. All BIA measurements were obtained with patients in a supine position on a hospital bed. To stabilize BIA values [±2 Ω for resistance (Rz) and ± 1 Ω for Xc (reactance)], patients remained in a supine position for 5 min before the BIA measurements were obtained, as fluid shifts occur after moving from standing to recumbence and directly affect the R and Z values. The BIA measurements of patients were standardized for sex and age using the data collected from healthy Italian adults (25). PA is expressed in degrees as arctan (Xc/Rz) × (180o/π). A standardized PA value (SPA) for each individual was determined from the sex- and age-matched reference population value by subtracting the reference PA value from the observed patient PA value and dividing the result by the respective age- and sex-reference standard deviation (SD). The technical accuracy of the BIA instrument was daily assessed using a precision circuit supplied by the BIA device manufacturer (AKERN, Florence, Italy). All measured Rz and Xc values were consistently ±1 Ω of the 385 Ohm reference value. In vivo reproducibility of the BIA measurements was also determined, with coefficients of variation (CV) of 1–2% for Rz and Xc.

Muscle ultrasonography of the quadriceps rectus femoris (QRF) of the lower extremity was performed using a 10–12 MHz probe and a multifrequency linear matrix (Mindray Z60, Madrid, Spain) in all subjects, who were positioned in a supine position. The evaluation was performed without compression at the level of the lower third from the superior pole of the patella and the anterior superior iliac spine, measuring the anteroposterior muscle thickness, circumference, and cross-sectional area (26). The ultrasonography was performed by a specific physician who was trained in this technique previously. The probe was aligned perpendicular to the longitudinal and transverse axes in the QRF, such as rectus femoris cross-sectional area (RF-CSA), rectus femoris circumference (RF-CIR), RF-axis (X-axis and Y-axis), and L-SAT (subcutaneous fat of the leg). For each parameter, three measurements were performed, and the mean was calculated. For the evaluation of the adipose tissue in the abdominal area, we measured the midpoint between the xiphoid appendix and the navel to capture the image, where total subcutaneous abdominal fat (T-SAT), superficial subcutaneous abdominal fat (S-SAT), and preperitoneal or visceral fat (VAT) data will be measured in centimeters (24). To calculate the global adipose tissue (GAT) and GAT index (GATi), the T-SAT, L-SAT, and VAT were added together, while for GATi, the sum of T-SAT, L-SAT, and VAT was divided by the individual’s height.

Functional assessment

Handgrip strength (HGS) was measured using a Jamar hand dynamometer (Asimow Engineering Co., Los Angeles, CA, USA). Grip strength was calculated in a seated position with the elbow flexed at 90° in the dominant hand. Patients were instructed to perform three maximal isometric contractions with brief pauses between measurements, and the maximum and mean values were recorded. The Timed Up and Go test was selected to evaluate functional capacity. It was performed with patients seated in a chair. The time taken to get up, walk 3 m, turn around, walk another 3 m, and sit back down was measured in seconds.

Assessment of nutritional status and definition of clinical outcomes

To diagnose malnutrition based on the GLIM criteria (7), weight loss, reduced BMI, and reduced muscle mass were categorized as phenotypic criteria, while reduced food intake/ assimilation and disease burden/inflammation were classified as etiologic criteria. For the diagnosis of malnutrition, the combination of at least one phenotypic criterion and one etiological criterion was required. To detect the reduced food intake, a dietary survey was conducted, and a 24-h reminder was collected for each patient. Through an online tool developed by the University of Valladolid,1 the percentage of macronutrientssuch as carbohydrates (g/day), proteins (g/day), lipids (g/day), fiber (g/day), and energy requirements (Kcal/day) were estimated. A low intake was considered when it represented less than 75% of the calculated requirements. Furthermore, in each clinical interview, the presence of gastrointestinal symptoms, such as dysphagia, nausea, vomiting, diarrhea, constipation, and abdominal pain, was investigated. These symptoms have been incorporated as supportive indicators into the GLIM consensus criterion to help identify poor food intake or assimilation. To detect inflammation, supportive proxy measures of inflammation (laboratory indicators) such as serum C-reactive protein (CRP), albumin, or prealbumin were also collected. The following diagnostic criteria for moderate and severe cases were used: Moderate cases included a BMI < 20 and age < 70 years, or FFMI <17 in male subjects with weight loss between 5 and 10%, and a BMI < 22 and age ≥ 70 years, or FFMI <15 in female subjects with weight loss between 5 and 10%. Severe cases were diagnosed with a BMI < 18.5 and age < 70 years with weight loss greater than 10%, and a BMI < 20 and age ≥ 70 years with weight loss greater than 10%. The European Working Group on Sarcopenia in Older People’s diagnostic criteria for sarcopenia were used to assess decreased muscle strength, with a cutoff of <27 kg in male patients and < 16 kg in female patients (10). The main clinical outcomes included a mortality rate at 600 days, complications derived by the RT treatment (dermatitis, dysphagia, mucositis, and asthenia), need for hospital admission during or after treatment (yes/no), or requirement of palliative care (yes/no). Palliative nature was defined as a disease in a terminal stage with no possibility of cure with conventional treatments. In these cases, medical care was offered so as to improve the quality of life and/or mitigate symptoms. The final status was classified as free or persistence/progression of disease. Disease-free status was defined as a remission of the tumor after the initial therapy without evidence of lesions in the imaging tools or clinical examination. Persistence of disease was used to identify patients whose cancer remained stable after initial therapy in the imaging tests and clinical examination. Progression of disease was defined as worsening of the tumor (local or distant), as evidenced by imaging scans or clinical examination. The disease status was established by the responsible physician. The clinical outcomes were extracted from the hospital’s medical record.

Statistical analysis

The results are presented as mean ± standard deviation (SD) for continuous variables and as numerical values (percentages) for categorical variables. Statistical tests such as the Student t-test or Wilcoxon test were employed based on the normality of the variables included. Pearson’s correlation coefficients among variables were calculated, and both linear and logistic regression analyses were conducted. A logistic regression analysis provided the odds ratio (OR) with 95% confidence intervals (CIs). The predictive capability of muscle mass variables was assessed using the receiver operating characteristic (ROC) curves and the area under the curve (AUC). Cox regression was used to analyze the association between mortality and body composition variables. Multivariate logistic regression of complications, hospital admission and palliative care and its relationship with nutritional assessment methods were performed. The results were adjusted by age, gender, BMI, hsCRP, and ECOG: [0 (comprising 217 patients) and a sum ranging from 1 to 5 (encompassing 239 patients)] (*p < 0.05; **p < 0.01; ***p < 0.001). A nomogram model predicting mortality risk was formulated and validated via multivariable COX regression analysis, utilizing the rms and foreign packages in R software. Decision tree analysis utilized the rpart package, and random forest analysis was conducted using the Randomforest package. All analyses and graphical representations were performed using R v.3.5.1 software (Integrated Development for R, RStudio, PBC, Boston, MA, USA). Statistical significance was set at a p-value of <0.05 or through the chi-squared test, when applicable.

Results

Characteristics of the patient

A total of 494 patients (386 men, 108 women) with a histologically confirmed HNC were included in the study. The median age of the patients was 63.9 years. Table 1 describes the baseline characteristics of the population. Demographic variables, BIA parameters, echography exploration, functional measurement, biochemical parameters, clinicopathological variables, and relevant clinical outcomes are shown in Table 1. Approximately 55% of the patients were in an advanced tumor stage (III:19.1% or IV:36.6%) at the time of diagnosis. Approximately half the patients had a normal performance status (ECOG 0), while 8.5% of the patients presented significant limitations in their health condition (ECOG≥2). The mean weight loss 3–6 months before diagnosis was 6.25%. A total of 54.6% received adjuvant chemotherapy. A total of 41% presented disease progression during the follow-up. Only 27% reached free disease status after initial treatment. Approximately 51.9% of the patients needed unintentional hospitalization, and in 63.3% of the patients, complications were revealed. Mortality at 600 days was confirmed to be at 24.9% (Table 1).

Table 1
www.frontiersin.org

Table 1. Baseline characteristics of the population of study.

Body composition measured by BIA and principal HNC outcomes

To understand the relationship between body composition and clinical outcome, Pearson’s correlation coefficient was used. The muscle mass determined by BIA and PA showed a negative correlation with the main clinical results (complications, palliative status, hospital admission, and overall survival, see Figure 1). This negative correlation was more prominent in male participants than in female participants, as can be seen comparatively in Supplementary Figure S1 (male participants) and Supplementary Figure S2 (female participants).

Figure 1
www.frontiersin.org

Figure 1. Correlation plots are presented to show the association between body composition (X-axis) and cancer complications related to head and neck cancer (Y-axis) of all participants. Pearson’s correlation coefficient was used, and an asterisk indicates a significant correlation between variables according to Pearson’s correlation test (*p < 0.05; **p < 0.01 and ***p < 0.001). BCM, body cell mass; BCMI, BCM index; FFMI, fat-free mass index; FM, fat mass; FMI, FM index; OR, odds ratio; PA, phase angle; SMI, skeletal muscle index; SPA, standardized PA.

The study presented the results of multivariate logistic regression analyzing the overall survival, complications, hospital admission, and palliative care in relation to BIA nutritional assessment, as detailed in Table 2. The outcomes were adjusted for age, sex, and BMI. PA, BCM, and BCMI emerged as significant prognostic factors for the overall survival. Elevated PA (as a continuous variable) exhibited a notably reduced risk for mortality (OR, 0.64; 95% CI, 0.43–0.94; p < 0.05), along with a decreased risk for hospital admission (OR, 0.75; 95% CI, 0.52–1.07; p < 0.05). High BCM was associated with a lower risk of mortality (OR, 0.88; 95% CI, 0.80–0.96; p < 0.01) and a reduced risk for hospital admission (OR, 0.91; 95% CI, 0.83–0.99; p < 0.05). On the other hand, SPA only revealed a significant low reduction for hospital admission (OR, 0.75; CI, 0.56–0.99; p < 0.05). Additionally, BCMI was identified as a significant prognostic factor for survival but not for other clinical outcomes. No statistically significant differences were observed among the remaining clinical variables analyzed (Table 2).

Table 2
www.frontiersin.org

Table 2. Multivariate logistic regression of nutritional assessment methods and principal head and neck cancer outcomes.

Decision tree and importance of BIA assessment to predict the risk for mortality

To assess the significance of BIA variables in predicting mortality, a random forest analysis was conducted (Figure 2A). The findings from this model underscore BCM as the most important predictor of survival and a determinant risk factor for mortality in HNC. Specifically, a BCM value above 17 kg correlates with a 91% chance of survival within the overall HNC population (p = 0.002). The accuracy of the model is shown in Figure 2B. Notably, it was confirmed that there is no evidence of multicollinearity within the model. Figure 2C illustrates the individual contribution of each variable included in the analysis. Among the variables, BCM, BCMI, FFMI, and PA demonstrate the highest discriminatory power for categorizing patients as either alive or deceased within the overall sample. Sex-based disparities are shown in Supplementary Figure S3 (male patients) and Supplementary Figure S4 (female patients). Furthermore, distinct random forest and decision tree models were developed based on sex. Figures 2D,E shows the predictive capacity of BCM on survival among HNC patients. Additionally, an ROC curve depicting the adjusted predictive ability of BCM for survival yielded an AUC of 0.73 (Figure 2D).

Figure 2
www.frontiersin.org

Figure 2. Random forest and decision tree of variables to predict mortality. (A) Decision tree performed with the most important variable in the model. Numbers inside the boxes indicate as follows: the left number represents the count of patients predicted as positive by the model; the right number represents the count of patients predicted as negative. The percentage indicates the proportion of individuals included in the model (B). Table to calculate the accuracy of the model (C). Quantitative contribution of each variable included in the analysis (D). The ROC curve of BCM and adjusted BCM to predict survival (E). The predictive model of BCM on survival in HNC patients. AUC, area under the curve; BCM, body cell mass; BCMI, BCM index; FFMI, fat-free mass index; FM, fat mass; FMI, FM index; OR, odds ratio; PA, phase angle; SMI, skeletal muscle index; SPA, standardized PA.

Predictive value and survival outcome

The predictive value of nutritional BIA assessment on survival in HNC patients divided by sex is shown in Table 3. The cut-point PA value that gave the most accurate prediction of survival was 5.1° (sensitivity 64% and specificity 63%; p < 0.001). This value was lower in female patients (PA = 4.8°; p < 0.01) than in male patients (5.1°, p < 0.001), as shown in Table 3. The cutoff point for BCM was set at 28.6 kg in male patients (sensitivity 46% and specificity 83%; p < 0.001) and 17 kg in female patients (sensitivity 86% and specificity 57%; p < 0.01; Table 3). Additionally, Figure 3 shows the AUC curves of the most prognostic variables adjusted for age, sex, and BMI. Supplementary Figure S5 shows the overall survival probability at 600 days based on significant predictive variable cutoff points. The low BCMI group (≤8 kg/m2) exhibited a significantly lower overall survival rate compared to the normal BCMI group (>8 kg/ m2), with a sensitivity of 77% and a specificity of 50% (p < 0.001). Similarly, patients with PA >5.1° demonstrated a higher survival probability when compared to patients with PA ≤5.1°, achieving 70.3% correct predictions (p < 0.001; Supplementary Figure S5).

Table 3
www.frontiersin.org

Table 3. Predictive value of nutritional assessment methods on survival in male and female patients with head and neck cancer.

Figure 3
www.frontiersin.org

Figure 3. The AUC showing the predictive model of the main BIA variables, adjusted by age, sex, and BMI. (A) AUC body cell mass index. (B) AUC skeletal mass index; (C) AUC fat-free mass index. (D) AUC phase angle. (E): AUC standardized phase angle. (F) AUC fat mass. (G). AUC fat mass index. AUC, area under the curve; BCM, body cell mass; BCMI, BCM index; FFMI, fat-free mass index; FM, fat mass; FMI, FM index; OR, odds ratio; PA, phase angle; SMI, skeletal muscle index; SPA, standardized PA.

Discussion

This study aimed to understand the prognostic significance of nutritional BIA assessment in patients with HNC. The findings from this study indicated that low BCM and PA were correlated with cancer outcomes, including increased risks for complications, hospital admissions, and mortality. Our data suggested that PA values below 5.1° in male patients and 4.8° in female patients served as useful reference prognostic values for Andalusian patients with HNC. Moreover, a BCM cutoff point exceeding 17 kg/m2 emerged as the most robust predictor for survival within the overall HNC population, which is associated with a survival probability rate exceeding 90%.

BIA, an easily accessible and non-invasive method, measures the electrical properties of a patient’s tissues and has been employed extensively for assessing body composition, encompassing measurements of FFM, FM, and TBW. In recent years, BIA has been largely utilized to evaluate PA and SPA, which are identified as prognostic factors for survival in severe diseases (13, 14), head and neck cancer (5, 23, 2733), and other tumors (1520, 34). However, varying values have been reported in the literature. For instance, the Swedish group led by Axelsson (33) proposed a slightly higher PA cutoff value of 5.95° for predicting 5-year survival in 128 subjects with HNC. Similar findings were provided by the German study group (32), concluding that patients with normal PA (>5.0°) exhibited notably better survival than malnourished patients (PA < 5.0°); however, the study’s limited sample size of only 42 patients could impact the robustness of conclusions.

Other studies suggested lower PA thresholds as significant markers of poorer prognosis in HNC. A study by Daniel Sat-Muñoz (21) validated that patients with PA <4.42° had reduced survival (19.8 months versus 34.4 months for those with PA >4.42°) and displayed compromised anthropometric, nutritional, and inflammatory status. Similar PA cutoff values have been reported in subsequent studies (27, 29, 35). Furthermore, the prospective randomized HEADNUT trial (29) revealed that pretherapeutic PA (4.7) served as a prognostic indicator for overall survival, while post-therapeutic PA did not influence survival probability. Additionally, a study conducted by Marie Lundberg et al. (35) demonstrated that low PA (<4.5°) correlated with prolonged hospital stays (19.1 versus 8.8 days in the normal PA group).

In our study, not only PA but also other bioelectrical nutritional parameters such as BCM and BCMI emerged as significant prognostic factors for survival. The ROC curve analysis, AUC calculations, and Kaplan–Meier curves were employed to explore their predictive value. Notably, the random forest analysis revealed BCM as the most critical variable in predicting mortality, marking a notable discovery in our study as it underscores the potential utility of BCM/BCMI as prognostic factors in HNC. Few studies have specifically focused on the association between BCM and survival in HNC patients, although correlations between body composition and prognostics have been more extensively explored in oropharyngeal cancer. For instance, a study by Bing Zhuang et al. (30) found that weight loss before RT and body composition changes during RT did not significantly impact the survival outcomes. However, patients with low appendicular skeletal muscle mass index (ASMI: male <7 kg/m2/female <5.7 kg/m2) before RT exhibited poorer overall survival. Additionally, evidence exists linking PA with body composition in the HNC population. The Polish group led by Teresa Małecka-Massalska (28) reported that PA before RT served as a useful marker for identifying HNC individuals at a high risk for unfavorable changes in body composition. Patients with low PA (4.36°) displayed over 9.3-fold higher odds of BMI reduction below 18.5 kg/m2 and 5.9-fold and 4.2-fold higher odds of lean mass and fat mass reduction post-therapy compared to patients with higher PA values (4.98°). Arman Arab et al. (36) conducted a systematic review and meta-analysis, including 14 studies covering 2,625 participants, to evaluate the PA’s predictive ability concerning survival in cancer patients (5 out of 14 studies involved HNC patients). Their analysis confirmed PA as a significant prognostic tool for predicting oncologic patients’ survival, although PA cutoff points varied from 3.8 to 5.9. Another recent systematic review (18) assessed SPA’s role in nutritional status and clinical outcomes in cancer patients. It concluded that, while there is a growing interest in SPA, definitive conclusions regarding its accuracy in predicting clinical outcomes (including complications and survival rates) cannot be drawn due to the absence of standardized SPA cutoff values. In our report, SPA emerged as a significant factor predicting survival, with a cutoff value of −0.75 in male patients and 0.3 in female patients.

The significance of our study lies in addressing the variability seen in PA and SPA cutoff points proposed in the existing literature. Notably, we not only presented the threshold values for the primary BIA variables but also stratified these specific cutoff points by sex. Another strength of our study involved controlling for various confounding factors through a comprehensive assessment of full body composition. Surprisingly, our findings highlighted BCM as the most robust predictor of survival among patients with HNC tumors, surpassing the predictive power of PA or SPA. BCM represents the metabolically active cell mass, encompassing FFM and immune function, and signifies the cell mass involved in critical physiological processes such as O2 consumption, CO2 production, and energy expenditure. Notably, it appears relatively less influenced by non-nutritional factors (8, 12).

Moreover, the multicenter nature of our study, involving a total of 12 tertiary hospitals and encompassing nearly 500 patients, reinforces the robustness of our findings. However, our study has several limitations. First, the relatively short follow-up time restricted our ability to assess the long-term prognosis. Second, inherent limitations of the BIA method, including potential interference from excess fluid and hydration status, could influence the parameters measured. Physiologically, AF values change depending on various factors such as sex (higher in male patients than in female patients), age (direct ratio relationship), BMI (direct correlation with extreme values, exhibiting an inverse correlation), and race. In our study, we adjusted the results for age, sex, BMI, ECOG, and hsCRP. Additionally, the smaller number of female participants might hinder the generalizability of the cutoff points to the female population, posing a limitation.

Conclusion

In conclusion, low values of BCM and PA at diagnosis were significantly associated with shorter overall survival in patients with HNC. A BMC threshold exceeding 17 kg was associated with a 91% probability of overall survival. Our study also suggested that PA values of <5.1° in male patients and <4.8° in female patients were useful prognostic tools for patients with HNC. To the best of our knowledge, this study is one of the most powerful studies, with a large number of patients, to provide reference values for BIA to predict poor survival in patients with HNC. Therefore, BIA could be considered an essential tool for nutritional assessment in this population that will allow clinicians to implement the appropriate procedures to improve HNC patients’ survival. Further studies will be necessary to demonstrate whether intensification of nutritional treatment in this population would enhance clinical outcomes.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by the Granada Biomedical Research Committee (CEI), Granada, Spain (Reference code: 2381-M1-22). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

IP-S: Writing – original draft, Writing – review & editing, Investigation. ADH-M: Writing – original draft, Writing – review & editing, Investigation. MTZ-M: Writing – review & editing, Investigation. RF-J: Writing – original draft, Writing – review & editing, Data curation, Supervision. MG-M: Investigation, Writing – review & editing. AM-G: Investigation, Writing – review & editing. FJV-L: Investigation, Writing – review & editing. MC-B: Investigation, Writing – review & editing. AV-P: Investigation, Writing – review & editing. CM-J: Investigation, Writing – review & editing. FPZ-R: Investigation, Writing – review & editing. JMR-R: Investigation, Writing – review & editing. LML-P: Investigation, Writing – review & editing. HB: Writing – review & editing. MJM-R: Funding acquisition, Resources, Supervision, Writing – review & editing, Investigation. JMG-A: Funding acquisition, Investigation, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by FRESENEIUS KABI. The APC was funded by FRESENIUS KABI. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Acknowledgments

The authors acknowledge the contributions of IBIMA (Malaga Biomedical Research Institute and BIONAND platform) that aided their efforts.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2024.1335052/full#supplementary-material

Footnotes

References

1. Head and Neck Guideline Steering Committee. Evidence-based practice guidelines for the nutritional management of adult patients with head and neck cancer. Sydney: Clinical Oncological Society of Australia. Available at: https://wiki.cancer.org.au/australi and https://wiki.cancer.org.au/australia/COSA:Head_and_neck_cancer_nutrition_guidelines

Google Scholar

2. Yamanaka, A, Yasui-Yamada, S, Furumoto, T, Kubo, M, Hayashi, H, Kitao, M, et al. (2022). Association of phase angle with muscle function and prognosis in patients with head and neck cancer undergoing chemoradiotherapy. Available at: https://pubmed.ncbi.nlm.nih.gov/36027816/

Google Scholar

3. Ganju, RG, Morse, R, Hoover, A, TenNapel, M, and Lominska, CE. The impact of sarcopenia on tolerance of radiation and outcome in patients with head and neck cancer receiving chemoradiation. Radiother Oncol. (2019) 137:117–24. doi: 10.1016/j.radonc.2019.04.023

PubMed Abstract | Crossref Full Text | Google Scholar

4. Sealy, MJ, Dechaphunkul, T, van der Schans, CP, Krijnen, WP, Roodenburg, JLN, Walker, J, et al. Low muscle mass is associated with early termination of chemotherapy related to toxicity in patients with head and neck cancer. Clin Nutr [Internet]. (2020) 39:501–9. doi: 10.1016/j.clnu.2019.02.029

PubMed Abstract | Crossref Full Text | Google Scholar

5. Bentahila, R, Giraud, P, Decazes, P, Kreps, S, Nay, P, Chatain, A, et al. The impact of sarcopenia on survival and treatment tolerance in patients with head and neck cancer treated with chemoradiotherapy. Cancer Med. (2023) 12:4170–83. doi: 10.1002/cam4.5278

PubMed Abstract | Crossref Full Text | Google Scholar

6. Álvarez-Hernández, J, Planas Vila, M, León-Sanz, M, García de Lorenzo, A, Celaya-Pérez, S, García-Lorda, P, et al. Prevalencia y costes de la malnutrición en pacientes hospitalizados; estudio PREDyCES®. Nutr Hosp. (2012) 27:1049–59. doi: 10.3305/nh.2012.27.4.5986

PubMed Abstract | Crossref Full Text | Google Scholar

7. Cederholm, T, Jensen, GL, Correia, MITD, Gonzalez, MC, Fukushima, R, Higashiguchi, T, et al. GLIM criteria for the diagnosis of malnutrition – a consensus report from the global clinical nutrition community. Clin Nutr. (2019) 38:1–9. doi: 10.1016/j.clnu.2018.08.002

PubMed Abstract | Crossref Full Text | Google Scholar

8. Norman, K, Stobäus, N, Pirlich, M, and Bosy-Westphal, A. Bioelectrical phase angle and impedance vector analysis – clinical relevance and applicability of impedance parameters. Clin Nutr [Internet]. (2012) 31:854–61. doi: 10.1016/j.clnu.2012.05.008

PubMed Abstract | Crossref Full Text | Google Scholar

9. Almeida, JMG, García, CG, Aguilar, IMV, Castañeda, VB, and Guerrero, DB. Morphofunctional assessment of patient’s nutritional status: a global approach. Nutr Hosp. (2021) 38:592–600. doi: 10.20960/nh.03378

Crossref Full Text | Google Scholar

10. Cruz-Jentoft, AJ, Baeyens, JP, Bauer, JM, Boirie, Y, Cederholm, T, Landi, F, et al. Sarcopenia: European consensus on definition and diagnosis: report of the European working group on sarcopenia in older people. Age Ageing. (2010) 39:412–23. doi: 10.1093/ageing/afq034

PubMed Abstract | Crossref Full Text | Google Scholar

11. Cruz-Jentoft, AJ, Bahat, G, Bauer, J, Boirie, Y, Bruyère, O, Cederholm, T, et al. Sarcopenia: revised European consensus on definition and diagnosis. Age Ageing. (2019) 48:16–31. doi: 10.1093/ageing/afy169

PubMed Abstract | Crossref Full Text | Google Scholar

12. HPFX, M, ARH, VZ, Kyle, UG, Bosaeus, I, De Lorenzo, AD, Deurenberg, P, et al. PCC of the EWGroup. Bioelectrical impedance analysis—part I: review of principles and methods. Clin Nutr. (2012) 23:1226–43. doi: 10.1016/j.clnu.2004.06.004

Crossref Full Text | Google Scholar

13. Hirose, S, Nakajima, T, Nozawa, N, Katayanagi, S, Ishizaka, H, Mizushima, Y, et al. Phase angle as an Indicator of sarcopenia, malnutrition, and Cachexia in inpatients with cardiovascular diseases. J Clin Med. (2020) 9:1–16. doi: 10.3390/jcm9082554

Crossref Full Text | Google Scholar

14. Do Amaral Paes, TC, KCC, DO, De Carvalho Padilha, P, and WAF, P. Phase angle assessment in critically ill cancer patients: relationship with the nutritional status, prognostic factors and death. J Crit Care. (2018) 44:430–5. doi: 10.1016/j.jcrc.2018.01.006

Crossref Full Text | Google Scholar

15. Souza, NC, Avesani, CM, Prado, CM, Martucci, RB, Rodrigues, VD, de Pinho, NB, et al. Phase angle as a marker for muscle abnormalities and function in patients with colorectal cancer. Clin Nutr. (2021) 40:4799–806. doi: 10.1016/j.clnu.2021.06.013

PubMed Abstract | Crossref Full Text | Google Scholar

16. Yasui-Yamada, S, Oiwa, Y, Saito, Y, Aotani, N, Matsubara, A, Matsuura, S, et al. Impact of phase angle on postoperative prognosis in patients with gastrointestinal and hepatobiliary-pancreatic cancer. Nutrition. (2020) 79-80:110891–80. doi: 10.1016/j.nut.2020.110891

Crossref Full Text | Google Scholar

17. Zhou, S, Yu, Z, Shi, X, Zhao, H, Dai, M, and Chen, W. The relationship between phase angle, nutrition status, and complications in patients with pancreatic head Cancer. Int J Environ Res Public Health. (2022) 19:6526. doi: 10.3390/ijerph19116426

Crossref Full Text | Google Scholar

18. Jiang, N, Zhang, J, Cheng, S, and Liang, B. The role of standardized phase angle in the assessment of nutritional status and clinical outcomes in Cancer patients: a systematic review of the literature. Nutrients. (2023) 15:50. doi: 10.3390/nu15010050

Crossref Full Text | Google Scholar

19. González-Barba, F, Balderas-Peña, LMA, Trujillo-Hernández, B, Cervantes-González, LM, González-Rodríguez, JA, Gutiérrez-Rodríguez, LX, et al. Phase angle and nutritional status: the impact on survival and health-related quality of life in locally advanced uterine cervical Cancer. Healthcare. (2023) 11:246. doi: 10.3390/healthcare11020246

Crossref Full Text | Google Scholar

20. Ji, W, Liu, XL, Zheng, KW, Yang, HM, Cui, JW, and Li, W. Correlation of phase angle with sarcopenia and its diagnostic value in elderly men with cancer. Nutrition. (2021) 84:111110. doi: 10.1016/j.nut.2020.111110

Crossref Full Text | Google Scholar

21. Sat-Muñoz, D, Martínez-Herrera, BE, González-Rodríguez, JA, Gutiérrez-Rodríguez, LX, Trujillo-Hernández, B, Quiroga-Morales, LA, et al. Phase angle, a cornerstone of outcome in head and neck Cancer. Nutrients. (2022) 14:3030. doi: 10.3390/nu14153030

Crossref Full Text | Google Scholar

22. Jung, AR, Roh, JL, Kim, JS, Kim, SB, Choi, SH, Nam, SY, et al. Prognostic value of body composition on recurrence and survival of advanced-stage head and neck cancer. Eur J Cancer. (2019) 116:98–106. doi: 10.1016/j.ejca.2019.05.006

PubMed Abstract | Crossref Full Text | Google Scholar

23. Grossberg, AJ, Chamchod, S, Fuller, CD, Mohamed, ASR, Heukelom, J, Eichelberger, H, et al. Association of Body Composition with Survival and Locoregional Control of radiotherapy-treated head and neck squamous cell carcinoma. JAMA Oncol. (2016) 2:782–9. doi: 10.1001/jamaoncol.2015.6339

PubMed Abstract | Crossref Full Text | Google Scholar

24. Nagtegaal, ID, Odze, RD, Klimstra, D, Paradis, V, Rugge, M, Schirmacher, P, et al. The 2019 WHO classification of tumours of the digestive system. Histopathology. (2020) 76:182–8. doi: 10.1111/his.13975

PubMed Abstract | Crossref Full Text | Google Scholar

25. Piccoli, A, Nigrelli, S, Caberlotto, A, Bottazzo, S, Rossi, B, and Pillon, LMQ. Bivariate normal values of the bioelectrical impedance vector in adult and elderly populations. Am J Clin Nutr. (1995) 61:269–70.

Google Scholar

26. García-Almeida, JM, García-García, C, Vegas-Aguilar, IM, Ballesteros Pomar, MD, Cornejo-Pareja, IM, Fernández Medina, B, et al. Nutritional ultrasound®: conceptualisation, technical considerations and standardisation. Endocrinol Diabetes Nutr. (2023) 70:74–84. doi: 10.1016/j.endinu.2022.03.008

Crossref Full Text | Google Scholar

27. Władysiuk, MS, Mlak, R, Morshed, K, Surtel, W, Brzozowska, A, and Małecka-Massalska, T. Bioelectrical impedance phase angle as a prognostic indicator of survival in head-and-neck cancer. Curr Oncol. (2016) 23:e481–7. doi: 10.3747/co.23.3181

PubMed Abstract | Crossref Full Text | Google Scholar

28. MaŁecka-Massalska, T, Powrózek, T, Prendecka, M, Mlak, R, Sobieszek, G, Brzozowski, W, et al. Phase angle as an objective and predictive factor of radiotherapy-induced changes in body composition of male patients with head and neck cancer. In Vivo (Brooklyn). (2019) 33:1645–51. doi: 10.21873/invivo.11650

PubMed Abstract | Crossref Full Text | Google Scholar

29. Löser, A, Abel, J, Kutz, LM, Krause, L, Finger, A, Greinert, F, et al. Head and neck cancer patients under (chemo-)radiotherapy undergoing nutritional intervention: results from the prospective randomized HEADNUT-trial. Radiother Oncol. (2021) 159:82–90. doi: 10.1016/j.radonc.2021.03.019

Crossref Full Text | Google Scholar

30. Zhuang, B, Xiao, S, Zhang, L, Wang, Y, Zhang, T, Jin, S, et al. Association between body composition and survival in head and neck cancer patients undergoing radiotherapy. Head Neck. (2022) 44:2046–54. doi: 10.1002/hed.27065

PubMed Abstract | Crossref Full Text | Google Scholar

31. van Rijn-Dekker, MI, van den Bosch, L, van den Hoek, JGM, Bijl, HP, van Aken, ESM, van der Hoorn, A, et al. Impact of sarcopenia on survival and late toxicity in head and neck cancer patients treated with radiotherapy. Radiother Oncol. (2020) 147:103–10. doi: 10.1016/j.radonc.2020.03.014

PubMed Abstract | Crossref Full Text | Google Scholar

32. Büntzel, J, Micke, O, Kisters, K, Büntzel, J, and Mücke, R. Malnutrition and survival – bioimpedance data in head neck cancer patients. In Vivo (Brooklyn). (2019) 33:979–82. doi: 10.21873/invivo.11567

PubMed Abstract | Crossref Full Text | Google Scholar

33. Axelsson, L, Silander, E, Bosaeus, I, and Hammerlid, E. Bioelectrical phase angle at diagnosis as a prognostic factor for survival in advanced head and neck cancer. Eur Arch Otorhinolaryngol. (2018) 275:2379–86. doi: 10.1007/s00405-018-5069-2

PubMed Abstract | Crossref Full Text | Google Scholar

34. Paiva, SI, Borges, LR, Halpern-Silveira, D, Assunção, MCF, Barros, AJD, and Gonzalez, MC. Standardized phase angle from bioelectrical impedance analysis as prognostic factor for survival in patients with cancer. Support Care Cancer. (2010) 19:187–92. doi: 10.1007/s00520-009-0798-9

Crossref Full Text | Google Scholar

35. Lundberg, M, Dickinson, A, Nikander, P, Orell, H, and Mäkitie, A. Low-phase angle in body composition measurements correlates with prolonged hospital stay in head and neck cancer patients. Acta Otolaryngol. (2019) 139:383–7. doi: 10.1080/00016489.2019.1566779

PubMed Abstract | Crossref Full Text | Google Scholar

36. Arab, A, Karimi, E, Vingrys, K, and Shirani, F. Is phase angle a valuable prognostic tool in cancer patients’ survival? A systematic review and meta-analysis of available literature. Clin Nutr. (2021) 40:3182–90. doi: 10.1016/j.clnu.2021.01.027

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: bioelectrical impedance analysis (BIA), nutritional assessment, phase angle (PA), body mass cell (BCM), head and neck cancer (HNC)

Citation: Prior-Sánchez I, Herrera-Martínez AD, Zarco-Martín MT, Fernández-Jiménez R, Gonzalo-Marín M, Muñoz-Garach A, Vilchez-López FJ, Cayón-Blanco M, Villarubia-Pozo A, Muñoz-Jiménez C, Zarco-Rodríguez FP, Rabat-Restrepo JM, Luengo-Pérez LM, Boughanem H, Martínez-Ramírez MJ and García-Almeida JM (2024) Prognostic value of bioelectrical impedance analysis in head and neck cancer patients undergoing radiotherapy: a VALOR® study. Front. Nutr. 11:1335052. doi: 10.3389/fnut.2024.1335052

Received: 08 November 2023; Accepted: 23 January 2024;
Published: 22 February 2024.

Edited by:

Arved Weimann, St. Georg Hospital, Germany

Reviewed by:

Nicole Kiss, Deakin University, Australia
Larissa Calixto-Lima, National Cancer Institute (INCA), Brazil

Copyright © 2024 Prior-Sánchez, Herrera-Martínez, Zarco-Martín, Fernández-Jiménez, Gonzalo-Marín, Muñoz-Garach, Vilchez-López, Cayón-Blanco, Villarubia-Pozo, Muñoz-Jiménez, Zarco-Rodríguez, Rabat-Restrepo, Luengo-Pérez, Boughanem, Martínez-Ramírez and García-Almeida. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Rocío Fernández-Jiménez, rociofernandeznutricion@gmail.com; Inmaculada Prior-Sánchez, priorsanchez@hotmail.com

These authors have contributed equally to this work and share first authorship

These authors have contributed equally to this work and share senior authorship

Download