<p><i>Trypanosoma cruzi</i>, the causative agent of Chagas´ disease, is an evolutionarily ancient species with distinct biological and immunological characteristics. A fundamental understanding of the basic biology of the parasite is necessary in order to develop reliable therapeutic and prophylactic agents against <i>T. cruzi</i>. We have, as a part of the <i>T. cruzi</i> genome project launched by the WHO, generated ESTs corresponding to about one third of the functional genes in the parasite. Only about 1/3 of the unique ESTs could be assigned a function upon sequence comparison to all publicly available data. Comparative analysis of the ESTs to functional genes in <i>S.</i> <i>cerevisiae</i> and <i>C. elegans</i> as well as to sequence data from all other kinetoplastids provided primary insights into the evolutionary divergence of <i>T. cruzi.</i> </p><p>A novel dispersed gene family (<i>DGC3</i>) was identified and shown to be present specifically on chromosome 3 and its homologue. Sequence analysis of ten isolated <i>DGC3</i> genes revealed a high sequence similarity of almost 98% among copies. The <i>DGC3</i> genes were transcribed, <i>trans</i>-spliced with the spliced leader and polyadenylated, but did not seem to have any protein-coding property. These data preliminary suggest that it encodes a novel family of functional RNA. </p><p>In the <i>T. cruzi</i> CL Brener strain, the two alleles of a single copy gene encoding the trypanothione synthetase (TcTRS) enzyme appeared to be highly polymorphic. The divergence of the deduced protein sequence was 4%, almost ten-fold higher than another protein, trypanothione reductase, involved in the same pathway. The observed allelic divergence might influence the TcTRS activity thereby having implications for drug design. Moreover, the <i>TcTRS</i> gene was found to be flanked by a number of genes involved in diverse functions and located to a pair of homologous chromosomes with a size difference of about 2 Mbp. </p><p>A gene potentially encoding the polypyrimidine-binding protein (TcPTB) was identified and characterised regarding its organisation and function. The deduced amino acid sequence was shown to comprise four RRM domains generally present in other PTBs. Interestingly, the <i>TcPTB</i> gene appeared to be expressed in a stage-specific manner implicating different functions during parasite development.</p>

Trypanosoma cruzi, the causative agent of Chagas´ disease, is an evolutionarily ancient species with distinct biological and immunological characteristics. A fundamental understanding of the basic biology of the parasite is necessary in order to develop reliable therapeutic and prophylactic agents against T. cruzi. We have, as a part of the T. cruzi genome project launched by the WHO, generated ESTs corresponding to about one third of the functional genes in the parasite. Only about 1/3 of the unique ESTs could be assigned a function upon sequence comparison to all publicly available data. Comparative analysis of the ESTs to functional genes in S. cerevisiae and C. elegans as well as to sequence data from all other kinetoplastids provided primary insights into the evolutionary divergence of T. cruzi. A novel dispersed gene family (DGC3) was identified and shown to be present specifically on chromosome 3 and its homologue. Sequence analysis of ten isolated DGC3 genes revealed a high sequence similarity of almost 98% among copies. The DGC3 genes were transcribed, trans-spliced with the spliced leader and polyadenylated, but did not seem to have any protein-coding property. These data preliminary suggest that it encodes a novel family of functional RNA. In the T. cruzi CL Brener strain, the two alleles of a single copy gene encoding the trypanothione synthetase (TcTRS) enzyme appeared to be highly polymorphic. The divergence of the deduced protein sequence was 4%, almost ten-fold higher than another protein, trypanothione reductase, involved in the same pathway. The observed allelic divergence might influence the TcTRS activity thereby having implications for drug design. Moreover, the TcTRS gene was found to be flanked by a number of genes involved in diverse functions and located to a pair of homologous chromosomes with a size difference of about 2 Mbp. A gene potentially encoding the polypyrimidine-binding protein (TcPTB) was identified and characterised regarding its organisation and function. The deduced amino acid sequence was shown to comprise four RRM domains generally present in other PTBs. Interestingly, the TcPTB gene appeared to be expressed in a stage-specific manner implicating different functions during parasite development.

