In this paper, we propose an automated white matter connectivity analysis method for machine learning classification and characterization of white matter abnormality via identification of ...discriminative fiber tracts. The proposed method uses diffusion MRI tractography and a data-driven approach to find fiber clusters corresponding to subdivisions of the white matter anatomy. Features extracted from each fiber cluster describe its diffusion properties and are used for machine learning. The method is demonstrated by application to a pediatric neuroimaging dataset from 149 individuals, including 70 children with autism spectrum disorder (ASD) and 79 typically developing controls (TDC). A classification accuracy of 78.33% is achieved in this cross-validation study. We investigate the discriminative diffusion features based on a two-tensor fiber tracking model. We observe that the mean fractional anisotropy from the second tensor (associated with crossing fibers) is most affected in ASD. We also find that local along-tract (central cores and endpoint regions) differences between ASD and TDC are helpful in differentiating the two groups. These altered diffusion properties in ASD are associated with multiple robustly discriminative fiber clusters, which belong to several major white matter tracts including the corpus callosum, arcuate fasciculus, uncinate fasciculus and aslant tract; and the white matter structures related to the cerebellum, brain stem, and ventral diencephalon. These discriminative fiber clusters, a small part of the whole brain tractography, represent the white matter connections that could be most affected in ASD. Our results indicate the potential of a machine learning pipeline based on white matter fiber clustering.
Radiomics is an emerging field which extracts quantitative radiology data from medical images and explores their correlation with clinical outcomes in a non-invasive manner. This review aims to ...assess whether radiomics is a useful and reproducible method for clinical management of hepatocellular carcinoma (HCC) by reviewing the strengths and weaknesses of current radiomics literature pertaining specifically to HCC. From an initial set of 48 articles recovered through database searches, 23 articles were retained to be included in this review after full screening. Among these 23 studies, 7 used a radiomics approach in magnetic resonance imaging (MRI). Only two studies applied radiomics to positron emission tomography–computed tomography (PET–CT). In the remaining 14 articles, a radiomics analysis was performed on computed tomography (CT). Eight studies dealt with the relationship between biological signatures and imaging findings, and can be classified as radiogenomic studies. For each study included in our review, we computed a Radiomics Quality Score (RQS) as proposed by Lambin et al. We found that the RQS (mean ± standard deviation) was 8.35 ± 5.38 (out of a possible maximum value of 36). Although these scores are fairly low, and radiomics has not yet reached clinical utility in HCC, it is important to underscore the fact that these early studies pave the way for the radiomics field with a focus on HCC. Radiomics is still a very young field, and is far from being mature, but it remains a very promising technology for the future for developing adequate personalized treatment as a non-invasive approach, for complementing or replacing tumor biopsies, as well as for developing novel prognostic biomarkers in HCC patients.
The mammalian heart must function as an efficient pump while simultaneously conducting electrical signals to drive the contraction process. In the ventricles, electrical activation begins at the ...insertion points of the Purkinje network in the endocardium. How does the diffusion component of the subsequent excitation wave propagate from the endocardium in a healthy heart wall without creating directional biases? We show that this is a consequence of the particular geometric organization of myocytes in the heart wall. Using a generalized helicoid to model fiber orientation, we treat the myocardium as a curved space via Riemannian geometry, and then use stochastic calculus to model local signal diffusion. Our analysis shows that the helicoidal arrangement of myocytes minimizes the directional biases that could lead to aberrant propagation, thereby explaining how electrophysiological principles are consistent with local measurements of cardiac fiber geometry. We discuss our results in the context of the need to balance electrical and mechanical requirements for heart function.
Heart wall myofibers wind as helices around the ventricles, strengthening them in a manner analogous to the reinforcement of concrete cylindrical columns by spiral steel cables Richart FE, et al. ...(1929) Univ of Illinois, Eng Exp Stn Bull 190. A multitude of such fibers, arranged smoothly and regularly, contract and relax as an integrated functional unit as the heart beats. To orchestrate this motion, fiber tangling must be avoided and pumping should be efficient. Current models of myofiber orientation across the heart wall suggest groupings into sheets or bands, but the precise geometry of bundles of myofibers is unknown. Here we show that this arrangement takes the form of a special minimal surface, the generalized helicoid Blair DE, Vanstone JR (1978) Minimal Submanifolds and Geodesics 13–16, closing the gap between individual myofibers and their collective wall structure. The model holds across species, with a smooth variation in its three curvature parameters within the myocardial wall providing tight fits to diffusion magnetic resonance images from the rat, the dog, and the human. Mathematically it explains how myofibers are bundled in the heart wall while economizing fiber length and optimizing ventricular ejection volume as they contract. The generalized helicoid provides a unique foundation for analyzing the fibrous composite of the heart wall and should therefore find applications in heart tissue engineering and in the study of heart muscle diseases.
The temporal pole (TP) is considered one of the major paralimbic cortical regions, and is involved in a variety of functions such as sensory perception, emotion, semantic processing, and social ...cognition. Based on differences in cytoarchitecture, the TP can be further subdivided into smaller regions (dorsal, ventrolateral and ventromedial), each forming key nodes of distinct functional networks. However, the brain structural connectivity profile of TP subregions is not fully clarified. Using diffusion MRI data in a set of 31 healthy subjects, we aimed to elucidate the comprehensive structural connectivity of three cytoarchitectonically distinct TP subregions. Diffusion tensor imaging (DTI) analysis suggested that major association fiber pathways such as the inferior longitudinal, middle longitudinal, arcuate, and uncinate fasciculi provide structural connectivity to the TP. Further analysis suggested partially overlapping yet still distinct structural connectivity patterns across the TP subregions. Specifically, the dorsal subregion is strongly connected with wide areas in the parietal lobe, the ventrolateral subregion with areas including constituents of the default-semantic network, and the ventromedial subregion with limbic and paralimbic areas. Our results suggest the involvement of the TP in a set of extensive but distinct networks of cortical regions, consistent with its functional roles.
Diffusion MRI (dMRI) data acquired on different scanners varies significantly in its content throughout the brain even if the acquisition parameters are nearly identical. Thus, proper harmonization ...of such data sets is necessary to increase the sample size and thereby the statistical power of neuroimaging studies. In this paper, we present a novel approach to harmonize dMRI data (the raw signal, instead of dMRI derived measures such as fractional anisotropy) using rotation invariant spherical harmonic (RISH) features embedded within a multi-modal image registration framework. All dMRI data sets from all sites are registered to a common template and voxel-wise differences in RISH features between sites at a group level are used to harmonize the signal in a subject-specific manner. We validate our method on diffusion data acquired from seven different sites (two GE, three Philips, and two Siemens scanners) on a group of age-matched healthy subjects. We demonstrate the efficacy of our method by statistically comparing diffusion measures such as fractional anisotropy, mean diffusivity and generalized fractional anisotropy across these sites before and after data harmonization. Validation was also done on a group oftest subjects, which were not used to “learn” the harmonization parameters. We also show results using TBSS before and after harmonization for independent validation of the proposed methodology. Using synthetic data, we show that any abnormality in diffusion measures due to disease is preserved during the harmonization process. Our experimental results demonstrate that, for nearly identical acquisition protocol across sites, scanner-specific differences in the signal can be removed using the proposed method in a model independent manner.
Dual-energy computed tomography (DECT) is an advanced CT computed tomography scanning technique enabling material characterization not possible with conventional CT scans. It allows the ...reconstruction of energy decay curves at each 3D image voxel, representing varied image attenuation at different effective scanning energy levels. In this paper, we develop novel unsupervised learning techniques based on mixture models and functional data analysis models to the clustering of DECT images. We design functional mixture models that integrate spatial image context in mixture weights, with mixture component densities being constructed upon the DECT energy decay curves as functional observations. We develop dedicated expectation-maximization algorithms for the maximum likelihood estimation of the model parameters. To our knowledge, this is the first article to develop statistical functional data analysis and model-based clustering techniques to take advantage of the full spectral information provided by DECT. We evaluate the application of DECT to head and neck squamous cell carcinoma. Current image-based evaluation of these tumors in clinical practice is largely qualitative, based on a visual assessment of tumor anatomic extent and basic one- or two-dimensional tumor size measurements. We evaluate our methods on 91 head and neck cancer DECT scans and compare our unsupervised clustering results to tumor contours traced manually by radiologists, as well as to several baseline algorithms. Given the inter-rater variability even among experts at delineating head and neck tumors, and given the potential importance of tissue reactions surrounding the tumor itself, our proposed methodology has the potential to add value in downstream machine learning applications for clinical outcome prediction based on DECT data in head and neck cancer.
The corticospinal tract (CST) is one of the most well studied tracts in human neuroanatomy. Its clinical significance can be demonstrated in many notable traumatic conditions and diseases such as ...stroke, spinal cord injury (SCI) or amyotrophic lateral sclerosis (ALS). With the advent of diffusion MRI and tractography the computational representation of the human CST in a 3D model became available. However, the representation of the entire CST and, specifically, the hand motor area has remained elusive. In this paper we propose a novel method, using manually drawn ROIs based on robustly identifiable neuroanatomic structures to delineate the entire CST and isolate its hand motor representation as well as to estimate their variability and generate a database of their volume, length and biophysical parameters. Using 37 healthy human subjects we performed a qualitative and quantitative analysis of the CST and the hand-related motor fiber tracts (HMFTs). Finally, we have created variability heat maps from 37 subjects for both the aforementioned tracts, which could be utilized as a reference for future studies with clinical focus to explore neuropathology in both trauma and disease states.
The cingulum bundle (CB) connects gray matter structures of the limbic system and as such has been implicated in the etiology of schizophrenia. There is growing evidence to suggest that the CB is ...actually comprised of a conglomeration of discrete sub-connections. The present study aimed to use Diffusion Tensor tractography to subdivide the CB into its constituent sub-connections, and to investigate the structural integrity of these sub-connections in patients with schizophrenia and matched healthy controls. Diffusion Tensor Imaging scans were acquired from 24 patients diagnosed with chronic schizophrenia and 26 matched healthy controls. Deterministic tractography was used in conjunction with FreeSurfer-based regions-of-interest to subdivide the CB into 5 sub-connections (I1 to I5). The patients with schizophrenia exhibited subnormal levels of FA in two cingulum sub-connections, specifically the fibers connecting the rostral and caudal anterior cingulate gyrus (I1) and the fibers connecting the isthmus of the cingulate with the parahippocampal cortex (I4). Furthermore, while FA in the I1 sub-connection was correlated with the severity of patients' positive symptoms (specifically hallucinations and delusions), FA in the I4 sub-connection was correlated with the severity of patients' negative symptoms (specifically affective flattening and anhedonia/asociality). These results support the notion that the CB is a conglomeration of structurally interconnected yet functionally distinct sub-connections, of which only a subset are abnormal in patients with schizophrenia. Furthermore, while acknowledging the fact that the present study only investigated the CB, these results suggest that the positive and negative symptoms of schizophrenia may have distinct neurobiological underpinnings.