Search

PurposeRecent advances in computational resources, including software libraries and hardware, have enabled the use of high-dimensional, multi-modal datasets to build Artificial Intelligence (AI) models and workflows for radiation therapy and image analysis. The purpose of Software Toolbox for RAdioTherapy and Imaging analysiS (STRATIS) is to provide cloud-enabled, easy-to-share software workflows to train and deploy AI models for transparency and multi-institutional collaboration.MethodSTRATIS leverages open source medical image informatics software for application-specific analysis. Jupyter notebooks for AI modeling workflows are provided with Python language as the base kernel. In addition to Python, workflows use software written in other languages, such as MATLAB, GNU-Octave, R, and C++, with the help of bridge libraries. The workflows can be run on a cloud platform, local workstation, or an institutional HPC cluster. Computational environments are provided in the form of publicly available docker images -and build scripts for local Anaconda environments. Utilities provided with STRATIS simplify bookkeeping of associations between imaging objects and allow chaining data processing operations defined via a setting file for AI models.ResultsWorkflows available on STRATIS can be broadly categorized into image segmentation, deformable image registration, and outcomes modeling for radiotherapy toxicity and tumor control using radiomics and dosimetry features. The STRATIS-forge GitHub organizationhttps://www.github.com/stratis-forgehosts build-scripts for Docker and Anaconda as well as Jupyter notebooks for analysis workflows. The software for building environments and workflow notebooks has open source-GNU-GPL copyright, and AI models retain the copyright chosen by their original developers.ConclusionSTRATIS enables researchers to deploy and share AI modeling workflows for radiotherapy and image analysis. STRATIS is publicly available on Terra.bio’s FireCloud platform with a pre-deployed computational environment and on GitHub organization for users pursuing local deployment.

BACKGROUND AND PURPOSE: To determine the ability of diffusion-weighted imaging (DWI) and dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) to predict long-term response of brain metastases prior to and within 72 hours of stereotactic radiosurgery (SRS). METHODS: In this prospective pilot study, multiple b-value DWI and T1-weighted DCE-MRI were performed in patients with brain metastases before and within 72 hours following SRS. Diffusion-weighted (DW) images were analyzed using the monoexponential and intravoxel incoherent motion (IVIM) models. DCE-MRI data were analyzed using the extended Tofts pharmacokinetic model. The parameters obtained with these methods were correlated with brain metastasis outcomes according to modified Response Assessment in Neuro-Oncology Brain Metastases (RANO-BM) criteria. RESULTS: We included 25 lesions from 16 patients; 16 patients underwent pre-SRS MRI and 12/16 patients underwent both pre- and early (within 72 hours) post-SRS MRI. The perfusion fraction (f) derived from IVIM early post-SRS was higher in lesions demonstrating progressive disease than in lesions demonstrating stable disease, partial response, or complete response (q = 0.041). Pre-SRS extracellular extravascular volume fraction, v(e), and volume transfer coefficient, k(trans), derived from DCE-MRI were higher in non-responders versus responders (q = 0.041). CONCLUSIONS: Quantitative DWI and DCE-MRI are feasible imaging methods in the pre- and early (within 72 hours) post-SRS evaluation of brain metastases. DWI- and DCE-MRI-derived parameters demonstrated physiologic changes (tumor cellularity and vascularity) and offer potentially useful biomarkers that can predict treatment response. This allows for initiation of alternate therapies within an effective time window that may help prevent disease progression.

