Multiple Kernel Learning of Environmental Data. Case study: Analysis and Mapping of Wind Fields

The paper presents the Multiple Kernel Learning (MKL) approach as a modelling and data exploratory tool and applies it to the problem of wind speed mapping. Support Vector Regression (SVR) is used to predict spatial variations of the mean wind speed from terrain features (slopes, terrain curvature, directional derivatives) generated at different spatial scales. Multiple Kernel Learning is applied to learn kernels for individual features and thematic feature subsets, both in the context of feature selection and optimal parameters determination. An empirical study on real-life data confirms the usefulness of MKL as a tool that enhances the interpretability of data-driven models.

Magnetic resonance imaging correlates of first-episode psychosis in young adult male patients: combined analysis of grey and white matter.

Background: Several patterns of grey and white matter changes have been separately described in young adults with first-episode psychosis. Concomitant investigation of grey and white matter densities in patients with first-episode psychosis without other psychiatric comorbidities that include all relevant imaging markers could provide clues to the neurodevelopmental hypothesis in schizophrenia. Methods: We recruited patients with first-episode psychosis diagnosed according to the DSM-IV-TR and matched controls. All participants underwent magnetic resonance imaging (MRI). Voxel-based morphometry (VBM) analysis and mean diffusivity voxel-based analysis (VBA) were used for grey matter data. Fractional anisotropy and axial, radial and mean diffusivity were analyzed using tract-based spatial statistics (TBSS) for white matter data. Results: We included 15 patients and 16 controls. The mean diffusivity VBA showed significantly greater mean diffusivity in the first-episode psychosis than in the control group in the lingual gyrus bilaterally, the occipital fusiform gyrus bilaterally, the right lateral occipital gyrus and the right inferior temporal gyrus. Moreover, the TBSS analysis revealed a lower fractional anisotropy in the first-episode psychosis than in the control group in the genu of the corpus callosum, minor forceps, corticospinal tract, right superior longitudinal fasciculus, left middle cerebellar peduncle, left inferior longitudinal fasciculus and the posterior part of the fronto-occipital fasciculus. This analysis also revealed greater radial diffusivity in the first-episode psychosis than in the control group in the right corticospinal tract, right superior longitudinal fasciculus and left middle cerebellar peduncle. Limitations: The modest sample size and the absence of women in our series could limit the impact of our results. Conclusion: Our results highlight the structural vulnerability of grey matter in posterior areas of the brain among young adult male patients with first-episode psychosis. Moreover, the concomitant greater radial diffusivity within several regions already revealed by the fractional anisotropy analysis supports the idea of a late myelination in patients with first-episode psychosis.

Individual prediction of cognitive decline in mild cognitive impairment using support vector machine-based analysis of diffusion tensor imaging data.

Although cross-sectional diffusion tensor imaging (DTI) studies revealed significant white matter changes in mild cognitive impairment (MCI), the utility of this technique in predicting further cognitive decline is debated. Thirty-five healthy controls (HC) and 67 MCI subjects with DTI baseline data were neuropsychologically assessed at one year. Among them, there were 40 stable (sMCI; 9 single domain amnestic, 7 single domain frontal, 24 multiple domain) and 27 were progressive (pMCI; 7 single domain amnestic, 4 single domain frontal, 16 multiple domain). Fractional anisotropy (FA) and longitudinal, radial, and mean diffusivity were measured using Tract-Based Spatial Statistics. Statistics included group comparisons and individual classification of MCI cases using support vector machines (SVM). FA was significantly higher in HC compared to MCI in a distributed network including the ventral part of the corpus callosum, right temporal and frontal pathways. There were no significant group-level differences between sMCI versus pMCI or between MCI subtypes after correction for multiple comparisons. However, SVM analysis allowed for an individual classification with accuracies up to 91.4% (HC versus MCI) and 98.4% (sMCI versus pMCI). When considering the MCI subgroups separately, the minimum SVM classification accuracy for stable versus progressive cognitive decline was 97.5% in the multiple domain MCI group. SVM analysis of DTI data provided highly accurate individual classification of stable versus progressive MCI regardless of MCI subtype, indicating that this method may become an easily applicable tool for early individual detection of MCI subjects evolving to dementia.

