Search

Objectives The clinical impact of brain microstructural abnormalities in multiple sclerosis (MS) remains elusive. We aimed to characterize the topography of longitudinal relaxation rate (R1) and quantitative susceptibility (χ) changes, as indices of iron and myelin, together with brain atrophy, and to clarify their contribution to cognitive and motor disability in MS. Methods In this cross-sectional study, voxel-based morphometry, and voxel-based quantification analyses of R1 and χ maps were conducted in gray matter (GM) and white matter (WM) of 117 MS patients and 53 healthy controls. Voxel-wise between-group differences were assessed with nonparametric permutation tests, while correlations between MRI metrics and clinical variables (global disability, cognitive and motor performance) were assessed both globally and voxel-wise within clusters emerging from the between-group comparisons. Results MS patients showed widespread R1 decrease associated with more limited modifications of χ, with atrophy mainly involving deep GM, posterior and infratentorial regions (p < 0.02). While R1 and χ showed a parallel reduction in several WM tracts (p < 0.001), reduced GM R1 values (p < 0.001) were associated with decreased thalamic χ (p < 0.001) and small clusters of increased χ in the caudate nucleus and prefrontal cortex (p < 0.02). In addition to the atrophy, χ values in the cingulum and corona radiata correlated with global disability and motor performance, while focal demyelination correlated with cognitive performance (p < 0.04). Conclusions We confirmed the presence of widespread R1 changes, involving both GM and WM, and atrophy in MS, with less extensive modifications of tissue χ. While atrophy and χ changes are related to global and motor disability, R1 changes are meaningful correlates of cognition. Key Points • Compared to healthy controls, multiple sclerosis patients showed R1 and χ changes suggestive of iron increase within the basal ganglia and reduced iron and myelin content within (subnuclei of) the thalamus. • Thalamic volume and χ changes significantly predicted clinical disability, as well as pulvinar R1 and χ changes, independently from atrophy. • Atrophy-independent R1 and χ changes, suggestive of thalamic iron and myelin depletion, may represent a sensitive marker of subclinical inflammation.

The ratio of T1-weighted/T2-weighted magnetic resonance images (T1w/T2w MRI) has been successfully applied at the cortical level since 2011 and is now one of the most used myelin mapping methods. However, no reports have explored the histological validity of T1w/T2w myelin mapping in white matter. Here we compare T1w/T2w with ex vivo postmortem histology and in vivo MRI methods, namely quantitative susceptibility mapping (QSM) and multi-echo T2 myelin water fraction (MWF) mapping techniques. We report a discrepancy between T1w/T2w myelin maps of the human corpus callosum and the histology and analyse the putative causes behind such discrepancy. T1w/T2w does not positively correlate with Luxol Fast Blue (LFB)-Optical Density but shows a weak to moderate, yet significant, negative correlation. On the contrary, MWF is strongly and positively correlated with LFB, whereas T1w/T2w and MWF maps are weakly negatively correlated. The discrepancy between T1w/T2w MRI maps, MWF and histological myelin maps suggests caution in using T1w/T2w as a white matter mapping method at the callosal level. While T1w/T2w imaging may correlate with myelin content at the cortical level, it is not a specific method to map myelin density in white matter.

Although the use of specific MRI criteria has significantly increased the diagnostic accuracy of multiple sclerosis (MS), reaching a correct neuroradiological diagnosis remains a challenging task, and therefore the search for new imaging biomarkers is crucial. This study aims to evaluate the incidence of one of the emerging neuroradiological signs highly suggestive of MS, the central vein sign (CVS), using data from Fabry disease (FD) patients as an index of microvascular disorder that could mimic MS.In this retrospective study, after the application of inclusion and exclusion criteria, MRI scans of 36 FD patients and 73 relapsing-remitting (RR) MS patients were evaluated. Among the RRMS participants, 32 subjects with a disease duration inferior to 5 years (early MS) were also analyzed. For all subjects, a Fazekas score (FS) was recorded, excluding patients with FS = 0. Different neuroradiological signs, including CVS, were evaluated on FLAIR T2-weighted and spoiled gradient recalled echo sequences.Among all the recorded neuroradiological signs, the most striking difference was found for the CVS, with a detectable prevalence of 78.1% (57/73) in RRMS and of 71.4% (25/32) in early MS patients, while this sign was absent in FD (0/36).Our results confirm the high incidence of CVS in MS, also in the early phases of the disease, while it seems to be absent in conditions with a different etiology. These results corroborate the possible role of CVS as a useful neuroradiological sign highly suggestive of MS.• The search for new imaging biomarkers is crucial to achieve a correct neuroradiological diagnosis of MS. • The CVS shows an incidence superior to 70% in MS patients, even in the early phases of the disease, while it appears to be absent in FD. • These findings further corroborate the possible future central role of CVS in distinguishing between MS and its mimickers.

