Axon mimicking hydrophilic hollow polycaprolactone microfibres for diffusion magnetic resonance imaging

Pore size estimation in axon-mimicking microfibers with diffusion-relaxation MRI

PURPOSE

This study aims to evaluate two distinct approaches for fiber radius estimation using diffusion-relaxation MRI data acquired in biomimetic microfiber phantoms that mimic hollow axons. The methods considered are the spherical mean power-law approach and a T2-based pore size estimation technique.

THEORY AND METHODS

A general diffusion-relaxation theoretical model for the spherical mean signal from water molecules within a distribution of cylinders with varying radii was introduced, encompassing the evaluated models as particular cases. Additionally, a new numerical approach was presented for estimating effective radii (i.e., MRI-visible mean radii) from the ground truth radii distributions, not reliant on previous theoretical approximations and adaptable to various acquisition sequences. The ground truth radii were obtained from scanning electron microscope images.

RESULTS

Both methods show a linear relationship between effective radii estimated from MRI data and ground-truth radii distributions, although some discrepancies were observed. The spherical mean power-law method overestimated fiber radii. Conversely, the T2-based method exhibited higher sensitivity to smaller fiber radii, but faced limitations in accurately estimating the radius in one particular phantom, possibly because of material-specific relaxation changes.

CONCLUSION

The study demonstrates the feasibility of both techniques to predict pore sizes of hollow microfibers. The T2-based technique, unlike the spherical mean power-law method, does not demand ultra-high diffusion gradients, but requires calibration with known radius distributions. This research contributes to the ongoing development and evaluation of neuroimaging techniques for fiber radius estimation, highlights the advantages and limitations of both methods, and provides datasets for reproducible research.

CACTUS: a computational framework for generating realistic white matter microstructure substrates

Monte-Carlo diffusion simulations are a powerful tool for validating tissue microstructure models by generating synthetic diffusion-weighted magnetic resonance images (DW-MRI) in controlled environments. This is fundamental for understanding the link between micrometre-scale tissue properties and DW-MRI signals measured at the millimetre-scale, optimizing acquisition protocols to target microstructure properties of interest, and exploring the robustness and accuracy of estimation methods. However, accurate simulations require substrates that reflect the main microstructural features of the studied tissue. To address this challenge, we introduce a novel computational workflow, CACTUS (Computational Axonal Configurator for Tailored and Ultradense Substrates), for generating synthetic white matter substrates. Our approach allows constructing substrates with higher packing density than existing methods, up to 95% intra-axonal volume fraction, and larger voxel sizes of up to 500μm3 with rich fibre complexity. CACTUS generates bundles with angular dispersion, bundle crossings, and variations along the fibres of their inner and outer radii and g-ratio. We achieve this by introducing a novel global cost function and a fibre radial growth approach that allows substrates to match predefined targeted characteristics and mirror those reported in histological studies. CACTUS improves the development of complex synthetic substrates, paving the way for future applications in microstructure imaging.

Direct jet coaxial electrospinning of axon-mimicking fibers for diffusion tensor imaging

Hollow polymer microfibers with variable microstructural and hydrophilic properties were proposed as building elements to create axon-mimicking phantoms for validation of diffusion tensor imaging (DTI). The axon-mimicking microfibers were fabricated in a mm-thick 3D anisotropic fiber strip, by direct jet coaxial electrospinning of PCL/polysiloxane-based surfactant (PSi) mixture as shell and polyethylene oxide (PEO) as core. Hydrophilic PCL-PSi fiber strips were first obtained by carefully selecting appropriate solvents for the core and appropriate fiber collector rotating and transverse speeds. The porous cross-section and anisotropic orientation of axon-mimicking fibers were then quantitatively evaluated using two ImageJ plugins-nearest distance (ND) and directionality based on their scanning electron microscopy (SEM) images. Third, axon-mimicking phantom was constructed from PCL-PSi fiber strips with variable porous-section and fiber orientation and tested on a 3T clinical MR scanner. The relationship between DTI measurements (mean diffusivity [MD] and fractional anisotropy [FA]) of phantom samples and their pore size and fiber orientation was investigated. Two key microstructural parameters of axon-mimicking phantoms including normalized pore distance and dispersion of fiber orientation could well interpret the variations in DTI measurements. Two PCL-PSi phantom samples made from different regions of the same fiber strips were found to have similar MD and FA values, indicating that the direct jet coaxial electrospun fiber strips had consistent microstructure. More importantly, the MD and FA values of the developed axon-mimicking phantoms were mostly in the biologically relevant range.

Cardiac q-space trajectory imaging by motion-compensated tensor-valued diffusion encoding in human heart in vivo

