Investigation of a potential upstream harmonization based on image appearance matching to improve radiomics features robustness: a phantom study

Objective. Radiomics is a promising valuable analysis tool consisting in extracting quantitative information from medical images. However, the extracted radiomics features are too sensitive to variations in used image acquisition and reconstruction parameters. This limited robustness hinders the generalizable validity of radiomics-assisted models. Our aim is to investigate a possible harmonization strategy based on matching image quality to improve feature robustness.Approach.We acquired CT scans of a phantom with two scanners across different dose levels and percentages of Iterative Reconstruction algorithms. The detectability index was used as a comprehensive task-based image quality metric. A statistical analysis based on the Intraclass Correlation Coefficient was performed to determine if matching image quality/appearance could enhance the robustness of radiomics features extracted from the phantom images. Additionally, an Artificial Neural Network was trained on these features to automatically classify the scanner used for image acquisition.Main results.We found that the ICC of the features across protocols providing a similar detectability index improves with respect to the ICC of the features across protocols providing a different detectability index. This improvement was particularly noticeable in features relevant for distinguishing between scanners.Significance.This preliminary study demonstrates that a harmonization based on image quality/appearance matching could improve radiomics features robustness and heterogeneous protocols can be used to obtain a similar image appearance in terms of the detectability index. Thus protocols with a lower dose level could be selected to reduce the amount of radiation dose delivered to the patient and simultaneously obtain a more robust quantitative analysis.

Radiomic- and dosiomic-based clustering development for radio-induced neurotoxicity in pediatric medulloblastoma

BACKGROUND

Texture analysis extracts many quantitative image features, offering a valuable, cost-effective, and non-invasive approach for individual medicine. Furthermore, multimodal machine learning could have a large impact for precision medicine, as texture biomarkers can underlie tissue microstructure. This study aims to investigate imaging-based biomarkers of radio-induced neurotoxicity in pediatric patients with metastatic medulloblastoma, using radiomic and dosiomic analysis.

METHODS

This single-center study retrospectively enrolled children diagnosed with metastatic medulloblastoma (MB) and treated with hyperfractionated craniospinal irradiation (CSI). Histological confirmation of medulloblastoma and baseline follow-up magnetic resonance imaging (MRI) were mandatory. Treatment involved helical tomotherapy (HT) delivering a dose of 39 Gray (Gy) to brain and spinal axis and a posterior fossa boost up to 60 Gy. Clinical outcomes, such as local and distant brain control and neurotoxicity, were recorded. Radiomic and dosiomic features were extracted from tumor regions on T1, T2, FLAIR (fluid-attenuated inversion recovery) MRI-maps, and radiotherapy dose distribution. Different machine learning feature selection and reduction approaches were performed for supervised and unsupervised clustering.

RESULTS

Forty-eight metastatic medulloblastoma patients (29 males and 19 females) with a mean age of 12 ± 6 years were enrolled. For each patient, 332 features were extracted. Greater level of abstraction of input data by combining selection of most performing features and dimensionality reduction returns the best performance. The resulting one-component radiomic signature yielded an accuracy of 0.73 with sensitivity, specificity, and precision of 0.83, 0.64, and 0.68, respectively.

CONCLUSIONS

Machine learning radiomic-dosiomic approach effectively stratified pediatric medulloblastoma patients who experienced radio-induced neurotoxicity. Strategy needs further validation in external dataset for its potential clinical use in ab initio management paradigms of medulloblastoma.

Using the gamma-index analysis for inter-fractional comparison of in-beam PET images for head-and-neck treatment monitoring in proton therapy: A Monte Carlo simulation study

GOAL

In-beam Positron Emission Tomography (PET) is a technique for in-vivo non-invasive treatment monitoring for proton therapy. To detect anatomical changes in patients with PET, various analysis methods exist, but their clinical interpretation is problematic. The goal of this work is to investigate whether the gamma-index analysis, widely used for dose comparisons, is an appropriate tool for comparing in-beam PET distributions. Focusing on a head-and-neck patient, we investigate whether the gamma-index map and the passing rate are sensitive to progressive anatomical changes.

METHODS/MATERIALS

We simulated a treatment course of a proton therapy patient using FLUKA Monte Carlo simulations. Gradual emptying of the sinonasal cavity was modeled through a series of artificially modified CT scans. The in-beam PET activity distributions from three fields were evaluated, simulating a planar dual head geometry. We applied the 3D-gamma evaluation method to compare the PET images with a reference image without changes. Various tolerance criteria and parameters were tested, and results were compared to the CT-scans.

RESULTS

Based on 210 MC simulations we identified appropriate parameters for the gamma-index analysis. Tolerance values of 3 mm/3% and 2 mm/2% were suited for comparison of simulated in-beam PET distributions. The gamma passing rate decreased with increasing volume change for all fields.

CONCLUSION

The gamma-index analysis was found to be a useful tool for comparing simulated in-beam PET images, sensitive to sinonasal cavity emptying. Monitoring the gamma passing rate behavior over the treatment course is useful to detect anatomical changes occurring during the treatment course.

Synthetic CT imaging for PET monitoring in proton therapy: a simulation study

Objective.This study addresses a fundamental limitation of in-beam positron emission tomography (IB-PET) in proton therapy: the lack of direct anatomical representation in the images it produces. We aim to overcome this shortcoming by pioneering the application of deep learning techniques to create synthetic control CT images (sCT) from combining IB-PET and planning CT scan data.Approach.We conducted simulations involving six patients who underwent irradiation with proton beams. Leveraging the architecture of a visual transformer (ViT) neural network, we developed a model to generate sCT images of these patients using the planning CT scans and the inter-fractional simulated PET activity maps during irradiation. To evaluate the model's performance, a comparison was conducted between the sCT images produced by the ViT model and the authentic control CT images-serving as the benchmark.Main results.The structural similarity index was computed at a mean value across all patients of 0.91, while the mean absolute error measured 22 Hounsfield Units (HU). Root mean squared error and peak signal-to-noise ratio values were 56 HU and 30 dB, respectively. The Dice similarity coefficient exhibited a value of 0.98. These values are comparable to or exceed those found in the literature. More than 70% of the synthetic morphological changes were found to be geometrically compatible with the ones reported in the real control CT scan.Significance.Our study presents an innovative approach to surface the hidden anatomical information of IB-PET in proton therapy. Our ViT-based model successfully generates sCT images from inter-fractional PET data and planning CT scans. Our model's performance stands on par with existing models relying on input from cone beam CT or magnetic resonance imaging, which contain more anatomical information than activity maps.

Deep learning based joint fusion approach to exploit anatomical and functional brain information in autism spectrum disorders

BACKGROUND

The integration of the information encoded in multiparametric MRI images can enhance the performance of machine-learning classifiers. In this study, we investigate whether the combination of structural and functional MRI might improve the performances of a deep learning (DL) model trained to discriminate subjects with Autism Spectrum Disorders (ASD) with respect to typically developing controls (TD).

MATERIAL AND METHODS

We analyzed both structural and functional MRI brain scans publicly available within the ABIDE I and II data collections. We considered 1383 male subjects with age between 5 and 40 years, including 680 subjects with ASD and 703 TD from 35 different acquisition sites. We extracted morphometric and functional brain features from MRI scans with the Freesurfer and the CPAC analysis packages, respectively. Then, due to the multisite nature of the dataset, we implemented a data harmonization protocol. The ASD vs. TD classification was carried out with a multiple-input DL model, consisting in a neural network which generates a fixed-length feature representation of the data of each modality (FR-NN), and a Dense Neural Network for classification (C-NN). Specifically, we implemented a joint fusion approach to multiple source data integration. The main advantage of the latter is that the loss is propagated back to the FR-NN during the training, thus creating informative feature representations for each data modality. Then, a C-NN, with a number of layers and neurons per layer to be optimized during the model training, performs the ASD-TD discrimination. The performance was evaluated by computing the Area under the Receiver Operating Characteristic curve within a nested 10-fold cross-validation. The brain features that drive the DL classification were identified by the SHAP explainability framework.

RESULTS

The AUC values of 0.66±0.05 and of 0.76±0.04 were obtained in the ASD vs. TD discrimination when only structural or functional features are considered, respectively. The joint fusion approach led to an AUC of 0.78±0.04. The set of structural and functional connectivity features identified as the most important for the two-class discrimination supports the idea that brain changes tend to occur in individuals with ASD in regions belonging to the Default Mode Network and to the Social Brain.

CONCLUSIONS

Our results demonstrate that the multimodal joint fusion approach outperforms the classification results obtained with data acquired by a single MRI modality as it efficiently exploits the complementarity of structural and functional brain information.

Effect of data harmonization of multicentric dataset in ASD/TD classification

