Visual ensemble selection of deep convolutional neural networks for 3D segmentation of breast tumors on dynamic contrast enhanced MRI.

OBJECTIVES: To develop a visual ensemble selection of deep convolutional neural networks (CNN) for 3D segmentation of breast tumors using T1-weighted dynamic contrast-enhanced (T1-DCE) MRI. METHODS: Multi-center 3D T1-DCE MRI (n = 141) were acquired for a cohort of patients diagnosed with locally advanced or aggressive breast cancer. Tumor lesions of 111 scans were equally divided between two radiologists and segmented for training. The additional 30 scans were segmented independently by both radiologists for testing. Three 3D U-Net models were trained using either post-contrast images or a combination of post-contrast and subtraction images fused at either the image or the feature level. Segmentation accuracy was evaluated quantitatively using the Dice similarity coefficient (DSC) and the Hausdorff distance (HD95) and scored qualitatively by a radiologist as excellent, useful, helpful, or unacceptable. Based on this score, a visual ensemble approach selecting the best segmentation among these three models was proposed. RESULTS: The mean and standard deviation of DSC and HD95 between the two radiologists were equal to 77.8 ± 10.0% and 5.2 ± 5.9 mm. Using the visual ensemble selection, a DSC and HD95 equal to 78.1 ± 16.2% and 14.1 ± 40.8 mm was reached. The qualitative assessment was excellent (resp. excellent or useful) in 50% (resp. 77%). CONCLUSION: Using subtraction images in addition to post-contrast images provided complementary information for 3D segmentation of breast lesions by CNN. A visual ensemble selection allowing the radiologist to select the most optimal segmentation obtained by the three 3D U-Net models achieved comparable results to inter-radiologist agreement, yielding 77% segmented volumes considered excellent or useful. KEY POINTS: • Deep convolutional neural networks were developed using T1-weighted post-contrast and subtraction MRI to perform automated 3D segmentation of breast tumors. • A visual ensemble selection allowing the radiologist to choose the best segmentation among the three 3D U-Net models outperformed each of the three models. • The visual ensemble selection provided clinically useful segmentations in 77% of cases, potentially allowing for a valuable reduction of the manual 3D segmentation workload for the radiologist and greatly facilitating quantitative studies on non-invasive biomarker in breast MRI.

Exploring the added value of pretherapeutic MR descriptors in predicting breast cancer pathologic complete response to neoadjuvant chemotherapy.

OBJECTIVES: To evaluate the association between pretreatment MRI descriptors and breast cancer (BC) pathological complete response (pCR) to neoadjuvant chemotherapy (NAC). MATERIALS AND METHODS: Patients with BC treated by NAC with a breast MRI between 2016 and 2020 were included in this retrospective observational single-center study. MR studies were described using the standardized BI-RADS and breast edema score on T2-weighted MRI. Univariable and multivariable logistic regression analyses were performed to assess variables association with pCR according to residual cancer burden. Random forest classifiers were trained to predict pCR on a random split including 70% of the database and were validated on the remaining cases. RESULTS: Among 129 BC, 59 (46%) achieved pCR after NAC (luminal (n = 7/37, 19%), triple negative (n = 30/55, 55%), HER2 + (n = 22/37, 59%)). Clinical and biological items associated with pCR were BC subtype (p < 0.001), T stage 0/I/II (p = 0.008), higher Ki67 (p = 0.005), and higher tumor-infiltrating lymphocytes levels (p = 0.016). Univariate analysis showed that the following MRI features, oval or round shape (p = 0.047), unifocality (p = 0.026), non-spiculated margins (p = 0.018), no associated non-mass enhancement (p = 0.024), and a lower MRI size (p = 0.031), were significantly associated with pCR. Unifocality and non-spiculated margins remained independently associated with pCR at multivariable analysis. Adding significant MRI features to clinicobiological variables in random forest classifiers significantly increased sensitivity (0.67 versus 0.62), specificity (0.69 versus 0.67), and precision (0.71 versus 0.67) for pCR prediction. CONCLUSION: Non-spiculated margins and unifocality are independently associated with pCR and can increase models performance to predict BC response to NAC. CLINICAL RELEVANCE STATEMENT: A multimodal approach integrating pretreatment MRI features with clinicobiological predictors, including tumor-infiltrating lymphocytes, could be employed to develop machine learning models for identifying patients at risk of non-response. This may enable consideration of alternative therapeutic strategies to optimize treatment outcomes. KEY POINTS: • Unifocality and non-spiculated margins are independently associated with pCR at multivariable logistic regression analysis. • Breast edema score is associated with MR tumor size and TIL expression, not only in TN BC as previously reported, but also in luminal BC. • Adding significant MRI features to clinicobiological variables in machine learning classifiers significantly increased sensitivity, specificity, and precision for pCR prediction.