PURPOSE

Tensor-valued diffusion encoding can probe more specific features of tissue microstructure than what is available by conventional diffusion weighting. In this work, we investigate the technical feasibility of tensor-valued diffusion encoding at high b-values with q-space trajectory imaging (QTI) analysis, in the human heart in vivo.

METHODS

Ten healthy volunteers were scanned on a 3T scanner. We designed time-optimal gradient waveforms for tensor-valued diffusion encoding (linear and planar) with second-order motion compensation. Data were analyzed with QTI. Normal values and repeatability were investigated for the mean diffusivity (MD), fractional anisotropy (FA), microscopic FA (μFA), isotropic, anisotropic and total mean kurtosis (MKi, MKa, and MKt), and orientation coherence (Cc ). A phantom, consisting of two fiber blocks at adjustable angles, was used to evaluate sensitivity of parameters to orientation dispersion and diffusion time.

RESULTS

QTI data in the left ventricular myocardium were MD = 1.62 ± 0.07 μm2 /ms, FA = 0.31 ± 0.03, μFA = 0.43 ± 0.07, MKa = 0.20 ± 0.07, MKi = 0.13 ± 0.03, MKt = 0.33 ± 0.09, and Cc  = 0.56 ± 0.22 (mean ± SD across subjects). Phantom experiments showed that FA depends on orientation dispersion, whereas μFA was insensitive to this effect.

CONCLUSION

We demonstrated the first tensor-valued diffusion encoding and QTI analysis in the heart in vivo, along with first measurements of myocardial μFA, MKi, MKa, and Cc . The methodology is technically feasible and provides promising novel biomarkers for myocardial tissue characterization.

Physical and digital phantoms for validating tractography and assessing artifacts

Fiber tractography is widely used to non-invasively map white-matter bundles in vivo using diffusion-weighted magnetic resonance imaging (dMRI). As it is the case for all scientific methods, proper validation is a key prerequisite for the successful application of fiber tractography, be it in the area of basic neuroscience or in a clinical setting. It is well-known that the indirect estimation of the fiber tracts from the local diffusion signal is highly ambiguous and extremely challenging. Furthermore, the validation of fiber tractography methods is hampered by the lack of a real ground truth, which is caused by the extremely complex brain microstructure that is not directly observable non-invasively and that is the basis of the huge network of long-range fiber connections in the brain that are the actual target of fiber tractography methods. As a substitute for in vivo data with a real ground truth that could be used for validation, a widely and successfully employed approach is the use of synthetic phantoms. In this work, we are providing an overview of the state-of-the-art in the area of physical and digital phantoms, answering the following guiding questions: "What are dMRI phantoms and what are they good for?", "What would the ideal phantom for validation fiber tractography look like?" and "What phantoms, phantom datasets and tools used for their creation are available to the research community?". We will further discuss the limitations and opportunities that come with the use of dMRI phantoms, and what future direction this field of research might take.

Validating pore size estimates in a complex microfiber environment on a human MRI system

PURPOSE

Recent advances in diffusion-weighted MRI provide "restricted diffusion signal fraction" and restricting pore size estimates. Materials based on co-electrospun oriented hollow cylinders have been introduced to provide validation for such methods. This study extends this work, exploring accuracy and repeatability using an extended acquisition on a 300 mT/m gradient human MRI scanner, in substrates closely mimicking tissue, that is, non-circular cross-sections, intra-voxel fiber crossing, intra-voxel distributions of pore-sizes, and smaller pore-sizes overall.

METHODS

In a single-blind experiment, diffusion-weighted data were collected from a biomimetic phantom on a 3T Connectom system using multiple gradient directions/diffusion times. Repeated scans established short-term and long-term repeatability. The total scan time (54 min) matched similar protocols used in human studies. The number of distinct fiber populations was estimated using spherical deconvolution, and median pore size estimated through the combination of CHARMED and AxCaliber3D framework. Diffusion-based estimates were compared with measurements derived from scanning electron microscopy.

RESULTS

The phantom contained substrates with different orientations, fiber configurations, and pore size distributions. Irrespective of one or two populations within the voxel, the pore-size estimates (~5 μm) and orientation-estimates showed excellent agreement with the median values of pore-size derived from scanning electron microscope and phantom configuration. Measurement repeatability depended on substrate complexity, with lower values seen in samples containing crossing-fibers. Sample-level repeatability was found to be good.

CONCLUSION

While no phantom mimics tissue completely, this study takes a step closer to validating diffusion microstructure measurements for use in vivo by demonstrating the ability to quantify microgeometry in relatively complex configurations.

