Continuous dynamic mapping of the corticospinal tract during surgery of motor eloquent brain tumors: evaluation of a new method

OBJECT The authors developed a new mapping technique to overcome the temporal and spatial limitations of classic subcortical mapping of the corticospinal tract (CST). The feasibility and safety of continuous (0.4-2 Hz) and dynamic (at the site of and synchronized with tissue resection) subcortical motor mapping was evaluated. METHODS The authors prospectively studied 69 patients who underwent tumor surgery adjacent to the CST (< 1 cm using diffusion tensor imaging and fiber tracking) with simultaneous subcortical monopolar motor mapping (short train, interstimulus interval 4 msec, pulse duration 500 μsec) and a new acoustic motor evoked potential alarm. Continuous (temporal coverage) and dynamic (spatial coverage) mapping was technically realized by integrating the mapping probe at the tip of a new suction device, with the concept that this device will be in contact with the tissue where the resection is performed. Motor function was assessed 1 day after surgery, at discharge, and at 3 months. RESULTS All procedures were technically successful. There was a 1:1 correlation of motor thresholds for stimulation sites simultaneously mapped with the new suction mapping device and the classic fingerstick probe (24 patients, 74 stimulation points; r(2) = 0.98, p < 0.001). The lowest individual motor thresholds were as follows: > 20 mA, 7 patients; 11-20 mA, 13 patients; 6-10 mA, 8 patients; 4-5 mA, 17 patients; and 1-3 mA, 24 patients. At 3 months, 2 patients (3%) had a persistent postoperative motor deficit, both of which were caused by a vascular injury. No patient had a permanent motor deficit caused by a mechanical injury of the CST. CONCLUSIONS Continuous dynamic mapping was found to be a feasible and ergonomic technique for localizing the exact site of the CST and distance to the motor fibers. The acoustic feedback and the ability to stimulate the tissue continuously and exactly at the site of tissue removal improves the accuracy of mapping, especially at low (< 5 mA) stimulation intensities. This new technique may increase the safety of motor eloquent tumor surgery.

Separation of small metabolites and lipids in spectra from biopsies by diffusion-weighted HR-MAS NMR: a feasibility study.

High Resolution Magic Angle Spinning (HR-MAS) NMR allows metabolic characterization of biopsies. HR-MAS spectra from tissues of most organs show strong lipid contributions that are overlapping metabolite regions, which hamper metabolite estimation. Metabolite quantification and analysis would benefit from a separation of lipids and small metabolites. Generally, a relaxation filter is used to reduce lipid contributions. However, the strong relaxation filter required to eliminate most of the lipids also reduces the signals for small metabolites. The aim of our study was therefore to investigate different diffusion editing techniques in order to employ diffusion differences for separating lipid and small metabolite contributions in the spectra from different organs for unbiased metabonomic analysis. Thus, 1D and 2D diffusion measurements were performed, and pure lipid spectra that were obtained at strong diffusion weighting (DW) were subtracted from those obtained at low DW, which include both small metabolites and lipids. This subtraction yielded almost lipid free small metabolite spectra from muscle tissue. Further improved separation was obtained by combining a 1D diffusion sequence with a T2-filter, with the subtraction method eliminating residual lipids from the spectra. Similar results obtained for biopsies of different organs suggest that this method is applicable in various tissue types. The elimination of lipids from HR-MAS spectra and the resulting less biased assessment of small metabolites have potential to remove ambiguities in the interpretation of metabonomic results. This is demonstrated in a reproducibility study on biopsies from human muscle.

Saddlepoint Approximations to the Probability of Ruin in Finite Time for the Compound Poisson Risk Process Perturbed by Diffusion

