02277nas a2200193 4500008004100000022001400041245011200055210006900167260000900236300001100245490000600256520164600262100002401908700002101932700002101953700002001974700002001994856006902014 2012 eng d a1176-934300aEvolutionary Genomics of Genes Involved in Olfactory Behavior in the Drosophila melanogaster Species Group.0 aEvolutionary Genomics of Genes Involved in Olfactory Behavior in c2012 a89-1040 v83 aPrevious comparative genomic studies of genes involved in olfactory behavior in Drosophila focused only on particular gene families such as odorant receptor and/or odorant binding proteins. However, olfactory behavior has a complex genetic architecture that is orchestrated by many interacting genes. In this paper, we present a comparative genomic study of olfactory behavior in Drosophila including an extended set of genes known to affect olfactory behavior. We took advantage of the recent burst of whole genome sequences and the development of powerful statistical tools to analyze genomic data and test evolutionary and functional hypotheses of olfactory genes in the six species of the Drosophila melanogaster species group for which whole genome sequences are available. Our study reveals widespread purifying selection and limited incidence of positive selection on olfactory genes. We show that the pace of evolution of olfactory genes is mostly independent of the life cycle stage, and of the number of life cycle stages, in which they participate in olfaction. However, we detected a relationship between evolutionary rates and the position that the gene products occupy in the olfactory system, genes occupying central positions tend to be more constrained than peripheral genes. Finally, we demonstrate that specialization to one host does not seem to be associated with bursts of adaptive evolution in olfactory genes in D. sechellia and D. erecta, the two specialists species analyzed, but rather different lineages have idiosyncratic evolutionary histories in which both historical and ecological factors have been involved.1 aLavagnino, Nicolás1 aSerra, François1 aArbiza, Leonardo1 aDopazo, Hernán1 aHasson, Esteban uhttp://www.ncbi.nlm.nih.gov/pmc/articles/PMC3273929/?tool=pubmed01377nas a2200313 4500008004100000245005700041210005500098260001300153300001300166490000600179520044900185100001700634700002200651700002700673700002200700700002000722700002200742700001800764700001900782700001700801700003300818700002800851700002000879700002600899700002500925700001600950700002100966856007600987 2011 eng d00aDiscovery of an ebolavirus-like filovirus in europe.0 aDiscovery of an ebolaviruslike filovirus in europe c2011 Oct ae10023040 v73 a

Filoviruses, amongst the most lethal of primate pathogens, have only been reported as natural infections in sub-Saharan Africa and the Philippines. Infections of bats with the ebolaviruses and marburgviruses do not appear to be associated with disease. Here we report identification in dead insectivorous bats of a genetically distinct filovirus, provisionally named Lloviu virus, after the site of detection, Cueva del Lloviu, in Spain.

1 aNegredo, Ana1 aPalacios, Gustavo1 aVázquez-Morón, Sonia1 aGonzález, Félix1 aDopazo, Hernán1 aMolero, Francisca1 aJuste, Javier1 aQuetglas, Juan1 aSavji, Nazir1 aMartínez, Maria, de la Cruz1 aHerrera, Jesus, Enrique1 aPizarro, Manuel1 aHutchison, Stephen, K1 aEchevarría, Juan, E1 aLipkin, Ian1 aTenorio, Antonio uhttp://clinbioinfosspa.es/content/discovery-ebolavirus-filovirus-europe02093nas a2200169 4500008004100000245011500041210006900156260001600225520142700241100002501668700001701693700002101710700002401731700002001755700001901775856012901794 2011 eng d00aEvolution of the biosynthesis of di-myo-inositol phosphate, a marker of adaptation to hot marine environments.0 aEvolution of the biosynthesis of dimyoinositol phosphate a marke c2011 Oct 263 a