We aimed to assess longitudinal changes in quantitative imaging metric values obtained from diffusion-weighted (DW-) and dynamic contrast-enhanced magnetic resonance imaging (DCE)-MRI at pre-treatment (TX[0]), immediately after the first fraction of stereotactic body radiotherapy (D1-TX[1]), and 6 weeks post-TX (Post-TX[2]) in patients with pancreatic ductal adenocarcinoma. Ten enrolled patients (n = 10) underwent DW- and DCE-MRI examinations on a 3.0 T scanner. The apparent diffusion coefficient, ADC (mm2/s), was derived from DW imaging data using a monoexponential model. The tissue relaxation rate, R1t, time-course data were fitted with a shutter-speed model, which provides estimates of the volume transfer constant, Ktrans (min−1), extravascular extracellular volume fraction, ve, and mean lifetime of intracellular water protons, τi (seconds). Wilcoxon rank-sum test compared the mean values, standard deviation, skewness, kurtosis, and relative percentage (r, %) changes (Δ) in ADC, Ktrans, ve, and τi values between the magnetic resonance examinations. rADCΔ2–0 values were significantly greater than rADCΔ1-0 values (P = .009). rKtransΔ2–0 values were significantly lower than rKtrans Δ1-0 values (P = .048). rveΔ2-1 and rveΔ2-0 values were significantly different (P = .016). rτiΔ2-1 values were significantly lower than rτiΔ2-0 values (P = .008). For group comparison, the pre-TX mean and kurtosis of ADC (P = .18 and P = .14), skewness and kurtosis of Ktrans values (P = .14 for both) showed a leaning toward significant difference between patients who experienced local control (n = 2) and failed early (n = 4). DW- and DCE-MRI-derived quantitative metrics could be useful biomarkers to evaluate longitudinal changes to stereotactic body radiotherapy in patients with pancreatic ductal adenocarcinoma.

We developed and tested the feasibility of computational fluid modeling (CFM) based on dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) for quantitative estimation of interstitial fluid pressure (IFP) and velocity (IFV) in patients with head and neck (HN) cancer with locoregional lymph node metastases. Twenty-two patients with HN cancer, with 38 lymph nodes, underwent pretreatment standard MRI, including DCE-MRI, on a 3-Tesla scanner. CFM simulation was performed with the finite element method in COMSOL Multiphysics software. The model consisted of a partial differential equation (PDE) module to generate 3D parametric IFP and IFV maps, using the Darcy equation and Ktrans values (min−1, estimated from the extended Tofts model) to reflect fluid influx into tissue from the capillary microvasculature. The Spearman correlation (ρ) was calculated between total tumor volumes and CFM estimates of mean tumor IFP and IFV. CFM-estimated tumor IFP and IFV mean ± standard deviation for the neck nodal metastases were 1.73 ± 0.39 (kPa) and 1.82 ± 0.9 × (10−7 m/s), respectively. High IFP estimates corresponds to very low IFV throughout the tumor core, but IFV rises rapidly near the tumor boundary where the drop in IFP is precipitous. A significant correlation was found between pretreatment total tumor volume and CFM estimates of mean tumor IFP (ρ = 0.50, P = 0.004). Future studies can validate these initial findings in larger patients with HN cancer cohorts using CFM of the tumor in concert with DCE characterization, which holds promise in radiation oncology and drug-therapy clinical trials.

BACKGROUND Diffusion-weighted imaging (DWI) is commonly used to detect prostate cancer, and a major clinical challenge is differentiating aggressive from indolent disease. PURPOSE To compare 14 site-specific parametric fitting implementations applied to the same dataset of whole-mount pathologically validated DWI to test the hypothesis that cancer differentiation varies with different fitting algorithms. STUDY TYPE Prospective. POPULATION Thirty-three patients prospectively imaged prior to prostatectomy. FIELD STRENGTH/SEQUENCE 3 T, field-of-view optimized and constrained undistorted single-shot DWI sequence. ASSESSMENT Datasets, including a noise-free digital reference object (DRO), were distributed to the 14 teams, where locally implemented DWI parameter maps were calculated, including mono-exponential apparent diffusion coefficient (MEADC), kurtosis (K), diffusion kurtosis (DK), bi-exponential diffusion (BID), pseudo-diffusion (BID*), and perfusion fraction (F). The resulting parametric maps were centrally analyzed, where differentiation of benign from cancerous tissue was compared between DWI parameters and the fitting algorithms with a receiver operating characteristic area under the curve (ROC AUC). STATISTICAL TEST Levene's test, P