A large deviations type approximation to the probability of ruin within a finite time for the compound Poisson risk process perturbed by diffusion is derived. This approximation is based on the saddlepoint method and generalizes the approximation for the non-perturbed risk process by Barndorff-Nielsen and Schmidli (Scand Actuar J 1995(2):169–186, 1995). An importance sampling approximation to this probability of ruin is also provided. Numerical illustrations assess the accuracy of the saddlepoint approximation using importance sampling as a benchmark. The relative deviations between saddlepoint approximation and importance sampling are very small, even for extremely small probabilities of ruin. The saddlepoint approximation is however substantially faster to compute.

Simultaneous multi-slice readout-segmented echo planar imaging for accelerated diffusion-weighted imaging of the breast

OBJECTIVES Readout-segmented echo planar imaging (rs-EPI) significantly reduces susceptibility artifacts in diffusion-weighted imaging (DWI) of the breast compared to single-shot EPI but is limited by longer scan times. To compensate for this, we tested a new simultaneous multi-slice (SMS) acquisition for accelerated rs-EPI. MATERIALS AND METHODS After approval by the local ethics committee, eight healthy female volunteers (age, 38.9±13.1 years) underwent breast MRI at 3T. Conventional as well as two-fold (2× SMS) and three-fold (3× SMS) slice-accelerated rs-EPI sequences were acquired at b-values of 50 and 800s/mm(2). Two independent readers analyzed the apparent diffusion coefficient (ADC) in fibroglandular breast parenchyma. The signal-to-noise ratio (SNR) was estimated based on the subtraction method. ADC and SNR were compared between sequences by using the Friedman test. RESULTS The acquisition time was 4:21min for conventional rs-EPI, 2:35min for 2× SMS rs-EPI and 1:44min for 3× SMS rs-EPI. ADC values were similar in all sequences (mean values 1.62×10(-3)mm(2)/s, p=0.99). Mean SNR was 27.7-29.6, and no significant differences were found among the sequences (p=0.83). CONCLUSION SMS rs-EPI yields similar ADC values and SNR compared to conventional rs-EPI at markedly reduced scan time. Thus, SMS excitation increases the clinical applicability of rs-EPI for DWI of the breast.

Optimized on-line enantioselective capillary electrophoretic method for kinetic and inhibition studies of drug metabolism mediated by cytochrome P450 enzymes.

Pharmacokinetic and pharmacodynamic properties of a chiral drug can significantly differ between application of the racemate and single enantiomers. During drug development, the characteristics of candidate compounds have to be assessed prior to clinical testing. Since biotransformation significantly influences drug actions in an organism, metabolism studies represent a crucial part of such tests. Hence, an optimized and economical capillary electrophoretic method for on-line studies of the enantioselective drug metabolism mediated by cytochrome P450 enzymes was developed. It comprises a diffusion-based procedure, which enables mixing of the enzyme with virtually any compound inside the nanoliter-scale capillary reactor and without the need of additional optimization of mixing conditions. For CYP3A4, ketamine as probe substrate and highly sulfated γ-cyclodextrin as chiral selector, improved separation conditions for ketamine and norketamine enantiomers compared to a previously published electrophoretically mediated microanalysis method were elucidated. The new approach was thoroughly validated for the CYP3A4-mediated N-demethylation pathway of ketamine and applied to the determination of its kinetic parameters and the inhibition characteristics in presence of ketoconazole and dexmedetomidine. The determined parameters were found to be comparable to literature data obtained with different techniques. The presented method constitutes a miniaturized and cost-effective tool, which should be suitable for the assessment of the stereoselective aspects of kinetic and inhibition studies of cytochrome P450-mediated metabolic steps within early stages of the development of a new drug.

Interface discontinuity factors in the modal eigenspace of the multigroup diffusion matrix

