The Draft Genome Sequences of 50 Salamander species (Caudata, Amphibia)

We present partial genome sequences of 50 salamander species (Urodela) from 10 genera and 4 families. These span nearly the entire range of genome sizes in salamanders, from ~14-130GB, the latter of which is among the largest of all animal genomes. Only three salamander genomes were available to this point, from Ambystomatidae (one species) and Salamandridae (two species from two genera), to which we have added Amphiumidae (one species), Plethodontidae (45 species from 6 genera), Proteidae (one species), and Sirenidae (three species from two genera). These span ~140 million years of evolutionary divergence, leaving only Cryptobranchidae, Hynobiidae, and Rhyacotritonidae as salamander families without genome assemblies. These data should facilitate additional future work on speciation and genome evolution, both within Urodela and across Animalia.

The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Commensal protists in reptiles display flexible host range and adaptation to ectothermic hosts

IMPORTANCE

Environmental factors like climate change and captive breeding can impact the gut microbiota and host health. Therefore, conservation efforts for threatened species may benefit from understanding how these factors influence animal microbiomes. Parabasalid protists are members of the mammalian microbiota that can modulate the immune system and impact susceptibility to infections. However, little is known about parabasalids in reptiles. Here, we profile reptile-associated parabasalids in wild and captive reptiles and find that captivity has minimal impact on parabasalid prevalence or diversity. However, because reptiles are cold-blooded (ectothermic), their microbiotas experience wider temperature fluctuation than microbes in warm-blooded animals. To investigate whether extreme weather patterns affect parabasalid-host interactions, we analyzed the gene expression in reptile-associated parabasalids and found that temperature differences significantly alter genes associated with host health. These results expand our understanding of parabasalids in this vulnerable vertebrate group and highlight important factors to be taken into consideration for conservation efforts.

Commensal protists in reptiles display flexible host range and adaptation to ectothermic hosts

Parabasalid protists recently emerged as keystone members of the mammalian microbiota with important effects on their host's health. However, the prevalence and diversity of parabasalids in wild reptiles and the consequences of captivity and other environmental factors on these symbiotic protists are unknown. Reptiles are ectothermic, and their microbiomes are subject to temperature fluctuations, such as those driven by climate change. Thus, conservation efforts for threatened reptile species may benefit from understanding how shifts in temperature and captive breeding influence the microbiota, including parabasalids, to impact host fitness and disease susceptibility. Here, we surveyed intestinal parabasalids in a cohort of wild reptiles across three continents and compared these to captive animals. Reptiles harbor surprisingly few species of parabasalids compared to mammals, but these protists exhibited a flexible host-range, suggesting specific adaptations to reptilian social structures and microbiota transmission. Furthermore, reptile-associated parabasalids are adapted to wide temperature ranges, although colder temperatures significantly altered the protist transcriptomes, with increased expression of genes associated with detrimental interactions with the host. Our findings establish that parabasalids are widely distributed in the microbiota of wild and captive reptiles and highlight how these protists respond to temperature swings encountered in their ectothermic hosts.

Gut microbial ecology of Philippine gekkonids: ecoevolutionary effects on microbiome compositions

Given the rapidly changing landscapes of habitats across the globe, a sound understanding of host-associated microbial communities and the ecoevolutionary forces that shape them is needed to assess general organismal adaptability. Knowledge of the symbiotic endogenous microbiomes of most reptilian species worldwide remains limited. We sampled gut microbiomes of geckos spanning nine species and four genera in the Philippines to (i) provide baseline data on gut microbiota in these host species, (ii) test for significant associations between host phylogenetic relationships and observed microbial assemblages, potentially indicative of phylosymbiosis, and (iii) identify correlations between multiple ecoevolutionary factors (e.g. species identity, habitat tendencies, range extents, and maximum body sizes) and gut microbiomes in Philippine gekkonids. We recovered no significant association between interspecific host genetic distances and observed gut microbiomes, providing limited evidence for phylosymbiosis in this group. Philippine gekkonid microbiomes were associated most heavily with host species identity, though marked variation among conspecifics at distinct sampling sites indicates that host locality influences gut microbiomes as well. Interestingly, individuals grouped as widespread and microendemic regardless of host species identity displayed significant differences in alpha and beta diversity metrics examined, likely driven by differences in rare OTU presence between groups. These results provide much needed insight in host-associated microbiomes in wild reptiles and the ecoevolutionary forces that structure such communities.

