Faculty Bibliography

BACKGROUND:Quantitative MRI can detect early changes in cartilage biochemical components, but its routine clinical implementation is challenging. PURPOSE/OBJECTIVE:along radial sections of the hip for accurate and reproducible multiparametric quantitative cartilage assessment in a clinically feasible scan time. STUDY TYPE/METHODS:Reproducibility, technical validation. SUBJECTS/PHANTOM/UNASSIGNED:A seven-compartment phantom and three healthy volunteers. FIELD STRENGTH/SEQUENCE/UNASSIGNED:at 3 T was developed. Automatic positioning and semiautomatic cartilage segmentation were implemented to improve consistency and simplify workflow. ASSESSMENT/RESULTS:Intra- and interscanner variability of our technique was assessed over multiple scans on three different MR scanners. STATISTICAL TESTS/UNASSIGNED:over six radial slices was calculated. Restricted maximum likelihood estimation of variance components was used to estimate intrasubject variances reflecting variation between results from the two scans using the same scanner (intrascanner variance) and variation among results from the three scanners (interscanner variance). RESULTS:. DATA CONCLUSION/UNASSIGNED:Our method, which includes slice positioning, model-based parameter estimation, and cartilage segmentation, is highly reproducible. It could enable employing quantitative hip cartilage evaluation for longitudinal and multicenter studies. LEVEL OF EVIDENCE/METHODS:1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018.

Autofluorescence is the emission of light by naturally occurring tissue components on the absorption of incident light. Autofluorescence within the eye is associated with several disorders, such as Age-related Macular Degeneration (AMD) which is a leading cause of central vision loss. Its pathogenesis is incompletely understood, but endogenous fluorophores in retinal tissue might play a role. Hyperspectral fluorescence microscopy of ex-vivo retinal tissue can be used to determine the fluorescence emission spectra of these fluorophores. Comparisons of spectra in healthy and diseased tissues can provide important insights into the pathogenesis of AMD. However, the spectrum from each pixel of the hyperspectral image is a superposition of spectra from multiple overlapping tissue components. As spectra cannot be negative, there is a need for a non-negative blind source separation model to isolate individual spectra. We propose a tensor formulation by leveraging multiple excitation wavelengths to excite the tissue sample. Arranging images from different excitation wavelengths as a tensor, a non-negative tensor decomposition can be performed to recover a provably unique low-rank model with factors representing emission and excitation spectra of these materials and corresponding abundance maps of autofluorescent substances in the tissue sample. We iteratively impute missing values common in fluorescence measurements using Expectation-Maximization and use L₂ regularization to reduce ill-posedness. Further, we present a framework for performing group hypothesis testing on hyperspectral images, finding significant differences in spectra between AMD and control groups in the peripheral macula. In the absence of ground truth, i.e. molecular identification of fluorophores, we provide a rigorous validation of chosen methods on both synthetic and real images where fluorescence spectra are known. These methodologies can be applied to the study of other pathologies presenting autofluorescence that can be captured by hyperspectral imaging.

The analysis of medical image time-series is becoming increasingly important as longitudinal imaging studies are maturing and large scale open imaging databases are becoming available. Image regression is widely used for several purposes: as a statistical representation for hypothesis testing, to bring clinical scores and images not acquired at the same time into temporal correspondence, or as a consistency filter to enforce temporal correlation. Geodesic image regression is the most prominent method, but the geodesic constraint limits the flexibility and therefore the application of the model, particularly when the observation time window is large or the anatomical changes are non-monotonic. In this paper, we propose to parameterize diffeomorphic flow by acceleration rather than velocity, as in the geodesic model. This results in a nonparametric image regression model which is completely flexible to capture complex change trajectories, while still constrained to be diffeomorphic and with a guarantee of temporal smoothness. We demonstrate the application of our model on synthetic 2D images as well as real 3D images of the cardiac cycle.

Spectral imaging is a ubiquitous tool in modern biochemistry. Despite acquiring dozens to thousands of spectral channels, existing technology cannot capture spectral images at the same spatial resolution as structural microscopy. Due to partial voluming and low light exposure, spectral images are often difficult to interpret and analyze. This highlights a need to upsample the low-resolution spectral image by using spatial information contained in the high-resolution image, thereby creating a fused representation with high specificity both spatially and spectrally. In this paper, we propose a framework for the fusion of co-registered structural and spectral microscopy images to create super-resolved representations of spectral images. As a first application, we super-resolve spectral images of retinal tissue imaged with confocal laser scanning microscopy, by using spatial information from structured illumination microscopy. Second, we super-resolve mass spectroscopic images of mouse brain tissue, by using spatial information from high-resolution histology images. We present a systematic validation of model assumptions crucial towards maintaining the original nature of spectra and the applicability of super-resolution. Goodness-of-fit for spectral predictions are evaluated through functional R² values, and the spatial quality of the super-resolved images are evaluated using normalized mutual information.