Preclinical and clinical evaluation of a new method to assess cardiac insulin resistance using nuclear imaging.

BACKGROUND: Myocardial insulin resistance (IR) could be a predictive factor of cardiovascular events. This study aimed to introduce a new method using 123I-6-deoxy-6-iodo-D-glucose (6DIG), a pure tracer of glucose transport, for the assessment of IR using cardiac dynamic nuclear imaging. METHODS: The protocol evaluated first in rat-models consisted in two 6DIG injections and one of insulin associated with planar imaging and blood sampling. Compartmental modeling was used to analyze 6DIG kinetics in basal and insulin conditions and to obtain an index of IR. As a part of a translational approach, a clinical study was then performed in 5 healthy and 6 diabetic volunteers. RESULTS: In rodent models, the method revealed reproducible when performed twice at 7 days apart in the same animal. Rosiglitazone, an insulin-sensitizing drug, induced a significant increase of myocardial IR index in obese Zucker rats from 0.96 ± 0.18 to 2.26 ± 0.44 (P<.05) after 7 days of an oral treatment, and 6DIG IR indexes correlated with the gold standard IR index obtained through the hyperinsulinemic-euglycemic clamp (r=.68, P<.02). In human, a factorial analysis was applied on images to obtain vascular and myocardial kinetics before compartmental modeling. 1.5-fold to 2.2-fold decreases in mean cardiac IR indexes from healthy to diabetic volunteers were observed without reaching statistical significance. CONCLUSIONS: These preclinical results demonstrate the reproducibility and sensibility of this novel imaging methodology. Although this first in-human study showed that this new method could be rapidly performed, larger studies need to be planned in order to confirm its performance.

How can we combat multicenter variability in MR radiomics? Validation of a correction procedure.

OBJECTIVE: Test a practical realignment approach to compensate the technical variability of MR radiomic features. METHODS: T1 phantom images acquired on 2 scanners, FLAIR and contrast-enhanced T1-weighted (CE-T1w) images of 18 brain tumor patients scanned on both 1.5-T and 3-T scanners, and 36 T2-weighted (T2w) images of prostate cancer patients scanned in one of two centers were investigated. The ComBat procedure was used for harmonizing radiomic features. Differences in statistical distributions in feature values between 1.5- and 3-T images were tested before and after harmonization. The prostate studies were used to determine the impact of harmonization to distinguish between Gleason grades (GGs). RESULTS: In the phantom data, 40 out of 42 radiomic feature values were significantly different between the 2 scanners before harmonization and none after. In white matter regions, the statistical distributions of features were significantly different (p < 0.05) between the 1.5- and 3-T images for 37 out of 42 features in both FLAIR and CE-T1w images. After harmonization, no statistically significant differences were observed. In brain tumors, 41 (FLAIR) or 36 (CE-T1w) out of 42 features were significantly different between the 1.5- and 3-T images without harmonization, against 1 (FLAIR) or none (CE-T1w) with harmonization. In prostate studies, 636 radiomic features were significantly different between GGs after harmonization against 461 before. The ability to distinguish between GGs using radiomic features was increased after harmonization. CONCLUSION: ComBat harmonization efficiently removes inter-center technical inconsistencies in radiomic feature values and increases the sensitivity of studies using data from several scanners. KEY POINTS: • Radiomic feature values obtained using different MR scanners or imaging protocols can be harmonized by combining off-the-shelf image standardization and feature realignment procedures. • Harmonized radiomic features enable one to pool data from different scanners and centers without a substantial loss of statistical power caused by intra- and inter-center variability. • The proposed realignment method is applicable to radiomic features from different MR sequences and tumor types and does not rely on any phantom acquisition.

A radiomics pipeline dedicated to Breast MRI: validation on a multi-scanner phantom study.

