Manifold Learning in MR spectroscopy using nonlinear dimensionality reduction and unsupervised clustering

Purpose To investigate whether nonlinear dimensionality reduction improves unsupervised classification of 1H MRS brain tumor data compared with a linear method. Methods In vivo single-voxel 1H magnetic resonance spectroscopy (55 patients) and 1H magnetic resonance spectroscopy imaging (MRSI) (29 patients) data were acquired from histopathologically diagnosed gliomas. Data reduction using Laplacian eigenmaps (LE) or independent component analysis (ICA) was followed by k-means clustering or agglomerative hierarchical clustering (AHC) for unsupervised learning to assess tumor grade and for tissue type segmentation of MRSI data. Results An accuracy of 93% in classification of glioma grade II and grade IV, with 100% accuracy in distinguishing tumor and normal spectra, was obtained by LE with unsupervised clustering, but not with the combination of k-means and ICA. With 1H MRSI data, LE provided a more linear distribution of data for cluster analysis and better cluster stability than ICA. LE combined with k-means or AHC provided 91% accuracy for classifying tumor grade and 100% accuracy for identifying normal tissue voxels. Color-coded visualization of normal brain, tumor core, and infiltration regions was achieved with LE combined with AHC. Conclusion Purpose To investigate whether nonlinear dimensionality reduction improves unsupervised classification of 1H MRS brain tumor data compared with a linear method. Methods In vivo single-voxel 1H magnetic resonance spectroscopy (55 patients) and 1H magnetic resonance spectroscopy imaging (MRSI) (29 patients) data were acquired from histopathologically diagnosed gliomas. Data reduction using Laplacian eigenmaps (LE) or independent component analysis (ICA) was followed by k-means clustering or agglomerative hierarchical clustering (AHC) for unsupervised learning to assess tumor grade and for tissue type segmentation of MRSI data. Results An accuracy of 93% in classification of glioma grade II and grade IV, with 100% accuracy in distinguishing tumor and normal spectra, was obtained by LE with unsupervised clustering, but not with the combination of k-means and ICA. With 1H MRSI data, LE provided a more linear distribution of data for cluster analysis and better cluster stability than ICA. LE combined with k-means or AHC provided 91% accuracy for classifying tumor grade and 100% accuracy for identifying normal tissue voxels. Color-coded visualization of normal brain, tumor core, and infiltration regions was achieved with LE combined with AHC. Conclusion The LE method is promising for unsupervised clustering to separate brain and tumor tissue with automated color-coding for visualization of 1H MRSI data after cluster analysis

Structural network efficiency is associated with cognitive impairment in small-vessel disease.

To characterize brain network connectivity impairment in cerebral small-vessel disease (SVD) and its relationship with MRI disease markers and cognitive impairment.METHODS: A cross-sectional design applied graph-based efficiency analysis to deterministic diffusion tensor tractography data from 115 patients with lacunar infarction and leukoaraiosis and 50 healthy individuals. Structural connectivity was estimated between 90 cortical and subcortical brain regions and efficiency measures of resulting graphs were analyzed. Networks were compared between SVD and control groups, and associations between efficiency measures, conventional MRI disease markers, and cognitive function were tested.RESULTS: Brain diffusion tensor tractography network connectivity was significantly reduced in SVD: networks were less dense, connection weights were lower, and measures of network efficiency were significantly disrupted. The degree of brain network disruption was associated with MRI measures of disease severity and cognitive function. In multiple regression models controlling for confounding variables, associations with cognition were stronger for network measures than other MRI measures including conventional diffusion tensor imaging measures. A total mediation effect was observed for the association between fractional anisotropy and mean diffusivity measures and executive function and processing speed.CONCLUSIONS: Brain network connectivity in SVD is disturbed, this disturbance is related to disease severity, and within a mediation framework fully or partly explains previously observed associations between MRI measures and SVD-related cognitive dysfunction. These cross-sectional results highlight the importance of network disruption in SVD and provide support for network measures as a disease marker in treatment studies

