Modular analysis of gene expression data with R.

SUMMARY: Large sets of data, such as expression profiles from many samples, require analytic tools to reduce their complexity. The Iterative Signature Algorithm (ISA) is a biclustering algorithm. It was designed to decompose a large set of data into so-called 'modules'. In the context of gene expression data, these modules consist of subsets of genes that exhibit a coherent expression profile only over a subset of microarray experiments. Genes and arrays may be attributed to multiple modules and the level of required coherence can be varied resulting in different 'resolutions' of the modular mapping. In this short note, we introduce two BioConductor software packages written in GNU R: The isa2 package includes an optimized implementation of the ISA and the eisa package provides a convenient interface to run the ISA, visualize its output and put the biclusters into biological context. Potential users of these packages are all R and BioConductor users dealing with tabular (e.g. gene expression) data. AVAILABILITY: http://www.unil.ch/cbg/ISA CONTACT: sven.bergmann@unil.ch

Analysis of angiotensin II- and ACTH-driven mineralocorticoid functions and omental adiposity in a non-genetic, hyperadipose female rat phenotype

The hypothalamic damage induced by neonatal treatment with monosodium l-glutamate (MSG) induces several metabolic abnormalities, resulting in a rat hyperleptinemic-hyperadipose phenotype. This study was conducted to explore the impact of the neonatal MSG treatment, in the adult (120 days old) female rat on: (a) the in vivo and in vitro mineralocorticoid responses to ACTH and angiotensin II (AII); (b) the effect of leptin on ACTH- and AII-stimulated mineralocorticoid secretions by isolated corticoadrenal cells; and (c) abdominal adiposity characteristics. Our data indicate that, compared with age-matched controls, MSG rats displayed: (1) enhanced and reduced mineralocorticoid responses to ACTH and AII treatments, respectively, effects observed in both in vivo and in vitro conditions; (2) adrenal refractoriness to the inhibitory effect of exogenous leptin on ACTH-stimulated aldosterone output by isolated adrenocortical cells; and (3) distorted omental adiposity morphology and function. This study supports that the adult hyperleptinemic MSG female rat is characterized by enhanced ACTH-driven mineralocorticoid function, impaired adrenal leptin sensitivity, and disrupted abdominal adiposity function. MSG rats could counteract undesirable effects of glucocorticoid excess, by developing a reduced AII-driven mineralocorticoid function. Thus, chronic hyperleptinemia could play a protective role against ACTH-mediated allostatic loads in the adrenal leptin resistant, MSG female rat phenotype.

Diseño e implementación de un framework de persistencia

Los Frameworks constituyen el nuevo paradigma en cuanto al desarrollo de software se refiere. Entre sus principales características se encuentran la facilidad para la reutilización de código. En este marco específico proporcionados por la tecnología usaremos la tecnología JAVA y su extensión en cuanto a la persistencia de datos. El Framework de persistencia es el responsable de gestionar la lógica de acceso a los datos en un SGBD (Sistema de Gestión de Bases de Datos), ya sea de entrada o salida, y ocultando los detalles más pesados relativos a la estructura propia de la Base de Datos utilizada, de manera completa y transparentemente. En conclusión, este proyecto se basa en un análisis del los Frameworks existentes, analizando sus características y profundizando en los detalles concretos de su actividad y manejo en cuanto a la persistencia.

CODAPACK3D. A new version of Compositional Data Package