ConnectomeViewer - Multi-Modal Multi-Level Network Visualization and Analysis.

Introduction: The field of Connectomic research is growing rapidly, resulting from methodological advances in structural neuroimaging on many spatial scales. Especially progress in Diffusion MRI data acquisition and processing made available macroscopic structural connectivity maps in vivo through Connectome Mapping Pipelines (Hagmann et al, 2008) into so-called Connectomes (Hagmann 2005, Sporns et al, 2005). They exhibit both spatial and topological information that constrain functional imaging studies and are relevant in their interpretation. The need for a special-purpose software tool for both clinical researchers and neuroscientists to support investigations of such connectome data has grown. Methods: We developed the ConnectomeViewer, a powerful, extensible software tool for visualization and analysis in connectomic research. It uses the novel defined container-like Connectome File Format, specifying networks (GraphML), surfaces (Gifti), volumes (Nifti), track data (TrackVis) and metadata. Usage of Python as programming language allows it to by cross-platform and have access to a multitude of scientific libraries. Results: Using a flexible plugin architecture, it is possible to enhance functionality for specific purposes easily. Following features are already implemented: * Ready usage of libraries, e.g. for complex network analysis (NetworkX) and data plotting (Matplotlib). More brain connectivity measures will be implemented in a future release (Rubinov et al, 2009). * 3D View of networks with node positioning based on corresponding ROI surface patch. Other layouts possible. * Picking functionality to select nodes, select edges, get more node information (ConnectomeWiki), toggle surface representations * Interactive thresholding and modality selection of edge properties using filters * Arbitrary metadata can be stored for networks, thereby allowing e.g. group-based analysis or meta-analysis. * Python Shell for scripting. Application data is exposed and can be modified or used for further post-processing. * Visualization pipelines using filters and modules can be composed with Mayavi (Ramachandran et al, 2008). * Interface to TrackVis to visualize track data. Selected nodes are converted to ROIs for fiber filtering The Connectome Mapping Pipeline (Hagmann et al, 2008) processed 20 healthy subjects into an average Connectome dataset. The Figures show the ConnectomeViewer user interface using this dataset. Connections are shown that occur in all 20 subjects. The dataset is freely available from the homepage (connectomeviewer.org). Conclusions: The ConnectomeViewer is a cross-platform, open-source software tool that provides extensive visualization and analysis capabilities for connectomic research. It has a modular architecture, integrates relevant datatypes and is completely scriptable. Visit www.connectomics.org to get involved as user or developer.

Predictive analysis and mapping of indoor radon concentrations in Switzerland

