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.

Measuring an effect size from dichotomized data : contrasted results whether using a correlation or an odds-ratio

It is well known that dichotomizing continuous data has the effect to decrease statistical power when the goal is to test for a statistical association between two variables. Modern researchers however are focusing not only on statistical significance but also on an estimation of the "effect size" (i.e., the strength of association between the variables) to judge whether a significant association is also clinically relevant. In this article, we are interested in the consequences of dichotomizing continuous data on the value of an effect size in some classical settings. It turns out that the conclusions will not be the same whether using a correlation or an odds ratio to summarize the strength of association between the variables: Whereas the value of a correlation is typically decreased by a factor pi/2 after each dichotomization, the value of an odds ratio is at the same time raised to the power 2. From a descriptive statistical point of view, it is thus not clear whether dichotomizing continuous data leads to a decrease or to an increase in the effect size, as illustrated using a data set to investigate the relationship between motor and intellectual functions in children and adolescents

The rationale of retinal endovascular fibrinolysis in the treatment of retinal vein occlusion : from experimental data to clinical application.

PURPOSE: : We describe a retinal endovascular fibrinolysis technique to directly reperfuse experimentally occluded retinal veins using a simple micropipette. METHODS: : Retinal vein occlusion was photochemically induced in 12 eyes of 12 minipigs: after intravenous injection of 10% fluorescein (1-mL bolus), the targeted retinal vein segment was exposed to thrombin (50 units) and to Argon laser (100-200 mW) through a pars plana approach. A beveled micropipette with a 30-μm-diameter sharp edge was used for micropuncture of the occluded vein and endovascular microinjection of tissue plasminogen activator (50 μg/mL) in 11 eyes. In one control eye, balanced salt solution was injected. The lesion site was examined histologically. RESULTS: : Retinal vein occlusion was achieved in all cases. Endovascular microinjection of tissue plasminogen activator or balanced salt solution led to reperfusion of the occluded retinal vein in all cases. Indicative of successful reperfusion were the following: continuous endovascular flow, unaffected collateral circulation, no optic disk ischemia, and no venous wall bleeding. However, balanced salt solution injection was accompanied by thrombus formation at the punctured site, whereas no thrombus was observed with tissue plasminogen activator injection. CONCLUSION: : Retinal endovascular fibrinolysis constitutes an efficient method of micropuncture and reperfusion of an experimentally occluded retinal vein. Thrombus formation at the punctured site can be prevented by injection of tissue plasminogen activator.

How quickly should we titrate antihypertensive medication? Systematic review modelling blood pressure response from trial data.

Context There are no evidence syntheses available to guide clinicians on when to titrate antihypertensive medication after initiation. Objective To model the blood pressure (BP) response after initiating antihypertensive medication. Data sources electronic databases including Medline, Embase, Cochrane Register and reference lists up to December 2009. Study selection Trials that initiated antihypertensive medication as single therapy in hypertensive patients who were either drug naive or had a placebo washout from previous drugs. Data extraction Office BP measurements at a minimum of two weekly intervals for a minimum of 4 weeks. An asymptotic approach model of BP response was assumed and non-linear mixed effects modelling used to calculate model parameters. Results and conclusions Eighteen trials that recruited 4168 patients met inclusion criteria. The time to reach 50% of the maximum estimated BP lowering effect was 1 week (systolic 0.91 weeks, 95% CI 0.74 to 1.10; diastolic 0.95, 0.75 to 1.15). Models incorporating drug class as a source of variability did not improve fit of the data. Incorporating the presence of a titration schedule improved model fit for both systolic and diastolic pressure. Titration increased both the predicted maximum effect and the time taken to reach 50% of the maximum (systolic 1.2 vs 0.7 weeks; diastolic 1.4 vs 0.7 weeks). Conclusions Estimates of the maximum efficacy of antihypertensive agents can be made early after starting therapy. This knowledge will guide clinicians in deciding when a newly started antihypertensive agent is likely to be effective or not at controlling BP.

Female-biased dispersal in the monogamous mammal Crocidura russula: evidence from field data and microsatellite patterns.

We investigated dispersal patterns in the monogamous Crocidura russula, based both on direct field observations (mark-recapture data) and on genetic analyses (microsatellite loci). Natal dispersal was found to be low. Most juveniles settled within their natal territory or one immediately adjacent. Migration rate was estimated to two individuals per year and per population. The correlation between genetic and geographical distances over a 16 km transect implies that migration occurs over short ranges. Natal dispersal was restricted to first-litter juveniles weaned in early May; this result suggests a direct dependence of dispersal on reproductive opportunities. Natal dispersal was highly female biased, a pattern unusual among mammals. Its association with monogamy provides support for the resource-competition model of dispersal. Our results demonstrate that a state-biased dispersal can be directly inferred from microsatellite genotype distributions, which opens new perspectives for empirical studies in this area.

The influence of supplementary health insurance on switching behaviour: evidence from Swiss data.