The duplication of genes followed by selection is perhaps the most prominent way in which molecular biological systems gain multiplicity, diversity and functional complexity in evolution. Whole genome duplications (WGDs) therefore have the potential of generating an extraordinary amount of evolutionary innovation. It is now accepted that the vertebrate lineage has gone through two rounds of WGD in its early stages, after the divergence of invertebrate chordates and before the emergence of jawed vertebrates. These basal vertebrate WGDs are called 2R for two rounds of whole genome duplication. An additional WGD called 3R occurred early in the evolution of teleost fishes, before the radiation of this species-rich group. This thesis describes the evolution of several endocrine and neuronal gene families in relation to the vertebrate WGDs, through a comparative genomic approach including both phylogenetic analyses and chromosomal location data across a wide range of vertebrate taxa. These results show that numerous endocrine gene families have expanded in 2R and in several cases also in 3R. These include the gene families of oxytocin and vasopressin receptors (OT/VP-R), somatostatin receptors (SSTR) and insulin-like growth factor binding proteins (IGFBP). For the OT/VP-R and SSTR families, previously undescribed subtypes were identified. The protein hormone family that includes growth hormone (GH), prolactin (PRL) and somatolactin (SL) acquired a new PRL gene in 2R, however the origins of GH, PRL and SL likely predate 2R. The corresponding family of receptors diversified during different time periods through a combination of local duplications and 3R. Neuronal gene families of the visual system have also expanded in 2R and 3R. The results presented here demonstrate that the vertebrate repertoire of visual opsin genes arose in 2R as part of chromosomal blocks that also include the OT/VP-R genes. The gene families including the transducin alpha, beta and gamma subunits also arose in 2R, hinting at the importance of these events in the diversification and specialization of phototransduction cascades for rods and cones. Thus, the whole genome duplications have been important contributors to the evolution of both vision and endocrine regulation in the vertebrates.

[Français] Une fraction importante des génomes eucaryotes est constituée de Gènes Répétés en Tandem (GRT). Un mécanisme fondamental dans l’évolution des GRT est la recombinaison inégale durant la méiose, entrainant la duplication locale (en tandem) de segments chromosomiques contenant un ou plusieurs gènes adjacents. Différents algorithmes ont été proposés pour inférer une histoire de duplication en tandem pour un cluster de GRT. Cependant, leur utilisation est limitée dans la pratique, car ils ne tiennent pas compte d’autres événements évolutifs pourtant fréquents, comme les inversions, les duplications inversées et les délétions. Cette thèse propose différentes approches algorithmiques permettant d’intégrer ces événements dans le modèle de duplication en tandem classique. Nos contributions sont les suivantes: • Intégrer les inversions dans un modèle de duplication en tandem simple (duplication d’un gène à la fois) et proposer un algorithme exact permettant de calculer le nombre minimal d’inversions s’étant produites dans l’évolution d’un cluster de GRT. • Généraliser ce modèle pour l’étude d’un ensemble de clusters orthologues dans plusieurs espèces. • Proposer un algorithme permettant d’inférer l’histoire évolutive d’un cluster de GRT en tenant compte des duplications en tandem, duplications inversées, inversions et délétions de segments chromosomiques contenant un ou plusieurs gènes adjacents. / [English] Tandemly arrayed genes (TAGs) represent an important fraction of most genomes. A fundamental mechanism at the origin of TAG clusters is unequal crossing-over during meiosis, leading to the duplication of chromosomal segments containing one or many adjacent genes. Such duplications are called tandem duplications, as the duplicated segment is placed next to the original one on the chromosome. Different algorithms have been proposed to infer the tandem duplication history of a TAG cluster. However, their applicability is limited in practice since they do not take into account other frequent evolutionary events such as inversion, inverted duplication and deletion. In this thesis, we propose different algorithmic approaches allowing to integrate these evolutionary events in the original tandem duplication model of evolution. Our contributions are summarized as follows: • We integrate inversion events in a tandem duplication model restricted to single gene duplications, and we propose an exact algorithm allowing to compute the minimum number of inversions explaining the evolution of a TAG cluster. • We generalize this model to the study of orthologous TAG clusters in different species. • We propose an algorithm allowing to infer the evolutionary history of a TAG cluster through tandem duplication, inverted duplication, inversion and deletion of chromosomal segments containing one or many adjacent genes.