The statistical analysis of compositional data should be treated using logratios of parts,which are difficult to use correctly in standard statistical packages. For this reason afreeware package, named CoDaPack was created. This software implements most of thebasic statistical methods suitable for compositional data.In this paper we describe the new version of the package that now is calledCoDaPack3D. It is developed in Visual Basic for applications (associated with Excel©),Visual Basic and Open GL, and it is oriented towards users with a minimum knowledgeof computers with the aim at being simple and easy to use.This new version includes new graphical output in 2D and 3D. These outputs could bezoomed and, in 3D, rotated. Also a customization menu is included and outputs couldbe saved in jpeg format. Also this new version includes an interactive help and alldialog windows have been improved in order to facilitate its use.To use CoDaPack one has to access Excel© and introduce the data in a standardspreadsheet. These should be organized as a matrix where Excel© rows correspond tothe observations and columns to the parts. The user executes macros that returnnumerical or graphical results. There are two kinds of numerical results: new variablesand descriptive statistics, and both appear on the same sheet. Graphical output appearsin independent windows. In the present version there are 8 menus, with a total of 38submenus which, after some dialogue, directly call the corresponding macro. Thedialogues ask the user to input variables and further parameters needed, as well aswhere to put these results. The web site http://ima.udg.es/CoDaPack contains thisfreeware package and only Microsoft Excel© under Microsoft Windows© is required torun the software.Kew words: Compositional data Analysis, Software

Analysis of GPCR properties of Frizzled receptors and establishment of the "in vitro" platform for screening of Frizzled agonists/antagonists as potential therapeutic agents

Morphogens of the Wnt protein family are the secreted lipoglycoprotein ligands which initiate several pathways heavily involved in the coordination of various developmental stages of organisms in the majority of animal species. Deregulation of these pathways in the adult leads to formation and sustaining of multiple types of cancer. The latter notion is reinforced by the fact that the very discovery of the first Wnt ligand was due to its role as the causative factor of carcinogenic transformation (Nusse and Varmus, 1982). Nowadays our knowledge on Wnt signaling has "moved with the times" and these pathways were identified to be often crucial for tumor formation, its interactions with the microenvironment, and promotion of the metastases (Huang and Du, 2008; Zerlin et al., 2008; Jessen, 2009). Thus the relevance of the pathway as the target for drug development has further increased in the light of modern paradigms of the complex cancer treatments which target also spreading and growth- promoting factors of tumors by specific and highly efficient substances (Pavet et al., 2010). Presently the field of the Wnt-targeting drug research is almost solely dominated by assays based on transcriptional activation induced by the signaling. This approach resulted in development of a number of promising substances (Lee et al., 2011). Despite its effectiveness, the method nevertheless suffers from several drawbacks. Among the major ones is the fact that this approach is prone to identify compounds targeting rather downstream effectors of the pathway, which are indiscriminately used by all the subtypes of the Wnt signaling. Additionally, proteins which are involved in several signaling cascades and not just the Wnt pathway turn out as targets of the new compounds. These issues increase risks of side effects due to off-target interactions and blockade of the pathway in healthy cells. In the present work we put forward a novel biochemical approach for drug development on the Wnt pathway. It targets Frizzleds (Fzs) - a family of 7-transmbembrane proteins which serve as receptors for Wnt ligands. They offer unique properties for the development of highly specific and effective drugs as they control all branches of the Wnt signaling. Recent advances in the understanding of the roles of heterotrimeric G proteins downstream from Fzs (Katanaev et al., 2005; Liu et al., 2005; Jernigan et al., 2010) suggest application of enzymatic properties of these effectors to monitor the receptor-mediated events. We have applied this knowledge in practice and established a specific and efficient method based on utilization of a novel high-throughput format of the GTP-binding assay to follow the activation of Fzs. This type of assay is a robust and well-established technology for the research and screenings on the GPCRs (Harrison and Traynor, 2003). The conventional method of detection involves the radioactively labeled non-hydrolysable GTP analog [35S]GTPyS. Its application in the large-scale screenings is however problematic which promoted development of the novel non-radioactive GTP analog GTP-Eu. The new molecule employs phenomenon of the time-resolved fluorescence to provide sensitivity comparable to the conventional radioactive substance. Initially GTP-Eu was tested only in one of many possible types of GTP-binding assays (Frang et al., 2003). In the present work we expand these limits by demonstrating the general comparability of the novel label with the radioactive method in various types of assays. We provide a biochemical characterization of GTP-Eu interactions with heterotrimeric and small GTPases and a comparative analysis of the behavior of the new label in the assays involving heterotrimeric G protein effectors. These developments in the GTP-binding assay were then applied to monitor G protein activation by the Fz receptors. The data obtained in mammalian cultured cell lines provides for the first time an unambiguous biochemical proof for direct coupling of Fzs with G proteins. The specificity of this interaction has been confirmed by the experiments with the antagonists of Fz and by the pertussis toxin-mediated deactivation. Additionally we have identified the specificity of Wnt3a towards several members of the Fz family and analyzed the properties of human Fz-1 which was found to be the receptor coupled to the Gi/o family of G proteins. Another process playing significant role in the functioning of every GPCR is endocytosis. This phenomenon can also be employed for drug screenings on GPCRs (Bickle, 2010). In the present work we have demonstrated that Drosophila Fz receptors are involved in an unusual for many GPCRs manifestation of the receptor-mediated internalization. Through combination of biochemical approaches and studies on Drosophila as the model organism we have shown that direct interactions of the Fzs and the α-subunit of the heterotrimeric G protein Go with the small GTPase Rab5 regulate internalization of the receptor in early endosomes. We provide data uncovering the decisive role of this self-promoted endocytosis in formation of a proper signaling output in the canonical as well as planar cell polarity (PCP) pathways regulated by Fz. The results of this work thus establish a platform for the high-throughput screening to identify substances active in the cancer-related Wnt pathways. This methodology has been adjusted and applied to provide the important insights in Fz functioning and will be instrumental for further investigations on the Wnt-mediated pathways.

