Imaging and quantifying salt-tracer transport in a riparian groundwater system by means of 3D ERT monitoring

Determining groundwater flow paths of infiltrated river water is necessary for studying biochemical processes in the riparian zone, but their characterization is complicated by strong temporal and spatial heterogeneity. We investigated to what extent repeat 3D surface electrical resistance tomography (ERT) can be used to monitor transport of a salt-tracer plume under close to natural gradient conditions. The aim is to estimate groundwater flow velocities and pathways at a site located within a riparian groundwater system adjacent to the perialpine Thur River in northeastern Switzerland. Our ERT time-lapse images provide constraints on the plume's shape, flow direction, and velocity. These images allow the movement of the plume to be followed for 35 m. Although the hydraulic gradient is only 1.43 parts per thousand, the ERT time-lapse images demonstrate that the plume's center of mass and its front propagate with velocities of 2x10(-4) m/s and 5x10(-4) m/s, respectively. These velocities are compatible with groundwater resistivity monitoring data in two observation wells 5 m from the injection well. Five additional sensors in the 5-30 m distance range did not detect the plume. Comparison of the ERT time-lapse images with a groundwater transport model and time-lapse inversions of synthetic ERT data indicate that the movement of the plume can be described for the first 6 h after injection by a uniform transport model. Subsurface heterogeneity causes a change of the plume's direction and velocity at later times. Our results demonstrate the effectiveness of using time-lapse 3D surface ERT to monitor flow pathways in a challenging perialpine environment over larger scales than is practically possible with crosshole 3D ERT.

Brain damage in methylmalonic aciduria: 2-methylcitrate induces cerebral ammonium accumulation and apoptosis in 3D organotypic brain cell cultures.

BACKGROUND: Methylmalonic aciduria is an inborn error of metabolism characterized by accumulation of methylmalonate (MMA), propionate and 2-methylcitrate (2-MCA) in body fluids. Early diagnosis and current treatment strategies aimed at limiting the production of these metabolites are only partially effective in preventing neurological damage. METHODS: To explore the metabolic consequences of methylmalonic aciduria on the brain, we used 3D organotypic brain cell cultures from rat embryos. We challenged the cultures at two different developmental stages with 1 mM MMA, propionate or 2-MCA applied 6 times every 12 h. In a dose-response experiment cultures were challenged with 0.01, 0.1, 0.33 and 1 mM 2-MCA. Immunohistochemical staining for different brain cell markers were used to assess cell viability, morphology and differentiation. Significant changes were validated by western blot analysis. Biochemical markers were analyzed in culture media. Apoptosis was studied by immunofluorescence staining and western blots for activated caspase-3. RESULTS: Among the three metabolites tested, 2-MCA consistently produced the most pronounced effects. Exposure to 2-MCA caused morphological changes in neuronal and glial cells already at 0.01 mM. At the biochemical level the most striking result was a significant ammonium increase in culture media with a concomitant glutamine decrease. Dose-response studies showed significant and parallel changes of ammonium and glutamine starting from 0.1 mM 2-MCA. An increased apoptosis rate was observed by activation of caspase-3 after exposure to at least 0.1 mM 2-MCA. CONCLUSION: Surprisingly, 2-MCA, and not MMA, seems to be the most toxic metabolite in our in vitro model leading to delayed axonal growth, apoptosis of glial cells and to unexpected ammonium increase. Morphological changes were already observed at 2-MCA concentrations as low as 0.01 mM. Increased apoptosis and ammonium accumulation started at 0.1 mM thus suggesting that ammonium accumulation is secondary to cell suffering and/or cell death. Local accumulation of ammonium in CNS, that may remain undetected in plasma and urine, may therefore play a key role in the neuropathogenesis of methylmalonic aciduria both during acute decompensations and in chronic phases. If confirmed in vivo, this finding might shift the current paradigm and result in novel therapeutic strategies.

A MODEL OF CREATINE DEFICIENCY SYNDROMES IN 3D BRAIN CELL CULTURES BY KNOCKDOWN OF GAMT AND SLC6A8 GENES