The present study aims to monitor longitudinal changes in simulated tumor interstitial fluid pressure (IFP) and velocity (IFV) values using dynamic contrast-enhanced (DCE)-MRI-based computational fluid modeling (CFM) in pancreatic ductal adenocarcinoma (PDAC) patients. Nine PDAC patients underwent MRI, including DCE-MRI, on a 3-Tesla MRI scanner at pre-treatment (TX (0)), after the first fraction of stereotactic body radiotherapy (SBRT, (D1-TX)), and six weeks post-TX (D2-TX). The partial differential equation of IFP formulated from the continuity equation, incorporating the Starling Principle of fluid exchange, Darcy velocity, and volume transfer constant (Ktrans), was solved in COMSOL Multiphysics software to generate IFP and IFV maps. Tumor volume (Vt), Ktrans, IFP, and IFV values were compared (Wilcoxon and Spearman) between the time- points. D2-TX Ktrans values were significantly different from pre-TX and D1-TX (p &lt, 0.05). The D1-TX and pre-TX mean IFV values exhibited a borderline significant difference (p = 0.08). The IFP values varying &lt, 3.0% between the three time-points were not significantly different (p &gt, 0.05). Vt and IFP values were strongly positively correlated at pre-TX (r = 0.90, p = 0.005), while IFV exhibited a strong negative correlation at D1-TX (r = −0.74, p = 0.045). Vt, Ktrans, IFP, and IFV hold promise as imaging biomarkers of early response to therapy in PDAC.

Current state-of-the-art models for estimating the pharmacokinetic parameters do not account for intervoxel movement of the contrast agent (CA). We introduce an optimal mass transport (OMT) formulation that naturally handles intervoxel CA movement and distinguishes between advective and diffusive flows. Ten patients with head and neck squamous cell carcinoma (HNSCC) were enrolled in the study between June 2014 and October 2015 and under-went DCE MRI imaging prior to beginning treatment. The CA tissue concentration information was taken as the input in the data-driven OMT model. The OMT approach was tested on HNSCC DCE data that provides quantitative information for forward flux (ΦF) and backward flux (ΦB). OMT-derived ΦF was compared with the volume transfer constant for CA, Ktrans, derived from the Extended Tofts Model (ETM). The OMT-derived flows showed a consistent jump in the CA diffusive behavior across the images in accordance with the known CA dynamics. The mean forward flux was 0.0082 ± 0.0091 (min-1) whereas the mean advective component was 0.0052±0.0086 (min-1) in the HNSCC patients. The diffusive percentages in forward and backward flux ranged from 8.67–18.76% and 12.76–30.36%, respectively. The OMT model accounts for intervoxel CA movement and results show that the forward flux (ΦF) is comparable with the ETM-derived Ktrans. This is a novel data-driven study based on optimal mass transport principles applied to patient DCE imaging to analyze CA flow in HNSCC.

Quantitative kurtosis phantoms are sought by multicenter clinical trials to establish accuracy and precision of quantitative imaging biomarkers on the basis of diffusion kurtosis imaging (DKI) parameters. We designed and evaluated precision, reproducibility, and long-term stability of a novel isotropic (i) DKI phantom fabricated using four families of chemicals based on vesicular and lamellar mesophases of liquid crystal materials. The constructed iDKI phantoms included negative control monoexponential diffusion materials to independently characterize noise and model-induced bias in quantitative kurtosis parameters. Ten test–retest DKI studies were performed on four scanners at three imaging centers over a six-month period. The tested prototype phantoms exhibited physiologically relevant apparent diffusion, Dapp, and kurtosis, Kapp, parameters ranging between 0.4 and 1.1 (×10−3 mm2/s) and 0.8 and 1.7 (unitless), respectively. Measured kurtosis phantom Kapp exceeded maximum fit model bias (0.1) detected for negative control (zero kurtosis) materials. The material-specific parameter precision [95% CI for Dapp: 0.013–0.022(×10−3 mm2/s) and for Kapp: 0.009–0.076] derived from the test–retest analysis was sufficient to characterize thermal and temporal stability of the prototype DKI phantom through correlation analysis of inter-scan variability. The present study confirms a promising chemical design for stable quantitative DKI phantom based on vesicular mesophase of liquid crystal materials. Improvements to phantom preparation and temperature monitoring procedures have potential to enhance precision and reproducibility for future multicenter iDKI phantom studies.

