Oxygen-enhanced MRI assessment of tumour hypoxia in head and neck cancer is feasible and well tolerated in the clinical setting.

BACKGROUND: Tumour hypoxia is a recognised cause of radiotherapy treatment resistance in head and neck squamous cell carcinoma (HNSCC). Current positron emission tomography-based hypoxia imaging techniques are not routinely available in many centres. We investigated if an alternative technique called oxygen-enhanced magnetic resonance imaging (OE-MRI) could be performed in HNSCC. METHODS: A volumetric OE-MRI protocol for dynamic T1 relaxation time mapping was implemented on 1.5-T clinical scanners. Participants were scanned breathing room air and during high-flow oxygen administration. Oxygen-induced changes in T1 times (ΔT1) and R2* rates (ΔR2*) were measured in malignant tissue and healthy organs. Unequal variance t-test was used. Patients were surveyed on their experience of the OE-MRI protocol. RESULTS: Fifteen patients with HNSCC (median age 59 years, range 38 to 76) and 10 non-HNSCC subjects (median age 46.5 years, range 32 to 62) were scanned; the OE-MRI acquisition took less than 10 min and was well tolerated. Fifteen histologically confirmed primary tumours and 41 malignant nodal masses were identified. Median (range) of ΔT1 times and hypoxic fraction estimates for primary tumours were -3.5% (-7.0 to -0.3%) and 30.7% (6.5 to 78.6%) respectively. Radiotherapy-responsive and radiotherapy-resistant primary tumours had mean estimated hypoxic fractions of 36.8% (95% confidence interval [CI] 17.4 to 56.2%) and 59.0% (95% CI 44.6 to 73.3%), respectively (p = 0.111). CONCLUSIONS: We present a well-tolerated implementation of dynamic, volumetric OE-MRI of the head and neck region allowing discernment of differing oxygen responses within biopsy-confirmed HNSCC. TRIAL REGISTRATION: ClinicalTrials.gov, NCT04724096 . Registered on 26 January 2021. RELEVANCE STATEMENT: MRI of tumour hypoxia in head and neck cancer using routine clinical equipment is feasible and well tolerated and allows estimates of tumour hypoxic fractions in less than ten minutes. KEY POINTS: • Oxygen-enhanced MRI (OE-MRI) can estimate tumour hypoxic fractions in ten-minute scanning. • OE-MRI may be incorporable into routine clinical tumour imaging. • OE-MRI has the potential to predict outcomes after radiotherapy treatment.

Methodology of generation of CFD meshes and 4D shape reconstruction of coronary arteries from patient-specific dynamic CT.

Due to the difficulties in retrieving both the time-dependent shapes of the vessels and the generation of numerical meshes for such cases, most of the simulations of blood flow in the cardiac arteries use static geometry. The article describes a methodology for generating a sequence of time-dependent 3D shapes based on images of different resolutions and qualities acquired from ECG-gated coronary artery CT angiography. The precision of the shape restoration method has been validated using an independent technique. The original proposed approach also generates for each of the retrieved vessel shapes a numerical mesh of the same topology (connectivity matrix), greatly simplifying the CFD blood flow simulations. This feature is of significant importance in practical CFD simulations, as it gives the possibility of using the mesh-morphing utility, minimizing the computation time and the need of interpolation between boundary meshes at subsequent time instants. The developed technique can be applied to generate numerical meshes in arteries and other organs whose shapes change over time. It is applicable to medical images produced by other than angio-CT modalities.

Feasibility of the J-PET to monitor the range of therapeutic proton beams.

PURPOSE: The aim of this work is to investigate the feasibility of the Jagiellonian Positron Emission Tomography (J-PET) scanner for intra-treatment proton beam range monitoring. METHODS: The Monte Carlo simulation studies with GATE and PET image reconstruction with CASToR were performed in order to compare six J-PET scanner geometries. We simulated proton irradiation of a PMMA phantom with a Single Pencil Beam (SPB) and Spread-Out Bragg Peak (SOBP) of various ranges. The sensitivity and precision of each scanner were calculated, and considering the setup's cost-effectiveness, we indicated potentially optimal geometries for the J-PET scanner prototype dedicated to the proton beam range assessment. RESULTS: The investigations indicate that the double-layer cylindrical and triple-layer double-head configurations are the most promising for clinical application. We found that the scanner sensitivity is of the order of 10-5 coincidences per primary proton, while the precision of the range assessment for both SPB and SOBP irradiation plans was found below 1 mm. Among the scanners with the same number of detector modules, the best results are found for the triple-layer dual-head geometry. The results indicate that the double-layer cylindrical and triple-layer double-head configurations are the most promising for the clinical application, CONCLUSIONS:: We performed simulation studies demonstrating that the feasibility of the J-PET detector for PET-based proton beam therapy range monitoring is possible with reasonable sensitivity and precision enabling its pre-clinical tests in the clinical proton therapy environment. Considering the sensitivity, precision and cost-effectiveness, the double-layer cylindrical and triple-layer dual-head J-PET geometry configurations seem promising for future clinical application.