It is estimated that around 230 people die each year due to radon (222Rn) exposure in Switzerland. 222Rn occurs mainly in closed environments like buildings and originates primarily from the subjacent ground. Therefore it depends strongly on geology and shows substantial regional variations. Correct identification of these regional variations would lead to substantial reduction of 222Rn exposure of the population based on appropriate construction of new and mitigation of already existing buildings. Prediction of indoor 222Rn concentrations (IRC) and identification of 222Rn prone areas is however difficult since IRC depend on a variety of different variables like building characteristics, meteorology, geology and anthropogenic factors. The present work aims at the development of predictive models and the understanding of IRC in Switzerland, taking into account a maximum of information in order to minimize the prediction uncertainty. The predictive maps will be used as a decision-support tool for 222Rn risk management. The construction of these models is based on different data-driven statistical methods, in combination with geographical information systems (GIS). In a first phase we performed univariate analysis of IRC for different variables, namely the detector type, building category, foundation, year of construction, the average outdoor temperature during measurement, altitude and lithology. All variables showed significant associations to IRC. Buildings constructed after 1900 showed significantly lower IRC compared to earlier constructions. We observed a further drop of IRC after 1970. In addition to that, we found an association of IRC with altitude. With regard to lithology, we observed the lowest IRC in sedimentary rocks (excluding carbonates) and sediments and the highest IRC in the Jura carbonates and igneous rock. The IRC data was systematically analyzed for potential bias due to spatially unbalanced sampling of measurements. In order to facilitate the modeling and the interpretation of the influence of geology on IRC, we developed an algorithm based on k-medoids clustering which permits to define coherent geological classes in terms of IRC. We performed a soil gas 222Rn concentration (SRC) measurement campaign in order to determine the predictive power of SRC with respect to IRC. We found that the use of SRC is limited for IRC prediction. The second part of the project was dedicated to predictive mapping of IRC using models which take into account the multidimensionality of the process of 222Rn entry into buildings. We used kernel regression and ensemble regression tree for this purpose. We could explain up to 33% of the variance of the log transformed IRC all over Switzerland. This is a good performance compared to former attempts of IRC modeling in Switzerland. As predictor variables we considered geographical coordinates, altitude, outdoor temperature, building type, foundation, year of construction and detector type. Ensemble regression trees like random forests allow to determine the role of each IRC predictor in a multidimensional setting. We found spatial information like geology, altitude and coordinates to have stronger influences on IRC than building related variables like foundation type, building type and year of construction. Based on kernel estimation we developed an approach to determine the local probability of IRC to exceed 300 Bq/m3. In addition to that we developed a confidence index in order to provide an estimate of uncertainty of the map. All methods allow an easy creation of tailor-made maps for different building characteristics. Our work is an essential step towards a 222Rn risk assessment which accounts at the same time for different architectural situations as well as geological and geographical conditions. For the communication of 222Rn hazard to the population we recommend to make use of the probability map based on kernel estimation. The communication of 222Rn hazard could for example be implemented via a web interface where the users specify the characteristics and coordinates of their home in order to obtain the probability to be above a given IRC with a corresponding index of confidence. Taking into account the health effects of 222Rn, our results have the potential to substantially improve the estimation of the effective dose from 222Rn delivered to the Swiss population.

Predicting current and future spatial community patterns of plant functional traits

Community-level patterns of functional traits relate to community assembly and ecosystem functioning. By modelling the changes of different indices describing such patterns - trait means, extremes and diversity in communities - as a function of abiotic gradients, we could understand their drivers and build projections of the impact of global change on the functional components of biodiversity. We used five plant functional traits (vegetative height, specific leaf area, leaf dry matter content, leaf nitrogen content and seed mass) and non-woody vegetation plots to model several indices depicting community-level patterns of functional traits from a set of abiotic environmental variables (topographic, climatic and edaphic) over contrasting environmental conditions in a mountainous landscape. We performed a variation partitioning analysis to assess the relative importance of these variables for predicting patterns of functional traits in communities, and projected the best models under several climate change scenarios to examine future potential changes in vegetation functional properties. Not all indices of trait patterns within communities could be modelled with the same level of accuracy: the models for mean and extreme values of functional traits provided substantially better predictive accuracy than the models calibrated for diversity indices. Topographic and climatic factors were more important predictors of functional trait patterns within communities than edaphic predictors. Overall, model projections forecast an increase in mean vegetation height and in mean specific leaf area following climate warming. This trend was important at mid elevation particularly between 1000 and 2000 m asl. With this study we showed that topographic, climatic and edaphic variables can successfully model descriptors of community-level patterns of plant functional traits such as mean and extreme trait values. However, which factors determine the diversity of functional traits in plant communities remains unclear and requires more investigations.

Sublimation of new matrix candidates for high spatial resolution imaging mass spectrometry of lipids: enhanced information in both positive and negative polarities after 1,5-diaminonapthalene deposition.