OBJECTIVE: Quantitative analysis in MRI is challenging due to variabilities in intensity distributions across patients, acquisitions and scanners and suffers from bias field inhomogeneity. Radiomic studies are impacted by these effects that affect radiomic feature values. This paper describes a dedicated pipeline to increase reproducibility in breast MRI radiomic studies. MATERIALS AND METHODS: T1, T2, and T1-DCE MR images of two breast phantoms were acquired using two scanners and three dual breast coils. Images were retrospectively corrected for bias field inhomogeneity and further normalised using Z score or histogram matching. Extracted radiomic features were harmonised between coils by the ComBat method. The whole pipeline was assessed qualitatively and quantitatively using statistical comparisons on two series of radiomic feature values computed in the gel mimicking the normal breast tissue or in dense lesions. RESULTS: Intra and inter-acquisition variabilities were strongly reduced by the standardisation pipeline. Harmonisation by ComBat lowered the percentage of radiomic features significantly different between the three coils from 87% after bias field correction and MR normalisation to 3% in the gel, while preserving or improving performance of lesion classification in the phantoms. DISCUSSION: A dedicated standardisation pipeline was developed to reduce variabilities in breast MRI, which paves the way for robust multi-scanner radiomic studies but needs to be assessed on patient data.

Validation of A Method to Compensate Multicenter Effects Affecting CT Radiomics.

Background Radiomics extracts features from medical images more precisely and more accurately than visual assessment. However, radiomics features are affected by CT scanner parameters such as reconstruction kernel or section thickness, thus obscuring underlying biologically important texture features. Purpose To investigate whether a compensation method could correct for the variations of radiomic feature values caused by using different CT protocols. Materials and Methods Phantom data involving 10 texture patterns and 74 patients in cohorts 1 (19 men; 42 patients; mean age, 60.4 years; September-October 2013) and 2 (16 men; 32 patients; mean age, 62.1 years; January-September 2007) scanned by using different CT protocols were retrospectively included. For any radiomic feature, the compensation approach identified a protocol-specific transformation to express all data in a common space that were devoid of protocol effects. The differences in statistical distributions between protocols were assessed by using Friedman tests before and after compensation. Principal component analyses were performed on the phantom data to evaluate the ability to distinguish between texture patterns after compensation. Results In the phantom data, the statistical distributions of features were different between protocols for all radiomic features and texture patterns (P < .05). After compensation, the protocol effect was no longer detectable (P > .05). Principal component analysis demonstrated that each texture pattern was no longer displayed as different clusters corresponding to different imaging protocols, unlike what was observed before compensation. The correction for scanner effect was confirmed in patient data with 100% (10 of 10 features for cohort 1) and 98% (87 of 89 features for cohort 2) of P values less than .05 before compensation, compared with 30% (three of 10) and 15% (13 of 89) after compensation. Conclusion Image compensation successfully realigned feature distributions computed from different CT imaging protocols and should facilitate multicenter radiomic studies. © RSNA, 2019 Online supplemental material is available for this article. See also the editorial by Steiger and Sood in this issue.

Transcutaneous Laryngeal Ultrasonography for Laryngeal Immobility Diagnosis in Patients with Voice Disorders After Thyroid/Parathyroid Surgery.

BACKGROUND: Transcutaneous laryngeal ultrasonography (TLUS) was recently developed to assess recurrent nerve palsy after thyroid/parathyroid surgery, with variable rates of efficiency. The aim of the current study was to evaluate this technique using subjective estimation and post-processing quantitative data. METHODS: Fifty subjects presenting with a recurrent nerve palsy and 50 "controls" presenting with voice, swallowing, or breathing disorders following thyroid/parathyroid surgery were prospectively included. All of them underwent a flexible laryngoscopy, considered the gold standard, and a ten-second TLUS clip within the 10 days following surgery. In addition to the subjective interpretation of vocal fold motion, two quantitative criteria taking into account motion symmetry (symmetry index, SI) and amplitude (mobility index) of the two hemi-larynges were defined on TLUS acquisitions in adduction and abduction. RESULTS: The subjective interpretation provided a sensitivity of 100% and a specificity of 96%, compared to the gold standard. The quantitative criteria provided a sensitivity and specificity of both 82%, when based on SI solely. When combining SI and mobility index, the sensitivity reached 94%, but the specificity fell to 66%. CONCLUSIONS: Visual assessment of recurrent nerve palsy using TLUS after thyroid/parathyroid surgery appeared a high sensitive and specific test compared to flexible laryngoscopy. Quantitative criteria are promising and need to be refined to better describe the whole TLUS video clip.