Distribution of snake species and snakebites in hotspots of Ethiopia

INTRODUCTION

In Sub-Saharan Africa, snakebites are a public health problem. In Ethiopia, clinical cases have been described, but little information exists on snakebites burden and its geographical distribution. The aim of this study was to document the spatial distribution of venomous snakes and snakebites in Ethiopia.

METHODOLOGY

In a cross-sectional observational study, venomous snakes were collected during snake catching activities in six Ethiopian hotspot areas between April 2015 and September 2020. Species and habitat were described. In the hotspot areas, routine health information data on reported snakebites was collected in 78 districts and subsequently used to map annual incidence per district.

RESULTS

A total of 333 snakes were collected and 14 species were identified. The most prevalent species were Bitis arietans, Bitis arietanus somalica, Echis pyramidum, known as vipers, and Naja pallida, known as cobra. The highest number of snakes (75) was observed in the Northwest and Eastern parts of Ethiopia, mostly in cultivation and man-made farm land, wooded and moist dry savanna. In each hotspot a wide variety of species was observed, although composition was different. The highest snakebite incidence overlapped with the high snakes densities in Northwest Ethiopia. The snakebite annual average incidence at district level was very heterogeneous and ranged from < 15 cases/100,000 inhabitants (44% of the districts) to 309.2 cases/100,000 inhabitants.

CONCLUSIONS

Snake diversity and distribution, linked to high incidence of snakebites in the hotspots, suggests a close interconnection between human, animal and environmental systems and could inform the need for antivenoms per geographical locality.

A comprehensive phylogeny of dwarf geckos of the genus Lygodactylus, with insights into their systematics and morphological variation

The 71 currently known species of dwarf geckos of the genus Lygodactylus are a clade of biogeographic interest due to their occurrence in continental Africa, Madagascar, and South America. Furthermore, because many species are morphologically cryptic, our knowledge of species-level diversity within this genus is incomplete, as indicated by numerous unnamed genetic lineages revealed in previous molecular studies. Here we provide an extensive multigene phylogeny covering 56 of the named Lygodactylus species, four named subspecies, and 34 candidate species of which 19 are newly identified in this study. Phylogenetic analyses, based on ∼10.1 kbp concatenated sequences of eight nuclear-encoded and five mitochondrial gene fragments, confirm the monophyly of 14 Lygodactylus species groups, arranged in four major clades. We recover two clades splitting from basal nodes, one comprising exclusively Malagasy species groups, and the other containing three clades. In the latter, there is a clade with only Madagascar species, which is followed by a clade containing three African and one South American species groups, and its sister clade containing six African and two Malagasy species groups. Relationships among species groups within these latter clades remain weakly supported. We reconstruct a Lygodactylus timetree based on a novel fossil-dated phylotranscriptomic tree of squamates, in which we included data from two newly sequenced Lygodactylus transcriptomes. We estimate the crown diversification of Lygodactylus started at 46 mya, and the dispersal of Lygodactylus among the main landmasses in the Oligocene and Miocene, 35-22 mya, but emphasize the wide confidence intervals of these estimates. The phylogeny suggests an initial out-of-Madagascar dispersal as most parsimonious, but accounting for poorly resolved nodes, an out-of-Africa scenario may only require one extra dispersal step. More accurate inferences into the biogeographic history of these geckos will likely require broader sampling of related genera and phylogenomic approaches to provide better topological support. A survey of morphological characters revealed that most of the major clades and species groups within Lygodactylus cannot be unambiguously characterized by external morphology alone, neither by unique character states nor by a diagnostic combination of character states. Thus, any future taxonomic work will likely benefit from integrative, phylogenomic approaches.

