Validation of Dynamic 3D MRI for Urodynamics Assessment Using an Anatomically Realistic In Vitro Model of the Bladder

Lowery urinary tract symptoms (LUTS) affect a large majority of the aging population. 3D Dynamic MRI shows promise as a noninvasive diagnostic tool that can assess bladder anatomy and function (urodynamics) while overcoming challenges associated with current urodynamic assessment methods. However, validation of this technique remains an unmet need. In this study, an anatomically realistic, bladder-mimicking in vitro flow model was created and used to systematically benchmark 3D dynamic MRI performance using a highly controllable syringe pump. Time-resolved volumes of the synthetic bladder model were obtained during simulated filling and voiding events and used to calculate volumetric flowrate. During MRI acquisitions, pressure during each event was recorded and used to create PV loops for work assessment. Error between control and MRI-derived volume for voiding and filling events exhibited 3.36% and 4.66% differences, respectively. A slight increase in average error was observed for MRI-derived flowrate when compared to the control flowrate (4.90% and 7.67% for voiding and filling, respectively). Overall, average error in segmented volumes increased with decreasing volume flowrate. Pressure drops were observed during voiding. Pressure increased during filling. Enhanced validation of novel 3D MRI urodynamics is achieved by using high-resolution PIV for visualizing and quantifying velocity inside the bladder model, which is not currently possible with 3D Dynamic MRI.

Movienet: Deep space-time-coil reconstruction network without k-space data consistency for fast motion-resolved 4D MRI

PURPOSE

To develop a novel deep learning approach for 4D-MRI reconstruction, named Movienet, which exploits space-time-coil correlations and motion preservation instead of k-space data consistency, to accelerate the acquisition of golden-angle radial data and enable subsecond reconstruction times in dynamic MRI.

METHODS

Movienet uses a U-net architecture with modified residual learning blocks that operate entirely in the image domain to remove aliasing artifacts and reconstruct an unaliased motion-resolved 4D image. Motion preservation is enforced by sorting the input image and reference for training in a linear motion order from expiration to inspiration. The input image was collected with a lower scan time than the reference XD-GRASP image used for training. Movienet is demonstrated for motion-resolved 4D MRI and motion-resistant 3D MRI of abdominal tumors on a therapeutic 1.5T MR-Linac (1.5-fold acquisition acceleration) and diagnostic 3T MRI scanners (2-fold and 2.25-fold acquisition acceleration for 4D and 3D, respectively). Image quality was evaluated quantitatively and qualitatively by expert clinical readers.

RESULTS

The reconstruction time of Movienet was 0.69 s (4 motion states) and 0.75 s (10 motion states), which is substantially lower than iterative XD-GRASP and unrolled reconstruction networks. Movienet enables faster acquisition than XD-GRASP with similar overall image quality and improved suppression of streaking artifacts.

CONCLUSION

Movienet accelerates data acquisition with respect to compressed sensing and reconstructs 4D images in less than 1 s, which would enable an efficient implementation of 4D MRI in a clinical setting for fast motion-resistant 3D anatomical imaging or motion-resolved 4D imaging.

Optimization of 3D dynamic speech MRI: Poisson-disc undersampling and locally higher-rank reconstruction through partial separability model with regional optimized temporal basis

PURPOSE

To improve the spatiotemporal qualities of images and dynamics of speech MRI through an improved data sampling and image reconstruction approach.

METHODS

For data acquisition, we used a Poisson-disc random under sampling scheme that reduced the undersampling coherence. For image reconstruction, we proposed a novel locally higher-rank partial separability model. This reconstruction model represented the oral and static regions using separate low-rank subspaces, therefore, preserving their distinct temporal signal characteristics. Regional optimized temporal basis was determined from the regional-optimized virtual coil approach. Overall, we achieved a better spatiotemporal image reconstruction quality with the potential of reducing total acquisition time by 50%.

RESULTS

The proposed method was demonstrated through several 2-mm isotropic, 64 mm total thickness, dynamic acquisitions with 40 frames per second and compared to the previous approach using a global subspace model along with other k-space sampling patterns. Individual timeframe images and temporal profiles of speech samples were shown to illustrate the ability of the Poisson-disc under sampling pattern in reducing total acquisition time. Temporal information of sagittal and coronal directions was also shown to illustrate the effectiveness of the locally higher-rank operator and regional optimized temporal basis. To compare the reconstruction qualities of different regions, voxel-wise temporal SNR analysis were performed.

CONCLUSION

Poisson-disc sampling combined with a locally higher-rank model and a regional-optimized temporal basis can drastically improve the spatiotemporal image quality and provide a 50% reduction in overall acquisition time.

Universal generative modeling in dual domains for dynamic MRI

Dynamic magnetic resonance image reconstruction from incomplete k-space data has generated great research interest due to its ability to reduce scan time. Nevertheless, the reconstruction problem remains a thorny issue due to its ill posed nature. Recently, diffusion models, especially score-based generative models, have demonstrated great potential in terms of algorithmic robustness and flexibility of utilization. Moreover, a unified framework through the variance exploding stochastic differential equation is proposed to enable new sampling methods and further extend the capabilities of score-based generative models. Therefore, by taking advantage of the unified framework, we propose a k-space and image dual-domain collaborative universal generative model (DD-UGM), which combines the score-based prior with a low-rank regularization penalty to reconstruct highly under-sampled measurements. More precisely, we extract prior components from both image and k-space domains via a universal generative model and adaptively handle these prior components for faster processing while maintaining good generation quality. Experimental comparisons demonstrate the noise reduction and detail preservation abilities of the proposed method. Moreover, DD-UGM can reconstruct data of different frames by only training a single frame image, which reflects the flexibility of the proposed model.

Super-Resolved Dynamic 3D Reconstruction of the Vocal Tract during Natural Speech

MRI is the gold standard modality for speech imaging. However, it remains relatively slow, which complicates imaging of fast movements. Thus, an MRI of the vocal tract is often performed in 2D. While 3D MRI provides more information, the quality of such images is often insufficient. The goal of this study was to test the applicability of super-resolution algorithms for dynamic vocal tract MRI. In total, 25 sagittal slices of 8 mm with an in-plane resolution of 1.6 × 1.6 mm2 were acquired consecutively using a highly-undersampled radial 2D FLASH sequence. The volunteers were reading a text in French with two different protocols. The slices were aligned using the simultaneously recorded sound. The super-resolution strategy was used to reconstruct 1.6 × 1.6 × 1.6 mm3 isotropic volumes. The resulting images were less sharp than the native 2D images but demonstrated a higher signal-to-noise ratio. It was also shown that the super-resolution allows for eliminating inconsistencies leading to regular transitions between the slices. Additionally, it was demonstrated that using visual stimuli and shorter text fragments improves the inter-slice consistency and the super-resolved image sharpness. Therefore, with a correct speech task choice, the proposed method allows for the reconstruction of high-quality dynamic 3D volumes of the vocal tract during natural speech.

Cervical Foraminal Changes in Patients with Intermittent Arm Radiculopathy Studied with a New MRI-Compatible Compression Device

Diagnosing cervical foraminal stenosis with intermittent arm radiculopathy is challenging due to discrepancies between MRI findings and symptoms. This can be attributed to the fact that MRI images are often obtained in a relaxed supine position. This study aims to evaluate the feasibility of the Dynamic MRI Compression System (DMRICS) and to assess possible changes in cervical foramina, with both quantitative measurements and qualitative grading systems, with MRI during a simulated Spurling test. Ten patients (five women and five men, ages 29-45) with previously confirmed cervical foraminal stenosis underwent MRI scans using DMRICS. MRI images were acquired in both relaxed and provoked states. A radiologist assessed 30 foramina (C4-C7) on the symptomatic side in both patient positions. Quantitative and qualitative measures were performed, including the numeric rating scale (NRS) and the Park and Kim grading systems. The provoked state induced concordant neck and arm pain in 9 of 10 patients. Significant shifts in Park and Kim foraminal gradings were noted: 13 of 27 Park gradings and 9 of 27 Kim gradings escalated post provocation. No quantitative changes were observed. This pilot study indicates that the DMRICS device has the potential to improve diagnostic accuracy for cervical radiculopathy, demonstrating induced cervical foraminal changes during a simulated Spurling test while performing MRI.

Time-Resolved Quantification of Patellofemoral Cartilage Deformation in Response to Loading and Unloading via Dynamic MRI With Prospective Motion Correction

BACKGROUND

In vivo cartilage deformation has been studied by static magnetic resonance imaging (MRI) with in situ loading, but knowledge about strain dynamics after load onset and release is scarce.

PURPOSE

To measure the dynamics of patellofemoral cartilage deformation and recovery in response to in situ loading and unloading by using MRI with prospective motion correction.

STUDY TYPE

Prospective.

SUBJECTS

Ten healthy male volunteers (age: [31.4 ± 3.2] years).

FIELD STRENGTH/SEQUENCE

T1-weighted RF-spoiled 2D gradient-echo sequence with a golden angle radial acquisition scheme, augmented with prospective motion correction, at 3 T.

ASSESSMENT

In situ knee loading was realized with a flexion angle of approximately 40° using an MR-compatible pneumatic loading device. The loading paradigm consisted of 2 minutes of unloaded baseline followed by a 5-minute loading bout with 50% body weight and an unloading period of 38 minutes. The cartilage strain was assessed as the mean distance between patellar and femoral bone-cartilage interfaces as a percentage of the initial (pre-load) distance.