Matrix sublimation has demonstrated to be a powerful approach for high-resolution matrix-assisted laser desorption ionization (MALDI) imaging of lipids, providing very homogeneous solvent-free deposition. This work presents a comprehensive study aiming to evaluate current and novel matrix candidates for high spatial resolution MALDI imaging mass spectrometry of lipids from tissue section after deposition by sublimation. For this purpose, 12 matrices including 2,5-dihydroxybenzoic acid (DHB), sinapinic acid (SA), α-cyano-4-hydroxycinnamic acid (CHCA), 2,6-dihydroxyacetphenone (DHA), 2',4',6'-trihydroxyacetophenone (THAP), 3-hydroxypicolinic acid (3-HPA), 1,8-bis(dimethylamino)naphthalene (DMAN), 1,8,9-anthracentriol (DIT), 1,5-diaminonapthalene (DAN), p-nitroaniline (NIT), 9-aminoacridine (9-AA), and 2-mercaptobenzothiazole (MBT) were investigated for lipid detection efficiency in both positive and negative ionization modes, matrix interferences, and stability under vacuum. For the most relevant matrices, ion maps of the different lipid species were obtained from tissue sections at high spatial resolution and the detected peaks were characterized by matrix-assisted laser desorption ionization time-of-flight/time-of-flight (MALDI-TOF/TOF) mass spectrometry. First proposed for imaging mass spectrometry (IMS) after sublimation, DAN has demonstrated to be of high efficiency providing rich lipid signatures in both positive and negative polarities with high vacuum stability and sub-20 μm resolution capacity. Ion images from adult mouse brain were generated with a 10 μm scanning resolution. Furthermore, ion images from adult mouse brain and whole-body fish tissue sections were also acquired in both polarity modes from the same tissue section at 100 μm spatial resolution. Sublimation of DAN represents an interesting approach to improve information with respect to currently employed matrices providing a deeper analysis of the lipidome by IMS.

Multiway Array Decomposition Analysis of EEGs in Alzheimer’s Disease

Methods for the extraction of features from physiological datasets are growing needs as clinical investigations of Alzheimer’s disease (AD) in large and heterogeneous population increase. General tools allowing diagnostic regardless of recording sites, such as different hospitals, are essential and if combined to inexpensive non-invasive methods could critically improve mass screening of subjects with AD. In this study, we applied three state of the art multiway array decomposition (MAD) methods to extract features from electroencephalograms (EEGs) of AD patients obtained from multiple sites. In comparison to MAD, spectral-spatial average filter (SSFs) of control and AD subjects were used as well as a common blind source separation method, algorithm for multiple unknown signal extraction (AMUSE). We trained a feed-forward multilayer perceptron (MLP) to validate and optimize AD classification from two independent databases. Using a third EEG dataset, we demonstrated that features extracted from MAD outperformed features obtained from SSFs AMUSE in terms of root mean squared error (RMSE) and reaching up to 100% of accuracy in test condition. We propose that MAD maybe a useful tool to extract features for AD diagnosis offering great generalization across multi-site databases and opening doors to the discovery of new characterization of the disease.

Spatial Polarization of the Ecological Footprint distribution

The international allocation of natural resources is determined, not by any ethical or ecological criteria, but by the dominance of market mechanisms. From a core-periphery perspective, this allocation may even be driven by historically determined structural patterns, with a core group of countries whose consumption appropriates most available natural resources, and another group, having low natural resource consumption, which plays a peripheral role. This article consists of an empirical distributional analysis of natural resource consumption (as measured by Ecological Footprints) whose purpose is to assess the extent to which the distribution of consumption responds to polarization (as opposed to mere inequality). To assess this, we estimate and decompose different polarization indices for a balanced sample of 119 countries over the period 1961 to 2007. Our results points toward a polarized distribution which is consistent with a core-periphery framework.

Predictive analysis and mapping of indoor radon concentrations in a complex environment using kernel estimation: an application to Switzerland.