Controlled contrast transcranial Doppler and arterial blood gas analysis to quantify shunt through patent foramen ovale.

BACKGROUND AND PURPOSE: A right-to-left shunt can be identified by contrast transcranial Doppler ultrasonography (c-TCD) at rest and/or after a Valsalva maneuver (VM) or by arterial blood gas (ABG) measurement. We assessed the influence of controlled strain pressures and durations during VM on the right-to-left passage of microbubbles, on which depends the shunt classification by c-TCD, and correlated it with the right-to-left shunt evaluation by ABG measurements in stroke patients with patent foramen ovale (PFO). METHODS: We evaluated 40 stroke patients with transesophageal echocardiography-documented PFO. The microbubbles were recorded with TCD at rest and after 4 different VM conditions with controlled duration and target strain pressures (duration in seconds and pressure in cm H2O, respectively): V5-20, V10-20, V5-40, and V10-40. The ABG analysis was performed after pure oxygen breathing in 34 patients, and the shunt was calculated as percentage of cardiac output. RESULTS: Among all VM conditions, V5-40 and V10-40 yielded the greatest median number of microbubbles (84 and 95, respectively; P<0.01). A significantly larger number of microbubbles were detected in V5-40 than in V5-20 (P<0.001) and in V10-40 than in V10-20 (P<0.01). ABG was not sensitive enough to detect a shunt in 31 patients. CONCLUSIONS: The increase of VM expiratory pressure magnifies the number of microbubbles irrespective of the strain duration. Because the right-to-left shunt classification in PFO is based on the number of microbubbles, a controlled VM pressure is advised for a reproducible shunt assessment. The ABG measurement is not sensitive enough for shunt assessment in stroke patients with PFO.

An Analysis of black-box optimization problems in reinsurance: evolutionary-based approaches