La créatine joue un rôle essentiel dans le métabolisme cellulaire par sa conversion, par la creatine kinase, en phosphocreatine permettant la régénération de l'ATP. La synthèse de créatine, chez les mammifères, s'effectue par une réaction en deux étapes impliquant Γ arginine: glycine amidinotransférase (AGAT) et la guanidinoacétate méthyltransférase (GAMT). L'entrée de créatine dans les cellules s'effectue par son transporteur, SLC6A8. Les déficiences en créatine, dues au déficit en GAMT, AGAT ou SLC6A8, sont fréquentes et caractérisées par une absence ou une forte baisse de créatine dans le système nerveux central. Alors qu'il est connu que AGAT, GAMT et SLC6A8 sont exprimés par le cerveau, les conséquences des déficiences en créatine sur les cellules nerveuses sont peu comprises. Le but de ce travail était de développer de nouveaux modèles expérimentaux des déficiences en Cr dans des cultures 3D de cellules nerveuses de rat en agrégats au moyen de l'interférence à l'ARN appliquée aux gènes GAMT et SLC6A8. Des séquences interférentes (shRNAs) pour les gènes GAMT et SLC6A8 ont été transduites par des vecteurs viraux AAV (virus adéno-associés), dans les cellules nerveuses en agrégats. Nous avons ainsi démontré une baisse de l'expression de GAMT au niveau protéique (mesuré par western blot), et ARN messager (mesuré par qPCR) ainsi qu'une variation caractérisitique de créatine et guanidinoacétate (mesuré par spectrométrie de masse). Après avoir validé nos modèles, nous avons montré que les knockdown de GAMT ou SLC6A8 affectent le développement des astrocytes et des neurones ou des oligodendrocytes et des astrocytes, respectivement, ainsi qu'une augmentation de la mort cellulaire et des modifications dans le pattern d'activation des voies de signalisation impliquant caspase 3 et p38 MAPK, ayant un rôle dans le processus d'apoptose. - Creatine plays essential roles in energy metabolism by the interconversion, by creatine kinase, to its phosphorylated analogue, phosphocreatine, allowing the regeneration of ATP. Creatine is synthesized in mammals by a two step mechanism involving arginine:glycine amidinotransferase (AGAT) and guanidinoacetate methyltransferase (GAMT). Creatine is taken up by cells by a specific transporter, SLC6A8. Creatine deficiency syndromes, due to defects in GAMT, AGAT and SLC6A8, are among the most frequent inborn errors of metabolism, and are characterized by an absence or a severe decrease of creatine in central nervous system, which is the main tissue affected. While it is known that AGAT, GAMT and SLC6A8 are expressed in CNS, many questions remain on the specific effects of AGAT, GAMT and SLC6A8 deficiencies on brain cells. Our aim was to develop new experimental models of creatine deficiencies by knockdown of GAMT and SLC6A8 genes by RNAi in 3D organotypic rat brain cell cultures in aggregates. Specific shRNAs for the GAMT and SLC6A8 genes were transduced in brain cell aggregates by adeno-associated viruses (AAV). The AAV-transduced shRNAs were able to efficiently knockdown the expression of our genes of interest, as shown by a strong decrease of protein by western blotting, a decrease of mRNA by qPCR or characteristic variations of creatine and guanidinoacetate by tandem mass spectrometry. After having validated our experimental models, we have also shown that GAMT and SLC6A8 knockdown affected the development of astrocytes and neurons or oligodendrocytes and astrocytes, respectively. We also observed an increase of cell death and variations in activation pattern of caspase 3 and p38 MAPK pathways, involved in apoptosis, in our experimental model.

Helical tomotherapy (HT) for the treatment of anal canal cancer: preliminary clinical results, and dosimetric comparison between HT and intensity-modulated or 3D conformal radiotherapy