STATISTICAL TESTS

Wilcoxon signed-rank tests (significance level: P < 0.05), Pearson correlation coefficient (r).

RESULTS

The cartilage compression and recovery behavior was characterized by a viscoelastic response. The elastic compression ([-12.5 ± 3.1]%) was significantly larger than the viscous compression ([-7.6 ± 1.5]%) and the elastic recovery ([10.5 ± 2.1]%) was significantly larger than the viscous recovery ([6.1 ± 1.8]%). There was a significant residual offset strain ([-3.6 ± 2.3]%) across the cohort. A significant negative correlation between elastic compression and elastic recovery was observed (r = -0.75).

DATA CONCLUSION

The in vivo cartilage compression and recovery time course in response to loading was successfully measured via dynamic MRI with prospective motion correction. The clinical relevance of the strain characteristics needs to be assessed in larger subject and patient cohorts.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 1.

Clinical profile and dynamic magnetic resonance imaging in Hirayama disease: a single-centered cross-sectional study in Nepal

Hirayama disease (HD) is juvenile monomelic amyotrophy of the distal upper limb first described by Hirayama in 1959 AD. HD is a benign condition with chronic microcirculatory changes. The hallmark of HD is necrosis of the anterior horns of the distal cervical spine.

Materials and Methods

Eighteen patients were assessed for clinical and radiological Hirayama disease. Clinical criteria included insidious onset nonprogressive chronic upper limb weakness and atrophy in teens or early twenties without sensory deficits and coarse tremors. MRI was done in a neutral position followed by neck flexion to evaluate cord atrophy and flattening, abnormal cervical curvature, loss of attachment between the posterior dural sac and subjacent lamina, anterior shifting of the posterior wall of the cervical dural canal, posterior epidural flow voids, and an enhancing epidural component with its dorsal extension.

Results

The mean age was 20.33 years, and the majority, 17 (94.4%), were male. Neutral-position MRI revealed loss of cervical lordosis in 5 (27.8%) patients, cord flattening in all patients with asymmetry in 10 (55.5%), and cord atrophy was observed in 13 (72.2%) patients with localized cervical cord atrophy in only 2 (11.1%) and extension of atrophy to dorsal cord in 11 (61.1%) patients. Intramedullary cord signal change was seen in 7 (38.9%) patients. Loss of attachment of posterior dura and subjacent lamina and anterior displacement of dorsal dura was seen in all patients. A crescent-shaped epidural intense enhancement was noted along the posterior aspect of the distal cervical canal in all patients, with dorsal level extension in 16 (88.89%) patients. The mean thickness of this epidural space was 4.38±2.26 (mean±2SD), and the mean extension was 5.5±4.6 vertebral levels (mean±2SD).

Conclusion

The high degree of clinical suspicion can guide additional contrast studies in flexion as a set MRI protocol for early detection and avoiding false negative diagnoses of HD.

Change in Descriptive Kinematic Parameters of Patients with Patellofemoral Instability When Compared to Individuals with Healthy Knees-A 3D MRI In Vivo Analysis

BACKGROUND

Patellofemoral instability (PFI) leads to chronic knee pain, reduced performance and chondromalacia patellae with consecutive osteoarthritis. Therefore, determining the exact patellofemoral contact mechanism, as well as the factors leading to PFI, is of great importance. The present study compares in vivo patellofemoral kinematic parameters and the contact mechanism of volunteers with healthy knees and patients with low flexion patellofemoral instability (PFI). The study was performed with a high-resolution dynamic MRI.

MATERIAL/METHODS

In a prospective cohort study, the patellar shift, patella rotation and the patellofemoral cartilage contact areas (CCA) of 17 patients with low flexion PFI were analyzed and compared with 17 healthy volunteers, matched via the TEA distance and sex, in unloaded and loaded conditions. MRI scans were carried out for 0°, 15° and 30° knee flexion in a custom-designed knee loading device. To suppress motion artifacts, motion correction was performed using a moiré phase tracking system with a tracking marker attached to the patella. The patellofemoral kinematic parameters and the CCA was calculated on the basis of semi-automated cartilage and bone segmentation and registrations.

RESULTS