Sonographic Dynamic Description of the Laryngeal Tract: Definition of Quantitative Measures to Characterize Vocal Fold Motion and Estimation of Their Normal Values.

Vocal fold motion was analyzed during free breathing using two-dimensional dynamic ultrasound imaging. Two cadavers were first analyzed to define easily identifiable landmarks. Motion of the laryngeal tract was then analyzed in an axial plane. Left and right arytenoids and thyroid cartilage were defined on images corresponding to abduction and adduction of the laryngeal tract. Associated area measurements were established for 50 healthy subjects. All area indices were significantly larger during abduction than adduction. Symmetry of motion was established by comparing each hemi-larynx, and mobility fractions were defined. Normal values of laryngeal motion during free breathing were thus established.

Multimodal MRI radiomic models to predict genomic mutations in diffuse intrinsic pontine glioma with missing imaging modalities.

Purpose: Predicting H3.1, TP53, and ACVR1 mutations in DIPG could aid in the selection of therapeutic options. The contribution of clinical data and multi-modal MRI were studied for these three predictive tasks. To keep the maximum number of subjects, which is essential for a rare disease, missing data were considered. A multi-modal model was proposed, collecting all available data for each patient, without performing any imputation. Methods: A retrospective cohort of 80 patients with confirmed DIPG and at least one of the four MR modalities (T1w, T1c, T2w, and FLAIR), acquired with two different MR scanners was built. A pipeline including standardization of MR data and extraction of radiomic features within the tumor was applied. The values of radiomic features between the two MR scanners were realigned using the ComBat method. For each prediction task, the most robust features were selected based on a recursive feature elimination with cross-validation. Five different models, one based on clinical data and one per MR modality, were developed using logistic regression classifiers. The prediction of the multi-modal model was defined as the average of all possible prediction results among five for each patient. The performances of the models were compared using a leave-one-out approach. Results: The percentage of missing modalities ranged from 6 to 11% across modalities and tasks. The performance of each individual model was dependent on each specific task, with an AUC of the ROC curve ranging from 0.63 to 0.80. The multi-modal model outperformed the clinical model for each prediction tasks, thus demonstrating the added value of MRI. Furthermore, regardless of performance criteria, the multi-modal model came in the first place or second place (very close to first). In the leave-one-out approach, the prediction of H3.1 (resp. ACVR1 and TP53) mutations achieved a balanced accuracy of 87.8% (resp. 82.1 and 78.3%). Conclusion: Compared with a single modality approach, the multi-modal model combining multiple MRI modalities and clinical features was the most powerful to predict H3.1, ACVR1, and TP53 mutations and provided prediction, even in the case of missing modality. It could be proposed in the absence of a conclusive biopsy.

Decrypting the information captured by MRI-radiomic features in predicting the response to neoadjuvant chemotherapy in breast cancer.

MRI-based radiomic models have shown promises in predicting the response to neoadjuvant chemotherapy in breast cancer. However, it is difficult to determine which information from the images contributes the most to the prediction: the distribution of gray-levels, the tumour heterogeneity, the shape of the lesions or the intensities of peritumoural regions. The purpose of this study is to dissociate the different sources of information to improve prediction results. Based on pre-treatment MR images from 103 patients, four types of 3D Volumes Of Interest were defined and arranged in multiple combinations. Combining features extracted from different regions proved to increase prediction performances. Clinical relevance- This study proposes a method based on analyses of MRI tumor heterogeneity, margins and peritumoral regions to improve the prediction of the response to neoadjuvant chemotherapy in breast cancer, which would help personalize patient treatment.

Automated estimation of aortic strain from steady-state free-precession and phase contrast MR images.

The strain values extracted from steady-state free-precession (SSFP) and phase contrast (PC) images acquired with a 1.5T scanner on a compliant flow phantom and within the thoracic aorta of 52 healthy subjects were compared. Aortic data were acquired perpendicular to the aorta at the level of the pulmonary artery bifurcation. Cross sectional areas were obtained by using an automatic and robust segmentation method. While a good correlation (r = 0.99) was found between the aortic areas extracted from SSFP and PC sequences, a lower correlation (r = 0.71) was found between the corresponding aortic strain values. Strain values estimated using SSFP and PC sequences were equally correlated with age. Interobserver reproducibility was better for SSFP than for PC. Strain values in the ascending and descending aorta were better correlated for SSFP (r = 0.8) than for PC (r = 0.65) and fitted with the expectation of a larger strain in the ascending aorta when using SSFP. The spatial and temporal resolutions of the acquisitions had a minor influence upon the estimated strain values. Thus, if PC acquisitions can be used to estimate both pulse wave velocity and aortic strain, an additional SSFP sequence may be useful to improve the accuracy in estimating the aortic strain.