Radiomic signature accurately predicts the risk of metastatic dissemination in late-stage non-small cell lung cancer.

Background: Non-small cell lung cancer (NSCLC) is the most common type of lung cancer, and the median overall survival (OS) is approximately 2-3 years among patients with stage III disease. Furthermore, it is one of the deadliest types of cancer globally due to non-specific symptoms and the lack of a biomarker for early detection. The most important decision that clinicians need to make after a lung cancer diagnosis is the selection of a treatment schedule. This decision is based on, among others factors, the risk of developing metastasis. Methods: A cohort of 115 NSCLC patients treated using chemotherapy and radiotherapy (RT) with curative intent was retrospectively collated and included patients for whom positron emission tomography/computed tomography (PET/CT) images, acquired before RT, were available. The PET/CT images were used to compute radiomic features extracted from a region of interest (ROI), the primary tumor. Radiomic and clinical features were then classified to stratify the patients into short and long time to metastasis, and regression analysis was used to predict the risk of metastasis. Results: Classification based on binarized metastasis-free survival (MFS) was applied with moderate success. Indeed, an accuracy of 0.73 was obtained for the selection of features based on the Wilcoxon test and logistic regression model. However, the Cox regression model for metastasis risk prediction performed very well, with a concordance index (C-index) score equal to 0.84. Conclusions: It is possible to accurately predict the risk of metastasis in NSCLC patients based on radiomic features. The results demonstrate the potential use of features extracted from cancer imaging in predicting the risk of metastasis.

Detection of range shifts in proton beam therapy using the J-PET scanner: a patient simulation study.

OBJECTIVE: The Jagiellonian PET (J-PET) technology, based on plastic scintillators, has been proposed as a cost effective tool for detecting range deviations during proton therapy. This study investigates the feasibility of using J-PET for range monitoring by means of a detailed Monte Carlo simulation study of 95 patients who underwent proton therapy at the Cyclotron Centre Bronowice (CCB) in Krakow, Poland. Approach: Discrepancies between prescribed and delivered treatments were artificially introduced in the simulations by means of shifts in patient positioning and in the Hounsfield unit to the relative proton stopping power calibration curve. A dual-layer, cylindrical J-PET geometry was simulated in an in-room monitoring scenario and a triple-layer, dual-head geometry in an in-beam protocol. The distribution of range shifts in reconstructed PET activity was visualised in the beam's eye view. Linear prediction models were constructed from all patients in the cohort, using the mean shift in reconstructed PET activity as a predictor of the mean proton range deviation. Main results: Maps of deviations in the range of reconstructed PET distributions showed agreement with those of deviations in dose range in most patients. The linear prediction model showed a good fit, with coefficient of determination r^2 = 0.84 (in-room) and 0.75 (in-beam). Residual standard error was below 1 mm: 0.33 mm (in-room) and 0.23 mm (in-beam). Significance: The precision of the proposed prediction models shows the sensitivity of the proposed J-PET scanners to shifts in proton range for a wide range of clinical treatment plans. Furthermore, it motivates the use of such models as a tool for predicting proton range deviations and opens up new prospects for investigations into the use of intra-treatment PET images for predicting clinical metrics that aid in the assessment of the quality of delivered treatment. .

GPU accelerated Monte Carlo scoring of positron emitting isotopes produced during proton therapy for PET verification.

Objective.Verification of delivered proton therapy treatments is essential for reaping the many benefits of the modality, with the most widely proposedin vivoverification technique being the imaging of positron emitting isotopes generated in the patient during treatment using positron emission tomography (PET). The purpose of this work is to reduce the computational resources and time required for simulation of patient activation during proton therapy using the GPU accelerated Monte Carlo code FRED, and to validate the predicted activity against the widely used Monte Carlo code GATE.Approach.We implement a continuous scoring approach for the production of positron emitting isotopes within FRED version 5.59.9. We simulate treatment plans delivered to 95 head and neck patients at Centrum Cyklotronowe Bronowice using this GPU implementation, and verify the accuracy using the Monte Carlo toolkit GATE version 9.0.Main results.We report an average reduction in computational time by a factor of 50 when using a local system with 2 GPUs as opposed to a large compute cluster utilising between 200 to 700 CPU threads, enabling simulation of patient activity within an average of 2.9 min as opposed to 146 min. All simulated plans are in good agreement across the two Monte Carlo codes. The two codes agree within a maximum of 0.95σon a voxel-by-voxel basis for the prediction of 7 different isotopes across 472 simulated fields delivered to 95 patients, with the average deviation over all fields being 6.4 × 10-3σ.Significance.The implementation of activation calculations in the GPU accelerated Monte Carlo code FRED provides fast and reliable simulation of patient activation following proton therapy, allowing for research and development of clinical applications of range verification for this treatment modality using PET to proceed at a rapid pace.

