Evolutionary biology is not a specialty, like genetics or development - it is an explanation of what is investigated by all biological specialties. Thus, the goal of this dissertation was to study both micro- and macroevolutionary processes in a multi-disciplinary framework. Population genetics, conservation, and phylogeny inference. The Jamaican boa (Epicrates subflavus) is an endemic species, whose natural populations greatly and constantly declined since the late 19th century, mainly due to predation by introduced species, human persecution, and habitat destruction. Using species-specific nuclear microsatellite loci and mitochondrial sequences, we investigated the population structure of this endangered reptile. All analyses pinpointed to an Eastern versus (Western+Central) pattern of differentiation in agreement with geological data and patterns of differentiation uncovered in other vertebrate and invertebrate Jamaican species. The same molecular markers were employed on 80 Jamaican boas of the European captive breeding program. This approach allowed us to (i) clarify all ambiguities in the studbook, (ii) correct parental allocation errors and (iii) assess the genetic diversity and the level of inbreeding of the current captive population. These results provide important insights for guiding the development of proper ex-situ and in-situ species survival and habitat management plans for this vulnerable snake. In the same framework of classical evolutionary genetics, we performed preliminary analyses of cytochrome b-like sequences in representatives of all cetacean families (but one), and revealed the presence of at least four nuclear mitochondrial pseudogenes that were independently inserted into the nuclear genome. Evo-Devo. The emergence of Evolutionary Developmental biology has caused a partial shift in the criteria for the selection of model species. Thus far, the main criterion was the relevance of a species for understanding human biology, whereas in the frame of the new discipline, it is the understanding of the generative mechanisms underlying biological diversity that is put forward. We discussed a few criteria and limitations of major relevance to the choice of model species for Evo-Devo studies, and applied a pragmatic approach to identify possible model species within Amniotes. Moreover, we developed MANTiS, an application pipeline that aims at integrating genomic, functional and expression data with evolutionary concepts, thus constituting the missing link between multi-species genome comparisons and functional analyses. Using MANTiS, we proceeded in the analysis of 35 metazoan full genomes for identifying all lineage-specific gene gains and losses. These results were combined with functional and expression analyses, and we demonstrated the much higher performance of MANTiS against popular databases of ortholog clusters (InParanoid, OrthoMCL, RoundUp). Finally, preliminary results of our attempt to adapt the new revolutionary technology of DNA sequencing in microfabricated high-density picoliter reactors (developed by 454/Roche) to the ultra-fast sequencing of brain full transcriptomes in multiple reptilian species are highly promising. As an example, the Crocodylus sample generated more than 72 Mbases (per run), which were successfully assembled in approximately 31,000 contigs. One third of the latter could be matched to known sequences in the transcriptome of related species. After fine-tuning of the in silico analyses, and incorporation of genomic sequence data, we expect our approach to provide important insights not only in the evolution of central nervous system novelties in vertebrates, but in transcriptomes in general as the brain transcriptome is one of the most complex among all organs.

<p>In the post-genomics era, the sheer volume of data is overwhelming without appropriate tools for data integration and analysis. Studying genomic sequences in the context of other related genomic sequences, i.e. comparative genomics, is a powerful technique enabling the identification of functionally interesting sequence regions based on the principal that similar sequences tend to be either homologous or provide similar functionality.</p> <p>Costs associated with full genome sequencing make it infeasible to sequence every genome of interest. Consequently, simple, smaller genomes are used as model organisms for more complex organisms, for instance, Mouse/Human. An annotated model organism provides a source of annotation for transcribed sequences and other gene regions of the more complex organism based on sequence homology. For example, the gene annotations from the model organism aid interpretation of expression studies in more complex organisms.</p> <p>To assist with comparative genomics research in the Arabidopsis/Brassica (Thale-cress/Canola) model-crop pair, a web-based, graphical genome browser (BioViz) was developed to display short Brassica genomic sequences in the context of the Arabidopsis model organism genome. This involved the development of graphical representations to integrate data from multiple sources and tools, and a novel user interface to provide the user with a more interactive web-based browsing experience. While BioViz was developed for the Arabidopsis/Brassica comparative genomics context, it could be applied to comparative browsing relative to other reference genomes.</p> <p>BioViz proved to be an valuable research support tool for Brassica / Arabidopsis comparative genomics. It provided convenient access to the underlying Arabidopsis annotation, allowed the user to view specific EST sequences in the context of the Arabidopsis genome and other related EST sequences. In addition, the limits to which the project pushed the SVG specification proved influential in the SVG community. The work done for BioViz inspired the definition of an opensource project to define standards for SVG based web applications and a standard framework for SVG based widget sets.</p>