Impact of 7-ketocholesterol and very long chain fatty acids on oligodendrocyte lipid membrane organization: evaluation via LAURDAN and FAMIS spectral image analysis.

In the context of multiple sclerosis and X-linked adrenoleukodystrophy, 7-ketocholesterol (7KC) and very long chain fatty acids (C24:0, C26:0) are supposed to induce side effects respectively on oligodendrocytes which are myelin (which is a lipoproteic complex) synthesizing cells. The effects of 7KC (25, 50 µM), C24:0 and C26:0 (10, 20 µM) on cell viability and lipid membrane organization were investigated on 158N murine oligodendrocytes. Concerning 7KC and fatty acids (at 20 µM only): 1) cell growth was strongly inhibited; 2) marked induction of cell death was revealed with propidium iodide (PI); 3) no apoptotic cells were found with C24:0 and C26:0 (absence of cells with condensed and/or fragmented nuclei, of FLICA positive cells and of PI negative/SYTO16 negative cells); 4) some apoptotic cells were detected with 7KC. Fatty acids (at 20 µM only) and 7KC also induced a disorganization of lipid membranes revealed with Merocyanine 540. So, to point out the effects of 7KC (25 µM), C24:0 and C26:0 (20 µM) on the lateral organization of lipid membranes, we used LAURDAN, which gives simultaneous information about morphology and phase state of lipid domains: its emission is blue in the ordered lipid phase, green in the disordered lipid phase. To overcome the qualitative filtering settings of blue and green emission colors, data obtained by mono- and bi-photon confocal microscopy were analyzed by spectral analysis. Sequences of emission images were obtained on both mono- and bi-photon confocal microscopes and processed by means of Factor Analysis of Medical Image Sequences (FAMIS), which is a relevant tool to unmix emission spectra and provide pure color images. Only 7KC was capable to induce a green emission with LAURDAN. Thus, at concentrations inducing oligodendrocyte cell death, 7KC (25 µM) is more efficient than C24:0 and C26:0 (20 µM), to trigger lateral lipid membrane disorganization.

Consistency of aortic distensibility and pulse wave velocity estimates with respect to the Bramwell-Hill theoretical model: a cardiovascular magnetic resonance study.

BACKGROUND: Arterial stiffness is considered as an independent predictor of cardiovascular mortality, and is increasingly used in clinical practice. This study aimed at evaluating the consistency of the automated estimation of regional and local aortic stiffness indices from cardiovascular magnetic resonance (CMR) data. RESULTS: Forty-six healthy subjects underwent carotid-femoral pulse wave velocity measurements (CF_PWV) by applanation tonometry and CMR with steady-state free-precession and phase contrast acquisitions at the level of the aortic arch. These data were used for the automated evaluation of the aortic arch pulse wave velocity (Arch_PWV), and the ascending aorta distensibility (AA_Distc, AA_Distb), which were estimated from ascending aorta strain (AA_Strain) combined with either carotid or brachial pulse pressure. The local ascending aorta pulse wave velocity AA_PWVc and AA_PWVb were estimated respectively from these carotid and brachial derived distensibility indices according to the Bramwell-Hill theoretical model, and were compared with the Arch_PWV. In addition, a reproducibility analysis of AA_PWV measurement and its comparison with the standard CF_PWV was performed. Characterization according to the Bramwell-Hill equation resulted in good correlations between Arch_PWV and both local distensibility indices AA_Distc (r = 0.71, p < 0.001) and AA_Distb (r = 0.60, p < 0.001); and between Arch_PWV and both theoretical local indices AA_PWVc (r = 0.78, p < 0.001) and AA_PWVb (r = 0.78, p < 0.001). Furthermore, the Arch_PWV was well related to CF_PWV (r = 0.69, p < 0.001) and its estimation was highly reproducible (inter-operator variability: 7.1%). CONCLUSIONS: The present work confirmed the consistency and robustness of the regional index Arch_PWV and the local indices AA_Distc and AA_Distb according to the theoretical model, as well as to the well established measurement of CF_PWV, demonstrating the relevance of the regional and local CMR indices.