Spatiotemporal shape models capture the dynamics of shape change over time and are an essential tool for monitoring and measuring anatomical growth or degeneration. In this paper we evaluate non-parametric shape regression on the challenging problem of modeling early childhood sub-cortical development starting from birth. Due to the flexibility of the model, it can be challenging to choose parameters which lead to a good model fit yet does not over fit. We systematically test a variety of parameter settings to evaluate model fit as well as the sensitivity of the method to specific parameters, and we explore the impact of missing data on model estimation.

BACKGROUND:Restricted and repetitive behaviors (RRBs), detectable by 12 months in many infants in whom autism spectrum disorder (ASD) is later diagnosed, may represent some of the earliest behavioral markers of ASD. However, brain function underlying the emergence of these key behaviors remains unknown. METHODS:Behavioral and resting-state functional connectivity (fc) magnetic resonance imaging data were collected from 167 children at high and low familial risk for ASD at 12 and 24 months (n = 38 at both time points). Twenty infants met criteria for ASD at 24 months. We divided RRBs into four subcategories (restricted, stereotyped, ritualistic/sameness, self-injurious) and used a data-driven approach to identify functional brain networks associated with the development of each RRB subcategory. RESULTS:Higher scores for ritualistic/sameness behavior were associated with less positive fc between visual and control networks at 12 and 24 months. Ritualistic/sameness and stereotyped behaviors were associated with less positive fc between visual and default mode networks at 12 months. At 24 months, stereotyped and restricted behaviors were associated with more positive fc between default mode and control networks. Additionally, at 24 months, stereotyped behavior was associated with more positive fc between dorsal attention and subcortical networks, whereas restricted behavior was associated with more positive fc between default mode and dorsal attention networks. No significant network-level associations were observed for self-injurious behavior. CONCLUSIONS:These observations mark the earliest known description of functional brain systems underlying RRBs, reinforce the construct validity of RRB subcategories in infants, and implicate specific neural substrates for future interventions targeting RRBs.

ITK-SNAP is an interactive software tool for manual and semi-automatic segmentation of 3D medical images. This paper summarizes major new features added to ITK-SNAP over the last decade. The main focus of the paper is on new features that support semi-automatic segmentation of multi-modality imaging datasets, such as MRI scans acquired using different contrast mechanisms (e.g., T1, T2, FLAIR). The new functionality uses decision forest classifiers trained interactively by the user to transform multiple input image volumes into a foreground/background probability map; this map is then input as the data term to the active contour evolution algorithm, which yields regularized surface representations of the segmented objects of interest. The new functionality is evaluated in the context of high-grade and low-grade glioma segmentation by three expert neuroradiogists and a non-expert on a reference dataset from the MICCAI 2013 Multi-Modal Brain Tumor Segmentation Challenge (BRATS). The accuracy of semi-automatic segmentation is competitive with the top specialized brain tumor segmentation methods evaluated in the BRATS challenge, with most results obtained in ITK-SNAP being more accurate, relative to the BRATS reference manual segmentation, than the second-best performer in the BRATS challenge; and all results being more accurate than the fourth-best performer. Segmentation time is reduced over manual segmentation by 2.5 and 5 times, depending on the rater. Additional experiments in interactive placenta segmentation in 3D fetal ultrasound illustrate the generalizability of the new functionality to a different problem domain.

Understanding of early brain changes has the potential to investigate imaging biomarkers for pre-symptomatic diagnosis and thus opportunity for optimal therapeutic intervention, for example in early diagnosis of infants at risk to autism or altered development of infants to drug exposure. In this paper, we propose a framework to analyze longitudinal changes of structural connectivity in the early developing infant brain by exploring underlying network components of brain structural connectivity and its changes with age. Structural connectivity is a non-negative sparse network. Projective non-negative matrix factorization (PNMF) offers benefits in sparsity and learning fewer parameters for non-negative sparse data. The number of matrix subcomponents was estimated by automatic relevance determination PNMF (ARDPNMF) for brain connectivity networks for the given data. We apply linear mixed effect modeling on the resulting loadings from ARDPNMF to model longitudinal network component changes over time. The proposed framework was validated on a synthetic example generated by known linear mixed effects on loadings of the known number of bases with different levels of additive noises. Feasibility of the framework on real data has been demonstrated by analysis of structural connectivity networks of high angular resonance diffusion imaging (HARDI) data from an ongoing neuroimaging study of autism. A total of 139 image data sets from high-risk and low-risk subjects acquired at multiple time points have been processed. Results demonstrate the feasibility of the framework to analyze connectivity network properties as a function of age and the potential to eventually explore differences associated with risk status