A fully adaptive numerical approximation for a two-dimensional epidemic model with nonlinear cross-diffusion

An epidemic model is formulated by a reactionâeuro"diffusion system where the spatial pattern formation is driven by cross-diffusion. The reaction terms describe the local dynamics of susceptible and infected species, whereas the diffusion terms account for the spatial distribution dynamics. For both self-diffusion and cross-diffusion, nonlinear constitutive assumptions are suggested. To simulate the pattern formation two finite volume formulations are proposed, which employ a conservative and a non-conservative discretization, respectively. An efficient simulation is obtained by a fully adaptive multiresolution strategy. Numerical examples illustrate the impact of the cross-diffusion on the pattern formation.

Hypertension is an independent predictor of mean platelet volume in patients with acute ischaemic stroke.

Introduction: Mean platelet volume (MPV) was shown to be significantly increased in patients with acute ischaemic stroke, especially in non-lacunar strokes. Moreover, some studies concluded that increased MPV is related to poor functional outcome after ischaemic stroke, although this association is still controversial. However, the determinants of MPV in patients with acute ischaemic stroke have never been investigated. Subjects and methods: We recorded the main demographic, clinical and laboratory data of consecutive patients with acute (admitted within 24 h after stroke onset) ischaemic stroke admitted in our Neurology Service between January 2003 and December 2008. MPV was generated at admission by the Sysmex XE-2100 automated cell counter (Sysmex Corporation, Kobe, Japan) from ethylenediaminetetraacetic acid blood samples stored at room temperature until measurement. The association of these parameters with MPV was investigated in univariate and multivariate analysis. Results: A total of 636 patients was included in our study. The median MPV was 10.4 ± 0.82 fL. In univariate analysis, glucose (β= 0.03, P= 0.05), serum creatinine (β= 0.002, P= 0.02), haemoglobin (β= 0.009, P < 0.001), platelet count (β=-0.002, P < 0.001) and history of arterial hypertension (β= 0.21, P= 0.005) were found to be significantly associated with MPV. In multivariate robust regression analysis, only hypertension and platelet count remained as independent determinants of MPV. Conclusions: In patients with acute ischaemic stroke, platelet count and history of hypertension are the only determinants of MPV.

Regional specificity of MRI contrast parameter changes in normal ageing revealed by voxel-based quantification (VBQ).

Normal ageing is associated with characteristic changes in brain microstructure. Although in vivo neuroimaging captures spatial and temporal patterns of age-related changes of anatomy at the macroscopic scale, our knowledge of the underlying (patho)physiological processes at cellular and molecular levels is still limited. The aim of this study is to explore brain tissue properties in normal ageing using quantitative magnetic resonance imaging (MRI) alongside conventional morphological assessment. Using a whole-brain approach in a cohort of 26 adults, aged 18-85years, we performed voxel-based morphometric (VBM) analysis and voxel-based quantification (VBQ) of diffusion tensor, magnetization transfer (MT), R1, and R2* relaxation parameters. We found age-related reductions in cortical and subcortical grey matter volume paralleled by changes in fractional anisotropy (FA), mean diffusivity (MD), MT and R2*. The latter were regionally specific depending on their differential sensitivity to microscopic tissue properties. VBQ of white matter revealed distinct anatomical patterns of age-related change in microstructure. Widespread and profound reduction in MT contrasted with local FA decreases paralleled by MD increases. R1 reductions and R2* increases were observed to a smaller extent in overlapping occipito-parietal white matter regions. We interpret our findings, based on current biophysical models, as a fingerprint of age-dependent brain atrophy and underlying microstructural changes in myelin, iron deposits and water. The VBQ approach we present allows for systematic unbiased exploration of the interaction between imaging parameters and extends current methods for detection of neurodegenerative processes in the brain. The demonstrated parameter-specific distribution patterns offer insights into age-related brain structure changes in vivo and provide essential baseline data for studying disease against a background of healthy ageing.

Thalamic connectivity in patients with essential tremor treated with MR imaging-guided focused ultrasound: in vivo fiber tracking by using diffusion-tensor MR imaging.

