Identification of signatures of positive selection has long been a major issue for understanding the unique features of any given species. However, only a fraction of human genes have been interrogated. Genes within segmental duplications are usually omitted due to the limitations of draft genome assemblies and the methodological reliance on accurate gene trees. In this work, we show the feasibility of a new method that does not need accurate gene trees or individual high-quality assemblies. We applied the concept to study exon evolution in the human genome, identifying 74 exons with evidence for rapid coding sequence evolution during human and Old World monkey evolution. Our results suggest abundant accelerated coding sequence evolution within duplicated regions of the genome and provide a more comprehensive view of the role of selection on the human genome. / La identificación de señales debidas a la acción de la selección positiva es de gran relevancia para desvelar características únicas de las especies. A pesar de ello, solo una fracción de genes humanos han sido analizados. Los genes incluidos en duplicaciones segmentarias son normalmente ignorados debido a limitaciones impuestas por la naturaleza preliminar de los genomas distintos al humano, así como por la dependencia en adecuados árboles filogenéticos. En este proyecto, demostramos la viabilidad de un nuevo método que no necesita árboles filogenéticos correctos ni ensamblajes de genomas de alta calidad. Hemos aplicado el concepto al genoma humano y hemos identificado 74 exones que muestran evidencia de haber evolucionado más rápidamente desde la separación de los humanos y los monos del viejo mundo. Nuestros resultados sugieren que ha habido abundante evolución acelerada dentro de las regiones duplicadas y ofrece una visión más esclarecedora del rol de la selección en la evolución del genoma humano.

Fibroblast growth factor (Fgf) signaling plays multiple inductive roles during development of vertebrates (Itoh 2007). Some Fgfs, such as Fgf8, are locally secreted and signal over a long range to provide positional information in the target tissue (Scholpp and Brand 2004). Fgf ligands signal in a receptor-dependent manner via tyrosine kinase receptors, four of which have been so far identified. Fgf8 signaling was shown to depend both on receptor activation as well as endocytosis. The specificity of Fgf ligands and receptors as well as the function of receptors in the control of the Fgf signaling range have been, however, largely unclear. In this study, we show that the putative Fgf8 receptor Fgfr1 is duplicated in zebrafish and that it acts redundantly in the formation of the posterior mesoderm. Also, in overexpression studies we confirm the notion that receptor endocytosis influences Fgf8 signaling range. Through TILLING mutant recovery and morpholino knockdown studies we also show that Fgfr2 is required for growth and skeletal development in zebrafish, whereas Fgfr4 is required for pectoral fin specification and growth.

Odorant-binding proteins (OBPs) are small soluble proteins expressed at high levels in the proximity of olfactory receptors. OBPs act as solubilizers and carriers of the lipophilic odorants in the aqueous mucus of mammals and other vertebrates. OBPs have now been studied nearly thirty years, but in comparison to the wealth of data available on their structural chemistry and molecular biology, our knowledge about gene expression and function of these proteins is still insufficient. This work provides new insights into the tissue specificity of OBP and presents several new sequences of genes governing these proteins in selected species of mice.

Duplications and deletions are a major source of genomic variation. Duplications, specifically, have a significant impact on gene genesis and dosage, and the malaria parasite P. falciparum has developed resistance to a growing number of anti-malarial drugs via gene duplication. It also contains highly duplicated families of antigenically variable allelic genes. While specific genes and families have been studied, a comprehensive analysis of duplications and deletions within the reference genome and population has not been performed. We analyzed the extent of segmental duplications (SD) in the reference genome for P. falciparum, primarily by a whole genome self alignment. We discovered that while 5% of the genome identified as SD, the distribution within the genome was partition clustered, with the vast majority localized to the subtelomeres. Within the SDs, we found an overrepresentation of genes encoding antigenically diverse proteins exposed to the extracellular membrane, specifically the var, rifin, and stevor gene families. To examine variation of duplications and deletions within the parasite populations, we designed a novel computational methodology to identify copy number variants (CNVs) from high throughput sequencing, using a read depth based approach refined with discordant read pairs. After validating the program against in vitro lab cultures, we analyzed isolates from Senegal for initial tests into clinical isolates. We then expanded our search to a global sample of 610 strains from Africa and South East Asia, identifying 68 CNV regions. Geographically, genic CNV were found on average in less than 10% of the population, indicating that CNV are rare. However, CNVs at high frequency were almost exclusively duplications associated with known drug resistant CNVs. We also identified the novel biallelic duplication of the crt gene – containing both the chloroquine resistant and sensitive allele. The synthesis of our SD and CNV analysis indicates a CNV conservative P. falciparum genome except where drug and human immune pressure select for gene duplication.