Machine Learning (ML) is nowadays an essential tool in the analysis of Magnetic Resonance Imaging (MRI) data, in particular in the identification of brain correlates in neurological and neurodevelopmental disorders. ML requires datasets of appropriate size for training, which in neuroimaging are typically obtained collecting data from multiple acquisition centers. However, analyzing large multicentric datasets can introduce bias due to differences between acquisition centers. ComBat harmonization is commonly used to address batch effects, but it can lead to data leakage when the entire dataset is used to estimate model parameters. In this study, structural and functional MRI data from the Autism Brain Imaging Data Exchange (ABIDE) collection were used to classify subjects with Autism Spectrum Disorders (ASD) compared to Typical Developing controls (TD). We compared the classical approach (external harmonization) in which harmonization is performed before train/test split, with an harmonization calculated only on the train set (internal harmonization), and with the dataset with no harmonization. The results showed that harmonization using the whole dataset achieved higher discrimination performance, while non-harmonized data and harmonization using only the train set showed similar results, for both structural and connectivity features. We also showed that the higher performances of the external harmonization are not due to larger size of the sample for the estimation of the model and hence these improved performance with the entire dataset may be ascribed to data leakage. In order to prevent this leakage, it is recommended to define the harmonization model solely using the train set.

Quality assessment, variability and reproducibility of anatomical measurements derived from T1-weighted brain imaging: The RIN-Neuroimaging Network case study

Initiatives for the collection of harmonized MRI datasets are growing continuously, opening questions on the reliability of results obtained in multi-site contexts. Here we present the assessment of the brain anatomical variability of MRI-derived measurements obtained from T1-weighted images, acquired according to the Standard Operating Procedures, promoted by the RIN-Neuroimaging Network. A multicentric dataset composed of 77 brain T1w acquisitions of young healthy volunteers (mean age = 29.7 ± 5.0 years), collected in 15 sites with MRI scanners of three different vendors, was considered. Parallelly, a dataset of 7 "traveling" subjects, each undergoing three acquisitions with scanners from different vendors, was also used. Intra-site, intra-vendor, and inter-site variabilities were evaluated in terms of the percentage standard deviation of volumetric and cortical thickness measures. Image quality metrics such as contrast-to-noise and signal-to-noise ratio in gray and white matter were also assessed for all sites and vendors. The results showed a measured global variability that ranges from 11% to 19% for subcortical volumes and from 3% to 10% for cortical thicknesses. Univariate distributions of the normalized volumes of subcortical regions, as well as the distributions of the thickness of cortical parcels appeared to be significantly different among sites in 8 subcortical (out of 17) and 21 cortical (out of 68) regions of i nterest in the multicentric study. The Bland-Altman analysis on "traveling" brain measurements did not detect systematic scanner biases even though a multivariate classification approach was able to classify the scanner vendor from brain measures with an accuracy of 0.60 ± 0.14 (chance level 0.33).

Quantification of pulmonary involvement in COVID-19 pneumonia: an upgrade of the LungQuant software for lung CT segmentation

Computed tomography (CT) scans are used to evaluate the severity of lung involvement in patients affected by COVID-19 pneumonia. Here, we present an improved version of the LungQuant automatic segmentation software (LungQuant v2), which implements a cascade of three deep neural networks (DNNs) to segment the lungs and the lung lesions associated with COVID-19 pneumonia. The first network (BB-net) defines a bounding box enclosing the lungs, the second one (U-net 1 ) outputs the mask of the lungs, and the final one (U-net 2 ) generates the mask of the COVID-19 lesions. With respect to the previous version (LungQuant v1), three main improvements are introduced: the BB-net, a new term in the loss function in the U-net for lesion segmentation and a post-processing procedure to separate the right and left lungs. The three DNNs were optimized, trained and tested on publicly available CT scans. We evaluated the system segmentation capability on an independent test set consisting of ten fully annotated CT scans, the COVID-19-CT-Seg benchmark dataset. The test performances are reported by means of the volumetric dice similarity coefficient (vDSC) and the surface dice similarity coefficient (sDSC) between the reference and the segmented objects. LungQuant v2 achieves a vDSC (sDSC) equal to 0.96 ± 0.01 (0.97 ± 0.01) and 0.69 ± 0.08 (0.83 ± 0.07) for the lung and lesion segmentations, respectively. The output of the segmentation software was then used to assess the percentage of infected lungs, obtaining a Mean Absolute Error (MAE) equal to 2%.

A multicenter evaluation of a deep learning software (LungQuant) for lung parenchyma characterization in COVID-19 pneumonia

BACKGROUND

The role of computed tomography (CT) in the diagnosis and characterization of coronavirus disease 2019 (COVID-19) pneumonia has been widely recognized. We evaluated the performance of a software for quantitative analysis of chest CT, the LungQuant system, by comparing its results with independent visual evaluations by a group of 14 clinical experts. The aim of this work is to evaluate the ability of the automated tool to extract quantitative information from lung CT, relevant for the design of a diagnosis support model.

METHODS

LungQuant segments both the lungs and lesions associated with COVID-19 pneumonia (ground-glass opacities and consolidations) and computes derived quantities corresponding to qualitative characteristics used to clinically assess COVID-19 lesions. The comparison was carried out on 120 publicly available CT scans of patients affected by COVID-19 pneumonia. Scans were scored for four qualitative metrics: percentage of lung involvement, type of lesion, and two disease distribution scores. We evaluated the agreement between the LungQuant output and the visual assessments through receiver operating characteristics area under the curve (AUC) analysis and by fitting a nonlinear regression model.

RESULTS

Despite the rather large heterogeneity in the qualitative labels assigned by the clinical experts for each metric, we found good agreement on the metrics compared to the LungQuant output. The AUC values obtained for the four qualitative metrics were 0.98, 0.85, 0.90, and 0.81.

CONCLUSIONS

Visual clinical evaluation could be complemented and supported by computer-aided quantification, whose values match the average evaluation of several independent clinical experts.

KEY POINTS

We conducted a multicenter evaluation of the deep learning-based LungQuant automated software. We translated qualitative assessments into quantifiable metrics to characterize coronavirus disease 2019 (COVID-19) pneumonia lesions. Comparing the software output to the clinical evaluations, results were satisfactory despite heterogeneity of the clinical evaluations. An automatic quantification tool may contribute to improve the clinical workflow of COVID-19 pneumonia.

Deriving quantitative information from multiparametric MRI via Radiomics: Evaluation of the robustness and predictive value of radiomic features in the discrimination of low-grade versus high-grade gliomas with machine learning

PURPOSE

Analysis pipelines based on the computation of radiomic features on medical images are widely used exploration tools across a large variety of image modalities. This study aims to define a robust processing pipeline based on Radiomics and Machine Learning (ML) to analyze multiparametric Magnetic Resonance Imaging (MRI) data to discriminate between high-grade (HGG) and low-grade (LGG) gliomas.

METHODS

The dataset consists of 158 multiparametric MRI of patients with brain tumor publicly available on The Cancer Imaging Archive, preprocessed by the BraTS organization committee. Three different types of image intensity normalization algorithms were applied and 107 features were extracted for each tumor region, setting the intensity values according to different discretization levels. The predictive power of radiomic features in the LGG versus HGG categorization was evaluated by using random forest classifiers. The impact of the normalization techniques and of the different settings in the image discretization was studied in terms of the classification performances. A set of MRI-reliable features was defined selecting the features extracted according to the most appropriate normalization and discretization settings.

RESULTS

The results show that using MRI-reliable features improves the performance in glioma grade classification (AUC=0.93±0.05) with respect to the use of raw (AUC=0.88±0.08) and robust features (AUC=0.83±0.08), defined as those not depending on image normalization and intensity discretization.

CONCLUSIONS

These results confirm that image normalization and intensity discretization strongly impact the performance of ML classifiers based on radiomic features. Thus, special attention should be provided in the image preprocessing step before typical radiomic and ML analysis are carried out.

In-vivo range verification analysis with in-beam PET data for patients treated with proton therapy at CNAO

Morphological changes that may arise through a treatment course are probably one of the most significant sources of range uncertainty in proton therapy. Non-invasive in-vivo treatment monitoring is useful to increase treatment quality. The INSIDE in-beam Positron Emission Tomography (PET) scanner performs in-vivo range monitoring in proton and carbon therapy treatments at the National Center of Oncological Hadrontherapy (CNAO). It is currently in a clinical trial (ID: NCT03662373) and has acquired in-beam PET data during the treatment of various patients. In this work we analyze the in-beam PET (IB-PET) data of eight patients treated with proton therapy at CNAO. The goal of the analysis is twofold. First, we assess the level of experimental fluctuations in inter-fractional range differences (sensitivity) of the INSIDE PET system by studying patients without morphological changes. Second, we use the obtained results to see whether we can observe anomalously large range variations in patients where morphological changes have occurred. The sensitivity of the INSIDE IB-PET scanner was quantified as the standard deviation of the range difference distributions observed for six patients that did not show morphological changes. Inter-fractional range variations with respect to a reference distribution were estimated using the Most-Likely-Shift (MLS) method. To establish the efficacy of this method, we made a comparison with the Beam’s Eye View (BEV) method. For patients showing no morphological changes in the control CT the average range variation standard deviation was found to be 2.5 mm with the MLS method and 2.3 mm with the BEV method. On the other hand, for patients where some small anatomical changes occurred, we found larger standard deviation values. In these patients we evaluated where anomalous range differences were found and compared them with the CT. We found that the identified regions were mostly in agreement with the morphological changes seen in the CT scan.