The aim of this study was to establish the repeatability measures of quantitative Gaussian and non-Gaussian diffusion metrics using diffusion-weighted imaging (DWI) data from phantoms and patients with head-and-neck and papillary thyroid cancers. The Quantitative Imaging Biomarker Alliance (QIBA) DWI phantom and a novel isotropic diffusion kurtosis imaging phantom were scanned at 3 different sites, on 1.5T and 3T magnetic resonance imaging systems, using standardized multiple b-value DWI acquisition protocol. In the clinical component of this study, a total of 60 multiple b-value DWI data sets were analyzed for test–retest, obtained from 14 patients (9 head-and-neck squamous cell carcinoma and 5 papillary thyroid cancers). Repeatability of quantitative DWI measurements was assessed by within-subject coefficient of variation (wCV%) and Bland–Altman analysis. In isotropic diffusion kurtosis imaging phantom vial with 2% ceteryl alcohol and behentrimonium chloride solution, the mean apparent diffusion (Dapp × 10−3 mm2/s) and kurtosis (Kapp, unitless) coefficient values were 1.02 and 1.68 respectively, capturing in vivo tumor cellularity and tissue microstructure. For the same vial, Dapp and Kapp mean wCVs (%) were ≤1.41% and ≤0.43% for 1.5T and 3T across 3 sites. For pretreatment head-and-neck squamous cell carcinoma, apparent diffusion coefficient, D, D*, K, and f mean wCVs (%) were 2.38%, 3.55%, 3.88%, 8.0%, and 9.92%, respectively, wCVs exhibited a higher trend for papillary thyroid cancers. Knowledge of technical precision and bias of quantitative imaging metrics enables investigators to properly design and power clinical trials and better discern between measurement variability versus biological change.

PurposeThis work presents a framework for deployment of deep learning image segmentation models for medical images across different operating systems and programming languages.MethodsComputational Environment for Radiological Research (CERR) platform was extended for deploying deep learning-based segmentation models to leverage CERR’s existing functionality for radiological data import, transformation, management, and visualization. The framework is compatible with MATLAB as well as GNU Octave and Python for license-free use. Pre and post processing configurations including parameters for pre-processing images, population of channels, and post-processing segmentations was standardized using JSON format. CPU and GPU implementations of pre-trained deep learning segmentation models were packaged using Singularity containers for use in Linux and Conda environment archives for Windows, macOS and Linux operating systems. The framework accepts images in various formats including DICOM and CERR’s planC and outputs segmentation in various formats including DICOM RTSTRUCT and planC objects. The ability to access the results readily in planC format enables visualization as well as radiomics and dosimetric analysis. The framework can be readily deployed in clinical software such as MIM via their extensions.ResultsThe open-source, GPL copyrighted framework developed in this work has been successfully used to deploy Deep Learning based segmentation models for five in-house developed and published models. These models span various treatment sites (H&N, Lung and Prostate) and modalities (CT, MR). Documentation for their usage and demo workflow is provided at https://github.com/cerr/CERR/wiki/Auto-Segmentation-models. The framework has also been used in clinical workflow for segmenting images for treatment planning and for segmenting publicly available large datasets for outcomes studies.ConclusionsThis work presented a comprehensive, open-source framework for deploying deep learning-based medical image segmentation models. The framework was used to translate the developed models to clinic as well as reproducible and consistent image segmentation across institutions, facilitating multi-institutional outcomes modeling studies.