Une réconciliation entre un arbre de gènes et un arbre d’espèces décrit une histoire d’évolution des gènes homologues en termes de duplications et pertes de gènes. Pour inférer une réconciliation pour un arbre de gènes et un arbre d’espèces, la parcimonie est généralement utilisée selon le nombre de duplications et/ou de pertes. Les modèles de réconciliation sont basés sur des critères probabilistes ou combinatoires. Le premier article définit un modèle combinatoire simple et général où les duplications et les pertes sont clairement identifiées et la réconciliation parcimonieuse n’est pas la seule considérée. Une architecture de toutes les réconciliations est définie et des algorithmes efficaces (soit de dénombrement, de génération aléatoire et d’exploration) sont développés pour étudier les propriétés combinatoires de l’espace de toutes les réconciliations ou seulement les plus parcimonieuses. Basée sur le processus classique nommé naissance-et-mort, un algorithme qui calcule la vraisemblance d’une réconciliation a récemment été proposé. Le deuxième article utilise cet algorithme avec les outils combinatoires décrits ci-haut pour calculer efficacement (soit approximativement ou exactement) les probabilités postérieures des réconciliations localisées dans le sous-espace considéré. Basé sur des taux réalistes (selon un modèle probabiliste) de duplication et de perte et sur des données réelles/simulées de familles de champignons, nos résultats suggèrent que la masse probabiliste de toute l’espace des réconciliations est principalement localisée autour des réconciliations parcimonieuses. Dans un contexte d’approximation de la probabilité d’une réconciliation, notre approche est une alternative intéressante face aux méthodes MCMC et peut être meilleure qu’une approche sophistiquée, efficace et exacte pour calculer la probabilité d’une réconciliation donnée. Le problème nommé Gene Tree Parsimony (GTP) est d’inférer un arbre d’espèces qui minimise le nombre de duplications et/ou de pertes pour un ensemble d’arbres de gènes. Basé sur une approche qui explore tout l’espace des arbres d’espèces pour les génomes considérés et un calcul efficace des coûts de réconciliation, le troisième article décrit un algorithme de Branch-and-Bound pour résoudre de façon exacte le problème GTP. Lorsque le nombre de taxa est trop grand, notre algorithme peut facilement considérer des relations prédéfinies entre ensembles de taxa. Nous avons testé notre algorithme sur des familles de gènes de 29 eucaryotes. / A reconciliation between a gene tree and a species tree depicts an evolutionary scenario of the homologous genes in terms of gene duplications and gene losses. To infer such a reconciliation given a gene tree and a species tree, parsimony is generally used according to the number of gene duplications and/or losses. The combinatorial models of reconciliation are based on probabilistic or combinatorial criteria. The first paper defines a simple and more general combinatorial model of reconciliation which clearly identifies duplication and loss events and does not only induce the most parsimonious reconciliation. An architecture of all possible reconciliations is developed together with efficient algorithms (that is counting, randomization, and exploration) to study combinatorial properties of the space of all reconciliations or only the most parsimonious ones. Based on the classical birth-death process, an algorithm that computes the likelihood of a reconciliation has recently been proposed. The second paper uses this algorithm together with the combinatorial tools described above to compute efficiently, either exactly or approximately, the posterior probability of the reconciliations located in the considered subspace. Based on realistic gene duplication and loss rates and on real/simulated datasets of fungal gene families, our results suggest that the probability mass of the whole space of reconciliations is mostly located around the most parsimonious ones. In the context of posterior probability approximation, our approach is a valuable alternative to a MCMC method and can competes against a sophisticated, efficient, and exact computation of the probability of a given reconciliation. The Gene Tree Parsimony (GTP) problem is to infer a species tree that minimizes the number of duplications and/or losses over a set of gene family trees. Based on a new approch that explores the whole species tree space for the considered taxa and an efficient computation of the reconciliation cost, the third paper describes a Branch-and- Bound algorithm that solves exactly the GTP problem. When the considered number of taxa is too large, our algorithm can naturally take into account predefined relationships between sets of taxa. We test our algorithm on a dataset of eukaryotic gene families spanning 29 taxa.