Disruption of white matter connectivity in chronic obstructive pulmonary disease.

BACKGROUND: Mild cognitive impairment is a common systemic manifestation of chronic obstructive pulmonary disease (COPD). However, its pathophysiological origins are not understood. Since, cognitive function relies on efficient communication between distributed cortical and subcortical regions, we investigated whether people with COPD have disruption in white matter connectivity. METHODS: Structural networks were constructed for 30 COPD patients (aged 54-84 years, 57% male, FEV1 52.5% pred.) and 23 controls (aged 51-81 years, 48% Male). Networks comprised 90 grey matter regions (nodes) interconnected by white mater fibre tracts traced using deterministic tractography (edges). Edges were weighted by the number of streamlines adjusted for a) streamline length and b) end-node volume. White matter connectivity was quantified using global and nodal graph metrics which characterised the networks connection density, connection strength, segregation, integration, nodal influence and small-worldness. Between-group differences in white matter connectivity and within-group associations with cognitive function and disease severity were tested. RESULTS: COPD patients' brain networks had significantly lower global connection strength (p = 0.03) and connection density (p = 0.04). There was a trend towards COPD patients having a reduction in nodal connection density and connection strength across the majority of network nodes but this only reached significance for connection density in the right superior temporal gyrus (p = 0.02) and did not survive correction for end-node volume. There were no other significant global or nodal network differences or within-group associations with disease severity or cognitive function. CONCLUSION: COPD brain networks show evidence of damage compared to controls with a reduced number and strength of connections. This loss of connectivity was not sufficient to disrupt the overall efficiency of network organisation, suggesting that it has redundant capacity that makes it resilient to damage, which may explain why cognitive dysfunction is not severe. This might also explain why no direct relationships could be found with cognitive measures. Smoking and hypertension are known to have deleterious effects on the brain. These confounding effects could not be excluded

White matter lesions characterise brain involvement in moderate to severe chronic obstructive pulmonary disease, but cerebral atrophy does not.

BACKGROUND: Brain pathology is relatively unexplored in chronic obstructive pulmonary disease (COPD). This study is a comprehensive investigation of grey matter (GM) and white matter (WM) changes and how these relate to disease severity and cognitive function. METHODS: T1-weighted and fluid-attenuated inversion recovery images were acquired for 31 stable COPD patients (FEV1 52.1% pred., PaO2 10.1 kPa) and 24 age, gender-matched controls. T1-weighted images were segmented into GM, WM and cerebrospinal fluid (CSF) tissue classes using a semi-automated procedure optimised for use with this cohort. This procedure allows, cohort-specific anatomical features to be captured, white matter lesions (WMLs) to be identified and includes a tissue repair step to correct for misclassification caused by WMLs. Tissue volumes and cortical thickness were calculated from the resulting segmentations. Additionally, a fully-automated pipeline was used to calculate localised cortical surface and gyrification. WM and GM tissue volumes, the tissue volume ratio (indicator of atrophy), average cortical thickness, and the number, size, and volume of white matter lesions (WMLs) were analysed across the whole-brain and regionally - for each anatomical lobe and the deep-GM. The hippocampus was investigated as a region-of-interest. Localised (voxel-wise and vertex-wise) variations in cortical gyrification, GM density and cortical thickness, were also investigated. Statistical models controlling for age and gender were used to test for between-group differences and within-group correlations. Robust statistical approaches ensured the family-wise error rate was controlled in regional and local analyses. RESULTS: There were no significant differences in global, regional, or local measures of GM between patients and controls, however, patients had an increased volume (p = 0.02) and size (p = 0.04) of WMLs. In patients, greater normalised hippocampal volume positively correlated with exacerbation frequency (p = 0.04), and greater WML volume was associated with worse episodic memory (p = 0.05). A negative relationship between WML and FEV1 % pred. approached significance (p = 0.06). CONCLUSIONS: There was no evidence of cerebral atrophy within this cohort of stable COPD patients, with moderate airflow obstruction. However, there were indications of WM damage consistent with an ischaemic pathology. It cannot be concluded whether this represents a specific COPD, or smoking-related, effect