Purpose To introduce a class of fast 3D quantitative MRI (qMRI) schemes (RESUMEN, for N = 1 , … , 4 ) that allow for a thorough characterization of microstructural properties of brain tissues. Methods An arbitrary multi-echo GRE acquisition optimized for quantitative susceptibility mapping (QSM) is complemented with an appropriate low flip-angle GRE sequence drawn from four possible choices. The acquired signals are processed to analytically derive the longitudinal relaxation ( R 1 ) and free induction decay ( R 2 ∗ ) rates, as well as the proton density (PD) and QSM. A comprehensive modeling of the excitation and B 1 - profiles and of the RF-spoiling is included in the acquisition and processing pipeline. Results The RESUMEN maps appear homogeneous throughout the field-of-view and exhibit comparable values and high SNR across the considered range of N values. Conclusions The introduced schemes represent a class of robust and flexible strategies to derive a thorough and fast qMRI study, suitable for a whole-brain acquisition with isotropic voxel resolution of 700 μ m in less than 15 min.

Precision medicine is an innovative and emerging approach to treatment that accounts for individual variability in genetic and environmental factors to identify and utilize the specific biomedical profile of a patient’s disease [...]

To investigate the interplay between spatial dose patterns and single nucleotide polymorphisms in the development of radiation-induced lymphopenia (RIL) in 186 non-small-cell lung cancer (NSCLC) patients undergoing chemo-radiotherapy (RT).This study included NSCLC patients enrolled in a randomized trial of protons vs. photons with available absolute lymphocyte counts at baseline and during RT and XRCC1-rs25487 genotyping data. After masking the GTV, planning CT scans and dose maps were spatially normalized to a common anatomical reference. A Voxel-Based Analysis (VBA) was performed to assess voxel-wise relationships of dosiomic and genomic explanatory variables with RIL. The underlying generalized linear model was designed to include both the explanatory variables (3D dose distributions and the XRCC1-rs25487 genotypes) and possible nuisance variables significantly correlated with RIL. The maps of model coefficients as well as their significance maps were generated.Measures for RIL definition during RT were characterized, including kinetic parameters for lymphocyte loss. The VBA generated three-dimensional maps of correlation between RIL and dose in lymphoid organs as well as organs with abundant blood pools. The identified voxel-wise relationships account for XRCC1-rs25487 polymorphism and demonstrate the variant AA genotype being detrimental to lymphocyte depletion (p = 0.03).The performed analyses blindly highlighted relevant anatomical regions that contributed most to lymphocyte depletion during RT and the interplay of the variant XRCC1-rs25487 AA genotype with the dose delivered to the primary lymphoid organs. These findings may help to guide the development of dosimetric RIL mitigation strategies for the application of effective individualized RT.

Recently, 2D or 3D methods for dose distribution analysis have been proposed as evolutions of the Dose Volume Histogram (DVH) approaches. Those methods, collectively referred to as pixel- or voxel-based (VB) methods, evaluate local dose response patterns and go beyond the organ-based philosophy of Normal Tissue Complication Probability (NTCP) modelling. VB methods have been introduced in the context of radiation oncology in the very last years following the virtuous example of neuroimaging experience. In radiation oncology setting, dose mapping is a suitable scheme to compare spatial patterns of local dose distributions between patients who develop toxicity and who do not. In this critical review, we present the methods that include spatial dose distribution information for evaluating different toxicity endpoints after radiation therapy. The review addresses two main topics. First, the critical aspects in dose map building, namely the spatial normalization of the dose distributions from different patients. Then, the issues related to the actual dose map comparison, i.e. the viable options for a robust VB statistical analysis and the potential pitfalls related to the adopted solutions. To elucidate the different theoretical and technical issues, the covered topics are illustrated in relation to practical applications found in the existing literature. We conclude the overview on the VB philosophy in radiation oncology by introducing new phenomenological approaches to NTCP modelling that accounts for inhomogeneous organ radiosensitivity.