Recent advances in DNA sequencing technology have enabled entire genomes to be sequenced quickly and accurately, resulting in an exponential increase in the number of organisms whose genome sequences have been elucidated. While the genome sequence of a given organism represents an important starting point in understanding its physiology, the functions of the protein products of many genes are still unknown; as such, computational methods for studying protein function are becoming increasingly important. In addition, this wealth of genomic information has created an unprecedented opportunity to compare the protein content of different organisms; among other applications, this can enable us to improve taxonomic classifications, to develop more accurate diagnostic tests for identifying particular bacteria, and to better understand protein content relationships in both closely-related and distantly-related organisms.<p> This thesis describes the design, evaluation, and use of a program called Proteome Subtraction and Intersection (PSI) that uses an idea called genome subtraction for discovering protein-phenotype relationships and for characterizing differences in protein content in groups of organisms. PSI takes as input a set of proteomes, as well as a partitioning of that set into a subset of "included" proteomes and a subset of "excluded" proteomes. Using reciprocal BLAST hits, PSI finds orthologous relationships among all the proteins in the proteomes from the original set, and then finds groups of orthologous proteins containing at least one orthologue from each of the proteomes in the "included" subset, and none from any of the proteomes in the "excluded" subset.<p> PSI is first applied to finding protein-phenotype relationships. By identifying proteins that are present in all sequenced isolates of the genus <i>Lactobacillus</i>, but not in the related bacterium <i>Pediococcus pentosaceus</i>, proteins are discovered that are likely to be responsible for the difference in cell shape between the lactobacilli and <i>P. pentosaceus</i>. In addition, proteins are identified that may be responsible for resistance to the antibiotic gatifloxacin in some lactic acid bacteria.<p> This thesis also explores the use of PSI for comparing protein content in groups of organisms. Based on the idea of genome subtraction, a novel metric is proposed for comparing the difference in protein content between two organisms. This metric is then used to create a phylogenetic tree for a large set of bacteria, which to the author's knowledge represents the largest phylogenetic tree created to date using protein content. In addition, PSI is used to find the proteomic cohesiveness of isolates of several bacterial species in order to support or refute their current taxonomic classifications.<p> Overall, PSI is a versatile tool with many interesting applications, and should become more and more valuable as additional genomic information becomes available.

Brachyspira species are bacteria that colonize the intestines of some mammalian and avian species with different degrees of pathogenicity. Brachyspira intermedia is a mild pig and bird pathogen with an unknown genomic sequence. In this project, we completed the genome of Brachyspira intermedia PWS/AT and did a comparative genomic analysis between B. intermedia PWS/AT and the already completed genomes of B. hyodysenteriae WA1, B. murdochii 56-150T and B. pilosicoli 95/1000. A table containing 15 classes of unique and shared genes was developed and analyzed in order to gain a better understanding of species-specific traits and clues behind the different degree of pathogenicity. Our result shows that genes are overall poorly annotated and further studies are of great importance for understanding different and shared properties. The largest number of unique features was found in B. intermedia and B. murdochii. B. hyodysenteriae and B. pilosicoli has most likely developed independently towards different biological niches and B. pilosicoli has undergone a major reductive evolution. One plasmid and six prophages were found in B. intermedia, where two of the phages appear to be capable of horizontal gene transfer. Further genome sequencing of more strains will probably increase the understanding of species-specific traits even more.