ProTheRaMon - a GATE simulation framework for proton therapy range monitoring using PET imaging.

OBJECTIVE: This paper reports on the implementation and shows examples of the use of the ProTheRaMon framework for simulating the delivery of proton therapy treatment plans and range monitoring using positron emission tomography (PET). ProTheRaMon offers complete processing of proton therapy treatment plans, patient CT geometries, and intra-treatment PET imaging, taking into account therapy and imaging coordinate systems and activity decay during the PET imaging protocol specific to a given proton therapy facility. We present the ProTheRaMon framework and illustrate its potential use case and data processing steps for a patient treated at the Cyclotron Centre Bronowice (CCB) proton therapy center in Krakow, Poland. APPROACH: The ProTheRaMon framework is based on GATE Monte Carlo software, the CASToR reconstruction package and in-house developed Python and bash scripts. The framework consists of five separated simulation and data processing steps, that can be further optimized according to the user's needs and specific settings of a given proton therapy facility and PET scanner design. MAIN RESULTS: ProTheRaMon is presented using example data from a patient treated at CCB and the J-PET scanner to demonstrate the application of the framework for proton therapy range monitoring. The output of each simulation and data processing stage is described and visualized. SIGNIFICANCE: We demonstrate that the ProTheRaMon simulation platform is a high-performance tool, capable of running on a computational cluster and suitable for multi-parameter studies, with databases consisting of large number of patients, as well as different PET scanner geometries and settings for range monitoring in a clinical environment. Due to its modular structure, the ProTheRaMon framework can be adjusted for different proton therapy centers and/or different PET detector geometries. It is available to the community via github.

The OpenGATE ecosystem for Monte Carlo simulation in medical physics.

This paper reviews the ecosystem of GATE, an open-source Monte Carlo toolkit for medical physics. Based on the shoulders of Geant4, the principal modules (geometry, physics, scorers) are described with brief descriptions of some key concepts (Volume, Actors, Digitizer). The main source code repositories are detailed together with the automated compilation and tests processes (Continuous Integration). We then described how the OpenGATE collaboration managed the collaborative development of about one hundred developers during almost 20 years. The impact of GATE on medical physics and cancer research is then summarized, and examples of a few key applications are given. Finally, future development perspectives are indicated.

Efficiency of 124I radioisotope production from natural and enriched tellurium dioxide using 124Te(p,xn)124I reaction.

BACKGROUND: 124I Iodine (T[Formula: see text] = 4.18 d) is the only long-life positron emitter radioisotope of iodine that may be used for both imaging and therapy as well as for 131I dosimetry. Its physical characteristics permits taking advantages of the higher Positron Emission Tomography (PET) image quality, whereas the availability of new molecules to be targeted with 124I makes it a novel innovative radiotracer probe for a specific molecular targeting. RESULTS: In this study Monte Carlo and SRIM/TRIM modelling was applied to predict the nuclear parameters of the 124I production process in a small medical cyclotron IBA 18/9 Cyclone. The simulation production yields for 124I and the polluting radioisotopes were  calculated for the natural and enriched 124TeO2  +  Al2O3  solid targets irradiated with 14.8 MeV protons. The proton beam was degraded energetically from 18 MeV with 0.2 mm Havar foil. The 124Te(p,xn)124I reactions were taken into account in the simulations. The optimal thickness of the target material was calculated using the SRIM/TRIM and Geant4 codes. The results of the simulations were compared with the experimental data obtained for the natural TeO2 +Al2O3 target. The dry distillation technique of the 124-iodine was applied. CONCLUSIONS: The experimental efficiency for the natural Te target was better than 41% with an average thick target (>0.8 mm) yield of 1.32 MBq/µAh. Joining the Monte Carlo and experimental approaches makes it possible to optimize the methodology for the 124I production from the expensive Te enriched targets.

Non-parametric MRI Brain Atlas for the Polish Population.