PURPOSE: The aim of this study was to develop models based on kernel regression and probability estimation in order to predict and map IRC in Switzerland by taking into account all of the following: architectural factors, spatial relationships between the measurements, as well as geological information. METHODS: We looked at about 240,000 IRC measurements carried out in about 150,000 houses. As predictor variables we included: building type, foundation type, year of construction, detector type, geographical coordinates, altitude, temperature and lithology into the kernel estimation models. We developed predictive maps as well as a map of the local probability to exceed 300 Bq/m(3). Additionally, we developed a map of a confidence index in order to estimate the reliability of the probability map. RESULTS: Our models were able to explain 28% of the variations of IRC data. All variables added information to the model. The model estimation revealed a bandwidth for each variable, making it possible to characterize the influence of each variable on the IRC estimation. Furthermore, we assessed the mapping characteristics of kernel estimation overall as well as by municipality. Overall, our model reproduces spatial IRC patterns which were already obtained earlier. On the municipal level, we could show that our model accounts well for IRC trends within municipal boundaries. Finally, we found that different building characteristics result in different IRC maps. Maps corresponding to detached houses with concrete foundations indicate systematically smaller IRC than maps corresponding to farms with earth foundation. CONCLUSIONS: IRC mapping based on kernel estimation is a powerful tool to predict and analyze IRC on a large-scale as well as on a local level. This approach enables to develop tailor-made maps for different architectural elements and measurement conditions and to account at the same time for geological information and spatial relations between IRC measurements.

Spatial and temporal regulation of " cdc16p " in the fission yeast " Schizosaccharomyces pombe "

Abstract: The fission yeast Schizosaccharomyces pombe has proven to be an excellent model system for the study of eukaryotic cell cycle control. S. pombe cells are rod-shaped and grow mainly by elongation at their tips. They divide by the means of a centrallyplaced division septum which provides two daughter cells of equal size. S. pombe cytokinesis begins at mitotic entry, when the division site is defined by formation of the contractile acto-myosin ring (CAR). Formation of the division septum is triggered at the end of mitosis by the spindle pole body (SPB) associated septation initiation network (SIN) proteins. SIN signalling requires activation of the GTPase spg1p, whose nucleotide status is regulated by the bipartite GAP byr4pcdc16p. Removal of cdc16p from the SPB during early mitosis is thought to allow priming of the SIN by association of cdc7p with both SPBs. During anaphase cdc7p is retained on the new SPB, which also recruits the kinase sid1 p and cdc14p, while the old SP8 reassembles the byr4-cdc16p GAP and is presumed not to signal; SPB asymmetry persists throughout anaphase. The trigger for inactivation of SIN signalling at the new SPB is unknown. This study has concentrated upon cdc16p. We have undertaken the analysis of the localisation of cdc16p using time-lapse microscopy. We have observed that the localisation of cdc16p is regulated at different transitions. We have shown that cdc16p is removed from the SPB prior to the onset of spindle formation and that it reappears asymmetrically at the beginning of anaphase B. We have also demonstrated that the resetting of the SIN at the new SPB is linked to completion of CAR contraction and septum formation. We propose the existence of a mechanism that monitors cytokinesis and that couples the activity of the SI N with the presence of the CAR. During the biochemical characterization of cdc16p, We have found that it is an unstable protein and that it is subjected to polyubiquitination by the SCF and proteasomal degradation. Together, these observations help to shed new light upon the mechanisms by which cytokinesis is regulated in S. pombe. Résumé: La levure Schizosaccharomyces pombe est un excellent organisme modèle pour l'étude du cycle cellulaire eucaryote. Les cellules S. pombe ont la forme de bâtonnets et croissent par l'allongement de leurs extrémités. Elles se divisent en formant, en leur milieu une paroi cellulaire, appelé septum, permettant ainsi l'obtention de deux cellules filles de même taille. Chez S. pombe, la cytokinèse commence en début de mitose lorsque le site de division est déterminé par la formation d'un anneau d'acto-myosine. Le septum, lui, est formé uniquement en fin de mitose par la contraction de l'anneau d'actomyosine. Cette contraction est sous le contrôle d'un réseau de signalisation cellulaire appelé le «réseau d'initiation de synthèse du septum » ou « septation initiation network » (SIN), qui se situe sur les pôles du fuseau mitotique. L'activation du SIN dépend d'une GTPase appelé spg1p dont le statut nucléotidique dépend des protéines cdclóp et byr4p qui forment un complexe qui favorise l'hydrolyse du GTP en GDP. En début de mitose, cdc16p ne se situe plus sur les poles du fuseau mitotique. La GTPase spg1p se retrouve donc principalement sous sa forme couplée au GTP, ce quí va permettre son interaction avec la kinase cdc7p. Cette protéine ainsi que deux autres kinases sid2p (avec mob1p) et sid1p (avec cdc14p) permettent la transmission du signal d'initiation de la contraction de l'anneau d'acto-myosine en fin d'anaphase. Pendant l'anaphase, cdc7p, sid1 p et cdc14p localisent sur un des deux pôles du fuseau mitotique. Il en est de même pour cdc1p et by14p et le pôle contenant cdc16p et byr4p est toujours différent de celui ou les régulateurs positifs du SIN se situent. En fin de cytokinèse, cdc16 et byr4p se retrouvent à nouveau sur chaque pôle des deux cellules filles. Dans cette étude, nous nous sommes concentrés sur l'analyse de la localisation de cdc16p pendant la mitose en utilisant une technique de microscopie en temps réel. Nous avons été en mesure de déterminer que le départ de cdc16p du pole s'effectue juste avant la formation du fuseau mitotique. Nous avons aussi découvert que la localisation asymétrique des composants du SIN dépend fortement de l'entrée en anaphase B. Finalement, Nous avons montré que distribution asymétrique des composants du SIN sur les pôles du fuseau mitotique dépendait aussi fortement de !a présence de l'anneau d'acto-myosine. Ceci nous permet donc de proposer l'existence d'un mécanisme cellulaire qui permet de s'assurer que la cytokinèse est achevée avant de diminuer la signalisation du SIN. Par ailleurs, des études biochimiques nous ont permis de montrer que cdc16p est dégradé par le proteosome. Ces travaux ont permis la découverte de nouveaux modes de régulation du SIN.