Virtual Ontogeny of Cortical Growth Preceding Mental Illness

BACKGROUND

Morphology of the human cerebral cortex differs across psychiatric disorders, with neurobiology and developmental origins mostly undetermined. Deviations in the tangential growth of the cerebral cortex during pre/perinatal periods may be reflected in individual variations in cortical surface area later in life.

METHODS

Interregional profiles of group differences in surface area between cases and controls were generated using T1-weighted magnetic resonance imaging from 27,359 individuals including those with attention-deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, major depressive disorder, schizophrenia, and high general psychopathology (through the Child Behavior Checklist). Similarity of interregional profiles of group differences in surface area and prenatal cell-specific gene expression was assessed.

RESULTS

Across the 11 cortical regions, group differences in cortical area for attention-deficit/hyperactivity disorder, schizophrenia, and Child Behavior Checklist were dominant in multimodal association cortices. The same interregional profiles were also associated with interregional profiles of (prenatal) gene expression specific to proliferative cells, namely radial glia and intermediate progenitor cells (greater expression, larger difference), as well as differentiated cells, namely excitatory neurons and endothelial and mural cells (greater expression, smaller difference). Finally, these cell types were implicated in known pre/perinatal risk factors for psychosis. Genes coexpressed with radial glia were enriched with genes implicated in congenital abnormalities, birth weight, hypoxia, and starvation. Genes coexpressed with endothelial and mural genes were enriched with genes associated with maternal hypertension and preterm birth.

CONCLUSIONS

Our findings support a neurodevelopmental model of vulnerability to mental illness whereby prenatal risk factors acting through cell-specific processes lead to deviations from typical brain development during pregnancy.

Subtly altered topological asymmetry of brain structural covariance networks in autism spectrum disorder across 43 datasets from the ENIGMA consortium

Small average differences in the left-right asymmetry of cerebral cortical thickness have been reported in individuals with autism spectrum disorder (ASD) compared to typically developing controls, affecting widespread cortical regions. The possible impacts of these regional alterations in terms of structural network effects have not previously been characterized. Inter-regional morphological covariance analysis can capture network connectivity between different cortical areas at the macroscale level. Here, we used cortical thickness data from 1455 individuals with ASD and 1560 controls, across 43 independent datasets of the ENIGMA consortium's ASD Working Group, to assess hemispheric asymmetries of intra-individual structural covariance networks, using graph theory-based topological metrics. Compared with typical features of small-world architecture in controls, the ASD sample showed significantly altered average asymmetry of networks involving the fusiform, rostral middle frontal, and medial orbitofrontal cortex, involving higher randomization of the corresponding right-hemispheric networks in ASD. A network involving the superior frontal cortex showed decreased right-hemisphere randomization. Based on comparisons with meta-analyzed functional neuroimaging data, the altered connectivity asymmetry particularly affected networks that subserve executive functions, language-related and sensorimotor processes. These findings provide a network-level characterization of altered left-right brain asymmetry in ASD, based on a large combined sample. Altered asymmetrical brain development in ASD may be partly propagated among spatially distant regions through structural connectivity.

Multi-site harmonization of MRI data uncovers machine-learning discrimination capability in barely separable populations: An example from the ABIDE dataset

•

Multi-site MRI data encode confounding information which may mask case-control differences.

•

The impact of the NeuroHarmonize method is evaluated in the ASD-control classification.

•

We verified the successful removal of the site effect by the harmonization protocol.

•

The increment in the classification performance is quantified after data harmonization.

•

We identified the anatomical features that contributed to the two-class separation.

Machine Learning (ML) techniques have been widely used in Neuroimaging studies of Autism Spectrum Disorders (ASD) both to identify possible brain alterations related to this condition and to evaluate the predictive power of brain imaging modalities. The collection and public sharing of large imaging samples has favored an even greater diffusion of the use of ML-based analyses. However, multi-center data collections may suffer the batch effect, which, especially in case of Magnetic Resonance Imaging (MRI) studies, should be curated to avoid confounding effects for ML classifiers and masking biases. This is particularly important in the study of barely separable populations according to MRI data, such as subjects with ASD compared to controls with typical development (TD). Here, we show how the implementation of a harmo- nization protocol on brain structural features unlocks the case-control ML separation capability in the analysis of a multi-center MRI dataset. This effect is demonstrated on the ABIDE data collection, involving subjects encompassing a wide age range. After data harmonization, the overall ASD vs. TD discrimination capability by a Random Forest (RF) classifier improves from a very low performance (AUC = 0.58 ± 0.04) to a still low, but reasonably significant AUC = 0.67 ± 0.03. The performances of the RF classifier have been evaluated also in the age-specific subgroups of children, adolescents and adults, obtaining AUC = 0.62 ± 0.02, AUC = 0.65 ± 0.03 and AUC = 0.69 ± 0.06, respectively. Specific and consistent patterns of anatomical differences related to the ASD condition have been identified for the three different age subgroups.

Consortium neuroscience of attention deficit/hyperactivity disorder and autism spectrum disorder: The ENIGMA adventure

Neuroimaging has been extensively used to study brain structure and function in individuals with attention deficit/hyperactivity disorder (ADHD) and autism spectrum disorder (ASD) over the past decades. Two of the main shortcomings of the neuroimaging literature of these disorders are the small sample sizes employed and the heterogeneity of methods used. In 2013 and 2014, the ENIGMA-ADHD and ENIGMA-ASD working groups were respectively, founded with a common goal to address these limitations. Here, we provide a narrative review of the thus far completed and still ongoing projects of these working groups. Due to an implicitly hierarchical psychiatric diagnostic classification system, the fields of ADHD and ASD have developed largely in isolation, despite the considerable overlap in the occurrence of the disorders. The collaboration between the ENIGMA-ADHD and -ASD working groups seeks to bring the neuroimaging efforts of the two disorders closer together. The outcomes of case-control studies of subcortical and cortical structures showed that subcortical volumes are similarly affected in ASD and ADHD, albeit with small effect sizes. Cortical analyses identified unique differences in each disorder, but also considerable overlap between the two, specifically in cortical thickness. Ongoing work is examining alternative research questions, such as brain laterality, prediction of case-control status, and anatomical heterogeneity. In brief, great strides have been made toward fulfilling the aims of the ENIGMA collaborations, while new ideas and follow-up analyses continue that include more imaging modalities (diffusion MRI and resting-state functional MRI), collaborations with other large databases, and samples with dual diagnoses.

Localization of anatomical changes in patients during proton therapy with in-beam PET monitoring: A voxel-based morphometry approach exploiting Monte Carlo simulations

Purpose

In‐beam positron emission tomography (PET) is one of the modalities that can be used for in vivo noninvasive treatment monitoring in proton therapy. Although PET monitoring has been frequently applied for this purpose, there is still no straightforward method to translate the information obtained from the PET images into easy‐to‐interpret information for clinical personnel. The purpose of this work is to propose a statistical method for analyzing in‐beam PET monitoring images that can be used to locate, quantify, and visualize regions with possible morphological changes occurring over the course of treatment.

Methods

We selected a patient treated for squamous cell carcinoma (SCC) with proton therapy, to perform multiple Monte Carlo (MC) simulations of the expected PET signal at the start of treatment, and to study how the PET signal may change along the treatment course due to morphological changes. We performed voxel‐wise two‐tailed statistical tests of the simulated PET images, resembling the voxel‐based morphometry (VBM) method commonly used in neuroimaging data analysis, to locate regions with significant morphological changes and to quantify the change.

Results

The VBM resembling method has been successfully applied to the simulated in‐beam PET images, despite the fact that such images suffer from image artifacts and limited statistics. Three dimensional probability maps were obtained, that allowed to identify interfractional morphological changes and to visualize them superimposed on the computed tomography (CT) scan. In particular, the characteristic color patterns resulting from the two‐tailed statistical tests lend themselves to trigger alarms in case of morphological changes along the course of treatment.

Conclusions

The statistical method presented in this work is a promising method to apply to PET monitoring data to reveal interfractional morphological changes in patients, occurring over the course of treatment. Based on simulated in‐beam PET treatment monitoring images, we showed that with our method it was possible to correctly identify the regions that changed. Moreover we could quantify the changes, and visualize them superimposed on the CT scan. The proposed method can possibly help clinical personnel in the replanning procedure in adaptive proton therapy treatments.

Enhancing the impact of Artificial Intelligence in Medicine: A joint AIFM-INFN Italian initiative for a dedicated cloud-based computing infrastructure

