Noise reduction by multiple path neural network using Attention mechanisms with an emphasis on robustness against Errors: A pilot study on brain Diffusion-Weighted images.

PURPOSE: In deep learning-based noise reduction, larger networks offer advanced and complex functionality by utilizing its greater degree of freedom, but come with increased unpredictability, raising the potential risk of unforeseen errors. Here, we introduce a novel denoising model for diffusion-weighted images that intentionally limits the network output freedom by incorporating multiple pathways with varying degrees of freedom, with the aim of minimizing the chance of unintended alterations to the input. The purpose of this pilot study is to assess the model's ability to perform effective denoising under the constraints. METHODS: Images from 10 healthy volunteers were used. Key innovations in our model development include: (1) neural network architecture that separated the function for calculating the specific output values from the function for adjusting the calculation for each pixel and (2) training that optimised the network based on both image and secondary obtained diffusion tensor. The generated images were compared with the original ones by measuring the deviation from ground truth images (averaged across eight acquisitions). RESULTS: The generated images demonstrated closer alignment with the ground truth images, both visually and statistically (Q < 0.05), compared to the original images. Furthermore, the advantage of the generated images over the original images was also found in the secondary obtained quantitative parameter maps with significance (Q < 0.05). CONCLUSION: The usefulness of the proposed method was suggested because it was successful in improving both the quality of the generated images and accuracy of the major diffusion parameter maps under the given restrictions.

Deep learning-based coronary computed tomography analysis to predict functionally significant coronary artery stenosis.

Fractional flow reserve derived from coronary CT (FFR-CT) is a noninvasive physiological technique that has shown a good correlation with invasive FFR. However, the use of FFR-CT is restricted by strict application standards, and the diagnostic accuracy of FFR-CT analysis may potentially be decreased by severely calcified coronary arteries because of blooming and beam hardening artifacts. The aim of this study was to evaluate the utility of deep learning (DL)-based coronary computed tomography (CT) data analysis in predicting invasive fractional flow reserve (FFR), especially in cases with severely calcified coronary arteries. We analyzed 184 consecutive cases (241 coronary arteries) which underwent coronary CT and invasive coronary angiography, including invasive FFR, within a three-month period. Mean coronary artery calcium scores were 963 ± 1226. We evaluated and compared the vessel-based diagnostic accuracy of our proposed DL model and a visual assessment to evaluate functionally significant coronary artery stenosis (invasive FFR < 0.80). A deep neural network was trained with consecutive short axial images of coronary arteries on coronary CT. Ninety-one coronary arteries of 89 cases (48%) had FFR-positive functionally significant stenosis. On receiver operating characteristics (ROC) analysis to predict FFR-positive stenosis using the trained DL model, average area under the curve (AUC) of the ROC curve was 0.756, which was superior to the AUC of visual assessment of significant (≥ 70%) coronary artery stenosis on CT (0.574, P = 0.011). The sensitivity, specificity, positive and negative predictive value (PPV and NPV), and accuracy of the DL model and visual assessment for detecting FFR-positive stenosis were 82 and 36%, 68 and 78%, 59 and 48%, 87 and 69%, and 73 and 63%, respectively. Sensitivity and NPV for the prediction of FFR-positive stenosis were significantly higher with our DL model than visual assessment (P = 0.0004, and P = 0.024). DL-based coronary CT data analysis has a higher diagnostic accuracy for functionally significant coronary artery stenosis than visual assessment.

Deep neural network-based synthetic image digital fluoroscopy using digitally reconstructed tomography.

We developed a deep neural network (DNN) to generate X-ray flat panel detector (FPD) images from digitally reconstructed radiographic (DRR) images. FPD and treatment planning CT images were acquired from patients with prostate and head and neck (H&N) malignancies. The DNN parameters were optimized for FPD image synthesis. The synthetic FPD images' features were evaluated to compare to the corresponding ground-truth FPD images using mean absolute error (MAE), peak signal-to-noise ratio (PSNR), and structural similarity index measure (SSIM). The image quality of the synthetic FPD image was also compared with that of the DRR image to understand the performance of our DNN. For the prostate cases, the MAE of the synthetic FPD image was improved (= 0.12 ± 0.02) from that of the input DRR image (= 0.35 ± 0.08). The synthetic FPD image showed higher PSNRs (= 16.81 ± 1.54 dB) than those of the DRR image (= 8.74 ± 1.56 dB), while SSIMs for both images (= 0.69) were almost the same. All metrics for the synthetic FPD images of the H&N cases were improved (MAE 0.08 ± 0.03, PSNR 19.40 ± 2.83 dB, and SSIM 0.80 ± 0.04) compared to those for the DRR image (MAE 0.48 ± 0.11, PSNR 5.74 ± 1.63 dB, and SSIM 0.52 ± 0.09). Our DNN successfully generated FPD images from DRR images. This technique would be useful to increase throughput when images from two different modalities are compared by visual inspection.