Gut microbial ecology of the Critically Endangered Fijian crested iguana (Brachylophus vitiensis): Effects of captivity status and host reintroduction on endogenous microbiomes

Animals often exhibit distinct microbial communities when maintained in captivity as compared to when in the wild. Such differentiation may be significant in headstart and reintroduction programs where individuals spend some time in captivity before release into native habitats. Using 16S rRNA gene sequencing, we (i) assessed differences in gut microbial communities between captive and wild Fijian crested iguanas (Brachylophus vitiensis) and (ii) resampled gut microbiota in captive iguanas released onto a native island to monitor microbiome restructuring in the wild. We used both cloacal swabs and fecal samples to further increase our understanding of gut microbial ecology in this IUCN Critically Endangered species. We found significant differentiation in gut microbial community composition and structure between captive and wild iguanas in both sampling schemes. Approximately two months postrelease, microbial communities in cloacal samples from formerly captive iguanas closely resembled wild counterparts. Interestingly, microbial communities in fecal samples from these individuals remained significantly distinct from wild conspecifics. Our results indicate that captive upbringings can lead to differences in microbial assemblages in headstart iguanas as compared to wild individuals even after host reintroduction into native conditions. This investigation highlights the necessity of continuous monitoring of reintroduced animals in the wild to ensure successful acclimatization and release.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

Venomous Snakes Reveal Ecological and Phylogenetic Factors Influencing Variation in Gut and Oral Microbiomes

The gastrointestinal tract (GIT) of vertebrates contains a series of organs beginning with the mouth and ending with the anus or cloacal opening. Each organ represents a unique environment for resident microorganisms. Due to their simple digestive anatomy, snakes are good models for studying microbiome variation along the GIT. Cloacal sampling captures the majority of the microbial diversity found in the GIT of snakes—yet little is known about the oral microbiota of snakes. Most research on the snake mouth and gut microbiota are limited to studies of a single species or captive-bred individuals. It therefore remains unclear how a host’s life history, diet, or evolutionary history correlate with differences in the microbial composition within the mouths and guts of wild snakes. We sampled the mouth and gut microbial communities from three species of Asian venomous snakes and utilized 16S rRNA microbial inventories to test if host phylogenetic and ecological differences correlate with distinct microbial compositions within the two body sites. These species occupy three disparate habitat types: marine, semi-arboreal, and arboreal, our results suggest that the diversity of snake mouth and gut microbial communities correlate with differences in both host ecology and phylogeny.

The Tiger Rattlesnake genome reveals a complex genotype underlying a simple venom phenotype

Variation in gene regulation is ubiquitous, yet identifying the mechanisms producing such variation, especially for complex traits, is challenging. Snake venoms provide a model system for studying the phenotypic impacts of regulatory variation in complex traits because of their genetic tractability. Here, we sequence the genome of the Tiger Rattlesnake, which possesses the simplest and most toxic venom of any rattlesnake species, to determine whether the simple venom phenotype is the result of a simple genotype through gene loss or a complex genotype mediated through regulatory mechanisms. We generate the most contiguous snake-genome assembly to date and use this genome to show that gene loss, chromatin accessibility, and methylation levels all contribute to the production of the simplest, most toxic rattlesnake venom. We provide the most complete characterization of the venom gene-regulatory network to date and identify key mechanisms mediating phenotypic variation across a polygenic regulatory network.

The scales of coevolution: comparative phylogeography and genetic demography of a locally adapted venomous predator and its prey

