Rgtsp: a generalized top scoring pairs package for class prediction.

SUMMARY: A top scoring pair (TSP) classifier consists of a pair of variables whose relative ordering can be used for accurately predicting the class label of a sample. This classification rule has the advantage of being easily interpretable and more robust against technical variations in data, as those due to different microarray platforms. Here we describe a parallel implementation of this classifier which significantly reduces the training time, and a number of extensions, including a multi-class approach, which has the potential of improving the classification performance. AVAILABILITY AND IMPLEMENTATION: Full C++ source code and R package Rgtsp are freely available from http://lausanne.isb-sib.ch/~vpopovic/research/. The implementation relies on existing OpenMP libraries.

Density-based hierarchical clustering of pyro-sequences on a large scale--the case of fungal ITS1.

MOTIVATION: Analysis of millions of pyro-sequences is currently playing a crucial role in the advance of environmental microbiology. Taxonomy-independent, i.e. unsupervised, clustering of these sequences is essential for the definition of Operational Taxonomic Units. For this application, reproducibility and robustness should be the most sought after qualities, but have thus far largely been overlooked. RESULTS: More than 1 million hyper-variable internal transcribed spacer 1 (ITS1) sequences of fungal origin have been analyzed. The ITS1 sequences were first properly extracted from 454 reads using generalized profiles. Then, otupipe, cd-hit-454, ESPRIT-Tree and DBC454, a new algorithm presented here, were used to analyze the sequences. A numerical assay was developed to measure the reproducibility and robustness of these algorithms. DBC454 was the most robust, closely followed by ESPRIT-Tree. DBC454 features density-based hierarchical clustering, which complements the other methods by providing insights into the structure of the data. AVAILABILITY: An executable is freely available for non-commercial users at ftp://ftp.vital-it.ch/tools/dbc454. It is designed to run under MPI on a cluster of 64-bit Linux machines running Red Hat 4.x, or on a multi-core OSX system. CONTACT: dbc454@vital-it.ch or nicolas.guex@isb-sib.ch.

Rapid and selective surveillance of Arabidopsis thaliana genome annotations with Centrifuge.

Centrifuge is a user-friendly system to simultaneously access Arabidopsis gene annotations and intra- and inter-organism sequence comparison data. The tool allows rapid retrieval of user-selected data for each annotated Arabidopsis gene providing, in any combination, data on the following features: predicted protein properties such as mass, pI, cellular location and transmembrane domains; SWISS-PROT annotations; Interpro domains; Gene Ontology records; verified transcription; BLAST matches to the proteomes of A.thaliana, Oryza sativa (rice), Caenorhabditis elegans, Drosophila melanogaster and Homo sapiens. The tool lends itself particularly well to the rapid analysis of contigs or of tens or hundreds of genes identified by high-throughput gene expression experiments. In these cases, a summary table of principal predicted protein features for all genes is given followed by more detailed reports for each individual gene. Centrifuge can also be used for single gene analysis or in a word search mode. AVAILABILITY: http://centrifuge.unil.ch/ CONTACT: edward.farmer@unil.ch.

Histone deacetylase inhibitors impair antibacterial defenses of macrophages.

Histone deacetylases (HDACs) control gene expression by deacetylating histones and nonhistone proteins. HDAC inhibitors (HDACi) are powerful anticancer drugs that exert anti-inflammatory and immunomodulatory activities. We recently reported a proof-of-concept study demonstrating that HDACi increase susceptibility to bacterial infections in vivo. Yet, still little is known about the effects of HDACi on antimicrobial innate immune defenses. Here we show that HDACi belonging to different chemical classes inhibit at multiple levels the response of macrophages to bacterial infection. HDACi reduce the phagocytosis and the killing of Escherichia coli and Staphylococcus aureus by macrophages. In line with these findings, HDACi decrease the expression of phagocytic receptors and inhibit bacteria-induced production of reactive oxygen and nitrogen species by macrophages. Consistently, HDACi impair the expression of nicotinamide adenine dinucleotide phosphate (NADPH) oxidase subunits and inducible nitric oxide synthase. These data indicate that HDACi have a strong impact on critical antimicrobial defense mechanisms in macrophages.

Linear relations between writhe and minimal crossing number in Conway families of ideal knots and links

We showed earlier how to predict the writhe of any rational knot or link in its ideal geometric configuration, or equivalently the average of the 3D writhe over statistical ensembles of random configurations of a given knot or link (Cerf and Stasiak 2000 Proc. Natl Acad. Sci. USA 97 3795). There is no general relation between the minimal crossing number of a knot and the writhe of its ideal geometric configuration. However, within individual families of knots linear relations between minimal crossing number and writhe were observed (Katritch et al 1996 Nature 384 142). Here we present a method that allows us to express the writhe as a linear function of the minimal crossing number within Conway families of knots and links in their ideal configuration. The slope of the lines and the shift between any two lines with the same

Manipulating sex ratio to increase population growth: the example of the Lesser Kestrel

Small or decreasing populations call for emergency actions like, for example, captive breeding programs. Such programs aim at rapidly increasing population sizes in order to reduce the loss of genetic variability and to avoid possible Allee effects. The Lesser Kestrel Falco naumanni is one of the species that is currently supported in several captive breeding programs at various locations. Here, we model the demographic and genetic consequences of potential management strategies that are based on offspring sex ratio manipulation. Increased population growth could be achieved by manipulating female conditions and/or male attractiveness in the captive breeders and consequently shifting the offspring sex ratio towards more female offspring, which are then used for reintroduction. Fragmenting populations into wild-breeding and captive-breeding demes and manipulating population sex ratio both immediately increase the inbreeding coefficient in the next generation (i.e. decrease N-e) but may, in the long term, reduce the loss of genetic variability if population growth is restricted by the number of females. We use the Lesser Kestrel and the wealth of information that is available on this species to predict the long-term consequences of various kinds of sex-ratio manipulation. We find that, in our example and possibly in many other cases, a sex-ratio manipulation that seems realistic could have a beneficial effect on the captive breeding program. However, the possible long-term costs and benefits of such measures need to be carefully optimized.