<p>In the post-genomics era, the sheer volume of data is overwhelming without appropriate tools for data integration and analysis. Studying genomic sequences in the context of other related genomic sequences, i.e. comparative genomics, is a powerful technique enabling the identification of functionally interesting sequence regions based on the principal that similar sequences tend to be either homologous or provide similar functionality.</p> <p>Costs associated with full genome sequencing make it infeasible to sequence every genome of interest. Consequently, simple, smaller genomes are used as model organisms for more complex organisms, for instance, Mouse/Human. An annotated model organism provides a source of annotation for transcribed sequences and other gene regions of the more complex organism based on sequence homology. For example, the gene annotations from the model organism aid interpretation of expression studies in more complex organisms.</p> <p>To assist with comparative genomics research in the Arabidopsis/Brassica (Thale-cress/Canola) model-crop pair, a web-based, graphical genome browser (BioViz) was developed to display short Brassica genomic sequences in the context of the Arabidopsis model organism genome. This involved the development of graphical representations to integrate data from multiple sources and tools, and a novel user interface to provide the user with a more interactive web-based browsing experience. While BioViz was developed for the Arabidopsis/Brassica comparative genomics context, it could be applied to comparative browsing relative to other reference genomes.</p> <p>BioViz proved to be an valuable research support tool for Brassica / Arabidopsis comparative genomics. It provided convenient access to the underlying Arabidopsis annotation, allowed the user to view specific EST sequences in the context of the Arabidopsis genome and other related EST sequences. In addition, the limits to which the project pushed the SVG specification proved influential in the SVG community. The work done for BioViz inspired the definition of an opensource project to define standards for SVG based web applications and a standard framework for SVG based widget sets.</p>

Recent advances in DNA sequencing technology have enabled entire genomes to be sequenced quickly and accurately, resulting in an exponential increase in the number of organisms whose genome sequences have been elucidated. While the genome sequence of a given organism represents an important starting point in understanding its physiology, the functions of the protein products of many genes are still unknown; as such, computational methods for studying protein function are becoming increasingly important. In addition, this wealth of genomic information has created an unprecedented opportunity to compare the protein content of different organisms; among other applications, this can enable us to improve taxonomic classifications, to develop more accurate diagnostic tests for identifying particular bacteria, and to better understand protein content relationships in both closely-related and distantly-related organisms.<p> This thesis describes the design, evaluation, and use of a program called Proteome Subtraction and Intersection (PSI) that uses an idea called genome subtraction for discovering protein-phenotype relationships and for characterizing differences in protein content in groups of organisms. PSI takes as input a set of proteomes, as well as a partitioning of that set into a subset of "included" proteomes and a subset of "excluded" proteomes. Using reciprocal BLAST hits, PSI finds orthologous relationships among all the proteins in the proteomes from the original set, and then finds groups of orthologous proteins containing at least one orthologue from each of the proteomes in the "included" subset, and none from any of the proteomes in the "excluded" subset.<p> PSI is first applied to finding protein-phenotype relationships. By identifying proteins that are present in all sequenced isolates of the genus <i>Lactobacillus</i>, but not in the related bacterium <i>Pediococcus pentosaceus</i>, proteins are discovered that are likely to be responsible for the difference in cell shape between the lactobacilli and <i>P. pentosaceus</i>. In addition, proteins are identified that may be responsible for resistance to the antibiotic gatifloxacin in some lactic acid bacteria.<p> This thesis also explores the use of PSI for comparing protein content in groups of organisms. Based on the idea of genome subtraction, a novel metric is proposed for comparing the difference in protein content between two organisms. This metric is then used to create a phylogenetic tree for a large set of bacteria, which to the author's knowledge represents the largest phylogenetic tree created to date using protein content. In addition, PSI is used to find the proteomic cohesiveness of isolates of several bacterial species in order to support or refute their current taxonomic classifications.<p> Overall, PSI is a versatile tool with many interesting applications, and should become more and more valuable as additional genomic information becomes available.

