Yin, Ziying

Abstract Glioblastoma progression is spatially heterogeneous, but conventional imaging provides limited information about where subsequent tumor progression is likely to occur. We developed a directional magnetic resonance elastography (MRE) framework to test whether local post-treatment tumor–brain interface mechanics are associated with later spatial tumor progression. In a secondary analysis of a prospectively acquired glioblastoma cohort, wedge-level viscoelastic instability features were extracted from the first post-treatment MRE scan and related to novel tumor burden on the second post-treatment scan after excluding tumor already present on pretreatment or first post-treatment imaging. Nine patients had longitudinal imaging suitable for spatial comparison; six lesions showed net interval growth and were included in the primary wedge-level directional analysis, while three non-growing lesions were retained for descriptive comparison. In growing lesions, several directional mechanical features were descriptively associated with later novel tumor burden. In cluster-aware models accounting for within-patient dependence among wedges, mean Δtanδ showed the most consistent association with later wedge-level novel tumor fraction across mixed-effects and generalized estimating equation analyses. Associations were directionally stable across wedge-width sensitivity analyses. These findings provide proof of principle that post-treatment glioblastoma interface mechanics contain spatially resolved information related to where later tumor emergence occurs, supporting further validation of directional MRE as a framework for longitudinal mapping of progression geometry.

Abstract Mechanical properties of biological tissues vary across spatial scales, yet radiomics typically relies on fixed, heuristic choices for neighbourhood size, kernel geometry, and spectral content - choices that can silently reshape the feature space before any modelling begins. We introduce a label-free, information-theoretic framework for selecting extraction parameters in multi-frequency MRE radiomics. For each configuration θ - neighbourhood radius r, kernel geometry k (sphere or shell), and frequency subset f - we extract a radiomics feature matrix and score it using an objective J(θ) that integrates distributional richness (Shannon entropy), cross-frequency coherence (canonical correlation), inter-feature redundancy (Spearman correlation), and bootstrap stability. We evaluate 121 configurations per tissue in multi-frequency MRE (30-60 Hz) of human brain, liver, and a calibrated phantom, and test robustness using 10,000 Dirichlet-sampled objective weightings. Across tissues, neighbourhood aggregation is consistently preferred over voxel-wise extraction, outperforming the no-neighbourhood baseline in 98.4-100% of weightings. External validation in 100 independent brain scans acquired with a different protocol and wider frequency range (20-90 Hz) confirms a reproducible mesoscopic plateau at r = 3-5 (9-15 mm), with a modal optimum at r = 4; omitting neighbourhood analysis reduces J(θ) by 38% relative to each subject’s optimum. Frequency-subset preferences replicate across datasets, with lower frequencies most frequently selected for brain. By turning ad hoc extraction choices into an outcome-free optimisation step, this framework improves reproducibility, reduces sensitivity to heuristic parameter choices, and generalises across acquisition protocols and imaging sites.

Abstract Background Conventional MRI cannot reliably distinguish radiation necrosis (RN) from recurrent metastasis after cranial radiotherapy, as both can show similar enhancement despite different biology. We tested whether these entities are mechanically non-equivalent in vivo and separable by MRE-derived viscoelastic metrics and perilesional interface-instability features. Methods In a prospective, histopathology-anchored cohort, 11 post-radiotherapy enhancing lesions were classified as RN (n=3) or recurrent/progressive tumor (n=8). MRE was acquired at 3.0 T with single-frequency 60-Hz excitation to derive storage modulus (G′), loss modulus (G″), and complex shear modulus magnitude (|G*|). Co-primary endpoints were median tumor G′ and |G*|, each tested one-sided (RN > tumor) with Holm correction across the two co-primary tests. Median tumor G″ was tested two-sided. A prespecified secondary 6-endpoint family (absolute and tumor/NAWM-normalized G′, G″, and |G*|) was analyzed with Benjamini–Hochberg FDR control. Exploratory instability mapping in a 0–6 mm peritumoral shell generated interface-topology metrics, including convexity. Results Absolute tumor-core medians were higher in RN than tumor for |G*| (1.79 vs 1.32 kPa; Cliff’s δ=0.67; q=0.10), G′ (1.62 vs 1.09 kPa; δ=0.50; q=0.14), and G″ (0.81 vs 0.46 kPa; δ=0.75; q=0.10). NAWM normalization improved separation: tumor/NAWM |G*| (2.26 vs 1.41; δ=0.92; q=0.04) and tumor/NAWM G″ (2.67 vs 0.87; δ=1.00; q=0.04) were FDR-significant. Convexity also differentiated RN from tumor (0.49 vs 0.36; δ=1.00; MWU p=0.01). Conclusions Tumor/NAWM G″, tumor/NAWM |G*|, convexity, and tumor G″ emerged as the strongest candidate features, indicating that RN is mechanically harder and more dissipative than recurrent metastasis. Signal strength was high (Cliff’s δ up to 1.00) but should be interpreted cautiously given sample size. Exploratory analyses further suggest that instability mapping captures biologically relevant interface behavior. These findings support a mechanics-based RN-versus-recurrence framework and justify prespecified, preregistered external validation.