This paper focuses on the switching behaviour of enrolees in the Swiss basic health insurance system. Even though the new Federal Law on Social Health Insurance (LAMal) was implemented in 1996 to promote competition among health insurers in basic insurance, there is limited evidence of premium convergence within cantons. This indicates that competition has not been effective so far, and reveals some inertia among consumers who seem reluctant to switch to less expensive funds. We investigate one possible barrier to switching behaviour, namely the influence of supplementary insurance. We use survey data on health plan choice (a sample of 1943 individuals whose switching behaviours were observed between 1997 and 2000) as well as administrative data relative to all insurance companies that operated in the 26 Swiss cantons between 1996 and 2005. The decision to switch and the decision to subscribe to a supplementary contract are jointly estimated.Our findings show that holding a supplementary insurance contract substantially decreases the propensity to switch. However, there is no negative impact of supplementary insurance on switching when the individual assesses his/her health as 'very good'. Our results give empirical support to one possible mechanism through which supplementary insurance might influence switching decisions: given that subscribing to basic and supplementary contracts with two different insurers may induce some administrative costs for the subscriber, holding supplementary insurance acts as a barrier to switch if customers who consider themselves 'bad risks' also believe that insurers reject applications for supplementary insurance on these grounds. In comparison with previous research, our main contribution is to offer a possible explanation for consumer inertia. Our analysis illustrates how consumer choice for one's basic health plan interacts with the decision to subscribe to supplementary insurance.

Novel methods improve prediction of species' distributions from occurrence data

Prediction of species' distributions is central to diverse applications in ecology, evolution and conservation science. There is increasing electronic access to vast sets of occurrence records in museums and herbaria, yet little effective guidance on how best to use this information in the context of numerous approaches for modelling distributions. To meet this need, we compared 16 modelling methods over 226 species from 6 regions of the world, creating the most comprehensive set of model comparisons to date. We used presence-only data to fit models, and independent presence-absence data to evaluate the predictions. Along with well-established modelling methods such as generalised additive models and GARP and BIOCLIM, we explored methods that either have been developed recently or have rarely been applied to modelling species' distributions. These include machine-learning methods and community models, both of which have features that may make them particularly well suited to noisy or sparse information, as is typical of species' occurrence data. Presence-only data were effective for modelling species' distributions for many species and regions. The novel methods consistently outperformed more established methods. The results of our analysis are promising for the use of data from museums and herbaria, especially as methods suited to the noise inherent in such data improve.

Comparison of Strategies to Detect Epistasis from eQTL Data.

Genome-wide association studies have been instrumental in identifying genetic variants associated with complex traits such as human disease or gene expression phenotypes. It has been proposed that extending existing analysis methods by considering interactions between pairs of loci may uncover additional genetic effects. However, the large number of possible two-marker tests presents significant computational and statistical challenges. Although several strategies to detect epistasis effects have been proposed and tested for specific phenotypes, so far there has been no systematic attempt to compare their performance using real data. We made use of thousands of gene expression traits from linkage and eQTL studies, to compare the performance of different strategies. We found that using information from marginal associations between markers and phenotypes to detect epistatic effects yielded a lower false discovery rate (FDR) than a strategy solely using biological annotation in yeast, whereas results from human data were inconclusive. For future studies whose aim is to discover epistatic effects, we recommend incorporating information about marginal associations between SNPs and phenotypes instead of relying solely on biological annotation. Improved methods to discover epistatic effects will result in a more complete understanding of complex genetic effects.

Reconstruction Of Past Climate Conditions Over Central Europe From Groundwater Data

Here we present a 30 000 years low-resolution climate record reconstructed from groundwater data. The investigated site is located in the Bohemian Cretaceous Basin, in the corridor between the Scandinavian ice sheet and the Alpine ice field. Noble gas temperatures (NGT), obtained from groundwater data, preserved multicentennial temperature variability and indicated a cooling of at least 5-7 °C during the last glacial maximum (LGM). This is further confirmed by the depleted δ18O and δ2H values at the LGM. High excess air (ΔNe) at the end of the Pleistocene is possibly related to abrupt changes in recharge dynamics due to progression and retreat of ice covers and permafrost. These results agree with the fact that during the LGM permafrost and small glaciers developed in the inner valleys of the Giant Mountains (located in the watershed of the aquifers). A temporal decrease of deuterium excess from the pre-industrial Holocene to present days is linked to an increase of the air temperatures, and probably also to an increase of water pressure at the source region of precipitation over the past few hundred years

Formica paralugubris (Hymenoptera, Formicidae) in the Italian Alps from new data and old data revisited.

We provide evidence that Formica paralugubris Seifert, 1996, a species of wood and recently described from Switzerland, is present in the Italian Alps. Until 1996, this species was confounded with F.lugubris Zetterstedt, 1838. We examine the wood and collection deposited at the University of Pavia (Italy) and collect new samples within the Italian Alps. Formica paralugubris seems to be more abundant than F.lugubris. Moreover, both species are found in sympatry in some localities

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.