Black-box optimization problems (BBOP) are de ned as those optimization problems in which the objective function does not have an algebraic expression, but it is the output of a system (usually a computer program). This paper is focussed on BBOPs that arise in the eld of insurance, and more speci cally in reinsurance problems. In this area, the complexity of the models and assumptions considered to de ne the reinsurance rules and conditions produces hard black-box optimization problems, that must be solved in order to obtain the optimal output of the reinsurance. The application of traditional optimization approaches is not possible in BBOP, so new computational paradigms must be applied to solve these problems. In this paper we show the performance of two evolutionary-based techniques (Evolutionary Programming and Particle Swarm Optimization). We provide an analysis in three BBOP in reinsurance, where the evolutionary-based approaches exhibit an excellent behaviour, nding the optimal solution within a fraction of the computational cost used by inspection or enumeration methods.

Time Series Input Selection using Multiple Kernel Learning

In this paper we study the relevance of multiple kernel learning (MKL) for the automatic selection of time series inputs. Recently, MKL has gained great attention in the machine learning community due to its flexibility in modelling complex patterns and performing feature selection. In general, MKL constructs the kernel as a weighted linear combination of basis kernels, exploiting different sources of information. An efficient algorithm wrapping a Support Vector Regression model for optimizing the MKL weights, named SimpleMKL, is used for the analysis. In this sense, MKL performs feature selection by discarding inputs/kernels with low or null weights. The approach proposed is tested with simulated linear and nonlinear time series (AutoRegressive, Henon and Lorenz series).

General trends in world biogeographic literature: a preliminary bibliometric analysis

Records with the search string biogeograph* were collected from the Science Citation Index (SCI). A total of 3456 records were downloaded for the 1945-2006 period from titles of articles and reviews, and 10,543 records were downloaded for 1991-2006, taking into consideration also abstracts and keywords. Temporal trends of publications, geographical and institutional distribution of the research output, authorship, and core journals were evaluated. There were as many as 122 countries carrying out biogeographic research; in the most recent period, USA is the top producing country, followed by the United Kingdom, Australia, France, Germany, Spain, and Canada. There were 17,493 authors contributing to the field. During 1991-2006 there were 4098 organizations with authors involved in biogeographic research; institutions with higher number of papers are the Natural History Museum (United Kingdom), the University of California, Berkeley (USA), the Museum National d'Histoire Naturelle (France), the Universidad Nacional Autónoma de México (Mexico), the American Museum of Natural History (USA) and the Russian Academy of Sciences (Russia). Research articles are spread over a variety of journals, with the Journal of Biogeography, Molecular Phylogenetics and Evolution, Molecular Ecology, and Biological Journal of the Linnean Society being the core journals. From 28,759 keywords retrieved those with the highest frequency were evolution, phylogeny, diversity, mitochondrial DNA, pattern(s), systematics, and population(s). We conclude that publications on biogeography have increased substantially during the last years, especially since 1998. The preferred journal for biogeographic papers is the Journal of Biogeography. Most frequent keywords seem to indicate that biogeography fits well within both evolutionary biology and ecology, with molecular biology and phylogenetics being important factors that drive their current development.

What explains the Great Moderation in the US? A structural analysis

This paper investigates what has caused output and inflation volatility to fall in the USusing a small scale structural model using Bayesian techniques and rolling samples. Thereare instabilities in the posterior of the parameters describing the private sector, the policyrule and the standard deviation of the shocks. Results are robust to the specification ofthe policy rule. Changes in the parameters describing the private sector are the largest,but those of the policy rule and the covariance matrix of the shocks explain the changes most.

On the robust effects of technology shocks on hours worked and output

We analyze the effects of neutral and investment-specific technology shockson hours and output. Long cycles in hours are captured in a variety of ways.Hours robustly fall in response to neutral shocks and robustly increase inresponse to investment specific shocks. The percentage of the variance ofhours (output) explained by neutral shocks is small (large); the opposite istrue for investment specific shocks. News shocks are uncorrelated with theestimated technology shocks.

Inflation dynamics: A structural econometric analysis

