Bioinformàtica : comparació de genomes d'eucariota

La cerca de similituds als codis genètics de dos espècies, ens permet obtenir molta informació de la evolució dels seus genomes. Aquesta informació afavoreix el descobriment de gens que es conserven amb la mateixa funcionalitat a diferents espècies. També té importants aplicacions mèdiques i ens permet entendre els processos evolutius que han portat a la diversitat d'espècies de l'actualitat. El present treball té l'objectiu d'automatitzar una sèrie de processos d'un servidor d'aplicacions web: http://platypus.uab.cat, que realitzin de forma òptima i eficient, la comparació dels genomes eucariotes, tots amb tots, conforme aquests genomes siguin seqüenciats. Així aquestes comparacions entre genomes de organismes superiors podran ser consultades via web.

Relaxation of selective constraints causes independent selenoprotein extinction in insect genomes

BACKGROUND: Selenoproteins are a diverse family of proteins notable for the presence of the 21st amino acid, selenocysteine. Until very recently, all metazoan genomes investigated encoded selenoproteins, and these proteins had therefore been believed to be essential for animal life. Challenging this assumption, recent comparative analyses of insect genomes have revealed that some insect genomes appear to have lost selenoprotein genes. METHODOLOGY/PRINCIPAL FINDINGS: In this paper we investigate in detail the fate of selenoproteins, and that of selenoprotein factors, in all available arthropod genomes. We use a variety of in silico comparative genomics approaches to look for known selenoprotein genes and factors involved in selenoprotein biosynthesis. We have found that five insect species have completely lost the ability to encode selenoproteins and that selenoprotein loss in these species, although so far confined to the Endopterygota infraclass, cannot be attributed to a single evolutionary event, but rather to multiple, independent events. Loss of selenoproteins and selenoprotein factors is usually coupled to the deletion of the entire no-longer functional genomic region, rather than to sequence degradation and consequent pseudogenisation. Such dynamics of gene extinction are consistent with the high rate of genome rearrangements observed in Drosophila. We have also found that, while many selenoprotein factors are concomitantly lost with the selenoproteins, others are present and conserved in all investigated genomes, irrespective of whether they code for selenoproteins or not, suggesting that they are involved in additional, non-selenoprotein related functions. CONCLUSIONS/SIGNIFICANCE: Selenoproteins have been independently lost in several insect species, possibly as a consequence of the relaxation in insects of the selective constraints acting across metazoans to maintain selenoproteins. The dispensability of selenoproteins in insects may be related to the fundamental differences in antioxidant defense between these animals and other metazoans.

Comparative analysis of cancer genes in the human and chimpanzee genomes

Background: Cancer is a major medical problem in modern societies. However, the incidence of this disease in non-human primates is very low. To study whether genetic differences between human and chimpanzee could contribute to their distinct cancer susceptibility, we have examined in the chimpanzee genome the orthologous genes of a set of 333 human cancer genes. Results: This analysis has revealed that all examined human cancer genes are present in chimpanzee, contain intact open reading frames and show a high degree of conservation between both species. However, detailed analysis of this set of genes has shown some differences in genes of special relevance for human cancer. Thus, the chimpanzee gene encoding p53 contains a Pro residue at codon 72, while this codon is polymorphic in humans and can code for Arg or Pro, generating isoforms with different ability to induce apoptosis or interact with p73. Moreover, sequencing of the BRCA1 gene has shown an 8 Kb deletion in the chimpanzee sequence that prematurely truncates the co-regulated NBR2 gene. Conclusion: These data suggest that small differences in cancer genes, as those found in tumor suppressor genes, might influence the differences in cancer susceptibility between human and chimpanzee. Nevertheless, further analysis will be required to determine the exact contribution of the genetic changes identified in this study to the different cancer incidence in non-human primates.

Sequencing primate genomes: what have we learned?

We summarize the progress in whole-genome sequencing and analyses of primate genomes. These emerging genome datasets have broadened our understanding of primate genome evolution revealing unexpected and complex patterns of evolutionary change. This includes the characterization of genome structural variation, episodic changes in the repeat landscape, differences in gene expression, new models regarding speciation, and the ephemeral nature of the recombination landscape. The functional characterization of genomic differences important in primate speciation and adaptation remains a significant challenge. Limited access to biological materials, the lack of detailed phenotypic data and the endangered status of many critical primate species have significantly attenuated research into the genetic basis of primate evolution. Next-generation sequencing technologies promise to greatly expand the number of available primate genome sequences; however, such draft genome sequences will likely miss critical genetic differences within complex genomic regions unless dedicated efforts are put forward to understand the full spectrum of genetic variation.