Introduction: The application of magnetic resonance imaging (MRI) to acquire detailed descriptions of the brain morphology in vivo is a driving force in brain mapping research. Most atlases are based on parametric statistics, however, the empirical results indicate that the population brain tissue distributions do not exhibit exactly a Gaussian shape. Our aim was to verify the population voxel-wise distribution of three main tissue classes: gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF), and to construct the brain templates for the Polish (Upper Silesian) healthy population with the associated non-parametric tissue probability maps (TPMs) taking into account the sex and age influence. Material and Methods: The voxel-wise distributions of these tissues were analyzed using the Shapiro-Wilk test. The non-parametric atlases were generated from 96 brains of the ethnically homogeneous, neurologically healthy, and radiologically verified group examined in a 3-Tesla MRI system. The standard parametric tissue proportion maps were also calculated for the sake of comparison. The maps were compared using the Wilcoxon signed-rank test and Kolmogorov-Smirnov test. The volumetric results segmented with the parametric and non-parametric templates were also analyzed. Results: The results confirmed that in each brain structure (regardless of the studied sub-population) the data distribution is skewed and apparently not Gaussian. The determined non-parametric and parametric templates were statistically compared, and significant differences were found between the maps obtained using both measures (the maps of GM, WM, and CSF). The impacts of applying the parametric and non-parametric TPMs on the segmentation process were also compared. The GM volumes are significantly greater when using the non-parametric atlas in the segmentation procedure, while the CSF volumes are smaller. Discussion and Conclusion: To determine the population atlases the parametric measures are uncritically and widely used. However, our findings suggest that the mean and parametric measures of such skewed distribution may not be the most appropriate summary statistic to find the best spatial representations of the structures in a standard space. The non-parametric methodology is more relevant and universal than the parametric approach in constructing the MRI brain atlases.

Circadian Profile of Salivary Melatonin Secretion in Hypoxic Ischemic Encephalopathy.

PURPOSE: In the present study, the salivary melatonin secretion in the hypoxic ischemic encephalopathy (HIE) children was measured. The logit model was fitted to the data to obtain the salivary dim light melatonin onsets (DLMOs), and the results were compared with the values estimated from the classic threshold method with a linear interpolation and those previously published for the blood measurements. MATERIALS AND METHODS: 9 patients suffering from HIE aged from 65 to 80 months were included in the study. The melatonin levels were assessed by a radioimmunoassay (RIA). The diurnal melatonin secretion was estimated using a nonlinear least squares method. Student's t-test and the Mann-Whitney U test were used for the comparisons of the obtained parameters. RESULTS: The circadian profiles of the melatonin secretion for both calculation methods do not differ statistically. The DLMO parameters obtained in the blood and saliva samples in children with hypoxic ischemic encephalopathy were similar.

Radiomics and artificial Intelligence for PET imaging analysis.

In recent years, processing of the imaging signal derived from CT, MR or positron emission has proven to be able to predict outcome parameters in cancer patients. The processing techniques of the signal constitute the discipline of radiomics. The quantitative analysis of medical images outperform the information that can be obtained through traditional visual analysis. The recognition of neoplasm molecular and genetic characteristics in a non-invasive way, based on routine radiological examinations, potentially allow complete tumor profiling and subsequent treatment customization at practically zero costs. This process is further boosted with the availability of increased computing power and development of artificial intelligence approaches.

Whole Brain and Cranial Size Adjustments in Volumetric Brain Analyses of Sex- and Age-Related Trends.

Our goal was to determine the influence of sex, age and the head/brain size on the compartmental brain volumes in the radiologically verified healthy population (96 subjects; 54 women and 42 men) from the Upper Silesia region in Poland. The MRI examinations were done using 3T Philips Achieva with the same T1-weighted and T2-weighted protocols. The image segmentation procedures were performed with SPM (Statistical Parameter Mapping) and FSL-FIRST software. The volumes of 14 subcortical structures for the left and right hemispheres and 4 overall volumes were calculated. The General Linear Models (GLM) analysis was used with and without the Total Brain Volume (TBV) and Intracranial Volume (ICV) parameters as the covariates to study the regional vs. global brain atrophy. After the ICV/TBV adjustments, the majority of sex differences in the specific volumes of interest (VOIs) revealed to be linked to the difference in the head/brain size parameters. The analysis also confirmed the significant effect of the aging process on the brain loss. After the TBV adjustment, the age- and sex-related volumetric trends for the gray and white matter volumes were observed: the negative age dependence of the gray matter volume is more pronounced in the males, while in case of the white matter the positive age-related trend in the female group is weaker. The local losses of the left caudate nucleus and the right thalamus are more advanced than the global brain atrophy. Different head-size correction strategies are not interchangeable and may yield various volumetric results, but when used together, facilitate studies on the regional dependencies inherent to a healthy, but aging, brain.