Genomic and radiomic data integration, namely radiogenomics, can provide meaningful knowledge in cancer diagnosis, prognosis and treatment. Despite several data structures based on multi-layer architecture proposed to combine multi-omic biological information, none of these has been designed and assessed to include radiomic data as well. To meet this need, we propose to use the MultiAssayExperiment (MAE), an R package that provides data structures and methods for manipulating and integrating multi-assay experiments, as a suitable tool to manage radiogenomic experiment data. To this aim, we first examine the role of radiogenomics in cancer phenotype definition, then the current state of radiogenomics data integration in public repository and, finally, challenges and limitations of including radiomics in MAE, designing an extended framework and showing its application on a case study from the TCGA-TCIA archives. Radiomic and genomic data from 91 patients have been successfully integrated in a single MAE object, demonstrating the suitability of the MAE data structure as container of radiogenomic data.

Purpose Radiation pneumonitis (RP) is commonly associated with thoracic radiation therapy, and its incidence is related to dose and volume of the normal lung in the path of radiation. Our aim was to investigate dose patterns associated with RP in patients enrolled in a randomized trial of intensity modulated radiation therapy (IMRT) versus passive scattering proton therapy (PSPT) for locally advanced non-small cell lung cancer. Methods We analyzed 178 patients prospectively treated with PSPT or IMRT for non-small cell lung cancer to a prescribed dose of 66 or 74 Gy in conventional daily fractionation with concurrent chemotherapy. Forty patients (22%) developed clinically symptomatic RP. Voxel-based analysis of local dose differences was done with a nonparametric permutation test accounting for multiple comparisons. From the obtained 3-dimensional significance maps, we derived clusters of voxels that exhibited dose differences between groups at a statistical significance level of 0.05. Results The voxel-based analysis highlighted that (1) significant dose differences between patients with and without RP were found in the lower part of the lungs and in the heart and (2) the anatomic regions significantly spared by PSPT and the clusters in which doses were significantly correlated with RP development were disjoint. Conclusions The analyzed trial data provide an unprecedented opportunity to substantiate previous hypotheses regarding the role of the heart and the lower lungs in the development of RP. Knowledge of this relationship between RP and thoracic regional radiosensitivity should be considered in clinical practice and in the design of future trials.