Real-time deep neural network-based automatic bowel gas segmentation on X-ray images for particle beam treatment.

Since particle beam distribution is vulnerable to change in bowel gas because of its low density, we developed a deep neural network (DNN) for bowel gas segmentation on X-ray images. We used 6688 image datasets from 209 cases as training data, 736 image datasets from 23 cases as validation data and 102 image datasets from 51 cases as test data (total 283 cases). For the training data, we prepared three types of digitally reconstructed radiographic (DRR) images (all-density, bone and gas) by projecting the treatment planning CT image data. However, the real X-ray images acquired in the treatment room showed low contrast that interfered with manual delineation of bowel gas. Therefore, we used synthetic X-ray images converted from DRR images in addition to real X-ray images.We evaluated DNN segmentation accuracy for the synthetic X-ray images using Intersection over Union, recall, precision, and the Dice coefficient, which measured 0.708 ± 0.208, 0.832 ± 0.170, 0.799 ± 0.191, and 0.807 ± 0.178, respectively. The evaluation metrics for the real X-images were less accurate than those for the synthetic X-ray images (0.408 ± 0237, 0.685 ± 0.326, 0.490 ± 0272, and 0.534 ± 0.271, respectively). Computation time was 29.7 ± 1.3 ms/image. Our DNN appears useful in increasing treatment accuracy in particle beam therapy.

Dual-Energy Computed Tomography-Based Iodine Concentration Estimation for Evaluating Choroidal Malignant Melanoma Response to Treatment: Optimization and Primary Validation.

Contrast-enhanced imaging for choroidal malignant melanoma (CMM) is mostly limited to detecting metastatic tumors, possibly due to difficulties in fixing the eye position. We aimed to (1) validate the appropriateness of estimating iodine concentration based on dual-energy computed tomography (DECT) for CMM and optimize the calculation parameters for estimation, and (2) perform a primary clinical validation by assessing the ability of this technique to show changes in CMM after charged-particle radiation therapy. The accuracy of the optimized estimate (eIC_optimized) was compared to an estimate obtained by commercial software (eIC_commercial) by determining the difference from the ground truth. Then, eIC_optimized, tumor volume, and CT values (80 kVp, 140 kVp, and synthesized 120 kVp) were measured at pre-treatment and 3 months and 1.5−2 years after treatment. The difference from the ground truth was significantly smaller in eIC_optimized than in eIC_commercial (p < 0.01). Tumor volume, CT values, and eIC_optimized all decreased significantly at 1.5−2 years after treatment, but only eIC_commercial showed a significant reduction at 3 months after treatment (p < 0.01). eIC_optimized can quantify contrast enhancement in primary CMM lesions and has high sensitivity for detecting the response to charged-particle radiation therapy, making it potentially useful for treatment monitoring.

Super-resolution generative adversarial networks with static T2*WI-based subject-specific learning to improve spatial difference sensitivity in fMRI activation.

The spatial resolution of fMRI is relatively poor and improvements are needed to indicate more specific locations for functional activities. Here, we propose a novel scheme, called Static T2*WI-based Subject-Specific Super Resolution fMRI (STSS-SRfMRI), to enhance the functional resolution, or ability to discriminate spatially adjacent but functionally different responses, of fMRI. The scheme is based on super-resolution generative adversarial networks (SRGAN) that utilize a T2*-weighted image (T2*WI) dataset as a training reference. The efficacy of the scheme was evaluated through comparison with the activation maps obtained from the raw unpreprocessed functional data (raw fMRI). MRI images were acquired from 30 healthy volunteers using a 3 Tesla scanner. The modified SRGAN reconstructs a high-resolution image series from the original low-resolution fMRI data. For quantitative comparison, several metrics were calculated for both the STSS-SRfMRI and the raw fMRI activation maps. The ability to distinguish between two different finger-tapping tasks was significantly higher [p = 0.00466] for the reconstructed STSS-SRfMRI images than for the raw fMRI images. The results indicate that the functional resolution of the STSS-SRfMRI scheme is superior, which suggests that the scheme is a potential solution to realizing higher functional resolution in fMRI images obtained using 3T MRI.

