Jacobo Reyes-Velasco

The Burmese python genome reveals the molecular basis for extreme adaptation in snakes

Significance The molecular basis of morphological and physiological adaptations in snakes is largely unknown. Here, we study these phenotypes using the genome of the Burmese python (Python molurus bivittatus), a model for extreme phenotypic plasticity and metabolic adaptation. We discovered massive rapid changes in gene expression that coordinate major changes in organ size and function after feeding. Many significantly responsive genes are associated with metabolism, development, and mammalian diseases. A striking number of genes experienced positive selection in ancestral snakes. Such genes were related to metabolism, development, lungs, eyes, heart, kidney, and skeletal structure—all highly modified features in snakes. Snake phenotypic novelty seems to be driven by the system-wide coordination of protein adaptation, gene expression, and changes in genome structure. Snakes possess many extreme morphological and physiological adaptations. Identification of the molecular basis of these traits can provide novel understanding for vertebrate biology and medicine. Here, we study snake biology using the genome sequence of the Burmese python (Python molurus bivittatus), a model of extreme physiological and metabolic adaptation. We compare the python and king cobra genomes along with genomic samples from other snakes and perform transcriptome analysis to gain insights into the extreme phenotypes of the python. We discovered rapid and massive transcriptional responses in multiple organ systems that occur on feeding and coordinate major changes in organ size and function. Intriguingly, the homologs of these genes in humans are associated with metabolism, development, and pathology. We also found that many snake metabolic genes have undergone positive selection, which together with the rapid evolution of mitochondrial proteins, provides evidence for extensive adaptive redesign of snake metabolic pathways. Additional evidence for molecular adaptation and gene family expansions and contractions is associated with major physiological and phenotypic adaptations in snakes; genes involved are related to cell cycle, development, lungs, eyes, heart, intestine, and skeletal structure, including GRB2-associated binding protein 1, SSH, WNT16, and bone morphogenetic protein 7. Finally, changes in repetitive DNA content, guanine-cytosine isochore structure, and nucleotide substitution rates indicate major shifts in the structure and evolution of snake genomes compared with other amniotes. Phenotypic and physiological novelty in snakes seems to be driven by system-wide coordination of protein adaptation, gene expression, and changes in the structure of the genome.

Incipient speciation with biased gene flow between two lineages of the Western Diamondback Rattlesnake (Crotalus atrox)

We used mitochondrial DNA sequence data from 151 individuals to estimate population genetic structure across the range of the Western Diamondback Rattlesnake (Crotalus atrox), a widely distributed North American pit viper. We also tested hypotheses of population structure using double-digest restriction site associated DNA (ddRADseq) data, incorporating thousands of nuclear genome-wide SNPs from 42 individuals. We found strong mitochondrial support for a deep divergence between eastern and western C. atrox populations, and subsequent intermixing of these populations in the Inter-Pecos region of the United States and Mexico. Our nuclear RADseq data also identify these two distinct lineages of C. atrox, and provide evidence for nuclear admixture of eastern and western alleles across a broad geographic region. We identified contrasting patterns of mitochondrial and nuclear genetic variation across this genetic fusion zone that indicate partially restricted patterns of gene flow, which may be due to either pre- or post-zygotic isolating mechanisms. The failure of these two lineages to maintain complete genetic isolation, and evidence for partially-restricted gene flow, imply that these lineages were in the early stages of speciation prior to secondary contact.

Phylogenetic relationships of the enigmatic longtailed rattlesnakes (Crotalus ericsmithi, C. lannomi, and C. stejnegeri).

The longtailed rattlesnakes of western Mexico represent an enigmatic group of poorly known venomous snake species: Crotalus ericsmithi, C. lannomi, and C. stejnegeri. In the 120 years since their discovery, fewer than twenty individuals have been deposited in natural history collections worldwide. These three species share similar morphological traits, including a particularly long tail that has been interpreted as either an ancestral condition among rattlesnakes or as derived within the longtailed group. An understanding of the phylogenetic distinctiveness and relationships among the longtailed rattlesnakes, and their relationships to other rattlesnake groups, has previously been hampered by a dearth of comparative material and tissues for collection of DNA sequence data. Facilitated by the recent availability of tissue samples from multiple individuals of each species, we estimate the phylogenetic relationships among the longtailed rattlesnakes and their placement among other rattlesnake groups, using DNA sequence data from three mitochondrial and three nuclear gene fragments. We explore phylogenetic signal in our data using Bayesian and maximum likelihood methods, species tree analyses and hypothesis testing. Our results strongly support the monophyly of longtailed rattlesnakes and suggest the three species diverged from each other during the mid to late Pliocene or early Pleistocene (~1.5-5.6 mya). Contrary to prevailing hypotheses, we find no evidence for an early or basal divergence of the longtailed clade within the rattlesnake tree, and instead estimate that it diverged relatively recently (~6.8 mya) from its sister lineage, composed of the diamondback rattlesnakes (C. atrox group) and the prairie rattlesnakes (C. viridis group). With our added sampling of lineages and identification of previously used problematic sequences, we provide a revised hypothesis for relationships among Crotalus species, yet underscore the need for future studies and new data to deliver a well-supported robust estimate of rattlesnake relationships.