ABSTRACT Purpose Active lesion expansion is clinically important because it can indicate disease progression and influence management and treatment planning. However, sequential imaging is often unavailable, and rapid growth can only be inferred from indirect findings such as active perilesional tissue injury. This study aimed to implement and test the feasibility of an MR elastography (MRE)‐based technique for assessing the expansile nature of mass lesions by detecting stiffness changes caused by mechanical strain manifestations in adjacent tissues. Methods A tissue strain mapping (TSM) algorithm was developed using spatial–temporal directional filtering to reveal the presence of perilesional latent strain that is not apparent in conventional MRE inversions. The technical feasibility was tested in tissue‐simulating phantoms containing a mass that was progressively expanded. Preliminary feasibility for clinical application was also demonstrated in representative in vivo examples. Resulting TSM maps were compared with conventional MRE stiffness maps. Results In the phantom experiments, TSM revealed the presence of local strain at the periphery of the expanding mass, with the effect increasing proportionally to the degree of expansion. Conventional MRE did not demonstrate comparable localized stiffening patterns. In vivo examples showed similar perilesional strain features that were not apparent on conventional MRE maps. Conclusion The feasibility of the proposed MRE‐based TSM technique was demonstrated in phantom and preliminary in vivo studies. These findings suggest that the technique is a promising tool for assessing the expansile nature of mass lesions and motivate further development of the TSM processing and exploration of potential clinical applications.