To introduce a class of fast 3D quantitative MRI (qMRI) schemes (RESUMEAn arbitrary multi-echo GRE acquisition optimized for quantitative susceptibility mapping (QSM) is complemented with an appropriate low flip-angle GRE sequence drawn from four possible choices. The acquired signals are processed to analytically derive the longitudinal relaxation (RThe RESUMEThe introduced schemes represent a class of robust and flexible strategies to derive a thorough and fast qMRI study, suitable for a whole-brain acquisition with isotropic voxel resolution of 700 μm in less than 15 min.

Purpose To investigate thoracic dose–response patterns for pericardial effusion (PCE) and mortality in patients treated for locally advanced Non-Small-Cell Lung Cancer (NSCLC) by Intensity Modulated RT (IMRT) or Passive-Scattering Proton Therapy (PSPT). Methods Among 178 patients, 43.5% developed grade ≥ 2 PCE. Clinical and dosimetric factors associated with PCE or overall survival (OS) were identified via multi-variable Cox proportional hazards modeling. The Voxel-Based Analyses (VBAs) of local dose differences between patients with and without PCE and mortality was performed. The robustness of VBA results was assessed by a novel characterization of spatial properties of dose distributions based on probabilistic independent component analysis (PICA) and connectograms. Results Several non-dosimetric variables were selected by the multivariable analysis for the considered outcomes, while the time-dependent PCE onset was uncorrelated with the OS (p = 0.34) at a multi-variable Cox analysis. Despite the significant PSPT dosimetric advantage, the RT technique did not affect the occurrence of PCE or OS. VBAs highlighted largely overlapping clusters significantly associated with PCE endpoints in heart and lungs. No significant dosimetric patterns related to mortality endpoints were found. PICA identified 43 components homogeneously scattered within thorax, while connectograms showed modest correlations between doses in main cardio-pulmonary substructures. Conclusions Spatially resolved analysis highlighted dose patterns related to radiation-induced cardiac toxiciy and the observed organ-based dose–response mismatch in PSPT and IMRT. Indeed, the thoracic regions spared by PSPT poorly overlapped with the areas involved in PCE development, as highlited by VBA. PICA and connectograms proved valuable tools for assessing the robusteness of obtained VBA inferences.

Simple Summary The pathophysiology of radiation pneumonitis (RP) after thoracic cancer radiation treatments is still not completely understood although the identification of underlying RP mechanisms may improve the therapeutic window of thoracic cancer patients. The aim of our retrospective study was to explore the dose–response patterns associated with RP by a multi-center voxel-based analysis. In a heterogeneously treated population of 382 thoracic cancer patients, we confirmed the previously described heart–lung interaction in the development of RP. The empowerment of VBA with a novel description of dose map spatial properties based on probabilistic independent component analysis (PICA) and connectograms provided valuable additional and independent information on the radiobiology of RP. Abstract This study investigates the dose–response patterns associated with radiation pneumonitis (RP) in patients treated for thoracic malignancies with different radiation modalities. To this end, voxel-based analysis (VBA) empowered by a novel strategy for the characterization of spatial properties of dose maps was applied. Data from 382 lung cancer and mediastinal lymphoma patients from three institutions treated with different radiation therapy (RT) techniques were analyzed. Each planning CT and biologically effective dose map (α/β = 3 Gy) was spatially normalized on a common anatomical reference. The VBA of local dose differences between patients with and without RP was performed and the clusters of voxels with dose differences that significantly correlated with RP at a p-level of 0.05 were generated accordingly. The robustness of VBA inference was evaluated by a novel characterization for spatial properties of dose maps based on probabilistic independent component analysis (PICA) and connectograms. This lays robust foundations to the obtained findings that the lower parts of the lungs and the heart play a prominent role in the development of RP. Connectograms showed that the dataset can support a radiobiological differentiation between the main heart and lung substructures.

Despite the key-role of the Prostate Imaging and Reporting and Data System (PI-RADS) in the diagnosis and characterization of prostate cancer (PCa), this system remains to be affected by several limitations, primarily associated with the interpretation of equivocal PI-RADS 3 lesions and with the debated role of Dynamic Contrast Enhanced-Magnetic Resonance Imaging (DCE-MRI), which is only used to upgrade peripheral PI-RADS category 3 lesions to PI-RADS category 4 if enhancement is focal. We aimed at investigating the usefulness of radiomics for detection of PCa lesions (Gleason Score ≥ 6) in PI-RADS 3 lesions and in peripheral PI-RADS 3 upgraded to PI-RADS 4 lesions (upPI-RADS 4). Multiparametric MRI (mpMRI) data of patients who underwent prostatic mpMRI between April 2013 and September 2018 were retrospectively evaluated. Biopsy results were used as gold standard. PI-RADS 3 and PI-RADS 4 lesions were re-scored according to the PI-RADS v2.1 before and after DCE-MRI evaluation. Radiomic features were extracted from T2-weighted MRI (T2), Apparent diffusion Coefficient (ADC) map and DCE-MRI subtracted images using PyRadiomics. Feature selection was performed using Wilcoxon-ranksum test and Minimum Redundancy Maximum Relevance (mRMR). Predictive models were constructed for PCa detection in PI-RADS 3 and upPI-RADS 4 lesions using at each step an imbalance-adjusted bootstrap resampling (IABR) on 1000 samples. 41 PI-RADS 3 and 32 upPI-RADS 4 lesions were analyzed. Among 293 radiomic features, the top selected features derived from T2 and ADC. For PI-RADS 3 stratification, second order model showed higher performances (Area Under the Receiver Operating Characteristic Curve—AUC— = 80%), while for upPI-RADS 4 stratification, first order model showed higher performances respect to superior order models (AUC = 89%). Our results support the significant role of T2 and ADC radiomic features for PCa detection in lesions scored as PI-RADS 3 and upPI-RADS 4. Radiomics models showed high diagnostic efficacy in classify PI-RADS 3 and upPI-RADS 4 lesions, outperforming PI-RADS v2.1 performance.

To develop and optimize a rapid magnetic resonance imaging (MRI) screening protocol for pancreatic cancer to be performed in conjunction with breast MRI screening in breast cancer susceptibility gene (BRCA)–positive individuals. An IRB-approved prospective study was conducted. The rapid screening pancreatic MR protocol was designed to be less than 10 min to be performed after a standard breast MRI protocol. Protocol consisted of coronal NT T2 SSFSE, axial NT T2 SSFSE and axial NT rFOV FOCUS DWI, and axial T1. Images were acquired with the patient in the same prone position of breast MRI using the built-in body coil. Image quality was qualitatively assessed by two radiologists with 12 and 13 years of MRI experience, respectively. The imaging protocol was modified until an endpoint of five consecutive patients with high-quality diagnostic images were achieved. Signal-to-noise ratio and contrast-to-noise ratio were assessed. The rapid pancreas MR protocol was successfully completed in all patients. Diagnostic image quality was achieved for all patients. Excellent image quality was achieved for low b values; however, image quality at higher b values was more variable. In one patient, a pancreatic neuroendocrine tumor was found and the patient was treated surgically. In four patients, small pancreatic cystic lesions were detected. In one subject, a hepatic mass was identified and confirmed as adenoma by liver MRI. Rapid MR protocol for pancreatic cancer screening is feasible and has the potential to play a role in screening BRCA patients undergoing breast MRI. • Develop and optimize a rapid magnetic resonance imaging (MRI) screening protocol for pancreatic cancer to be performed in conjunction with breast MRI screening in BRCA mutation positive individuals.

The aim of our study is to characterize the risk of radiation-induced esophagitis (RE) in a cohort of Non-Small-Cell Lung Cancer (NSCLC) patients treated with concurrent chemotherapy and photon/proton therapy. For each patient, the RE was graded according to the CTCAE v.3. The esophageal dose-volume histograms (DVHs) were extracted. Voxel-based analyses (VBAs) were performed to assess the spatial patterns of the dose differences between patients with and without RE of grade ≥ 2. Two hierarchical NTCP models were developed by multivariable stepwise logistic regression based on non-dosimetric factors and on the DVH metrics for the whole esophagus and its anatomical subsites identified by the VBA. In the 173 analyzed patients, 76 (44%) developed RE of grade ≥ 2 at a median follow-up time of 31 days. The VBA identified regions of significant association between dose and RE in a region encompassing the thoracic esophagus. We developed two NTCP models, including the RT modality and a dosimetric factor: V55Gy for the model related to the whole esophagus, and the mean dose for the model designed on the thoracic esophagus. The cross-validated performance showed good predictions for both models (ROC-AUC of 0.70 and 0.73, respectively). The only slight improvement provided by the analysis of the thoracic esophageal subsites might be due to the relevant sparing of cervical and lower thoracic esophagus in the analyzed cohort. Further studies on larger cohorts and a more heterogeneous set of dose distributions are needed to validate these preliminary findings and shed further light on the spatial patterns of RE development.

Introduction Chemotherapy prolongs the survival of patients with advanced and metastatic tumors. Since the liver plays an active role in the metabolism of chemotherapy agents, hepatic injury is a common adverse effect. The purpose of this study is to compare a novel quantitative chemical shift encoded magnetic resonance imaging (CSE-MRI) method with conventional T1-weighted In and Out of phase (T1 IOP) MR for evaluating the reproducibility of the methods in an oncologic population exposed to chemotherapy. Materials and methods This retrospective study was approved by the institutional review board with a waiver for informed consent. The study included patients who underwent chemotherapy, no suspected liver iron overload, and underwent upper abdomen MRI. Two radiologists independently draw circular ROIsin the liver parenchyma. The fat fraction was calculated from IOP imaging and measured from IDEAL-IQ fat fraction maps. Two different equations were used to estimate fat with IOP sequences. Intra-class correlation coefficient and repeatability coefficient were estimated to evaluate agreement between two readers on iron level and fat fraction measurement. Results CSE-MRI showed a higher reliability in fat quantification compared with both IOP methods, with a substantially higher inter-reader agreement (0.961 vs 0.372). This has important clinical implications. Conclusion The novel CSE-MRI method described here provides increased reproducibility and confidence in diagnosing hepatic steatosis in a oncologic clinical setting. IDEAL-IQ has been proved to be more reproducible than conventional IOP imaging.

Background Liver cancer is the sixth most common cancer worldwide and the second leading cause of cancer mortality. Chronic liver disease caused by viral infection, alcohol abuse, or other factors can lead to cirrhosis. Cirrhosis is the most important clinical risk factor for hepatocellular carcinoma (HCC) whereby the normal hepatic architecture is replaced by fibrous septa and a spectrum of nodules ranging from benign regenerative nodules to HCC, each one of them with different imaging features. Conclusions Furthermore, advanced techniques including the quantification of hepatic and intralesional fat and iron, magnetic resonance elastography, radiomics, radiogenomics, and positron emission tomography (PET)-MRI are highly promising for the extraction of new imaging biomarkers that reflect the tumor microenvironment and, in the future, may add decision-making value in the management of patients with HCC.

PURPOSE: To investigate the value of second-opinion interpretation of cross-sectional images by subspecialized radiologists to diagnose recurrent pancreatic cancer after surgery. METHODS: The IRB approved and issued a waiver of informed consent for this retrospective study. Initial and second-opinion interpretations of 69 consecutive submitted MRI or CT follow up after pancreatic cancer resection between January 1, 2009 and December 31, 2013 were evaluated by one oncologic imaging radiologist, who was blinded to patient’s clinical details and histopathologic data. The reviewer was asked to classify each interpretation in reference of the diagnosis of PDAC recurrence. It was also recorded if the radiologic interpretation recommended additional imaging studies to confirm recurrence. The diagnosis of recurrence was determined by pathology when available, otherwise by imaging follow up, clinical or laboratory assessments. Cohen’s kappa statistic was used to assess agreement between initial and second-opinion interpretations. The differences between the initial and second-opinion interpretations were examined using McNemar test or Bowker’s test of symmetry. RESULTS: Disagreement on recurrence between the initial report and the second-opinion interpretation was observed in 32% of cases (22/69; k=0.44). Second-opinion interpretations had a higher sensitivity and a higher specificity on recurrence compared to the initial interpretations (0.93 vs 0.75 and 0.90 vs 0.68, respectively), and the difference in specificity was significant (p=0.016). Additional imaging studies were recommended more frequently in the initial interpretation (22% vs 6%, p=0.006). CONCLUSIONS: Our study shows the second-opinion interpretation by subspecialized radiologists improves the detection of pancreatic cancer recurrence after surgical resection.

Objective: The aims of this study were to: a) evaluate tumor response rates using modified-Response-evaluation-criteria-in-solid-tumors (mRecist) criteria, b) evaluate safety of Degradable Starch Microspheres Trans-arterial-chemo-embolization (DSMs-TACE) for unresectable hepatocellular-carcinoma (HCC) treatment. Materials and Methods: We prospectively enrolled 24 HCC cirrhotic patients (21/3 M/F, mean age 66.3 years) to be treated with repeated DSMs-TACE procedures, performed at 4-6 week intervals on the basis of tumor response and patients tolerance. Clinical and biochemical evaluations were performed before and after each procedure. Treatment response was also assessed by Computed-tomography (CT) or Magnetic-resonance-imaging (MRI)-scan 4-6 weeks following each procedure. Results: In our experience, DSMs-TACE was both safe and effective. A total of 53 DSMs-TACE procedures were performed (2.2 per patient). No procedure-related death was observed. Complete Response (CR) was observed in 5/24 (20.8%), 4/17 (23.5%) and 5/12 (41.6%) patients after the first, second and third procedure, respectively. At the end of each treatment, all patients experienced at least a partial response. At the end of the repeated procedures, no differences between mono- or bi-lobar disease were observed in patients with CR (64.2% vs 50%; p=ns). In most cases, treatment discontinuation was due to worsening liver function. Conclusion: DSMs-TACE is a valid, well-tolerated alternative treatment to Lipiodol-TACE or DEB-TACE, as it has demonstrated to achieve a relatively high percentage of complete tumor necrosis. CR rates were similar between patients with mono- or bi-lobar disease indicating the possibility of carrying-out repeated procedure in a safe and effective way in both types of patients.