Background: To report a single-center experience in 19 patients (pts) with anal canal cancer treated with helical tomotherapy (HT) and concurrent chemotherapy, and compare the dosimetric results with fixed-field intensitymodulated radiotherapy (IMRT) and 3D conformal radiotherapy (3D RT). Materials and Methods: Between 2007 and 2008, 19 consecutive pts were treated with HT and concurrent CT for anal canal cancer. Median age was 59 years (range, 38−83), and female/male ratio was 14/5. The majority of the pts had T2 or T3 tumours (68.4%), and 52.6% had positive lymph nodes. In all 19 pts, pelvic and inguinal nodes, and tumour irradiation was given using HT upto a median dose of 36 Gy (1.8 Gy/fr) followed by a 1-week gap. A boost dose of 23.4 Gy (1.8 Gy/fr) was delivered to the tumour and involved nodes using 3DRT (n = 12), HT (n = 6), or IMRT (n = 1). Simultaneous integrated boost was used in none of the pts. All but one patient with a T1N0 tumour received concomitant mitomycin/5- fluorouracil (n = 12) or mitomycin/capecitabin (n = 7) CT. Toxicity was scored according to the Common Terminology Criteria for Adverse Events (NCICTCAE v3.0). HT plans and treatments were generated using Tomotherapy, Inc., software and hardware; and 3D or IMRT boost plans with the CMS treatment planning system (TPS), using 6−18 MV photons from a Siemens Primus accelerator. For dosimetric comparison, computed tomography data sets of 10 pts were imported into the TPS, and 3D and 5-field step-andshoot IMRT plans were generated for each case. Plans were optimized with the aim of assessing organs at risk (OAR) and healthy-tissue sparing while enforcing highly conformal target coverage, and evaluated by dose-volume histograms (DVH) of planning target volumes (PTV) and OAR. Results: With a median follow-up of 13 months (range, 3−18), all pts are alive and well; except one patient developing local recurrence at 12 months. No patient developed grade 3 or more acute toxicity. No unplanned treatment interruption was necessary because of toxicity. With 360-degree-of-freedom beam projection, HT showed an advantage over 3D or IMRT plans in terms of dose conformity around the PTV, and dose gradients were steeper outside the PTV, resulting in reduced doses to OARs. Using HT, acute toxicity was acceptable, and seemed to be better than historical standards. Conclusion: We conclude that HT combined with concurrent chemotherapy for anal canal cancer is effective and tolerable. Compared to 3DRT or 5-field IMRT, there is better conformity around the PTV, and OAR sparing.

Segmentation automatique de la matière grise et de la matière blanche par un algorithme génétique : une étude de validation

Résumé Objectif: l'observation des variations de volume de la matière grise (MG), de la matière blanche (MB), et du liquide céphalo-rachidien (LCR) est particulièrement utile dans l'étude de nombreux processus physiopathologiques, la mesure quantitative 'in vivo' de ces volumes présente un intérêt considérable tant en recherche qu'en pratique clinique. Cette étude présente et valide une méthode de segmentation automatique du cerveau avec mesure des volumes de MG et MB sur des images de résonance magnétique. Matériel et Méthode: nous utilisons un algorithme génétique automatique pour segmenter le cerveau en MG, MB et LCR à partir d'images tri-dimensionnelles de résonance magnétique en pondération Ti. Une étude morphométrique a été conduite sur 136 sujets hommes et femmes de 15 à 74 ans. L'algorithme a ensuite été validé par 5 approches différentes: I. Comparaison de mesures de volume sur un cerveau de cadavre par méthode automatique et par mesure de déplacement d'eau selon la méthode d'Archimède. 2. Comparaison de mesures surfaces sur des images bidimensionnelles segmentées soit par un traçage manuel soit par la méthode automatique. 3. Evaluation de la fiabilité de la segmentation par acquisitions et segmentations itératives du même cerveau. 4. Les volumes de MG, MB et LCR ont été utilisés pour une étude du vieillissement normal de la population. 5. Comparaison avec les données existantes de la littérature. Résultats: nous avons pu observer une variation de la mesure de 4.17% supplémentaire entre le volume d'un cerveau de cadavre mesuré par la méthode d'Archimède, en majeure partie due à la persistance de tissus après dissection_ La comparaison des méthodes de comptage manuel de surface avec la méthode automatique n'a pas montré de variation significative. L'épreuve du repositionnement du même sujet à diverses reprises montre une très bonne fiabilité avec une déviation standard de 0.46% pour la MG, 1.02% pour la MB et 3.59% pour le LCR, soit 0.19% pour le volume intracrânien total (VICT). L'étude morphométrique corrobore les résultats des études anatomiques et radiologiques existantes. Conclusion: la segmentation du cerveau par un algorithme génétique permet une mesure 100% automatique, fiable et rapide des volumes cérébraux in vivo chez l'individu normal.

