Search

Motivation: Protein sequence alignments are essential to structural, evolutionary and functional analysis, but their accuracy is often limited by sequence similarity unless molecular structures are available. Protein structures predicted at experimental grade accuracy, as achieved by AlphaFold2, could therefore have a major impact on sequence analysis., Results: Here, we find that multiple sequence alignments estimated on AlphaFold2 predictions are almost as accurate as alignments estimated on experimental structures and significantly closer to the structural reference than sequence-based alignments. We also show that AlphaFold2 structural models of relatively low quality can be used to obtain highly accurate alignments. These results suggest that, besides structure modeling, AlphaFold2 encodes higher-order dependencies that can be exploited for sequence analysis., Availability and Implementation: All data, analyses and results are available on Zenodo (https://doi.org/10.5281/zenodo.7031286). The code and scripts have been deposited in GitHub (https://github.com/cbcrg/msa-af2-nf) and the various containers in (https://cloud.sylabs.io/library/athbaltzis/af2/alphafold, https://hub.docker.com/r/athbaltzis/pred)., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2022. Published by Oxford University Press.)

Aneuploidy causes system-wide disruptions in the stochiometric balances of transcripts, proteins, and metabolites, often resulting in detrimental effects for the organism. The protozoan parasite Leishmania has an unusually high tolerance for aneuploidy, but the molecular and functional consequences for the pathogen remain poorly understood. Here, we addressed this question in vitro and present the first integrated analysis of the genome, transcriptome, proteome, and metabolome of highly aneuploid Leishmania donovani strains. Our analyses unambiguously establish that aneuploidy in Leishmania proportionally impacts the average transcript- and protein abundance levels of affected chromosomes, ultimately correlating with the degree of metabolic differences between closely related aneuploid strains. This proportionality was present in both proliferative and non-proliferative in vitro promastigotes. However, as in other Eukaryotes, we observed attenuation of dosage effects for protein complex subunits and in addition, non-cytoplasmic proteins. Differentially expressed transcripts and proteins between aneuploid Leishmania strains also originated from non-aneuploid chromosomes. At protein level, these were enriched for proteins involved in protein metabolism, such as chaperones and chaperonins, peptidases, and heat-shock proteins. In conclusion, our results further support the view that aneuploidy in Leishmania can be adaptive. Additionally, we believe that the high karyotype diversity in vitro and absence of classical transcriptional regulation make Leishmania an attractive model to study processes of protein homeostasis in the context of aneuploidy and beyond., Competing Interests: The authors have declared that no competing interests exist.

Many fields of biology rely on the inference of accurate multiple sequence alignments (MSA) of biological sequences. Unfortunately, the problem of assembling an MSA is NP-complete thus limiting computation to approximate solutions using heuristics solutions. The progressive algorithm is one of the most popular frameworks for the computation of MSAs. It involves pre-clustering the sequences and aligning them starting with the most similar ones. The scalability of this framework is limited, especially with respect to accuracy. We present here an alternative approach named regressive algorithm. In this framework, sequences are first clustered and then aligned starting with the most distantly related ones. This approach has been shown to greatly improve accuracy during scale-up, especially on datasets featuring 10,000 sequences or more. Another benefit is the possibility to integrate third-party clustering methods and third-party MSA aligners. The regressive algorithm has been tested on up to 1.5 million sequences, its implementation is available in the T-Coffee package.

Many next-generation sequencing datasets contain only relative information because of biological and technical factors that limit the total number of transcripts observed for a given sample. It is not possible to interpret any one component in isolation. The field of compositional data analysis has emerged with alternative methods for relative data based on log-ratio transforms. However, these data often contain many more features than samples, and thus require creative new ways to reduce the dimensionality of the data. The summation of parts, called amalgamation, is a practical way of reducing dimensionality, but can introduce a non-linear distortion to the data. We exploit this non-linearity to propose a powerful yet interpretable dimension method called data-driven amalgamation. Our new method, implemented in the user-friendly R package amalgam, can reduce the dimensionality of compositional data by finding amalgamations that optimally (i) preserve the distance between samples, or (ii) classify samples as diseased or not. Our benchmark on 13 real datasets confirm that these amalgamations compete with state-of-the-art methods in terms of performance, but result in new features that are easily understood: they are groups of parts added together., (© The Author(s) 2019. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

Since the turn of the century, technological advances have made it possible to obtain the molecular profile of any tissue in a cost-effective manner. Among these advances are sophisticated high-throughput assays that measure the relative abundances of microorganisms, RNA molecules, and metabolites. While these data are most often collected to gain new insights into biological systems, they can also be used as biomarkers to create clinically useful diagnostic classifiers. How best to classify high-dimensional -omics data remains an area of active research. However, few explicitly model the relative nature of these data and instead rely on cumbersome normalizations. This report (i) emphasizes the relative nature of health biomarkers, (ii) discusses the literature surrounding the classification of relative data, and (iii) benchmarks how different transformations perform for regularized logistic regression across multiple biomarker types. We show how an interpretable set of log contrasts, called balances, can prepare data for classification. We propose a simple procedure, called discriminative balance analysis, to select groups of 2 and 3 bacteria that can together discriminate between experimental conditions. Discriminative balance analysis is a fast, accurate, and interpretable alternative to data normalization. IMPORTANCE High-throughput sequencing provides an easy and cost-effective way to measure the relative abundance of bacteria in any environmental or biological sample. When these samples come from humans, the microbiome signatures can act as biomarkers for disease prediction. However, because bacterial abundance is measured as a composition, the data have unique properties that make conventional analyses inappropriate. To overcome this, analysts often use cumbersome normalizations. This article proposes an alternative method that identifies pairs and trios of bacteria whose stoichiometric presence can differentiate between diseased and nondiseased samples. By using interpretable log contrasts called balances, we developed an entirely normalization-free classification procedure that reduces the feature space and improves the interpretability, without sacrificing classifier performance., (Copyright © 2020 Quinn and Erb.)

Multiple sequence alignments (MSAs) are used for structural 1,2 and evolutionary predictions 1,2 , but the complexity of aligning large datasets requires the use of approximate solutions 3 , including the progressive algorithm 4 . Progressive MSA methods start by aligning the most similar sequences and subsequently incorporate the remaining sequences, from leaf to root, based on a guide tree. Their accuracy declines substantially as the number of sequences is scaled up 5 . We introduce a regressive algorithm that enables MSA of up to 1.4 million sequences on a standard workstation and substantially improves accuracy on datasets larger than 10,000 sequences. Our regressive algorithm works the other way around from the progressive algorithm and begins by aligning the most dissimilar sequences. It uses an efficient divide-and-conquer strategy to run third-party alignment methods in linear time, regardless of their original complexity. Our approach will enable analyses of extremely large genomic datasets such as the recently announced Earth BioGenome Project, which comprises 1.5 million eukaryotic genomes 6 .

Background: Next-generation sequencing (NGS) has made it possible to determine the sequence and relative abundance of all nucleotides in a biological or environmental sample. A cornerstone of NGS is the quantification of RNA or DNA presence as counts. However, these counts are not counts per se: their magnitude is determined arbitrarily by the sequencing depth, not by the input material. Consequently, counts must undergo normalization prior to use. Conventional normalization methods require a set of assumptions: they assume that the majority of features are unchanged and that all environments under study have the same carrying capacity for nucleotide synthesis. These assumptions are often untestable and may not hold when heterogeneous samples are compared., Results: Methods developed within the field of compositional data analysis offer a general solution that is assumption-free and valid for all data. Herein, we synthesize the extant literature to provide a concise guide on how to apply compositional data analysis to NGS count data., Conclusions: In highlighting the limitations of total library size, effective library size, and spike-in normalizations, we propose the log-ratio transformation as a general solution to answer the question, "Relative to some important activity of the cell, what is changing?", (© The Author(s) 2019. Published by Oxford University Press.)

Obesity is an important health problem with a strong environmental component that is acquiring pandemic proportion. The high availability of caloric dense foods promotes overeating potentially causing obesity. Animal models are key to validate novel therapeutic strategies, but researchers must carefully select the appropriate model to draw the right conclusions. Obesity is defined by an increased body mass index greater than 30 and characterized by an excess of adipose tissue. However, the regulation of food intake involves a close interrelationship between homeostatic and non-homeostatic factors. Studies in animal models have shown that intermittent access to sweetened or calorie-dense foods induces changes in feeding behavior. However, these studies are focused mainly on the final outcome (obesity) rather than on the primary dysfunction underlying the overeating of palatable foods. We describe a protocol to study overeating in mice using diet-induced obesity (DIO). This method can be applied to free choice between palatable food and a standard rodent chow or to forced intake of calorie-dense and/or palatable diets. Exposure to such diets is sufficient to promote changes in meal pattern that we register and analyze during the period of weight gain allowing the longitudinal characterization of feeding behavior in mice. Abnormal eating behaviors such as binge eating or snacking, behavioral alterations commonly observed in obese humans, can be detected using our protocol. In the free-choice procedure, mice develop a preference for the rewarding palatable food showing the reinforcing effect of this diet. Compulsive components of feeding are reflected by maintenance of feeding despite an adverse bitter taste caused by adulteration with quinine and by the negligence of standard chow when access to palatable food is ceased or temporally limited. Our strategy also enables to identify compulsive overeating in mice under a high-caloric regime by using limited food access and finally, we propose complementary behavioral tests to confirm the non-homeostatic food-taking triggered by these foods. Finally, we describe how to computationally explore large longitudinal behavioral datasets., Competing Interests: Competing interestsThe authors have no conflicts of interest to declare., (Copyright © 2019 The Authors; exclusive licensee Bio-protocol LLC.)

We present a new web application to query and visualize time-series behavioral data: the Pergola web-server. This server provides a user-friendly interface for exploring longitudinal behavioral data taking advantage of the Pergola Python library. Using the server, users can process the data applying some basic operations, such as binning or grouping, while formatting the data into existing genomic formats. Thanks to this repurposing of genomics standards, the application automatically renders an interactive data visualization based on sophisticated genome visualization tools. Our tool allows behavioral scientists to share, display and navigate complex behavioral data comprising multiple individuals and multiple data types, in a scalable and flexible manner. A download option allows for further analysis using genomic tools. The server can be a great resource for the field in a time where behavioral science is entering a data-intensive cycle thanks to high-throughput behavioral phenotyping platforms. Pergola is publicly available at http://pergola.crg.eu/., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Autism spectrum disorders are early onset neurodevelopmental disorders characterized by deficits in social communication and restricted repetitive behaviors, yet they are quite heterogeneous in terms of their genetic basis and phenotypic manifestations. Recently, de novo pathogenic mutations in DYRK1A, a chromosome 21 gene associated to neuropathological traits of Down syndrome, have been identified in patients presenting a recognizable syndrome included in the autism spectrum. These mutations produce DYRK1A kinases with partial or complete absence of the catalytic domain, or they represent missense mutations located within this domain. Here, we undertook an extensive biochemical characterization of the DYRK1A missense mutations reported to date and show that most of them, but not all, result in enzymatically dead DYRK1A proteins. We also show that haploinsufficient Dyrk1a +/- mutant mice mirror the neurological traits associated with the human pathology, such as defective social interactions, stereotypic behaviors and epileptic activity. These mutant mice present altered proportions of excitatory and inhibitory neocortical neurons and synapses. Moreover, we provide evidence that alterations in the production of cortical excitatory neurons are contributing to these defects. Indeed, by the end of the neurogenic period, the expression of developmental regulated genes involved in neuron differentiation and/or activity is altered. Therefore, our data indicate that altered neocortical neurogenesis could critically affect the formation of cortical circuits, thereby contributing to the neuropathological changes in DYRK1A haploinsufficiency syndrome., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

The growing appetite of behavioral neuroscience for automated data production is prompting the need for new computational standards allowing improved interoperability, reproducibility, and shareability. We show here how these issues can be solved by repurposing existing genomic formats whose structure perfectly supports the handling of time series. This allows existing genomic analysis and visualization tools to be deployed onto behavioral data. As a proof of principle, we implemented the conversion procedure in Pergola, an open source software, and used genomics tools to reproduce results obtained in mouse, fly, and worm. We also show how common genomics techniques such as principal component analysis, hidden Markov modeling, and volcano plots can be deployed on the reformatted behavioral data. These analyses are easy to share because they depend on the scripting of public software. They are also easy to reproduce thanks to their integration within Nextflow, a workflow manager using containerized software., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

Motivation: Although seldom acknowledged explicitly, count data generated by sequencing platforms exist as compositions for which the abundance of each component (e.g. gene or transcript) is only coherently interpretable relative to other components within that sample. This property arises from the assay technology itself, whereby the number of counts recorded for each sample is constrained by an arbitrary total sum (i.e. library size). Consequently, sequencing data, as compositional data, exist in a non-Euclidean space that, without normalization or transformation, renders invalid many conventional analyses, including distance measures, correlation coefficients and multivariate statistical models., Results: The purpose of this review is to summarize the principles of compositional data analysis (CoDA), provide evidence for why sequencing data are compositional, discuss compositionally valid methods available for analyzing sequencing data, and highlight future directions with regard to this field of study., Supplementary Information: Supplementary data are available at Bioinformatics online.

Obesity represents an important risk factor contributing to the global burden of disease. The current obesogenic environment with easy access to calorie-dense foods is fueling this obesity epidemic. However, how these foods contribute to the progression of feeding behavior changes that lead to overeating is not well understood and needs systematic assessment. Using novel automated methods for the high-throughput screening of behavior, we here examine mice meal pattern upon long-term exposure to a free-choice chocolate-mixture diet and a high-fat diet with face validity for a rapid development of obesity induced by unhealthy food regularly consumed in our societies. We identified rapid diet-specific behavioral changes after exposure to those high-caloric diets. Mice fed with high-fat chow, showed long-lasting meal pattern disturbances, which initiate with a stable loss of circadian feeding rhythmicity. Mice receiving a chocolate-mixture showed qualitatively similar changes, though less marked, consisting in a transient disruption of the feeding behavior and the circadian feeding rhytmicity. Strikingly, compulsive-like eating behavior is triggered immediately after exposure to both high-fat food and chocolate-mixture diet, well before any changes in body weight could be observed. We propose these changes as behavioral biomarkers of prodromal states of obesity that could allow early intervention., (© 2018 Society for the Study of Addiction.)

A major problem in treating obesity is the high rate of relapse to abnormal food-taking habits after maintaining an energy balanced diet. Alterations of eating behavior such as compulsive-like behavior and lack of self-control over food intake play a critical role in relapse. In this study, we used an operant paradigm of food-seeking behavior on two different diet-induced obesity models, a free-choice chocolate-mixture diet and a high-fat diet with face validity for a rapid development of obesity or for unhealthy food regularly consumed in our societies. A reduced operant performance and motivation for the hedonic value of palatable chocolate pellets was revealed in both obesity mouse models. However, only mice exposed to high-fat diet showed an increased compulsive-like behavior in the absence of the reinforcer further characterized by impaired operant learning, enhanced impulsivity and intensified inflexibility. We used principal component analysis to globally identify the specific behaviors responsible for the differences among diet groups. Learning impairment and inflexible behaviors contributed to a first principal component, explaining the largest proportion of the variance in the high-fat diet mice phenotype. Reinforcement, impulsion and compulsion were the main contributors to the second principal component explaining the differences in the chocolate-mixture mice behavioral phenotype. These behaviors were not exclusive of chocolate group because some high-fat individuals showed similar values on this component. These data indicate that extended access to hypercaloric diets differentially modifies operant behavior learning, behavioral flexibility, impulsive-like and compulsive-like behavior, and these effects were dependent on the exposure to each specific diet., (© 2017 Society for the Study of Addiction.)

Background: Genomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction., Results: We generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species., Conclusions: The genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.

Intellectual disability in Down syndrome (DS) is accompanied by altered neuro-architecture, deficient synaptic plasticity, and excitation-inhibition imbalance in critical brain regions for learning and memory. Recently, we have demonstrated beneficial effects of a combined treatment with green tea extract containing (-)-epigallocatechin-3-gallate (EGCG) and cognitive stimulation in young adult DS individuals. Although we could reproduce the cognitive-enhancing effects in mouse models, the underlying mechanisms of these beneficial effects are unknown. Here, we explored the effects of a combined therapy with environmental enrichment (EE) and EGCG in the Ts65Dn mouse model of DS at young age. Our results show that combined EE-EGCG treatment improved corticohippocampal-dependent learning and memory. Cognitive improvements were accompanied by a rescue of cornu ammonis 1 (CA1) dendritic spine density and a normalization of the proportion of excitatory and inhibitory synaptic markers in CA1 and dentate gyrus., Competing Interests: Authors report no conflict of interest.

This review provides an overview on the development of Multiple sequence alignment (MSA) methods and their main applications. It is focused on progress made over the past decade. The three first sections review recent algorithmic developments for protein, RNA/DNA and genomic alignments. The fourth section deals with benchmarks and explores the relationship between empirical and simulated data, along with the impact on method developments. The last part of the review gives an overview on available MSA local reliability estimators and their dependence on various algorithmic properties of available methods., (© The Author 2015. Published by Oxford University Press. For Permissions, please email: journals.permissions@oup.com.)

Correlation is ubiquitously used in gene expression analysis although its validity as an objective criterion is often questionable. If no normalization reflecting the original mRNA counts in the cells is available, correlation between genes becomes spurious. Yet the need for normalization can be bypassed using a relative analysis approach called log-ratio analysis. This approach can be used to identify proportional gene pairs, i.e. a subset of pairs whose correlation can be inferred correctly from unnormalized data due to their vanishing log-ratio variance. To interpret the size of non-zero log-ratio variances, a proposal for a scaling with respect to the variance of one member of the gene pair was recently made by Lovell et al. Here we derive analytically how spurious proportionality is introduced when using a scaling. We base our analysis on a symmetric proportionality coefficient (briefly mentioned in Lovell et al.) that has a number of advantages over their statistic. We show in detail how the choice of reference needed for the scaling determines which gene pairs are identified as proportional. We demonstrate that using an unchanged gene as a reference has huge advantages in terms of sensitivity. We also explore the link between proportionality and partial correlation and derive expressions for a partial proportionality coefficient. A brief data-analysis part puts the discussed concepts into practice.

Background: Legumes are the third largest family of angiosperms and the second most important crop class. Legume genomes have been shaped by extensive large-scale gene duplications, including an approximately 58 million year old whole genome duplication shared by most crop legumes., Results: We report the genome and the transcription atlas of coding and non-coding genes of a Mesoamerican genotype of common bean (Phaseolus vulgaris L., BAT93). Using a comprehensive phylogenomics analysis, we assessed the past and recent evolution of common bean, and traced the diversification of patterns of gene expression following duplication. We find that successive rounds of gene duplications in legumes have shaped tissue and developmental expression, leading to increased levels of specialization in larger gene families. We also find that many long non-coding RNAs are preferentially expressed in germ-line-related tissues (pods and seeds), suggesting that they play a significant role in fruit development. Our results also suggest that most bean-specific gene family expansions, including resistance gene clusters, predate the split of the Mesoamerican and Andean gene pools., Conclusions: The genome and transcriptome data herein generated for a Mesoamerican genotype represent a counterpart to the genomic resources already available for the Andean gene pool. Altogether, this information will allow the genetic dissection of the characters involved in the domestication and adaptation of the crop, and their further implementation in breeding strategies for this important crop.

Down syndrome (DS) individuals present increased risk for Alzheimer's disease (AD) neuropathology and AD-type dementia. Here, we investigated the use of green tea extracts containing (-)-epigallocatechin-3-gallate (EGCG), as co-adjuvant to enhance the effects of environmental enrichment (EE) in Ts65Dn mice, a segmental trisomy model of DS that partially mimics DS/AD pathology, at the age of initiation of cognitive decline. Classical repeated measures ANOVA showed that combined EE-EGCG treatment was more efficient than EE or EGCG alone to improve specific spatial learning related variables. Using principal component analysis (PCA) we found that several spatial learning parameters contributed similarly to a first PC and explained a large proportion of the variance among groups, thus representing a composite learning measure. This PC1 revealed that EGCG or EE alone had no significant effect. However, combined EE-EGCG significantly ameliorated learning alterations of middle age Ts65Dn mice. Interestingly, PCA revealed an increased variability along learning sessions with good and poor learners in Ts65Dn, and this stratification did not disappear upon treatments. Our results suggest that combining EE and EGCG represents a viable therapeutic approach for amelioration of age-related cognitive decline in DS, although its efficacy may vary across individuals.

Multiple sequence alignments (MSAs) are a prerequisite for a wide variety of evolutionary analyses. Published assessments and benchmark data sets for protein and, to a lesser extent, global nucleotide MSAs are available, but less effort has been made to establish benchmarks in the more general problem of whole-genome alignment (WGA). Using the same model as the successful Assemblathon competitions, we organized a competitive evaluation in which teams submitted their alignments and then assessments were performed collectively after all the submissions were received. Three data sets were used: Two were simulated and based on primate and mammalian phylogenies, and one was comprised of 20 real fly genomes. In total, 35 submissions were assessed, submitted by 10 teams using 12 different alignment pipelines. We found agreement between independent simulation-based and statistical assessments, indicating that there are substantial accuracy differences between contemporary alignment tools. We saw considerable differences in the alignment quality of differently annotated regions and found that few tools aligned the duplications analyzed. We found that many tools worked well at shorter evolutionary distances, but fewer performed competitively at longer distances. We provide all data sets, submissions, and assessment programs for further study and provide, as a resource for future benchmarking, a convenient repository of code and data for reproducing the simulation assessments., (© 2014 Earl et al.; Published by Cold Spring Harbor Laboratory Press.)

T-Coffee, for Tree-based consistency objective function for alignment evaluation, is a versatile multiple sequence alignment (MSA) method suitable for aligning virtually any type of biological sequences. T-Coffee provides more than a simple sequence aligner; rather it is a framework in which alternative alignment methods and/or extra information (i.e., structural, evolutionary, or experimental information) can be combined to reach more accurate and more meaningful MSAs. T-Coffee can be used either by running input data via the Web server ( http://tcoffee.crg.cat/apps/tcoffee/index.html ) or by downloading the T-Coffee package. Here, we present how the package can be used in its command line mode to carry out the most common tasks and multiply align proteins, DNA, and RNA sequences. This chapter particularly emphasizes on the description of T-Coffee special flavors also called "modes," designed to address particular biological problems.

Predicting protein functional classes such as localization sites and modifications plays a crucial role in function annotation. Given a tremendous amount of sequence data yielded from high-throughput sequencing experiments, the need of efficient and interpretable prediction strategies has been rapidly amplified. Our previous approach for subcellular localization prediction, PSLDoc, archives high overall accuracy for Gram-negative bacteria. However, PSLDoc is computational intensive due to incorporation of homology extension in feature extraction and probabilistic latent semantic analysis in feature reduction. Besides, prediction results generated by support vector machines are accurate but generally difficult to interpret. In this work, we incorporate three new techniques to improve efficiency and interpretability. First, homology extension is performed against a compact non-redundant database using a fast search model to reduce running time. Second, correspondence analysis (CA) is incorporated as an efficient feature reduction to generate a clear visual separation of different protein classes. Finally, functional classes are predicted by a combination of accurate compact set (CS) relation and interpretable one-nearest neighbor (1-NN) algorithm. Besides localization data sets, we also apply a human protein kinase set to validate generality of our proposed method. Experiment results demonstrate that our method make accurate prediction in a more efficient and interpretable manner. First, homology extension using a fast search on a compact database can greatly accelerate traditional running time up to twenty-five times faster without sacrificing prediction performance. This suggests that computational costs of many other predictors that also incorporate homology information can be largely reduced. In addition, CA can not only efficiently identify discriminative features but also provide a clear visualization of different functional classes. Moreover, predictions based on CS achieve 100% precision. When combined with 1-NN on unpredicted targets by CS, our method attains slightly better or comparable performance compared with the state-of-the-art systems.

We address the challenge of regulatory sequence alignment with a new method, Pro-Coffee, a multiple aligner specifically designed for homologous promoter regions. Pro-Coffee uses a dinucleotide substitution matrix estimated on alignments of functional binding sites from TRANSFAC. We designed a validation framework using several thousand families of orthologous promoters. This dataset was used to evaluate the accuracy for predicting true human orthologs among their paralogs. We found that whereas other methods achieve on average 73.5% accuracy, and 77.6% when trained on that same dataset, the figure goes up to 80.4% for Pro-Coffee. We then applied a novel validation procedure based on multi-species ChIP-seq data. Trained and untrained methods were tested for their capacity to correctly align experimentally detected binding sites. Whereas the average number of correctly aligned sites for two transcription factors is 284 for default methods and 316 for trained methods, Pro-Coffee achieves 331, 16.5% above the default average. We find a high correlation between a method's performance when classifying orthologs and its ability to correctly align proven binding sites. Not only has this interesting biological consequences, it also allows us to conclude that any method that is trained on the ortholog data set will result in functionally more informative alignments.

Motivation: Probabilistic approaches for inferring transcription factor binding sites (TFBSs) and regulatory motifs from DNA sequences have been developed for over two decades. Previous work has shown that prediction accuracy can be significantly improved by incorporating features such as the competition of multiple transcription factors (TFs) for binding to nearby sites, the tendency of TFBSs for co-regulated TFs to cluster and form cis-regulatory modules and explicit evolutionary modeling of conservation of TFBSs across orthologous sequences. However, currently available tools only incorporate some of these features, and significant methodological hurdles hampered their synthesis into a single consistent probabilistic framework., Results: We present MotEvo, a integrated suite of Bayesian probabilistic methods for the prediction of TFBSs and inference of regulatory motifs from multiple alignments of phylogenetically related DNA sequences, which incorporates all features just mentioned. In addition, MotEvo incorporates a novel model for detecting unknown functional elements that are under evolutionary constraint, and a new robust model for treating gain and loss of TFBSs along a phylogeny. Rigorous benchmarking tests on ChIP-seq datasets show that MotEvo's novel features significantly improve the accuracy of TFBS prediction, motif inference and enhancer prediction., Availability: Source code, a user manual and files with several example applications are available at www.swissregulon.unibas.ch.

T-Coffee (Tree-based consistency objective function for alignment evaluation) is a versatile multiple sequence alignment (MSA) method suitable for aligning most types of biological sequences. The main strength of T-Coffee is its ability to combine third party aligners and to integrate structural (or homology) information when building MSAs. The series of protocols presented here show how the package can be used to multiply align proteins, RNA and DNA sequences. The protein section shows how users can select the most suitable T-Coffee mode for their data set. Detailed protocols include T-Coffee, the default mode, M-Coffee, a meta version able to combine several third party aligners into one, PSI (position-specific iterated)-Coffee, the homology extended mode suitable for remote homologs and Expresso, the structure-based multiple aligner. We then also show how the T-RMSD (tree based on root mean square deviation) option can be used to produce a functionally informative structure-based clustering. RNA alignment procedures are described for using R-Coffee, a mode able to use predicted RNA secondary structures when aligning RNA sequences. DNA alignments are illustrated with Pro-Coffee, a multiple aligner specific of promoter regions. We also present some of the many reformatting utilities bundled with T-Coffee. The package is an open-source freeware available from http://www.tcoffee.org/.

We present and validate BlastR, a method for efficiently and accurately searching non-coding RNAs. Our approach relies on the comparison of di-nucleotides using BlosumR, a new log-odd substitution matrix. In order to use BlosumR for comparison, we recoded RNA sequences into protein-like sequences. We then showed that BlosumR can be used along with the BlastP algorithm in order to search non-coding RNA sequences. Using Rfam as a gold standard, we benchmarked this approach and show BlastR to be more sensitive than BlastN. We also show that BlastR is both faster and more sensitive than BlastP used with a single nucleotide log-odd substitution matrix. BlastR, when used in combination with WU-BlastP, is about 5% more accurate than WU-BlastN and about 50 times slower. The approach shown here is equally effective when combined with the NCBI-Blast package. The software is an open source freeware available from www.tcoffee.org/blastr.html.

The availability of sequence specificities for a substantial fraction of yeast's transcription factors and comparative genomic algorithms for binding site prediction has made it possible to comprehensively annotate transcription factor binding sites genome-wide. Here we use such a genome-wide annotation for comprehensively studying promoter architecture in yeast, focusing on the distribution of transcription factor binding sites relative to transcription start sites, and the architecture of TATA and TATA-less promoters. For most transcription factors, binding sites are positioned further upstream and vary over a wider range in TATA promoters than in TATA-less promoters. In contrast, a group of 6 'proximal promoter motifs' (GAT1/GLN3/DAL80, FKH1/2, PBF1/2, RPN4, NDT80, and ROX1) occur preferentially in TATA-less promoters and show a strong preference for binding close to the transcription start site in these promoters. We provide evidence that suggests that pre-initiation complexes are recruited at TATA sites in TATA promoters and at the sites of the other proximal promoter motifs in TATA-less promoters. TATA-less promoters can generally be classified by the proximal promoter motif they contain, with different classes of TATA-less promoters showing different patterns of transcription factor binding site positioning and nucleosome coverage. These observations suggest that different modes of regulation of transcription initiation may be operating in the different promoter classes. In addition we show that, across all promoter classes, there is a close match between nucleosome free regions and regions of highest transcription factor binding site density. This close agreement between transcription factor binding site density and nucleosome depletion suggests a direct and general competition between transcription factors and nucleosomes for binding to promoters.

Budding yeast telomeres and cryptic mating-type loci are enriched at the nuclear envelope, forming foci that sequester silent information regulators (SIR factors), much as heterochromatic chromocenters in higher eukaryotes sequester HP1. Here we examine the impact of such subcompartments for regulating transcription genome-wide. We show that the efficiency of subtelomeric reporter gene repression depends not only on the strength of SIR factor recruitment by cis-acting elements, but also on the accumulation of SIRs in such perinuclear foci. To monitor the effects of disrupting this subnuclear compartment, we performed microarray analyses under conditions that eliminate telomere anchoring, while preserving SIR complex integrity. We found 60 genes reproducibly misregulated. Among those with increased expression, 22% were within 20 kb of a telomere, confirming that the nuclear envelope (NE) association of telomeres helps repress natural subtelomeric genes. In contrast, loci that were down-regulated were distributed over all chromosomes. Half of this ectopic repression was SIR complex dependent. We conclude that released SIR factors can promiscuously repress transcription at nontelomeric genes despite the presence of "anti-silencing" mechanisms. Bioinformatic analysis revealed that promoters bearing the PAC (RNA Polymerase A and C promoters) or Abf1 binding consenses are consistently down-regulated by mislocalization of SIR factors. Thus, the normal telomeric sequestration of SIRs both favors subtelomeric repression and prevents promiscuous effects at a distinct subset of promoters. This demonstrates that patterns of gene expression can be regulated by changing the spatial distribution of repetitive DNA sequences that bind repressive factors.

The autonomously replicating sequence binding factor 1 (Abf1) was initially identified as an essential DNA replication factor and later shown to be a component of the regulatory network controlling mitotic and meiotic cell cycle progression in budding yeast. The protein is thought to exert its functions via specific interaction with its target site as part of distinct protein complexes, but its roles during mitotic growth and meiotic development are only partially understood. Here, we report a comprehensive approach aiming at the identification of direct Abf1-target genes expressed during fermentation, respiration, and sporulation. Computational prediction of the protein's target sites was integrated with a genome-wide DNA binding assay in growing and sporulating cells. The resulting data were combined with the output of expression profiling studies using wild-type versus temperature-sensitive alleles. This work identified 434 protein-coding loci as being transcriptionally dependent on Abf1. More than 60% of their putative promoter regions contained a computationally predicted Abf1 binding site and/or were bound by Abf1 in vivo, identifying them as direct targets. The present study revealed numerous loci previously unknown to be under Abf1 control, and it yielded evidence for the protein's variable DNA binding pattern during mitotic growth and meiotic development.