Defining thalamic nuclei and topographic connectivity gradients in vivo.

The thalamus consists of multiple nuclei that have been previously defined by their chemoarchitectual and cytoarchitectual properties ex vivo. These form discrete, functionally specialized, territories with topographically arranged graduated patterns of connectivity. However, previous in vivo thalamic parcellation with MRI has been hindered by substantial inter-individual variability or discrepancies between MRI derived segmentations and histological sections. Here, we use the Euclidean distance to characterize probabilistic tractography distributions derived from diffusion MRI. We generate 12 feature maps by performing voxel-wise parameterization of the distance histograms (6 feature maps) and the distribution of three-dimensional distance transition gradients generated by applying a Sobel kernel to the distance metrics. We use these 12 feature maps to delineate individual thalamic nuclei, then extract the tractography profiles for each and calculate the voxel-wise tractography gradients. Within each thalamic nucleus, the tractography gradients were topographically arranged as distinct non-overlapping cortical networks with transitory overlapping mid-zones. This work significantly advances quantitative segmentation of the thalamus in vivo using 3T MRI. At an individual subject level, the thalamic segmentations consistently achieve a close relationship with a priori histological atlas information, and resolve in vivo topographic gradients within each thalamic nucleus for the first time. Additionally, these techniques allow individual thalamic nuclei to be closely aligned across large populations and generate measures of inter-individual variability that can be used to study both basic function and pathological processes in vivo

Quasi-diffusion magnetic resonance imaging (QDI): A fast, high b-value diffusion imaging technique.

To enable application of non-Gaussian diffusion magnetic resonance imaging (dMRI) techniques in large-scale clinical trials and facilitate translation to clinical practice there is a requirement for fast, high contrast, techniques that are sensitive to changes in tissue structure which provide diagnostic signatures at the early stages of disease. Here we describe a new way to compress the acquisition of multi-shell b-value diffusion data, Quasi-Diffusion MRI (QDI), which provides a probe of subvoxel tissue complexity using short acquisition times (1-4 min). We also describe a coherent framework for multi-directional diffusion gradient acquisition and data processing that allows computation of rotationally invariant quasi-diffusion tensor imaging (QDTI) maps. QDI is a quantitative technique that is based on a special case of the Continuous Time Random Walk model of diffusion dynamics and assumes the presence of non-Gaussian diffusion properties within tissue microstructure. QDI parameterises the diffusion signal attenuation according to the rate of decay (i.e. diffusion coefficient, D in mm2 s-1) and the shape of the power law tail (i.e. the fractional exponent, α). QDI provides analogous tissue contrast to Diffusional Kurtosis Imaging (DKI) by calculation of normalised entropy of the parameterised diffusion signal decay curve, Hn, but does so without the limitations of a maximum b-value. We show that QDI generates images with superior tissue contrast to conventional diffusion imaging within clinically acceptable acquisition times of between 84 and 228 s. We show that QDI provides clinically meaningful images in cerebral small vessel disease and brain tumour case studies. Our initial findings suggest that QDI may be added to routine conventional dMRI acquisitions allowing simple application in clinical trials and translation to the clinical arena

Brain tumor classification using the diffusion tensor image segmentation (D-SEG) technique.