Using 3D surface datasets to understand landslide evolution: From analogue models to real case study

Early detection of landslide surface deformation with 3D remote sensing techniques, as TLS, has become a great challenge during last decade. To improve our understanding of landslide deformation, a series of analogue simulation have been carried out on non-rigid bodies coupled with 3D digitizer. All these experiments have been carried out under controlled conditions, as water level and slope angle inclination. We were able to follow 3D surface deformation suffered by complex landslide bodies from precursory deformation still larger failures. These experiments were the basis for the development of a new algorithm for the quantification of surface deformation using automatic tracking method on discrete points of the slope surface. To validate the algorithm, comparisons were made between manually obtained results and algorithm surface displacement results. Outputs will help in understanding 3D deformation during pre-failure stages and failure mechanisms, which are fundamental aspects for future implementation of 3D remote sensing techniques in early warning systems.

High throughput 3D super-resolution microscopy reveals Caulobacter crescentus in vivo Z-ring organization.

We created a high-throughput modality of photoactivated localization microscopy (PALM) that enables automated 3D PALM imaging of hundreds of synchronized bacteria during all stages of the cell cycle. We used high-throughput PALM to investigate the nanoscale organization of the bacterial cell division protein FtsZ in live Caulobacter crescentus. We observed that FtsZ predominantly localizes as a patchy midcell band, and only rarely as a continuous ring, supporting a model of "Z-ring" organization whereby FtsZ protofilaments are randomly distributed within the band and interact only weakly. We found evidence for a previously unidentified period of rapid ring contraction in the final stages of the cell cycle. We also found that DNA damage resulted in production of high-density continuous Z-rings, which may obstruct cytokinesis. Our results provide a detailed quantitative picture of in vivo Z-ring organization.

Multimodal evaluation for medical image segmentation

This paper is a joint effort between five institutionsthat introduces several novel similarity measures andcombines them to carry out a multimodal segmentationevaluation. The new similarity measures proposed arebased on the location and the intensity values of themisclassified voxels as well as on the connectivity andthe boundaries of the segmented data. We showexperimentally that the combination of these measuresimprove the quality of the evaluation. The study that weshow here has been carried out using four differentsegmentation methods from four different labs applied toa MRI simulated dataset of the brain. We claim that ournew measures improve the robustness of the evaluation andprovides better understanding about the differencebetween segmentation methods.

Influence of the ankle arthropathies on the 3D kinematics of the lower limbs during gait: analysis with an ambulatory system.

Introduction: Ankle arthropathy is associated with a decreased motion of the ankle-hindfoot during ambulation. Ankle arthrodesis was shown to result in degeneration of the neighbour joints of the foot. Inversely, total ankle arthroplasty conceptually preserves the adjacent joints because of the residual mobility of the ankle but this has not been demonstrated yet in vivo. It has also been reported that degenerative ankle diseases, and even arthrodesis, do not result in alteration of the knee and hip joints. We present the preliminary results of a new approach of this problem based on ambulatory gait analysis. Patients and Methods: Motion analysis of the lower limbs was performed using a Physilog® (BioAGM, CH) system consisting of three-dimensional (3D) accelerometer and gyroscope, coupled to a magnetic system (Liberty©, Polhemus, USA). Both systems have been validated. Three groups of two patients were included into this pilot study and compared to healthy subjects (controls) during level walking: patients with ankle osteoarthritis (group 1), patients treated by ankle arthrodesis (group 2), patients treated by total ankle prosthesis (group 3). Results: Motion patterns of all analyzed joints over more than 20 gait cycles in each subject were highly repeatable. Motion amplitude of the ankle-hindfoot in control patients was similar to recently reported results. Ankle arthrodesis limited the motion of the ankle-hindfoot in the sagittal and horizontal planes. The prosthetic ankle allowed a more physiologic movement in the sagittal plane only. Ankle arthritis and its treatments did not influence the range of motion of the knee and hip joint during stance phase, excepted for a slight decrease of the hip flexion in groups 1 and 2. Conclusion: The reliability of the system was shown by the repeatability of the consecutive measurements. The results of this preliminary study were similar to those obtained through laboratory gait analysis. However, our system has the advantage to allow ambulatory analysis of 3D kinematics of the lower limbs outside of a gait laboratory and in real life conditions. To our knowledge this is a new concept in the analysis of ankle arthropathy and its treatments. Therefore, there is a potential to address specific questions like the difficult comparison of the benefits of ankle arthroplasty versus arthrodesis. The encouraging results of this pilot study offer the perspective to analyze the consequences of ankle arthropathy and its treatments on the biomechanics of the lower limbs ambulatory, in vivo and in daily life conditions.