OpenDose: Open-Access Resource for Nuclear Medicine Dosimetry.

Radiopharmaceutical dosimetry depends on the localization in space and time of radioactive sources and requires the estimation of the amount of energy emitted by the sources deposited within targets. In particular, when computing resources are not accessible, this task can be performed using precomputed tables of specific absorbed fractions (SAFs) or S values based on dosimetric models. The aim of the OpenDose collaboration is to generate and make freely available a range of dosimetric data and tools. Methods: OpenDose brings together resources and expertise from 18 international teams to produce and compare traceable dosimetric data using 6 of the most popular Monte Carlo codes in radiation transport (EGSnrc/EGS++, FLUKA, GATE, Geant4, MCNP/MCNPX, and PENELOPE). SAFs are uploaded, together with their associated statistical uncertainties, in a relational database. S values are then calculated from monoenergetic SAFs on the basis of the radioisotope decay data presented in International Commission on Radiological Protection Publication 107. Results: The OpenDose collaboration produced SAFs for all source region and target combinations of the 2 International Commission on Radiological Protection Publication 110 adult reference models. SAFs computed from the different Monte Carlo codes were in good agreement at all energies, with SDs below individual statistical uncertainties. Calculated S values were in good agreement with OLINDA/EXM 2.0 (commercial) and IDAC-Dose 2.1 (free) software. A dedicated website (www.opendose.org) has been developed to provide easy and open access to all data. Conclusion: The OpenDose website allows the display and downloading of SAFs and the corresponding S values for 1,252 radionuclides. The OpenDose collaboration, open to new research teams, will extend data production to other dosimetric models and implement new free features, such as online dosimetric tools and patient-specific absorbed dose calculation software, together with educational resources.

Mathematical evaluation of melatonin secretion in hypoxic ischemic encephalopathy.

OBJECTIVES: The aims of the study were to assess the kinematics of the lower limbs and pelvis during normal walking in professional ballet dancers and to investigate relationships between movements of segments of the lower limbs and pelvis. METHODS: Thirty one professional ballet dancers and twenty eight controls completed five walking trials at their preferred speed. Kinematic data in the basic anatomical planes for ankle, knee, and hip joints as well as for the pelvis were collected with an optoelectronic motion system. RESULTS: The female ballet dancers had in comparison with the controls significantly larger (p < 0.01) knee flexion in the swing phase and hip abduction in the preswing phase. Compared to the control group, the male ballet dancers had significantly larger dorsiflexion in the final stance and the total pelvic tilt range of motion. The number of significant correlations between kinematic parameters was higher in the female ballet dancers. CONCLUSIONS: It can be concluded that specific movement techniques and compensatory strategies used in ballet dance can alter relationships between movements of segments of the lower limbs during normal walking. The relationships between movements in the joints of the lower limbs and pelvis are stronger in women.

SPECT-CT fusion: a new diagnostic tool for endocrinology.

INTRODUCTION: The introduction of the DICOM format in all diagnostic imaging devices allowed co-registering SPECT, CT, MR and other types of biomedical imaging. Fusion can be performed by dedicated hybrid devices or by means of software. The fusion algorithm consists of two steps: coregistration and simultaneous visualization. Our center gradually implemented SPECT-CT fusion in clinical diagnostic work-up of several endocrinologic and oncologic diseases more than 2 years ago. MATERIAL AND METHODS: An easy and fast algorithm in terms of computational complexity of image fusion was presented and applied to 81 consecutive cases. Thirty-two patients were scheduled to SPECT-CT fusion after thyroidectomy and (131)I treatment for thyroid cancer, twelve after somatostatin receptor scintigraphy, seven after 131I MIBG therapy, six after diagnostic MIBG scintigraphy with (123)I or (131)I, three after parathyroid scintigraphy and two after bone scan. The most common indication to the fusion was the need of metabolic characterization of suspected lesions detected on CT scan. The anatomic localization of a focal uptake seen on SPECT and the evaluation of the radiometabolic therapy effect folloved. RESULTS: A variance of error level observed was a result of human factor, decision on marker's placement, respiratory movements and marker's displacement between acquisitions. However, 74% of patients in our series have fusion results classified as "very good" or "good". CONCLUSIONS: The selection of patients, the training of the personel and the cooperation with radiologists are the most important factors for a correct application and interpretation of the SPECT-CT image fusion.