Purpose Contouring inconsistencies are known but understudied in clinical radiation therapy trials. We applied auto-contouring to the Radiation Therapy Oncology Group (RTOG) 0617 dose escalation trial data. We hypothesized that the trial heart doses were higher than reported due to inconsistent and insufficient heart segmentation. We tested our hypothesis by comparing doses between deep-learning (DL) segmented hearts and trial hearts. Methods and Materials The RTOG 0617 data were downloaded from The Cancer Imaging Archive; the 442 patients with trial hearts and dose distributions were included. All hearts were resegmented using our DL pipeline and quality assured to meet the requirements for clinical implementation. Dose (V5%, V30%, and mean heart dose) was compared between the 2 sets of hearts (Wilcoxon signed-rank test). Each dose metric was associated with overall survival (Cox proportional hazards). Lastly, 18 volume similarity metrics were assessed for the hearts and correlated with |DoseDL – DoseRTOG0617| (linear regression; significance: P ≤ .0028; corrected for 18 tests). Results Dose metrics were significantly higher for DL hearts compared with trial hearts (eg, mean heart dose: 15 Gy vs 12 Gy; P = 5.8E-16). All 3 DL heart dose metrics were stronger overall survival predictors than those of the trial hearts (median, P = 2.8E-5 vs 2.0E-4). Thirteen similarity metrics explained |DoseDL – DoseRTOG0617|; the axial distance between the 2 centers of mass was the strongest predictor (CENTAxial; median, R2 = 0.47; P = 6.1E-62). CENTAxial agreed with the qualitatively identified inconsistencies in the superior direction. The trial’s qualitative heart contouring score was not correlated with |DoseDL – DoseRTOG0617| (median, R2 = 0.01; P = .02) or with any of the similarity metrics (median, Rs = 0.13 [range, –0.22 to 0.31]). Conclusions Using a coherent heart definition, as enabled through our open-source DL algorithm, the trial heart doses in RTOG 0617 were found to be significantly higher than previously reported, which may have led to an even more rapid mortality accumulation. Auto-segmentation is likely to reduce contouring and dose inconsistencies and increase the quality of clinical RT trials.

Objective. Delineating swallowing and chewing structures aids in radiotherapy (RT) treatment planning to limit dysphagia, trismus, and speech dysfunction. We aim to develop an accurate and efficient method to automate this process. Approach. CT scans of 242 head and neck (H&N) cancer patients acquired from 2004 to 2009 at our institution were used to develop auto-segmentation models for the masseters, medial pterygoids, larynx, and pharyngeal constrictor muscle using DeepLabV3+. A cascaded framework was used, wherein models were trained sequentially to spatially constrain each structure group based on prior segmentations. Additionally, an ensemble of models, combining contextual information from axial, coronal, and sagittal views was used to improve segmentation accuracy. Prospective evaluation was conducted by measuring the amount of manual editing required in 91 H&N CT scans acquired February-May 2021. Main results. Medians and inter-quartile ranges of Dice similarity coefficients (DSC) computed on the retrospective testing set (N = 24) were 0.87 (0.85–0.89) for the masseters, 0.80 (0.79–0.81) for the medial pterygoids, 0.81 (0.79–0.84) for the larynx, and 0.69 (0.67–0.71) for the constrictor. Auto-segmentations, when compared to two sets of manual segmentations in 10 randomly selected scans, showed better agreement (DSC) with each observer than inter-observer DSC. Prospective analysis showed most manual modifications needed for clinical use were minor, suggesting auto-contouring could increase clinical efficiency. Trained segmentation models are available for research use upon request via https://github.com/cerr/CERR/wiki/Auto-Segmentation-models. Significance. We developed deep learning-based auto-segmentation models for swallowing and chewing structures in CT and demonstrated its potential for use in treatment planning to limit complications post-RT. To the best of our knowledge, this is the only prospectively-validated deep learning-based model for segmenting chewing and swallowing structures in CT. Segmentation models have been made open-source to facilitate reproducibility and multi-institutional research.