Histology-driven data mining of lipid signatures from multiple imaging mass spectrometry analyses: application to human colorectal cancer liver metastasis biopsies.

Imaging mass spectrometry (IMS) represents an innovative tool in the cancer research pipeline, which is increasingly being used in clinical and pharmaceutical applications. The unique properties of the technique, especially the amount of data generated, make the handling of data from multiple IMS acquisitions challenging. This work presents a histology-driven IMS approach aiming to identify discriminant lipid signatures from the simultaneous mining of IMS data sets from multiple samples. The feasibility of the developed workflow is evaluated on a set of three human colorectal cancer liver metastasis (CRCLM) tissue sections. Lipid IMS on tissue sections was performed using MALDI-TOF/TOF MS in both negative and positive ionization modes after 1,5-diaminonaphthalene matrix deposition by sublimation. The combination of both positive and negative acquisition results was performed during data mining to simplify the process and interrogate a larger lipidome into a single analysis. To reduce the complexity of the IMS data sets, a sub data set was generated by randomly selecting a fixed number of spectra from a histologically defined region of interest, resulting in a 10-fold data reduction. Principal component analysis confirmed that the molecular selectivity of the regions of interest is maintained after data reduction. Partial least-squares and heat map analyses demonstrated a selective signature of the CRCLM, revealing lipids that are significantly up- and down-regulated in the tumor region. This comprehensive approach is thus of interest for defining disease signatures directly from IMS data sets by the use of combinatory data mining, opening novel routes of investigation for addressing the demands of the clinical setting.

Hyperhidrosis: a new and often early symptom in Fabry disease. International experience and data from the Fabry Outcome Survey.

Hypohidrosis is a classic feature of Fabry disease; in contrast, hyperhidrosis has only been rarely described. The aim of the study is to characterise the baseline descriptive data on hyperhidrosis (frequency, age at onset, sex ratio and outcome with and without enzyme replacement therapy) in hemizygous male and heterozygous female patients with Fabry disease. We describe case histories of five patients with Fabry disease and hyperhidrosis seen at three different centres. We have also analysed a cohort of 21 paediatric patients in the UK and a large European cohort of patients enrolled in the Fabry Outcome Survey (FOS). Five patients (three female, two male) with hyperhidrosis were originally identified, although each had additional symptoms related to Fabry disease. The age at onset of hyperhidrosis was less than 18 years in four cases. In the cohort of 21 paediatric patients (12 female, nine male), one female had hyperhidrosis; the age at onset of this symptom was 11 years. In the FOS cohort, 66 of 714 patients with Fabry disease had hyperhidrosis (44 of 369 females, 11.9%; 22 of 345 males, 6.4%). The female predominance was observed in seven of nine countries from which data were analysed. Hyperhidrosis is an increasingly recognised feature of the Fabry disease phenotype. It is more prevalent in females than in males and often appears in childhood or adolescence. The efficacy of enzyme replacement therapy on this recently recognised symptom should be assessed.

Are invasive marsh frogs (Pelophylax ridibundus) replacing the native P. lessonae /P. esculentus hybridogenetic complex in Western Europe? Genetic evidence from a field study.

The water-frog L-E system, widespread in Western Europe, comprises the pool frog Pelophylax lessonae and the hybridogenetic edible frog P. esculentus, which originated from hybridization between pool frogs and marsh frogs (P. ridibundus). In P. esculentus, the lessonae (L) genome is eliminated during meiosis and has to be gained anew each generation from a P. lessonae partner, while the ridibundus (R') genome is transmitted clonally. It therefore accumulates deleterious mutations, so that R'R' offspring from P. esculentus x P. esculentus crosses are normally unviable. This system is now threatened by invasive P. ridibundus (RR) imported from Eastern Europe and the Balkans. We investigated the genetic interactions between invasive marsh frogs and native water frogs in a Swiss wetland area, and used genetic data collected in the field to validate several components of a recently postulated mechanism of species replacement. We identified neo-ridibundus individuals derived from crosses between invasive ridibundus and native esculentus, as well as newly arisen hybridogenetic esculentus lineages stemming from crosses between invasive ridibundus (RR) and native lessonae (LL). As their ridibundus genomes are likely to carry less deleterious mutations, such lineages are expected to produce viable ridibundus offspring, contributing to species replacement. However, such crosses with invasive ridibundus only occurred at a limited scale; moreover, RR x LL crosses did not induce any introgression from the ridibundus to the lessonae genome. We did not find any ridibundus stemming from crosses between ancient esculentus lineages. Despite several decades of presence on the site, introduced ridibundus individuals only represent 15% of sampled frogs, and their spatial distribution seems shaped by specific ecological requirements rather than history of colonization. We therefore expect the three taxa to coexist stably in this area.