The synthesis of di-myo-inositol phosphate (DIP), a common compatible solute in hyperthermophiles, involves the consecutive actions of inositol-1-phosphate cytidylyltransferase (IPCT) and di-myo-inositol phosphate phosphate synthase (DIPPS). In most cases, both activities are present in a single gene product, but separate genes are also found in a few organisms. Genes for IPCT and DIPPS were found in the genomes of 33 organisms, all with thermophilic/hyperthermophilic lifestyles. Phylogeny of IPCT/DIPPS revealed an incongruent topology with 16S RNA phylogeny, thus suggesting horizontal gene transfer. The phylogenetic tree of the DIPPS domain was rooted by using phosphatidylinositol phosphate synthase sequences as out-group. The root locates at the separation of genomes with fused and split genes. We propose that the gene encoding DIPPS was recruited from the biosynthesis of phosphatidylinositol. The last DIP-synthesizing ancestor harboured separated genes for IPCT and DIPPS and this architecture was maintained in a crenarchaeal lineage, and transferred by horizontal gene transfer to hyperthermophilic marine Thermotoga species. It is plausible that the driving force for the assembly of those two genes in the early ancestor is related to the acquired advantage of DIP producers to cope with high temperature. This work corroborates the view that Archaea were the first hyperthermophilic organisms.

1 aGonçalves, Luís, G1 aBorges, Nuno1 aSerra, François1 aFernandes, Pedro, L1 aDopazo, Hernán1 aSantos, Helena uhttp://clinbioinfosspa.es/content/evolution-biosynthesis-di-myo-inositol-phosphate-marker-adaptation-hot-marine-environments02524nas a2200169 4500008004100000245006800041210006600109260000900175300001200184490000600196520197400202100002102176700002102197700002002218700001802238856009802256 2011 eng d00aGenome-wide heterogeneity of nucleotide substitution model fit.0 aGenomewide heterogeneity of nucleotide substitution model fit c2011 a896-9080 v33 a

At a genomic scale, the patterns that have shaped molecular evolution are believed to be largely heterogeneous. Consequently, comparative analyses should use appropriate probabilistic substitution models that capture the main features under which different genomic regions have evolved. While efforts have concentrated in the development and understanding of model selection techniques, no descriptions of overall relative substitution model fit at the genome level have been reported. Here, we provide a characterization of best-fit substitution models across three genomic data sets including coding regions from mammals, vertebrates, and Drosophila (24,000 alignments). According to the Akaike Information Criterion (AIC), 82 of 88 models considered were selected as best-fit models at least in one occasion, although with very different frequencies. Most parameter estimates also varied broadly among genes. Patterns found for vertebrates and Drosophila were quite similar and often more complex than those found in mammals. Phylogenetic trees derived from models in the 95% confidence interval set showed much less variance and were significantly closer to the tree estimated under the best-fit model than trees derived from models outside this interval. Although alternative criteria selected simpler models than the AIC, they suggested similar patterns. All together our results show that at a genomic scale, different gene alignments for the same set of taxa are best explained by a large variety of different substitution models and that model choice has implications on different parameter estimates including the inferred phylogenetic trees. After taking into account the differences related to sample size, our results suggest a noticeable diversity in the underlying evolutionary process. All together, we conclude that the use of model selection techniques is important to obtain consistent phylogenetic estimates from real data at a genomic scale.

1 aArbiza, Leonardo1 aPatricio, Mateus1 aDopazo, Hernán1 aPosada, David uhttp://clinbioinfosspa.es/content/genome-wide-heterogeneity-nucleotide-substitution-model-fit02573nas a2200301 4500008004100000022001400041245006900055210006600124260001300190300001300203490000600216520167900222653001201901653002301913653001501936653001901951653003401970653001302004653001202017653001402029653002302043653002702066100002102093700002102114700002002135700002002155856009602175 2011 eng d a1553-735800aNatural selection on functional modules, a genome-wide analysis.0 aNatural selection on functional modules a genomewide analysis c2011 Mar ae10010930 v73 a

