Prognostic accuracy of cerebral blood flow measurement by perfusion computed tomography, at the time of emergency room admission, in acute stroke patients.

The purpose of this study was to determine the prognostic accuracy of perfusion computed tomography (CT), performed at the time of emergency room admission, in acute stroke patients. Accuracy was determined by comparison of perfusion CT with delayed magnetic resonance (MR) and by monitoring the evolution of each patient's clinical condition. Twenty-two acute stroke patients underwent perfusion CT covering four contiguous 10mm slices on admission, as well as delayed MR, performed after a median interval of 3 days after emergency room admission. Eight were treated with thrombolytic agents. Infarct size on the admission perfusion CT was compared with that on the delayed diffusion-weighted (DWI)-MR, chosen as the gold standard. Delayed magnetic resonance angiography and perfusion-weighted MR were used to detect recanalization. A potential recuperation ratio, defined as PRR = penumbra size/(penumbra size + infarct size) on the admission perfusion CT, was compared with the evolution in each patient's clinical condition, defined by the National Institutes of Health Stroke Scale (NIHSS). In the 8 cases with arterial recanalization, the size of the cerebral infarct on the delayed DWI-MR was larger than or equal to that of the infarct on the admission perfusion CT, but smaller than or equal to that of the ischemic lesion on the admission perfusion CT; and the observed improvement in the NIHSS correlated with the PRR (correlation coefficient = 0.833). In the 14 cases with persistent arterial occlusion, infarct size on the delayed DWI-MR correlated with ischemic lesion size on the admission perfusion CT (r = 0.958). In all 22 patients, the admission NIHSS correlated with the size of the ischemic area on the admission perfusion CT (r = 0.627). Based on these findings, we conclude that perfusion CT allows the accurate prediction of the final infarct size and the evaluation of clinical prognosis for acute stroke patients at the time of emergency evaluation. It may also provide information about the extent of the penumbra. Perfusion CT could therefore be a valuable tool in the early management of acute stroke patients.

A stable and 40Ar/39Ar isotope study of a major thrust in the Helvetic nappes (Swiss Alps) - Evidence for fluid-flow and constraints on nappe kinematics

Stable isotope and Ar-40/Ar-39 measurements,were made on samples associated with a major tectonic discontinuity in the Helvetic Alps, the basal thrust of the Diablerets nappe (external zone of the Alpine Belt) in order to determine both the importance of fluids in this thrust zone and the timing of thrusting. A systematic decrease in the delta(18)O values (up to 6 parts per thousand) of calcite, quartz, and white mica exists within a 10- to 70-m-wide zone over a distance of 37 km along the thrust, and they become more pronounced toward the root of the nappe. A similar decrease in the delta(13)C values of calcite is observed only in the deepest sections (up to 3 parts per thousand). The delta D-SMOW (SMOW = standard mean ocean water) values of white mica are -54 parts per thousand +/- 8 parts per thousand (n = 22) and are independent of the distance from the thrust. These variations are interpreted to reflect syntectonic solution reprecipitation during fluid passage along the thrust. The calculated delta(18)O and delta D values (versus SMOW) for the fluid in equilibrium with the analyzed minerals is 12 parts per thousand to 16 parts per thousand and -30 parts per thousand to +5 parts per thousand, respectively, for assumed temperatures of 250 to 450 degrees C. The isotopic and structural data are consistent with fluids derived from the deep-seated roots of the Helvetic nappes where large volumes of Mesozoic sediments were metamorphosed to the amphibolite facies, It is suggested that connate and metamorphic waters, overpressured by rapid tectonic burial in a subductive system escaped by upward infiltration along moderately dipping pathways until they reached the main shear zone at the base of the moving pile, where they were channeled toward the surface, This model also explains the mechanism by which large amounts of fluid were removed from the Mesozoic sediments during Alpine metamorphism. White mica Ar-49/Ar-39 ages vary from 27 Ma far from the Diablerets thrust to 15 Ma along the thrust. An older component is observed in micas far from the thrust, interpreted as a detrital signature, and indicates that regional metamorphic temperatures were less than about 350 degrees C. The;plateau and near plateau ages nearest the thrust are consistent with either neocrystallization of white mica or argon loss by recrystallization during thrusting, which may have been enhanced in the zones of highest fluid flow. The 15 Ma Ar-40/Ar-39 age plateau measured on white mica sampled exactly on the thrust surface dates the end of both fluid flow and tectonic transport.