Impact of ComBat and a Multi-Model approach to deal with multi-scanner and missing MRI data in a small cohort study. Application to H3K27M mutation prediction in patients with DIPG.

Radiomics was proposed to identify tumor phenotypes non-invasively from quantitative imaging features. Calculating a large amount of information on images, allows the development of reliable classification models. In multi-modal imaging protocols, the question arises of adding an imaging modality to improve model performance. In addition, in the implementation of clinical protocols, some modalities are not acquired or are of insufficient quality and cannot be reliably taken into account. Furthermore, multi-scanner studies generate some variability in the acquisition and data. Some methodological solutions using ComBat and a multi-model approach were tested to take these two issues into account. It was applied to a cohort of 88 patients with Diffuse Intrinsic Pontine Glioma (DIPG). Sixteen models using radiomic features computed using 0, 1, 2, 3 or 4 MRI modalities were proposed. Based on Leave-One-Out Cross-Validation, F1 weighted scores ranged from 0.66 to 0.85. A model of majority voting using the prediction of all the models available for one given patient was finally applied, reducing drastically the number of unclassified patients.Clinical relevance- In case of patients with DIPG, the prediction of H3 mutation is of prime importance in case of inconclusive biopsy or in the absence of it. It could suggest orientations for new chemotherapy drugs associated with the radiation therapy.

Object Detection Improves Tumour Segmentation in MR Images of Rare Brain Tumours.

Tumour lesion segmentation is a key step to study and characterise cancer from MR neuroradiological images. Presently, numerous deep learning segmentation architectures have been shown to perform well on the specific tumour type they are trained on (e.g., glioblastoma in brain hemispheres). However, a high performing network heavily trained on a given tumour type may perform poorly on a rare tumour type for which no labelled cases allows training or transfer learning. Yet, because some visual similarities exist nevertheless between common and rare tumours, in the lesion and around it, one may split the problem into two steps: object detection and segmentation. For each step, trained networks on common lesions could be used on rare ones following a domain adaptation scheme without extra fine-tuning. This work proposes a resilient tumour lesion delineation strategy, based on the combination of established elementary networks that achieve detection and segmentation. Our strategy allowed us to achieve robust segmentation inference on a rare tumour located in an unseen tumour context region during training. As an example of a rare tumour, Diffuse Intrinsic Pontine Glioma (DIPG), we achieve an average dice score of 0.62 without further training or network architecture adaptation.

Sex moderations in the relationship between aortic stiffness, cognition, and cerebrovascular reactivity in healthy older adults.

It is well established that sex differences exist in the manifestation of vascular diseases. Arterial stiffness (AS) has been associated with changes in cerebrovascular reactivity (CVR) and cognitive decline in aging. Specifically, older adults with increased AS show a decline on executive function (EF) tasks. Interestingly, the relationship between AS and CVR is more complex, where some studies show decreased CVR with increased AS, and others demonstrate preserved CVR despite higher AS. Here, we investigated the possible role of sex on these hemodynamic relationships. Acquisitions were completed in 48 older adults. Pseudo-continuous arterial spin labeling (pCASL) data were collected during a hypercapnia challenge. Aortic pulse wave velocity (PWV) data was acquired using cine phase contrast velocity series. Cognitive function was assessed with a comprehensive neuropsychological battery, and a composite score for EF was calculated using four cognitive tests from the neuropsychological battery. A moderation model test revealed that sex moderated the relationship between PWV and CVR and PWV and EF, but not between CVR and EF. Together, our results indicate that the relationships between central stiffness, cerebral hemodynamics and cognition are in part mediated by sex.

Automated segmentation of the aorta from phase contrast MR images: validation against expert tracing in healthy volunteers and in patients with a dilated aorta.

PURPOSE: To assess if segmentation of the aorta can be accurately achieved using the modulus image of phase contrast (PC) magnetic resonance (MR) acquisitions. MATERIALS AND METHODS: PC image sequences containing both the ascending and descending aorta of 52 subjects were acquired using three different MR scanners. An automated segmentation technique, based on a 2D+t deformable surface that takes into account the features of PC aortic images, such as flow-related effects, was developed. The study was designed to: 1) assess the variability of our approach and its robustness to the type of MR scanner, and 2) determine its sensitivity to aortic dilation and its accuracy against an expert manual tracing. RESULTS: Interobserver variability in the lumen area was 0.59 +/- 0.92% for the automated approach versus 10.09 +/- 8.29% for manual segmentation. The mean Dice overlap measure was 0.945 +/- 0.014. The method was robust to the aortic size and highly correlated (r = 0.99) with the manual tracing in terms of aortic area and diameter. CONCLUSION: A fast and robust automated segmentation of the aortic lumen was developed and successfully tested on images provided by various MR scanners and acquired on healthy volunteers as well as on patients with a dilated aorta.