Classically, the functional consequences of natural selection over genomes have been analyzed as the compound effects of individual genes. The current paradigm for large-scale analysis of adaptation is based on the observed significant deviations of rates of individual genes from neutral evolutionary expectation. This approach, which assumed independence among genes, has not been able to identify biological functions significantly enriched in positively selected genes in individual species. Alternatively, pooling related species has enhanced the search for signatures of selection. However, grouping signatures does not allow testing for adaptive differences between species. Here we introduce the Gene-Set Selection Analysis (GSSA), a new genome-wide approach to test for evidences of natural selection on functional modules. GSSA is able to detect lineage specific evolutionary rate changes in a notable number of functional modules. For example, in nine mammal and Drosophilae genomes GSSA identifies hundreds of functional modules with significant associations to high and low rates of evolution. Many of the detected functional modules with high evolutionary rates have been previously identified as biological functions under positive selection. Notably, GSSA identifies conserved functional modules with many positively selected genes, which questions whether they are exclusively selected for fitting genomes to environmental changes. Our results agree with previous studies suggesting that adaptation requires positive selection, but not every mutation under positive selection contributes to the adaptive dynamical process of the evolution of species.

10aAnimals10aDatabases, Genetic10aDrosophila10aGenome, Insect10aGenome-Wide Association Study10aGenomics10aMammals10aPhylogeny10aSelection, Genetic10aSequence Analysis, DNA1 aSerra, François1 aArbiza, Leonardo1 aDopazo, Joaquin1 aDopazo, Hernán uhttp://clinbioinfosspa.es/content/natural-selection-functional-modules-genome-wide-analysis02020nas a2200349 4500008004100000022001400041245011700055210006900172260001300241300001100254490000700265520090400272653002501176653001301201653001301214653001401227653002301241653001301264100002101277700002101298700002301319700002001342700002101362700001701383700002401400700003301424700002401457700002001481700002001501700002001521856012901541 2011 eng d a1362-496200aPhylemon 2.0: a suite of web-tools for molecular evolution, phylogenetics, phylogenomics and hypotheses testing.0 aPhylemon 20 a suite of webtools for molecular evolution phylogen c2011 Jul aW470-40 v393 a

Phylemon 2.0 is a new release of the suite of web tools for molecular evolution, phylogenetics, phylogenomics and hypotheses testing. It has been designed as a response to the increasing demand of molecular sequence analyses for experts and non-expert users. Phylemon 2.0 has several unique features that differentiates it from other similar web resources: (i) it offers an integrated environment that enables evolutionary analyses, format conversion, file storage and edition of results; (ii) it suggests further analyses, thereby guiding the users through the web server; and (iii) it allows users to design and save phylogenetic pipelines to be used over multiple genes (phylogenomics). Altogether, Phylemon 2.0 integrates a suite of 30 tools covering sequence alignment reconstruction and trimming; tree reconstruction, visualization and manipulation; and evolutionary hypotheses testing.

10aEvolution, Molecular10aGenomics10aInternet10aPhylogeny10aSequence Alignment10aSoftware1 aSánchez, Rubén1 aSerra, François1 aTárraga, Joaquín1 aMedina, Ignacio1 aCarbonell, José1 aPulido, Luis1 aDe Maria, Alejandro1 aCapella-Gutíerrez, Salvador1 aHuerta-Cepas, Jaime1 aGabaldón, Toni1 aDopazo, Joaquin1 aDopazo, Hernán uhttp://clinbioinfosspa.es/content/phylemon-20-suite-web-tools-molecular-evolution-phylogenetics-phylogenomics-and-hypotheses01509nas a2200169 4500008004100000245009800041210006900139260001300208300001100221490000700232520088300239100002201122700003201144700002001176700002201196856012101218 2011 eng d00aPhylogenetic and in silico structural analysis of the Parkinson disease-related kinase PINK1.0 aPhylogenetic and in silico structural analysis of the Parkinson  c2011 Apr a369-780 v323 a

Parkinson disease (PD) is the second most common neurodegenerative disorder and is characterized by the loss of dopaminergic neurons in the substantia nigra. Mutations in PINK1 were shown to cause recessive familial PD, and today are proposed to be associated with the disease via mitochondrial dysfunction and oxidative damage. The PINK1 gene comprises eight exons, which encode a ubiquitously expressed 581 amino acid protein that contains an N-terminal mitochondrial targeting domain and a serine/threonine protein kinase. To better understand the relationship between PINK1 and PD we have first analyzed the evolutionary history of the gene showing its late emergence in evolution. In addition, we have modeled the three-dimensional structure of PINK1 and found some evidences that help to explain the effect of some PD-related mutations in this protein’s function.