PURPOSE: To use diffusion-tensor (DT) magnetic resonance (MR) imaging in patients with essential tremor who were treated with transcranial MR imaging-guided focused ultrasound lesion inducement to identify the structural connectivity of the ventralis intermedius nucleus of the thalamus and determine how DT imaging changes correlated with tremor changes after lesion inducement. MATERIALS AND METHODS: With institutional review board approval, and with prospective informed consent, 15 patients with medication-refractory essential tremor were enrolled in a HIPAA-compliant pilot study and were treated with transcranial MR imaging-guided focused ultrasound surgery targeting the ventralis intermedius nucleus of the thalamus contralateral to their dominant hand. Fourteen patients were ultimately included. DT MR imaging studies at 3.0 T were performed preoperatively and 24 hours, 1 week, 1 month, and 3 months after the procedure. Fractional anisotropy (FA) maps were calculated from the DT imaging data sets for all time points in all patients. Voxels where FA consistently decreased over time were identified, and FA change in these voxels was correlated with clinical changes in tremor over the same period by using Pearson correlation. RESULTS: Ipsilateral brain structures that showed prespecified negative correlation values of FA over time of -0.5 or less included the pre- and postcentral subcortical white matter in the hand knob area; the region of the corticospinal tract in the centrum semiovale, in the posterior limb of the internal capsule, and in the cerebral peduncle; the thalamus; the region of the red nucleus; the location of the central tegmental tract; and the region of the inferior olive. The contralateral middle cerebellar peduncle and bilateral portions of the superior vermis also showed persistent decrease in FA over time. There was strong correlation between decrease in FA and clinical improvement in hand tremor 3 months after lesion inducement (P < .001). CONCLUSION: DT MR imaging after MR imaging-guided focused ultrasound thalamotomy depicts changes in specific brain structures. The magnitude of the DT imaging changes after thalamic lesion inducement correlates with the degree of clinical improvement in essential tremor.

Double-tagged fluorescent bacterial bioreporter for the study of polycyclic aromatic hydrocarbon diffusion and bioavailability.

Bacterial degradation of polycyclic aromatic hydrocarbons (PAHs), ubiquitous contaminants from oil and coal, is typically limited by poor accessibility of the contaminant to the bacteria. In order to measure PAH availability in complex systems, we designed a number of diffusion-based assays with a double-tagged bacterial reporter strain Burkholderia sartisoli RP037-mChe. The reporter strain is capable of mineralizing phenanthrene (PHE) and induces the expression of enhanced green fluorescent protein (eGFP) as a function of the PAH flux to the cell. At the same time, it produces a second autofluorescent protein (mCherry) in constitutive manner. Quantitative epifluorescence imaging was deployed in order to record reporter signals as a function of PAH availability. The reporter strain expressed eGFP proportionally to dosages of naphthalene or PHE in batch liquid cultures. To detect PAH diffusion from solid materials the reporter cells were embedded in 2 cm-sized agarose gel patches, and fluorescence was recorded over time for both markers as a function of distance to the PAH source. eGFP fluorescence gradients measured on known amounts of naphthalene or PHE served as calibration for quantifying PAH availability from contaminated soils. To detect reporter gene expression at even smaller diffusion distances, we mixed and immobilized cells with contaminated soils in an agarose gel. eGFP fluorescence measurements confirmed gel patch diffusion results that exposure to 2-3 mg lampblack soil gave four times higher expression than to material contaminated with 10 or 1 (mg PHE) g(-1).

Le problème de l'interdépendance dans la comparaison des politiques publiques infranationales : une méta-analyse des mécanismes de diffusion

Un des problèmes majeurs de la méthode comparative, et notamment de la comparaison entre politiques publiques infranationales, est le fait que les phénomènes à comparer, au lieu de consister en unités analytiques indépendantes les unes des autres, sont au contraire susceptibles de s'influencer mutuellement. Cette contribution aborde la problématique de l'interdépendance, avec une attention particulière aux différences entre politiques publiques internationales et infranationales. Elle synthétisera et discutera les apports des travaux empiriques sur la diffusion, à travers une métaanalyse des « conditions nécessaires » pour les différents mécanismes de diffusion des politiques publiques. Dans ce but, elle analysera, un échantillon de la littérature scientifique avec la technique des ensembles flous, et illustré à l'aide d'exemples concrets.

High b-value diffusion-weighted imaging: A sensitive method to reveal white matter differences in schizophrenia.

Over the last 10 years, diffusion-weighted imaging (DWI) has become an important tool to investigate white matter (WM) anomalies in schizophrenia. Despite technological improvement and the exponential use of this technique, discrepancies remain and little is known about optimal parameters to apply for diffusion weighting during image acquisition. Specifically, high b-value diffusion-weighted imaging known to be more sensitive to slow diffusion is not widely used, even though subtle myelin alterations as thought to happen in schizophrenia are likely to affect slow-diffusing protons. Schizophrenia patients and healthy controls were scanned with a high b-value (4000s/mm(2)) protocol. Apparent diffusion coefficient (ADC) measures turned out to be very sensitive in detecting differences between schizophrenia patients and healthy volunteers even in a relatively small sample. We speculate that this is related to the sensitivity of high b-value imaging to the slow-diffusing compartment believed to reflect mainly the intra-axonal and myelin bound water pool. We also compared these results to a low b-value imaging experiment performed on the same population in the same scanning session. Even though the acquisition protocols are not strictly comparable, we noticed important differences in sensitivities in the favor of high b-value imaging, warranting further exploration.