An open-source library of implementations for deep-learning based image segmentation and radiotherapy outcomes models is presented in this work. As oncology treatment planning becomes increasingly driven by automation, such a library of model implementations is crucial to (i) validate existing models on datasets collected at different institutions, (ii) automate segmentation, (iii) create ensembles for improving performance and (iv) incorporate validated models in the clinical workflow. The library was developed with the Computational Environment for Radiological Research (CERR) software platform. Centralizing model implementations in CERR builds upon its rich set of radiotherapy and radiomics tools and caters to the world-wide user base. CERR provides well-validated feature extraction pipelines for radiotherapy dosimetry and radiomics with fine control over the calculation settings, allowing users to select the appropriate parameters used in model derivation. Models for automatic image segmentation are distributed via containers, with seamless i/o to and from CERR. Container technology allow segmentation models to be deployed with a variety of scientific computing architectures. The library includes implementations of popular radiotherapy models outlined in the Quantitative Analysis of Normal Tissue Effects in the Clinic effort and recently published literature. Radiomics models include features from the Image Biomarker Standardization Initiative and application-specific features found to be relevant across multiple sites and image modalities. Deep learning-based image segmentation models include state-of-the-art networks such as DeepLab and other problem-specific architectures. The source code is distributed at https://www.github.com/cerr/CERR under the GNU-GPL copyleft with additional restrictions on clinical and commercial use and provision to dual license in future.

Patients with asymptomatic kidney stones have a high rate of progression to becoming symptomatic kidney stones when followed for several years. Small kidney stones are often found incidentally on imaging when evaluating patients for kidney donation, and there is a concern that after nephrectomy, the donor may become symptomatic and incur damage to the remaining kidney. We reviewed kidney donors at our institution with asymptomatic stones and surveyed them several years after donation to see if the stones became clinically active.

Current hypotheses stipulate core symptoms of schizophrenia (SZ) result from the brain's incapacity to integrate neural processes. Converging diffusion magnetic resonance imaging and graph theory studies provide evidence of macrostructural alterations in SZ. However, age-related topological changes within and between white matter (WM) networks and its relationship to gene expression with disease progression remain incompletely understood. This cross-sectional study uses network modeling to investigate changes in WM network organization with disease progression in chronic SZ as well its relationship with gene expression in healthy brains. First, we replicate prior findings demonstrating altered global WM network topology in SZ. Novel results show significantly altered age-related network degradation patterns in patients compared with controls. Specifically, controls show stereotyped, linear global network decline with age. In contrast, patients show nonlinear network decline with age. Further analysis reveals lack of significant topological decline in younger adult patients, which is subsequently followed by stereotyped linear decline in older adult patients. Node-specific analyses show significant topological differences in frontal and limbic regions of younger adult patients compared with age-matched controls, which become less pronounced with age in older adult patients compared with age-matched controls. Lastly, we show several gene expression profiles, including DISC1, are associated with age-related changes in WM disconnectivity. Together, these findings provide novel WM topological and genetic evidence supporting neurodevelopmental models of SZ, suggesting that network remodeling continues throughout the third decade of life before stabilizing.