Artificial Intelligence (AI) techniques have been implemented in the field of Medical Imaging for more than forty years. Medical Physicists, Clinicians and Computer Scientists have been collaborating since the beginning to realize software solutions to enhance the informative content of medical images, including AI-based support systems for image interpretation. Despite the recent massive progress in this field due to the current emphasis on Radiomics, Machine Learning and Deep Learning, there are still some barriers to overcome before these tools are fully integrated into the clinical workflows to finally enable a precision medicine approach to patients' care. Nowadays, as Medical Imaging has entered the Big Data era, innovative solutions to efficiently deal with huge amounts of data and to exploit large and distributed computing resources are urgently needed. In the framework of a collaboration agreement between the Italian Association of Medical Physicists (AIFM) and the National Institute for Nuclear Physics (INFN), we propose a model of an intensive computing infrastructure, especially suited for training AI models, equipped with secure storage systems, compliant with data protection regulation, which will accelerate the development and extensive validation of AI-based solutions in the Medical Imaging field of research. This solution can be developed and made operational by Physicists and Computer Scientists working on complementary fields of research in Physics, such as High Energy Physics and Medical Physics, who have all the necessary skills to tailor the AI-technology to the needs of the Medical Imaging community and to shorten the pathway towards the clinical applicability of AI-based decision support systems.

Quantification of pulmonary involvement in COVID-19 pneumonia by means of a cascade of two U-nets: training and assessment on multiple datasets using different annotation criteria

Purpose

This study aims at exploiting artificial intelligence (AI) for the identification, segmentation and quantification of COVID-19 pulmonary lesions. The limited data availability and the annotation quality are relevant factors in training AI-methods. We investigated the effects of using multiple datasets, heterogeneously populated and annotated according to different criteria.

Methods

We developed an automated analysis pipeline, the LungQuant system, based on a cascade of two U-nets. The first one (U-net\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_1$$\end{document}1) is devoted to the identification of the lung parenchyma; the second one (U-net\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2) acts on a bounding box enclosing the segmented lungs to identify the areas affected by COVID-19 lesions. Different public datasets were used to train the U-nets and to evaluate their segmentation performances, which have been quantified in terms of the Dice Similarity Coefficients. The accuracy in predicting the CT-Severity Score (CT-SS) of the LungQuant system has been also evaluated.

Results

Both the volumetric DSC (vDSC) and the accuracy showed a dependency on the annotation quality of the released data samples. On an independent dataset (COVID-19-CT-Seg), both the vDSC and the surface DSC (sDSC) were measured between the masks predicted by LungQuant system and the reference ones. The vDSC (sDSC) values of 0.95±0.01 and 0.66±0.13 (0.95±0.02 and 0.76±0.18, with 5 mm tolerance) were obtained for the segmentation of lungs and COVID-19 lesions, respectively. The system achieved an accuracy of 90% in CT-SS identification on this benchmark dataset.

Conclusion

We analysed the impact of using data samples with different annotation criteria in training an AI-based quantification system for pulmonary involvement in COVID-19 pneumonia. In terms of vDSC measures, the U-net segmentation strongly depends on the quality of the lesion annotations. Nevertheless, the CT-SS can be accurately predicted on independent test sets, demonstrating the satisfactory generalization ability of the LungQuant.

Supplementary Information

The online version supplementary material available at 10.1007/s11548-021-02501-2.

Strategies to develop radiomics and machine learning models for lung cancer stage and histology prediction using small data samples

Predictive models based on radiomics and machine-learning (ML) need large and annotated datasets for training, often difficult to collect. We designed an operative pipeline for model training to exploit data already available to the scientific community. The aim of this work was to explore the capability of radiomic features in predicting tumor histology and stage in patients with non-small cell lung cancer (NSCLC). We analyzed the radiotherapy planning thoracic CT scans of a proprietary sample of 47 subjects (L-RT) and integrated this dataset with a publicly available set of 130 patients from the MAASTRO NSCLC collection (Lung1). We implemented intra- and inter-sample cross-validation strategies (CV) for evaluating the ML predictive model performances with not so large datasets. We carried out two classification tasks: histology classification (3 classes) and overall stage classification (two classes: stage I and II). In the first task, the best performance was obtained by a Random Forest classifier, once the analysis has been restricted to stage I and II tumors of the Lung1 and L-RT merged dataset (AUC = 0.72 ± 0.11). For the overall stage classification, the best results were obtained when training on Lung1 and testing of L-RT dataset (AUC = 0.72 ± 0.04 for Random Forest and AUC = 0.84 ± 0.03 for linear-kernel Support Vector Machine). According to the classification task to be accomplished and to the heterogeneity of the available dataset(s), different CV strategies have to be explored and compared to make a robust assessment of the potential of a predictive model based on radiomics and ML.

Three dimensional MRF obtains highly repeatable and reproducible multi-parametric estimations in the healthy human brain at 1.5T and 3T

Magnetic resonance fingerprinting (MRF) is highly promising as a quantitative MRI technique due to its accuracy, robustness, and efficiency. Previous studies have found high repeatability and reproducibility of 2D MRF acquisitions in the brain. Here, we have extended our investigations to 3D MRF acquisitions covering the whole brain using spiral projection k-space trajectories. Our travelling head study acquired test/retest data from the brains of 12 healthy volunteers and 8 MRI systems (3 systems at 3 T and 5 at 1.5 T, all from a single vendor), using a study design not requiring all subjects to be scanned at all sites. The pulse sequence and reconstruction algorithm were the same for all acquisitions. After registration of the MRF-derived PD T1 and T2 maps to an anatomical atlas, coefficients of variation (CVs) were computed to assess test/retest repeatability and inter-site reproducibility in each voxel, while a General Linear Model (GLM) was used to determine the voxel-wise variability between all confounders, which included test/retest, subject, field strength and site. Our analysis demonstrated a high repeatability (CVs 0.7-1.3% for T1, 2.0-7.8% for T2, 1.4-2.5% for normalized PD) and reproducibility (CVs of 2.0-5.8% for T1, 7.4-10.2% for T2, 5.2-9.2% for normalized PD) in gray and white matter. Both repeatability and reproducibility improved when compared to similar experiments using 2D acquisitions. Three-dimensional MRF obtains highly repeatable and reproducible estimations of T1 and T2, supporting the translation of MRF-based fast quantitative imaging into clinical applications.

Virtual Histology of Cortical Thickness and Shared Neurobiology in 6 Psychiatric Disorders

IMPORTANCE

Large-scale neuroimaging studies have revealed group differences in cortical thickness across many psychiatric disorders. The underlying neurobiology behind these differences is not well understood.

OBJECTIVE

To determine neurobiologic correlates of group differences in cortical thickness between cases and controls in 6 disorders: attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder (BD), major depressive disorder (MDD), obsessive-compulsive disorder (OCD), and schizophrenia.

DESIGN, SETTING, AND PARTICIPANTS

Profiles of group differences in cortical thickness between cases and controls were generated using T1-weighted magnetic resonance images. Similarity between interregional profiles of cell-specific gene expression and those in the group differences in cortical thickness were investigated in each disorder. Next, principal component analysis was used to reveal a shared profile of group difference in thickness across the disorders. Analysis for gene coexpression, clustering, and enrichment for genes associated with these disorders were conducted. Data analysis was conducted between June and December 2019. The analysis included 145 cohorts across 6 psychiatric disorders drawn from the ENIGMA consortium. The numbers of cases and controls in each of the 6 disorders were as follows: ADHD: 1814 and 1602; ASD: 1748 and 1770; BD: 1547 and 3405; MDD: 2658 and 3572; OCD: 2266 and 2007; and schizophrenia: 2688 and 3244.

MAIN OUTCOMES AND MEASURES

Interregional profiles of group difference in cortical thickness between cases and controls.

RESULTS

A total of 12 721 cases and 15 600 controls, ranging from ages 2 to 89 years, were included in this study. Interregional profiles of group differences in cortical thickness for each of the 6 psychiatric disorders were associated with profiles of gene expression specific to pyramidal (CA1) cells, astrocytes (except for BD), and microglia (except for OCD); collectively, gene-expression profiles of the 3 cell types explain between 25% and 54% of variance in interregional profiles of group differences in cortical thickness. Principal component analysis revealed a shared profile of difference in cortical thickness across the 6 disorders (48% variance explained); interregional profile of this principal component 1 was associated with that of the pyramidal-cell gene expression (explaining 56% of interregional variation). Coexpression analyses of these genes revealed 2 clusters: (1) a prenatal cluster enriched with genes involved in neurodevelopmental (axon guidance) processes and (2) a postnatal cluster enriched with genes involved in synaptic activity and plasticity-related processes. These clusters were enriched with genes associated with all 6 psychiatric disorders.

CONCLUSIONS AND RELEVANCE

In this study, shared neurobiologic processes were associated with differences in cortical thickness across multiple psychiatric disorders. These processes implicate a common role of prenatal development and postnatal functioning of the cerebral cortex in these disorders.

Autism Spectrum Disorder and Childhood Apraxia of Speech: Early Language-Related Hallmarks across Structural MRI Study

