Multilayer perceptron with local constraint as an emerging method in spatial data analysis

The use of Geographic Information Systems has revolutionalized the handling and the visualization of geo-referenced data and has underlined the critic role of spatial analysis. The usual tools for such a purpose are geostatistics which are widely used in Earth science. Geostatistics are based upon several hypothesis which are not always verified in practice. On the other hand, Artificial Neural Network (ANN) a priori can be used without special assumptions and are known to be flexible. This paper proposes to discuss the application of ANN in the case of the interpolation of a geo-referenced variable.

IQRray, a new method for Affymetrix microarray quality control, and the homologous organ conservation score, a new benchmark method for quality control metrics.

MOTIVATION: Microarray results accumulated in public repositories are widely reused in meta-analytical studies and secondary databases. The quality of the data obtained with this technology varies from experiment to experiment, and an efficient method for quality assessment is necessary to ensure their reliability. RESULTS: The lack of a good benchmark has hampered evaluation of existing methods for quality control. In this study, we propose a new independent quality metric that is based on evolutionary conservation of expression profiles. We show, using 11 large organ-specific datasets, that IQRray, a new quality metrics developed by us, exhibits the highest correlation with this reference metric, among 14 metrics tested. IQRray outperforms other methods in identification of poor quality arrays in datasets composed of arrays from many independent experiments. In contrast, the performance of methods designed for detecting outliers in a single experiment like Normalized Unscaled Standard Error and Relative Log Expression was low because of the inability of these methods to detect datasets containing only low-quality arrays and because the scores cannot be directly compared between experiments. AVAILABILITY AND IMPLEMENTATION: The R implementation of IQRray is available at: ftp://lausanne.isb-sib.ch/pub/databases/Bgee/general/IQRray.R. CONTACT: Marta.Rosikiewicz@unil.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Quantifying rates of landscape evolution and tectonic processes by thermochronology and numerical modeling of crustal heat transport using PECUBE

PECUBE is a three-dimensional thermal-kinematic code capable of solving the heat production-diffusion-advection equation under a temporally varying surface boundary condition. It was initially developed to assess the effects of time-varying surface topography (relief) on low-temperature thermochronological datasets. Thermochronometric ages are predicted by tracking the time-temperature histories of rock-particles ending up at the surface and by combining these with various age-prediction models. In the decade since its inception, the PECUBE code has been under continuous development as its use became wider and addressed different tectonic-geomorphic problems. This paper describes several major recent improvements in the code, including its integration with an inverse-modeling package based on the Neighborhood Algorithm, the incorporation of fault-controlled kinematics, several different ways to address topographic and drainage change through time, the ability to predict subsurface (tunnel or borehole) data, prediction of detrital thermochronology data and a method to compare these with observations, and the coupling with landscape-evolution (or surface-process) models. Each new development is described together with one or several applications, so that the reader and potential user can clearly assess and make use of the capabilities of PECUBE. We end with describing some developments that are currently underway or should take place in the foreseeable future. (C) 2012 Elsevier B.V. All rights reserved.

Multiplex ultra-performance liquid chromatography-tandem mass spectrometry method for simultaneous quantification in human plasma of fluconazole, itraconazole, hydroxyitraconazole, posaconazole, voriconazole, voriconazole-N-oxide, anidulafungin, and caspofungin.