The non-LTR retrotransposons in Ciona intestinalis: new insights into the evolution of chordate genomes

Background: Non-long terminal repeat (non-LTR) retrotransposons have contributed to shaping the structure and function of genomes. In silico and experimental approaches have been used to identify the non-LTR elements of the urochordate Ciona intestinalis. Knowledge of the types and abundance of non-LTR elements in urochordates is a key step in understanding their contribution to the structure and function of vertebrate genomes. Results: Consensus elements phylogenetically related to the I, LINE1, LINE2, LOA and R2 elements of the 14 eukaryotic non-LTR clades are described from C. intestinalis. The ascidian elements showed conservation of both the reverse transcriptase coding sequence and the overall structural organization seen in each clade. The apurinic/apyrimidinic endonuclease and nucleic-acid-binding domains encoded upstream of the reverse transcriptase, and the RNase H and the restriction enzyme-like endonuclease motifs encoded downstream of the reverse transcriptase were identified in the corresponding Ciona families. Conclusions: The genome of C. intestinalis harbors representatives of at least five clades of non-LTR retrotransposons. The copy number per haploid genome of each element is low, less than 100, far below the values reported for vertebrate counterparts but within the range for protostomes. Genomic and sequence analysis shows that the ascidian non-LTR elements are unmethylated and flanked by genomic segments with a gene density lower than average for the genome. The analysis provides valuable data for understanding the evolution of early chordate genomes and enlarges the view on the distribution of the non-LTR retrotransposons in eukaryotes.

Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution

Background: Chemoreception is a widespread mechanism that is involved in critical biologic processes, including individual and social behavior. The insect peripheral olfactory system comprises three major multigene families: the olfactory receptor (Or), the gustatory receptor (Gr), and the odorant-binding protein (OBP) families. Members of the latter family establish the first contact with the odorants, and thus constitute the first step in the chemosensory transduction pathway.Results: Comparative analysis of the OBP family in 12 Drosophila genomes allowed the identification of 595 genes that encode putative functional and nonfunctional members in extant species, with 43 gene gains and 28 gene losses (15 deletions and 13 pseudogenization events). The evolution of this family shows tandem gene duplication events, progressive divergence in DNA and amino acid sequence, and prevalence of pseudogenization events in external branches of the phylogenetic tree. We observed that the OBP arrangement in clusters is maintained across the Drosophila species and that purifying selection governs the evolution of the family; nevertheless, OBP genes differ in their functional constraints levels. Finally, we detect that the OBP repertoire evolves more rapidly in the specialist lineages of the Drosophila melanogaster group (D. sechellia and D. erecta) than in their closest generalists.Conclusion: Overall, the evolution of the OBP multigene family is consistent with the birth-and-death model. We also found that members of this family exhibit different functional constraints, which is indicative of some functional divergence, and that they might be involved in some of the specialization processes that occurred through the diversification of the Drosophila genus.

CBS: an open platform that integrates predictive methods and epigenetics information to characterize conserved regulatory features in multiple Drosophila genomes.

