Impact of Collared Peccaries Dycotiles tajacu (Artiodactyla: Tayassuidae) on understory vegetation in the tropical rainforest of the Nogal-La Selva Biological Corridor, Costa Rica / Impacto del pecarí de collar, Dycotiles tajacu (Artiodactyla: Tayassuidae) en la vegetación del sotobosque del bosque tropical húmedo del Corredor Biológico Local Nogal-La Selva, Costa Rica

Abstract Introduction: Evidence suggests that herbivores, such as peccaries, shape vegetation structure and diversity through predation, trampling, dispersal, and rooting behavior. Objective: To evaluate the impact of peccaries (Dycotiles tajacu) on the understory vegetation of the tropical rainforest in the Nogal-La Selva Local Biological Corridor, Costa Rica, comparing a site with the absence of peccaries to another with the presence of these animals. Methodology: From June to November 2021, 20 experimental exclusions and 20 free access plots, each measuring 2 m2 were used to quantify herbivory, the number of leaf blades, damaged leaves, healthy leaves, sapling height, and fallen biomass at both sites. Results: A higher sapling density was found in the Nogal Reserve, but a lower sapling diversity, while in La Selva there was a higher sapling diversity, but a lower density of seedlings. Herbivory and sapling height in La Selva exceeded those in Nogal. The exclusion of peccaries reduced seedling damage but did not affect the dynamics of fallen biomass. Conclusion: For the design, implementation, and evaluation of the effectiveness of biological corridors, it is crucial to consider plant-animal interactions to enhance the flow of ecological processes through functional and structural connectivity, analyzed from interactions such as those presented in this paper.

Resumen Introducción: Existe evidencia que herbívoros, como los saínos, dan forma a la estructura y diversidad de la vegetación a través del comportamiento de depredación, pisoteo, dispersión y enraizamiento. Objetivo: Evaluar el impacto de los saínos (Dycotiles tajacu) en la vegetación del sotobosque del bosque tropical húmedo en el Corredor Biológico Local Nogal-La Selva, Costa Rica, en un sitio con ausencia y en otro con presencia de saínos. Métodos: De junio a noviembre de 2021 se utilizaron 20 exclusiones experimentales y 20 parcelas de acceso libre de 2 m2, se cuantifico la herbivoría, número de láminas foliares, hojas dañadas, hojas sanas, altura de brinzales y biomasa caída en ambos sitios. Resultados: Se encontró una mayor densidad de brinzales en Reserva Nogal pero una menor diversidad, contrario en La Selva donde se encontró una mayor diversidad de brinzales, pero una menor densidad de plántulas. La herbivoría y la altura de brinzales en La Selva fue mayor que en Nogal. La exclusión de los saínos disminuyó el daño a las plántulas, pero no afectó la dinámica de la biomasa caída. Conclusión: Es necesario contemplar para el diseño, implementación y evaluación de la efectividad de corredores biológicos, las interacciones planta-animal, para potencializar el flujo de procesos ecológicos mediante la conectividad funcional y estructural, analizada a partir de interacciones como las presentadas en este trabajo.

Intrapatient variation in PD-L1 expression and tumor mutational burden and the impact on outcomes to immune checkpoint inhibitor therapy in patients with non-small cell lung cancer.

BACKGROUND: PD-L1 tumor proportion score (TPS) and tumor mutational burden (TMB) are key predictive biomarkers for immune checkpoint inhibitors (ICI) efficacy in non-small cell lung cancer (NSCLC). Data on their variation across multiple samples are limited. METHODS: Patients with NSCLC and multiple PD-L1 TPS and/or TMB assessments were included. Clinicopathologic and genomic data were analyzed according to PD-L1 and TMB variation. RESULTS: In total, 402 PD-L1 sample pairs and 413 TMB sample pairs were included. Concordance between pairs was moderate for PD-L1 (ρ=0.53, P<0.0001) and high for TMB (ρ=0.80, P<0.0001). Shorter time between biopsies correlated with higher concordance in PD-L1, but not in TMB. Major increases (ΔTPS≥+50%) and decreases (ΔTPS≤-50%) in PD-L1 were observed in 9.7% and 8.0% of cases, respectively. PD-L1, but not TMB, decreased with intervening ICI (P=0.02). Acquired copy number loss of CD274, PDCD1LG2, and JAK2 were associated with major decrease in PD-L1 (q<0.05). Among patients with multiple PD-L1 assessments before ICI, cases where all samples had a PD-L1 ≥1%, compared to cases with at least one sample with PD-L1 <1% and another with PD-L1 ≥1%, achieved improved objective response rate and progression-free survival (PFS). Among patients with at least one PD-L1 <1% and one ≥1% before ICI, cases where the most proximal sample was PD-L1 ≥1% had longer median PFS compared to cases where the most proximal PD-L1 was <1%. Among patients with multiple TMB assessments before ICI, patients with a TMB ≥10 mut/Mb based on the most recent assessment, as compared to those with a TMB <10 mut/Mb, achieved improved PFS and OS to ICI; instead, no differences were observed when patients were categorized using the oldest TMB assessment. CONCLUSION: Despite intrapatient concordance in PD-L1 and TMB, variation in these biomarkers can influence ICI outcomes, warranting consideration for reassessment prior to ICI initiation when feasible.