Health technology assessment in evolution - focal therapy in localised prostate cancer.

Focal therapy in prostate cancer aims to treat only the part of the gland harboring clinically significant disease while preserving the rest of the tissue. This approach may substantially reduce treatment-related toxicity without compromising disease control outcomes. Short- to medium-term functional and oncological results in prospective interventional studies are promising, but comparative effectiveness research against standard of care is required to incorporate focal therapy among standard options. In this review, we discuss the actual stage of assessment and results of sources of energy commonly used to deliver focal therapy. We also provide our viewpoint on how the field will evolve in the near future.

FastEpistasis: a high performance computing solution for quantitative trait epistasis.

Motivation: Genome-wide association studies have become widely used tools to study effects of genetic variants on complex diseases. While it is of great interest to extend existing analysis methods by considering interaction effects between pairs of loci, the large number of possible tests presents a significant computational challenge. The number of computations is further multiplied in the study of gene expression quantitative trait mapping, in which tests are performed for thousands of gene phenotypes simultaneously. Results: We present FastEpistasis, an efficient parallel solution extending the PLINK epistasis module, designed to test for epistasis effects when analyzing continuous phenotypes. Our results show that the algorithm scales with the number of processors and offers a reduction in computation time when several phenotypes are analyzed simultaneously. FastEpistasis is capable of testing the association of a continuous trait with all single nucleotide polymorphism ( SNP) pairs from 500 000 SNPs, totaling 125 billion tests, in a population of 5000 individuals in 29, 4 or 0.5 days using 8, 64 or 512 processors.

Long-distance wolf recolonization of France and Switzerland inferred from non-invasive genetic sampling over a period of 10 years

In the early 1900s, the wolf (Canis lupus) was extirpated from France and Switzerland. There is growing evidence that the species is presently recolonizing these countries in the western Alps. By sequencing the mitochondrial DNA (mtDNA) control region of various samples mainly collected in the field (scats, hairs, regurgitates, blood or tissue; n = 292), we could (1) develop a non-invasive method enabling the unambiguous attribution of these samples to wolf, fox (Vulpes vulpes) or dog (Canis familiaris), among others; (2) demonstrate that Italian, French and Swiss wolves share the same mtDNA haplotype, a haplotype that has never been found in any other wolf population world-wide. Combined together, field and genetic data collected over 10 years corroborate the scenario of a natural expansion of wolves from the Italian source population. Furthermore, such a genetic approach is of conservation significance, since it has important consequences for management decisions. This first long-term report using non-invasive sampling demonstrates that long-distance dispersers are common, supporting the hypothesis that individuals may often attempt to colonize far from their native pack, even in the absence of suitable corridors across habitats characterized by intense human activities.

Peroxisome proliferator-activated receptor agonists.

The peroxisome proliferator-activated receptors are a family of three ligand-activated transcription factors. Fibrate antihyperlipidemic drugs and thiazolidinedione antihyperglycemic drugs were recently identified as synthetic ligands for these receptors. In addition, certain unsaturated fatty acids and eicosanoids were shown to bind the receptors, and thus represent naturally occurring PPAR ligands. The synthetic and natural ligands have proven to be powerful tools in dissecting the biology of these orphan receptors.

Evolutionary footprint of coevolving positions in genes.

MOTIVATION: The analysis of molecular coevolution provides information on the potential functional and structural implication of positions along DNA sequences, and several methods are available to identify coevolving positions using probabilistic or combinatorial approaches. The specific nucleotide or amino acid profile associated with the coevolution process is, however, not estimated, but only known profiles, such as the Watson-Crick constraint, are usually considered a priori in current measures of coevolution. RESULTS: Here, we propose a new probabilistic model, Coev, to identify coevolving positions and their associated profile in DNA sequences while incorporating the underlying phylogenetic relationships. The process of coevolution is modeled by a 16 × 16 instantaneous rate matrix that includes rates of transition as well as a profile of coevolution. We used simulated, empirical and illustrative data to evaluate our model and to compare it with a model of 'independent' evolution using Akaike Information Criterion. We showed that the Coev model is able to discriminate between coevolving and non-coevolving positions and provides better specificity and specificity than other available approaches. We further demonstrate that the identification of the profile of coevolution can shed new light on the process of dependent substitution during lineage evolution.

A qualitative continuous model of cellular auxin and brassinosteroid signaling and their crosstalk.

Motivation: Hormone pathway interactions are crucial in shaping plant development, such as synergism between the auxin and brassinosteroid pathways in cell elongation. Both hormone pathways have been characterized in detail, revealing several feedback loops. The complexity of this network, combined with a shortage of kinetic data, renders its quantitative analysis virtually impossible at present.Results: As a first step towards overcoming these obstacles, we analyzed the network using a Boolean logic approach to build models of auxin and brassinosteroid signaling, and their interaction. To compare these discrete dynamic models across conditions, we transformed them into qualitative continuous systems, which predict network component states more accurately and can accommodate kinetic data as they become available. To this end, we developed an extension for the SQUAD software, allowing semi-quantitative analysis of network states. Contrasting the developmental output depending on cell type-specific modulators enabled us to identify a most parsimonious model, which explains initially paradoxical mutant phenotypes and revealed a novel physiological feature.