To devise a novel Spatial Normalization framework for Voxel-based analysis (VBA) in brain radiotherapy. VBAs rely on accurate spatial normalization of different patients' planning CTs on a common coordinate system (CCS). The cerebral anatomy, well characterized by MRI, shows instead poor contrast in CT, resulting in potential inaccuracies in VBAs based on CT alone.We analyzed 50 meningioma patients treated with proton-therapy, undergoing planning CT and T1-weighted (T1w) MRI. The spatial normalization pipeline based on MR and CT images consisted in: intra-patient registration of CT to T1w, inter-patient registration of T1w to MNI space chosen as CCS, doses propagation to MNI. The registration quality was compared with that obtained by Statistical Parametric Mapping software (SPM), used as benchmark. To evaluate the accuracy of dose normalization, the dose organ overlap (DOO) score was computed on gray matter, white matter and cerebrospinal fluid before and after normalization. In addition, the trends in the DOOs distribution were investigated by means of cluster analysis.The registration quality was higher for the proposed method compared to SPM (p 0.001). The DOO scores showed a significant improvement after normalization (p 0.001). The cluster analysis highlighted 2 clusters, with one of them including the majority of data and exhibiting acceptable DOOs.Our study presents a robust tool for spatial normalization, specifically tailored for brain dose VBAs. Furthermore, the cluster analysis provides a formal criterion for patient exclusion in case of non-acceptable normalization results. The implemented framework lays the groundwork for future reliable VBAs in brain irradiation studies.