Coevolutionary theory predicts that differences in the genetic demography of interacting species can influence patterns of local adaptation by affecting the potential of local populations to respond to selection. We conducted a comparative phylogeographical study of venomous rattlesnakes and their venom-resistant ground squirrel prey across California, and assessed how effective population size (Ne) estimates correspond with a previously documented pattern of rattlesnake local adaptation. Using RAD sequencing markers, we detected lineage relationships among both the rattlesnakes (Crotalus oreganus ssp.) and ground squirrels (Otospermophilus sp.) that are incongruent with previous phylogenetic hypotheses. Both rattlesnakes and squirrels share a deep divergence at the Sacramento–San Joaquin River Delta. At this broad phylogeographical scale, we found that the locally adapted rattlesnakes had higher Ne than squirrels. At the population scale, snakes also had larger Ne accompanied by larger values of several metrics of population genetic diversity. However, the specific magnitude of local adaptation of venom activity to ground squirrel venom resistance was not significantly correlated with local differences in Ne or other diversity statistics between predator and prey populations, suggesting that other factors in the geographic mosaic of coevolution contribute to the specific local-scale outcomes of this interaction. These results suggest an evolutionary mechanism that may explain some (but clearly not all) of rattlesnake local adaptation in this coevolutionary interaction – larger population sizes raise the adaptive potential of rattlesnakes compared to ground squirrels.

Gradual and Discrete Ontogenetic Shifts in Rattlesnake Venom Composition and Assessment of Hormonal and Ecological Correlates

Ontogenetic shifts in venom occur in many snakes but establishing their nature as gradual or discrete processes required additional study. We profiled shifts in venom expression from the neonate to adult sizes of two rattlesnake species, the eastern diamondback and the timber rattlesnake. We used serial sampling and venom chromatographic profiling to test if ontogenetic change occurs gradually or discretely. We found evidence for gradual shifts in overall venom composition in six of eight snakes, which sometimes spanned more than two years. Most chromatographic peaks shift gradually, but one quarter shift in a discrete fashion. Analysis of published diet data showed gradual shifts in overall diet composition across the range of body sizes attained by our eight study animals, while the shifts in abundance of different prey classes varied in form from gradual to discrete. Testosterone concentrations were correlated with the change in venom protein composition, but the relationship is not strong enough to suggest causation. Venom research employing simple juvenile versus adult size thresholds may be failing to account for continuous variation in venom composition lifespan. Our results imply that venom shifts represent adaptive matches to dietary shifts and highlight venom for studies of alternative gene regulatory mechanisms.

Congruence and Conflict in the Higher-Level Phylogenetics of Squamate Reptiles: An Expanded Phylogenomic Perspective

Genome-scale data have the potential to clarify phylogenetic relationships across the tree of life but have also revealed extensive gene tree conflict. This seeming paradox, whereby larger data sets both increase statistical confidence and uncover significant discordance, suggests that understanding sources of conflict is important for accurate reconstruction of evolutionary history. We explore this paradox in squamate reptiles, the vertebrate clade comprising lizards, snakes, and amphisbaenians. We collected an average of 5103 loci for 91 species of squamates that span higher-level diversity within the clade, which we augmented with publicly available sequences for an additional 17 taxa. Using a locus-by-locus approach, we evaluated support for alternative topologies at 17 contentious nodes in the phylogeny. We identified shared properties of conflicting loci, finding that rate and compositional heterogeneity drives discordance between gene trees and species tree and that conflicting loci rarely overlap across contentious nodes. Finally, by comparing our tests of nodal conflict to previous phylogenomic studies, we confidently resolve 9 of the 17 problematic nodes. We suggest this locus-by-locus and node-by-node approach can build consensus on which topological resolutions remain uncertain in phylogenomic studies of other contentious groups. [Anchored hybrid enrichment (AHE); gene tree conflict; molecular evolution; phylogenomic concordance; target capture; ultraconserved elements (UCE).].

Phylogenetic and spatial distribution of evolutionary diversification, isolation, and threat in turtles and crocodilians (non-avian archosauromorphs)

Background

The origin of turtles and crocodiles and their easily recognized body forms dates to the Triassic and Jurassic. Despite their long-term success, extant species diversity is low, and endangerment is extremely high compared to other terrestrial vertebrate groups, with ~ 65% of ~ 25 crocodilian and ~ 360 turtle species now threatened by exploitation and habitat loss. Here, we combine available molecular and morphological evidence with statistical and machine learning algorithms to present a phylogenetically informed, comprehensive assessment of diversification, threat status, and evolutionary distinctiveness of all extant species.