Tc-99m-MAA SPECT-derived tumor-to-normal liver ratios for partition model dosimetry in selective internal radiation therapy (SIRT)

Aim: When planning SIRT using 90Y microspheres, the partition model is used to refine the activity calculated by the body surface area (BSA) method to potentially improve the safety and efficacy of treatment. For this partition model dosimetry, accurate determination of mean tumor-to-normal liver ratio (TNR) is critical since it directly impacts absorbed dose estimates. This work aimed at developing and assessing a reliable methodology for the calculation of 99mTc-MAA SPECT/CT-derived TNR ratios based on phantom studies. Materials and methods: IQ NEMA (6 hot spheres) and Kyoto liver phantoms with different hot/background activity concentration ratios were imaged on a SPECT/CT (GE Infinia Hawkeye 4). For each reconstruction with the IQ phantom, TNR quantification was assessed in terms of relative recovery coefficients (RC) and image noise was evaluated in terms of coefficient of variation (COV) in the filled background. RCs were compared using OSEM with Hann, Butterworth and Gaussian filters, as well as FBP reconstruction algorithms. Regarding OSEM, RCs were assessed by varying different parameters independently, such as the number of iterations (i) and subsets (s) and the cut-off frequency of the filter (fc). The influence of the attenuation and diffusion corrections was also investigated. Furthermore, both 2D-ROIs and 3D-VOIs contouring were compared. For this purpose, dedicated Matlab© routines were developed in-house for automatic 2D-ROI/3D-VOI determination to reduce intra-user and intra-slice variability. Best reconstruction parameters and RCs obtained with the IQ phantom were used to recover corrected TNR in case of the Kyoto phantom for arbitrary hot-lesion size. In addition, we computed TNR volume histograms to better assess uptake heterogeneityResults: The highest RCs were obtained with OSEM (i=2, s=10) coupled with the Butterworth filter (fc=0.8). Indeed, we observed a global 20% RC improvement over other OSEM settings and a 50% increase as compared to the best FBP reconstruction. In any case, both attenuation and diffusion corrections must be applied, thus improving RC while preserving good image noise (COV<10%). Both 2D-ROI and 3D-VOI analysis lead to similar results. Nevertheless, we recommend using 3D-VOI since tumor uptake regions are intrinsically 3D. RC-corrected TNR values lie within 17% around the true value, substantially improving the evaluation of small volume (<15 mL) regions. Conclusions: This study reports the multi-parameter optimization of 99mTc MAA SPECT/CT images reconstruction in planning 90Y dosimetry for SIRT. In phantoms, accurate quantification of TNR was obtained using OSEM coupled with Butterworth and RC correction.

l0-deconvolution for compressive diffusion MRI

Diffusion MRI is a well established imaging modality providing a powerful way to non-invasively probe the structure of the white matter. Despite the potential of the technique, the intrinsic long scan times of these sequences have hampered their use in clinical practice. For this reason, a wide variety of methods have been proposed to shorten acquisition times. [...] We here review a recent work where we propose to further exploit the versatility of compressed sensing and convex optimization with the aim to characterize the fiber orientation distribution sparsity more optimally. We re-formulate the spherical deconvolution problem as a constrained l0 minimization.

Diffusion tensor imaging shows structural remodeling of stroke mirror region: results from a pilot study.

Background: The role of the non-injured hemisphere in stroke recovery is poorly understood. In this pilot study, we sought to explore the presence of structural changes detectable by diffusion tensor imaging (DTI) in the contralesional hemispheres of patients who recovered well from ischemic stroke. Methods: We analyzed serial DTI data from 16 stroke patients who had moderate initial neurological deficits (NIHSS scores 3-12) and good functional outcome at 3-6 months (NIHSS score 0 or modified Rankin Score ≤1). We segmented the brain tissue in gray and white matter (GM and WM) and measured the apparent diffusion coefficient (ADC) and fractional anisotropy in the infarct, in the contralesional infarct mirror region as well as in concentrically expanding regions around them. Results: We found that GM and WM ADC significantly increased in the infarct region (p < 0.01) from acute to chronic time points, whereas in the infarct mirror region, GM and WM ADC increased (p < 0.01) and WM fractional anisotropy decreased (p < 0.05). No significant changes were detected in other regions. Conclusion: DTI-based metrics are sensitive to regional structural changes in the contralesional hemisphere during stroke recovery. Prospective studies in larger cohorts with varying levels of recovery are needed to confirm our findings.