Nutrient sensing kinases lie at the heart of cellular health and homeostasis, allowing cells to quickly adapt to changing environments. Target of Rapamycin (TOR) and PAS kinase (PASK, or PASKIN) are two such nutrient kinases, conserved from yeast to man. In yeast, these kinases each have paralogs. The two TOR paralogs in yeast mimic the mammalian TORC1 and TORC2 complexes, except both Tor1 and Tor2 may contribute to TORC1 or TORC2 function. The two PAS kinase paralogs are paired with the TOR paralogs, meaning that both Psk1 and Psk2 regulate TORC1, while Psk2 suppresses a temperature-sensitive allele of Tor2. Herein we review the evolutionary models for these paralogs, their function in yeast and mammalian cells, as well as the overlapping function of PAS kinase and TOR. We also use Rice University’s Direct Coupling Analysis algorithms to analyze co-evolutionary relationships and identify potential interaction sites between PAS kinase and several of its substrates.

In mono- und eudikotylen Pflanzen kodiert eine Genfamilie (RpoT, RNA-Polymerase des T3/T7-Typs) mitochondriale und plastidäre RNA-Polymerasen (RNAP), die den ungeraden T-Phagen-Polymerasen ähneln. RpoT-Gene von Angiospermen sind gut charakterisiert, während aus tiefer abzweigenden Pflanzenspecies bisher lediglich die Gene aus dem Moos Physcomitrella beschrieben wurden. Um einen Beitrag zur Aufklärung der molekularen Evolution der RpoT-Polymerasen im Pflanzenreich zu liefern und um Erkenntnisse über die potentielle Bedeutung von multiplen Phagen-Typ (RNAP) in Pflanzen zu gewinnen, wurden die RpoT-Gene aus dem Lycophyten Selaginella moellendorffii und aus dem basalen Angiosperm Nuphar advena identifiziert und charakterisiert. Selaginella moellendorffii (Moosfarn)-Trace-Sequenzdaten mit hoher Ähnlichkeit zu RpoT-Sequenzen von Angiospermen wurden benutzt, um das full-length SmRpoT-Gen und die entsprechende cDNA zu isolieren. Die SmRpoT-mRNA ist 3542 nt lang und weist einen offenen Leserahmen von 3006 nt auf, der für ein putatives Protein aus 1002 Aminosäuren mit einer molekularen Masse von 113 kDa kodiert. Das SmRpoT-Gen besteht aus 19 Exons und 18 Introns, die in ihren Positionen mit denen aus den Angiosperm- und Physcomitrella-Genen konserviert sind. Mittels Southernblot-Analyse wurde nachgewiesen, dass S. moellendorffii ein single-copy RpoT-Gen kodiert. Für das N-terminale Transitpeptid von SmRpoT konnte gezeigt werden, dass es bei transienter Expression in Arabidopsis- und Selaginella-Protoplasten den Transport von GFP (green fluorescent protein) exclusiv in Mitochondrien vermittelt. In N. advena wurden mittels Screening einer BAC-Bibliothek drei RpoT-Gene identifiziert. Sowohl die genomischen als auch die cDNA-Sequenzen wurden aufgeklärt. Die NaRpoT-mRNAs kodieren putative Polypeptide von 996, 990 und 985 Aminosären. Alle drei Gene besitzen 19 Exons und 18 Introns, die in ihren Positionen mit denen der RpoT-Gene aus Selaginella und allen anderen Landpflanzen konserviert sind. Die kodierten Proteine weisen auf Aminosäureebene einen hohen Konservierungsgrad auf, einschließlich aller essentiellen Regionen und Aminosäurereste, die für die T7-RNAP bekannt sind. Die N-terminalen Transitpeptide zweier der kodierten RNAP, NaRpoTm1 und NaRpoTm2, vermittelten den Import von GFP exclusiv in Mitochondrien, während die dritte Polymerase, NaRpoTp, in Chloroplasten importiert wurde. Interessanterweise muß die Translation der NaRpoTp-mRNA an einem CUG-Codon initiiert werden, um ein funktionelles Protein mit plastidärem Transitpeptid zu erhalten. Die N. advena RpoTp-RNAP ist somit neben AGAMOUS aus Arabidopsis und der RpoTp-RNAP aus Nicotiana, ein weiteres Beispiel für jene selten vorkommenden pflanzlichen mRNAs, deren Translation exclusiv an nicht-AUG-Codons initiiert wird. Die Rekonstruktion von phylogenetischen Bäumen resultierte in unterschiedlichen Positionen für die Selaginella- und Nuphar-Polymerasen: Im Gegensatz zu der RpoT-Polymerase aus S. moellendorffii und denen aus Physcomitrella, die in den phylogenetischen Analysen Schwesterpositionen zu allen anderen Phagentyp-RNAP der Angiospermen einnehmen, clusterten die Nuphar-RpoTs zusammen mit den deutlich separierten mitochondrialen (NaRpoTm1 und NaRpoTm2) und plastidären (NaRpoTp) Polymerasen. Selaginella kodiert eine einzige mitochondriale RNAP, während Nuphar zwei mitochondriale und eine plastidäre RNAP besitzt. Die Identifizierung einer Plastiden-lokalisierten Phagentyp-RNAP in diesem basalen Eudikotylen, die ortholog zu allen anderen RpoT-Enzymen der Blütenpflanzen ist, läßt darauf schließen, daß die Acquisition einer nukleär kodierten plastidären RNAP, die noch in den Lycopoden fehlt, nach der Trennung der Leucopoden von allen anderen Tracheophyten erfolgte. Eine “dual-targeting” RNAP (mitochondrial und plastidär lokalisiert), wie sie in Eudikotylen, nicht jedoch in Monokotylen vorkommt, wurde weder in Selaginella noch in Nuphar nachgewiesen, vermutlich ist sie ein evolutionäres Novum von eudikotylen Pflanzen wie Arabidopsis. / In mono- and eudicot plants, a small nuclear gene family (RpoT, RNA polymerase of the T3/T7 type) encodes mitochondrial as well as chloroplast RNA polymerases homologous to the T-odd bacteriophage enzymes. RpoT genes from angiosperms are well characterized, whereas data from deeper branching plant species until recently were limited to the moss Physcomitrella. To elucidate the molecular evolution of the RpoT polymerases in the plant kingdom and to get more insight into the potential importance of having more than one phage-type RNA polymerase (RNAP) available, we identified and characterized RpoT genes in the lycophyte Selaginella moellendorffii and the basal eudicot Nuphar advena. Selaginella moellendorffii (spikemoss) sequence trace data encoding a polypeptide highly similar to angiosperm and moss phage-type organelle RNA polymerases were used to isolate a BAC clone containing the full-length gene SmRpoT as well as the corresponding cDNA. The SmRpoT mRNA comprises 3452 nt with an open reading frame of 3,006 nt, encoding a putative protein of 1,002 amino acids with a molecular mass of 113 kDa. The SmRpoT gene comprises 19 exons and 18 introns, conserved in their position with those of the angiosperm and Physcomitrella RpoT genes. Using Southern blot analysis, it was shown that S. moellendorffii encodes a single RpoT gene. The N-terminal transit peptide of SmRpoT was shown to confer targeting of green fluorescent protein (GFP) exclusively to mitochondria after transient expression in Arabidopsis and Selaginella protoplasts. In Nuphar advena three RpoT genes were identified by BAC library screening. Both genomic gene sequences and full-length cDNAs were determined. The NaRpoT mRNAs specify putative polypeptides of 996, 990 and 985 amino acids, respectively. All three genes comprise 19 exons and 18 introns, conserved in their positions with those from S. moellendorffii and the RpoT genes of other land plants. The encoded proteins show a high degree of conservation at the amino acid sequence level, including all functional crucial regions and residues known from the phage T7 RNAP. The N-terminal transit peptides of two of the encoded polymerases, NaRpoTm1 and NaRpoTm2, conferred targeting of GFP exclusively to mitochondria, whereas the third polymerase, NaRpoTp, was targeted to chloroplasts. Remarkably, translation of NaRpoTp mRNA has to be initiated at a CUG codon to generate a functional plastid transit peptide. Thus, besides AGAMOUS in Arabidopsis and the Nicotiana RpoTp polymerase, N. advena RpoTp provides another example for a plant mRNA that is exclusively translated from a non-AUG codon. Reconstruction of phylogenetic trees revealed different positions of the RpoTs from the lycophyte Selaginella and the basal eudicot Nuphar. In contrast to the RpoTs of S. moellendorffii and those of the moss Physcomitrella, which are according to the phylogenetic analyses in sister positions to all other phage-type polymerases of angiosperms, the Nuphar RpoTs clustered with the well separated clades of mitochondrial (NaRpoTm1 and NaRpoTm2) and plastid (NaRpoTp) polymerases. Selaginella encodes a single mitochondrial RNAP, whereas Nuphar harbors two mitochondrial and one plastid phage-type polymerases. Identification of a plastid localized phage-type RNAP in this basal eudicot, orthologous to all other RpoTp enzymes of flowering plants, suggests that the acquisition of a nuclear encoded plastid RNA polymerase, not present in lycopods, took place after the split of lycopods from all other tracheophytes. A dual-targeted mitochondrial and plastid RNA polymerase (RpoTmp), as present in eudicots but not monocots, was not detected in Nuphar or Selaginella suggesting that its occurrence is an evolutionary novelty of eudicotyledoneous plants like Arabidopsis.

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.