Interface discontinuity factors based on the Generalized Equivalence Theory are commonly used in nodal homogenized diffusion calculations so that diffusion average values approximate heterogeneous higher order solutions. In this paper, an additional form of interface correction factors is presented in the frame of the Analytic Coarse Mesh Finite Difference Method (ACMFD), based on a correction of the modal fluxes instead of the physical fluxes. In the ACMFD formulation, implemented in COBAYA3 code, the coupled multigroup diffusion equations inside a homogenized region are reduced to a set of uncoupled modal equations through diagonalization of the multigroup diffusion matrix. Then, physical fluxes are transformed into modal fluxes in the eigenspace of the diffusion matrix. It is possible to introduce interface flux discontinuity jumps as the difference of heterogeneous and homogeneous modal fluxes instead of introducing interface discontinuity factors as the ratio of heterogeneous and homogeneous physical fluxes. The formulation in the modal space has been implemented in COBAYA3 code and assessed by comparison with solutions using classical interface discontinuity factors in the physical space

On the Extension of the Analytic Nodal Diffusion Solver ANDES to Sodium Fast Reactors

Within the framework of the Collaborative Project for a European Sodium Fast Reactor, the reactor physics group at UPM is working on the extension of its in-house multi-scale advanced deterministic code COBAYA3 to Sodium Fast Reactors (SFR). COBAYA3 is a 3D multigroup neutron kinetics diffusion code that can be used either as a pin-by-pin code or as a stand-alone nodal code by using the analytic nodal diffusion solver ANDES. It is coupled with thermalhydraulics codes such as COBRA-TF and FLICA, allowing transient analysis of LWR at both fine-mesh and coarse-mesh scales. In order to enable also 3D pin-by-pin and nodal coupled NK-TH simulations of SFR, different developments are in progress. This paper presents the first steps towards the application of COBAYA3 to this type of reactors. ANDES solver, already extended to triangular-Z geometry, has been applied to fast reactor steady-state calculations. The required cross section libraries were generated with ERANOS code for several configurations. The limitations encountered in the application of the Analytic Coarse Mesh Finite Difference (ACMFD) method –implemented inside ANDES– to fast reactors are presented and the sensitivity of the method when using a high number of energy groups is studied. ANDES performance is assessed by comparison with the results provided by ERANOS, using a mini-core model in 33 energy groups. Furthermore, a benchmark from the NEA for a small 3D FBR in hexagonal-Z geometry and 4 energy groups is also employed to verify the behavior of the code with few energy groups.

Neighborhood-corrected interface discontinuity factors for multi-group pin-by-pin diffusion calculations for LWR

Performing three-dimensional pin-by-pin full core calculations based on an improved solution of the multi-group diffusion equation is an affordable option nowadays to compute accurate local safety parameters for light water reactors. Since a transport approximation is solved, appropriate correction factors, such as interface discontinuity factors, are required to nearly reproduce the fully heterogeneous transport solution. Calculating exact pin-by-pin discontinuity factors requires the knowledge of the heterogeneous neutron flux distribution, which depends on the boundary conditions of the pin-cell as well as the local variables along the nuclear reactor operation. As a consequence, it is impractical to compute them for each possible configuration; however, inaccurate correction factors are one major source of error in core analysis when using multi-group diffusion theory. An alternative to generate accurate pin-by-pin interface discontinuity factors is to build a functional-fitting that allows incorporating the environment dependence in the computed values. This paper suggests a methodology to consider the neighborhood effect based on the Analytic Coarse-Mesh Finite Difference method for the multi-group diffusion equation. It has been applied to both definitions of interface discontinuity factors, the one based on the Generalized Equivalence Theory and the one based on Black-Box Homogenization, and for different few energy groups structures. Conclusions are drawn over the optimal functional-fitting and demonstrative results are obtained with the multi-group pin-by-pin diffusion code COBAYA3 for representative PWR configurations.

A Semi-Lagrangian Particle Level Set Finite Element Method for Interface Problems