BACKGROUND: There is an increasing demand for noninvasive brain tumor biomarkers to guide surgery and subsequent oncotherapy. We present a novel whole-brain diffusion tensor imaging (DTI) segmentation (D-SEG) to delineate tumor volumes of interest (VOIs) for subsequent classification of tumor type. D-SEG uses isotropic (p) and anisotropic (q) components of the diffusion tensor to segment regions with similar diffusion characteristics. METHODS: DTI scans were acquired from 95 patients with low- and high-grade glioma, metastases, and meningioma and from 29 healthy subjects. D-SEG uses k-means clustering of the 2D (p,q) space to generate segments with different isotropic and anisotropic diffusion characteristics. RESULTS: Our results are visualized using a novel RGB color scheme incorporating p, q and T2-weighted information within each segment. The volumetric contribution of each segment to gray matter, white matter, and cerebrospinal fluid spaces was used to generate healthy tissue D-SEG spectra. Tumor VOIs were extracted using a semiautomated flood-filling technique and D-SEG spectra were computed within the VOI. Classification of tumor type using D-SEG spectra was performed using support vector machines. D-SEG was computationally fast and stable and delineated regions of healthy tissue from tumor and edema. D-SEG spectra were consistent for each tumor type, with constituent diffusion characteristics potentially reflecting regional differences in tissue microstructure. Support vector machines classified tumor type with an overall accuracy of 94.7%, providing better classification than previously reported. CONCLUSIONS: D-SEG presents a user-friendly, semiautomated biomarker that may provide a valuable adjunct in noninvasive brain tumor diagnosis and treatment planning

Optimization of quasi-diffusion magnetic resonance imaging for quantitative accuracy and time-efficient acquisition.

Purpose Quasi-diffusion MRI (QDI) is a novel quantitative technique based on the continuous time random walk model of diffusion dynamics. QDI provides estimates of the diffusion coefficient D1,2, in mm2 s−1 and a fractional exponent, α, defining the non-Gaussianity of the diffusion signal decay. Here, the b-value selection for rapid clinical acquisition of QDI tensor imaging (QDTI) data is optimized. Methods Clinically appropriate QDTI acquisitions were optimized in healthy volunteers with respect to a multi-b-value reference (MbR) dataset comprising 29 diffusion-sensitized images arrayed between b = 0 and 5000 s mm−2. The effects of varying maximum b-value (bmax), number of b-value shells, and the effects of Rician noise were investigated. Results QDTI measures showed bmax dependence, most significantly for α in white matter, which monotonically decreased with higher bmax leading to improved tissue contrast. Optimized 2 b-value shell acquisitions showed small systematic differences in QDTI measures relative to MbR values, with overestimation of D1,2 and underestimation of α in white matter, and overestimation of D1,2 and α anisotropies in gray and white matter. Additional shells improved the accuracy, precision, and reliability of QDTI estimates with 3 and 4 shells at bmax = 5000 s mm−2, and 4 b-value shells at bmax = 3960 s mm−2, providing minimal bias in D1,2 and α compared to the MbR. Conclusion A highly detailed optimization of non-Gaussian dMRI for in vivo brain imaging was performed. QDI provided robust parameterization of non-Gaussian diffusion signal decay in clinically feasible imaging times with high reliability, accuracy, and precision of QDTI measures

The effect of phosphodiesterase-5 inhibitors on cerebral blood flow in humans: A systematic review.

Agents that augment cerebral blood flow (CBF) could be potential treatments for vascular cognitive impairment. Phosphodiesterase-5 inhibitors are vasodilating drugs established in the treatment of erectile dysfunction (ED) and pulmonary hypertension. We reviewed published data on the effects of phosphodiesterase-5 inhibitors on CBF in adult humans. A systematic review according to PRISMA guidelines was performed. Embase, Medline and Cochrane Library Trials databases were searched. Sixteen studies with 353 participants in total were retrieved. Studies included healthy volunteers and patients with migraine, ED, type 2 diabetes, stroke, pulmonary hypertension, Becker muscular dystrophy and subarachnoid haemorrhage. Most studies used middle cerebral artery flow velocity to estimate CBF. Few studies employed direct measurements of tissue perfusion. Resting CBF velocity was unaffected by phosphodiesterase-5 inhibitors, but cerebrovascular regulation was improved in ED, pulmonary hypertension, diabetes, Becker's and a group of healthy volunteers. This evidence suggests that phosphodiesterase-5 inhibitors improve responsiveness of the cerebral vasculature, particularly in disease states associated with an impaired endothelial dilatory response. This supports the potential therapeutic use of phosphodiesterase-5 inhibitors in vascular cognitive impairment where CBF is reduced. Further studies with better resolution of deep CBF are warranted. The review is registered on the PROSPERO database (registration number CRD42016029668)