Crossover from gas-like to liquid-like molecular diffusion in a simple supercritical fluid.

According to textbooks, no physical observable can be discerned allowing to distinguish a liquid from a gas beyond the critical point. Yet, several proposals have been put forward challenging this view and various transition boundaries between a gas-like and a liquid-like behaviour, including the so-called Widom and Frenkel lines, and percolation line, have been suggested to delineate the supercritical state space. Here we report observation of a crossover from gas-like (Gaussian) to liquid-like (Lorentzian) self-dynamic structure factor by incoherent quasi-elastic neutron scattering measurements on supercritical fluid methane as a function of pressure, along the 200 K isotherm. The molecular self-diffusion coefficient was derived from the best Gaussian (at low pressures) or Lorentzian (at high pressures) fits to the neutron spectra. The Gaussian-to-Lorentzian crossover is progressive and takes place at about the Widom line intercept (59 bar). At considerably higher pressures, a liquid-like jump diffusion mechanism properly describes the supercritical fluid on both sides of the Frenkel line. The present observation of a gas-like to liquid-like crossover in the self dynamics of a simple supercritical fluid confirms emerging views on the unexpectedly complex physics of the supercritical state, and could have planet-wide implications and possible industrial applications in green chemistry.

Epidemiological and clinical description of patients with oropharyngeal cancer treated in a public oncology referral hospital in Chile.

Introduction: The incidence of squamous carcinoma of the oropharynx (OPSCC) has presented an increase worldwide, a fact that occurs along with a phenomenon of epidemiological transition, whose pathogenesis is linked to human papilloma virus (HPV) in a significant part of the cases. Published evidence at the Latin American level is scarce. The present study aims to evaluate the epidemiological and clinical characteristics of patients with oropharyngeal cancer treated in a public oncology reference centre in Chile. Methodology: A cross-sectional study was carried out. Patients with histological confirmation of OPSCC aged 18 years or older, referred to the National Cancer Institute of Chile between 2012 and 2023 were included. The association with HPV was determined by immunohistochemistry for p16. Results: 178 patients were analysed, most of them in locoregionally advanced stages involving the palatine tonsil. Seventy-seven percent were male, with a median age of 60 years. Sixty-seven percent of patients were positive for p16, with a progressive increase to 85% in the last 2 years of the study. The p16(+) patients were younger and had fewer classical risk factors. Primary treatment was radiotherapy in 94% of patients. Conclusion: The epidemiological profile of patients with OPSCC treated in a Chilean public oncology referral centre reflects the epidemiological transition observed in developed countries. This change justifies the need to adapt health policies and conduct research that considers the characteristics of this new epidemiological profile.

Tumor mutational burden for the prediction of PD-(L)1 blockade efficacy in cancer: challenges and opportunities.

Tumor mutational burden (TMB) is a biomarker that measures the number of somatic mutations in a tumor's genome. TMB has emerged as a predictor of response to immune checkpoint inhibitors (ICIs) in various cancer types, and several studies have shown that patients with high TMB have better outcomes when treated with programmed death-ligand 1-based therapies. Recently, the Food and Drug Administration has approved TMB as a companion diagnostic for the use of pembrolizumab in solid tumors. However, despite its potential, the use of TMB as a biomarker for immunotherapy efficacy is limited by several factors. Here we review the limitations of TMB in predicting immunotherapy outcomes in patients with cancer and discuss potential strategies to optimize its use in the clinic.

A Monte Carlo study on the secondary neutron generation by oxygen ion beams for radiotherapy and its comparison to lighter ions.