Inversion prepared coronary MR angiography: direct visualization of coronary blood flow.

PURPOSE: Visualization of coronary blood flow by means of a slice-selective inversion pre-pulse in concert with bright-blood coronary MRA. MATERIALS AND METHODS: Coronary magnetic resonance angiography (MRA) of the right coronary artery (RCA) was performed in eight healthy adult subjects on a 1.5 Tesla MR system (Gyroscan ACS-NT, Philips Medical Systems, Best, NL) using a free-breathing navigator-gated and cardiac-triggered 3D steady-state free-precession (SSFP) sequence with radial k-space sampling. Imaging was performed with and without a slice-selective inversion pre-pulse, which was positioned along the main axis of the coronary artery but perpendicular to the imaging volume. Objective image quality parameters such as SNR, CNR, maximal visible vessel length, and vessel border definition were analyzed. RESULTS: In contrast to conventional bright-blood 3D coronary MRA, the selective inversion pre-pulse provided a direct measure of coronary blood flow. In addition, CNR between the RCA and right ventricular blood pool was increased and the vessels had a tendency towards better delineation. Blood SNR and CNR between right coronary blood and epicardial fat were comparable in both sequences. CONCLUSION: The combination of a free-breathing navigator-gated and cardiac-triggered 3D SSFP sequence with a slice-selective inversion pre-pulse allows for direct and directional visualization of coronary blood flow with the additional benefit of improved contrast between coronary and right ventricular blood pool.

Ambient vertical flow in long-screen wells: a case study in the Fontainebleau Sands Aquifer (France)

A tritium (H-3) profile was constructed in a long-screened well (LSW) of the Fontainebleau Sands Aquifer (France), and the data were combined with temperature logs to gain insight into the potential effects of the ambient vertical flow (AVF) of water through the well on the natural aquifer stratification. AVF is commonly taken into account in wells located in fracture aquifers or intercepting two different aquifers with distinct hydraulic heads. However, due to the vertical hydraulic gradient of the flow lines intercepted by wells, AVF of groundwater is a common process within any type of aquifer. The detection of 3H in the deeper parts of the studied well ( approximate depth 50m), where H-3-free groundwater is expected, indicates that shallow young water is being transported downwards through the well itself. The temperature logs show a nearly zero gradient with depth, far below the mean geothermal gradient in sedimentary basins. The results show that the age distribution of groundwater samples might be biased in relation to the age distribution in the surroundings of the well. The use of environmental tracers to investigate aquifer properties, particularly in LSWs, is then limited by the effects of the AVF of water that naturally occurs through the well.

Study of gene flow through a hybrid zone in the common shrew (Sorex araneus) using microsatellites

The common shrew (Sorer araneus) is subdivided into several chromosomal races. As hybrid zones between them have been characterized, this organism is of particular interest in studying the role of chromosomes in speciation. Six microsatellite loci were used to evaluate the level of gene how in the S. araneus hybrid zone between the Cordon and Valais races. Most of these loci were very polymorphic, the total number of alleles detected per locus ranging from 3 to 20. Using Mantel tests, we showed that the effect of rivers as barriers to gene flow is less important at this sampling scale. The effect of the chromosomal race is of particular importantance in diminishing gene flow.

Volumetric quantification of atherosclerotic plaque in CT considering partial volume effect.

Coronary artery calcification (CAC) is quantified based on a computed tomography (CT) scan image. A calcified region is identified. Modified expectation maximization (MEM) of a statistical model for the calcified and background material is used to estimate the partial calcium content of the voxels. The algorithm limits the region over which MEM is performed. By using MEM, the statistical properties of the model are iteratively updated based on the calculated resultant calcium distribution from the previous iteration. The estimated statistical properties are used to generate a map of the partial calcium content in the calcified region. The volume of calcium in the calcified region is determined based on the map. The experimental results on a cardiac phantom, scanned 90 times using 15 different protocols, demonstrate that the proposed method is less sensitive to partial volume effect and noise, with average error of 9.5% (standard deviation (SD) of 5-7mm(3)) compared with 67% (SD of 3-20mm(3)) for conventional techniques. The high reproducibility of the proposed method for 35 patients, scanned twice using the same protocol at a minimum interval of 10 min, shows that the method provides 2-3 times lower interscan variation than conventional techniques.