Autism Spectrum Disorder (ASD) and Childhood Apraxia of Speech (CAS) are developmental disorders with distinct diagnostic criteria and different epidemiology. However, a common genetic background as well as overlapping clinical features between ASD and CAS have been recently reported. To date, brain structural language-related abnormalities have been detected in both the conditions, but no study directly compared young children with ASD, CAS and typical development (TD). In the current work, we aim: (i) to test the hypothesis that ASD and CAS display neurostructural differences in comparison with TD through morphometric Magnetic Resonance Imaging (MRI)-based measures (ASD vs. TD and CAS vs. TD); (ii) to investigate early possible disease-specific brain structural patterns in the two clinical groups (ASD vs. CAS); (iii) to evaluate predictive power of machine-learning (ML) techniques in differentiating the three samples (ASD, CAS, TD). We retrospectively analyzed the T1-weighted brain MRI scans of 68 children (age range: 34-74 months) grouped into three cohorts: (1) 26 children with ASD (mean age ± standard deviation: 56 ± 11 months); (2) 24 children with CAS (57 ± 10 months); (3) 18 children with TD (55 ± 13 months). Furthermore, a ML analysis based on a linear-kernel Support Vector Machine (SVM) was performed. All but one brain structures displayed significant higher volumes in both ASD and CAS children than TD peers. Specifically, ASD alterations involved fronto-temporal regions together with basal ganglia and cerebellum, while CAS alterations are more focused and shifted to frontal regions, suggesting a possible speech-related anomalies distribution. Caudate, superior temporal and hippocampus volumes directly distinguished the two conditions in terms of greater values in ASD compared to CAS. The ML analysis identified significant differences in brain features between ASD and TD children, whereas only some trends in the ML classification capability were detected in CAS as compared to TD peers. Similarly, the MRI structural underpinnings of two clinical groups were not significantly different when evaluated with linear-kernel SVM. Our results may represent the first step towards understanding shared and specific neural substrate in ASD and CAS conditions, which subsequently may contribute to early differential diagnosis and tailoring specific early intervention.

Numerical and Workbench Design of 2.35 T Double-Tuned (¹H/²³Na) Nested RF Birdcage Coils Suitable for Animal Size MRI

The birdcage Radio Frequency (RF) coil is one of the most used configurations in Magnetic Resonance Imaging (MRI) scanners for the detection of the proton (¹H) signal over a large homogeneous volume. More recently, birdcage RF coils have been successfully used also in the field of X-nuclei MRI, where the signal of a second nucleus (e.g. ¹³C, ²³Na, ³¹P, and many others) needs to be detected with high sensitivity and spatial homogeneity. To this purpose several technical solutions have been adopted to design Double Tuned (DT) volume RF coils, including the recent configuration of the nested birdcage RF coils. One of the main problems in the design of DT RF coils is the decoupling between the ¹H and X channels, and a number of solutions have been adopted over the years. In this work, based on numerical and workbench methods, we report the decoupling optimization of DT (¹H/²³Na) nested RF birdcage coils suitable for 2.35 T MRI scanners encompassing an inner Low-Pass (LP) birdcage used for X-nuclei, an outer High-Pass (HP) birdcage for ¹H and an external cylindrical RF shield. We show that a suitable geometrical selection of the two coaxial RF birdcage coils (relative angular orientation, diameters and lengths) and RF shield (diameter, length) allows a significant decoupling optimization. We also provide valuable information about the RF B₁⁺ field homogeneity and efficiency. Our approach was validated both with numerical simulations and workbench testing using DT nested RF coil prototypes.

Subcortical Brain Volume, Regional Cortical Thickness, and Cortical Surface Area Across Disorders: Findings From the ENIGMA ADHD, ASD, and OCD Working Groups

OBJECTIVE

Attention deficit hyperactivity disorder (ADHD), autism spectrum disorder (ASD), and obsessive-compulsive disorder (OCD) are common neurodevelopmental disorders that frequently co-occur. The authors sought to directly compare these disorders using structural brain imaging data from ENIGMA consortium data.

METHODS

Structural T1-weighted whole-brain MRI data from healthy control subjects (N=5,827) and from patients with ADHD (N=2,271), ASD (N=1,777), and OCD (N=2,323) from 151 cohorts worldwide were analyzed using standardized processing protocols. The authors examined subcortical volume, cortical thickness, and cortical surface area differences within a mega-analytical framework, pooling measures extracted from each cohort. Analyses were performed separately for children, adolescents, and adults, using linear mixed-effects models adjusting for age, sex, and site (and intracranial volume for subcortical and surface area measures).

RESULTS

No shared differences were found among all three disorders, and shared differences between any two disorders did not survive correction for multiple comparisons. Children with ADHD compared with those with OCD had smaller hippocampal volumes, possibly influenced by IQ. Children and adolescents with ADHD also had smaller intracranial volume than control subjects and those with OCD or ASD. Adults with ASD showed thicker frontal cortices compared with adult control subjects and other clinical groups. No OCD-specific differences were observed across different age groups and surface area differences among all disorders in childhood and adulthood.

CONCLUSIONS

The study findings suggest robust but subtle differences across different age groups among ADHD, ASD, and OCD. ADHD-specific intracranial volume and hippocampal differences in children and adolescents, and ASD-specific cortical thickness differences in the frontal cortex in adults, support previous work emphasizing structural brain differences in these disorders.

Retrospective rigid motion correction of three-dimensional magnetic resonance fingerprinting of the human brain

PURPOSE

To obtain three-dimensional (3D), quantitative and motion-robust imaging with magnetic resonance fingerprinting (MRF).

METHODS

Our acquisition is based on a 3D spiral projection k-space scheme. We compared different orderings of trajectory interleaves in terms of rigid motion-correction robustness. In all tested orderings, we considered the whole dataset as a sum of 56 segments of 7-s duration, acquired sequentially with the same flip angle schedule. We performed a separate image reconstruction for each segment, producing whole-brain navigators that were aligned to the first segment using normalized correlation. The estimated rigid motion was used to correct the k-space data, and the aligned data were matched with the dictionary to obtain motion-corrected maps.

RESULTS

A significant improvement on the motion-affected maps after motion correction is evident with the suppression of motion artifacts. Correlation with the motionless baseline improved by 20% on average for both T₁ and T₂ estimations after motion correction. In addition, the average motion-induced quantification bias of 70 ms for T₁ and 18 ms for T₂ values was reduced to 12 ms and 6 ms, respectively, improving the reliability of quantitative estimations.

CONCLUSION

We established a method that allows correcting 3D rigid motion on a 7-s timescale during the reconstruction of MRF data using self-navigators, improving the image quality and the quantification robustness.

Lower gray matter volumes of frontal lobes and insula in adolescents with anorexia nervosa restricting type: Findings from a Brain Morphometry Study

BACKGROUND

Brain atrophy in anorexia nervosa (AN) is one of the most marked structural brain changes observed in mental disorders. In this study, we propose a whole brain analysis approach to characterize global and regional cerebral volumes in adolescents with restricting-type anorexia nervosa (AN-r).

METHODS

A total of 48 adolescent females (age range 13-18 years) were enrolled in the study (24 right-handed AN-r in the early stages of the illness and treated in the same clinical setting and 24 age-matched healthy controls [HC]). High-resolution T1-weighted magnetic resonance images were acquired. Cerebral volumes, including the total amounts of gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF) were obtained with the Statistical Parametric Mapping software (SPM8); specific cortical regional volumes were computed by applying an atlas-based cortical parcellation to the SPM8 GM segments. Analysis of variance (ANOVA) was performed to identify any significant between-group differences in global and regional brain volumes.

RESULTS

The analyses revealed reduced total GM volumes (p = 0.02) and increased CSF (p = 0.05) in AN-r, compared with HC. No significant between-group difference was found in WM volumes. At the regional level, significantly lower GM volumes in both frontal lobes (p = 0.006) and in the left insula (p = 0.016) were detected. No significant relationships were found between cerebral volumes and duration of illness, psychiatric comorbidities, psychopharmacological treatment, prepubertal phase, or presence of amenorrhea.

CONCLUSIONS

The topographic distribution of GM reduction in a homogenous group of AN-r involves regions responsible for the emotional and cognitive deficits associated with the illness. These findings are discussed in relation to the roles of the insular cortex and the frontal lobes.

Measuring the effects of confounders in medical supervised classification problems: the Confounding Index (CI)

Over the years, there has been growing interest in using machine learning techniques for biomedical data processing. When tackling these tasks, one needs to bear in mind that biomedical data depends on a variety of characteristics, such as demographic aspects (age, gender, etc.) or the acquisition technology, which might be unrelated with the target of the analysis. In supervised tasks, failing to match the ground truth targets with respect to such characteristics, called confounders, may lead to very misleading estimates of the predictive performance. Many strategies have been proposed to handle confounders, ranging from data selection, to normalization techniques, up to the use of training algorithm for learning with imbalanced data. However, all these solutions require the confounders to be known a priori. To this aim, we introduce a novel index that is able to measure the confounding effect of a data attribute in a bias-agnostic way. This index can be used to quantitatively compare the confounding effects of different variables and to inform correction methods such as normalization procedures or ad-hoc-prepared learning algorithms. The effectiveness of this index is validated on both simulated data and real-world neuroimaging data.