Risk analysis of climate change on coffee nematodes and leaf miner in Brazil

The objective of this work was to assess the potential impact of climate change on the spatial distribution of coffee nematodes (races of Meloidogyne incognita) and leaf miner (Leucoptera coffeella), using a Geographic Information System. Assessment of the impacts of climate change on pest infestations and disease epidemics in crops is needed as a basis for revising management practices to minimize crop losses as climatic conditions shift. Future scenarios focused on the decades of the 2020's, 2050's, and 2080's (scenarios A2 and B2) were obtained from five General Circulation Models available on Data Distribution Centre from Intergovernmental Panel on Climate Change. Geographic distribution maps were prepared using models to predict the number of generations of the nematodes and leaf miner. Maps obtained in scenario A2 allowed prediction of an increased infestation of the nematode and of the pest, due to greater number of generations per month, than occurred under the climatological normal from 1961-1990. The number of generations also increased in the B2 scenario, but was lower than in the A2 scenario for both organisms.

Classification of very high spatial resolution imagery using mathematical morphology and support vector machines

We investigate the relevance of morphological operators for the classification of land use in urban scenes using submetric panchromatic imagery. A support vector machine is used for the classification. Six types of filters have been employed: opening and closing, opening and closing by reconstruction, and opening and closing top hat. The type and scale of the filters are discussed, and a feature selection algorithm called recursive feature elimination is applied to decrease the dimensionality of the input data. The analysis performed on two QuickBird panchromatic images showed that simple opening and closing operators are the most relevant for classification at such a high spatial resolution. Moreover, mixed sets combining simple and reconstruction filters provided the best performance. Tests performed on both images, having areas characterized by different architectural styles, yielded similar results for both feature selection and classification accuracy, suggesting the generalization of the feature sets highlighted.

Spatial visualization of DNA in solution.

An image analysis method is presented which allows for the reconstruction of the three-dimensional path of filamentous objects from two of their projections. Starting with stereo pairs, this method is used to trace the trajectory of DNA molecules embedded in vitreous ice and leads to a faithful representation of their three-dimensional shape in solution. This computer-aided reconstruction is superior to the subjective three-dimensional impression generated by observation of stereo pairs of micrographs because it enables one to look at the reconstructed molecules from any chosen direction and distance and allows quantitative analysis such as determination of distances, curvature, persistence length, and writhe of DNA molecules in solution.