We develop and estimate a structural model of inflation that allowsfor a fraction of firms that use a backward looking rule to setprices. The model nests the purely forward looking New KeynesianPhillips curve as a particular case. We use measures of marginalcosts as the relevant determinant of inflation, as the theorysuggests, instead of an ad-hoc output gap. Real marginal costsare a significant and quantitatively important determinant ofinflation. Backward looking price setting, while statisticallysignificant, is not quantitatively important. Thus, we concludethat the New Keynesian Phillips curve provides a good firstapproximation to the dynamics of inflation.

GPS analysis of human locomotion: further evidence for long-range correlations in stride-to-stride fluctuations of gait parameters.

During free walking, gait is automatically adjusted to provide optimal mechanical output and minimal energy expenditure; gait parameters, such as cadence, fluctuate from one stride to the next around average values. It was described that this fluctuation exhibited long-range correlations and fractal-like patterns. In addition, it was suggested that these long-range correlations disappeared if the participant followed the beep of metronome to regulate his or her pace. Until now, these fractal fluctuations were only observed for stride interval, because no technique existed to adequately analyze an extended time of free walking. The aim of the present study was to measure walking speed (WS), step frequency (SF) and step length (SL) with high accuracy (<1 cm) satellite positioning method (global positioning system or GPS) in order to detect long-range correlations in the stride-to-stride fluctuations. Eight participants walked 30 min under free and constrained (metronome) conditions. Under free walking conditions, DFA (detrended fluctuation analysis) and surrogate data tests showed that the fluctuation of WS, SL and SF exhibited a fractal pattern (i.e., scaling exponent alpha: 0.5 < alpha < 1) in a large majority of participants (7/8). Under constrained conditions (metronome), SF fluctuations became significantly anti-correlated (alpha < 0.5) in all participants. However, the scaling exponent of SL and WS was not modified. We conclude that, when the walking pace is controlled by an auditory signal, the feedback loop between the planned movement (at supraspinal level) and the sensory inputs induces a continual shifting of SF around the mean (persistent anti-correlation), but with no effect on the fluctuation dynamics of the other parameters (SL, WS).

Post-detection analysis for grating-based ultra-small angle X-ray scattering.

Until recently, the hard X-ray, phase-sensitive imaging technique called grating interferometry was thought to provide information only in real space. However, by utilizing an alternative approach to data analysis we demonstrated that the angular resolved ultra-small angle X-ray scattering distribution can be retrieved from experimental data. Thus, reciprocal space information is accessible by grating interferometry in addition to real space. Naturally, the quality of the retrieved data strongly depends on the performance of the employed analysis procedure, which involves deconvolution of periodic and noisy data in this context. The aim of this article is to compare several deconvolution algorithms to retrieve the ultra-small angle X-ray scattering distribution in grating interferometry. We quantitatively compare the performance of three deconvolution procedures (i.e., Wiener, iterative Wiener and Lucy-Richardson) in case of realistically modeled, noisy and periodic input data. The simulations showed that the algorithm of Lucy-Richardson is the more reliable and more efficient as a function of the characteristics of the signals in the given context. The availability of a reliable data analysis procedure is essential for future developments in grating interferometry.

10

This contribution introduces Data Envelopment Analysis (DEA), a performance measurement technique. DEA helps decision makers for the following reasons: (1) By calculating an efficiency score, it indicates if a firm is efficient or has capacity for improvement; (2) By setting target values for input and output, it calculates how much input must be decreased or output increased in order to become efficient; (3) By identifying the nature of returns to scale, it indicates if a firm has to decrease or increase its scale (or size) in order to minimise the average total cost; (4) By identifying a set of benchmarks, it specifies which other firms' processes need to be analysed in order to improve its own practices. This contribution presents the essentials about DEA, alongside a case study to intuitively understand its application. It also introduces Win4DEAP, a software package that conducts efficiency analysis based on DEA methodology. The methodical background of DEA is presented for more demanding readers. Finally, four advanced topics of DEA are treated: adjustment to the environment, preferences, sensitivity analysis and time series data.