Automatic worm detection to solve overlapping problems using a convolutional neural network.

The nematode Caenorhabditis elegans is a powerful experimental model to investigate vital functions of higher organisms. We recently established a novel method, named "pond assay for the sensory systems (PASS)", that dramatically improves both the evaluation accuracy of sensory response of worms and the efficiency of experiments. This method uses many worms in numbers that are impractical to count manually. Although several automated detection systems have been introduced, detection of overlapped worms remains difficult. To overcome this problem, we developed an automated worm detection system based on a deep neural network (DNN). Our DNN was based on a "YOLOv4″ one-stage detector with one-class classification (OCC) and multi-class classification (MCC). The OCC defined a single class for worms, while the MCC defined four classes for the number of overlapped worms. For the training data, a total of 2000 model sub-images were prepared by manually drawing square worm bounding boxes from 150 images. To make simulated images, a total of 10-80 model images for each class were randomly selected and randomly placed on a simulated microscope field. A total of 19,000 training datasets and 1000 validation datasets with a ground-truth bounding-box were prepared. We evaluated detection accuracy using 150 images, which were different from the training data. Evaluation metrics were detection error, precision, recall, and average precision (AP). Precision values were 0.91 for both OCC and MCC. However, the recall value for MCC (= 0.93) was higher than that for OCC (= 0.79). The number of detection errors for OCC increased with increasing the ground truth; however, that for MCC was independent of the ground truth. AP values were 0.78 and 0.90 for the OCC and the MCC, respectively. Our worm detection system with MCC provided better detection accuracy for large numbers of worms with overlapping positions than that with the OCC.

Quantitative measurement of diffusion-weighted imaging signal using expression-controlled aquaporin-4 cells: Comparative study of 2-compartment and diffusion kurtosis imaging models.

The purpose of this study was to compare parameter estimates for the 2-compartment and diffusion kurtosis imaging models obtained from diffusion-weighted imaging (DWI) of aquaporin-4 (AQP4) expression-controlled cells, and to look for biomarkers that indicate differences in the cell membrane water permeability. DWI was performed on AQP4-expressing and non-expressing cells and the signal was analyzed with the 2-compartment and diffusion kurtosis imaging models. For the 2-compartment model, the diffusion coefficients (Df, Ds) and volume fractions (Ff, Fs, Ff = 1-Fs) of the fast and slow compartments were estimated. For the diffusion kurtosis imaging model, estimates of the diffusion kurtosis (K) and corrected diffusion coefficient (D) were obtained. For the 2-compartment model, Ds and Fs showed clear differences between AQP4-expressing and non-expressing cells. Fs was also sensitive to cell density. There was no clear relationship with the cell type for the diffusion kurtosis imaging model parameters. Changes to cell membrane water permeability due to AQP4 expression affected DWI of cell suspensions. For the 2-compartment and diffusion kurtosis imaging models, Ds was the parameter most sensitive to differences in AQP4 expression.

Distribution of intraperitoneally administered deuterium-labeled water in aquaporin-4-knockout mouse brain after middle cerebral artery occlusion.

Introduction: As the movement of water in the brain is known to be involved in neural activity and various brain pathologies, the ability to assess water dynamics in the brain will be important for the understanding of brain function and the diagnosis and treatment of brain diseases. Aquaporin-4 (AQP4) is a membrane channel protein that is highly expressed in brain astrocytes and is important for the movement of water molecules in the brain. Methods: In this study, we investigated the contribution of AQP4 to brain water dynamics by administering deuterium-labeled water (D2O) intraperitoneally to wild-type and AQP4 knockout (AQP4-ko) mice that had undergone surgical occlusion of the middle cerebral artery (MCA). Water dynamics in the infarct region and on either side of the anterior cerebral artery (ACA) was monitored with proton-density-weighted imaging (PDWI) performed on a 7T animal MRI. Results: D2O caused a negative signal change quickly after administration. The AQP4-ko mice showed a delay of the time-to-minimum in both the contralateral and ipsilateral ACA regions compared to wild-type mice. Also, only the AQP4- ko mice showed a delay of the time-to-minimum in the ipsilateral ACA region compared to the contralateral side. In only the wild-type mice, the signal minimum in the ipsilateral ACA region was higher than that in the contralateral ACA region. In the infarct region, the signal attenuation was slower for the AQP4-ko mice in comparison to the wild-type mice. Discussion: These results suggest that AQP4 loss affects water dynamics in the ACA region not only in the infarct region. Dynamic PDWI after D2O administration may be a useful tool for showing the effects of AQP4 in vivo.