ABSTRACT The xyloglucan endo-transglycosylase/hydrolase (XTH) genes are found in all vascular and some nonvascular plants. The XTH genes encode proteins which comprise a subfamily of glycoside hydrolase (GH) family 16 in the Carbohydrate-Active enZYmes (CAZY) classification. The XTH gene products are believed to play intrinsic role in cell wall modification during growth and development throughout the lifetime of the plant. In the present investigation, biochemical and reverse genetic approaches were used to better understand the functions of individual members of the XTH gene family of two important plants: the model organism Arabidopsis thaliana and the grain crop barley (Hordeum vulgare). A phylogenetic tree of the xyloglucan-active enzymes of GH16 has previously been constructed, where enzymes with similar activities have been shown to cluster together. Several members of phylogenetic Group I/II and III-B, predicted to exhibit xyloglucan endo-transglycosylase activity (EC 2.4.1.207) and members of Group III-A, predicted to exhibit xyloglucan endo-hydrolase activity (EC 3.2.1.151), were included to analyze the functional diversity of XTH gene products. A heterologous expression system using the yeast Pichia pastoris was found to be effective for recombinant protein production with a success rate of ca. 50%. XTH gene products were obtained in soluble and active forms for subsequent biochemical characterization. In order to be able to screen larger numbers of protein producing clones, a fast and easy method is required to identify clones expressing active protein in high enough amounts. Thus, a miniaturized XET/XEH assay for high-throughput analysis was developed, which was able to identify activities with good precision and with a reduced time and materials consumption and a reduced work load. Enzyme kinetic analysis indicated that the XET or XEH activity of all XTH gene products characterized in the present study corresponded to predictions based on the previously revised phylogenetic clustering. To gain insight into the biological function of the predominant XEHs AtXTH31 and AtXTH32, which are highly expressed in rapidly developing tissues, a reverse genetic approach was employed using T-DNA insertion lines of the A. thaliana Columbia ecotype. Genotypic and phenotypic characterization, together with in situ assays of XET and XEH activities, in single- and double-knock-out mutants indicated that these Group III-A enzymes are active in expanding tissues of the A. thaliana roots and hypocotyl.  Although suppression of in muro XEH activity was clearly observed in the double-knock-out, no significant growth phenotype was observed, with the exception that radicle emergence appeared to be faster than in the wild type plants. Keywords: Arabidopis thaliana, Hordeum vulgare, plant cell wall, xyloglucan, glycoside hydrolase family 16, xyloglucan endo-transglycosylase/hydrolase gene family, xyloglucan endo-transglycosylase, xyloglucan endo-hydrolase, heterologous protein expression, Pichia pastoris, T-DNA insertion, in situ XET/XEH assay, high-throughput screening / QC 20110114

Une réconciliation entre un arbre de gènes et un arbre d’espèces décrit une histoire d’évolution des gènes homologues en termes de duplications et pertes de gènes. Pour inférer une réconciliation pour un arbre de gènes et un arbre d’espèces, la parcimonie est généralement utilisée selon le nombre de duplications et/ou de pertes. Les modèles de réconciliation sont basés sur des critères probabilistes ou combinatoires. Le premier article définit un modèle combinatoire simple et général où les duplications et les pertes sont clairement identifiées et la réconciliation parcimonieuse n’est pas la seule considérée. Une architecture de toutes les réconciliations est définie et des algorithmes efficaces (soit de dénombrement, de génération aléatoire et d’exploration) sont développés pour étudier les propriétés combinatoires de l’espace de toutes les réconciliations ou seulement les plus parcimonieuses. Basée sur le processus classique nommé naissance-et-mort, un algorithme qui calcule la vraisemblance d’une réconciliation a récemment été proposé. Le deuxième article utilise cet algorithme avec les outils combinatoires décrits ci-haut pour calculer efficacement (soit approximativement ou exactement) les probabilités postérieures des réconciliations localisées dans le sous-espace considéré. Basé sur des taux réalistes (selon un modèle probabiliste) de duplication et de perte et sur des données réelles/simulées de familles de champignons, nos résultats suggèrent que la masse probabiliste de toute l’espace des réconciliations est principalement localisée autour des réconciliations parcimonieuses. Dans un contexte d’approximation de la probabilité d’une réconciliation, notre approche est une alternative intéressante face aux méthodes MCMC et peut être meilleure qu’une approche sophistiquée, efficace et exacte pour calculer la probabilité d’une réconciliation donnée. Le problème nommé Gene Tree Parsimony (GTP) est d’inférer un arbre d’espèces qui minimise le nombre de duplications et/ou de pertes pour un ensemble d’arbres de gènes. Basé sur une approche qui explore tout l’espace des arbres d’espèces pour les génomes considérés et un calcul efficace des coûts de réconciliation, le troisième article décrit un algorithme de Branch-and-Bound pour résoudre de façon exacte le problème GTP. Lorsque le nombre de taxa est trop grand, notre algorithme peut facilement considérer des relations prédéfinies entre ensembles de taxa. Nous avons testé notre algorithme sur des familles de gènes de 29 eucaryotes. / A reconciliation between a gene tree and a species tree depicts an evolutionary scenario of the homologous genes in terms of gene duplications and gene losses. To infer such a reconciliation given a gene tree and a species tree, parsimony is generally used according to the number of gene duplications and/or losses. The combinatorial models of reconciliation are based on probabilistic or combinatorial criteria. The first paper defines a simple and more general combinatorial model of reconciliation which clearly identifies duplication and loss events and does not only induce the most parsimonious reconciliation. An architecture of all possible reconciliations is developed together with efficient algorithms (that is counting, randomization, and exploration) to study combinatorial properties of the space of all reconciliations or only the most parsimonious ones. Based on the classical birth-death process, an algorithm that computes the likelihood of a reconciliation has recently been proposed. The second paper uses this algorithm together with the combinatorial tools described above to compute efficiently, either exactly or approximately, the posterior probability of the reconciliations located in the considered subspace. Based on realistic gene duplication and loss rates and on real/simulated datasets of fungal gene families, our results suggest that the probability mass of the whole space of reconciliations is mostly located around the most parsimonious ones. In the context of posterior probability approximation, our approach is a valuable alternative to a MCMC method and can competes against a sophisticated, efficient, and exact computation of the probability of a given reconciliation. The Gene Tree Parsimony (GTP) problem is to infer a species tree that minimizes the number of duplications and/or losses over a set of gene family trees. Based on a new approch that explores the whole species tree space for the considered taxa and an efficient computation of the reconciliation cost, the third paper describes a Branch-and- Bound algorithm that solves exactly the GTP problem. When the considered number of taxa is too large, our algorithm can naturally take into account predefined relationships between sets of taxa. We test our algorithm on a dataset of eukaryotic gene families spanning 29 taxa.

