Publications

2021

Godino-Moya A, on-Lara R-M \ia M, andez MM \in-F, Prieto C, opez CA-L. Elastic AlignedSENSE for Dynamic MR Reconstruction: A Proof of Concept in Cardiac Cine. Entropy (Basel). 2021;23(5):555. doi:10.3390/e23050555
Numerous methods in the extensive literature on magnetic resonance imaging (MRI) reconstruction exploit temporal redundancy to accelerate cardiac cine. Some of them include motion compensation, which involves high computational costs and long runtimes. In this work, we proposed a method-elastic alignedSENSE (EAS)-for the direct reconstruction of a motion-free image plus a set of nonrigid deformations to reconstruct a 2D cardiac sequence. The feasibility of the proposed approach was tested in 2D Cartesian and golden radial multi-coil breath-hold cardiac cine acquisitions. The proposed approach was compared against parallel imaging compressed sense (sPICS) and group-wise motion corrected compressed sense (GWCS) reconstructions. EAS provides better results on objective measures with considerable less runtime when an acceleration factor is higher than 10×. Subjective assessment of an expert, however, invited proposing the combination of EAS and GWCS as a preferable alternative to GWCS or EAS in isolation.
Li S, Zheng Y, Sun W, Lasič S, Szczepankiewicz F, Wei Q, Han S, Zhang S, Zhong X, Wang L, et al. Glioma Grading, Molecular Feature Classification, and Microstructural Characterization Using MR Diffusional Variance Decomposition (DIVIDE) Imaging. Eur Radiol. 2021;31(11):8197–8207. doi:10.1007/s00330-021-07959-x
OBJECTIVE: To evaluate the potential of diffusional variance decomposition (DIVIDE) for grading, molecular feature classification, and microstructural characterization of gliomas. MATERIALS AND METHODS: Participants with suspected gliomas underwent DIVIDE imaging, yielding parameter maps of fractional anisotropy (FA), mean diffusivity (MD), anisotropic mean kurtosis (MKA), isotropic mean kurtosis (MKI), total mean kurtosis (MKT), MKA/MKT, and microscopic fractional anisotropy (μFA). Tumor type and grade, isocitrate dehydrogenase (IDH) 1/2 mutant status, and the Ki-67 labeling index (Ki-67 LI) were determined after surgery. Statistical analysis included 33 high-grade gliomas (HGG) and 17 low-grade gliomas (LGG). Tumor diffusion metrics were compared between HGG and LGG, among grades, and between wild and mutated IDH types using appropriate tests according to normality assessment results. Receiver operating characteristic and Spearman correlation analysis were also used for statistical evaluations. RESULTS: FA, MD, MKA, MKI, MKT, μFA, and MKA/MKT differed between HGG and LGG (FA: p = 0.047; MD: p = 0.037, others p < 0.001), and among glioma grade II, III, and IV (FA: p = 0.048; MD: p = 0.038, others p < 0.001). All diffusion metrics differed between wild-type and mutated IDH tumors (MKI: p = 0.003; others: p < 0.001). The metrics that best discriminated between HGG and LGGs and between wild-type and mutated IDH tumors were MKT and FA respectively (area under the curve 0.866 and 0.881). All diffusion metrics except FA showed significant correlation with Ki-67 LI, and MKI had the highest correlation coefficient (rs = 0.618). CONCLUSION: DIVIDE is a promising technique for glioma characterization and diagnosis. KEY POINTS: • DIVIDE metrics MKI is related to cell density heterogeneity while MKA and μFA are related to cell eccentricity. • DIVIDE metrics can effectively differentiate LGG from HGG and IDH mutation from wild-type tumor, and showed significant correlation with the Ki-67 labeling index. • MKI was larger than MKA which indicates predominant cell density heterogeneity in gliomas. • MKA and MKI increased with grade or degree of malignancy, however with a relatively larger increase in the cell eccentricity metric MKA in relation to the cell density heterogeneity metric MKI.
Eriksson J, Nelson D, Holst A, Hellgren E, Friman O, Oldner A. Temporal patterns of organ dysfunction after severe trauma. Crit Care. 2021;25(1):165. doi:10.1186/s13054-021-03586-6
BACKGROUND: Understanding temporal patterns of organ dysfunction (OD) may aid early recognition of complications after trauma and assist timing and modality of treatment strategies. Our aim was to analyse and characterise temporal patterns of OD in intensive care unit-admitted trauma patients. METHODS: We used group-based trajectory modelling to identify temporal trajectories of OD after trauma. Modelling was based on the joint development of all six subdomains comprising the sequential organ failure assessment score measured daily during the first two weeks post trauma. Further, the time for trajectories to stabilise and transition to final group assignments were evaluated.
Afzali M, Nilsson M, Palombo M, Jones DK. SPHERIOUSLY? The challenges of estimating sphere radius non-invasively in the human brain from diffusion MRI. Neuroimage. 2021;237:118183. doi:10.1016/j.neuroimage.2021.118183
The Soma and Neurite Density Imaging (SANDI) three-compartment model was recently proposed to disentangle cylindrical and spherical geometries, attributed to neurite and soma compartments, respectively, in brain tissue. There are some recent advances in diffusion-weighted MRI signal encoding and analysis (including the use of multiple so-called ’b-tensor’ encodings and analysing the signal in the frequency-domain) that have not yet been applied in the context of SANDI. In this work, using: (i) ultra-strong gradients; (ii) a combination of linear, planar, and spherical b-tensor encodings; and (iii) analysing the signal in the frequency domain, three main challenges to robust estimation of sphere size were identified: First, the Rician noise floor in magnitude-reconstructed data biases estimates of sphere properties in a non-uniform fashion. It may cause overestimation or underestimation of the spherical compartment size and density. This can be partly ameliorated by accounting for the noise floor in the estimation routine. Second, even when using the strongest diffusion-encoding gradient strengths available for human MRI, there is an empirical lower bound on the spherical signal fraction and radius that can be detected and estimated robustly. For the experimental setup used here, the lower bound on the sphere signal fraction was approximately 10%. We employed two different ways of establishing the lower bound for spherical radius estimates in white matter. The first, examining power-law relationships between the DW-signal and diffusion weighting in empirical data, yielded a lower bound of 7μm, while the second, pure Monte Carlo simulations, yielded a lower limit of 3μm and in this low radii domain, there is little differentiation in signal attenuation. Third, if there is sensitivity to the transverse intra-cellular diffusivity in cylindrical structures, e.g., axons and cellular projections, then trying to disentangle two diffusion-time-dependencies using one experimental parameter (i.e., change in frequency-content of the encoding waveform) makes spherical radii estimates particularly challenging. We conclude that due to the aforementioned challenges spherical radii estimates may be biased when the corresponding sphere signal fraction is low, which must be considered.
Di Biase MA, Cetin-Karayumak S, Lyall AE, Zalesky A, Cho KIK, Zhang F, Kubicki M, Rathi Y, Lyons MG, Bouix S, et al. White Matter Changes in Psychosis Risk Relate to Development and Are Not Impacted by the Transition to Psychosis. Mol Psychiatry. 2021;26(11):6833–44. doi:10.1038/s41380-021-01128-8
Subtle alterations in white matter microstructure are observed in youth at clinical high risk (CHR) for psychosis. However, the timing of these changes and their relationships to the emergence of psychosis remain unclear. Here, we track the evolution of white matter abnormalities in a large, longitudinal cohort of CHR individuals comprising the North American Prodrome Longitudinal Study (NAPLS-3). Multi-shell diffusion magnetic resonance imaging data were collected across multiple timepoints (1-5 over 1 year) in 286 subjects (aged 12-32 years): 25 CHR individuals who transitioned to psychosis (CHR-P; 61 scans), 205 CHR subjects with unknown transition outcome after the 1-year follow-up period (CHR-U; 596 scans), and 56 healthy controls (195 scans). Linear mixed effects models were fitted to infer the impact of age and illness-onset on variation in the fractional anisotropy of cellular tissue (FAT) and the volume fraction of extracellular free water (FW). Baseline measures of white matter microstructure did not differentiate between HC, CHR-U and CHR-P individuals. However, age trajectories differed between the three groups in line with a developmental effect: CHR-P and CHR-U groups displayed higher FAT in adolescence, and 4% lower FAT by 30 years of age compared to controls. Furthermore, older CHR-P subjects (20+ years) displayed 4% higher FW in the forceps major (p < 0.05). Prospective analysis in CHR-P did not reveal a significant impact of illness onset on regional FAT or FW, suggesting that transition to psychosis is not marked by dramatic change in white matter microstructure. Instead, clinical high risk for psychosis-regardless of transition outcome-is characterized by subtle age-related white matter changes that occur in tandem with development.
e IO, Makris N, Densmore M, McKinnon MC, Lanius RA. Altered basal forebrain BOLD signal variability at rest in posttraumatic stress disorder: A potential candidate vulnerability mechanism for neurodegeneration in PTSD. Hum Brain Mapp. 2021;42(11):3561–75. doi:10.1002/hbm.25454
Individuals with posttraumatic stress disorder (PTSD) are at increased risk for the development of various forms of dementia. Nevertheless, the neuropathological link between PTSD and neurodegeneration remains unclear. Degeneration of the human basal forebrain constitutes a pathological hallmark of neurodegenerative diseases, such as Alzheimer’s and Parkinson’s disease. In this seed-based resting-state (rs-)fMRI study identifying as outcome measure the temporal BOLD signal fluctuation magnitude, a seed-to-voxel analyses assessed temporal correlations between the average BOLD signal within a bilateral whole basal forebrain region-of-interest and each whole-brain voxel among individuals with PTSD (n = 65), its dissociative subtype (PTSD+DS) (n = 38) and healthy controls (n = 46). We found that compared both with the PTSD and healthy controls groups, the PTSD+DS group exhibited increased BOLD signal variability within two nuclei of the seed region, specifically in its extended amygdaloid region: the nucleus accumbens and the sublenticular extended amygdala. This finding is provocative, because it mimics staging models of neurodegenerative diseases reporting allocation of neuropathology in early disease stages circumscribed to the basal forebrain. Here, underlying candidate etiopathogenetic mechanisms are neurovascular uncoupling, decreased connectivity in local- and large-scale neural networks, or disrupted mesolimbic dopaminergic circuitry, acting indirectly upon the basal forebrain cholinergic pathways. These abnormalities may underpin reward-related deficits representing a putative link between persistent traumatic memory in PTSD and anterograde memory deficits in neurodegeneration. Observed alterations of the basal forebrain in the dissociative subtype of PTSD point towards the urgent need for further exploration of this region as a potential candidate vulnerability mechanism for neurodegeneration in PTSD.
He J, Zhang F, Xie G, Yao S, Feng Y, Bastos DCA, Rathi Y, Makris N, Kikinis R, Golby AJ, et al. Comparison of Multiple Tractography Methods for Reconstruction of the Retinogeniculate Visual Pathway Using Diffusion MRI. Hum Brain Mapp. 2021;42(12):3887–3904. doi:10.1002/hbm.25472
The retinogeniculate visual pathway (RGVP) conveys visual information from the retina to the lateral geniculate nucleus. The RGVP has four subdivisions, including two decussating and two nondecussating pathways that cannot be identified on conventional structural magnetic resonance imaging (MRI). Diffusion MRI tractography has the potential to trace these subdivisions and is increasingly used to study the RGVP. However, it is not yet known which fiber tracking strategy is most suitable for RGVP reconstruction. In this study, four tractography methods are compared, including constrained spherical deconvolution (CSD) based probabilistic (iFOD1) and deterministic (SD-Stream) methods, and multi-fiber (UKF-2T) and single-fiber (UKF-1T) unscented Kalman filter (UKF) methods. Experiments use diffusion MRI data from 57 subjects in the Human Connectome Project. The RGVP is identified using regions of interest created by two clinical experts. Quantitative anatomical measurements and expert anatomical judgment are used to assess the advantages and limitations of the four tractography methods. Overall, we conclude that UKF-2T and iFOD1 produce the best RGVP reconstruction results. The iFOD1 method can better quantitatively estimate the percentage of decussating fibers, while the UKF-2T method produces reconstructed RGVPs that are judged to better correspond to the known anatomy and have the highest spatial overlap across subjects. Overall, we find that it is challenging for current tractography methods to both accurately track RGVP fibers that correspond to known anatomy and produce an approximately correct percentage of decussating fibers. We suggest that future algorithm development for RGVP tractography should take consideration of both of these two points.
Langbein BJ, Szczepankiewicz F, Westin C-F, Bay C, Maier SE, Kibel AS, Tempany CM, Fennessy FM. A Pilot Study of Multidimensional Diffusion MRI for Assessment of Tissue Heterogeneity in Prostate Cancer. Invest Radiol. 2021;56(12):845–53. doi:10.1097/RLI.0000000000000796
OBJECTIVES: The objectives of this exploratory study were to investigate the feasibility of multidimensional diffusion magnetic resonance imaging (MddMRI) in assessing diffusion heterogeneity at both a macroscopic and microscopic level in prostate cancer (PCa). MATERIALS AND METHODS: Informed consent was obtained from 46 subjects who underwent 3.0-T prostate multiparametric MRI, complemented with a prototype spin echo-based MddMRI sequence in this institutional review board-approved study. Prostate cancer tumors and comparative normal tissue from each patient were contoured on both apparent diffusion coefficient and MddMRI-derived mean diffusivity (MD) maps (from which microscopic diffusion heterogeneity [MKi] and microscopic diffusion anisotropy were derived) using 3D Slicer. The discriminative ability of MddMRI-derived parameters to differentiate PCa from normal tissue was determined using the Friedman test. To determine if tumor diffusion heterogeneity is similar on macroscopic and microscopic scales, the linear association between SD of MD and mean MKi was estimated using robust regression (bisquare weighting). Hypothesis testing was 2 tailed; P values less than 0.05 were considered statistically significant. RESULTS: All MddMRI-derived parameters could distinguish tumor from normal tissue in the fixed-effects analysis (P < 0.0001). Tumor MKi was higher (P < 0.05) compared with normal tissue (median, 0.40; interquartile range, 0.29-0.52 vs 0.20-0.18; 0.25), as was tumor microscopic diffusion anisotropy (0.55; 0.36-0.81 vs 0.20-0.15; 0.28). The MKi could not be predicted (no significant association) by SD of MD. There was a significant correlation between tumor volume and SD of MD (R2 = 0.50, slope = 0.008 μm2/ms per millimeter, P < 0.001) but not between tumor volume and MKi. CONCLUSIONS: This explorative study demonstrates that MddMRI provides novel information on MKi and microscopic anisotropy, which differ from measures at the macroscopic level. MddMRI has the potential to characterize tumor tissue heterogeneity at different spatial scales.
Herberthson M, Boito D, Haije TD, Feragen A, Westin C-F, Özarslan E. Q-space Trajectory Imaging with Positivity Constraints (QTI+). Neuroimage. 2021;238:118198. doi:10.1016/j.neuroimage.2021.118198
Q-space trajectory imaging (QTI) enables the estimation of useful scalar measures indicative of the local tissue structure. This is accomplished by employing generalized gradient waveforms for diffusion sensitization alongside a diffusion tensor distribution (DTD) model. The first two moments of the underlying DTD are made available by acquisitions at low diffusion sensitivity (b-values). Here, we show that three independent conditions have to be fulfilled by the mean and covariance tensors associated with distributions of symmetric positive semidefinite tensors. We introduce an estimation framework utilizing semi-definite programming (SDP) to guarantee that these conditions are met. Applying the framework on simulated signal profiles for diffusion tensors distributed according to non-central Wishart distributions demonstrates the improved noise resilience of QTI+ over the commonly employed estimation methods. Our findings on a human brain data set also reveal pronounced improvements, especially so for acquisition protocols featuring few number of volumes. Our method’s robustness to noise is expected to not only improve the accuracy of the estimates, but also enable a meaningful interpretation of contrast in the derived scalar maps. The technique’s performance on shorter acquisitions could make it feasible in routine clinical practice.
Fan X, Wang C, Han J, Ding X, Tang S, Ning L. Role of TRPV4-P2X7 Pathway in Neuropathic Pain in Rats with Chronic Compression of the Dorsal Root Ganglion. Neurochem Res. 2021;46(8):2143–53. doi:10.1007/s11064-021-03352-8
Transient receptor potential vanilloid 4 (TRPV4) is a Ca2+-permeable non-selective cation channel that is involved in the development of neuropathic pain. P2X7 receptor (P2X7) belongs to a class of ATP-gated nonselective cation channels that plays an important role in neuropathic pain. Nevertheless, little is known about the interaction between them for neuropathic pain. In this paper, we investigated role of TRPV4-P2X7 pathway in neuropathic pain. We evaluated the effect of TRPV4-P2X7 pathway on neuropathic pain in a chronic compression of the dorsal root ganglion (DRG) (hereafter termed CCD) model. We analyzed the effect of P2X7 on mechanical and thermal hyperalgesia mediated by TRPV4 in CCD. Furthermore, we assessed the effect of TRPV4 on the expression of P2X7 and the release of IL-1β and IL-6 in DRG after CCD. We found that intraperitoneal injection of TRPV4 agonist GSK-1016790A led to a significant increase of mechanical and thermal hyperalgesia in CCD, which was partially suppressed by P2X7 blockade with antagonist Brilliant Blue G (BBG). Then, we further noticed that GSK-1016790A injection increased the P2X7 expression of CCD, which was decreased by TRPV4 blockade with antagonist RN-1734 and HC-067047. Furthermore, we also discovered that the expressions of IL-1β and IL-6 were upregulated by GSK-1016790A injection but reduced by RN-1734 and HC-067047. Our results provide evidence that P2X7 contributes to development of neuropathic pain mediated by TRPV4 in the CCD model, which may be the basis for treatment of neuropathic pain relief.