Exploring cell membrane water exchange in aquaporin-4-deficient ischemic mouse brain using diffusion-weighted MRI.

BACKGROUND: Aquaporin-4 is a membrane channel protein that is highly expressed in brain astrocytes and facilitates the transport of water molecules. It has been suggested that suppression of aquaporin-4 function may be an effective treatment for reducing cellular edema after cerebral infarction. It is therefore important to develop clinically applicable measurement systems to evaluate and better understand the effects of aquaporin-4 suppression on the living body. METHODS: Animal models of focal cerebral ischemia were created by surgically occluding the middle cerebral artery of wild-type and aquaporin-4 knockout mice, after which multi-b-value multi-diffusion-time diffusion-weighted imaging measurements were performed. Data were analyzed with both the apparent diffusion coefficient (ADC) model and a compartmental water-exchange model. RESULTS: ADCs were estimated for five different b value ranges. The ADC of aquaporin-4 knockout mice in the contralateral region was significantly higher than that of wild-type mice for each range. In contrast, aquaporin-4 knockout mice had significantly lower ADC than wild-type mice in ischemic tissue for each b-value range. Genotype-dependent differences in the ADC were particularly significant for the lowest ranges in normal tissue and for the highest ranges in ischemic tissue. The ADCs measured at different diffusion times were significantly different for both genotypes. Fitting of the water-exchange model to the ischemic region data found that the water-exchange time in aquaporin-4 knockout mice was approximately 2.5 times longer than that in wild-type mice. CONCLUSIONS: Multi-b-value multi-diffusion-time diffusion-weighted imaging may be useful for in vivo research and clinical diagnosis of aquaporin-4-related diseases.

Monte Carlo Simulator for Diffusion-weighted Imaging Sequences.

We developed a Monte Carlo simulator for diffusion-weighted imaging sequences which displays the motion of water molecules and computes the dynamic phase dispersion due to the applied motion probing gradients. This simulator can be used to validate the analytical equations of diffusion models and understand their limitations due to their approximations. Here, we introduce the software and some specific use cases. The software can be downloaded from the following website: https://www.nirs.qst.go.jp/amr_diag.

New algorithm using L1 regularization for measuring electron energy spectra.

Retrieving the spectrum of physical radiation from experimental measurements typically involves using a mathematical algorithm to deconvolve the instrument response function from the measured signal. However, in the field of signal processing known as "Source Separation" (SS), which refers to the process of computationally retrieving the separate source components that generate an overlapping signal on the detector, the deconvolution process can become an ill-posed problem and crosstalk complicates the separation of the individual sources. To overcome this problem, we have designed a magnetic spectrometer for inline electron energy spectrum diagnosis and developed an analysis algorithm using techniques applicable to the problem of SS. An unknown polychromatic electron spectrum is calculated by sparse coding using a Gaussian basis function and an L1 regularization algorithm with a sparsity constraint. This technique is verified by using a specially designed magnetic field electron spectrometer. We use Monte Carlo simulations of the detector response to Maxwellian input energy distributions with electron temperatures of 5.0 MeV, 10.0 MeV, and 15.0 MeV to show that the calculated sparse spectrum can reproduce the input spectrum with an optimum energy bin width automatically selected by the L1 regularization. The spectra are reproduced with a high accuracy of less than 4.0% error, without an initial value. The technique is then applied to experimental measurements of intense laser accelerated electron beams from solid targets. Our analysis concept of spectral retrieval and automatic optimization of energy bin width by sparse coding could form the basis of a novel diagnostic method for spectroscopy.

Multiple sclerosis lesions in motor tracts from brain to cervical cord: spatial distribution and correlation with disability.

