The Discovery of XY Sex Chromosomes in a \u3cem\u3eBoa\u3c/em\u3e and \u3cem\u3ePython\u3c/em\u3e

For over 50 years, biologists have accepted that all extant snakes share the same ZW sex chromosomes derived from a common ancestor [1, 2, 3], with different species exhibiting sex chromosomes at varying stages of differentiation. Accordingly, snakes have been a well-studied model for sex chromosome evolution in animals [1, 4]. A review of the literature, however, reveals no compelling support that boas and pythons possess ZW sex chromosomes [2, 5]. Furthermore, phylogenetic patterns of facultative parthenogenesis in snakes and a sex-linked color mutation in the ball python (Python regius) are best explained by boas and pythons possessing an XY sex chromosome system [6, 7]. Here we demonstrate that a boa (Boa imperator) and python (Python bivittatus) indeed possess XY sex chromosomes, based on the discovery of male-specific genetic markers in both species. We use these markers, along with transcriptomic and genomic data, to identify distinct sex chromosomes in boas and pythons, demonstrating that XY systems evolved independently in each lineage. This discovery highlights the dynamic evolution of vertebrate sex chromosomes and further enhances the value of snakes as a model for studying sex chromosome evolution

An ancient adaptive episode of convergent molecular evolution confounds phylogenetic inference

Convergence can mislead phylogenetic inference by mimicking shared ancestry, but has been detected only rarely in molecular evolution. Here, we show that significant convergence occurred in snake and agamid lizard mitochondrial genomes. Most evidence, and most of the mitochondrial genome, supports one phylogenetic tree, but a subset of mostly amino acid-altering mitochondrial sites strongly support a radically different phylogeny. These sites are convergent, probably selected, and overwhelm the signal from other sites. This suggests that convergent molecular evolution can seriously mislead phylogenetics, even with large data sets. Radical phylogenies inconsistent with previous evidence should be treated cautiously

Genome-wide signatures of convergent evolution in echolocating mammals

Evolution is typically thought to proceed through divergence of genes, proteins, and ultimately phenotypes(1-3). However, similar traits might also evolve convergently in unrelated taxa due to similar selection pressures(4,5). Adaptive phenotypic convergence is widespread in nature, and recent results from a handful of genes have suggested that this phenomenon is powerful enough to also drive recurrent evolution at the sequence level(6-9). Where homoplasious substitutions do occur these have long been considered the result of neutral processes. However, recent studies have demonstrated that adaptive convergent sequence evolution can be detected in vertebrates using statistical methods that model parallel evolution(9,10) although the extent to which sequence convergence between genera occurs across genomes is unknown. Here we analyse genomic sequence data in mammals that have independently evolved echolocation and show for the first time that convergence is not a rare process restricted to a handful of loci but is instead widespread, continuously distributed and commonly driven by natural selection acting on a small number of sites per locus. Systematic analyses of convergent sequence evolution in 805,053 amino acids within 2,326 orthologous coding gene sequences compared across 22 mammals (including four new bat genomes) revealed signatures consistent with convergence in nearly 200 loci. Strong and significant support for convergence among bats and the dolphin was seen in numerous genes linked to hearing or deafness, consistent with an involvement in echolocation. Surprisingly we also found convergence in many genes linked to vision: the convergent signal of many sensory genes was robustly correlated with the strength of natural selection. This first attempt to detect genome-wide convergent sequence evolution across divergent taxa reveals the phenomenon to be much more pervasive than previously recognised

A Proposal to Sequence the Genome of a Garter Snake (Thamnophis sirtalis)

Here we develop an argument in support of sequencing a garter snake (Thamnophis sirtalis) genome, and outline a plan to accomplish this. This snake is a common, widespread, nonvenomous North American species that has served as a model for diverse studies in evolutionary biology, physiology, genomics, behavior and coevolution. The anole lizard is currently the only genome sequence available for a non-avian reptile. Thus, the garter snake at this time would be the first available snake genome sequence and as such would provide much needed comparative representation of non-avian reptilian genomes, and would also allow critical new insights for vertebrate comparative genomic studies. We outline the major areas of discovery that the availability of the garter snake genome would enable, and describe a plan for whole-genome sequencing.Organismic and Evolutionary Biolog

Microsatellite discovery in an insular amphibian (Grandisonia alternans) with comments on cross-species utility and the accuracy of locus identification from unassembled Illumina data

The Seychelles archipelago is unique among isolated oceanic islands because it features an endemic radiation of caecilian amphibians (Gymnophiona). In order to develop population genetics resources for this system, we identified microsatellite loci using unassembled Illumina MiSeq data generated from a genomic library of Grandisonia alternans, a species that occurs on multiple islands in the archipelago. Applying a recently described method (PALFINDER) we identified 8001 microsatellite loci that were potentially informative for population genetics analyses. Of these markers, we screened 60 loci using five individuals, directly sequenced several amplicons to confirm their identity, and then used eight loci to score allele sizes in 64 G. alternans individuals originating from five islands. A number of these individuals were sampled using non-lethal methods, demonstrating the efficacy of non-destructive molecular sampling in amphibian research. Although two loci satisfied our criteria as diploid, neutrally evolving loci with the statistical power to detect population structure, our success in identifying reliable loci was very low. Additionally, we discovered some issues with primer redundancy and differences between Illumina and Sanger sequences that suggest some Illumina-inferred loci are invalid. We investigated cross-species utility for eight loci and found most could be successfully amplified, sequenced and aligned across other species and genera of caecilians from the Seychelles. Thus, our study in part supported the validity of using PALFINDER with unassembled reads for microsatellite discovery within and across species, but importantly identified major limitations to applying this approach to small datasets (ca. 1 million reads) and loci with small tandem repeat sizes