Background Platelet-activating factor acetylhydrolase 1B1 (LIS1), a critical mediator of neuronal migration in developing brain, is expressed throughout life. However, relatively little is known about LIS1 function in the mature brain. We previously demonstrated that LIS1 involvement in the formation and turnover of synaptic protrusions and synapses of young brain after neuronal migration is complete. Here we examine the requirement for LIS1 to maintain hippocampal circuit function in adulthood. Methods Effects of conditional Lis1 inactivation in excitatory pyramidal neurons, starting in juvenile mouse brain, were probed using high-resolution approaches combining mouse genetics, designer receptor exclusively activated by designer drug technology to specifically manipulate CA1 pyramidal neuron excitatory activity, electrophysiology, hippocampus-selective behavioral testing, and magnetic resonance imaging tractography to examine the connectivity of LIS1-deficient neurons. Results We found progressive excitatory and inhibitory postsynaptic dysfunction as soon as 10 days after conditional inactivation of Lis1 targeting CA1 pyramidal neurons. Surprisingly, by postnatal day 60 it also caused CA1 histological disorganization, with a selective decline in parvalbumin-expressing interneurons and further reduction in inhibitory neurotransmission. Accompanying these changes were behavioral and cognitive deficits that could be rescued by either designer receptor exclusively activated by designer drug–directed specific increases in CA1 excitatory transmission or pharmacological enhancement of gamma-aminobutyric acid transmission. Lagging behind electrophysiological changes was a progressive, selective decline in neural connectivity, affecting hippocampal efferent pathways documented by magnetic resonance imaging tractography. Conclusions LIS1 supports synaptic function and plasticity of mature CA1 neurons. Postjuvenile loss of LIS1 disrupts the structure and cellular composition of the hippocampus, its connectivity with other brain regions, and cognition dependent on hippocampal circuits.

Unraveling the relationship between molecular signatures in the brain and their functional, architectonic, and anatomic correlates is an important neuroscientific goal. It is still not well understood whether the diversity demonstrated by histological studies in the human brain is reflected in the spatial patterning of whole brain transcriptional profiles. Using genome-wide maps of transcriptional distribution of the human brain by the Allen Brain Institute, we test the hypothesis that gene expression profiles are specific to anatomically described brain regions. In this work, we demonstrate that this is indeed the case by showing that gene similarity clusters appear to respect conventional basal-cortical and caudal-rostral gradients. To fully investigate the causes of this observed spatial clustering, we test a connectionist hypothesis that states that the spatial patterning of gene expression in the brain is simply reflective of the fiber tract connectivity between brain regions. We find that although gene expression and structural connectivity are not determined by each other, they do influence each other with a high statistical significance. This implies that spatial diversity of gene expressions is a result of mainly location-specific features but is influenced by neuronal connectivity, such that like cellular species preferentially connects with like cells.

The brain's myelin content can be mapped by T2-relaxometry, which resolves multiple differentially relaxing T2 pools from multi-echo MRI. Unfortunately, the conventional fitting procedure is a hard and numerically ill-posed problem. Consequently, the T2 distributions and myelin maps become very sensitive to noise and are frequently difficult to interpret diagnostically. Although regularization can improve stability, it is generally not adequate, particularly at relatively low signal to noise ratio (SNR) of around 100-200. The purpose of this study was to obtain a fitting algorithm which is able to overcome these difficulties and generate usable myelin maps from noisy acquisitions in a realistic scan time. To this end, we restrict the T2 distribution to only 3 distinct resolvable tissue compartments, modeled as Gaussians: myelin water, intra/extra-cellular water and a slow relaxing cerebrospinal fluid compartment. We also impose spatial smoothness expectation that volume fractions and T2 relaxation times of tissue compartments change smoothly within coherent brain regions. The method greatly improves robustness to noise, reduces spatial variations, improves definition of white matter fibers, and enhances detection of demyelinating lesions. Due to efficient design, the additional spatial aspect does not cause an increase in processing time. The proposed method was applied to fast spiral acquisitions on which conventional fitting gives uninterpretable results. While these fast acquisitions suffer from noise and inhomogeneity artifacts, our preliminary results indicate the potential of spatially constrained 3-pool T2 relaxometry.