1 aCardona, Fernando1 aSánchez-Mut, Jose, Vicente1 aDopazo, Hernán1 aPérez-Tur, Jordi uhttp://clinbioinfosspa.es/content/phylogenetic-and-silico-structural-analysis-parkinson-disease-related-kinase-pink101456nas a2200253 4500008004100000245009300041210006900134260000900203300000700212490000700219520063100226100002300857700001500880700001900895700002300914700002800937700001500965700001700980700002100997700002401018700002001042700002001062856012001082 2011 eng d00aRecent human evolution has shaped geographical differences in susceptibility to disease.0 aRecent human evolution has shaped geographical differences in su c2011 a550 v123 a

Searching for associations between genetic variants and complex diseases has been a very active area of research for over two decades. More than 51,000 potential associations have been studied and published, a figure that keeps increasing, especially with the recent explosion of array-based Genome-Wide Association Studies. Even if the number of true associations described so far is high, many of the putative risk variants detected so far have failed to be consistently replicated and are widely considered false positives. Here, we focus on the world-wide patterns of replicability of published association studies.

1 aMarigorta, Urko, M1 aLao, Oscar1 aCasals, Ferran1 aCalafell, Francesc1 aMorcillo-Suarez, Carlos1 aFaria, Rui1 aBosch, Elena1 aSerra, François1 aBertranpetit, Jaume1 aDopazo, Hernán1 aNavarro, Arcadi uhttp://clinbioinfosspa.es/content/recent-human-evolution-has-shaped-geographical-differences-susceptibility-disease02322nas a2200241 4500008004100000245007500041210006900116260001300185300001100198490000700209520154000216100002401756700003301780700002101813700001901834700003501853700002001888700001801908700002701926700002001953700001701973856009001990 2011 eng d00aRole of tomato BRANCHED1-like genes in the control of shoot branching.0 aRole of tomato BRANCHED1like genes in the control of shoot branc c2011 Aug a701-140 v673 a

In angiosperms, shoot branching greatly determines overall plant architecture and affects fundamental aspects of plant life. Branching patterns are determined by genetic pathways conserved widely across angiosperms. In Arabidopsis thaliana (Brassicaceae, Rosidae) BRANCHED1 (BRC1) plays a central role in this process, acting locally to arrest axillary bud growth. In tomato (Solanum lycopersicum, Solanaceae, Asteridae) we have identified two BRC1-like paralogues, SlBRC1a and SlBRC1b. These genes are expressed in arrested axillary buds and both are down-regulated upon bud activation, although SlBRC1a is transcribed at much lower levels than SlBRC1b. Alternative splicing of SlBRC1a renders two transcripts that encode two BRC1-like proteins with different C-t domains due to a 3’-terminal frameshift. The phenotype of loss-of-function lines suggests that SlBRC1b has retained the ancestral role of BRC1 in shoot branch suppression. We have isolated the BRC1a and BRC1b genes of other Solanum species and have studied their evolution rates across the lineages. These studies indicate that, after duplication of an ancestral BRC1-like gene, BRC1b genes continued to evolve under a strong purifying selection that was consistent with the conserved function of SlBRC1b in shoot branching control. In contrast, the coding sequences of Solanum BRC1a genes have evolved at a higher evolution rate. Branch-site tests indicate that this difference does not reflect relaxation but rather positive selective pressure for adaptation.