Trabecular bone score (TBS) the new parameter of 2D texture analysis for the evaluation of 3D bone micro architecture status

X-ray is a technology that is used for numerous applications in the medical field. The process of X-ray projection gives a 2-dimension (2D) grey-level texture from a 3- dimension (3D) object. Until now no clear demonstration or correlation has positioned the 2D texture analysis as a valid indirect evaluation of the 3D microarchitecture. TBS is a new texture parameter based on the measure of the experimental variogram. TBS evaluates the variation between 2D image grey-levels. The aim of this study was to evaluate existing correlations between 3D bone microarchitecture parameters - evaluated from μCT reconstructions - and the TBS value, calculated on 2D projected images. 30 dried human cadaveric vertebrae were acquired on a micro-scanner (eXplorer Locus, GE) at isotropic resolution of 93 μm. 3D vertebral body models were used. The following 3D microarchitecture parameters were used: Bone volume fraction (BV/TV), Trabecular thickness (TbTh), trabecular space (TbSp), trabecular number (TbN) and connectivity density (ConnD). 3D/2D projections has been done by taking into account the Beer-Lambert Law at X-ray energy of 50, 100, 150 KeV. TBS was assessed on 2D projected images. Correlations between TBS and the 3D microarchitecture parameters were evaluated using a linear regression analysis. Paired T-test is used to assess the X-ray energy effects on TBS. Multiple linear regressions (backward) were used to evaluate relationships between TBS and 3D microarchitecture parameters using a bootstrap process. BV/TV of the sample ranged from 18.5 to 37.6% with an average value at 28.8%. Correlations' analysis showedthat TBSwere strongly correlatedwith ConnD(0.856≤r≤0.862; p<0.001),with TbN (0.805≤r≤0.810; p<0.001) and negatively with TbSp (−0.714≤r≤−0.726; p<0.001), regardless X-ray energy. Results show that lower TBS values are related to "degraded" microarchitecture, with low ConnD, low TbN and a high TbSp. The opposite is also true. X-ray energy has no effect onTBS neither on the correlations betweenTBS and the 3Dmicroarchitecture parameters. In this study, we demonstrated that TBS was significantly correlated with 3D microarchitecture parameters ConnD and TbN, and negatively with TbSp, no matter what X-ray energy has been used. This article is part of a Special Issue entitled ECTS 2011. Disclosure of interest: None declared.

JULIDE: a software tool for 3D reconstruction and statistical analysis of autoradiographic mouse brain sections.

In this article we introduce JULIDE, a software toolkit developed to perform the 3D reconstruction, intensity normalization, volume standardization by 3D image registration and voxel-wise statistical analysis of autoradiographs of mouse brain sections. This software tool has been developed in the open-source ITK software framework and is freely available under a GPL license. The article presents the complete image processing chain from raw data acquisition to 3D statistical group analysis. Results of the group comparison in the context of a study on spatial learning are shown as an illustration of the data that can be obtained with this tool.

Coronary MR imaging using free-breathing 3D steady-state free precession with radial k-space sampling.