Prostate cancer (PCa) is a disease affecting an increasing number of men worldwide. Several efforts have been made to identify imaging biomarkers to non-invasively detect and characterize PCa, with substantial improvements thanks to multiparametric Magnetic Resonance Imaging (mpMRI). In recent years, diffusion kurtosis imaging (DKI) was proposed to be directly related to tissue physiological and pathological characteristic, while the radiomic approach was proven to be a key method to study cancer imaging phenotypes. Our aim was to compare a standard radiomic model for PCa detection, built using T2-weighted (T2W) and Apparent Diffusion Coefficient (ADC), with an advanced one, including DKI and quantitative Dynamic Contrast Enhanced (DCE), while also evaluating differences in prediction performance when using 2D or 3D lesion segmentation. The obtained results in terms of diagnostic accuracy were high for all of the performed comparisons, reaching values up to 0.99 for the area under a receiver operating characteristic curve (AUC), and 0.98 for both sensitivity and specificity. In comparison, the radiomic model based on standard features led to prediction performances higher than those of the advanced model, while greater accuracy was achieved by the model extracted from 3D segmentation. These results provide new insights into active topics of discussion, such as choosing the most convenient acquisition protocol and the most appropriate postprocessing pipeline to accurately detect and characterize PCa.

The importance of Diffusion Weighted Imaging (DWI) in prostate cancer (PCa) diagnosis have been widely handled in literature. In the last decade, due to the mono-exponential model limitations, several studies investigated non-Gaussian DWI models and their utility in PCa diagnosis. Since their results were often inconsistent and conflicting, we performed a systematic review of studies from 2012 examining the most commonly used Non-Gaussian DWI models for PCa detection and characterization. A meta-analysis was conducted to assess the ability of each Non-Gaussian model to detect PCa lesions and distinguish between low and intermediate/high grade lesions. Weighted mean differences and 95% confidence intervals were calculated and the heterogeneity was estimated using the I2 statistic. 29 studies were selected for the systematic review, whose results showed inconsistence and an unclear idea about the actual usefulness and the added value of the Non-Gaussian model parameters. 12 studies were considered in the meta-analyses, which showed statistical significance for several non-Gaussian parameters for PCa detection, and to a lesser extent for PCa characterization. Our findings showed that Non-Gaussian model parameters may potentially play a role in the detection and characterization of PCa but further studies are required to identify a standardized DWI acquisition protocol for PCa diagnosis.