Results

In contrast to other terrestrial vertebrates and their own diversity in the fossil record, the recent extant lineages of turtles and crocodilians have not experienced any global mass extinctions or lineage-wide shifts in diversification rate or body-size evolution over time. We predict threat statuses for 114 as-yet unassessed or data-deficient species and identify a concentration of threatened turtles and crocodilians in South and Southeast Asia, western Africa, and the eastern Amazon. We find that unlike other terrestrial vertebrate groups, extinction risk increases with evolutionary distinctiveness: a disproportionate amount of phylogenetic diversity is concentrated in evolutionarily isolated, at-risk taxa, particularly those with small geographic ranges. Our findings highlight the important role of geographic determinants of extinction risk, particularly those resulting from anthropogenic habitat-disturbance, which affect species across body sizes and ecologies.

Conclusions

Extant turtles and crocodilians maintain unique, conserved morphologies which make them globally recognizable. Many species are threatened due to exploitation and global change. We use taxonomically complete, dated molecular phylogenies and various approaches to produce a comprehensive assessment of threat status and evolutionary distinctiveness of both groups. Neither group exhibits significant overall shifts in diversification rate or body-size evolution, or any signature of global mass extinctions in recent, extant lineages. However, the most evolutionarily distinct species tend to be the most threatened, and species richness and extinction risk are centered in areas of high anthropogenic disturbance, particularly South and Southeast Asia. Range size is the strongest predictor of threat, and a disproportionate amount of evolutionary diversity is at risk of imminent extinction.

Availability of Recently Described Ptychadena (Anura: Ptychadenidae) Nomina from Ethiopia

Frogs of the genus Ptychadena (Boulenger, 1917) have long been identified as harboring cryptic diversity and comprising numerous species-complexes (Largen 1997; Zimkus et al. 2017), and many authors have recognized the particularly high hidden richness in the Ethiopian highlands (Largen 1997; see Largen Spawls 2010 and refs. within). This cryptic diversity was confirmed by recent molecular studies (Freilich et al. 2014; Smith et al. 2017a, Reyes-Velasco et al. 2018). Those authors identified a congruent set of evolutionarily distinct candidate species using both mitochondrial and nuclear DNA, and described the geographic and ecological isolation of these species in detail (Freilich et al. 2014; Smith et al. 2017a).

A new species of African Snake-Eyed Skink (Scincidae: Panaspis) from Ethiopia

A recent molecular phylogenetic revision of the snake-eyed skinks (genus Panaspis Cope, 1868) uncovered extensive cryptic diversity, including several new species from throughout sub-Saharan Africa. Here, we describe one of these from Ethiopia as Panaspis annettesabinae sp. nov. This description is based on a previous molecular phylogeny and morphological, scalation, and coloration data collected from the type specimen. Phylogenetic analyses place the species alone in what we term the P. annettesabinae species group from Ethiopia. This group forms the sister lineage to a large southern African radiation and suggests a potential northern origin for much of the extant diversity of Panaspis. Many new taxa have recently been discovered in the genus and region, and there are several historical records of Panaspis from elsewhere in Ethiopia. Thus, we suggest that the range of this new species (known only from a single specimen at present) may be much larger, and that additional undescribed species may exist in northern sub-Saharan Africa.

Sexual Dichromatism Drives Diversification within a Major Radiation of African Amphibians