1 aMartín-Trillo, Mar1 aGrandío, Eduardo, González1 aSerra, François1 aMarcel, Fabien1 aRodríguez-Buey, María, Luisa1 aSchmitz, Gregor1 aTheres, Klaus1 aBendahmane, Abdelhafid1 aDopazo, Hernán1 aCubas, Pilar uhttp://clinbioinfosspa.es/content/role-tomato-branched1-genes-control-shoot-branching01980nas a2200217 4500008004100000022001400041245007200055210006900127260000900196300001100205490000600216520130500222100001601527700002001543700002101563700002901584700002001613700002501633700002501658856007901683 2011 eng d a1932-620300aSexual selection halts the relaxation of protamine 2 among rodents.0 aSexual selection halts the relaxation of protamine 2 among roden c2011 ae292470 v63 aSexual selection has been proposed as the driving force promoting the rapid evolutionary changes observed in some reproductive genes including protamines. We test this hypothesis in a group of rodents which show marked differences in the intensity of sexual selection. Levels of sperm competition were not associated with the evolutionary rates of protamine 1 but, contrary to expectations, were negatively related to the evolutionary rate of cleaved- and mature-protamine 2. Since both domains were found to be under relaxation, our findings reveal an unforeseen role of sexual selection: to halt the degree of degeneration that proteins within families may experience due to functional redundancy. The degree of relaxation of protamine 2 in this group of rodents is such that in some species it has become dysfunctional and it is not expressed in mature spermatozoa. In contrast, protamine 1 is functionally conserved but shows directed positive selection on specific sites which are functionally relevant such as DNA-anchoring domains and phosphorylation sites. We conclude that in rodents protamine 2 is under relaxation and that sexual selection removes deleterious mutations among species with high levels of sperm competition to maintain the protein functional and the spermatozoa competitive.1 aLüke, Lena1 aVicens, Alberto1 aSerra, François1 aLuque-Larena, Juan, Jose1 aDopazo, Hernán1 aRoldan, Eduardo, R S1 aGomendio, Montserrat uhttp://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.002924701951nas a2200205 4500008004100000245007400041210006900115260001500184300001300199490000700212520125700219100002601476700003001502700002501532700002001557700001601577700002101593700003101614856010001645 2011 eng d00aSUS1 introns are required for efficient mRNA nuclear export in yeast.0 aSUS1 introns are required for efficient mRNA nuclear export in y c2011 Oct 1 a8599-6110 v393 a

Efficient coupling between mRNA synthesis and export is essential for gene expression. Sus1/ENY2, a component of the SAGA and TREX-2 complexes, is involved in both transcription and mRNA export. While most yeast genes lack introns, we previously reported that yeast SUS1 bears two. Here we show that this feature is evolutionarily conserved and critical for Sus1 function. We determine that while SUS1 splicing is inefficient, it responds to cellular conditions, and intronic mutations either promoting or blocking splicing lead to defects in mRNA export and cell growth. Consistent with this, we find that an intron-less SUS1 only partially rescues sus1Δ phenotypes. Remarkably, splicing of each SUS1 intron is also affected by the presence of the other and by SUS1 exonic sequences. Moreover, by following SUS1 RNA and protein levels we establish that nonsense-mediated decay (NMD) pathway and the splicing factor Mud2 both play a role in SUS1 expression. Our data (and those of the accompanying work by Hossain et al.) provide evidence of the involvement of splicing, translation, and decay in the regulation of early events in mRNP biogenesis; and imply the additional requirement for a balance in splicing isoforms from a single gene.

1 aCuenca-Bono, Bernardo1 aGarcía-Molinero, Varinia1 aPascual-García, Pau1 aDopazo, Hernán1 aLlopis, Ana1 aVilardell, Josep1 aRodríguez-Navarro, Susana uhttp://clinbioinfosspa.es/content/sus1-introns-are-required-efficient-mrna-nuclear-export-yeast03089nas a2200385 4500008004100000022001400041245013400055210006900189260001300258300001200271490000700283520182800290653001502118653001002133653002602143653002302169653002802192653002502220653003802245653002202283653001802305653001102323653001802334653001902352653003602371653001302407653003302420100002202453700002102475700001802496700002002514700002002534700002502554856012402579 2008 eng d a1098-100400aUse of estimated evolutionary strength at the codon level improves the prediction of disease-related protein mutations in humans.0 aUse of estimated evolutionary strength at the codon level improv c2008 Jan a198-2040 v293 a