Thoracic radiation therapy (RT) is often associated with lung side effects, whose etiology is still controversial. Our aim was to explore correlations between local dose in the thoracic anatomy and the radiation-induced lung damage (RILD). To this end, we designed a robust scheme for voxel-based analysis (VBA) to explore dose patterns associated with RILD in non-small-cell lung cancer (NSCLC) patients receiving stereotactic body RT (SBRT). We analyzed 106 NSCLC SBRT patients (median prescription dose: 50 Gy; range: [40-54] Gy) in 4 fractions (range: [3-5]) with clinical and dosimetric records suitable for the analysis. The incidence of acute G1 RILD (RTOG grade ⩾ 1) was 68%. Each planning CT and dose map was spatially normalized to a common anatomical reference using a B-spline inter-patient registration algorithm after masking the gross tumor volume. The tumor-subtracted dose maps were converted into biologically effective dose maps (α/β = 3 Gy). VBA was performed according to a non-parametric permutation test accounting for multiple comparison, based on a cluster analysis method. The underlying general linear model of RILD was designed to include dose maps and each non-dosimetric variable significantly correlated with RILD. The clusters of voxels with dose differences significantly correlated with RILD at a given p -level (S p ) were generated. The only non-dosimetric variable significantly correlated with RILD was the chronic obstructive pulmonary disease (p = 0.034). Patients with G1 RILD received significantly (p ⩽ 0.05) higher doses in two voxel clusters S 0.05 in the lower-left lung (14 cm3) and in an area (64 cm3) largely included within the ventricles. The applied VBA represents a powerful tool to probe the dose susceptibility of inhomogeneous organs in clinical radiobiology studies. The identified subregions with dose differences associated with G1 RILD in both the heart and lower lungs endorse a trend of previously reported hypotheses on lung toxicity radiobiology.