Therapeutic drug monitoring (TDM) may contribute to optimizing the efficacy and safety of antifungal therapy because of the large variability in drug pharmacokinetics. Rapid, sensitive, and selective laboratory methods are needed for efficient TDM. Quantification of several antifungals in a single analytical run may best fulfill these requirements. We therefore developed a multiplex ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) method requiring 100 μl of plasma for simultaneous quantification within 7 min of fluconazole, itraconazole, hydroxyitraconazole, posaconazole, voriconazole, voriconazole-N-oxide, caspofungin, and anidulafungin. Protein precipitation with acetonitrile was used in a single extraction procedure for eight analytes. After reverse-phase chromatographic separation, antifungals were quantified by electrospray ionization-triple-quadrupole mass spectrometry by selected reaction monitoring detection using the positive mode. Deuterated isotopic compounds of azole antifungals were used as internal standards. The method was validated based on FDA recommendations, including assessment of extraction yields, matrix effect variability (<9.2%), and analytical recovery (80.1 to 107%). The method is sensitive (lower limits of azole quantification, 0.01 to 0.1 μg/ml; those of echinocandin quantification, 0.06 to 0.1 μg/ml), accurate (intra- and interassay biases of -9.9 to +5% and -4.0 to +8.8%, respectively), and precise (intra- and interassay coefficients of variation of 1.2 to 11.1% and 1.2 to 8.9%, respectively) over clinical concentration ranges (upper limits of quantification, 5 to 50 μg/ml). Thus, we developed a simple, rapid, and robust multiplex UPLC-MS/MS assay for simultaneous quantification of plasma concentrations of six antifungals and two metabolites. This offers, by optimized and cost-effective lab resource utilization, an efficient tool for daily routine TDM aimed at maximizing the real-time efficacy and safety of different recommended single-drug antifungal regimens and combination salvage therapies, as well as a tool for clinical research.

An optimized method for establishing high purity murine CD8+ T cell cultures.

Establishing CD8(+) T cell cultures has been empirical and the published methods have been largely individual laboratory based. In this study, we optimized culturing conditions and show that IL-2 concentration is the most critical factor for the success of establishing CD8(+) T cell cultures. High IL-2 concentration encouraged T cells to non-specifically proliferate, express a B cell marker, B220, and undergo apoptosis. These cells also lose typical irregular T cell morphology and are incapable of sustaining long-term cultures. Using tetramer and intracellular cytokine assessments, we further demonstrated that many antigen-specific T cells have been rendered nonfunctional when expanded under high IL-2 concentration. When IL-2 is used in the correct range, B220-mediated cell depletion greatly enhanced the success rate of such T cell cultures.

A rapid, reliable method of evaluating growth and viability of intraerythrocytic protozoan hemoparasites using fluorescence flow cytometry

Fluorescence flow cytometry was employed to assess the potential of a vital dye, hydroethiedine, for use in the detection and monitoring of the viability of hemoparasites in infected erythrocytes, using Babesia bovis as a model parasite. The studies demonstrated that hydroethidine is taken up by B. bovis and metabolically converted to the DNA binding fluorochrone, ethidium. Following uptake of the dye, erythrocytes contamine viable parasites were readily distinguished and quantitated. Timed studies with the parasiticidal drug, Ganaseg, showed that it is possible to use the fluorochrome assay to monitor the effects of the drug on the rate of replication and viability of B. bovis in culture. The assay provides a rapid method for evaluation of the in vitro effect of drugs on hemoparasites and for analysis of the effect of various components of the immune response, such as lymphokines, monocyte products, antibodies, and effector cells (T, NK, LAK, ADCC) on the growth and viability of intraerythrocytic parasites.

Combining food web and species distribution models for improved community projections

The ability to model biodiversity patterns is of prime importance in this era of severe environmental crisis. Species assemblage along environmental gradient is subject to the interplay of biotic interactions in complement to abiotic environmental filtering. Accounting for complex biotic interactions for a wide array of species remains so far challenging. Here, we propose to use food web models that can infer the potential interaction links between species as a constraint in species distribution models. Using a plant-herbivore (butterfly) interaction dataset, we demonstrate that this combined approach is able to improve both species distribution and community forecasts. Most importantly, this combined approach is very useful in rendering models of more generalist species that have multiple potential interaction links, where gap in the literature may be recurrent. Our combined approach points a promising direction forward to model the spatial variation of entire species interaction networks. Our work has implications for studies of range shifting species and invasive species biology where it may be unknown how a given biota might interact with a potential invader or in future climate.

Colonoscopy is the preferred colorectal cancer screening method in a population-based program.