Despite important efforts to solve the clinico-radiological paradox, correlation between lesion load and physical disability in patients with multiple sclerosis remains modest. One hypothesis could be that lesion location in corticospinal tracts plays a key role in explaining motor impairment. In this study, we describe the distribution of lesions along the corticospinal tracts from the cortex to the cervical spinal cord in patients with various disease phenotypes and disability status. We also assess the link between lesion load and location within corticospinal tracts, and disability at baseline and 2-year follow-up. We retrospectively included 290 patients (22 clinically isolated syndrome, 198 relapsing remitting, 39 secondary progressive, 31 primary progressive multiple sclerosis) from eight sites. Lesions were segmented on both brain (T2-FLAIR or T2-weighted) and cervical (axial T2- or T2*-weighted) MRI scans. Data were processed using an automated and publicly available pipeline. Brain, brainstem and spinal cord portions of the corticospinal tracts were identified using probabilistic atlases to measure the lesion volume fraction. Lesion frequency maps were produced for each phenotype and disability scores assessed with Expanded Disability Status Scale score and pyramidal functional system score. Results show that lesions were not homogeneously distributed along the corticospinal tracts, with the highest lesion frequency in the corona radiata and between C2 and C4 vertebral levels. The lesion volume fraction in the corticospinal tracts was higher in secondary and primary progressive patients (mean = 3.6 ± 2.7% and 2.9 ± 2.4%), compared to relapsing-remitting patients (1.6 ± 2.1%, both P < 0.0001). Voxel-wise analyses confirmed that lesion frequency was higher in progressive compared to relapsing-remitting patients, with significant bilateral clusters in the spinal cord corticospinal tracts (P < 0.01). The baseline Expanded Disability Status Scale score was associated with lesion volume fraction within the brain (r = 0.31, P < 0.0001), brainstem (r = 0.45, P < 0.0001) and spinal cord (r = 0.57, P < 0.0001) corticospinal tracts. The spinal cord corticospinal tracts lesion volume fraction remained the strongest factor in the multiple linear regression model, independently from cord atrophy. Baseline spinal cord corticospinal tracts lesion volume fraction was also associated with disability progression at 2-year follow-up (P = 0.003). Our results suggest a cumulative effect of lesions within the corticospinal tracts along the brain, brainstem and spinal cord portions to explain physical disability in multiple sclerosis patients, with a predominant impact of intramedullary lesions.

The Utility of a Convolutional Neural Network for Generating a Myelin Volume Index Map from Rapid Simultaneous Relaxometry Imaging.

PURPOSE: A current algorithm to obtain a synthetic myelin volume fraction map (SyMVF) from rapid simultaneous relaxometry imaging (RSRI) has a potential problem, that it does not incorporate information from surrounding pixels. The purpose of this study was to develop a method that utilizes a convolutional neural network (CNN) to overcome this problem. METHODS: RSRI and magnetization transfer images from 20 healthy volunteers were included. A CNN was trained to reconstruct RSRI-related metric maps into a myelin volume-related index (generated myelin volume index: GenMVI) map using the MVI map calculated from magnetization transfer images (MTMVI) as reference. The SyMVF and GenMVI maps were statistically compared by testing how well they correlated with the MTMVI map. The correlations were evaluated based on: (i) averaged values obtained from 164 atlas-based ROIs, and (ii) pixel-based comparison for ROIs defined in four different tissue types (cortical and subcortical gray matter, white matter, and whole brain). RESULTS: For atlas-based ROIs, the overall correlation with the MTMVI map was higher for the GenMVI map than for the SyMVF map. In the pixel-based comparison, correlation with the MTMVI map was stronger for the GenMVI map than for the SyMVF map, and the difference in the distribution for the volunteers was significant (Wilcoxon sign-rank test, P < 0.001) in all tissue types. CONCLUSION: The proposed method is useful, as it can incorporate more specific information about local tissue properties than the existing method. However, clinical validation is necessary.

Imaging of Hypoxic Tumor: Correlation between Diffusion-weighted MR Imaging and 18F-fluoroazomycin Arabinoside Positron Emission Tomography in Head and Neck Carcinoma.