Objective.To study the secondary neutrons generated by primary oxygen beams for cancer treatment and compare the results to those from primary protons, helium, and carbon ions. This information can provide useful insight into the positioning of neutron detectors in phantom for future experimental dose assessments.Approach.Mono-energetic oxygen beams and spread-out Bragg peaks were simulated using the Monte Carlo particle transport codesFLUktuierende KAskade, tool for particle simulation, and Monte Carlo N-Particle, with energies within the therapeutic range. The energy and angular distribution of the secondary neutrons were quantified.Main results.The secondary neutron spectra generated by primary oxygen beams present the same qualitative trend as for other primary ions. The energy distributions resemble continuous spectra with one peak in the thermal/epithermal region, and one other peak in the fast/relativistic region, with the most probable energy ranging from 94 up to 277 MeV and maximum energies exceeding 500 MeV. The angular distribution of the secondary neutrons is mainly downstream-directed for the fast/relativistic energies, whereas the thermal/epithermal neutrons present a more isotropic propagation. When comparing the four different primary ions, there is a significant increase in the most probable energy as well as the number of secondary neutrons per primary particle when increasing the mass of the primaries.Significance.Most previous studies have only presented results of secondary neutrons generated by primary proton beams. In this work, secondary neutrons generated by primary oxygen beams are presented, and the obtained energy and angular spectra are added as supplementary material. Furthermore, a comparison of the secondary neutron generation by the different primary ions is given, which can be used as the starting point for future studies on treatment plan comparison and secondary neutron dose optimisation. The distal penumbra after the maximum dose deposition appears to be a suitable location for in-phantom dose assessments.

Recent developments at the Calibration Laboratory for radiation protection instruments and dosemeters of the Paul Scherrer Institute.

The Calibration Laboratory of the Paul Scherrer Institute is responsible for the calibration and verification of radiation protection instruments and dosemeters for internal as well as external customers, such as nuclear power plants and hospitals. The Laboratory is authorised by the Federal Institute of Metrology METAS to perform legal verifications and is accredited as an inspection body, according to ISO 17020, and as a calibration laboratory, according to ISO 17025 by the Swiss Accreditation Service. In addition to routine activities (e.g. calibration of dosemeters, irradiation of passive and active dosemeters on phantoms), the Laboratory is involved in various research and development projects with the goal of supporting the increasingly complex calibration needs and to comply with new regulations, particularly those related to clearance limits in free release measurements. Here we present an overview of the currently on-going projects, which include the investigation of Monte Carlo methods for the calibration of clearance monitors and the development of new calibration procedures for wipe test counters. We also discuss the challenges encountered in these projects, since different stakeholders, ranging from legal authorities to equipment manufacturers and users, are involved in the calibration process.

Evolutionary analysis of JC polyomavirus in Misiones' population yields insight into the population dynamics of the early human dispersal in the Americas.

BACKGROUND: JC polyomavirus (JCV) has an ethno-geographical distribution across human populations. OBJECTIVE: Study the origins of the population of Misiones (Argentina) by using JCV as genetic marker. METHODS: Viral detection and characterization was conducted by PCR amplification and evolutionary analysis of the intergenic region sequences. RESULTS: 22 out of 121 samples were positive for JCV, including 5 viral lineages: MY (n = 8), Eu-a (n = 7), B1-c (n = 4), B1-b (n = 2) and Af2 (n = 1). MY sequences clustered within a branch of Native American origin that diverged from its Asian counterpart about 21,914 years ago (HPD 95% interval 15,383-30,177), followed by a sustained demographic expansion around 5000 years ago. CONCLUSIONS: JCV in Misiones reflects the multiethnic origin of the current population, with an important Amerindian contribution. Analysis of the MY viral lineage shows a pattern consistent with the arrival of early human migrations to the Americas and a population expansion by the pre-Columbian native societies.

Transcriptomic profile comparison of monocytes from rheumatoid arthritis patients in treatment with methotrexate, anti-TNFa, abatacept or tocilizumab.

It is well documented that patients affected by rheumatoid arthritis (RA) have distinct susceptibility to the different biologic DMARDs available on the market, probably because of the many facets of the disease. Monocytes are deeply involved in the pathogenesis of RA and we therefore evaluated and compared the transcriptomic profile of monocytes isolated from patients on treatment with methotrexate alone or in combination with tocilizumab, anti-TNFα or abatacept and from healthy donors. Whole-genome transcriptomics yielded a list of regulated genes by Rank Product statistics and DAVID was then used for functional annotation enrichment analysis. Last, data were validated by qRT-PCR. Abatacept, tocilizumab and anti-TNFa cohorts were separately compared with methotrexate, leading to the identification of 78, 6, and 436 differentially expressed genes, respectively. The upper-most ranked genes were related to inflammatory processes and immune responses. Such an approach draws the genomic profile of monocytes in treated RA patients and lays the basis for finding gene signature for tailored therapeutic choices.