Background: Accurate characterization of tumor stiffness and adhesion is crucial when surgically resecting meningiomas. Currently, such determinations are difficult to ascertain through standard preoperative imaging diagnostics and clinical experience. Herein, we explored magnetic resonance elastography (MRE) as a means of assessing these parameters and thereby improving neurosurgical resections. Methods: This prospective study of 193 patients with meningiomas involved three distinct patient subsets (cohort 1, 88 [2017–2021]; cohort 2, 55 [2022–2023]; cohort 3, 50 [2023–2024]. We analyzed preoperatively obtained MRE data on meningioma stiffness (cohorts 1 and 2) and adhesion (cohort 2) for accuracy, while matching subjects with or without MRE studies (cohort 3, 1:1 ratio) to define the role of preoperative MRE in improving operative outcomes. Results: Areas under the curve (AUCs) for MRE in gauging meningioma stiffness and adhesion were 0.87 and 0.89, respectively. Both parameters significantly surpassed empiric predictions of stiffness and adhesion rendered by surgeons and radiologists, based on routine magnetic resonance imaging (MRI). MRE-determined stiffness measures independently correlated with prolonged resection duration, whereas MRE-assessed adhesion showed independent associations with increased blood loss, longer duration of tumor resection, and higher risk of pial injury. Ultimately, results based on cohort 3 demonstrated that MRE has great potential in this setting to improve neurosurgical procedures by significantly reducing intraoperative bleeding, resection duration, and risks of pial injury. Conclusion: MRE is a promising, noninvasive tool for accurately assessing meningioma stiffness and adhesion, standing to improve the precision and safety of needed surgical intervention.

ABSTRACT Purpose Abnormal restrictions at myofascial interfaces are a key factor in chronic low back pain (CLBP) with myofascial pain syndrome (MPS), but no established imaging method exists for quantitative assessment. This pilot study evaluated the clinical feasibility of MR elastography (MRE)‐based slip interface imaging (SII) for visualizing myofascial mobility in CLBP‐MPS patients versus age−/sex‐/BMI‐matched healthy participants. A new SII biomarker, maximal normalized displacement discontinuity (), was introduced, and the previously developed normalized octahedral shear strain (NOSS) map was used for comparison. Preliminary assessments of vibration‐frequency effects, repeatability, and group differences in inter‐muscular mobility were conducted. Methods MRE was conducted at 30, 60, and 90 Hz to compare with NOSS for visualizing inter‐muscular slip interfaces at the quadratus lumborum‐erector spinae (QL‐ES) and erector spinae‐multifidus (ES‐M) muscle boundaries at L3‐L4. Repeatability of and NOSS was assessed at 30 Hz (identified as the best‐performing frequency) using within‐subject coefficients of variation and a linear mixed‐effects model from within‐ and between‐day scans. Clinical feasibility was then assessed via Wilcoxon rank‐sum tests in six CLBP‐MPS patients and matched controls. Results at 30 Hz exhibited superior delineation of slip interfaces compared to higher frequencies and NOSS. demonstrated significantly lower variation and better repeatability than NOSS. Healthy participants showed continuous, well‐defined slip interfaces, whereas CLBP‐MPS patients displayed disrupted interfaces with significant mobility reduction at L4 ES‐M. Conclusion SII with lower vibration frequency and map enhances visualization of lower back myofascial mobility and shows potential for assessing differences between healthy individuals and CLBP‐MPS patients.

AbstractPurposeWhile standard MR elastography (MRE) uses harmonic mechanical waves, there are some applications in which imaging and analysis of transient mechanical motion are of interest. Wavelet MRE has been developed for detecting broadband motion from transient excitation. The goal of this study was to evaluate the accuracy and efficiency of wavelet MRE for transient displacement detection in brain MRE applications.MethodsTransient motion was induced in a gel phantom, while MRE images were acquired using bipolar motion‐encoding gradient (MEG) at multiple frequencies (20–200 Hz). Displacements were estimated using (1) combinations of multiple MEGs forming the wavelet MRE and (2) deconvolution from a single MEG. These estimated displacements were used to calculate the MRE phase for each MEG. Correlation (r2) between the calculated and acquired phases was evaluated.Three healthy volunteers were scanned twice in a clinical scanner using wavelet MRE with an occipital impact. Time‐resolved brain translation, rotation, and maximal principal strain were calculated. Repeatability was assessed both qualitatively and through Pearson correlation.ResultsWavelet MRE outperformed standard MRE in displacement estimation, showing higher correlations between calculated and acquired phase, even with fewer phase offsets. In the volunteer study, consistent temporal motion dynamics and spatial maximal principal strain distributions across repeated scans demonstrated the repeatability of wavelet MRE.ConclusionThis study validated the accuracy and efficiency of wavelet MRE for broadband motion detection and demonstrated its feasibility and repeatability in vivo. This technique shows promise for advancing our understanding of the injury risks and mechanisms associated with sports‐related head trauma.

ABSTRACTBackgroundMagnetic resonance elastography (MRE) quantifies tissue viscoelasticity, offering insights into intracranial meningiomas. MRE‐derived parameters, including the storage modulus (G') and the magnitude of the complex modulus (|G*|), may aid neurosurgical planning, but their repeatability remains unexplored in brain tumors.Purpose/HypothesisTo evaluate the repeatability of G' and |G*| in meningiomas, the influence of tumor volume and spatial location on variability, and ranking stability across paired elastograms.Study TypeProspective.Population/SubjectsSeventeen paired MRE scans from 16 patients with meningiomas (mean age 64 ± 12 years, 10 females) scheduled for resection.Field Strength/SequenceMRE was performed on a 3 T system using a single‐shot spin‐echo echo‐planar imaging (SE‐EPI) sequence.AssessmentTumors were segmented on contrast‐enhanced T1‐weighted images; distance and volume metrics were extracted. Regions of interest were applied to G' and |G*| elastograms. Repeatability was evaluated using the coefficient of variation (CV%). Ranking stability was assessed by comparing median G' or |G*| values.Statistical TestsThe difference in CV% was assessed using the Wilcoxon signed rank test. Linear regression assessed the effect of tumor volume and distance from the external occipital protuberance (EOP) on CV%. The Wilcoxon signed‐rank test evaluated ranking stability and scan‐rescan CV%. Statistical significance was defined as p < 0.05.ResultsThe difference in CV% G' (14.89 ± 12.32) and CV% |G*| (13.88 ± 12.60) was not significant (p = 0.61). Tumor volume inversely correlated with CV% (G': β = −0.75 |G*|: β = −0.59). Tumor‐EOP distance had no significant effect (G': p = 0.58, |G*|: p = 0.23). Rankings remained stable (G': p = 0.82, |G*|: p = 0.71).Data ConclusionMRE‐derived parameters are repeatable in intracranial meningiomas, with G' and |G*| being equally stable. Larger tumors yielded more consistent estimates. Preserved ranking stability supports the use of MRE in preoperative planning.Evidence Level2.Technical EfficacyStage 2.

Abstract Purpose To evaluate the test–retest repeatability of a rapid, free‐breathing two‐dimensional (2D) MR elastography (MRE) technique and to assess the reliability of the liver‐stiffness measurements compared with conventional breath‐hold MREs. Methods Fifteen and 115 participants were enrolled in the technical repeatability and measurement equivalence assessment cohorts, respectively. All participants underwent rapid free‐breathing and conventional breath‐hold 2D MRE (twice in repeatability cohort) on 1.5T scanners. Both methods have four phase offsets over one harmonic motion cycle at 60 Hz; one and 10 cycles were collected and processed in breath‐hold and free‐breathing MREs, respectively. The liver stiffness measurements of free‐breathing (LS F ) and breath‐hold MRE (LS B ) were calculated from manually drawn regions of interest. The repeatability coefficients and Spearman correlation were used to assess technical repeatability and measurement agreement between LS F and LS B . Univariable and multivariable linear regressions were performed to identify potential influencing factors, including age, sex, body mass index, and fat fraction, in the measurement agreement between LS B and LS F . Results The repeatability coefficient of free‐breathing 2D MRE is comparable to breath‐hold MRE (LS F : 20.8%, LS B : 20.4%). LS F showed a strong agreement and significant correlation with LS B (LS F = 1.01 × LS B , ρ = 0.94, p < 0.001). The measurement agreement between LS F and LS B was only significantly affected by sex ( p = 0.047) after adjusting for confounding factors of age, body mass index, and fat fraction. Conclusions The nongated, free‐breathing, multislice 2D MRE technique can provide reliable liver stiffness measurements compared with conventional breath‐hold 2D MRE. It could provide a comfortable alternative method with reliable liver stiffness measurements for patients who have difficulty in suspending respiration.