The effects of sevoflurane and halothane anesthesia on cerebral blood flow velocity in children.

We compared cerebral blood flow velocity during anesthesia with sevoflurane and halothane in 23 children admitted for elective surgery (age, 0.4-9.7 yr; median age, 1.9 yr; ASA physical status I-II). Inhaled induction was performed in a randomized sequence with sevoflurane or halothane. Under steady-state conditions, cerebral blood flow velocity (systolic [V(s)], mean [V(mn)], and diastolic [VD]) were measured by a blinded investigator using transcranial pulsed Doppler ultrasonography. The anesthetic was then changed. CBFV measurements were repeated after washout of the first anesthetic and after steady-state of the second (equivalent minimal alveolar concentration to first anesthetic). The resistance index was calculated. VD and V(mn) were significantly lower during sevoflurane (V(mn) 1.35 m/s) than during halothane (V(mn) 1.50 m/s; P = 0.001), whereas V(s) was unchanged. The resistance index was lower during halothane (P < 0.001). Our results indicate lower vessel resistance and higher mean velocity during halothane than during sevoflurane. IMPLICATIONS: The mean cerebral blood flow velocity is significantly decreased in children during inhaled anesthesia with sevoflurane than during halothane. This might be relevant for the choice of anesthetic in children with risk of increased intracranial pressure, neurosurgery, or craniofacial osteotomies.

Quantification of dose nonuniformities by voxel-based dosimetry in patients receiving 90Y-ibritumomab-tiuxetan.

OBJECTIVE: To assess the impact of nonuniform dose distribution within lesions and tumor-involved organs of patients receiving Zevalin, and to discuss possible implications of equivalent uniform biological effective doses (EU-BED) on treatment efficacy and toxicity. MATLAB? -based software for voxel-based dosimetry was adopted for this purpose. METHODS: Eleven lesions from seven patients with either indolent or aggressive non-Hodgkin lymphoma were analyzed, along with four organs with disease. Absorbed doses were estimated by a direct integration of single-voxel kinetic data from serial tomographic images. After proper corrections, differential BED distributions and surviving cell fractions were estimated, allowing for the calculation of EU-BED. To quantify dose uniformity in each target area, a heterogeneity index was defined. RESULTS: Average doses were below those prescribed by conventional radiotherapy to eradicate lymphoma lesions. Dose heterogeneity and effect on tumor control varied among lesions, with no apparent relation to tumor mass. Although radiation doses to involved organs were safe, unexpected liver toxicity occurred in one patient who presented with a pattern of diffuse infiltration. CONCLUSION: Voxel-based dosimetry and radiobiologic modeling can be successfully applied to lesions and tumor-involved organs, representing a methodological advance over estimation of mean absorbed doses. However, effects on tumor control and organ toxicity still cannot be easily predicted.

Oligo-Miocene extensional tectonics and fluid flow across the Northern Snake Range detachment system, Nevada

The Northern Snake Range (Nevada) represents a spectacular example of a metamorphic core complex and exposes a complete section from the mylonitic footwall into the hanging wall of a fossil detachment system. Paired geochronological and stable isotopic data of mylonitic quartzite within the detachment footwall reveal that ductile deformation and infiltration of meteoric fluids occurred between 27 and 23 Ma. Ar-40/Ar-39 ages display complex recrystallization-cooling relationships but decrease systematically from 26.9 +/- 0.2 Ma at the top to 21.3 +/- 0.2 Ma at the bottom of footwall mylonite. Hydrogen isotope (delta D) values in white mica are very low (-150 to -145 %) within the top 80-90 m of detachment footwall, in contrast to values obtained from the deeper part of the section where values range from -77 to -64 %, suggesting that time-integrated interaction between rock and meteoric fluid was restricted to the uppermost part of the mylonitic footwall. Pervasive mica-water hydrogen isotope exchange is difficult to reconcile with models of Ar-40 loss during mylonitization solely by volume diffusion. Rather, we interpret the Ar-40/Ar-39 ages of white mica with low-delta D values to date syn-mylonitic hydrogen and argon isotope exchange, and we conclude that the hydrothermal system of the Northern Snake Range was active during late Oligocene (27-23 Ma) and has been exhumed by the combined effects of ductile strain, extensional detachment faulting, and erosion.