Ecdysozoan mitogenomics: evidence for a common origin of the legged invertebrates, the Panarthropoda

Ecdysozoa is the recently recognized clade of molting animals that comprises the vast majority of extant animal species and the most important invertebrate model organisms—the fruit fly and the nematode worm. Evolutionary relationships within the ecdysozoans remain, however, unresolved, impairing the correct interpretation of comparative genomic studies. In particular, the affinities of the three Panarthropoda phyla (Arthropoda, Onychophora, and Tardigrada) and the position of Myriapoda within Arthropoda (Mandibulata vs. Myriochelata hypothesis) are among the most contentious issues in animal phylogenetics. To elucidate these relationships, we have determined and analyzed complete or nearly complete mitochondrial genome sequences of two Tardigrada, Hypsibius dujardini and Thulinia sp. (the first genomes to date for this phylum); one Priapulida, Halicryptus spinulosus; and two Onychophora, Peripatoides sp. and Epiperipatus biolleyi; and a partial mitochondrial genome sequence of the Onychophora Euperipatoides kanagrensis. Tardigrada mitochondrial genomes resemble those of the arthropods in term of the gene order and strand asymmetry, whereas Onychophora genomes are characterized by numerous gene order rearrangements and strand asymmetry variations. In addition, Onychophora genomes are extremely enriched in A and T nucleotides, whereas Priapulida and Tardigrada are more balanced. Phylogenetic analyses based on concatenated amino acid coding sequences support a monophyletic origin of the Ecdysozoa and the position of Priapulida as the sister group of a monophyletic Panarthropoda (Tardigrada plus Onychophora plus Arthropoda). The position of Tardigrada is more problematic, most likely because of long branch attraction (LBA). However, experiments designed to reduce LBA suggest that the most likely placement of Tardigrada is as a sister group of Onychophora. The same analyses also recover monophyly of traditionally recognized arthropod lineages such as Arachnida and of the highly debated clade Mandibulata

Modeling And Partitioning The Nucleotide Evolutionary Process For Phylogenetic And Comparative Genomic Inference

The transformation of genomic data into functionally relevant information about the composition of biological systems hinges critically on the field of computational genome biology, at the core of which lies comparative genomics. The aim of comparative genomics is to extract meaningful functional information from the differences and similarities observed across genomes of different organisms. We develop and test a novel framework for applying complex models of nucleotide evolution to solve phylogenetic and comparative genomic problems, and demonstrate that these techniques are crucial for accurate comparative evolutionary inferences. Additionally, we conduct an exploratory study using vertebrate mitochondrial genomes as a model to identify the reciprocal influences that genome structure, nucleotide evolution, and multi-level molecular function may have on one another. Collectively this work represents a significant and novel contribution to accurately modeling and characterizing patterns of nucleotide evolution, a contribution that enables the enhanced detection of patterns of genealogical relationships, selection, and function in comparative genomic datasets. Our work with entire mitochondrial genomes highlights a coordinated evolutionary shift that simultaneously altered genome architecture, replication, nucleotide evolution and molecular function (of proteins, RNAs, and the genome itself). Current research in computational biology, including the advances included in this dissertation, continue to close the gap that impedes the transformation of genomic data into powerful tools for the analysis and understanding of biological systems function

Historic genetic structuring and paraphyly within the Great-tailed Grackle

The Great-tailed Grackle (Quiscalus mexicanus) and Boat-tailed Grackle (Q. major) are sister species that have expanded their ranges during historical times. This expansion has created an area of sympatry between these species in Texas and Louisiana, and between distinctive Great-tailed Grackle subspecies in the southwestern United States and northern Mexico. We investigated the evolutionary histories of both species using mitochondrial DNA sequence data and modern phylogenetic methods. Our results reveal genetic structure within Great-tailed, but not Boat-tailed Grackles. Great-tailed Grackles are separated into two clades, but range expansion in the north has led to secondary contact between them. Boat-tailed Grackles are monophyletic and are embedded within the Great-tailed Grackle assemblage, rendering the latter paraphyletic. These results reveal a complex phylogeographic pattern caused by recent range expansion and secondary contact of once allopatric units

Report from the First Snake Genomics and Integrative Biology Meeting

This report summarizes the proceedings of the 1st Snake Genomics and Integrative Biology Meeting held in Vail, CO USA, 5-8 October 2011. The meeting had over twenty registered participants, and was conducted as a single session of presentations. Goals of the meeting included coordination of genomic data collection and fostering collaborative interactions among researchers using snakes as model systems.Organismic and Evolutionary Biolog

Reorganisation ofHoxdregulatory landscapes during the evolution of a snake-like body plan

Within land vertebrate species, snakes display extreme variations in their body plan, characterized by the absence of limbs and an elongated morphology. Such a particular interpretation of the basic vertebrate body architecture has often been associated with changes in the function or regulation of Hox genes. Here, we use an interspecies comparative approach to investigate different regulatory aspects at the snake HoxD locus. We report that, unlike in other vertebrates, snake mesoderm-specific enhancers are mostly located within the HoxD cluster itself rather than outside. In addition, despite both the absence of limbs and an altered Hoxd gene regulation in external genitalia, the limb-associated bimodal HoxD chromatin structure is maintained at the snake locus. Finally, we show that snake and mouse orthologous enhancer sequences can display distinct expression specificities. These results show that vertebrate morphological evolution likely involved extensive reorganisation at Hox loci, yet within a generally conserved regulatory framework.Fundação para a Ciência e Tecnologia grant: (PTDB/BEX-BID/0899/2014); Schweizerischer Nationalfonds zur Forderung der Wissenschaftlichen Forschung grant: (310030B_138662); Claraz Foundation; Université de Genève; Instituto Federal de Educação Ciência e Tecnologia do Espírito Santo