Adaptive evolution of Toll-like receptors in birds Hana Velová, PhD thesis 6 Abstract Toll-like receptors (TLRs) are one of the key and presumably also evolutionary most original components of animal immune system. As Pattern recognition receptors they form the first line of innate immune defence against various pathogens. The proper receptor binding of pathogenic ligands is crucial for their correct recognition and for subsequent triggering of an appropriate immune response. Because there exists a direct interaction between the receptor surface and the pathogenic ligand, host-pathogen coevolution on molecular level can be predicted. Thus, through variability of their ligands, TLRs are exposed to extensive selective pressures that may be detected on both genetic and protein levels. Surprisingly, the variability we revealed in birds is even higher than previously expected based on the reports from other vertebrates, mainly mammals. In my doctoral thesis I summarise the results of my contribution to the avian TLR research. We were the first who experimentally verify the absence of functional TLR5 in several avian species and duplication of TLR7 in others. We finally resolved the origin of duplication in TLR1 and in TLR2 family. An important part of my research project focused on the prediction of potentially...

Gene tree-species tree reconciliation is the problem of mapping each node in a gene tree to a position in a species tree. Several methods have been used to address this problem. Variational inference is a method for finding the best approximation to the true distribution in a family of distributions. In this project, we investigated whether variational inference is a useful method to address the gene tree-species tree reconciliation problem. The distribution of trees is modeled by a so-called Subsplit Bayesian Network (SBN), and the evolution process is modeled by a birth-death process with constant duplication- and loss rate. We implemented the method in Python and compared it with A Variational Approach to Bayesian Phylogenetic Inference [1] (VBPI) [1] using synthetic data. The result showed that our method outperformed VBPI in most test cases. / Genträd-artträdsförsoning är problemet med att kartlägga varje nod i ett genträd till en position i ett artträd. Flera metoder har använts för att lösa detta problem. Variationsinferens är en metod för att hitta den bästa approximationen till den sanna fördelningen i en familj av sannolikhetsfördelningar. I det här projektet undersökte vi om variationsinferens är en användbar metod för att lösa Genträd-artträdsförsoningproblemet. Fördelningen av träd modelleras av ett så kallat subsplit Bayesian-nätverk (SBN), och evolutionsprocessen är modellerad av en födelse-dödsprocess med konstant duplicering- och förlusthastighet. Vi implementerade metoden i Python och jämförde den med VBPI [1] med syntetisk data. Resultatet visade att vår metod överträffade VBPI i de flesta testfallen.

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.