Hybrid Multiscale Finite Volume method for two-phase flow in porous media

We present a novel hybrid (or multiphysics) algorithm, which couples pore-scale and Darcy descriptions of two-phase flow in porous media. The flow at the pore-scale is described by the Navier?Stokes equations, and the Volume of Fluid (VOF) method is used to model the evolution of the fluid?fluid interface. An extension of the Multiscale Finite Volume (MsFV) method is employed to construct the Darcy-scale problem. First, a set of local interpolators for pressure and velocity is constructed by solving the Navier?Stokes equations; then, a coarse mass-conservation problem is constructed by averaging the pore-scale velocity over the cells of a coarse grid, which act as control volumes; finally, a conservative pore-scale velocity field is reconstructed and used to advect the fluid?fluid interface. The method relies on the localization assumptions used to compute the interpolators (which are quite straightforward extensions of the standard MsFV) and on the postulate that the coarse-scale fluxes are proportional to the coarse-pressure differences. By numerical simulations of two-phase problems, we demonstrate that these assumptions provide hybrid solutions that are in good agreement with reference pore-scale solutions and are able to model the transition from stable to unstable flow regimes. Our hybrid method can naturally take advantage of several adaptive strategies and allows considering pore-scale fluxes only in some regions, while Darcy fluxes are used in the rest of the domain. Moreover, since the method relies on the assumption that the relationship between coarse-scale fluxes and pressure differences is local, it can be used as a numerical tool to investigate the limits of validity of Darcy's law and to understand the link between pore-scale quantities and their corresponding Darcy-scale variables.

Heterogeneity of human monocytes: an optimized four-color flow cytometry protocol for analysis of monocyte subsets.

Monocytes are central mediators in the development of atherosclerotic plaques. They circulate in blood and eventually migrate into tissue including the vessel wall where they give rise to macrophages and dendritic cells. The existence of monocyte subsets with distinct roles in homeostasis and inflammation suggests specialization of function. These subsets are identified based on expression of the CD14 and CD16 markers. Routinely applicable protocols remain elusive, however. Here, we present an optimized four-color flow cytometry protocol for analysis of human blood monocyte subsets using a specific PE-Cy5-conjugated monoclonal antibody (mAb) to HLA-DR, a PE-Cy7-conjugated mAb to CD14, a FITC-conjugated mAb to CD16, and PE-conjugated mAbs to additional markers relevant to monocyte function. Classical CD14(+)CD16(-) monocytes (here termed "Mo1" subset) expressed high CCR2, CD36, CD64, and CD62L, but low CX(3)CR1, whereas "nonclassical" CD14(lo)CD16(+) monocytes (Mo3) essentially showed the inverse expression pattern. CD14(+)CD16(+) monocytes (Mo2) expressed high HLA-DR, CD36, and CD64. In patients with stable coronary artery disease (n = 13), classical monocytes were decreased, whereas "nonclassical" monocytes were increased 90% compared with healthy subjects with angiographically normal coronary arteries (n = 14). Classical monocytes from CAD patients expressed higher CX(3)CR1 and CCR2 than controls. Thus, stable CAD is associated with expansion of the nonclassical monocyte subset and increased expression of inflammatory markers on monocytes. Flow cytometric analysis of monocyte subsets and marker expression may provide valuable information on vascular inflammation. This may translate into the identification of monocyte subsets as selective therapeutic targets, thus avoiding adverse events associated with indiscriminate monocyte inhibition.

Quantitative integration of geophysical and hydraulic data : from the local towards the regional scale range