Two low coverage bird genomes and a comparison of reference-guided versus de novo genome assemblies

As a greater number and diversity of high-quality vertebrate reference genomes become available, it is increasingly feasible to use these references to guide new draft assemblies for related species. Reference-guided assembly approaches may substantially increase the contiguity and completeness of a new genome using only low levels of genome coverage that might otherwise be insufficient for de novo genome assembly. We used low-coverage (∼3.5–5.5x) Illumina paired-end sequencing to assemble draft genomes of two bird species (the Gunnison Sage-Grouse, Centrocercus minimus, and the Clarks Nutcracker, Nucifraga columbiana). We used these data to estimate de novo genome assemblies and reference-guided assemblies, and compared the information content and completeness of these assemblies by comparing CEGMA gene set representation, repeat element content, simple sequence repeat content, and GC isochore structure among assemblies. Our results demonstrate that even lower-coverage genome sequencing projects are capable of producing informative and useful genomic resources, particularly through the use of reference-guided assemblies.

Molecular systematics of the genus Sonora (Squamata: Colubridae) in central and western Mexico

Mexico possesses high levels of endemic biodiversity, especially for squamate reptiles. However, the evolutionary relationships among many reptiles in this region are not well known. The closely related genera of Sonora Baird and Girard 1853 and Procinura Cope 1879 are coralsnake mimics found from the central and western United States to southwestern Mexico and Baja California. Although species delimitation in this group has historically relied upon colour pattern and other morphological characters, many populations of these species display colour pattern polymorphism, which may confound taxonomy. We used molecular phylogenetics to assess the evolutionary relationships and delimit species within Sonora, focusing on the phylogenetic position of Procinura and the validity of S. mutabilis and aequalis. We sequenced two mitochondrial (ND4 and cytb) and two nuclear (c-mos and RAG-1) genes for the single species of Procinura and each of the four species of Sonora. We analysed these sequences using maximum parsimony, maximum likelihood and Bayesian phylogenetic analyses on separately concatenated mitochondrial and nuclear datasets. Additionally, we used Bayesian coalescent methods to build a species tree (Bayesian species tree analysis) and delimit species boundaries (Bayesian species delimitation). All methods indicated that Procinura is deeply nested within Sonora, and most individual species are well supported. However, we found that one taxon (S. aequalis) is paraphyletic with regard to another (S. mutabilis). We recommend that the genus Procinura be synonymised with Sonora and that S. aequalis be synonymised with S. mutabilis. Additionally, the phylogenetic patterns that we document are broadly congruent with a Miocene or Pliocene divergence between S. michoacanensis and S. mutabilis along the Trans-Mexican Volcanic Belt. Finally, our data are consistent with the early evolution of coralsnake mimicry and colour pattern polymorphism within the genus Sonora.

Microsatellite landscape evolutionary dynamics across 450 million years of vertebrate genome evolution

The evolutionary dynamics of simple sequence repeats (SSRs or microsatellites) across the vertebrate tree of life remain largely undocumented and poorly understood. In this study, we analyzed patterns of genomic microsatellite abundance and evolution across 71 vertebrate genomes. The highest abundances of microsatellites exist in the genomes of ray-finned fishes, squamate reptiles, and mammals, while crocodilian, turtle, and avian genomes exhibit reduced microsatellite landscapes. We used comparative methods to infer evolutionary rates of change in microsatellite abundance across vertebrates and to highlight particular lineages that have experienced unusually high or low rates of change in genomic microsatellite abundance. Overall, most variation in microsatellite content, abundance, and evolutionary rate is observed among major lineages of reptiles, yet we found that several deeply divergent clades (i.e., squamate reptiles and mammals) contained relatively similar genomic microsatellite compositions. Archosauromorph reptiles (turtles, crocodilians, and birds) exhibit reduced genomic microsatellite content and the slowest rates of microsatellite evolution, in contrast to squamate reptile genomes that have among the highest rates of microsatellite evolution. Substantial branch-specific shifts in SSR content in primates, monotremes, rodents, snakes, and fish are also evident. Collectively, our results support multiple major shifts in microsatellite genomic landscapes among vertebrates.

Revisiting the phylogeography, demography and taxonomy of the frog genus Ptychadena in the Ethiopian highlands with the use of genome-wide SNP data