We present a quasi-monotone semi-Lagrangian particle level set (QMSL-PLS) method for moving interfaces. The QMSL method is a blend of first order monotone and second order semi-Lagrangian methods. The QMSL-PLS method is easy to implement, efficient, and well adapted for unstructured, either simplicial or hexahedral, meshes. We prove that it is unconditionally stable in the maximum discrete norm, � · �h,∞, and the error analysis shows that when the level set solution u(t) is in the Sobolev space Wr+1,∞(D), r ≥ 0, the convergence in the maximum norm is of the form (KT/Δt)min(1,Δt � v �h,∞ /h)((1 − α)hp + hq), p = min(2, r + 1), and q = min(3, r + 1),where v is a velocity. This means that at high CFL numbers, that is, when Δt > h, the error is O( (1−α)hp+hq) Δt ), whereas at CFL numbers less than 1, the error is O((1 − α)hp−1 + hq−1)). We have tested our method with satisfactory results in benchmark problems such as the Zalesak’s slotted disk, the single vortex flow, and the rising bubble.

Stability of liquid bridges between an elliptical and a circular supporting disk

A numerical method has been developed to determine the stability limits for liquid bridges held between noncircular supporting disks and the application to a configuration with a circular and an elliptical disk subjected to axial acceleration has been made. The numerical method led to results very different from the available analytical solution which has been revisited and a better approximation has been obtained. It has been found that just retaining one more term in the asymptotic analysis the solution reproduces the real behavior of the configuration and the numerical results.

Image processing methods for human brain connectivity analysis from in-vivo diffusion MRI

The structural connectivity of the brain is considered to encode species-wise and subject-wise patterns that will unlock large areas of understanding of the human brain. Currently, diffusion MRI of the living brain enables to map the microstructure of tissue, allowing to track the pathways of fiber bundles connecting the cortical regions across the brain. These bundles are summarized in a network representation called connectome that is analyzed using graph theory. The extraction of the connectome from diffusion MRI requires a large processing flow including image enhancement, reconstruction, segmentation, registration, diffusion tracking, etc. Although a concerted effort has been devoted to the definition of standard pipelines for the connectome extraction, it is still crucial to define quality assessment protocols of these workflows. The definition of quality control protocols is hindered by the complexity of the pipelines under test and the absolute lack of gold-standards for diffusion MRI data. Here we characterize the impact on structural connectivity workflows of the geometrical deformation typically shown by diffusion MRI data due to the inhomogeneity of magnetic susceptibility across the imaged object. We propose an evaluation framework to compare the existing methodologies to correct for these artifacts including whole-brain realistic phantoms. Additionally, we design and implement an image segmentation and registration method to avoid performing the correction task and to enable processing in the native space of diffusion data. We release PySDCev, an evaluation framework for the quality control of connectivity pipelines, specialized in the study of susceptibility-derived distortions. In this context, we propose Diffantom, a whole-brain phantom that provides a solution to the lack of gold-standard data. The three correction methodologies under comparison performed reasonably, and it is difficult to determine which method is more advisable. We demonstrate that susceptibility-derived correction is necessary to increase the sensitivity of connectivity pipelines, at the cost of specificity. Finally, with the registration and segmentation tool called regseg we demonstrate how the problem of susceptibility-derived distortion can be overcome allowing data to be used in their original coordinates. This is crucial to increase the sensitivity of the whole pipeline without any loss in specificity.

Design of a Semi-analytical Method to Describe Plasma-Wall Interactions