Patients with low flexion PFI showed a significant reduction in patellofemoral CCA for 0° (unloaded: p = 0.002, loaded: p = 0.004), 15° (unloaded: p = 0.014, loaded: p = 0.001) and 30° (unloaded: p = 0.008; loaded: p = 0.001) flexion compared to healthy subjects. Additionally, patients with PFI revealed a significantly increased patellar shift when compared to volunteers with healthy knees at 0° (unloaded: p = 0.033; loaded: p = 0.031), 15° (unloaded: p = 0.025; loaded: p = 0.014) and 30° flexion (unloaded: p = 0.030; loaded: p = 0.034) There were no significant differences for patella rotation between patients with PFI and the volunteers, except when, under load at 0° flexion, PFI patients showed increased patellar rotation (p = 0.005. The influence of quadriceps activation on the patellofemoral CCA is reduced in patients with low flexion PFI.

CONCLUSION

Patients with PFI showed different patellofemoral kinematics at low flexion angles in both unloaded and loaded conditions compared to volunteers with healthy knees. Increased patellar shifts and decreased patellofemoral CCAs were observed in low flexion angles. The influence of the quadriceps muscle is diminished in patients with low flexion PFI. Therefore, the goal of patellofemoral stabilizing therapy should be to restore a physiologic contact mechanism and improve patellofemoral congruity for low flexion angles.

Quantifying Velopharyngeal Motion Variation in Speech Sound Production Using an Audio-Informed Dynamic MRI Atlas

New developments in dynamic magnetic resonance imaging (MRI) facilitate high-quality data acquisition of human velopharyngeal deformations in real-time speech. With recently established speech motion atlases, group analysis is made possible via spatially and temporally aligned datasets in the atlas space from a desired population of interest. In practice, when analyzing motion characteristics from various subjects performing a designated speech task, it is observed that different subjects' velopharyngeal deformation patterns could vary during the pronunciation of the same utterance, regardless of the spatial and temporal alignment of their MRI. Since such variation can be subtle, identification and extraction of unique patterns out of these high-dimensional datasets is a challenging task. In this work, we present a method that computes and visualizes subtle deformation variation patterns as principal components of a subject group's dynamic motion fields in the atlas space. Coupled with the real-time speech audio recordings during image acquisition, the key time frames that contain maximum speech variations are identified by the principal components of temporally aligned audio waveforms, which in turn inform the temporal location of the maximum spatial deformation variation. Henceforth, the motion fields between the key frames and the reference frame for each subject are computed and warped into the common atlas space, enabling a direct extraction of motion variation patterns via quantitative analysis. The method was evaluated on a dataset of twelve healthy subjects. Subtle velopharyngeal motion differences were visualized quantitatively to reveal pronunciation-specific patterns among different subjects.

4D Golden-Angle Radial MRI at Subsecond Temporal Resolution

Intraframe motion blurring, as a major challenge in free-breathing dynamic MRI, can be reduced if high temporal resolution can be achieved. To address this challenge, this work proposes a highly accelerated 4D (3D + time) dynamic MRI framework with subsecond temporal resolution that does not require explicit motion compensation. The method combines standard stack-of-stars golden-angle radial sampling and tailored GRASP-Pro (Golden-angle RAdial Sparse Parallel imaging with imProved performance) reconstruction. Specifically, 4D dynamic MRI acquisition is performed continuously without motion gating or sorting. The k-space centers in stack-of-stars radial data are organized to guide estimation of a temporal basis, with which GRASP-Pro reconstruction is employed to enforce joint low-rank subspace and sparsity constraints. This new basis estimation strategy is the new feature proposed for subspace-based reconstruction in this work to achieve high temporal resolution (e.g., subsecond/3D volume). It does not require sequence modification to acquire additional navigation data, it is compatible with commercially available stack-of-stars sequences, and it does not need an intermediate reconstruction step. The proposed 4D dynamic MRI approach was tested in abdominal motion phantom, free-breathing abdominal MRI, and dynamic contrast-enhanced MRI (DCE-MRI). Our results have shown that GRASP-Pro reconstruction with the new basis estimation strategy enables highly-accelerated 4D dynamic imaging at subsecond temporal resolution (with five spokes or less for each dynamic frame per image slice) for both free-breathing non-DCE-MRI and DCE-MRI. In the abdominal phantom, better image quality with lower root mean square error and higher structural similarity index was achieved using GRASP-Pro compared with standard GRASP. With the ability to acquire each 3D image in less than 1 s, intraframe respiratory blurring can be intrinsically reduced for body applications with our approach, which eliminates the need for explicit motion detection and motion compensation.

The Influence of Surgical Realignment Procedures on Dynamic Patellar Tracking: A Dynamic Magnetic Resonance Imaging-Controlled Feasibility Study

Persisting patellar maltracking following surgical realignment often remains unseen. The aim of this study was to analyze the effects of realignment procedures on patellofemoral kinematics in patients with patellofemoral instability (PFI) and patellofemoral maltracking (PM) by using dynamic magnetic resonance imaging (MRI). Patients planned for surgical patellar realignment due to PFI and a clinically and radiologically apparent PM between December 2019 and May 2022 were included. Patients without PM, limited range of motion, joint effusion, or concomitant injuries were excluded. Dynamic mediolateral translation (dMPT) and patella tilt (dPT) were measured preoperatively and three months postoperatively. In 24 patients (7 men, 17 women; mean age 23.0 years), 10 tibial tubercle transfers, 5 soft tissue patella tendon transfers, 6 trochleoplasties, 3 lateral lengthenings, 1 varizating distal femoral osteotomy (DFO), and 1 torsional DFO were performed. At final follow-up, dMPT (from 10.95 ± 5.93 mm to 4.89 ± 0.40 mm, p < 0.001) and dPT (from 14.50° ± 10.33° to 8.44° ± 7.46°, p = 0.026) were significantly improved. All static radiological parameters were corrected to physiological values. Surgical patellar realignment contributed to the significant improvement of patellofemoral kinematics, with an approximation to normal values. The postoperative application of dynamic MRI allowed for a quantification of the performed correction, allowing for a postoperative control of success.

Magnetic resonance imaging assessed enteric motility and luminal content analysis in patients with severe bloating and visible distension

BACKGROUND

Gastrointestinal symptoms in functional gut disorders occur without any discernible structural gut abnormality. Preliminary observations on enteric MRI suggest possible abnormal content and motility of the terminal ileum (TI) in constipation-predominant IBS (IBS-C) with severe bloating, and in functional bloating and distension (FABD) patients. We investigated whether MRI can quantify differences in small bowel (SB) content and motility between patients and healthy controls (HCs).

METHODS

11 IBS-C (mean age 40 [21-52] years; 10 women) and 7 FABD (36 [21-56]; all women) patients with bloating and 20 HCs (28 [22-48]; 6 women) underwent enteric MRI, including dynamic motility and anatomical sequences. Three texture analysis (TA) parameters assessed the homogeneity of the luminal content, with ratios calculated between the TI and (1) the SB and (2) the ascending colon. Four TI motility metrics were derived. Ascending colon diameter (ACD) was measured. A comparison between HCs and patients was performed independently for: (1) three TA parameters, (2) four TI motility metrics, and (3) ACD.

KEY RESULTS

Compared with HCs, patients had TI:colon ratios higher for TA contrast (p < 0.001), decreased TI motility (lower mean motility [p = 0.04], spatial motility variation [p = 0.03], and area of motile TI [p = 0.03]), and increased ACD (p = 0.001).

CONCLUSIONS AND INFERENCES

IBS-C and FABD patients show reduced TI motility and differences in luminal content compared with HCs. This potentially indicates reflux of colonic contents or delayed clearance of the TI, which alongside increased ACD may contribute to symptoms of constipation and bloating.

A novel method for measuring bowel motility and velocity with dynamic magnetic resonance imaging in two and three dimensions

Increasingly, dynamic magnetic resonance imaging (MRI) has potential as a noninvasive and accessible tool for diagnosing and monitoring gastrointestinal motility in healthy and diseased bowel. However, current MRI methods of measuring bowel motility have limitations: requiring bowel preparation or long acquisition times; providing mainly surrogate measures of motion; and estimating bowel-wall movement in just two dimensions. In this proof-of-concept study we apply a method that provides a quantitative measure of motion within the bowel, in both two and three dimensions, using existing, vendor-implemented MRI pulse sequences with minimal bowel preparation. This method uses a minimised cost function to fit linear vectors in the spatial and temporal domains. It is sensitised to the spatial scale of the bowel and aims to address issues relating to the low signal-to-noise in high-temporal resolution dynamic MRI scans, previously compensated for by performing thick-slice (10-mm) two-dimensional (2D) coronal scans. We applied both 2D and three-dimensional (3D) scanning protocols in two healthy volunteers. For 2D scanning, analysis yielded bi-modal velocity peaks, with a mean antegrade motion of 5.5 mm/s and an additional peak at ~9 mm/s corresponding to longitudinal peristalsis, as supported by intraoperative data from the literature. Furthermore, 3D scans indicated a mean forward motion of 4.7 mm/s, and degrees of antegrade and retrograde motion were also established. These measures show promise for the noninvasive assessment of bowel motility, and have the potential to be tuned to particular regions of interest and behaviours within the bowel.

Dynamic Mediolateral Patellar Translation Is a Sex- and Size-Independent Parameter of Adult Proximal Patellar Tracking Using Dynamic 3 Tesla Magnetic Resonance Imaging

PURPOSE

To provide normal values for physiological patellofemoral tracking in a representative group of healthy individuals, as well as sex differences, using real-time 3T-magnetic resonance imaging (MRI) and to test for the reliability of the presented technique.

METHODS

One hundred knees of healthy individuals with no history of patellofemoral symptoms were scanned with dynamic MRI sequences, during repetitive cycles of flexion (40°) and full extension. Within a 30-seconds time-frame, three simultaneous, transverse slices were acquired. Dynamic mediolateral patellar translation (dMPT) and dynamic patellar tilt (dPT) were measured on two occasions by two independent examiners. Common radiological parameters were measured using static MRI, and correlations were calculated.

RESULTS

100 knees (53 right, 47 left; age: 26.7 ± 4.4 years; BMI: 22.5 ± 3.1) of 57 individuals (27 females, 30 males) were included. Mean height was 170.1 ± 7.7 cm in women and 181.8 ± 6.4 cm in men. Average patella diameter was 37.9 ± 2.7 (95% CI 37.1-38.7) mm in women and 42.4 ± 3.2 (95% CI 41.5-43.3) mm in men. In females, the patellar diameters and intercondylar distances were significantly smaller than in males (P < .001). Radiological parameters for patellar maltracking were within the normal range. During the range of motion, mean dMPT was 1.7 ± 2.4 (95% CI .9-2.5) mm in females and 1.8 ± 2.7 (95% CI 1.1-2.6) mm in males (P = .766). Mean dPT was 1.3 ± 2.9° (95% CI .4-2.1°) in females and -0.2 ± 3.8° (95% CI -1.2-.9°) in males (P = .036). Neither dMPT nor dPT was correlated with height, BMI, or patellar diameter. Intercondylar distance correlated weakly with dPT (r = -.241; P = .041). Intra- and interrater reliability were excellent for dMPT and dPT.

CONCLUSION

Dynamic mediolateral patellar translation is a size- and sex-independent parameter for proximal patellar tracking. In healthy individuals without patellofemoral abnormalities normal dMPT proximal to the trochlea groove was 1.7 ± 2.5 (1.2-2.2) mm, independent of size or sex. Normal dPT showed a dependency on sex and was 1.3 ± 2.9 (.4-2.1)° in women and -0.2 ± 3.8 (-1.2-0.9)° in men.

LEVEL OF EVIDENCE

Level II, diagnostic study.

Characterization of surface motion patterns in highly deformable soft tissue organs from dynamic MRI: An application to assess 4D bladder motion

BACKGROUND AND OBJECTIVES

Dynamic Magnetic Resonance Imaging (MRI) may capture temporal anatomical changes in soft tissue organs with high-contrast but the obtained sequences usually suffer from limited volume coverage which makes the high-resolution reconstruction of organ shape trajectories a major challenge in temporal studies. Because of the variability of abdominal organ shapes across time and subjects, the objective of the present study is to go towards 3D dense velocity measurements to fully cover the entire surface and to extract meaningful features characterizing the observed organ deformations and enabling clinical action or decision.

METHODS

We present a pipeline for characterization of bladder surface dynamics during deep respiratory movements. For a compact shape representation, the reconstructed temporal volumes were first used to establish subject-specific dynamical 4D mesh sequences using the large deformation diffeomorphic metric mapping (LDDMM) framework. Then, we performed a statistical characterization of organ dynamics from mechanical parameters such as mesh elongations and distortions. Since we refer to organs as non-flat surfaces, we have also used the mean curvature change as metric to quantify surface evolution. However, the numerical computation of curvature is strongly dependant on the surface parameterization (i.e. the mesh resolution). To cope with this dependency, we employed a non-parametric method for surface deformation analysis. Independent of parameterization and minimizing the length of the geodesic curves, it stretches smoothly the surface curves towards a sphere by minimizing a Dirichlet energy. An Eulerian PDE approach is used to derive a shape descriptor from the curve-shortening flow. Intercorrelations between individuals' motion patterns are computed using the Laplace-Beltrami Operator (LBO) eigenfunctions for spherical mapping.

RESULTS

Application to extracting characterization correlation curves for locally-controlled simulated shape trajectories demonstrates the stability of the proposed shape descriptor. Its usability was shown on MRI acquired for seven healthy participants for which the bladder was highly deformed by maximum of inspiration. As expected, the study showed that deformations occured essentially on the top lateral regions.

CONCLUSION

Promising results were obtained, showing the organ in its 3D complexity during deformation due to strain conditions. Smooth genus-0 manifold reconstruction from sparse dynamic MRI data is employed to perform a statistical shape analysis for the determination of bladder deformation.

QdMRI: A system for comprehensive analysis of thoracic dynamics via dynamic MRI

Quantitative thoracic dynamic magnetic resonance imaging (QdMRI), a recently developed technique, provides a potential solution for evaluating treatment effects in thoracic insufficiency syndrome (TIS). In this paper, we integrate all related algorithms and modules during our work from the past 10 years on TIS into one system, named QdMRI, to address the following questions: (1) How to effectively acquire dynamic images? For many TIS patients, subjects are unable to cooperate with breathing instructions during image acquisition. Image acquisition can only be implemented under free-breathing conditions, and it is not feasible to use a surrogate device for tracing breathing signals. (2) How to assess the thoracic structures from the acquired image, such as lungs, left and right, separately? (3) How to depict the dynamics of thoracic structures due to respiration motion? (4) How to use the structural and functional information for the quantitative evaluation of surgical TIS treatment and for the design of the surgery plan? The QdMRI system includes 4 major modules: dynamic MRI (dMRI) acquisition, 4D image construction, image segmentation (from 4D image), and visualization of segmentation results, dynamic measurements, and comparisons of measurements from TIS patients with those from normal children. Scanning/image acquisition time for one subject is ~20 minutes, 4D image construction time is ~5 minutes, image segmentation of lungs via deep learning is 70 seconds for all time points (with the average DICE 0.96 in healthy children), and measurement computation time is 2 seconds.

Regional diaphragm motion analysis via dynamic MRI

Breathing-related movement analysis is important in the study of many disease processes. The analysis of diaphragmatic motion via thoracic imaging in particular is important in a variety of disorders. Compared to computed tomography (CT) and fluoroscopy, dynamic magnetic resonance imaging (dMRI) has several advantages, such as better soft tissue contrast, no ionizing radiation, and greater flexibility in selecting scanning planes. In this paper, we propose a novel method for full diaphragmatic motion analysis via free-breathing dMRI. Firstly, after 4D dMRI image construction in a cohort of 51 normal children, we manually delineated the diaphragm on sagittal plane dMRI images at end-inspiration and end-expiration. Then, 25 points were selected uniformly and homologously on each hemi-diaphragm surface. Based on the inferior-superior displacements of these 25 points between end-expiration (EE) and end-inspiration (EI) time points, we obtained their velocities. We then summarized 13 parameters from these velocities for each hemi-diaphragm to provide a quantitative regional analysis of diaphragmatic motion. We observed that the regional velocities of the right hemi-diaphragm were almost always statistically significantly greater than those of the left hemi-diaphragm in homologous locations. There was a significant difference for sagittal curvatures but not for coronal curvatures between the two hemi-diaphragms. Using this methodology, future larger scale prospective studies may be considered to confirm our findings in the normal state and to quantitatively assess regional diaphragmatic dysfunction when various disease conditions are present.

Translation and rotation analysis based on stress MRI for the diagnosis of anterior cruciate ligament tears

BACKGROUND

Due to the increasing need for a detailed biomechanical analysis of anterior cruciate ligament (ACL) lesions, the aim of the study was to develop a method of direct measurement of the three-dimensional tibial translation and rotation based on stress MRI.

METHODS

For the purpose of the study, thirty patients with acute ACL rupture and 17 healthy control subjects were selected. Based on clinical examination, they were qualified for MRI examination using the Arthroholder Device prototype to perform anterior tibial translation. Each examination was performed at 30° of knee flexion, initially without tibia translation and then using the force applied to the calf of 80 N. The femur and tibia were separately registered using rigid local SimpleITK landmark refinement; translation and rotation parameters were then calculated using the 3D transformation algorithms. The significance level was set at 0.05.

RESULTS

Initially, the device and method for obtaining the parameters of the 3D translation and rotation were validated. The pooled Standard Deviation for translation parameters was 0.81 mm and for rotation parameters 0.87°. Compared to the control group, statistically significant differences were found in parameters such as Anterior Shift [(median ± interquartile range) 3.89 mm ±6.55 vs. 0.90 mm ±2.78, P=0.002238] and External Rotation (-0.55° ±3.88 vs. -2.87° ±2.40, P=0.005074). Statistically significant correlations were observed in combined groups between Anterior Shift and parameters such as External Rotation (P=0.001611), PCL Tibial Attachment Point (pPCL) Anterior Shift (<0.000001), Rolimeter Measurement (P=0.000016), and Side-to-Side Difference (SSD) (P=0.000383). A significant statistical correlation was also observed between External Rotation and parameters such as Rolimeter (P=0.02261) and SSD (P=0.03458).

CONCLUSIONS

The analysis of the anterior tibia translation using stress MRI and the proposed three-dimensional calculation method allows for a detailed analysis of the tibial translation and rotation parameters. The correlations showed the importance of external rotation during anterior tibial translation.

Automatic upper airway segmentation in static and dynamic MRI via anatomy-guided convolutional neural networks

PURPOSE

Upper airway segmentation on MR images is a prerequisite step for quantitatively studying the anatomical structure and function of the upper airway and surrounding tissues. However, the complex variability of intensity and shape of anatomical structures and different modes of image acquisition commonly used in this application makes automatic upper airway segmentation challenging. In this paper, we develop and test a comprehensive deep learning-based segmentation system for use on MR images to address this problem.

MATERIALS AND METHODS

In our study, both static and dynamic MRI data sets are utilized, including 58 axial static 3D MRI studies, 22 mid-retropalatal dynamic 2D MRI studies, 21 mid-retroglossal dynamic 2D MRI studies, 36 mid-sagittal dynamic 2D MRI studies, and 23 isotropic dynamic 3D MRI studies, involving a total of 160 subjects and over 20 000 MRI slices. Samples of static and 2D dynamic MRI data sets were randomly divided into training, validation, and test sets by an approximate ratio of 5:2:3. Considering that the variability of annotation data among 3D dynamic MRIs was greater than for other MRI data sets, we increased the ratio of training data for these data to improve the robustness of the model. We designed a unified framework consisting of the following procedures. For static MRI, a generalized region-of-interest (GROI) strategy is applied to localize the partitions of nasal cavity and other portions of upper airway in axial data sets as two separate subobjects. Subsequently, the two subobjects are segmented by two separate 2D U-Nets. The two segmentation results are combined as the whole upper airway structure. The GROI strategy is also applied to other MRI modes. To minimize false-positive and false-negative rates in the segmentation results, we employed a novel loss function based explicitly on these rates to train the segmentation networks. An inter-reader study is conducted to test the performance of our system in comparison to human variability in ground truth (GT) segmentation of these challenging structures.

RESULTS

The proposed approach yielded mean Dice coefficients of 0.84±0.03, 0.89±0.13, 0.84±0.07, and 0.86±0.05 for static 3D MRI, mid-retropalatal/mid-retroglossal 2D dynamic MRI, mid-sagittal 2D dynamic MRI, and isotropic dynamic 3D MRI, respectively. The quantitative results show excellent agreement with manual delineation results. The inter-reader study results demonstrate that the segmentation performance of our approach is statistically indistinguishable from manual segmentations considering the inter-reader variability in GT.

CONCLUSIONS

The proposed method can be utilized for routine upper airway segmentation from static and dynamic MR images with high accuracy and efficiency. The proposed approach has the potential to be employed in other dynamic MRI-related applications, such as lung or heart segmentation.

A case report of Kuttner tumor mimicking a malignant tumor, leading to overtreatment

Preoperative diagnosis with multimodal approaches might lead to overtreatment. Cautious understanding of cytology and dynamic contrast-enhanced magnetic resonance imaging is required when a Kuttner tumor is cited as differential diagnosis.

Effect of sleep on upper airway dynamics in obese adolescents with obstructive sleep apnea syndrome

STUDY OBJECTIVES

The biomechanical basis of obstructive sleep apnea syndrome (OSAS) may influence upper airway dynamics. In this study, we investigate dynamic changes during respiration in wakefulness and sleep in obese adolescents with and without OSAS.

METHODS

Respiratory-gated dynamic magnetic resonance imaging (MRI) at the retropalatal and retroglossal regions was performed with simultaneous measurement of SpO2 and nasal-oral mask airflow and pressure. Airway cross-sectional area (CSA) was determined using AMIRA. Percent change in CSA was calculated from five continuous tidal breaths in states of wakefulness and sleep. Mixed effects models were used to evaluate interactions between group (OSAS/control), site (retropalatal/retroglossal), and stage (wake/sleep).

RESULTS

We studied 24 children with OSAS (mean age 15.49 ± 2.00 years, mean apnea-hypopnea index [AHI] 16.53 ± 8.72 events/h) and 19 controls (mean age 14.86 ± 1.75 years, mean AHI 2.12 ± 1.69 events/h). Groups were similar in age, sex, height, weight, and BMI Z-score. Participants with OSAS had a 48.17% greater increase in percent change of airway CSA during sleep than controls (p < 0.0001), while there was no difference between groups during wakefulness (p = 0.6589). Additionally, participants with OSAS had a 48.80% increase in percent change of airway CSA during sleep as compared with wakefulness (p < 0.0001), whereas no such relationship was observed in controls (p = 0.5513).

CONCLUSIONS

This study demonstrates significant effects of sleep on upper airway dynamics in obese children with OSAS. Dynamic MRI with physiological data can potentially provide further insight into the biomechanical basis of OSAS and assist in more effective management.

Data-Driven Regularization Parameter Selection in Dynamic MRI

In dynamic MRI, sufficient temporal resolution can often only be obtained using imaging protocols which produce undersampled data for each image in the time series. This has led to the popularity of compressed sensing (CS) based reconstructions. One problem in CS approaches is determining the regularization parameters, which control the balance between data fidelity and regularization. We propose a data-driven approach for the total variation regularization parameter selection, where reconstructions yield expected sparsity levels in the regularization domains. The expected sparsity levels are obtained from the measurement data for temporal regularization and from a reference image for spatial regularization. Two formulations are proposed. Simultaneous search for a parameter pair yielding expected sparsity in both domains (S-surface), and a sequential parameter selection using the S-curve method (Sequential S-curve). The approaches are evaluated using simulated and experimental DCE-MRI. In the simulated test case, both methods produce a parameter pair and reconstruction that is close to the root mean square error (RMSE) optimal pair and reconstruction. In the experimental test case, the methods produce almost equal parameter selection, and the reconstructions are of high perceived quality. Both methods lead to a highly feasible selection of the regularization parameters in both test cases while the sequential method is computationally more efficient.

Lung parenchymal characterization via thoracic dynamic MRI in normal children and pediatric patients with TIS

Quantitative thoracic dynamic magnetic resonance imaging (QdMRI), a recently developed technique, provides a potential solution for evaluating treatment effects in thoracic insufficiency syndrome (TIS). In this paper, we demonstrate how lung parenchymal characteristics can be assessed via intensity properties in lung dynamic MRI, a modality suitable for use in pediatric patients. The QdMRI-based approach includes dynamic MR image acquisition, 4D image construction, image pre-processing with non-uniformity correction and intensity standardization, and lung segmentation from the 4D constructed image via a deep learning approach, as well as extraction of image parenchymal intensity properties from the segmented lungs and statistical comparisons among different clinical scenarios. We include 22 dMRI scans from 11 TIS patients (each with both pre-operative and post-operative scans) and 23 dMRI scans from healthy children. Two-sided paired t-testing is performed to compare lung intensity properties between end of expiration (EE) and end of inspiration (EI) within TIS patients (pre-operative and post-operative, separately) and normal children. We also compare the lung intensity properties at EE and EI among pre-operative TIS patients, post-operative TIS patients, and normal children. Experimental results show that lung (T2) intensity at EI is significantly lower than that at EE and lung intensity of post-operative TIS patients is significantly lower than that in pre-operative TIS patients and closer to that of normal children than to that of pre-operative TIS patients, indicating improvement in lung aeration. To our knowledge, this is the first study to provide a quantitative dynamic functional method to analyze lung parenchyma during tidal breathing on dynamic MRI in both healthy children and pediatric patients with TIS.

Estimation of the dynamic volume of each lung via rapid limited-slice dynamic MRI

Dynamic lung volumetric parameters are useful for clinical assessment of many thoracic disorders, given that respiration is a dynamic process. Estimation of such parameters based on imaging and analysis is an important goal to achieve if implementation in routine clinical practice is to become a reality. Compared to CT, dynamic thoracic MRI has several advantages including better soft tissue contrast, lack of ionizing radiation, and flexibility in selecting scanning planes. 4D dynamic MRI seems to be the best choice for some clinical applications, notwithstanding the major limitation of a long image acquisition time (~45 minutes). Therefore, approaches to acquire images and estimate volumetric parameters rapidly is highly desirable in dynamic MRI-based clinical applications. In this paper, we present a technique for estimating lung volumetric parameters from limited-slices dynamic thoracic MRI, greatly reducing the number of slices to be scanned and therefore also the time required for image acquisition. We demonstrate a relative RMS error of predicted lung volumes of less than 5% by utilizing only 5 sagittal MRI slices through each lung compared to the current full scan involving about 20 slices per lung. As such, this approach can lead to time-saving during scan acquisition and therefore increased patient comfort and convenience for practical real-world clinical applications. This may potentially also improve image quality and usability due to the reduction of patient motion, abnormal breathing patterns, etc. ensuing from improved patient comfort and scan duration.

Cine magnetic resonance imaging provides novel predictors of early continence recovery after radical prostatectomy: Assessment of the dynamics of pelvic floor muscles

AIMS

Postprostatectomy incontinence is a major complication of prostatectomy. Although pelvic floor muscle training can successfully treat postprostatectomy incontinence, evidence for how muscle movement affects continence recovery is lacking. We evaluated dynamic factors of prostatectomy patients using cine magnetic resonance imaging to identify risk factors for postprostatectomy incontinence and reveal the contribution of pelvic floor muscles to continence recovery.

METHODS

A total of 128 prostate cancer patients who underwent robot-assisted laparoscopic surgery were enrolled. Cine magnetic resonance imaging was performed preoperatively and 6 months after surgery. Continence was defined as pad-free or use of safety pads. We defined the bladder neck elevation distance during pelvic floor muscle training as the bladder elevation distance. Patients with continence recovery within 1 month comprised the continence group (n = 48); other patients comprised the incontinence group (n = 80).

RESULTS

The preoperative bladder elevation distance was significantly longer in the continence group than in the incontinence group (10.4 vs. 8.2 mm; p < .001). The postoperative bladder elevation distance of the continence group tended to be longer (9.9 vs. 8.9 mm; p = .057). Multivariate analysis showed that the preoperative bladder elevation distance significantly contributed to continence recovery (p = .016). Patients with a longer preoperative bladder elevation distance (>8.5 mm) experienced continence recovery significantly faster than patients with a shorter distance (<8.5mm) (p = .038).

CONCLUSIONS

Bladder elevation distance, a novel dynamic parameter, was strongly associated with early continence recovery. Cine magnetic resonance imaging can assess a patient's risk of postprostatectomy incontinence and guide pelvic floor muscle training.

4D image construction from free-breathing MRI slice acquisitions of the thorax based on a concept of flux

Retrospective 4D image construction from continuously acquired 2D slices is a necessary step to achieve high-quality 4D images. Self-gating methods, which extract breathing signals only from image information without any external gating technology, have much potential, such as in pediatric patients with thoracic insufficiency syndrome (TIS) who suffer from extreme malformations of the chest wall, diaphragm, and spine, leading to breathing that is very complex with lots of abnormal respiration cycles, including very deep or shallow cycles. Existing methods do not work well in this clinical scenario and most are not fully automatic, requiring some manual interactive operations. In this paper, we propose a fully automatic 4D dMRI construction method based on the concept of flux to address the 4D image construction from 2D slices of subjects with complex respiration. Firstly, we extract the breathing signal for each location based on the flux of the optical flow vector field of the body region from the image series. Then, we give a full analysis for all cycles and extract several normal ones and map them to one cosine respiration model for each location. After that, we re-sample one normal cycle from the respiration model for each location independently. All of these resampled normal cycles form the final constructed 4D image. Qualitative and quantitative evaluations on 25 subjects show that the proposed method can handle datasets from subjects with more complex respiration and achieves good self-consistency results while maintaining time and space continuity.

Characteristics of dynamic magnetic resonance imaging of symptomatic chronic calcifying tendinitis: preliminary case reports

Background

The symptoms of chronic calcifying tendinitis consist of shoulder contracture and impingement sign. However, there have been no reports about the use of imaging studies to differentiate these 2 clinical symptoms. A "burning sign" caused by abnormal blood flow was previously reported in the shoulder joint in patients with frozen shoulder by dynamic magnetic resonance imaging. This burning sign was related to pain. The purpose of this study was to investigate the dynamic magnetic resonance imaging findings in patients with symptomatic chronic calcifying tendinitis and to examine the relationship between the location of the burning sign and the physical findings.

Methods

We retrospectively analyzed data for 6 patients with symptomatic chronic calcifying tendinitis (mean age, 55.5 ± 9.3 years; 4 women). The range of shoulder motion, impingement sign, and location of the burning sign were assessed.

Results

Four patients had an impingement sign without shoulder contracture, and the other 2 patients had shoulder contracture. All the patients with an impingement sign also had a burning sign around the calcium deposit and no enhancement in the rotator interval and axillary pouch. Conversely, all the patients with contracture had a burning sign in the rotator interval and axillary pouch and no enhancement around the calcium deposit.

Conclusions

Dynamic magnetic resonance imaging identified 2 types of findings in patients with symptomatic chronic calcifying tendinitis: a burning sign in the rotator interval and axillary pouch or around the calcium deposit. The former pattern may be related to shoulder contracture, whereas the latter may be related to impingement sign.

Thirty years of translational research in Mobility Medicine: Collection of abstracts of the 2020 Padua Muscle Days

More than half a century of skeletal muscle research is continuing at Padua University (Italy) under the auspices of the Interdepartmental Research Centre of Myology (CIR-Myo), the European Journal of Translational Myology (EJTM) and recently also with the support of the A&CM-C Foundation for Translational Myology, Padova, Italy. The Volume 30(1), 2020 of the EJTM opens with the collection of abstracts for the conference "2020 Padua Muscle Days: Mobility Medicine 30 years of Translational Research". This is an international conference that will be held between March 18-21, 2020 in Euganei Hills and Padova in Italy. The abstracts are excellent examples of translational research and of the multidimensional approaches that are needed to classify and manage (in both the acute and chronic phases) diseases of Mobility that span from neurologic, metabolic and traumatic syndromes to the biological process of aging. One of the typical aim of Physical Medicine and Rehabilitation is indeed to reduce pain and increase mobility enough to enable impaired persons to walk freely, garden, and drive again. The excellent contents of this Collection of Abstracts reflect the high scientific caliber of researchers and clinicians who are eager to present their results at the PaduaMuscleDays. A series of EJTM Communications will also add to this preliminary evidence.

Assessment of fasted and fed gastrointestinal contraction frequencies in healthy subjects using continuously tagged MRI

BACKGROUND

Continuously tagged MRI during free breathing can assess bowel motility at frequencies as low as the slow wave, motility pattern range. This study aimed to evaluate noninvasive gastrointestinal-tagged MRI for small bowel motility assessment and to observe the physiological response to a 300-kcal meal challenge in healthy, overnight-fasted volunteers.

METHODS

After overnight fasting, 16 healthy subjects (7 women, mean age 25.5, range 19-37 years) underwent a free breathing, tagged MRI scan to capture small bowel motility. Each subject underwent a (a) baseline motility scan, (b) food challenge, (c) postchallenge scan, and (d) second postchallenge scan (after 20 minutes). Motility was quantified using a frequency analysis technique for measuring the spectral power of the strain, referred to as motility score. Motility score was assessed in 20 frequency intervals between 1 and 20 contractions per minute (cpm), and the data were analyzed with linear mixed-effect models.

KEY RESULT

The stimulation protocol demonstrated an immediate, food-induced, motility response in the low-frequency range (2-10 cpm), which is consistent with the stomach and small bowel frequency range (3-12 cpm).

CONCLUSIONS AND INFERENCES

This study shows that this MRI tagging technique is able to quantify the fasted-to-fed response to a 300-kcal meal challenge within the specific small bowel motility frequency range in healthy subjects. The food provocation MRI protocol provides a tool to explore the gut's response to a stimulus in specific motility frequency ranges in patients with gastrointestinal dysmotility and functional disorders.

Application of radial GRAPPA techniques to single- and multislice dynamic speech MRI using a 16-channel neurovascular coil

PURPOSE

To investigate: (1) the feasibility of using through-time radial GeneRalized Autocalibrating Partially Parallel Acquisitions (rGRAPPA) and hybrid radial GRAPPA (h-rGRAPPA) in single- and multislice dynamic speech MRI; (2) whether single-slice dynamic speech MRI at a rate of 15 frames per second (fps) or higher and with adequate image quality can be achieved using these radial GRAPPA techniques.

METHODS

Seven healthy adult volunteers were imaged at 3T using a 16-channel neurovascular coil and 2 spoiled gradient echo sequences (radial trajectory, field of view = 192 × 192 mm² , acquired pixel size = 2.4 × 2.4 mm² ). One sequence imaged a single slice at 16.8 fps, the other imaged 2 interleaved slices at 7.8 fps per slice. Image sets were reconstructed using rGRAPPA and h-rGRAPPA, and their image quality was compared using the root mean square error, structural similarity index, and visual assessments.

RESULTS

Image quality deteriorated when fewer than 170 calibration frames were used in the rGRAPPA reconstruction. rGRAPPA image sets demonstrated: (1) in 97% of cases, a similar image quality to h-rGRAPPA image sets reconstructed using a k-space segment size of 4, (2) in 98% of cases, a better image quality than h-rGRAPPA image sets reconstructed using a k-space segment size of 32.

CONCLUSION

This study confirmed: (1) the feasibility of using rGRAPPA and h-rGRAPPA in single- and multislice dynamic speech MRI, (2) that single-slice speech imaging at a frame rate higher than 15 fps and with adequate image quality can be achieved using these radial GRAPPA techniques.

Efficient Dynamic Parallel MRI Reconstruction for the Low-Rank Plus Sparse Model

The low-rank plus sparse (L+S) decomposition model enables the reconstruction of under-sampled dynamic parallel magnetic resonance imaging (MRI) data. Solving for the low-rank and the sparse components involves non-smooth composite convex optimization, and algorithms for this problem can be categorized into proximal gradient methods and variable splitting methods. This paper investigates new efficient algorithms for both schemes. While current proximal gradient techniques for the L+S model involve the classical iterative soft thresholding algorithm (ISTA), this paper considers two accelerated alternatives, one based on the fast iterative shrinkage-thresholding algorithm (FISTA), and the other with the recent proximal optimized gradient method (POGM). In the augmented Lagrangian (AL) framework, we propose an efficient variable splitting scheme based on the form of the data acquisition operator, leading to simpler computation than the conjugate gradient (CG) approach required by existing AL methods. Numerical results suggest faster convergence of the efficient implementations for both frameworks, with POGM providing the fastest convergence overall and the practical benefit of being free of algorithm tuning parameters.

Detecting the effects of a standardized meal challenge on small bowel motility with MRI in prepared and unprepared bowel

OBJECTIVE

MRI is increasingly used to evaluate small bowel contractility. The objective of this study was to validate a clinically practical stimulation test (300-kcal meal) for small bowel motility.

METHODS

Thirty-one healthy subjects underwent dynamic MRI to capture global small bowel motility after ±10h fasting, of which 15 underwent bowel preparation consisting of 1 L 2.5% mannitol solution and 16 did not. Each subject underwent (1) a baseline motility scan (2) a food challenge (3) a post-challenge scan, and (4) second post-challenge scan (after ±20 minutes). This protocol was repeated within 2 weeks. Motility was quantified using a validated motility assessment technique.

KEY RESULTS

Motility in prepared subjects at baseline was significantly higher than motility in unprepared subjects (0.36 AU vs 0.18 AU, P < 0.001). In the prepared group, the food challenge produced an 8% increase in motility (P = 0.33) while in the unprepared subjects a significant increase of 30% was observed (P < 0.001). Responses to food remained insignificant (P = 0.21) and significant (P = 0.003), for the prepared and unprepared subjects, respectively, ±20 minutes post food challenge. These results were confirmed in the repeated scan session.

CONCLUSION & INFERENCES

A significant response to a 300-kcal meal was measured within 10 minutes in unprepared bowel, supporting the clinical use of this challenge to provoke and assess motility changes. A caloric challenge did not produce an observable increase in motility in mannitol prepared subjects.

Increased perfusion in dynamic gadolinium-enhanced MRI correlates with areas of bone repair and of bone necrosis in patients with Kienböck's disease

BACKGROUND

Osteonecrosis of the lunate, Kienböck's disease, can lead to fragmentation of the lunate, carpal collapse, and severe osteoarthritis. Since the etiology of Kienböck's disease is impaired circulation, a diagnostic method capable of assessing perfusion would be valuable. Recent studies have suggested that dynamic contrast-enhanced (DCE) MR examinations at 3 T can assess perfusion in healthy carpal bones.

PURPOSE

To evaluate the use of DCE-MR for assessing perfusion in the lunate bone in patients with Kienböck's disease. Furthermore, to compare perfusion with histopathology with a focus on bone viability.

STUDY TYPE

Prospective case-control study.

POPULATION

Fourteen patients with Kienböck's disease and a control group of 19 healthy subjects. Field Strength: 3 T with T₁ -weighted fat-saturated contrast-enhanced gradient echo series.

ASSESSMENT

Features of the enhancement curves from the DCE-MR examinations, time to peak (TTP), maximum slope (MS), and maximum enhancement (ME) assessed by a radiologist. Six of 14 patients were surgerized with lunate excision, allowing comparison between features of the enhancement curves and histopathology.

STATISTICAL TESTS

Mann-Whitney U-test. P < 0.05 was considered a statistically significant difference.

RESULTS

Patients with Kienböck's disease showed significantly higher and faster perfusion parameters compared with the control group, the mean value of the TTP in patients was 126.73 sec, in controls 189.79 sec (P = 0.024), ME in patients 173.55 AU, in controls 28.46 AU (P < 0.001), and MS in patients 5.04 AU, in controls 1.06 AU (P < 0.001). When compared with histopathology, increased perfusion was seen in areas of bone formation but also in necrosis. Areas of normal bone showed low perfusion.

DATA CONCLUSION

DCE-MRI at 3 T can diagnose altered perfusion in patients with Kienböck's disease. Increased perfusion cannot definitely be used as a marker of bone viability.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2019;50:481-489.

3D dynamic MRI of the vocal tract during natural speech

PURPOSE

To develop and evaluate a technique for 3D dynamic MRI of the full vocal tract at high temporal resolution during natural speech.

METHODS

We demonstrate 2.4 × 2.4 × 5.8 mm³ spatial resolution, 61-ms temporal resolution, and a 200 × 200 × 70 mm³ FOV. The proposed method uses 3D gradient-echo imaging with a custom upper-airway coil, a minimum-phase slab excitation, stack-of-spirals readout, pseudo golden-angle view order in k_x -k_y , linear Cartesian order along k_z , and spatiotemporal finite difference constrained reconstruction, with 13-fold acceleration. This technique is evaluated using in vivo vocal tract airway data from 2 healthy subjects acquired at 1.5T scanner, 1 with synchronized audio, with 2 tasks during production of natural speech, and via comparison with interleaved multislice 2D dynamic MRI.

RESULTS

This technique captured known dynamics of vocal tract articulators during natural speech tasks including tongue gestures during the production of consonants "s" and "l" and of consonant-vowel syllables, and was additionally consistent with 2D dynamic MRI. Coordination of lingual (tongue) movements for consonants is demonstrated via volume-of-interest analysis. Vocal tract area function dynamics revealed critical lingual constriction events along the length of the vocal tract for consonants and vowels.

CONCLUSION

We demonstrate feasibility of 3D dynamic MRI of the full vocal tract, with spatiotemporal resolution adequate to visualize lingual movements for consonants and vocal tact shaping during natural productions of consonant-vowel syllables, without requiring multiple repetitions.

RECOVERY OF POINT CLOUDS ON SURFACES: APPLICATION TO IMAGE RECONSTRUCTION

We introduce a framework for the recovery of points on a smooth surface in high-dimensional space, with application to dynamic imaging. We assume the surface to be the zero-level set of a bandlimited function. We show that the exponential maps of the points on the surface satisfy annihilation relations, implying that they lie in a finite dimensional subspace. We rely on nuclear norm minimization of the maps to recover the points from noisy and undersampled measurements. Since this direct approach suffers from the curse of dimensionality, we introduce an iterative reweighted algorithm that uses the "kernel trick". The resulting algorithm has similarities to iterative algorithms used in graph signal processing (GSP); this framework can be seen as a continuous domain alternative to discrete GSP theory. The use of the algorithm in recovering free breathing and ungated cardiac data shows the potential of this framework in practical applications.

Strain Map of the Tongue in Normal and ALS Speech Patterns from Tagged and Diffusion MRI

Amyotrophic Lateral Sclerosis (ALS) is a neurological disease that causes death of neurons controlling muscle movements. Loss of speech and swallowing functions is a major impact due to degeneration of the tongue muscles. In speech studies using magnetic resonance (MR) techniques, diffusion tensor imaging (DTI) is used to capture internal tongue muscle fiber structures in three-dimensions (3D) in a non-invasive manner. Tagged magnetic resonance images (tMRI) are used to record tongue motion during speech. In this work, we aim to combine information obtained with both MR imaging techniques to compare the functionality characteristics of the tongue between normal and ALS subjects. We first extracted 3D motion of the tongue using tMRI from fourteen normal subjects in speech. The estimated motion sequences were then warped using diffeomorphic registration into the b0 spaces of the DTI data of two normal subjects and an ALS patient. We then constructed motion atlases by averaging all warped motion fields in each b0 space, and computed strain in the line of action along the muscle fiber directions provided by tractography. Strain in line with the fiber directions provides a quantitative map of the potential active region of the tongue during speech. Comparison between normal and ALS subjects explores the changing volume of compressing tongue tissues in speech facing the situation of muscle degradation. The proposed framework provides for the first time a dynamic map of contracting fibers in ALS speech patterns, and has the potential to provide more insight into the detrimental effects of ALS on speech.

Evaluation of gastrointestinal motility with MRI: Advances, challenges and opportunities

Dynamic magnetic resonance imaging (MRI) of gastrointestinal motility has developed rapidly over the past few years. The non-invasive and non-ionizing character of MRI is an important advantage together with the fact that it is fast and can visualize the entire gastrointestinal tract. Advances in imaging and quantification techniques have facilitated assessment of gastric, small intestinal, and colonic motility in a clinical setting. Automated quantitative motility assessment using dynamic MRI meets the need for non-invasive techniques. Recently, studies have begun to examine this technique in patients, including those with IBD, pseudo-obstruction and functional bowel disorders. Remaining challenges for clinical implementation are processing the large amount of data, standardization and validation of the numerous MRI metrics and subsequently assessment of the potential role of dynamic MRI. This review examines the methods, advances, and remaining challenges of evaluation of gastrointestinal motility with MRI. It accompanies an article by Khalaf et al. in this journal that describes a new protocol for assessment of pan-intestinal motility in fasted and fed state in a single MRI session.

Accelerated 4D self-gated MRI of tibiofemoral kinematics

Anatomical (static) magnetic resonance imaging (MRI) is the most useful imaging technique for the evaluation and assessment of internal derangement of the knee, but does not provide dynamic information and does not allow the study of the interaction of the different tissues during motion. As knee pain is often only experienced during dynamic tasks, the ability to obtain four-dimensional (4D) images of the knee during motion could improve the diagnosis and provide a deeper understanding of the knee joint. In this work, we present a novel approach for dynamic, high-resolution, 4D imaging of the freely moving knee without the need for external triggering. The dominant knee of five healthy volunteers was scanned during a flexion/extension task. To evaluate the effects of non-uniform motion and poor coordination skills on the quality of the reconstructed images, we performed a comparison between fully free movement and movement instructed by a visual cue. The trigger signal for self-gating was extracted using principal component analysis (PCA), and the images were reconstructed using a parallel imaging and compressed sensing reconstruction pipeline. The reconstructed 4D movies were scored for image quality and used to derive bone kinematics through image registration. Using our method, we were able to obtain 4D high-resolution movies of the knee without the need for external triggering hardware. The movies obtained with and without instruction did not differ significantly in terms of image scoring and quantitative values for tibiofemoral kinematics. Our method showed to be robust for the extraction of the self-gating signal even for uninstructed motion. This can make the technique suitable for patients who, as a result of pain, may find it difficult to comply exactly with instructions. Furthermore, bone kinematics can be derived from accelerated MRI without the need for additional hardware for triggering.

Fast dynamic ventilation MRI of hyperpolarized 129 Xe using spiral imaging

Purpose

To develop and optimize a rapid dynamic hyperpolarized ¹²⁹Xe ventilation (DXeV) MRI protocol and investigate the feasibility of capturing pulmonary signal‐time curves in human lungs.

Theory and Methods

Spiral k‐space trajectories were designed with the number of interleaves N_int = 1, 2, 4, and 8 corresponding to voxel sizes of 8 mm, 5 mm, 4 mm, and 2.5 mm, respectively, for field of view = 15 cm. DXeV images were acquired from a gas‐flow phantom to investigate the ability of N_int = 1, 2, 4, and 8 to capture signal‐time curves. A finite element model was constructed to investigate gas‐flow dynamics corroborating the experimental signal‐time curves. DXeV images were also carried out in six subjects (three healthy and three chronic obstructive pulmonary disease subjects).

Results

DXeV images and numerical modelling of signal‐time curves permitted the quantification of temporal and spatial resolutions for different numbers of spiral interleaves. The two‐interleaved spiral (N_int = 2) was found to be the most time‐efficient to obtain DXeV images and signal‐time curves of whole lungs with a temporal resolution of 624 ms for 13 slices. Signal‐time curves were well matched in three healthy volunteers. The Spearman's correlations of chronic obstructive pulmonary disease subjects were statistically different from three healthy subjects (P < 0.05).

Conclusion

The N_int = 2 spiral demonstrates the successful acquisition of DXeV images and signal‐time curves in healthy subjects and chronic obstructive pulmonary disease patients. Magn Reson Med 79:2597–2606, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

A Four-dimensional Motion Field Atlas of the Tongue from Tagged and Cine Magnetic Resonance Imaging

Representation of human tongue motion using three-dimensional vector fields over time can be used to better understand tongue function during speech, swallowing, and other lingual behaviors. To characterize the inter-subject variability of the tongue's shape and motion of a population carrying out one of these functions it is desirable to build a statistical model of the four-dimensional (4D) tongue. In this paper, we propose a method to construct a spatio-temporal atlas of tongue motion using magnetic resonance (MR) images acquired from fourteen healthy human subjects. First, cine MR images revealing the anatomical features of the tongue are used to construct a 4D intensity image atlas. Second, tagged MR images acquired to capture internal motion are used to compute a dense motion field at each time frame using a phase-based motion tracking method. Third, motion fields from each subject are pulled back to the cine atlas space using the deformation fields computed during the cine atlas construction. Finally, a spatio-temporal motion field atlas is created to show a sequence of mean motion fields and their inter-subject variation. The quality of the atlas was evaluated by deforming cine images in the atlas space. Comparison between deformed and original cine images showed high correspondence. The proposed method provides a quantitative representation to observe the commonality and variability of the tongue motion field for the first time, and shows potential in evaluation of common properties such as strains and other tensors based on motion fields.

Parotid malignancies are not homogenous in terms of their presentation on dynamic MRI - a preliminary report

CONCLUSION

At least two groups of parotid malignancies exist, including one whose parameters of dynamic MRI closely resemble those of pleomorphic adenomas. Also tumors with long time to peak enhancement after administration of a contrast agent (T_peak) and low washout rate (WR) should be considered malignant, especially single masses with concomitant lymphadenopathy located within the deep parotid lobe.

OBJECTIVE

To verify if malignant tumors of the parotid are homogeneous in terms of parameters of preoperative dynamic MRI: T_peak and WR.

METHODS

The retrospective analysis included 221 surgical patients with parotid tumors. Aside from fine needle biopsy, pre-operative examination included dynamic and diffusion-weighted MRI. Final diagnosis was based on histopathological examination of the surgical specimen.

RESULTS

Twenty-four of 221 (10.8%) malignant lesions were identified. Using k-means clustering based on T_peak and WR values, two distinct clusters of parotid malignancies were identified. The cut-off value for T_peak optimally differentiating between the clusters was 140 s; the cut-off value for WR could not be identified. The two clusters did not differ in terms of dynamic and diffusion-weighted MRI parameters, patient age, sex and prevalence of lymphadenopathy. Significant inter-cluster differences were found in the prevalence of deep parotid lobe involvement and presence of a single mass.

Characteristics of dynamic magnetic resonance imaging of idiopathic severe frozen shoulder

BACKGROUND

The purpose of this study was to evaluate the effectiveness of current techniques for dynamic 3-dimensional magnetic resonance imaging (MRI) in the diagnosis of idiopathic severe frozen shoulder (FS).

MATERIALS AND METHODS

Subjects consisted of 5 healthy volunteers and 16 patients with idiopathic severe FS. We defined severe idiopathic FS as follows: range of motion ≤100° in forward flexion, ≤10° in external rotation, and ≤L5 in internal rotation. All patients suffered from continued global range of motion loss for at least 6 months. We evaluated the diagnostic characteristics of 3-dimensional dynamic MRI in FS patients compared with those in healthy volunteers.

RESULTS

MRI of all FS patients displayed an abnormal intake of blood flow from the acromial arterial network and the branches of circumflex humeral arteries into the axillary pouch and the rotator interval. We named this finding "burning sign." The burning sign was present at all phases of the condition. In the FS group, the patients with enhanced deposition of contrast medium in the axillary pouch in the delayed phase (n = 11) had a statistically significant score for pain during exercise, higher than that of patients with reduced deposition of contrast medium at the same site (n = 5; P = .027).

CONCLUSION

Burning sign is an abnormal finding that appears in dynamic MRI of severe FS. Hence, the burning sign may be associated with pain and inflammation in idiopathic FS.

A continuous-infusion dynamic MRI model at 3.0 Tesla for the serial quantitative evaluation of microvascular proliferation in an animal model of glioblastoma multiforme

PURPOSE

To develop a continuous-infusion dynamic MRI technique to characterize tumor-associated microvascular proliferation (MVP) in a rat brain model of glioblastoma multiforme.

METHODS

The proposed model assumes effects due to tumor-associated MVP (eg, vascular permeability, K^trans ; intravascular plasma fraction, v_p ) cannot be individually separated and solves for a single parameter (k^vasc ) that quantifies the T₁ -weighted contrast enhancement from dynamic images acquired during continuous contrast agent (CA) infusion. Untreated C6 tumor-bearing animals (N = 6) were serially imaged on postoperative days (PODs) 14 and 18 with a 3 Tesla clinical scanner utilizing a dynamic spatial and temporal resolution of 0.38 × 0.38 × 1.5 mm³ and 3.47 s, respectively.

RESULTS

An association was present between PODs 14 and 18 for median tumor k^vasc (Pearson's r = 0.94, P = 0.0052) and CA concentration ([CA], derived from pre- and postcontrast R₁ maps; r = 0.94, P = 0.0054). On POD 18, there was a voxel-based association between k^vasc and [CA] within each tumor (0.45 < r < 0.82, P < 0.001). However, voxel-based subregions demonstrated a reduced association between k^vasc and [CA] (N = 5; -0.08 < r < 0.22, P > 0.05) or an inverse association (N = 1; r = -0.28, P = 0.001), indicating differences between locations of vascular permeability and subsequent CA pooling in tumors.

CONCLUSION

The continuous-infusion method may provide a quantitative measure for characterizing and monitoring tumor-associated MVP. Magn Reson Med 78:1824-1838, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Prevalence of dynamic magnetic resonance imaging-identified pelvic organ prolapse in pre- and postmenopausal women without clinically evident pelvic organ descent

Background Dynamic magnetic resonance imaging (dMRI) is an imaging tool that can be used to evaluate and stage pelvic organ prolapse (POP). Greater understanding of the incidental detection of POP in asymptomatic patients is needed. Purpose To evaluate the prevalence of dMRI-detected POP in pre-and postmenopausal women who were imaged for reasons unrelated to pelvic floor dysfunction. Material and Methods A total of 227 women who had diagnoses that did not include POP underwent abdominal/pelvic dMRI. Patients with a positive gynecological examination for or a clinical history of POP ( n = 11), hysterectomy ( n = 4), or gynecologic-oncology surgery ( n = 2) were excluded, as well as patients who were unable to strain during MRI ( n = 11). A total of 199 patients without visible prolapse were enrolled in the study. An H-line, M-line, pubococcygeal line (PCL), and mid-pubic line (MPL) were used to detect and grade prolapse. Results The prevalence of dMRI-identified POP was higher in postmenopausal subjects. The PCL led to a greater frequency of prolapse detection than the MPL. The frequency of middle compartment descent was similar regardless of whether the PCL or MPL was used as a reference line. There was a higher incidence of prolapse in the posterior compartment. Using an H-line and PCL as references, the anterior and posterior compartments were found to significantly differ between pre- and postmenopausal subjects. The MRI parameters that were used to define POP were not correlated with parity, vaginal birth, BMI, or fetal birth weight. With respect to the MPL, age was correlated with both the presence of an elongated H-line and with descent. Conclusion Dynamic MRI identified incidental pelvic organ prolapse in asymptomatic patients. The prevalence of dMRI-detected POP was higher in postmenopausal women without visible prolapse. These findings suggest the need for further studies to identify how to modify the currently used dMRI thresholds for postmenopausal women.

ACCELERATED DYNAMIC MRI USING STRUCTURED LOW RANK MATRIX COMPLETION

We introduce a fast structured low-rank matrix completion algorithm with low memory & computational demand to recover the dynamic MRI data from undersampled measurements. The 3-D dataset is modeled as a piecewise smooth signal, whose discontinuities are localized to the zero sets of a bandlimited function. We show that a structured matrix corresponding to convolution with the Fourier coefficients of the signal derivatives is highly low-rank. This property enables us to recover the signal from undersampled measurements. The application of this scheme in dynamic MRI shows significant improvement over state of the art methods.

Computational fluid dynamics upper airway effective compliance, critical closing pressure, and obstructive sleep apnea severity in obese adolescent girls

Obstructive sleep apnea syndrome (OSAS) is associated with anatomical abnormalities restricting upper airway size and functional factors decreasing pharyngeal dilator activity in sleep. In this study we hypothesized that OSAS is also associated with altered pharyngeal mechanical compliance during wakefulness. Five OSAS and six control obese girls between 14 and 18 years of age were studied. All underwent polysomnography, critical closing pressure (P_crit) studies, and dynamic MRI of the upper airway during awake tidal breathing. Effective airway compliance was defined as the slope of cross-sectional area vs. average pressure between maximum inspiration and maximum expiration along the pharyngeal airway. Pharyngeal pressure fields were calculated by using image-based computational fluid dynamics and nasal resistance. Spearman correlations were calculated to test associations between apnea-hypopnea index (AHI), P_crit, and airway compliance. Effective compliances in the nasopharynx (C_NP) and velopharynx (C_VP) were lower and negative in OSAS compared with controls: -4.4 vs. 1.9 (mm²/cmH₂O, P = 0.012) and -2.1 vs. 3.9 (mm²/cmH₂O, P = 0.021), respectively, suggesting a strong phasic pharyngeal dilator activity during inspiration in OSAS compared with controls. For all subjects, C_NP and AHI correlated negatively (r_S = -0.69, P = 0.02), and passive P_crit correlated with C_NP (r_S = -0.76, P = 0.006) and with AHI (r_S = 0.86, P = 0.0006). Pharyngeal mechanics obtained during wakefulness could be used to characterize subjects with OSAS. Moreover, negative effective compliance during wakefulness and its correlation to AHI and P_crit suggest that phasic dilator activity of the upper pharynx compensates for negative pressure loads in these subjects.

Evaluation of transcatheter arterial embolization therapy on hepatocellular carcinomas using contrast-enhanced harmonic power Doppler sonography: comparison with CT, power Doppler sonography, and dynamic MRI

PURPOSE

The aim of this study was to assess and compare the sensitivity of power Doppler sonography, contrast-enhanced sonography, plain computed tomography (CT), and dynamic magnetic resonance imaging (MRI) for detecting hepatocellular carcinoma (HCC) nodules incompletely treated with transcatheter arterial embolization (TAE).

METHODS

A total of 63 unresectable HCC nodules were examined in this study. The HCCs were treated with TAE. All patients underwent plain CT, power Doppler sonography, contrast-enhanced harmonic power Doppler sonography, and dynamic MRI 1 week after TAE. The sensitivity of each modality to incompletely treated HCC nodules was compared. Detection of the residual viable HCC on angiography or tumor biopsy was regarded as the gold standard for the diagnosis of incomplete treatment.

RESULTS

Twenty-four nodules (38%) were diagnosed as incompletely treated. The sensitivities of plain CT, power Doppler sonography, contrast-enhanced harmonic power Doppler sonography, and dynamic MRI to these incompletely treated nodules were 42% (10/24), 46% (11/24), 88% (21/24), and 79% (19/24), respectively. Eighty percent (19 nodules) of the 24 incompletely treated nodules were located within a depth of less than 8 cm. The sensitivities of plain CT, power Doppler sonography, contrast-enhanced harmonic power Doppler sonography, and dynamic MRI to these superficial incompletely treated nodules were 37% (7/19), 53% (10/19), 100% (19/19), and 74% (14/19), respectively. In contrast, the sensitivities of each modality to deeply located nodules were 60% (3/5), 20% (1/5), 40% (2/5), and 100% (5/5), respectively.

CONCLUSION

Plain CT and power Doppler sonography had a low sensitivity to HCC nodules incompletely treated with TAE. Except for those that were deeply located, contrast-enhanced harmonic sonography showed the highest sensitivity in detecting incompletely treated HCC nodules.