Clonal KEAP1 mutations with loss of heterozygosity share reduced immunotherapy efficacy and low immune cell infiltration in lung adenocarcinoma.

BACKGROUND: KEAP1 mutations have been associated with reduced survival in lung adenocarcinoma (LUAD) patients treated with immune checkpoint inhibitors (ICIs), particularly in the presence of STK11/KRAS alterations. We hypothesized that, beyond co-occurring genomic events, clonality prediction may help identify deleterious KEAP1 mutations and their counterparts with retained sensitivity to ICIs. PATIENTS AND METHODS: Beta-binomial modelling of sequencing read counts was used to infer KEAP1 clonal inactivation by combined somatic mutation and loss of heterozygosity (KEAP1 C-LOH) versus partial inactivation [KEAP1 clonal diploid-subclonal (KEAP1 CD-SC)] in the Memorial Sloan Kettering Cancer Center (MSK) MetTropism cohort (N = 2550). Clonality/LOH prediction was compared to a streamlined clinical classifier that relies on variant allele frequencies (VAFs) and tumor purity (TP) (VAF/TP ratio). The impact of this classification on survival outcomes was tested in two independent cohorts of LUAD patients treated with immunotherapy (MSK/Rome N = 237; DFCI N = 461). Immune-related features were studied by exploiting RNA-sequencing data (TCGA) and multiplexed immunofluorescence (DFCI mIF cohort). RESULTS: Clonality/LOH inference in the MSK MetTropism cohort overlapped with a clinical classification model defined by the VAF/TP ratio. In the ICI-treated MSK/Rome discovery cohort, predicted KEAP1 C-LOH mutations were associated with shorter progression-free survival (PFS) and overall survival (OS) compared to KEAP1 wild-type cases (PFS log-rank P = 0.001; OS log-rank P < 0.001). Similar results were obtained in the DFCI validation cohort (PFS log-rank P = 0.006; OS log-rank P = 0.014). In both cohorts, we did not observe any significant difference in survival outcomes when comparing KEAP1 CD-SC and wild-type tumors. Immune deconvolution and multiplexed immunofluorescence revealed that KEAP1 C-LOH and KEAP1 CD-SC differed for immune-related features. CONCLUSIONS: KEAP1 C-LOH mutations are associated with an immune-excluded phenotype and worse clinical outcomes among advanced LUAD patients treated with ICIs. By contrast, survival outcomes of patients whose tumors harbored KEAP1 CD-SC mutations were similar to those with KEAP1 wild-type LUADs.

Single-cell immune profiling reveals thymus-seeding populations, T cell commitment, and multilineage development in the human thymus.

T cell development in the mouse thymus has been studied extensively, but less is known regarding T cell development in the human thymus. We used a combination of single-cell techniques and functional assays to perform deep immune profiling of human T cell development, focusing on the initial stages of prelineage commitment. We identified three thymus-seeding progenitor populations that also have counterparts in the bone marrow. In addition, we found that the human thymus physiologically supports the development of monocytes, dendritic cells, and NK cells, as well as limited development of B cells. These results are an important step toward monitoring and guiding regenerative therapies in patients after hematopoietic stem cell transplantation.

Urinary volatile Organic compounds in non-alcoholic fatty liver disease (NAFLD), type two diabetes mellitus (T2DM) and NAFLD-T2DM coexistence.