Electric probes are objects immersed in the plasma with sharp boundaries which collect of emit charged particles. Consequently, the nearby plasma evolves under abrupt imposed and/or naturally emerging conditions. There could be localized currents, different time scales for plasma species evolution, charge separation and absorbing-emitting walls. The traditional numerical schemes based on differences often transform these disparate boundary conditions into computational singularities. This is the case of models using advection-diffusion differential equations with source-sink terms (also called Fokker-Planck equations). These equations are used in both, fluid and kinetic descriptions, to obtain the distribution functions or the density for each plasma species close to the boundaries. We present a resolution method grounded on an integral advancing scheme by using approximate Green's functions, also called short-time propagators. All the integrals, as a path integration process, are numerically calculated, what states a robust grid-free computational integral method, which is unconditionally stable for any time step. Hence, the sharp boundary conditions, as the current emission from a wall, can be treated during the short-time regime providing solutions that works as if they were known for each time step analytically. The form of the propagator (typically a multivariate Gaussian) is not unique and it can be adjusted during the advancing scheme to preserve the conserved quantities of the problem. The effects of the electric or magnetic fields can be incorporated into the iterative algorithm. The method allows smooth transitions of the evolving solutions even when abrupt discontinuities are present. In this work it is proposed a procedure to incorporate, for the very first time, the boundary conditions in the numerical integral scheme. This numerical scheme is applied to model the plasma bulk interaction with a charge-emitting electrode, dealing with fluid diffusion equations combined with Poisson equation self-consistently. It has been checked the stability of this computational method under any number of iterations, even for advancing in time electrons and ions having different time scales. This work establishes the basis to deal in future work with problems related to plasma thrusters or emissive probes in electromagnetic fields.

Sorting by diffusion: An asymmetric obstacle course for continuous molecular separation

A separation technique employing a microfabricated sieve has been demonstrated by observing the motion of DNA molecules of different size. The sieve consists of a two-dimensional lattice of obstacles whose asymmetric disposition rectifies the Brownian motion of molecules driven through the device, causing them to follow paths that depend on their diffusion coefficient. A nominal 6% resolution by length of DNA molecules in the size range 15–30 kbp may be achieved in a 4-inch (10-cm) silicon wafer. The advantage of this method is that samples can be loaded and sorted continuously, in contrast to the batch mode commonly used in gel electrophoresis.

“Coarse” stability and bifurcation analysis using time-steppers: A reaction-diffusion example

Evolutionary, pattern forming partial differential equations (PDEs) are often derived as limiting descriptions of microscopic, kinetic theory-based models of molecular processes (e.g., reaction and diffusion). The PDE dynamic behavior can be probed through direct simulation (time integration) or, more systematically, through stability/bifurcation calculations; time-stepper-based approaches, like the Recursive Projection Method [Shroff, G. M. & Keller, H. B. (1993) SIAM J. Numer. Anal. 30, 1099–1120] provide an attractive framework for the latter. We demonstrate an adaptation of this approach that allows for a direct, effective (“coarse”) bifurcation analysis of microscopic, kinetic-based models; this is illustrated through a comparative study of the FitzHugh-Nagumo PDE and of a corresponding Lattice–Boltzmann model.

Co-Permeability of 3H-Labeled Water and 14C-Labeled Organic Acids across Isolated Plant Cuticles1 : Investigating Cuticular Paths of Diffusion and Predicting Cuticular Transpiration

Penetration of 3H-labeled water (3H2O) and the 14C-labeled organic acids benzoic acid ([14C]BA), salicylic acid ([14C]SA), and 2,4-dichlorophenoxyacetic acid ([14C]2,4-D) were measured simultaneously in isolated cuticular membranes of Prunus laurocerasus L., Ginkgo biloba L., and Juglans regia L. For each of the three pairs of compounds (3H2O/[14C]BA, 3H2O/[14C]SA, and 3H2O/[14C]2,4-D) rates of cuticular water penetration were highly correlated with the rates of penetration of the organic acids. Therefore, water and organic acids penetrated the cuticles by the same routes. With the combination 3H2O/[14C]BA, co-permeability was measured with isolated cuticles of nine other plant species. Permeances of 3H2O of all 12 investigated species were highly correlated with the permeances of [14C]BA (r2 = 0.95). Thus, cuticular transpiration can be predicted from BA permeance. The application of this experimental method, together with the established prediction equation, offers the opportunity to answer several important questions about cuticular transport physiology in future investigations.