Theory predicts that sexually dimorphic traits under strong sexual selection, particularly those involved with intersexual signaling, can accelerate speciation and produce bursts of diversification. Sexual dichromatism (sexual dimorphism in color) is widely used as a proxy for sexual selection and is associated with rapid diversification in several animal groups, yet studies using phylogenetic comparative methods to explicitly test for an association between sexual dichromatism and diversification have produced conflicting results. Sexual dichromatism is rare in frogs, but it is both striking and prevalent in African reed frogs, a major component of the diverse frog radiation termed Afrobatrachia. In contrast to most other vertebrates, reed frogs display female-biased dichromatism in which females undergo color transformation, often resulting in more ornate coloration in females than in males. We produce a robust phylogeny of Afrobatrachia to investigate the evolutionary origins of sexual dichromatism in this radiation and examine whether the presence of dichromatism is associated with increased rates of net diversification. We find that sexual dichromatism evolved once within hyperoliids and was followed by numerous independent reversals to monochromatism. We detect significant diversification rate heterogeneity in Afrobatrachia and find that sexually dichromatic lineages have double the average net diversification rate of monochromatic lineages. By conducting trait simulations on our empirical phylogeny, we demonstrate that our inference of trait-dependent diversification is robust. Although sexual dichromatism in hyperoliid frogs is linked to their rapid diversification and supports macroevolutionary predictions of speciation by sexual selection, the function of dichromatism in reed frogs remains unclear. We propose that reed frogs are a compelling system for studying the roles of natural and sexual selection on the evolution of sexual dichromatism across micro- and macroevolutionary timescales.

The emerging field of venom-microbiomics for exploring venom as a microenvironment, and the corresponding Initiative for Venom Associated Microbes and Parasites (iVAMP)