INTRODUCTION: Accumulating evidence have shown a significant correlation between urinary volatile organic compounds (VOCs) profile and the manifestation of several physiological and pathological states, including liver diseases. Previous studies have investigated the urinary metabolic signature as a non-invasive tool for the early discrimination between non-alcoholic fatty liver (NAFL) and non-alcoholic steatohepatitis (NASH), which nowadays represents one of the most important challenges in this context, feasible only by carrying out liver biopsy. OBJECTIVES: The aim of the study was to investigate the differences in the urinary VOCs profiles of non-alcoholic fatty liver disease (NAFLD) patients, diabetes mellitus (T2DM) subjects and NAFLD/T2DM patients. METHODS: Headspace solid-phase microextraction (HS-SPME) coupled with gas chromatography-mass spectrometry (GC-MS) was applied to profile the urinary VOCs. Urine samples were analysed both under acid and alkaline conditions, to obtain a range of urinary volatiles with different physicochemical properties. RESULTS: Urinary VOCs profiles of 13 NAFLD patients, 13 T2DM subjects and 13 NAFLD/T2DM patients were investigated by multivariate and univariate data analysis techniques which allowed to identify 21 volatiles under alkaline conditions able to describe the NAFLD/T2DM group concerning the other two groups. CONCLUSION: Our results suggest that VOCs signatures can improve the knowledge of the pathological condition where NAFLD coexists with T2DM and discovering new features that are not simply the sum of the two diseases. These preliminary findings may be considered as hypothesis-generating, to be clearly confirmed by larger prospective investigations.

Use of cryotherapy to treat obstructing papilloma of an accessory tracheal bronchus: case report.

BACKGROUND: Tracheal papillomatosis is a relatively rare condition with limited data on successful treatment modalities. To our knowledge, this is the first report to describe a papilloma arising from an accessory bronchus. Furthermore, this case report demonstrates successful treatment with clinical and patient-centered improvements after use of Spray Cryotherapy. CASE PRESENTATION: A 71-year-old woman presented with one year history of recurrent fevers and intermittent hemoptysis. Imaging and video bronchoscopy revealed an obstructing papilloma of an accessory tracheal bronchus to the right upper lobe. She was treated with debridement followed by multiple cryotherapy treatments resulting in complete clinical and radiographic resolution of her post-obstructive pneumonia. CONCLUSIONS: This case report not only supports existing literature on the use of cryotherapy for airway diseases but also presents a unique form of obstructing papilloma confined to an accessory bronchus, the only report of its kind based on extensive literature review.

Biological Response Dynamics to Cochlear Implantation: Modeling and Modulating the Electrode-Tissue Interface.

OBJECTIVE: The present study primarily aims to study and model the impedance dynamics following cochlear implant (CI) surgery in humans. Secondarily, to observe how the modulation of a single-dose topic dexamethasone alters this response. DESIGN: CI impedance and impedance subcomponents were measured in a day-by-day basis between CI surgery and its activation ( 1 month). Impedance dynamics were mathematically modeled. Furthermore, we conducted a randomized double-blind placebo-controlled clinical trial to evaluate the effect of dexamethasone application. RESULTS: The model showed good fitting for all impedance's measures and both modulations over time. The clinical impedance and impedance subcomponents showed different patterns between groups during the first week postimplantation. Single-dose topical dexamethasone had a transient effect on Impedances, postponing (for 3 days) but not preventing the rise in impedance and its application was mainly effective at the base of the cochlea. CONCLUSIONS: The proposed mathematical fitting properly resembles the impedance dynamics in humans and their modulation due to the dexamethasone. The impedance subcomponents' dynamic pattern showed a time-course consistent with the biological processes of the foreign body reaction. Single-dose topic dexamethasone has a temporal reduction effect in impedance values on the basal turn of the cochlea but does not prevent the rise on the long-term.

Analytical validity of DecisionDx-SCC, a gene expression profile test to identify risk of metastasis in cutaneous squamous cell carcinoma (SCC) patients.

BACKGROUND: To improve identification of patients with cutaneous squamous cell carcinoma (SCC) at high risk for metastatic disease, the DecisionDx-SCC assay, a prognostic 40-gene expression profile (40-GEP) test, was developed and validated. The 40-GEP assay utilizes RT-PCR gene expression analysis on primary tumor biopsy tissue to evaluate the expression of 34 signature gene targets and 6 normalization genes. The test provides classifications of low risk (Class 1), moderate risk (Class 2A), and high risk (Class 2B) of metastasis within 3 years of diagnosis. The primary objective of this study was to validate the analytical performance of the 40-gene expression signature. METHODS: The repeatability and reproducibility of the 40-GEP test was evaluated by performance of inter-assay, intra-assay, and inter-operator precision experiments along with monitoring the reliability of sample and reagent stability for class call concordance. The technical performance of clinical orders from September 2020 through July 2021 for the 40-GEP test was assessed. RESULTS: Patient hematoxylin and eosin (H&E) stained slides were reviewed by a board-certified pathologist to assess minimum acceptable tumor content. Class specific controls (Class 1 and Class 2B) were evaluated with Levey-Jennings analysis and demonstrated consistent and reproducible results. Inter-assay, inter-operator and intra-assay concordance were all ≥90%, with short-term and long-term RNA stability also meeting minimum concordance requirements. Of the 2586 orders received, 93.5% remained eligible for testing, with 97.1% of all tested samples demonstrating actionable class call results. CONCLUSION: DecisionDx-SCC demonstrates a high degree of analytical precision, yielding high concordance rates across multiple performance experiments, along with exhibiting robust technical reliability on clinical samples.