Des progrès significatifs ont été réalisés dans le domaine de l'intégration quantitative des données géophysique et hydrologique l'échelle locale. Cependant, l'extension à de plus grandes échelles des approches correspondantes constitue encore un défi majeur. Il est néanmoins extrêmement important de relever ce défi pour développer des modèles fiables de flux des eaux souterraines et de transport de contaminant. Pour résoudre ce problème, j'ai développé une technique d'intégration des données hydrogéophysiques basée sur une procédure bayésienne de simulation séquentielle en deux étapes. Cette procédure vise des problèmes à plus grande échelle. L'objectif est de simuler la distribution d'un paramètre hydraulique cible à partir, d'une part, de mesures d'un paramètre géophysique pertinent qui couvrent l'espace de manière exhaustive, mais avec une faible résolution (spatiale) et, d'autre part, de mesures locales de très haute résolution des mêmes paramètres géophysique et hydraulique. Pour cela, mon algorithme lie dans un premier temps les données géophysiques de faible et de haute résolution à travers une procédure de réduction déchelle. Les données géophysiques régionales réduites sont ensuite reliées au champ du paramètre hydraulique à haute résolution. J'illustre d'abord l'application de cette nouvelle approche dintégration des données à une base de données synthétiques réaliste. Celle-ci est constituée de mesures de conductivité hydraulique et électrique de haute résolution réalisées dans les mêmes forages ainsi que destimations des conductivités électriques obtenues à partir de mesures de tomographic de résistivité électrique (ERT) sur l'ensemble de l'espace. Ces dernières mesures ont une faible résolution spatiale. La viabilité globale de cette méthode est testée en effectuant les simulations de flux et de transport au travers du modèle original du champ de conductivité hydraulique ainsi que du modèle simulé. Les simulations sont alors comparées. Les résultats obtenus indiquent que la procédure dintégration des données proposée permet d'obtenir des estimations de la conductivité en adéquation avec la structure à grande échelle ainsi que des predictions fiables des caractéristiques de transports sur des distances de moyenne à grande échelle. Les résultats correspondant au scénario de terrain indiquent que l'approche d'intégration des données nouvellement mise au point est capable d'appréhender correctement les hétérogénéitées à petite échelle aussi bien que les tendances à gande échelle du champ hydraulique prévalent. Les résultats montrent également une flexibilté remarquable et une robustesse de cette nouvelle approche dintégration des données. De ce fait, elle est susceptible d'être appliquée à un large éventail de données géophysiques et hydrologiques, à toutes les gammes déchelles. Dans la deuxième partie de ma thèse, j'évalue en détail la viabilité du réechantillonnage geostatique séquentiel comme mécanisme de proposition pour les méthodes Markov Chain Monte Carlo (MCMC) appliquées à des probmes inverses géophysiques et hydrologiques de grande dimension . L'objectif est de permettre une quantification plus précise et plus réaliste des incertitudes associées aux modèles obtenus. En considérant une série dexemples de tomographic radar puits à puits, j'étudie deux classes de stratégies de rééchantillonnage spatial en considérant leur habilité à générer efficacement et précisément des réalisations de la distribution postérieure bayésienne. Les résultats obtenus montrent que, malgré sa popularité, le réechantillonnage séquentiel est plutôt inefficace à générer des échantillons postérieurs indépendants pour des études de cas synthétiques réalistes, notamment pour le cas assez communs et importants où il existe de fortes corrélations spatiales entre le modèle et les paramètres. Pour résoudre ce problème, j'ai développé un nouvelle approche de perturbation basée sur une déformation progressive. Cette approche est flexible en ce qui concerne le nombre de paramètres du modèle et lintensité de la perturbation. Par rapport au rééchantillonage séquentiel, cette nouvelle approche s'avère être très efficace pour diminuer le nombre requis d'itérations pour générer des échantillons indépendants à partir de la distribution postérieure bayésienne. - Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending corresponding approaches beyond the local scale still represents a major challenge, yet is critically important for the development of reliable groundwater flow and contaminant transport models. To address this issue, I have developed a hydrogeophysical data integration technique based on a two-step Bayesian sequential simulation procedure that is specifically targeted towards larger-scale problems. The objective is to simulate the distribution of a target hydraulic parameter based on spatially exhaustive, but poorly resolved, measurements of a pertinent geophysical parameter and locally highly resolved, but spatially sparse, measurements of the considered geophysical and hydraulic parameters. To this end, my algorithm links the low- and high-resolution geophysical data via a downscaling procedure before relating the downscaled regional-scale geophysical data to the high-resolution hydraulic parameter field. I first illustrate the application of this novel data integration approach to a realistic synthetic database consisting of collocated high-resolution borehole measurements of the hydraulic and electrical conductivities and spatially exhaustive, low-resolution electrical conductivity estimates obtained from electrical resistivity tomography (ERT). The overall viability of this method is tested and verified by performing and comparing flow and transport simulations through the original and simulated hydraulic conductivity fields. The corresponding results indicate that the proposed data integration procedure does indeed allow for obtaining faithful estimates of the larger-scale hydraulic conductivity structure and reliable predictions of the transport characteristics over medium- to regional-scale distances. The approach is then applied to a corresponding field scenario consisting of collocated high- resolution measurements of the electrical conductivity, as measured using a cone penetrometer testing (CPT) system, and the hydraulic conductivity, as estimated from electromagnetic flowmeter and slug test measurements, in combination with spatially exhaustive low-resolution electrical conductivity estimates obtained from surface-based electrical resistivity tomography (ERT). The corresponding results indicate that the newly developed data integration approach is indeed capable of adequately capturing both the small-scale heterogeneity as well as the larger-scale trend of the prevailing hydraulic conductivity field. The results also indicate that this novel data integration approach is remarkably flexible and robust and hence can be expected to be applicable to a wide range of geophysical and hydrological data at all scale ranges. In the second part of my thesis, I evaluate in detail the viability of sequential geostatistical resampling as a proposal mechanism for Markov Chain Monte Carlo (MCMC) methods applied to high-dimensional geophysical and hydrological inverse problems in order to allow for a more accurate and realistic quantification of the uncertainty associated with the thus inferred models. Focusing on a series of pertinent crosshole georadar tomographic examples, I investigated two classes of geostatistical resampling strategies with regard to their ability to efficiently and accurately generate independent realizations from the Bayesian posterior distribution. The corresponding results indicate that, despite its popularity, sequential resampling is rather inefficient at drawing independent posterior samples for realistic synthetic case studies, notably for the practically common and important scenario of pronounced spatial correlation between model parameters. To address this issue, I have developed a new gradual-deformation-based perturbation approach, which is flexible with regard to the number of model parameters as well as the perturbation strength. Compared to sequential resampling, this newly proposed approach was proven to be highly effective in decreasing the number of iterations required for drawing independent samples from the Bayesian posterior distribution.