Evaluation of the intra- and inter-method agreement of brain MRI segmentation software packages: A comparison between SPM12 and FreeSurfer v6.0

PURPOSE

The lack of inter-method agreement can produce inconsistent results in neuroimaging studies. We evaluated the intra-method repeatability and the inter-method reproducibility of two widely-used automatic segmentation methods for brain MRI: the FreeSurfer (FS) and the Statistical Parametric Mapping (SPM) software packages.

METHODS

We segmented the gray matter (GM), the white matter (WM) and subcortical structures in test-retest MRI data of healthy volunteers from Kirby-21 and OASIS datasets. We used Pearson's correlation (r), Bland-Altman plot and Dice index to study intra-method repeatability and inter-method reproducibility. In order to test whether different processing methods affect the results of a neuroimaging-based group study, we carried out a statistical comparison between male and female volume measures.

RESULTS

A high correlation was found between test-retest volume measures for both SPM (r in the 0.98-0.99 range) and FS (r in the 0.95-0.99 range). A non-null bias between test-retest FS volumes was detected for GM and WM in the OASIS dataset. The inter-method reproducibility analysis measured volume correlation values in the 0.72-0.98 range and the overlap between the segmented structures assessed by the Dice index was in the 0.76-0.83 range. SPM systematically provided significantly greater GM volumes and lower WM and subcortical volumes with respect to FS. In the male vs. female brain volume comparisons, inconsistencies arose for the OASIS dataset, where the gender-related differences appear subtler with respect to the Kirby dataset.

CONCLUSIONS

The inter-method reproducibility should be evaluated before interpreting the results of neuroimaging studies.

Multi-site repeatability and reproducibility of MR fingerprinting of the healthy brain at 1.5 and 3.0 T

Fully-quantitative MR imaging methods are useful for longitudinal characterization of disease and assessment of treatment efficacy. However, current quantitative MRI protocols have not been widely adopted in the clinic, mostly due to lengthy scan times. Magnetic Resonance Fingerprinting (MRF) is a new technique that can reconstruct multiple parametric maps from a single fast acquisition in the transient state of the MR signal. Due to the relative novelty of this technique, the repeatability and reproducibility of quantitative measurements obtained using MRF has not been extensively studied. Our study acquired test/retest data from the brains of nine healthy volunteers, each scanned on five MRI systems (two at 3.0 T and three at 1.5 T, all from a single vendor) located at two different centers. The pulse sequence and reconstruction algorithm were the same for all acquisitions. After registration of the MRF-derived M0, T1 and T2 maps to an anatomical atlas, coefficients-of-variation (CVs) were computed to assess test/retest repeatability and inter-site reproducibility in each voxel, while a General Linear Model (GLM) was used to determine the voxel-wise variability between all confounders, which included test/retest, subject, field strength and site. Our analysis demonstrated an excellent repeatability (CVs of 2-3% for T1, 5-8% for T2, 3% for normalized-M0) and a good reproducibility (CVs of 3-8% for T1, 8-14% for T2, 5% for normalized-M0) in grey and white matter.

Evaluation of Altered Functional Connections in Male Children With Autism Spectrum Disorders on Multiple-Site Data Optimized With Machine Learning

No univocal and reliable brain-based biomarkers have been detected to date in Autism Spectrum Disorders (ASD). Neuroimaging studies have consistently revealed alterations in brain structure and function of individuals with ASD; however, it remains difficult to ascertain the extent and localization of affected brain networks. In this context, the application of Machine Learning (ML) classification methods to neuroimaging data has the potential to contribute to a better distinction between subjects with ASD and typical development controls (TD). This study is focused on the analysis of resting-state fMRI data of individuals with ASD and matched TD, available within the ABIDE collection. To reduce the multiple sources of heterogeneity that impact on understanding the neural underpinnings of autistic condition, we selected a subgroup of 190 subjects (102 with ASD and 88 TD) according to the following criteria: male children (age range: 6.5-13 years); rs-fMRI data acquired with open eyes; data from the University sites that provided the largest number of scans (KKI, NYU, UCLA, UM). Connectivity values were evaluated as the linear correlation between pairs of time series of brain areas; then, a Linear kernel Support Vector Machine (L-SVM) classification, with an inter-site cross-validation scheme, was carried out. A permutation test was conducted to identify over-connectivity and under-connectivity alterations in the ASD group. The mean L-SVM classification performance, in terms of the area under the ROC curve (AUC), was 0.75 ± 0.05. The highest performance was obtained using data from KKI, NYU and UCLA sites in training and data from UM as testing set (AUC = 0.83). Specifically, stronger functional connectivity (FC) in ASD with respect to TD involve (p < 0.001) the angular gyrus with the precuneus in the right (R) hemisphere, and the R frontal operculum cortex with the pars opercularis of the left (L) inferior frontal gyrus. Weaker connections in ASD group with respect to TD are the intra-hemispheric R temporal fusiform cortex with the R hippocampus, and the L supramarginal gyrus with L planum polare. The results indicate that both under- and over-FC occurred in a selected cohort of ASD children relative to TD controls, and that these functional alterations are spread in different brain networks.

Brainstem enlargement in preschool children with autism: Results from an intermethod agreement study of segmentation algorithms

The intermethod agreement between automated algorithms for brainstem segmentation is investigated, focusing on the potential involvement of this structure in Autism Spectrum Disorders (ASD). Inconsistencies highlighted in previous studies on brainstem in the population with ASD may in part be a result of poor agreement in the extraction of structural features between different methods. A sample of 76 children with ASD and 76 age-, gender-, and intelligence-matched controls was considered. Volumetric analyses were performed using common tools for brain structures segmentation, namely FSL-FIRST, FreeSurfer (FS), and Advanced Normalization Tools (ANTs). For shape analysis SPHARM-MAT was employed. Intermethod agreement was quantified in terms of Pearson correlations between pairs of volumes obtained by the different methods. The degree of overlap between segmented masks was quantified in terms of the Dice index. Both Pearson correlations and Dice indices, showed poor agreement between FSL-FIRST and the other methods (ANTs and FS), which by contrast, yielded Pearson correlations greater than 0.93 and average Dice indices greater than 0.76 when compared with each other. As with volume, shape analyses exhibited discrepancies between segmentation methods, with particular differences noted between FSL-FIRST and the others (ANT and FS), with under- and over-segmentation in specific brainstem regions. These data suggest that research on brain structure alterations should cross-validate findings across multiple methods. We consistently detected an enlargement of brainstem volume in the whole sample and in the male cohort across multiple segmentation methods, a feature particularly driven by the subgroup of children with idiopathic intellectual disability associated with ASD.

Vascular Function Is Improved After an Environmental Enrichment Program

Environmental enrichment may slow cognitive decay possibly acting through an improvement in vascular function. Aim of the study was to assess the effects of a 7-month cognitive, social, and physical training program on cognitive and vascular function in patients with mild cognitive impairment. In a single-center, randomized, parallel-group study, 113 patients (age, 65–89 years) were randomized to multidomain training (n=55) or usual care (n=58). All participants underwent neuropsychological tests and vascular evaluation, including brachial artery flow-mediated dilation, carotid–femoral pulse wave velocity, carotid distensibility, and assessment of circulating hematopoietic CD34+ and endothelial progenitor cells. At study entry, an age-matched control group (n=45) was also studied. Compared with controls, patients had at study entry a reduced flow-mediated dilation (2.97±2.14% versus 3.73±2.06%; P =0.03) and hyperemic stimulus (shear rate area under the curve, 19.1±15.7 versus 25.7±15.1×10 −3 ; P =0.009); only the latter remained significant after adjustment for confounders ( P =0.03). Training improved Alzheimer disease assessment scale cognitive (training, 14.0±4.8 to 13.1±5.5; nontraining, 12.1±3.9 to 13.2±4.8; P for interaction visit×training=0.02), flow-mediated dilation (2.82±2.19% to 3.40±1.81%, 3.05±2.08% to 2.24±1.59%; P =0.006; P =0.023 after adjustment for diameter and shear rate area under the curve), and circulating hematopoietic CD34 + cells and prevented the decline in carotid distensibility (18.4±5.3 to 20.0±6.6, 23.9±11.0 to 19.5±7.1 Pa −1 ; P =0.005). The only clinical predictor of improvement of cognitive function after training was established hypertension. There was no correlation between changes in measures of cognitive and vascular function. In conclusion, a multidomain training program slows cognitive decline, especially in hypertensive individuals. This effect is accompanied by improved systemic endothelial function, mobilization of progenitor CD34 + cells, and preserved carotid distensibility.

Clinical Trial Registration—

URL: http://www.clinicaltrials.gov . Unique identifier: NCT01725178.