Les duplicacions segmentàries (DSs), o també anomenades duplicons o Low copy Repeats (LCRs), són regions de coma mínim 1 kb amb un alt nivell d'identitat (>90%), que estan presents almenys dues vegades en el genoma. La regió 8p23.1 consta de 6.5 Mb a la part distal del braç curt del cromosoma 8 i està flanquejada per duplicacions segmentàries. Degut a la seva arquitectura genòmica aquesta regió és susceptible a patir reordenaments mediats per recombinació homòloga no al·lèlica entre les DSs, com per exemple la inversió polimòrfica de 8p23.1 [inv(8)(p23)], present en un de cada quatre individus de la població general europea i japonesa, així com d'altres reorganitzacions menys corrents.El treball realitzat en aquesta tesi doctoral pretén aprofundir en la caracterització de la complexa arquitectura genòmica d'aquesta regió. En la nostra primera aproximació a l'estudi de les DSs que flanquegen la regió cromosòmica 8p23.1, es va identificar una nova família gènica específica de primats, la família gènica FAM90A.Així, bona part d'aquesta tesi doctoral està centrada en l'anàlisi de l'origen, formació, evolució i expansió de FAM90A en els homínids. Per altra banda també s'ha analitzant en detall la variabilitat de FAM90A com a variant en número de còpia (CNV) en diferents poblacions humanes.Finalment, s'ha establert la freqüència de la inversió que afecta a 8p23.1 en població espanyola. També s'ha procedit a genotipar diversos individus homozigots per la inversió i s'ha predit l' estatus de la inversió en 150 individus del projecte HapMap i s'ha analitzat l'efecte que té aquesta reorganització sobre els nivells d'expressió dels gens de la regió. / Segmental duplications (SDs), also known as duplicons or Low Copy Repeats (LCRs), are regions of a minimum of 1 kb with a high sequence identity level (>90%), which are present at least two times in the genome. The 8p23.1 region extends 6.5 Mb at the distal part of the short arm of chromosome 8 and it is flanked by segmental duplications. Due to its genomic architecture the region is prone to suffer rearrangements mediated by non-allelic homologous recombination between these SDs, such as the polymorphic inversion of 8p23.1 [inv(8)(p23)], which is present in one out of every four of European and Japanese general population individuals, as well as other less frequent rearrangements.The aim of the work presented in this doctoral thesis is to get insights in the characterization of the genomic architecture of this complex region. Our first approach to study the SDs flanking 8p23.1 region resulted in the identification of a novel gene family which is primate specific, the FAM90A gene family. Thus, this doctoral thesis is mainly focused on the analysis of the origins, formation, evolution and expansion of FAM90A in hominoids. It has also been analyzed in detail the variability of FAM90A as a copy number variant (CNV) in different human populations.Finally, it has been established the frequency of the inversion affecting 8p23.1 region in the Spanish population. Several homozygous inverted individuals have been genotyped and the status for the inversion has been predicted for 150 HapMap individuals, as well as the effect of this rearrangement on the gene expression levels of the genes contained in the region.