Quantification, self-renewal, and genetic tracing of FL1+ tumor-initiating cells in a large cohort of human gliomas.

Evidence has emerged that the initiation and growth of gliomas is sustained by a subpopulation of cancer-initiating cells (CICs). Because of the difficulty of using markers to tag CICs in gliomas, we have previously exploited more robust phenotypic characteristics, including a specific morphology and intrincic autofluorescence, to identify and isolate a subpopulation of glioma CICs, called FL1(+). The objective of this study was to further validate our method in a large cohort of human glioma and a mouse model of glioma. Seventy-four human gliomas of all grades and the GFAP-V(12)HA-ras B8 mouse model were analyzed for in vitro self-renewal capacity and their content of FL1(+). Nonneoplastic brain tissue and embryonic mouse brain were used as control. Genetic traceability along passages was assessed with microsatellite analysis. We found that FL1(+) cells from low-grade gliomas and from control nonneoplasic brain tissue show a lower level of autofluorescence and undergo a restricted number of cell divisions before dying in culture. In contrast, we found that FL1(+) cells derived from many but not all high-grade gliomas acquire high levels of autofluorescence and can be propagated in long-term cultures. Moreover, FL1(+) cells show a remarkable traceability over time in vitro and in vivo. Our results show that FL1(+) cells can be found in all specimens of a large cohort of human gliomas of different grades and in a model of genetically induced mouse glioma as well as nonneoplastic brain. However, their self-renewal capacity is variable and seems to be dependent on the tumor grade.

Error estimate and control in the MSFV method for multiphase flow in porous media