Cortical and Subcortical Brain Morphometry Differences Between Patients With Autism Spectrum Disorder and Healthy Individuals Across the Lifespan: Results From the ENIGMA ASD Working Group

OBJECTIVE

Neuroimaging studies show structural differences in both cortical and subcortical brain regions in children and adults with autism spectrum disorder (ASD) compared with healthy subjects. Findings are inconsistent, however, and it is unclear how differences develop across the lifespan. The authors investigated brain morphometry differences between individuals with ASD and healthy subjects, cross-sectionally across the lifespan, in a large multinational sample from the Enhancing Neuroimaging Genetics Through Meta-Analysis (ENIGMA) ASD working group.

METHOD

The sample comprised 1,571 patients with ASD and 1,651 healthy control subjects (age range, 2-64 years) from 49 participating sites. MRI scans were preprocessed at individual sites with a harmonized protocol based on a validated automated-segmentation software program. Mega-analyses were used to test for case-control differences in subcortical volumes, cortical thickness, and surface area. Development of brain morphometry over the lifespan was modeled using a fractional polynomial approach.

RESULTS

The case-control mega-analysis demonstrated that ASD was associated with smaller subcortical volumes of the pallidum, putamen, amygdala, and nucleus accumbens (effect sizes [Cohen's d], 0.13 to -0.13), as well as increased cortical thickness in the frontal cortex and decreased thickness in the temporal cortex (effect sizes, -0.21 to 0.20). Analyses of age effects indicate that the development of cortical thickness is altered in ASD, with the largest differences occurring around adolescence. No age-by-ASD interactions were observed in the subcortical partitions.

CONCLUSIONS

The ENIGMA ASD working group provides the largest study of brain morphometry differences in ASD to date, using a well-established, validated, publicly available analysis pipeline. ASD patients showed altered morphometry in the cognitive and affective parts of the striatum, frontal cortex, and temporal cortex. Complex developmental trajectories were observed for the different regions, with a developmental peak around adolescence. These findings suggest an interplay in the abnormal development of the striatal, frontal, and temporal regions in ASD across the lifespan.

GPCALMA: a Grid-based Tool for Mammographic Screening

Objectives: The next generation of high energy physics (HEP) experiments requires a GRID approach to a distributed computing system: the key concept is the Virtual Organisation (VO), a group of distributed users with a common goal and the will to share their resources.

Methods: A similar approach, applied to a group of hospitals that joined the GPCALMA project (Grid Platform for Computer Assisted Library for MAmmography), will allow common screening programs for early diagnosis of breast and, in the future, lung cancer. The application code makes use of neural networks for the image analysis and is useful in improving the radiologists' diagnostic performance. GRID services allow remote image analysis and interactive online diagnosis, with a potential for a relevant reduction of the delays presently associated with screening programs.

Results and Conclusions: A prototype of the system, based on AliEn GRID Services [1], is already available, with a central server running common services [2] and several clients connecting to it. Mammograms can be acquired in any location; the related informatio required to select and access them at any time is stored in a common service called Data Catalogue, which can be queried by any client. Thanks to the PROOF facility [3], the result of a query can be used as input for analy-sis algorithms, which are executed on the nodes where the input images are stored,. The selected approach avoids data transfers for all the images with a negative diagnosis and allows an almost real time diagnosis for the set of images with high cancer probability.

Semiautomated Evaluation of the Primary Motor Cortex in Patients with Amyotrophic Lateral Sclerosis at 3T

BACKGROUND AND PURPOSE

Amyotrophic lateral sclerosis is a neurodegenerative disease involving the upper and lower motor neurons. In amyotrophic lateral sclerosis, pathologic changes in the primary motor cortex include Betz cell depletion and the presence of reactive iron-loaded microglia, detectable on 7T MR images as atrophy and T2*-hypointensity. Our purposes were the following: 1) to investigate the signal hypointensity-to-thickness ratio of the primary motor cortex as a radiologic marker of upper motor neuron involvement in amyotrophic lateral sclerosis with a semiautomated method at 3T, 2) to compare 3T and 7T results, and 3) to evaluate whether semiautomated measurement outperforms visual image assessment.

MATERIALS AND METHODS

We investigated 27 patients and 13 healthy subjects at 3T, and 19 patients and 18 healthy subjects at 7T, performing a high-resolution 3D multiecho T2*-weighted sequence targeting the primary motor cortex. The signal hypointensity-to-thickness ratio of the primary motor cortex was calculated with a semiautomated method depicting signal intensity profiles of the cortex. Images were also visually classified as "pathologic" or "nonpathologic" based on the primary motor cortex signal intensity and thickness.

RESULTS

The signal hypointensity-to-thickness ratio of the primary motor cortex was greater in patients than in controls (P < .001), and it correlated with upper motor neuron impairment in patients (ρ = 0.57, P < .001). The diagnostic accuracy of the signal hypointensity-to-thickness ratio was high at 3T (area under the curve = 0.89) and even higher at 7T (area under the curve = 0.94). The sensitivity of the semiautomated method (0.81) outperformed the sensitivity of the visual assessment (0.56-0.63) at 3T.

CONCLUSIONS

The signal hypointensity-to-thickness ratio of the primary motor cortex calculated with a semiautomated method is suggested as a radiologic marker of upper motor neuron burden in patients with amyotrophic lateral sclerosis. This semiautomated method may be useful for improving the subjective radiologic evaluation of upper motor neuron pathology in patients suspected of having amyotrophic lateral sclerosis.

ARIANNA: A research environment for neuroimaging studies in autism spectrum disorders