Ultra-High Field (7 Tesla) MRI for Gamma Knife Surgery Targeting of the Ventro-Intermediate Nucleus: A Pilot in Vivo Study

Aim: Gamma Knife surgery (GKS) is a non-invasive neurosurgical stereotactic procedure, increasingly used as an alternative to open functional procedures. This includes the targeting of the ventro-intermediate (Vim) nucleus of the thalamus for tremor. We currently perform an indirect targeting, using the "quadrilatere of Guyot," as the Vim nucleus is not visible on current 3 Tesla (T) MRI acquisitions. The primary objective of the current study was to enhance anatomic imaging for Vim GKS using high-field (7 T) MRI, with the aim of refining the visualization and precision of anatomical targeting. Method: Five young healthy subjects (mean age 23 years) were scanned both on 3 and 7 T MRI in Lausanne University Hospital (CHUV) and Center for Biomedical Imaging (CIBM). Classical T1-weighted MPRAGE, T2 CISS sequences (replacing former ventriculography) and diffusion tensor imaging were acquired at 3T. We obtained high-resolution susceptibility weighted images (SWI) at 7T for the visualization of thalamic subparts. SWI was further integrated for the first time into Leksell Gamma Plan® (LGP) software and co-registered with the 3T images. A simulation of targeting of the Vim was done using the "quadrilatere of Guyot" methodology on the 3T images. Furthermore, a correlation with the position of the found target on SWI was performed. The atlas of Morel et al. was used to confirm the findings on a detailed computer analysis outside LGP. Also, 3T and 7T MRI of one patient undergoing GKS Vim thalamotomy, were obtained before and 2 years after the procedure, and studied similarly. Results: The use of SWI provided a superior resolution and improved image contrast within the central gray matter. This allowed visualization and direct delineation of groups of thalamic nuclei in vivo, including the Vim. The position of the target, as assessed with the "quadrilatere of Guyot" method on 3 T, perfectly matched with the supposed one of the Vim on the SWI. Furthermore, a 3-dimensional model of the Vim target area was created on the basis of 3T and 7T images. Conclusion: This is the first report of the integration of SWI high-field MRI into the LGP in healthy subjects and in one patient treated GKS Vim thalamotomy. This approach aims at the improvement of targeting validation and further direct targeting of the Vim in tremor. The anatomical correlation between the direct visualization on 7T and the current targeting methods on 3T seems to show a very good anatomical matching.

RNA interference against aquaporin-4 decreases the apparent diffusion coefficient in the normal brain

Aim: Diffusion weighted magnetic resonance imaging (MRI) is now widely used in human brain diagnosis.1 To date molecular mechanisms underlying changes in Apparent Diffusion Coefficient (ADC) signals remain poorly understood. AQP4, localized to astrocytes, is one of the most highly expressed cerebral AQPs.2 AQP4 is involved in water movement within the cell membrane of cultured astrocytes.3 We hypothesize that AQP4 contributes to water diffusion and underlying ADC values in normal brain. Methods: We used an RNA interference (RNAi) protocol in vivo, to acutely knockdown expression of AQP4 in rat brain and to determine whether this was associated with changes in brain ADC values using MRI protocols as previously described.4 RNAi was performed using specific small interference RNA (siRNA) against AQP4 (siAQP4) and a non-targeted-siRNA (siGLO) as a control. The specificity and efficiency of the siAQP4 were first tested in vitro in astrocyte and hippocampal slice cultures. In vivo, siRNAs were injected into the rat cortex 3d prior to MRI acquisition and AQP4 was assessed by western blot (n=4) and immunohistochemistry (n=6). Histology was performed on adjacent slices. Results: siAQP4 application on primary astrocyte cultures induced a 76% decrease in AQP4 expression after 4 days. In hippocampal slice cultures; we also found a significant decrease in AQP4 expression in astrocytes after siAQP4. In vivo, injection of non-targeted siRNA (siGLO) tagged with CY3 allowed us to show that GFAP positive cells (astrocytes) were positively stained with CY3-siGLO, showing efficient transfection. Western blot and immunohistochemical analysis showed that siAQP4 induced a ~30% decrease in AQP4 expression without modification of tissue properties or cell death. After siAQP4 treatment, a significant decrease in ADC values (~50%) were observed without altered of T2 values. Conclusions: Together these results suggest that AQP4 reduces water diffusion through the astrocytic plasma membrane and decreases ADC values. Our findings demonstrate for the first time that astrocytic AQP4 contributes significantly to brain water diffusion and ADC values in normal brain. These results open new avenues to interpretation of ADC values under normal physiological conditions and in acute and chronic brain injuries.