ConFiG: Contextual Fibre Growth to generate realistic axonal packing for diffusion MRI simulation.
Ontology highlight
ABSTRACT:
SUBMITTER: Callaghan R
PROVIDER: S-EPMC7903162 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BACKGROUND: Models of cellular molecular systems are built from components such as biochemical reactions (including interactions between ligands and membrane-bound proteins), conformational changes and active and passive transport. A discrete, stochastic description of the kinetics is often essential to capture the behavior of the system accurately. Where spatial effects play a prominent role the complex morphology of cells may have to be represented, along with aspects such as chemical localization and diffusion. This high level of detail makes efficiency a particularly important consideration for software that is designed to simulate such systems. RESULTS: We describe STEPS, a stochastic reaction-diffusion simulator developed with an emphasis on simulating biochemical signaling pathways accurately and efficiently. STEPS supports all the above-mentioned features, and well-validated support for SBML allows many existing biochemical models to be imported reliably. Complex boundaries can be represented accurately in externally generated 3D tetrahedral meshes imported by STEPS. The powerful Python interface facilitates model construction and simulation control. STEPS implements the composition and rejection method, a variation of the Gillespie SSA, supporting diffusion between tetrahedral elements within an efficient search and update engine. Additional support for well-mixed conditions and for deterministic model solution is implemented. Solver accuracy is confirmed with an original and extensive validation set consisting of isolated reaction, diffusion and reaction-diffusion systems. Accuracy imposes upper and lower limits on tetrahedron sizes, which are described in detail. By comparing to Smoldyn, we show how the voxel-based approach in STEPS is often faster than particle-based methods, with increasing advantage in larger systems, and by comparing to MesoRD we show the efficiency of the STEPS implementation. CONCLUSION: STEPS simulates models of cellular reaction-diffusion systems with complex boundaries with high accuracy and high performance in C/C++, controlled by a powerful and user-friendly Python interface. STEPS is free for use and is available at http://steps.sourceforge.net/
| S-EPMC3472240 | biostudies-literature
Project description:Diffusion MRI and tractography allow for investigation of the architectural configuration of white matter in vivo, offering new avenues for applications like presurgical planning. Despite the promising outlook, there are many pitfalls that complicate its use for (clinical) application. Amongst these are inaccuracies in the geometry of the diffusion profiles on which tractography is based, and poor alignment with neighboring profiles. Recently developed contextual processing techniques, including enhancement and well-posed geometric sharpening, have shown to result in sharper and better aligned diffusion profiles. However, the research that has been conducted up to now is mainly of theoretical nature, and so far these techniques have only been evaluated by visual inspection of the diffusion profiles. In this work, the method is evaluated in a clinically relevant application: the reconstruction of the optic radiation for epilepsy surgery. For this evaluation we have developed a framework in which we incorporate a novel scoring procedure for individual pathways. We demonstrate that, using enhancement and sharpening, the extraction of an anatomically plausible reconstruction of the optic radiation from a large amount of probabilistic pathways is greatly improved in three healthy controls, where currently used methods fail to do so. Furthermore, challenging reconstructions of the optic radiation in three epilepsy surgery candidates with extensive brain lesions demonstrate that it is beneficial to integrate these methods in surgical planning.
| S-EPMC4117467 | biostudies-literature
Project description:Diffusion MRI combined with biophysical modeling allows for the description of a white matter (WM) fiber bundle in terms of compartment specific white matter tract integrity (WMTI) metrics, which include intra-axonal diffusivity (Daxon), extra-axonal axial diffusivity (De||), extra-axonal radial diffusivity (De?), axonal water fraction (AWF), and tortuosity (?) of extra-axonal space. Here we derive these parameters from diffusion kurtosis imaging to examine their relationship to concentrations of global WM N-acetyl-aspartate (NAA), creatine (Cr), choline (Cho) and myo-Inositol (mI), as measured with proton MR spectroscopy ((1)H-MRS), in a cohort of 25 patients with mild traumatic brain injury (MTBI). We found statistically significant (p<0.05) positive correlations between NAA and Daxon, AWF, ?, and fractional anisotropy; negative correlations between NAA and De,? and the overall radial diffusivity (D?). These correlations were supported by similar findings in regional analysis of the genu and splenium of the corpus callosum. Furthermore, a positive correlation in global WM was noted between Daxon and Cr, as well as a positive correlation between De|| and Cho, and a positive trend between De|| and mI. The specific correlations between NAA, an endogenous probe of the neuronal intracellular space, and WMTI metrics related to the intra-axonal space, combined with the specific correlations of De|| with mI and Cho, both predominantly present extra-axonally, corroborate the overarching assumption of many advanced modeling approaches that diffusion imaging can disentangle between the intra- and extra-axonal compartments in WM fiber bundles. Our findings are also generally consistent with what is known about the pathophysiology of MTBI, which appears to involve both intra-axonal injury (as reflected by a positive trend between NAA and Daxon) as well as axonal shrinkage, demyelination, degeneration, and/or loss (as reflected by correlations between NAA and De?, AWF, and ?).
| S-EPMC4651014 | biostudies-literature
Project description:The cingulum is a prominent white matter tract that supports prefrontal, parietal, and temporal lobe interactions. Despite being composed of both short and long association fibres, many MRI-based reconstructions (tractography) of the cingulum depict an essentially uniform tract that almost encircles the corpus callosum. The present study tested the validity of dividing this tract into subdivisions corresponding to the 'parahippocampal', 'retrosplenial', and 'subgenual' portions of the cingulum. These three cingulum subdivisions occupied different medial-lateral locations, producing a topographic arrangement of cingulum fibres. Other comparisons based on these different reconstructions indicate that only a small proportion of the total white matter in the cingulum traverses the length of the tract. In addition, both the radial diffusivity and fractional anisotropy of the subgenual subdivision differed from that of the retrosplenial subdivision which, in turn, differed from that of the parahippocampal subdivision. The extent to which the radial diffusivity scores and the fractional anisotropy scores correlated between the various cingulum subdivisions proved variable, illustrating how one subdivision may not act as a proxy for other cingulum subdivisions. Attempts to relate the status of the cingulum, as measured by MRI-based fibre tracking, with cognitive or affective measures will, therefore, depend greatly on how and where the cingulum is reconstructed. The present study provides a new framework for subdividing the cingulum, based both on its known connectivity and MRI-based properties.
| S-EPMC3611599 | biostudies-literature
Project description:PurposeAnalytical phantoms have closed form Fourier transform expressions and are used to simulate MRI acquisitions. Existing three-dimensional (3D) analytical phantoms are unable to accurately model shapes of biomedical interest. The goal of this study was to demonstrate that polyhedral analytical phantoms have closed form Fourier transform expressions and can accurately represent 3D biomedical shapes.MethodsThe Fourier transform of a polyhedron was implemented and its accuracy in representing faceted and smooth surfaces was characterized. Realistic anthropomorphic polyhedral brain and torso phantoms were constructed and their use in simulated 3D and two-dimensional (2D) MRI acquisitions was described.ResultsUsing polyhedra, the Fourier transform of faceted shapes can be computed to within machine precision. Smooth surfaces can be approximated with increasing accuracy by increasing the number of facets in the polyhedron; the additional accumulated numerical imprecision of the Fourier transform of polyhedra with many faces remained small. Simulations of 3D and 2D brain and 2D torso cine acquisitions produced realistic reconstructions free of high frequency edge aliasing compared with equivalent voxelized/rasterized phantoms.ConclusionAnalytical polyhedral phantoms are easy to construct and can accurately simulate shapes of biomedical interest. Magn Reson Med 76:663-678, 2016. © 2015 Wiley Periodicals, Inc.
| S-EPMC4837112 | biostudies-literature
Project description:Analyzing the fluid behavior in complex porous media like gas diffusion layers (GDLs) in polymer electrolyte fuel cells (PEFCs) can be accurately done using the lattice Boltzmann method (LBM). This article shows the data obtained from a study in which diffusion parameters such as porosity, gas phase tortuosity and diffusibility are computed considering simulated porous media [1]. The data were computed when a water drop obstacle is placed inside the GDL domain and the size of the water-drop is varied. Additionally, figures showing the evolution of the flow velocity field are presented alongside graphics that presents the change in local and bulk porosity for each obstacle size. Finally, there is a detailed method explanation concerning the implementation of the lattice Boltzmann method and a general description of computational codes for the domain and obstacle generation as well as the boundary conditions simulation. Data and processes in this article can be exploited in new attempts to solve real case problems in complex mesoscale media.
| S-EPMC6838375 | biostudies-literature
Project description:PURPOSE:Quantitative evaluation of dynamic contrast enhanced MRI (DCE-MRI) allows for estimating perfusion, vessel permeability, and tissue volume fractions by fitting signal intensity curves to pharmacokinetic models. These compart mental models assume rapid equilibration of contrast agent within each voxel. However, there is increasing evidence that this assumption is violated for small molecular weight gadolinium chelates. To evaluate the error introduced by this invalid assumption, we simulated DCE-MRI experiments with volume fractions computed from entire histological tumor cross-sections obtained from murine studies. METHODS:A 2D finite element model of a diffusion-compensated Tofts-Kety model was developed to simulate dynamic T1 signal intensity data. Digitized histology slices were segmented into vascular (vp ), cellular and extravascular extracellular (ve ) volume fractions. Within this domain, Ktrans (the volume transfer constant) was assigned values from 0 to 0.5?min-1 . A representative signal enhancement curve was then calculated for each imaging voxel and the resulting simulated DCE-MRI data analyzed by the extended Tofts-Kety model. RESULTS:Results indicated parameterization errors of -19.1%?±?10.6% in Ktrans , -4.92%?±?3.86% in ve , and 79.5%?±?16.8% in vp for use of Gd-DTPA over 4 tumor domains. CONCLUSION:These results indicate a need for revising the standard model of DCE-MRI to incorporate a correction for slow diffusion of contrast agent. Magn Reson Med 80:330-340, 2018. © 2017 International Society for Magnetic Resonance in Medicine.
| S-EPMC5876107 | biostudies-literature
Project description:Cerebral white matter exhibits age-related degenerative changes during the course of normal aging, including decreases in axon density and alterations in axonal structure. Noninvasive approaches to measure these microstructural alterations throughout the lifespan would be invaluable for understanding the substrate and regional variability of age-related white matter degeneration. Recent advances in diffusion magnetic resonance imaging (MRI) have leveraged high gradient strengths to increase sensitivity toward axonal size and density in the living human brain. Here, we examined the relationship between age and indices of axon diameter and packing density using high-gradient strength diffusion MRI in 36 healthy adults (aged 22-72) in well-defined central white matter tracts in the brain. A recently validated method for inferring the effective axonal compartment size and packing density from diffusion MRI measurements acquired with 300 mT/m maximum gradient strength was applied to the in vivo human brain to obtain indices of axon diameter and density in the corpus callosum, its sub-regions, and adjacent anterior and posterior fibers in the forceps minor and forceps major. The relationships between the axonal metrics, corpus callosum area and regional gray matter volume were also explored. Results revealed a significant increase in axon diameter index with advancing age in the whole corpus callosum. Similar analyses in sub-regions of the corpus callosum showed that age-related alterations in axon diameter index and axon density were most pronounced in the genu of the corpus callosum and relatively absent in the splenium, in keeping with findings from previous histological studies. The significance of these correlations was mirrored in the forceps minor and forceps major, consistent with previously reported decreases in FA in the forceps minor but not in the forceps major with age. Alterations in the axonal imaging metrics paralleled decreases in corpus callosum area and regional gray matter volume with age. Among older adults, results from cognitive testing suggested an association between larger effective compartment size in the corpus callosum, particularly within the genu of the corpus callosum, and lower scores on the Montreal Cognitive Assessment, largely driven by deficits in short-term memory. The current study suggests that high-gradient diffusion MRI may be sensitive to the axonal substrate of age-related white matter degeneration reflected in traditional DTI metrics and provides further evidence for regionally selective alterations in white matter microstructure with advancing age.
| S-EPMC6467051 | biostudies-literature
Project description:BackgroundDeep Brain Stimulation (DBS) of the Globus pallidus internus (GPi) is gold standard treatment in medically refractory dystonia. Recent evidence indicates that stimulation effects are also due to axonal modulation and affection of a fibre network. For the GPi, the pallidothalamic tracts are known to be the major motor efferent pathways. The aim of this study is to explore the anatomic vicinity of these tracts and DBS electrodes in dystonia applying diffusion tractography.MethodsDiffusion MRI was acquired in ten patients presenting for DBS for dystonia. We applied both a conventionally used probabilistic tractography algorithm (FSL) as well as a probabilistic streamline tracking approach, based on constrained spherical deconvolution and particle filtering with anatomic priors, to the datasets. DBS electrodes were coregistered to the diffusion datasets.ResultsWe were able to delineate the pallidothalamic tracts in all patients. Using the streamline approach, we were able to distinguish between the two sub-components of the tracts, the ansa lenticularis and the fasciculus lenticularis. Clinically efficient DBS electrodes displayed a close anatomic vicinity pathway of the pallidothalamic tracts, and their course was consistent with previous tracer labelling studies. Although we present only anatomic data, we interpret these findings as evidence of the possible involvement of fibre tracts to the clinical effect in DBS. Electrophysiological intraoperative recordings would be needed to complement our findings. In the future, a clear and individual delineation of the pallidothalamic tracts could optimize the stereotactic process of optimal electrode localization. Hum Brain Mapp 38:1224-1232, 2017. © 2016 Wiley Periodicals, Inc.
| S-EPMC6867131 | biostudies-literature
Project description:BackgroundMicroarray technologies have become common tools in biological research. As a result, a need for effective computational methods for data analysis has emerged. Numerous different algorithms have been proposed for analyzing the data. However, an objective evaluation of the proposed algorithms is not possible due to the lack of biological ground truth information. To overcome this fundamental problem, the use of simulated microarray data for algorithm validation has been proposed.ResultsWe present a microarray simulation model which can be used to validate different kinds of data analysis algorithms. The proposed model is unique in the sense that it includes all the steps that affect the quality of real microarray data. These steps include the simulation of biological ground truth data, applying biological and measurement technology specific error models, and finally simulating the microarray slide manufacturing and hybridization. After all these steps are taken into account, the simulated data has realistic biological and statistical characteristics. The applicability of the proposed model is demonstrated by several examples.ConclusionThe proposed microarray simulation model is modular and can be used in different kinds of applications. It includes several error models that have been proposed earlier and it can be used with different types of input data. The model can be used to simulate both spotted two-channel and oligonucleotide based single-channel microarrays. All this makes the model a valuable tool for example in validation of data analysis algorithms.
| S-EPMC1574357 | biostudies-literature