Advances in high-throughput DNA sequencing technologies have determined an explosion in the number of sequenced bacterial genomes. Comparative sequence analysis frequently reveals evidences of homologous recombination occurring with different mechanisms and rates in different species, but the large-scale use of computational methods to identify recombination events is hampered by their high computational costs. Here, we propose a new method to identify recombination events in large datasets of whole genome sequences. Using a filtering procedure of the gene conservation profiles of a test genome against a panel of strains, this algorithm identifies sets of contiguous genes acquired by homologous recombination. The locations of the recombination breakpoints are determined using a statistical test that is able to account for the differences in the natural rate of evolution between different genes. The algorithm was tested on a dataset of 75 genomes of Staphylococcus aureus and 50 genomes comprising different streptococcal species, and was able to detect intra-species recombination events in S. aureus and in Streptococcus pneumoniae. Furthermore, we found evidences of an inter-species exchange of genetic material between S. pneumoniae and Streptococcus mitis, a closely related commensal species that colonizes the same ecological niche. The method has been implemented in an R package, Reco, which is freely available from supplementary material, and provides a rapid screening tool to investigate recombination on a genome-wide scale from sequence data.

Abstract Motivation: Eu.Gene Analyzer is an easy-to-use, stand-alone application that allows rapid and powerful microarray data analysis in the context of biological pathways. Its intuitive graphical user interface makes it an easy and flexible tool, even for the first-time user. Eu.Gene supports a variety of array platforms, organisms and pathway ontologies, transparently deals with multiple nomenclature systems and seamlessly integrates data from different sources. Two different statistical methods, the Fisher Exact Test and the Gene Set Enrichment Analysis (GSEA), are implemented to identify biological pathways transcriptionally affected under experimental conditions. A suite of tools is offered to define, visualize and share custom non-redundant pathway sets. In conclusion, Eu.Gene Analyzer is a new software application that takes advantage of information from multiple pathway databases to build a comprehensive interpretation of experimental results in a simple, intuitive environment. Availability: Download of Eu.Gene Analyzer Java version is available free of charge for academic users. Please visit the web page: http://www.ducciocavalieri.org/bio/Eugene.htm Contact: [email protected] Supplementary information: http://www.ducciocavalieri.org/bio/Eugene/Suppl_Inf

MOTIVATION: Microarray studies permit to quantify expression levels on a global scale by measuring transcript abundance of thousands of genes simultaneously. A difficulty when analysing expression measures is how to model variability for the whole set of genes. It is usually unrealistic to assume a common variance for each gene. Several approaches to model gene-specific variances are proposed. We take advantage of calibration experiments, in which the probes hybridized on the two channels come from the same population (self-self experiment). In this case it is possible to estimate the gene-specific variance, to be incorporated in comparative experiments on the same tissue, cellular line or species. RESULTS: We present two approaches to introduce prior information on gene-specific variability from a calibration experiment: an empirical Bayes model and a full Bayesian hierarchical model. We apply the methods in the analysis of human lipopolysaccharide-stimulated leukocyte experiments. AVAILABILITY: The calculations are implemented in WinBugs. The codes are available on request from the authors.