BACKGROUND AND STUDY AIMS: Various screening methods for colorectal cancer (CRC) are promoted by professional societies; however, few data are available about the factors that determine patient participation in screening, which is crucial to the success of population-based programs. This study aimed (i) to identify factors that determine acceptance of screening and preference of screening method, and (ii) to evaluate procedure success, detection of colorectal neoplasia, and patient satisfaction with screening colonoscopy. PATIENTS AND METHODS: Following a public awareness campaign, the population aged 50 - 80 years was offered CRC screening in the form of annual fecal occult blood tests, flexible sigmoidoscopy, a combination of both, or colonoscopy. RESULTS: 2731 asymptomatic persons (12.0 % of the target population) registered with and were eligible to take part in the screening program. Access to information and a positive attitude to screening were major determinants of participation. Colonoscopy was the method preferred by 74.8 % of participants. Advanced colorectal neoplasia was present in 8.5 %; its prevalence was higher in males and increased with age. Significant complications occurred in 0.5 % of those undergoing colonoscopy and were associated with polypectomy or sedation. Most patients were satisfied with colonoscopy and over 90 % would choose it again for CRC screening. CONCLUSIONS: In this population-based study, only a small proportion of the target population underwent CRC screening despite an extensive information campaign. Colonoscopy was the preferred method and was safe. The determinants of participation in screening and preference of screening method, together with the distribution of colorectal neoplasia in different demographic categories, provide a rationale for improving screening procedures.

A mixed finite element method for nonlinear diffusion equations

We propose a mixed finite element method for a class of nonlinear diffusion equations, which is based on their interpretation as gradient flows in optimal transportation metrics. We introduce an appropriate linearization of the optimal transport problem, which leads to a mixed symmetric formulation. This formulation preserves the maximum principle in case of the semi-discrete scheme as well as the fully discrete scheme for a certain class of problems. In addition solutions of the mixed formulation maintain exponential convergence in the relative entropy towards the steady state in case of a nonlinear Fokker-Planck equation with uniformly convex potential. We demonstrate the behavior of the proposed scheme with 2D simulations of the porous medium equations and blow-up questions in the Patlak-Keller-Segel model.

Segmentation of brain structures in presence of a space-occupying lesion.

Brain deformations induced by space-occupying lesions may result in unpredictable position and shape of functionally important brain structures. The aim of this study is to propose a method for segmentation of brain structures by deformation of a segmented brain atlas in presence of a space-occupying lesion. Our approach is based on an a priori model of lesion growth (MLG) that assumes radial expansion from a seeding point and involves three steps: first, an affine registration bringing the atlas and the patient into global correspondence; then, the seeding of a synthetic tumor into the brain atlas providing a template for the lesion; finally, the deformation of the seeded atlas, combining a method derived from optical flow principles and a model of lesion growth. The method was applied on two meningiomas inducing a pure displacement of the underlying brain structures, and segmentation accuracy of ventricles and basal ganglia was assessed. Results show that the segmented structures were consistent with the patient's anatomy and that the deformation accuracy of surrounding brain structures was highly dependent on the accurate placement of the tumor seeding point. Further improvements of the method will optimize the segmentation accuracy. Visualization of brain structures provides useful information for therapeutic consideration of space-occupying lesions, including surgical, radiosurgical, and radiotherapeutic planning, in order to increase treatment efficiency and prevent neurological damage.

In silico prediction of human carboxylesterase-1 (hCES1) metabolism combining docking analyses and MD simulations.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted in developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas little has been done to predict the hydrolytic activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES1. The study involves both docking analyses of known substrates to develop predictive models, and molecular dynamics (MD) simulations to reveal the in situ behavior of substrates and products, with particular attention being paid to the influence of their ionization state. The results emphasize some crucial properties of the hCES1 catalytic cavity, confirming that as a trend with several exceptions, hCES1 prefers substrates with relatively smaller and somewhat polar alkyl/aryl groups and larger hydrophobic acyl moieties. The docking results underline the usefulness of the hydrophobic interaction score proposed here, which allows a robust prediction of hCES1 catalysis, while the MD simulations show the different behavior of substrates and products in the enzyme cavity, suggesting in particular that basic substrates interact with the enzyme in their unprotonated form.