We investigated the usefulness of diffusion-weighted imaging (DWI) for detecting changes in the structure of hypoxic cells by evaluating the correlation between 18F-fluoroazomycin arabinoside (FAZA) positron emission tomography activity and DWI parameters in head and neck carcinoma. The diffusion coefficient corresponding to the slow compartment of a two-compartment model had a significant positive correlation with FAZA activity (ρ = 0.58, P = 0.016), whereas the diffusional kurtosis from diffusion kurtosis imaging had a significant negative correlation (ρ = -0.62, P = 0.008), which suggests that those DWI parameters might be useful as indicators for changes in cell structure.

The Utility of Applying Various Image Preprocessing Strategies to Reduce the Ambiguity in Deep Learning-based Clinical Image Diagnosis.

PURPOSE: A general problem of machine-learning algorithms based on the convolutional neural network (CNN) technique is that the reason for the output judgement is unclear. The purpose of this study was to introduce a strategy that may facilitate better understanding of how and why a specific judgement was made by the algorithm. The strategy is to preprocess the input image data in different ways to highlight the most important aspects of the images for reaching the output judgement. MATERIALS AND METHODS: T2-weighted brain image series falling into two age-ranges were used. Classifying each series into one of the two age-ranges was the given task for the CNN model. The images from each series were preprocessed in five different ways to generate five different image sets: (1) subimages from the inner area of the brain, (2) subimages from the periphery of the brain, (3-5) subimages of brain parenchyma, gray matter area, and white matter area, respectively, extracted from the subimages of (2). The CNN model was trained and tested in five different ways using one of these image sets. The network architecture and all the parameters for training and testing remained unchanged. RESULTS: The judgement accuracy achieved by training was different when the image set used for training was different. Some of the differences was statistically significant. The judgement accuracy decreased significantly when either extra-parenchymal or gray matter area was removed from the periphery of the brain (P < 0.05). CONCLUSION: The proposed strategy may help visualize what features of the images were important for the algorithm to reach correct judgement, helping humans to understand how and why a particular judgement was made by a CNN.

Reference region extraction by clustering for the pharmacokinetic analysis of dynamic contrast-enhanced MRI in prostate cancer.

PURPOSE: Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) measures changes in the concentration of an administered contrast agent to quantitatively evaluate blood circulation in a tumor or normal tissues. This method uses a pharmacokinetic analysis based on the time course of a reference region, such as muscle, rather than arterial input function. However, it is difficult to manually define a homogeneous reference region. In the present study, we developed a method for automatic extraction of the reference region using a clustering algorithm based on a time course pattern for DCE-MRI studies of patients with prostate cancer. METHODS: Two feature values related to the shape of the time course were extracted from the time course of all voxels in the DCE-MRI images. Each voxel value of T1-weighted images acquired before administration were also added as anatomical data. Using this three-dimensional feature vector, all voxels were segmented into five clusters by the Gaussian mixture model, and one of these clusters that included the gluteus muscle was selected as the reference region. RESULTS: Each region of arterial vessel, muscle, and fat was segmented as a different cluster from the tumor and normal tissues in the prostate. In the extracted reference region, other tissue elements including scattered fat and blood vessels were removed from the muscle region. CONCLUSIONS: Our proposed method can automatically extract the reference region using the clustering algorithm with three types of features based on the time course pattern and anatomical data. This method may be useful for evaluating tumor circulatory function in DCE-MRI studies.

Comparison of magnetization transfer contrast of conventional and simultaneous multislice turbo spin echo acquisitions focusing on excitation time interval.

PURPOSE: Image contrast differs between conventional multislice turbo spin echo (conventional TSE) and multiband turbo spin echo (SMS-TSE). Difference in time interval between excitations for adjacent slices (SETI) might cause this difference. This study aimed to evaluate the influence of SETI on MT effect for conventional TSE and compare conventional TSE with SMS-TSE in this respect. MATERIALS AND METHODS: Three different agar concentration phantoms were scanned with conventional TSE by adjusting SETI and TR. Signal change for different SETI was evaluated using Pearson's correlation analysis. SMS-TSE was acquired by changing TR similarly. Three human volunteers were scanned with similar settings to evaluate reproducibility of the phantom results in human brain. RESULTS: In conventional TSE, shorter SETI induced larger signal reduction. Longer TR and higher agar concentration emphasized this characteristic. Significant linear correlation (P < 0.05) was found in the major cases. The SMS-TSE signal intensity in each TR and phantom was smaller than the assumable levels in conventional TSE when the slices were simultaneously excited. Similar characteristic was observed in human brain. CONCLUSION: Shorter SETI results in larger MT effect in conventional TSE. The contrast change in SMS-TSE was larger than the supposable level from simultaneous excitation, which needs consideration in clinics.