The rapid development of large-scale genomic sequencing has dramatically changed the field of genetics, in part through the development of comparative genomics. Fungal comparative genomics is particularly powerful given the large number of genomes currently available, their compact architecture, and their relative ease of genetic manipulation. Fungal comparative genomics was employed in this work to address two related questions. First, it was used along with computational thermodynamic methods to predict non-coding RNA (ncRNA) in Saccharomyces cerevisiae. Sets of positive and negative control genes were evaluated to determine the effect of window sizes and step sizes on the sensitivity of ncRNA identification. The approach was then applied to predict ncRNA genes on chromosome 6 of S. cerevisiae and S. bayanus. Northern blot analysis, rapid amplification of cDNA ends (RACE), and publicly available cDNA library data were used to test the predictions. Strong experimental evidence was accumulated for four new ncRNA genes. Potential structural elements in the 5' and 3' untranslated regions (UTRs) of six annotated protein-coding genes were also identified. This work shows that thermodynamic approaches, coupled with comparative genomics, are powerful tools for predicting structural ncRNA. Second, comparative genomic approaches were employed to identify a non-reciprocal transfer event from Cryptococcus neoformans var. grubii to var. neoformans ~2 million years ago involving a 14 gene (~40 kb) region. The majority of clinical and environmental var. neoformans strains from around the world contain this sequence obtained from var. grubii. The introgression event likely occurred via an incomplete inter-varietal sexual cycle creating a hybrid intermediate where mobile elements common to both lineages mediated the exchange. The subsequent duplication in laboratory strains of a fragment of this same genomic region supports evolutionary theories that instabilities in subtelomeric regions promote adaptive evolution through gene amplification and subsequent adaptation. These data indicate that DNA exchange between closely related sympatric varieties or species may be a recurrent theme in the evolution of fungal species. It further suggests that while evolutionary divergence is the primary force driving speciation, rare introgression events also play a potentially important role. / Dissertation

MicroRNAs are a large family of 21-26 nucleotide non-coding RNAs with a role in the post-transcriptional regulation of gene expression. In recent years, microRNAs have been proposed to play a significant role in the expansion of organism complexity. MicroRNAs are expressed in a cell or tissue-specific manner during embryonic development, suggesting a role in cellular differentiation. For example, Let-7 is a metazoan microRNA that acts as developmental timer between larval stages in C. elegans. We conducted a comparative study that determined the distribution of microRNA families among metazoans, including the identification of new family members for several species. MicroRNA families appear to have evolved in bursts of evolution that correlate with the advent of major metazoan groups such as vertebrates, eutherians, primates and hominids. Most microRNA families identified in these organisms appeared with or after the advent of vertebrates. Only a few of them appear to be shared between vertebrates and invertebrates. The distribution of these microRNA families supports the idea that at least one whole genome duplication event (WGS) predates the advent of vertebrates. Gene ontology analyses of the genes these microRNA families regulate show enrichments for functions related to cell differentiation and morphogenesis. MicroRNA genes appear to be under great selective constraints. Identification of conserved regions by comparative genomics allows for the computational identification of microRNAs. We have identified and characterized ultraconserved regions between the genomes of the honey bee (Apis mellifera) and the parasitic wasp (Nasonia vitripennis), and developed a strategy for the identification of microRNAs based on regions of ultraconservation. Ultraconserved regions preferentially localize within introns and intergenic regions, and are enriched in functions related to neural development. Introns harboring ultraconserved elements appear to be under negative selection and under a level of constraint that is higher than in their exonic counterparts. This level of constraint suggests functional roles yet to be discovered and suggests that introns are major players in the regulation of biological processes. Our computational strategy was able to identify new microRNA genes shared between honey bee and wasp. We recovered 41 of 45 previously validated microRNAs for these organisms, and we identified several new ones. A significant fraction of these microRNA candidates are located in introns and intergenic regions and are organized in genomic clusters. Expression of 13 of these new candidates was verified by 454 sequencing.