Predicting the functional impact of protein variation is one of the most challenging problems in bioinformatics. A rapidly growing number of genome-scale studies provide large amounts of experimental data, allowing the application of rigorous statistical approaches for predicting whether a given single point mutation has an impact on human health. Up until now, existing methods have limited their source data to either protein or gene information. Novel in this work, we take advantage of both and focus on protein evolutionary information by using estimated selective pressures at the codon level. Here we introduce a new method (SeqProfCod) to predict the likelihood that a given protein variant is associated with human disease or not. Our method relies on a support vector machine (SVM) classifier trained using three sources of information: protein sequence, multiple protein sequence alignments, and the estimation of selective pressure at the codon level. SeqProfCod has been benchmarked with a large dataset of 8,987 single point mutations from 1,434 human proteins from SWISS-PROT. It achieves 82% overall accuracy and a correlation coefficient of 0.59, indicating that the estimation of the selective pressure helps in predicting the functional impact of single-point mutations. Moreover, this study demonstrates the synergic effect of combining two sources of information for predicting the functional effects of protein variants: protein sequence/profile-based information and the evolutionary estimation of the selective pressures at the codon level. The results of large-scale application of SeqProfCod over all annotated point mutations in SWISS-PROT (available for download at http://sgu.bioinfo.cipf.es/services/Omidios/; last accessed: 24 August 2007), could be used to support clinical studies.

10aAlgorithms10aCodon10aComputational Biology10aDatabases, Protein10aDNA Mutational Analysis10aEvolution, Molecular10aGenetic Predisposition to Disease10aGenetic Variation10aGenome, Human10aHumans10aIduronic Acid10aPoint Mutation10aPolymorphism, Single Nucleotide10aProteins10aTumor Suppressor Protein p531 aCapriotti, Emidio1 aArbiza, Leonardo1 aCasadio, Rita1 aDopazo, Joaquin1 aDopazo, Hernán1 aMarti-Renom, Marc, A uhttp://clinbioinfosspa.es/content/use-estimated-evolutionary-strength-codon-level-improves-prediction-disease-related-002514nas a2200337 4500008004100000022001400041245007300055210006900128260001600197300000800213490000600221520148900227653001501716653002301731653002401754653003001778653002301808653002101831653002401852653001301876653002001889653002801909100002501937700002101962700002001983700002402003700001902027700002002046700002002066856009002086 2007 eng d a1471-210500aFrom genes to functional classes in the study of biological systems.0 aFrom genes to functional classes in the study of biological syst c2007 Apr 03 a1140 v83 a

BACKGROUND: With the popularization of high-throughput techniques, the need for procedures that help in the biological interpretation of results has increased enormously. Recently, new procedures inspired in systems biology criteria have started to be developed.

RESULTS: Here we present FatiScan, a web-based program which implements a threshold-independent test for the functional interpretation of large-scale experiments that does not depend on the pre-selection of genes based on the multiple application of independent tests to each gene. The test implemented aims to directly test the behaviour of blocks of functionally related genes, instead of focusing on single genes. In addition, the test does not depend on the type of the data used for obtaining significance values, and consequently different types of biologically informative terms (gene ontology, pathways, functional motifs, transcription factor binding sites or regulatory sites from CisRed) can be applied to different classes of genome-scale studies. We exemplify its application in microarray gene expression, evolution and interactomics.

CONCLUSION: Methods for gene set enrichment which, in addition, are independent from the original data and experimental design constitute a promising alternative for the functional profiling of genome-scale experiments. A web server that performs the test described and other similar ones can be found at: http://www.babelomics.org.

10aAlgorithms10aChromosome Mapping10aComputer Simulation10aGene Expression Profiling10aModels, Biological10aMultigene Family10aSignal Transduction10aSoftware10aSystems biology10aUser-Computer Interface1 aAl-Shahrour, Fátima1 aArbiza, Leonardo1 aDopazo, Hernán1 aHuerta-Cepas, Jaime1 aMinguez, Pablo1 aMontaner, David1 aDopazo, Joaquin uhttp://clinbioinfosspa.es/content/genes-functional-classes-study-biological-systems-0