n this paper the iterative MSFV method is extended to include the sequential implicit simulation of time dependent problems involving the solution of a system of pressure-saturation equations. To control numerical errors in simulation results, an error estimate, based on the residual of the MSFV approximate pressure field, is introduced. In the initial time steps in simulation iterations are employed until a specified accuracy in pressure is achieved. This initial solution is then used to improve the localization assumption at later time steps. Additional iterations in pressure solution are employed only when the pressure residual becomes larger than a specified threshold value. Efficiency of the strategy and the error control criteria are numerically investigated. This paper also shows that it is possible to derive an a-priori estimate and control based on the allowed pressure-equation residual to guarantee the desired accuracy in saturation calculation.

Lesion Recognition in Multiple Sclerosis: A Sequence Comparison and Quantification Study at 3T

Introduction Lesion detection in multiple sclerosis (MS) is an essential part of its clinical diagnosis. In addition, radiological characterisation of MS lesions is an important research field that aims at distinguishing different MS types, monitoring drug response and prognosis. To date, various MR protocols have been proposed to obtain optimal lesion contrast for early and comprehensive diagnosis of the MS disease. In this study, we compare the sensitivity of five different MR contrasts for lesion detection: (i) the DIR sequence (Double Inversion Recovery, [4]), (ii) the Dark-fluid SPACE acquisition schemes, a 3D variant of a 2D FLAIR sequence [1], (iii) the MP2RAGE [2], an MP-RAGE variant that provides homogeneous T1 contrast and quantitative T1-values, and the sequences currently used for clinical MS diagnosis (2D FLAIR, MP-RAGE). Furthermore, we investigate the T1 relaxation times of cortical and sub-cortical regions in the brain hemispheres and the cerebellum at 3T. Methods 10 early-stage female MS patients (age: 31.64.7y; disease duration: 3.81.9y; disability score, EDSS: 1.80.4) and 10 healthy controls (age and gender-matched: 31.25.8y) were included in the study after obtaining informed written consent according to the local ethic protocol. All experiments were performed at 3T (Magnetom Trio a Tim System, Siemens, Germany) using a 32-channel head coil [5]. The imaging protocol included the following sequences, (all except for axial FLAIR 2D with 1x1x1.2 mm3 voxel and 256x256x160 matrix): DIR (TI1/TI2/TR XX/3652/10000 ms, iPAT=2, TA 12:02 min), MP-RAGE (TI/TR 900/2300 ms, iPAT=3, TA 3:47 min); MP2RAGE (TI1/TI2/TR 700/2500/5000 ms, iPAT=3, TA 8:22 min, cf. [2]); 3D FLAIR SPACE (only for patient 4-6, TI/TR 1800/5000 ms, iPAT=2, TA=5;52 min, cf. [1]); Axial FLAIR (0.9x0.9x2.5 mm3, 256x256x44 matrix, TI/TR 2500/9000 ms, iPAT=2, TA 4:05 min). Lesions were identified by two experienced neurologist and radiologist, manually contoured and assigned to regional locations (s. table 1). Regional lesion masks (RLM) from each contrast were compared for number and volumes of lesions. In addition, RLM were merged in a single "master" mask, which represented the sum of the lesions of all contrasts. T1 values were derived for each location from this mask for patients 5-10 (3D FLAIR contrast was missing for patient 1-4). Results & Discussion The DIR sequence appears the most sensitive for total lesions count, followed by the MP2RAGE (table 1). The 3D FLAIR SPACE sequence turns out to be more sensitive than the 2D FLAIR, presumably due to reduced partial volume effects. Looking for sub-cortical hemispheric lesions, the DIR contrast appears to be equally sensitive to the MP2RAGE and SPACE, but most sensitive for cerebellar MS plaques. The DIR sequence is also the one that reveals cortical hemispheric lesions best. T1 relaxation times at 3T in the WM and GM of the hemispheres and the cerebellum, as obtained with the MP2RAGE sequence, are shown in table 2. Extending previous studies, we confirm overall longer T1-values in lesion tissue and higher standard deviations compared to the non-lesion tissue and control tissue in healthy controls. We hypothesize a biological (different degree of axonal loss and demyelination) rather than technical origin. Conclusion In this study, we applied 5 MR contrasts including two novel sequences to investigate the contrast of highest sensitivity for early MS diagnosis. In addition, we characterized for the first time the T1 relaxation time in cortical and sub-cortical regions of the hemispheres and the cerebellum. Results are in agreement with previous publications and meaningful biological interpretation of the data.