The major histocompatibility complex (MHC) is a genomic region containing genes of immunomodulatory importance. MHC class I genes encode cell-surface glycoproteins that present peptides to circulating T cells, playing a key role in recognition of self and non-self. Studies of MHC loci in vertebrates have examined levels of polymorphism and molecular evolutionary processes generating diversity. The bovine MHC (BoLA) has been associated with disease susceptibility, resistance and progression. To delineate mechanisms by which MHC class I genes evolved to function optimally in a species like cattle, it is necessary to study genomic organization of BoLA to define gene content, and investigate characteristics of expressed class I molecules. This study describes development of a physical map of BoLA class I region derived from screening two BAC libraries, isolating positive clones and confirming gene content, order and chromosomal location through PCR, novel BAC end sequencing techniques, and selected BAC shotgun cloning and/or sequencing and FISH analysis. To date, this is the most complete ordered BAC array encompassing the BoLA class I region from the class III boundary to the extended class I region. Characterization of a frameshift allele exhibiting trans-species polymorphism in Bos and Bison by flow cytometry, real-time RT-PCR, 1D and 2D gel analysis is also described. This frameshift allele encodes an early termination signal within the antigen recognition site (ARS) of exon 3 of the BoLA BSA-Adel class I gene predicting a truncated class I protein that is soluble. An ability to assess MHC diversity in populations and provision of animals with defined MHC haplotypes and genetic content for experimental research is necessary in developing a basis upon which to build functional studies to elucidate associations between haplotype and disease in bovids. The BoLA class I region is immunologically important for disease association studies in an economically important species. This study provides knowledge of gene content and organization within the class I MHC region in cattle, providing a template for more detailed analysis and elucidation of complex disease associations through functional genomics and comparative analysis, as well as evolution of the MHC in bovids to optimize a population’s immune response.

Chemokines triggered by Toll-like receptors (TLRs) are small chemoattractant proteins, which mainly regulate leukocyte trafficking in inflammatory reactions via interaction with G protein-coupled receptors. Forty-two chemokines and 19 cognate receptors have been found in the human genome. Prior to this study, only 11 chicken chemokines and 7 receptors had been reported. The objectives of this study were to identify systematically chicken chemokines and their cognate receptor genes in the chicken genome and to annotate these genes and ligand-receptor binding by a comparative genomics approach. Twenty-three chemokine and 14 chemokine receptor genes were identified in the chicken genome. The number of coding exons in these genes and the syntenies are highly conserved between human, mouse, and chicken although the amino acid sequence homologies are generally low between mammalian and chicken chemokines. Chicken genes were named with the systematic nomenclature used in humans and mice based on phylogeny, synteny, and sequence homology. The independent nomenclature of chicken chemokines and chemokine receptors suggests that the chicken may have ligand-receptor pairings similar to mammals. The TLR family represents evolutionarily conserved components of the patternrecognizing receptors (PRRs) of the innate immune system that recognize specific pathogen-associated molecular patterns (PAMPs) through their ectodomains (ECDs). TLR's ECDs contain 19 to 25 tandem copies of leucine-rich repeat (LRR) motifs. TLRs play important roles in the activation of pro-inflammatory cytokines, chemokines and modulation of antigen-specific adaptive immune responses. To date, nine TLRs have been reported in chicken, along with a non-functional TLR8. Two non-mammalian TLRs, TLR21 and TLR22, have been identified in pufferfish and zebrafish. The objectives of this study were to determine if there is the existence of chicken genes homologous to fish-specific TLRs, and if possible ligands of these receptors exist. After searching the chicken genome sequence and EST database, a novel chicken TLR homologous to fish TLR21 was identified. Phylogenetic analysis indicated that the identified chicken TLR is the orthologue of TLR21 in fish. Bioinformatic analysis of potential PAMP binding sites within LRR insertions showed that CpG DNA is the putative ligand of this receptor.