Background: Information about the composition of regulatory regions is of great value for designing experiments to functionally characterize gene expression. The multiplicity of available applications to predict transcription factor binding sites in a particular locus contrasts with the substantial computational expertise that is demanded to manipulate them, which may constitute a potential barrier for the experimental community. Results: CBS (Conserved regulatory Binding Sites, http://compfly.bio.ub.es/CBS) is a public platform of evolutionarily conserved binding sites and enhancers predicted in multiple Drosophila genomes that is furnished with published chromatin signatures associated to transcriptionally active regions and other experimental sources of information. The rapid access to this novel body of knowledge through a user-friendly web interface enables non-expert users to identify the binding sequences available for any particular gene, transcription factor, or genome region. Conclusions: The CBS platform is a powerful resource that provides tools for data mining individual sequences and groups of co-expressed genes with epigenomics information to conduct regulatory screenings in Drosophila.

Consequences of lineage-specific gene loss on functional evolution of surviving paralogs: ALDH1A and retinoic acid signaling in vertebrate genomes

Genome duplications increase genetic diversity and may facilitate the evolution of gene subfunctions. Little attention, however, has focused on the evolutionary impact of lineage-specific gene loss. Here, we show that identifying lineage-specific gene loss after genome duplication is important for understanding the evolution of gene subfunctions in surviving paralogs and for improving functional connectivity among human and model organism genomes. We examine the general principles of gene loss following duplication, coupled with expression analysis of the retinaldehyde dehydrogenase Aldh1a gene family during retinoic acid signaling in eye development as a case study. Humans have three ALDH1A genes, but teleosts have just one or two. We used comparative genomics and conserved syntenies to identify loss of ohnologs (paralogs derived from genome duplication) and to clarify uncertain phylogenies. Analysis showed that Aldh1a1 and Aldh1a2 form a clade that is sister to Aldh1a3-related genes. Genome comparisons showed secondarily loss of aldh1a1 in teleosts, revealing that Aldh1a1 is not a tetrapod innovation and that aldh1a3 was recently lost in medaka, making it the first known vertebrate with a single aldh1a gene. Interestingly, results revealed asymmetric distribution of surviving ohnologs between co-orthologous teleost chromosome segments, suggesting that local genome architecture can influence ohnolog survival. We propose a model that reconstructs the chromosomal history of the Aldh1a family in the ancestral vertebrate genome, coupled with the evolution of gene functions in surviving Aldh1a ohnologs after R1, R2, and R3 genome duplications. Results provide evidence for early subfunctionalization and late subfunction-partitioning and suggest a mechanistic model based on altered regulation leading to heterochronic gene expression to explain the acquisition or modification of subfunctions by surviving ohnologs that preserve unaltered ancestral developmental programs in the face of gene loss.

Manipulació genètica de vectors virals per a estratègies de teràpia gènica en malalties autoimmunes

Els tractaments actuals per a la Malaltia de Crohn es basen en l’administració de proteïnes recombinants amb capacitat immunosupressora. Aquests tractaments són sistèmics, crònics i causen el desenvolupament d’efectes secundaris severs com infeccions amb patògens oportunistes i una major incidència de tumors. Passar d’una administració sistèmica dels immunosupressors a una administració local sembla clau per evitar la severitat i el nombre de complicacions secundàries. Ha sigut descrit recentment que el subconjunt cel•lular Th-17 juga un paper central en moltes malalties autoimmunes en humans, com ara la Malaltia de Crohn, Escleroderma, Artritis reumatoide o Psoriasis. Per determinar com la via del Th-17 es veu afectada durant el desenvolupament i els estadis aguts de la malaltia, analitzarem una gran bateria de interleuquines/citoquines i els mecanismes intracel•lulars d’activació/inhibició de l’activitat de les cèl•lules del Sistema Immune. En un segon pas, clonarem els gens seleccionats relacionats amb la immunomodulació de la via del Th-17 en vectors recentment desenvolupats (adenovirus quimera 5/40 o diferents pseudotips d’AAVs ), ja que hipotitzem que una producció local d’immunomoduladors s’assolirà amb una eficiència major fent servir aquests vectors, evitant així els efectes secundaris. A aquests efectes, hem començat el projecte clonant els gens de la IL10 i del receptor transmembrenal de la IL23 en adenovirus i en genomes d’AAV.

Bioinformática: interfaz gráfica para comparación de genomas vía web

El proyecto consiste en un entorno gráfico cuyo fin es el de visualizar, estudiar e interpretar la conservación de código genético existente entre los diferentes genomas. Una interface que permite cargar hasta ocho genomas para ser comparados en detalle, por pares o entre todos ellos a la vez. El gráfico que se muestra en la interfaz, representa los Maximal Unique Matchings entre cada par de genomas, lo que significa coincidencias de la mayor longitud posible no repetidas, en las secuencias de ADN de las especies comparadas. La finalidad es el estudio de las evoluciones que han ido apareciendo entre diferentes organismos o los genes que comparten unas especies con otras.

Desarrollo de una colección de Líneas Casi Isogénicas (NIL) para el estudio de caracteres de interés agronómico en especies del género Fragaria

Fragaria vesca posee varias características que la hacen interesante como especie modelo en la familia Rosaceae. Además, su naturaleza diploide permite sortear la complejidad genética de la fresa cultivada. Considerando que los genomas de las especies diploides y octoploides de Fragaria mantienen una relación de sintenia y colinearidad, estamos desarrollando una colección de Líneas Casi Isogénicas (NIL) en Fragaria diploide, usando Fragaria bucharica, una variedad asiática, como donante de introgresiones y Fragaria vesca var. Reine des Vallées, una variedad francesa comúnmente cultivada en España para usos industriales, como parental recurrente. Para obtener introgresiones de F. bucharica en un fondo genético homogéneo de F. vesca, se desarrolló un retrocruzamiento y se obtuvo una población BC1 que fue analizada para la presencia de alelos de F. bucharica a lo largo de los 7 grupos de ligamiento (LG) del mapa de referencia de Fragaria. Los individuos con bajo número de introgresiones de gran tamaño, que en conjunto abarcaron todo el genoma de F. bucharica, fueron seleccionados y retrocruzados con el parental recurrente. La descendencia fue analizada para los loci segregantes y los individuos BC2 con 1 ó 2 introgresiones fueron seleccionados como líneas donadoras de introgresiones de las NIL. Hasta el momento se han descrito tres fenotipos dominantes diferentes bajo condiciones de invernadero en varias NIL heterozigóticas: Las plantas con introgresiones en LG2 presentan estolones. Las plantas con introgresiones en LG1 y LG3 presentan floración temprana. Las plantas con introgresiones en LG6 presentan floración estacional. Actualmente se está trabajando en la selección y caracterización de introgresiones de pequeño tamaño mediante la autofecundación de las líneas seleccionadas. La caracterización fenotípica de los recombinantes seleccionados permitirá localizar y estimar los patrones de herencia de los genes implicados en los caracteres descritos, además de la identificación de nuevos caracteres recesivos que aparecerán en homozigosis.

Origen de la multicel·lularitat : una aproximació genòmica i funcional

En el present estudi s'analitza l'origen i evolució de 2 molècules claus pera entendre la multicel·lularitat dels animals: les molècules d'adhesió integrines i els factors de transcripció T-box. S’utilitzen els genomes recentment publicats de protists unicel•lulars parents propers dels animals. S’analitza l’origen i evolució d’aquests gens mitjançant anàlisi filogènic, determinació de motius funcionals i també tècniques de biologia molecular. A més, es documenta un cas de transferència gènica horitzontal des d'un eucariota cap a un procariota, fenomen poc habitual. Les principals conclusions són que tant l’adhesoma d'integrina com els gens T-box tenen un origen molt anterior als animals, en un context unicel•lular, i que després foren cooptats pel llinatge multicel•lular dels animals.

Bioinformática: identificar genes en una interfaz gráfica vía web para la comparación de genomas

Este trabajo desarrolla el proceso de diseño e implementación de una interfaz web que permite la exploración en detalle de las relaciones entre genomas completos. La interfaz permite la comparación simultánea de nueve genomas, representando en cada gráfica las relaciones entre cada par de genomas junto los genes identificados de cada uno de ellos. Es capaz de trabajar con genomas del dominio Eukaryota y se adapta a la capacidad de cómputo de la máquina cliente. La información representada son MUMs (Maximal Unique Matching, secuencia máxima y única encontrada en ambos genomas) y SuperMUMs (agrupación de MUMs mediante Approximate String Matching). Los datos son previamente calculados y accesibles desde un servidor web.

Parallelization of whole genome alignment

With the advent of High performance computing, it is now possible to achieve orders of magnitude performance and computation e ciency gains over conventional computer architectures. This thesis explores the potential of using high performance computing to accelerate whole genome alignment. A parallel technique is applied to an algorithm for whole genome alignment, this technique is explained and some experiments were carried out to test it. This technique is based in a fair usage of the available resource to execute genome alignment and how this can be used in HPC clusters. This work is a rst approximation to whole genome alignment and it shows the advantages of parallelism and some of the drawbacks that our technique has. This work describes the resource limitations of current WGA applications when dealing with large quantities of sequences. It proposes a parallel heuristic to distribute the load and to assure that alignment quality is mantained.

Análisis bioinformáticos sobre la tecnología Hadoop

Desde el inicio del proyecto del genoma humano y su éxito en el año 2001 se han secuenciado genomas de multitud de especies. La mejora en las tecnologías de secuenciación ha generado volúmenes de datos con un crecimiento exponencial. El proyecto Análisis bioinformáticos sobre la tecnología Hadoop abarca la computación paralela de datos biológicos como son las secuencias de ADN. El estudio ha sido encauzado por la naturaleza del problema a resolver. El alineamiento de secuencias genéticas con el paradigma MapReduce.