Intra- and inter-operator reproducibility of US point shear-wave elastography in various organs: evaluation in phantoms and healthy volunteers.

PURPOSE: This study was conducted in order to assess the intra- and interoperator reproducibility of shear-wave speed (SWS) measurement on elasticity phantoms and healthy volunteers using ultrasound-based point shear-wave elastography. MATERIALS AND METHODS: This study was approved by the institutional review board. Two operators measured the SWS of five elasticity phantoms and seven organs (thyroid, lymph node, muscle, spleen, kidney, pancreas, and liver) of 30 healthy volunteers with 1.0-4.5 MHz convex (4C1) and 4.0-9.0 MHz linear (9L4) transducers. The phantom measurements were repeated ten times, while the volunteer measurements were performed five times each. Intra- and interoperator reproducibility was assessed. Interoperator reproducibility was also evaluated with the 95% Bland-Altman limits of agreement (LOA). RESULTS: In phantoms, all intraclass correlation coefficients (ICCs) were above 0.90 and the 95% LOA between the two operators were less than ± 18%. In volunteers, intraoperator ICCs were > 0.75 for all regions except the pancreas. Interoperator ICC was above 0.75 for the right lobe of the liver (depth 4 cm) and the kidney, but the 95% LOA was less than ± 25% only for the liver. CONCLUSION: Although excellent in phantoms, interoperator reproducibility was insufficient for all regions in the volunteers other than the right hepatic lobe at a depth of 4 cm. Clinicians should be aware of the 95% LOA when using SWS in patients. KEY POINTS: • Our phantom study indicated a high reproducibility for shear-wave speed (SWS) measurements with point shear-wave elastography (pSWE). • In volunteers, intraoperator reproducibility was generally high, but the interoperator reproducibility was not high enough except for the right hepatic lobe at 4 cm depth. • To evaluate interoperator reproducibility, the 95% limits of agreement (LOA) between operators should be considered in addition to the intraclass correlation coefficient (ICC).

Spatial distribution of multiple sclerosis lesions in the cervical spinal cord.

Spinal cord lesions detected on MRI hold important diagnostic and prognostic value for multiple sclerosis. Previous attempts to correlate lesion burden with clinical status have had limited success, however, suggesting that lesion location may be a contributor. Our aim was to explore the spatial distribution of multiple sclerosis lesions in the cervical spinal cord, with respect to clinical status. We included 642 suspected or confirmed multiple sclerosis patients (31 clinically isolated syndrome, and 416 relapsing-remitting, 84 secondary progressive, and 73 primary progressive multiple sclerosis) from 13 clinical sites. Cervical spine lesions were manually delineated on T2- and T2*-weighted axial and sagittal MRI scans acquired at 3 or 7 T. With an automatic publicly-available analysis pipeline we produced voxelwise lesion frequency maps to identify predilection sites in various patient groups characterized by clinical subtype, Expanded Disability Status Scale score and disease duration. We also measured absolute and normalized lesion volumes in several regions of interest using an atlas-based approach, and evaluated differences within and between groups. The lateral funiculi were more frequently affected by lesions in progressive subtypes than in relapsing in voxelwise analysis (P < 0.001), which was further confirmed by absolute and normalized lesion volumes (P < 0.01). The central cord area was more often affected by lesions in primary progressive than relapse-remitting patients (P < 0.001). Between white and grey matter, the absolute lesion volume in the white matter was greater than in the grey matter in all phenotypes (P < 0.001); however when normalizing by each region, normalized lesion volumes were comparable between white and grey matter in primary progressive patients. Lesions appearing in the lateral funiculi and central cord area were significantly correlated with Expanded Disability Status Scale score (P < 0.001). High lesion frequencies were observed in patients with a more aggressive disease course, rather than long disease duration. Lesions located in the lateral funiculi and central cord area of the cervical spine may influence clinical status in multiple sclerosis. This work shows the added value of cervical spine lesions, and provides an avenue for evaluating the distribution of spinal cord lesions in various patient groups.