Comparative analysis of signal models for microscopic fractional anisotropy estimation using q-space trajectory encoding

Microscopic diffusion anisotropy imaging using diffusion-weighted MRI and multidimensional diffusion encoding is a promising method for quantifying clinically and scientifically relevant microstructural properties of neural tissue. Several methods for estimating microscopic fractional anisotropy (µFA), a normalized measure of microscopic diffusion anisotropy, have been introduced but the differences between the methods have received little attention thus far. In this study, the accuracy and precision of µFA estimation using q-space trajectory encoding and different signal models were assessed using imaging experiments and simulations. Three healthy volunteers and a microfibre phantom were imaged with five non-zero b-values and gradient waveforms encoding linear and spherical b-tensors. Since the ground-truth µFA was unknown in the imaging experiments, Monte Carlo random walk simulations were performed using axon-mimicking fibres for which the ground truth was known. Furthermore, parameter bias due to time-dependent diffusion was quantified by repeating the simulations with tuned waveforms, which have similar power spectra, and with triple diffusion encoding, which, unlike q-space trajectory encoding, is not based on the assumption of time-independent diffusion. The truncated cumulant expansion of the powder-averaged signal, gamma-distributed diffusivities assumption, and q-space trajectory imaging, a generalization of the truncated cumulant expansion to individual signals, were used to estimate µFA. The gamma-distributed diffusivities assumption consistently resulted in greater µFA values than the second order cumulant expansion, 0.1 greater when averaged over the whole brain. In the simulations, the generalized cumulant expansion provided the most accurate estimates. Importantly, although time-dependent diffusion caused significant overestimation of µFA using all the studied methods, the simulations suggest that the resulting bias in µFA is less than 0.1 in human white matter.

Coaxial electrospun biomimetic copolymer fibres for application in diffusion magnetic resonance imaging

Objective. The use of diffusion magnetic resonance imaging (dMRI) opens the door to characterizing brain microstructure because water diffusion is anisotropic in axonal fibres in brain white matter and is sensitive to tissue microstructural changes. As dMRI becomes more sophisticated and microstructurally informative, it has become increasingly important to use a reference object (usually called an imaging phantom) for validation of dMRI. This study aims to develop axon-mimicking physical phantoms from biocopolymers and assess their feasibility for validating dMRI measurements.Approach. We employed a simple and one-step method-coaxial electrospinning-to prepare axon-mimicking hollow microfibres from polycaprolactone-b-polyethylene glycol (PCL-b-PEG) and poly(D, L-lactide-co-glycolic) acid (PLGA), and used them as building elements to create axon-mimicking phantoms. Electrospinning was firstly conducted using two types of PCL-b-PEG and two types of PLGA with different molecular weights in various solvents, with different polymer concentrations, for determining their spinnability. Polymer/solvent concentration combinations with good fibre spinnability were used as the shell material in the following co-electrospinning process in which the polyethylene oxide polymer was used as the core material. Following the microstructural characterization of both electrospun and co-electrospun fibres using optical and electron microscopy, two prototype phantoms were constructed from co-electrospun anisotropic hollow microfibres after inserting them into water-filled test tubes.Main results. Hollow microfibres that mimic the axon microstructure were successfully prepared from the appropriate core and shell material combinations. dMRI measurements of two phantoms on a 7 tesla (T) pre-clinical scanner revealed that diffusivity and anisotropy measurements are in the range of brain white matter.Significance. This feasibility study showed that co-electrospun PCL-b-PEG and PLGA microfibre-based axon-mimicking phantoms could be used in the validation of dMRI methods which seek to characterize white matter microstructure.

Stability and reproducibility of co-electrospun brain-mimicking phantoms for quality assurance of diffusion MRI sequences

Grey and white matter mimicking phantoms are important for assessing variations in diffusion MR measures at a single time point and over an extended period of time. This work investigates the stability of brain-mimicking microfibre phantoms and reproducibility of their MR derived diffusion parameters. The microfibres were produced by co-electrospinning and characterized by scanning electron microscopy (SEM). Grey matter and white matter phantoms were constructed from random and aligned microfibres, respectively. MR data were acquired from these phantoms over a period of 33 months. SEM images revealed that only small changes in fibre microstructure occurred over 30 months. The coefficient of variation in MR measurements across all time-points was between 1.6% and 3.4% for MD across all phantoms and FA in white matter phantoms. This was within the limits expected for intra-scanner variability, thereby confirming phantom stability over 33 months. These specialised diffusion phantoms may be used in a clinical environment for intra and inter-site quality assurance purposes, and for validation of quantitative diffusion biomarkers.