Understanding the diversification of biological lineages is central to evolutionary studies. To properly study the process of speciation, it is necessary to link micro-evolutionary studies with macro-evolutionary mechanisms. Micro-evolutionary studies require proper sampling across a taxon’s range to adequately infer genetic diversity. Here we use the grass frogs of the genus Ptychadena from the Ethiopian highlands as a model to study the process of lineage diversification in this unique biodiversity hotspot. We used thousands of genome-wide SNPs obtained from double digest restriction site associated DNA sequencing (ddRAD-seq) in populations of the Ptychadena neumanni species complex from the Ethiopian highlands in order to infer their phylogenetic relationships and genetic structure, as well as to study their demographic history. Our genome-wide phylogenetic study supports the existence of approximately 13 lineages clustered into 3 species groups. Our phylogenetic and phylogeographic reconstructions suggest that those endemic lineages diversified in allopatry, and subsequently specialized to different habitats and elevations. Demographic analyses point to a continuous decrease in the population size across the majority of lineages and populations during the Pleistocene, which is consistent with a continuous period of aridification that East Africa experienced since the Pliocene. We discuss the taxonomic implications of our analyses and, in particular, we warn against the recent practice to solely use Bayesian species delimitation methods when proposing taxonomic changes.

Synopsis and taxonomic revision of three genera in the snake tribe Sonorini

ABSTRACT Delimiting species is a crucial goal of integrative biology, and yet can be misled by homoplasy and high levels of morphological variation. The snake tribe Sonorini contains three genera that have long confounded taxonomists: Chilomeniscus, Chionactis and Sonora. Dynamic colour evolution in this group, including rampant geographic variation in colour and colour polymorphism, has led to a chaotic taxonomy. We used mitochondrial and high-throughput nuclear data (ddRADseq) and complete taxonomic sampling of each genus to reconstruct phylogenetic relationships and systematically revise the genus. Our research revealed that Sonora is paraphyletic with regards to Chilomeniscus and Chionactis and that at least one species (S. semiannulata) is paraphyletic with respect to at least one other recognized species. Additionally, we found substantial undescribed genetic diversity within multiple species which is incongruent with morphological variation in coloration. Accordingly, we proposed synonymizing Chionactis and Chilomeniscus with Sonora, which has taxonomic priority over both genera. As we found genetic evidence that supported some of the historically delimited diversity within multiple taxa, we revised species-level taxonomy accordingly. This new taxonomy recognizes a revised genus of Sonora that contains 15 species of diminutive and often brightly coloured snakes that are distributed from central Mexico to north-western USA. http://www.zoobank.org/urn:lsid:zoobank.org:pub:45A553D8-6435-4E0A-84ED-DF31E2CCD872

Two new species of Eleutherodactylus (subgenus Syrrhophus) from western Mexico.

We describe two new species of Eleutherodactylus, subgenus Syrrhophus, from two separate mountain ranges in western Mexico. Eleutherodactylus grunwaldi sp. nov. inhabits the Sierra de Manantlán in Colima and Jalisco from 1300 to 2200 m, whereas E. wixarika sp. nov. is known from a single locality in the Sierra Huichola of northern Jalisco at 2400 m, but is probably more widespread. Eleutherodactylus grunwaldi is readily distinguishable from most members of mainland Syrrhophus by a combination of its large size, broad, truncate digital pads more than three times the narrowest part of the digit, and a black and green marbled color pattern. This species is saxicolous, inhabiting limestone outcrops, and has been found in caves during the dry season. Eleutherodactylus wixarika is a moderate sized species, most similar to E. teretistes, E. pallidus and E. modestus. It is distinguished from all other members of the subgenus by the combination a tuberculate, reddish dorsum, lack of compact lumbar glands, and expanded digital pads less than twice the width of the narrowest part of the digit. This species inhabits areas with secondary vegetation in pine forest. Males of both species call at night during the rainy season. The advertisement call of both species consists of a short, narrow band, pure-tone note organized into a discrete train at a rate of about six times per minute. Spectral and temporal acoustic properties differ between species. The subgenus Syrrhophus of the genus Eleutherodactylus is one of the most poorly studied groups of frogs in Mexico but probably one of the most diverse.

The Conservation Status of Several Endemic Mexican Kinosternid Turtles

Abstract The goal of this project was to carry out surveys in northern and central México to locate extant populations of Kinosternon hirtipes chapalaense, K. h. magdalense, and K. h. megacephalum, evaluate their conservation status, help establish captive breeding colonies of these taxa in the near future, and collect tissues for phylogenetic studies. During 2010 and 2011, we were able to locate extant populations of 2 of the subspecies (K. h. chapalaense and K. h. magdalense), but K. h. megacephalum could not be found. We also failed to locate natural springs in a radius of approximately 60 km from the type locality (and only known locality) of K. h. megacephalum and believe that this turtle is extinct.