Effect of UV-B Radiation on Flavonoids and Phenols Accumulation in Tempisque (Sideroxylon capiri Pittier) Callus.

Tempisque (Sideroxylon capiri Pittier) is classified as a threatened species and has been reported with a high content of phenols and flavonoids in the leaves. The use of abiotic elicitors such as radiation has been reported due to the changes it produces in the metabolism of plants by activating their defense mechanisms and increasing the biosynthesis of bioactive compounds with antioxidant capacity such as phenols and flavonoids. Therefore, the aim of this work was to evaluate the effect of UV-B radiation on growth parameters and the synthesis of bioactive compounds in in vitro culture of tempisque callus. For the callus induction, we used thidiazuron (TDZ) and 2,4-dichlorophenoxyacetic acid (2,4-D) at 0, 0.5 and 1 mg/L. Calluses were exposed to UV-B radiation (0, 1, 2, 3 and 4 h/day) for two and four weeks. The highest callus formation index was obtained with TDZ and 2,4-D at 1 mg/mL. The greatest increase in the concentration of phenols and flavonoids was detected in the fourth week with 4 h of exposure per day. The highest concentrations of quercetin (230 µg/g dry weight), kaempferol (235 µg/g dry weight) and gallic acid (240 µg/g dry weight) were found in callus obtained from leaves explants.

Sentinel lymph node biopsy in melanoma: beyond histologic factors.

Sentinel lymph node (SLN) biopsy should be performed with the technical expertise required to correctly identify the sentinel node, in the context of understanding both the likelihood of positivity in a given patient and the prognostic significance of a positive or negative result. National Comprehensive Cancer Network guidelines recommend SLN biopsy for all cutaneous melanoma patients with primary tumor thickness greater than 1 mm and in select patients with thickness between 0.8 and 1 mm, yet admit a lack of consistent clarity in its utility for prognosis and therapeutic value in tumors < 1 mm and leave the decision for undergoing the procedure up to the patient and treating physician. Recent studies have evaluated specific patient populations, tumor histopathologic characteristics, and gene expression profiling and their use in predicting SLN positivity. These data have given insight into improving the physician's ability to potentially predict SLN positivity, shedding light on if and when omission of SLN biopsy in specific patients based on clinicopathological characteristics might be appropriate. This review provides discussion and insight into these recent advancements.

Assessment of secondary neutrons in particle therapy by Monte Carlo simulations.

Objective. The purpose of this study is to estimate the energy and angular distribution of secondary neutrons inside a phantom in hadron therapy, which will support decisions on detector choice and experimental setup design for in-phantom secondary neutron measurements.Approach. Dedicated Monte Carlo simulations were implemented, considering clinically relevant energies of protons, helium and carbon ions. Since scored quantities can vary from different radiation transport models, the codes FLUKA, TOPAS and MCNP were used. The geometry of an active scanning beam delivery system for heavy ion treatment was implemented, and simulations of pristine and spread-out Bragg peaks were carried out. Previous studies, focused on specific ion types or single energies, are qualitatively in agreement with the obtained results.Main results. The secondary neutrons energy distributions present a continuous spectrum with two peaks, one centred on the thermal/epithermal region, and one on the high-energy region, with the most probable energy ranging from 19 up to 240 MeV, depending on the ion type and its initial energy. The simulations show that the secondary neutron energies may exceed 400 MeV and, therefore, suitable neutron detectors for this energy range shall be needed. Additionally, the angular distribution of the low energy neutrons is quite isotropic, whereas the fast/relativistic neutrons are mainly scattered in the down-stream direction.Significance. It would be possible to minimize the influence of the heavy ions when measuring the neutron-generated recoil protons by selecting appropriate measurement positions within the phantom. Although there are discrepancies among the three Monte Carlo codes, the results agree qualitatively and in order of magnitude, being sufficient to support further investigations with the ultimate goal of mapping the secondary neutron doses both in- and out-of-field in hadrontherapy. The obtained secondary neutron spectra are available as supplementary material.