The complexity and heterogeneity of Autism Spectrum Disorders (ASD) require the implementation of dedicated analysis techniques to obtain the maximum from the interrelationship among many variables that describe affected individuals, spanning from clinical phenotypic characterization and genetic profile to structural and functional brain images. The ARIANNA project has developed a collaborative interdisciplinary research environment that is easily accessible to the community of researchers working on ASD (https://arianna.pi.infn.it). The main goals of the project are: to analyze neuroimaging data acquired in multiple sites with multivariate approaches based on machine learning; to detect structural and functional brain characteristics that allow the distinguishing of individuals with ASD from control subjects; to identify neuroimaging-based criteria to stratify the population with ASD to support the future development of personalized treatments. Secure data handling and storage are guaranteed within the project, as well as the access to fast grid/cloud-based computational resources. This paper outlines the web-based architecture, the computing infrastructure and the collaborative analysis workflows at the basis of the ARIANNA interdisciplinary working environment. It also demonstrates the full functionality of the research platform. The availability of this innovative working environment for analyzing clinical and neuroimaging information of individuals with ASD is expected to support researchers in disentangling complex data thus facilitating their interpretation.

Hippocampal subfields at ultra high field MRI: An overview of segmentation and measurement methods

The hippocampus is one of the most interesting and studied brain regions because of its involvement in memory functions and its vulnerability in pathological conditions, such as neurodegenerative processes. In the recent years, the increasing availability of Magnetic Resonance Imaging (MRI) scanners that operate at ultra‐high field (UHF), that is, with static magnetic field strength ≥7T, has opened new research perspectives. Compared to conventional high‐field scanners, these systems can provide new contrasts, increased signal‐to‐noise ratio and higher spatial resolution, thus they may improve the visualization of very small structures of the brain, such as the hippocampal subfields. Studying the morphometry of the hippocampus is crucial in neuroimaging research because changes in volume and thickness of hippocampal subregions may be relevant in the early assessment of pathological cognitive decline and Alzheimer's Disease (AD). The present review provides an overview of the manual, semi‐automated and fully automated methods that allow the assessment of hippocampal subfield morphometry at UHF MRI, focusing on the different hippocampal segmentation produced. © 2017 The Authors Hippocampus Published by Wiley Periodicals, Inc.

The effect of age, sex and clinical features on the volume of Corpus Callosum in pre-schoolers with Autism Spectrum Disorder: a case-control study

A growing body of literature has identified volume alterations of the corpus callosum (CC) in subjects with autism spectrum disorders (ASD). However, to date very few investigations have been conducted on pre-school-age ASD children. This study aims to compare the volume of CC and its sub-regions between pre-schoolers with ASD and controls (CON) and to examine their relationship to demographic and clinical variables (sex, age, non-verbal IQ -NVIQ-, expressive non-echolalic language, emotional and behavioural problems, and autism severity). The volume of CC of 40 pre-schoolers with ASD (20 males and 20 females; mean age: 49 ± 12 months; mean NVIQ: 73 ± 22) and 40 sex-, age-, and NVIQ-matched CON subjects (20 M and 20 F; mean age: 49 ± 14 months; mean NVIQ: 73 ± 23) were quantified applying the FreeSurfer automated parcellation software on Magnetic Resonance images. No significant volumetric differences in CC total volume and in its sub-regions between ASD and CON were found using total brain volume as a covariate. Analogously, absence of CC volumetric differences was evident when boys and girls with ASD were compared with their matched controls. The CC total volume of younger ASD male subjects was found significantly larger with respect to matched CON, which is consistent with the atypical growth trajectory widely reported in these young children. The CC total volume was negatively correlated with autism severity, whereas no association between CC volume and other clinical variables was detected. If replicated, the indirect relationship between CC volume and autism severity suggests the involvement of CC in core ASD symptoms.

One-Class Support Vector Machines Identify the Language and Default Mode Regions As Common Patterns of Structural Alterations in Young Children with Autism Spectrum Disorders

The identification of reliable brain endophenotypes of autism spectrum disorders (ASD) has been hampered to date by the heterogeneity in the neuroanatomical abnormalities detected in this condition. To handle the complexity of neuroimaging data and to convert brain images in informative biomarkers of pathology, multivariate analysis techniques based on Support Vector Machines (SVM) have been widely used in several disease conditions. They are usually trained to distinguish patients from healthy control subjects by making a binary classification. Here, we propose the use of the One-Class Classification (OCC) or Data Description method that, in contrast to two-class classification, is based on a description of one class of objects only. This approach, by defining a multivariate normative rule on one class of subjects, allows recognizing examples from a different category as outliers. We applied the OCC to 314 regional features extracted from brain structural Magnetic Resonance Imaging (MRI) scans of young children with ASD (21 males and 20 females) and control subjects (20 males and 20 females), matched on age [range: 22-72 months of age; mean = 49 months] and non-verbal intelligence quotient (NVIQ) [range: 31-123; mean = 73]. We demonstrated that a common pattern of features characterize the ASD population. The OCC SVM trained on the group of ASD subjects showed the following performances in the ASD vs. controls separation: the area under the receiver operating characteristic curve (AUC) was 0.74 for the male and 0.68 for the female population, respectively. Notably, the ASD vs. controls discrimination results were maximized when evaluated on the subsamples of subjects with NVIQ ≥ 70, leading to AUC = 0.81 for the male and AUC = 0.72 for the female populations, respectively. Language regions and regions from the default mode network-posterior cingulate cortex, pars opercularis and pars triangularis of the inferior frontal gyrus, and transverse temporal gyrus-contributed most to distinguishing individuals with ASD from controls, arguing for the crucial role of these areas in the ASD pathophysiology. The observed brain patterns associate preschoolers with ASD independently of their age, gender and NVIQ and therefore they are expected to constitute part of the ASD brain endophenotype.

Rehabilitative Interventions and Brain Plasticity in Autism Spectrum Disorders: Focus on MRI-Based Studies

Clinical and research evidence supports the efficacy of rehabilitative intervention for improving targeted skills or global outcomes in individuals with autism spectrum disorder (ASD). However, putative mechanisms of structural and functional brain changes are poorly understood. This review aims to investigate the research literature on the neural circuit modifications after non-pharmacological intervention. For this purpose, longitudinal studies that used magnetic resonance imaging (MRI)-based techniques at the start and at the end of the trial to evaluate the neural effects of rehabilitative treatment in subjects with ASD were identified. The six included studies involved a limited number of patients in the active group (from 2 to 16), and differed by acquisition method (task-related and resting-state functional MRI) as well as by functional MRI tasks. Overall, the results produced by the selected investigations demonstrated brain plasticity during the treatment interval that results in an activation/functional connectivity more similar to those of subjects with typical development (TD). Repeated MRI evaluation may represent a promising tool for the detection of neural changes in response to treatment in patients with ASD. However, large-scale randomized controlled trials after standardized rehabilitative intervention are required before translating these preliminary results into clinical use.

The effect of gender on the neuroanatomy of children with autism spectrum disorders: a support vector machine case-control study

Background

Genetic, hormonal, and environmental factors contribute since infancy to sexual dimorphism in regional brain structures of subjects with typical development. However, the neuroanatomical differences between male and female children with autism spectrum disorders (ASD) are an intriguing and still poorly investigated issue. This study aims to evaluate whether the brain of young children with ASD exhibits sex-related structural differences and if a correlation exists between clinical ASD features and neuroanatomical underpinnings.

Methods

A total of 152 structural MRI scans were analysed. Specifically, 76 young children with ASD (38 males and 38 females; 2–7 years of age; mean = 53 months, standard deviation = 17 months) were evaluated employing a support vector machine (SVM)-based analysis of the grey matter (GM). Group comparisons consisted of 76 age-, gender- and non-verbal-intelligence quotient-matched children with typical development or idiopathic developmental delay without autism.

Results

For both genders combined, SVM showed a significantly increased GM volume in young children with ASD with respect to control subjects, predominantly in the bilateral superior frontal gyrus (Brodmann area –BA– 10), bilateral precuneus (BA 31), bilateral superior temporal gyrus (BA 20/22), whereas less GM in patients with ASD was found in right inferior temporal gyrus (BA 37). For the within gender comparisons (i.e., females with ASD vs. controls and males with ASD vs. controls), two overlapping regions in bilateral precuneus (BA 31) and left superior frontal gyrus (BA 9/10) were detected. Sex-by-group analyses revealed in males with ASD compared to matched controls two male-specific regions of increased GM volume (left middle occipital gyrus—BA 19—and right superior temporal gyrus—BA 22). Comparisons in females with and without ASD demonstrated increased GM volumes predominantly in the bilateral frontal regions. Additional regions of significantly increased GM volume in the right anterior cingulate cortex (BA 32) and right cerebellum were typical only of females with ASD.

Conclusions

Despite the specific behavioural correlates of sex-dimorphism in ASD, brain morphology as yet remains unclear and requires future dedicated investigations. This study provides evidence of structural brain gender differences in young children with ASD that possibly contribute to the different phenotypic disease manifestations in males and females.

Electronic supplementary material

The online version of this article (doi:10.1186/s13229-015-0067-3) contains supplementary material, which is available to authorized users.

Temporal lobe connects regression and macrocephaly to autism spectrum disorders

Interictal electroencephalogram (EEG) abnormalities are frequently associated with autism spectrum disorders (ASD), although their relationship with the clinical features of ASD, particularly the regressive onset, remains controversial. The aim of this study was to investigate whether the characteristics of interictal EEG abnormalities might help to distinguish and predict definite phenotypes within the heterogeneity of ASD. We reviewed the awake and sleep interictal EEGs of 220 individuals with idiopathic ASD, either with or without a history of seizures. EEG findings were analyzed with respect to a set of clinical variables to explore significant associations. A brain morphometry study was also carried out on a subgroup of patients. EEG abnormalities were seen in 154/220 individuals (70%) and were mostly focal (p < 0.01) with an anterior localization (p < 0.001). They were detected more frequently during sleep (p < 0.01), and were associated with a regressive onset of ASD (p < 0.05), particularly in individuals with focal temporal localization (p < 0.05). This association was also stronger in regressive patients with concurrent macrocephaly, together with a relative volumetric reduction of the right temporal cortex (p < 0.05). Indeed, concurrence of temporal EEG abnormalities, regression and macrocephaly might possibly define a distinct endophenotype of ASD. EEG-based endophenotypes could be useful to untangle the complexity of ASD, helping to establish anatomic or pathophysiologic subtypes of the disorder.

Evaluation of 3D radio-frequency electromagnetic fields for any matching and coupling conditions by the use of basis functions

A procedure for evaluating radio-frequency electromagnetic fields in anatomical human models for any matching and coupling conditions is introduced. The procedure resorts to the extraction of basis functions: such basis functions, which represent the fields produced by each individual port without any residual coupling, are derived through an algebraic procedure which uses the S parameter matrix and the fields calculated in one (only) full-wave simulation. The basis functions are then used as building-blocks for calculating the fields for any other S parameter matrix. The proposed approach can be used both for volume coil driven in quadrature and for parallel transmission configuration.

Gray Matter Alterations in Young Children with Autism Spectrum Disorders: Comparing Morphometry at the Voxel and Regional Level

BACKGROUND AND PURPOSE

Sophisticated algorithms to infer disease diagnosis, pathology progression and patient outcome are increasingly being developed to analyze brain MRI data. They have been successfully implemented in a variety of diseases and are currently investigated in the field of neuropsychiatric disorders, including autism spectrum disorder (ASD). We aim to test the ability to predict ASD from subtle morphological changes in structural magnetic resonance imaging (sMRI).

METHODS

The analysis of sMRI of a cohort of male ASD children and controls matched for age and nonverbal intelligence quotient (NVIQ) has been carried out with two widely used preprocessing software packages (SPM and Freesurfer) to extract brain morphometric information at different spatial scales. Then, support vector machines have been implemented to classify the brain features and to localize which brain regions contribute most to the ASD-control separation.

RESULTS

The features extracted from the gray matter subregions provide the best classification performance, reaching an area under the receiver operating characteristic curve (AUC) of 74%. This value is enhanced to 80% when considering only subjects with NVIQ over 70.

CONCLUSIONS

Despite the subtle impact of ASD on brain morphology and a limited cohort size, results from sMRI-based classifiers suggest a consistent network of altered brain regions.