Venom is a known source of novel antimicrobial natural products. The substantial, increasing number of these discoveries have unintentionally culminated in the misconception that venom and venom-producing glands are largely sterile environments. Culture-dependent and -independent studies on the microbial communities in venom microenvironments reveal the presence of archaea, algae, bacteria, fungi, protozoa, and viruses. Venom-centric microbiome studies are relatively sparse to date with the adaptive advantages that venom-associated microbes might offer to their hosts, or that hosts might provide to venom-associated microbes, remaining largely unknown. We highlight the potential for the discovery of venom microbiomes within the adaptive landscape of venom systems. The considerable number of convergently evolved venomous animals, juxtaposed with the comparatively few known studies to identify microbial communities in venom, provides new possibilities for both biodiversity and therapeutic discoveries. We present an evidence-based argument for integrating microbiology as part of venomics (i.e., venom-microbiomics) and introduce iVAMP, the Initiative for Venom Associated Microbes and Parasites (https://ivamp-consortium.github.io/), as a growing collaborative consortium. We express commitment to the diversity, inclusion and scientific collaboration among researchers interested in this emerging subdiscipline through expansion of the iVAMP consortium.

•

Venom-microbiome studies as an integrative field of venomics and microbiology.

•

Argument for multi-omics-based discovery through a microenvironment framework.

•

Introduction of a venom-microbiome research consortium (iVAMP).

Systematics of South American snail-eating snakes (Serpentes, Dipsadini), with the description of five new species from Ecuador and Peru

A molecular phylogeny of the Neotropical snail-eating snakes (tribe Dipsadini) is presented including 43 (24 for the first time) of the 77 species, sampled for both nuclear and mitochondrial genes. Morphological and phylogenetic support was found for four new species of Dipsas and one of Sibon, which are described here based on their unique combination of molecular, meristic, and color pattern characteristics. Sibynomorphus is designated as a junior subjective synonym of Dipsas. Dipsas latifrontalis and D. palmeri are resurrected from the synonymy of D. peruana. Dipsas latifasciata is transferred from the synonymy of D. peruana to the synonymy of D. palmeri. A new name, D. jamespetersi, is erected for the taxon currently known as Sibynomorphus petersi. Re-descriptions of D. latifrontalis and D. peruana are presented, as well as the first photographic voucher of an adult specimen of D. latifrontalis, along with photographs of all known Ecuadorian Dipsadini species. The first country record of D. variegata in Ecuador is provided and D. oligozonata removed from the list of Peruvian herpetofauna. With these changes, the number of Dipsadini reported in Ecuador increases to 22, 18 species of Dipsas and four of Sibon.

Correction to 'The role of climatic and geological events in generating diversity in Ethiopian grass frogs (genus Ptychadena)'

[This corrects the article DOI: 10.1098/rsos.170021.].

The role of climatic and geological events in generating diversity in Ethiopian grass frogs (genus Ptychadena)

Ethiopia is a world biodiversity hotspot and harbours levels of biotic endemism unmatched in the Horn of Africa, largely due to topographic-and thus habitat-complexity, which results from a very active geological and climatic history. Among Ethiopian vertebrate fauna, amphibians harbour the highest levels of endemism, making amphibians a compelling system for the exploration of the impacts of Ethiopia's complex abiotic history on biotic diversification. Grass frogs of the genus Ptychadena are notably diverse in Ethiopia, where they have undergone an evolutionary radiation. We used molecular data and expanded taxon sampling to test for cryptic diversity and to explore diversification patterns in both the highland radiation and two widespread lowland Ptychadena. Species delimitation results support the presence of nine highland species and four lowland species in our dataset, and divergence dating suggests that both geologic events and climatic fluctuations played a complex and confounded role in the diversification of Ptychadena in Ethiopia. We rectify the taxonomy of the endemic P. neumanni species complex, elevating one formally synonymized name and describing three novel taxa. Finally, we describe two novel lowland Ptychadena species that occur in Ethiopia and may be more broadly distributed.

A new large species of Bitis Gray, 1842 (Serpentes: Viperidae) from the Bale Mountains of Ethiopia

A new species of viperine viperid snake is described, Bitis harenna sp. nov. The new species is a member of the subgenus Macrocerastes based on it having three scales separating the nasal and rostral shields, and on the combination of 'divisions' of dorsal scale rows on the upper flanks and 'fusions' of rows on the lower flanks. Bitis harenna sp. nov. is distinguished from other members of the subgenus by its unique colour pattern, posterior parietal flange on the lateral wall of the braincase, and possibly by differences in scalation and head proportions. Only a single museum specimen is known, a female collected from 'Dodola' in Ethiopia probably in the late 1960s and previously identified as a possibly unusually coloured and patterned B. parviocula. A live, presumably male, specimen very closely resembling the holotype of Bitis harenna sp. nov. was photographed on the Harenna escarpment of the Bale Mountains National Park, Ethiopia in 2013, providing secure occurrence data and evidence that the holotype is not a uniquely aberrant specimen. A revised key to the species of Bitis in Ethiopia is presented. Aspects of body scalation are compared among species of the subgenus Macrocerastes and between species of Macrocerastes and Bitis, and several systematic characters are highlighted and clarified.

Amphibians and reptiles of the Calakmul Biosphere Reserve, México, with new records

We provide a list of amphibians and reptiles of the Calakmul Biosphere Reserve in the southern half of the Mexican Yucatan, in the state of Campeche. The study area was sampled through opportunistic, transect and pitfall trap surveys conducted for three successive years. These surveys resulted in a total of 2,359 amphibian and reptile encounters, belonging to 20 amphibian and 69 reptile species from 24 total families. We present herein the records for one snake, one chelonian and two salamander species not previously recorded in the area.

Phylogenetic Analysis of Bacterial Communities in Different Regions of the Gastrointestinal Tract of Agkistrodon piscivorus, the Cottonmouth Snake

Vertebrates are metagenomic organisms in that they are composed not only of their own genes but also those of their associated microbial cells. The majority of these associated microorganisms are found in the gastrointestinal tract (GIT) and presumably assist in processes such as energy and nutrient acquisition. Few studies have investigated the associated gut bacterial communities of non-mammalian vertebrates, and most rely on captive animals and/or fecal samples only. Here we investigate the gut bacterial community composition of a squamate reptile, the cottonmouth snake, Agkistrodon piscivorus through pyrosequencing of the bacterial 16S rRNA gene. We characterize the bacterial communities present in the small intestine, large intestine and cloaca. Many bacterial lineages present have been reported by other vertebrate gut community studies, but we also recovered unexpected bacteria that may be unique to squamate gut communities. Bacterial communities were not phylogenetically clustered according to GIT region, but there were statistically significant differences in community composition between regions. Additionally we demonstrate the utility of using cloacal swabs as a method for sampling snake gut bacterial communities.