Correction for Magnetic Field Inhomogeneities and Normalization of Voxel Values Are Needed to Better Reveal the Potential of MR Radiomic Features in Lung Cancer.

Purpose: To design and validate a preprocessing procedure dedicated to T2-weighted MR images of lung cancers so as to improve the ability of radiomic features to distinguish between adenocarcinoma and other histological types. Materials and Methods: A discovery set of 52 patients with advanced lung cancer who underwent T2-weighted MR imaging at 3 Tesla in a single center study from August 2017 to May 2019 was used. Findings were then validated using a validation set of 19 additional patients included from May to October 2019. Tumor type was obtained from the pathology report after trans-thoracic needle biopsy, metastatic lymph node or metastasis samples, or surgical excisions. MR images were preprocessed using N4ITK bias field correction and by normalizing voxel intensities with fat as a reference region. Segmentation and extraction of radiomic features were performed with LIFEx software on the raw images, on the N4ITK-corrected images and on the fully preprocessed images. Two analyses were conducted where radiomic features were extracted: (1) from the whole tumor volume (3D analysis); (2) from all slices encompassing the tumor (2D analysis). Receiver operating characteristic (ROC) analysis was used to identify features that could distinguish between adenocarcinoma and other histological types. Sham experiments were also designed to control the number of false positive findings. Results: There were 31 (12) adenocarcinomas and 21 (7) other histological types in the discovery (validation) set. In 2D, preprocessing increased the number of discriminant radiomic features from 8 without preprocessing to 22 with preprocessing. 2D analysis yielded more features able to identify adenocarcinoma than 3D analysis (12 discriminant radiomic features after preprocessing in 3D). Preprocessing did not increase false positive findings as no discriminant features were identified in any of the sham experiments. The greatest sensitivity of the 2D analysis applied to preprocessed data was confirmed in the validation set. Conclusion: Correction for magnetic field inhomogeneities and normalization of voxel values are essential to reveal the full potential of radiomic features to identify the tumor histological type from MR T2-weighted images, with classification performance similar to those reported in PET/CT and in multiphase CT in lung cancers.

A downsampling strategy to assess the predictive value of radiomic features.

Many studies are devoted to the design of radiomic models for a prediction task. When no effective model is found, it is often difficult to know whether the radiomic features do not include information relevant to the task or because of insufficient data. We propose a downsampling method to answer that question when considering a classification task into two groups. Using two large patient cohorts, several experimental configurations involving different numbers of patients were created. Univariate or multivariate radiomic models were designed from each configuration. Their performance as reflected by the Youden index (YI) and Area Under the receiver operating characteristic Curve (AUC) was compared to the stable performance obtained with the highest number of patients. A downsampling method is described to predict the YI and AUC achievable with a large number of patients. Using the multivariate models involving machine learning, YI and AUC increased with the number of patients while they decreased for univariate models. The downsampling method better estimated YI and AUC obtained with the largest number of patients than the YI and AUC obtained using the number of available patients and identifies the lack of information relevant to the classification task when no such information exists.

Regularized estimation of contrast agent attenuation to improve the imaging of microbubbles in small animal studies.

Quantitative analysis of tissue perfusion using contrast-enhanced ultrasound is still limited by shadowing, which is caused by inadequate compensation for microbubble contrast agent attenuation. Many previous methods have been developed for attenuation correction in soft tissues. However, no method has been proposed to correct for microbubble attenuation in vivo. In this article, a model to estimate microbubble attenuation is presented, using the time-intensity variation in a highly echogenic distal area without contrast uptake. This model is based on the assumption that a linear relationship holds between local microbubble attenuation and local backscatter. The model was applied to 12 murine renal perfusion studies. Parametric images of microbubble attenuation were generated, corresponding to dynamic contrast agent-specific sequences without shadowing. Contrast uptake kinetics consistent with the physiology were retrieved in all perfused areas. This method therefore proved to be of potential interest in the quantification of tissue perfusion in small animal studies.