PURPOSE: To investigate the potential of free-breathing 3D steady-state free precession (SSFP) imaging with radial k-space sampling for coronary MR-angiography (MRA), coronary projection MR-angiography and coronary vessel wall imaging. MATERIALS AND METHODS: A navigator-gated free-breathing T2-prepared 3D SSFP sequence (TR = 6.1 ms, TE = 3.0 ms, flip angle = 120 degrees, field-of-view = 360 mm(2)) with radial k-space sampling (384 radials) was implemented for coronary MRA. For projection coronary MRA, this sequence was combined with a 2D selective aortic spin tagging pulse. Coronary vessel wall imaging was performed using a high-resolution inversion-recovery black-blood 3D radial SSFP sequence (384 radials, TR = 5.3 ms, TE = 2.7 ms, flip angle = 55 degrees, reconstructed resolution 0.35 x 0.35 x 1.2 mm(3)) and a local re-inversion pulse. Six healthy volunteers (two for each sequence) were investigated. Motion artifact level was assessed by two radiologists. Results: In coronary MRA, the coronary lumen was displayed with a high signal and high contrast to the surrounding lumen. Projection coronary MRA demonstrated selective visualization of the coronary lumen while surrounding tissue was almost completely suppressed. In coronary vessel wall imaging, the vessel wall was displayed with a high signal when compared to the blood pool and the surrounding tissue. No visible motion artifacts were seen. Conclusion: 3D radial SSFP imaging enables coronary MRA, coronary projection MRA and coronary vessel wall imaging with a low motion artifact level.

Segmentation and grid generation for numerical simulations of vad connections with patient-specific data

Objectives: We are interested in the numerical simulation of the anastomotic region comprised between outflow canula of LVAD and the aorta. Segmenta¬tion, geometry reconstruction and grid generation from patient-specific data remain an issue because of the variable quality of DICOM images, in particular CT-scan (e.g. metallic noise of the device, non-aortic contrast phase). We pro¬pose a general framework to overcome this problem and create suitable grids for numerical simulations.Methods: Preliminary treatment of images is performed by reducing the level window and enhancing the contrast of the greyscale image using contrast-limited adaptive histogram equalization. A gradient anisotropic diffusion filter is applied to reduce the noise. Then, watershed segmentation algorithms and mathematical morphology filters allow reconstructing the patient geometry. This is done using the InsightToolKit library (www.itk.org). Finally the Vascular Model¬ing ToolKit (www.vmtk.org) and gmsh (www.geuz.org/gmsh) are used to create the meshes for the fluid (blood) and structure (arterial wall, outflow canula) and to a priori identify the boundary layers. The method is tested on five different patients with left ventricular assistance and who underwent a CT-scan exam.Results: This method produced good results in four patients. The anastomosis area is recovered and the generated grids are suitable for numerical simulations. In one patient the method failed to produce a good segmentation because of the small dimension of the aortic arch with respect to the image resolution.Conclusions: The described framework allows the use of data that could not be otherwise segmented by standard automatic segmentation tools. In particular the computational grids that have been generated are suitable for simulations that take into account fluid-structure interactions. Finally the presented method features a good reproducibility and fast application.

A multidimensional segmentation evaluation for medical image data.

Evaluation of segmentation methods is a crucial aspect in image processing, especially in the medical imaging field, where small differences between segmented regions in the anatomy can be of paramount importance. Usually, segmentation evaluation is based on a measure that depends on the number of segmented voxels inside and outside of some reference regions that are called gold standards. Although some other measures have been also used, in this work we propose a set of new similarity measures, based on different features, such as the location and intensity values of the misclassified voxels, and the connectivity and the boundaries of the segmented data. Using the multidimensional information provided by these measures, we propose a new evaluation method whose results are visualized applying a Principal Component Analysis of the data, obtaining a simplified graphical method to compare different segmentation results. We have carried out an intensive study using several classic segmentation methods applied to a set of MRI simulated data of the brain with several noise and RF inhomogeneity levels, and also to real data, showing that the new measures proposed here and the results that we have obtained from the multidimensional evaluation, improve the robustness of the evaluation and provides better understanding about the difference between segmentation methods.