INTRODUCTION Chronic liver disease resulting in fibrosis, and ultimately cirrhosis, is a significant cause of morbidity and mortality worldwide. None of the conventional imaging techniques are able to detect early fibrosis and compare its grade with the histopathologic scale. Liver biopsy, as the diagnostic standard for liver fibrosis, also has limitations and is not well accepted by patients. Magnetic resonance elastography is a well-established technique for evaluating liver stiffness and may replace invasive procedures. Detection of liver fibrosis in its early stages, however, requires a detailed knowledge of normal liver stiffness. OBJECTIVES This study aimed to determine normal liver stiffness values in healthy volunteers. PATIENTS AND METHODS A total of 102 volunteers (mean age, 21.6 years; range, 20–28 years) with no history of gastrointestinal, hepatobiliary, or cardiovascular disease were enrolled in the study. Liver stiffness was evaluated by magnetic resonance elastography with a 1.5T clinical magnetic resonance scanner. Images of the induced transverse wave propagation were obtained and converted to tissue stiffness maps (elastograms). RESULTS The mean (SD) liver stiffness for the entire group of volunteers was 2.14 (0.28) kPa (range, 1.37–2.66 kPa). For women, the mean (SD) stiffness value was 2.14 (0.30) kPa (range, 1.37–2.66 kPa), and for men, 2.14 (0.25) kPa (range, 1.54–2.54 kPa). CONCLUSIONS Liver stiffness in a healthy adult cohort did not exceed 2.7 kPa and is not influenced by sex, body mass index, or fat content.

Mathematical models of normal tissue complication probability (NTCP) able to robustly predict radiation-induced morbidities (RIM) play an essential role in the identification of a personalized optimal plan, and represent the key to maximizing the benefits of technological advances in radiation therapy (RT). Most modern RT techniques pose, however, new challenges in estimating the risk of RIM. The aim of this report is to schematically review NTCP models in the framework of advanced radiation therapy techniques. Issues relevant to hypofractionated stereotactic body RT and ion beam therapy are critically reviewed. Reirradiation scenarios for new or recurrent malignances and NTCP are also illustrated. A new phenomenological approach to predict RIM is suggested.

In radiation oncology, the need for a modern Normal Tissue Complication Probability (NTCP) philosophy to include voxel-based evidence on organ radio-sensitivity (RS) has been acknowledged. Here a new formalism (Probabilistic Atlas for Complication Estimation, PACE) to predict radiation-induced morbidity (RIM) is presented. The adopted strategy basically consists in keeping the structure of a classical, phenomenological NTCP model, such as the Lyman-Kutcher-Burman (LKB), and replacing the dose distribution with a collection of RIM odds, including also significant non-dosimetric covariates, as input of the model framework. The theory was first demonstrated in silico on synthetic dose maps, classified according to synthetic outcomes. PACE was then applied to a clinical dataset of thoracic cancer patients classified for lung fibrosis. LKB models were trained for comparison. Overall, the obtained learning curves showed that the PACE model outperformed the LKB and predicted synthetic outcomes with an accuracy >0.8. On the real patients, PACE performance, evaluated by both discrimination and calibration, was significantly higher than LKB. This trend was confirmed by cross-validation. Furthermore, the capability to infer the spatial pattern of underlying RS map for the analyzed RIM was successfully demonstrated, thus paving the way to new perspectives of NTCP models as learning tools.