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 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.

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.

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).].

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.

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.

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.

A species-level phylogeny of Trachylepis (Scincidae: Mabuyinae) provides insight into their reproductive mode evolution.

Trachylepis (Mabuyinae) includes ∼80 species of fully-limbed skinks found primarily in Africa and Madagascar, but a robust species-level phylogeny for this genus is lacking and this impedes studies on a wide-range of topics from biogeography to character evolution. Trachylepis and its close relatives (which together form the Mabuya group or Mabuyinae) are notable in that they have undergone multiple transitions and remarkable specializations in their reproductive modes. A Trachylepis phylogeny will be particularly useful for investigating reproductive evolution, because it includes species that exhibit oviparity, viviparity, and bimodal parity (species with both oviparous and viviparous populations). We sequenced DNA at four mitochondrial and five nuclear loci for 67 (∼84% of) Trachylepis species to infer a phylogeny for this genus. We performed stochastic character mapping of parity mode under a variety of parity mode transition models to infer ancestral parity mode states and the number and type of parity mode transitions. We recovered a strongly supported phylogeny of Trachylepis that is generally consistent with earlier phylogenetic studies. The best-fit model of reproductive mode evolution supports an oviparous ancestor for Trachylepis, and supports at least three viviparity to oviparity transitions. We compared parity mode evolution under the overall best-fit model (no constraints on parity mode transitions) to the best-fit model among the subset of models that assume viviparity to oviparity transitions are impossible. Our results support a model of reproductive evolution that allows for reversibility from viviparity to oviparity, a process that is not generally accepted. Alternatively, the best-fit model of evolution among the set of models that eliminate reversals from viviparity to oviparity suggests that bimodal reproduction may have persisted for millions of years within multiple lineages.

Hidden diversity within the depauperate genera of the snake tribe Lampropeltini (Serpentes, Colubridae).

Accurate representation of lineage diversity through complete taxon sampling is crucial to understanding the evolution of biodiversity, particularly when using molecular phylogenetics to estimate evolutionary relationships. In this interest, taxonomic diversity is often used as a proxy for lineage diversity even though the two concepts are not synonymous. We explore this within the snake tribe Lampropeltini which includes some of the most conspicuous and heavily studied snakes in North America. Both the taxonomy and hypothesized relationships within this tribe have been in flux. The number of species has increased from 23 to 51 over the last thirty years, predominately within three of the nine genera (Lampropeltis, Pantherophis, Pituophis). The remaining six depauperate genera (Arizona, Bogertophis, Cemophora, Pseudelaphe, Rhinocheilus, and Senticolis) have been poorly represented in phylogenetic studies. To estimate evolutionary relationships and determine if the dichotomy in depauperate and speciose genera within Lampropeltini is a function of taxon sampling or truly represents the lineage diversity, we estimated the phylogeny of this group using nuclear and mitochondrial loci in a concatenated and coalescent framework with the largest sampling of the six depauperate genera to date. In addition, we estimated the divergence dates among the genera to assess whether the instability of Lampropeltini phylogenetic relationships is due to an adaptive radiation. While some nodes still remain unresolved, the generic-level relationships we recovered agree with those of a recent next-generation study that used a much larger set of loci for fewer individuals. We also tested two putative species, Arizona pacata and Pseudelaphe phaescens, for the first time phylogenetically and find evidence that they are distinct lineages. Overall, we find that the taxonomic and genetic diversity are not correlated in Lampropeltini and that representing putative diversity in phylogenies will lead to a better estimate of evolutionary histories, especially in groups with complex radiations.

Microbiome evolution along divergent branches of the vertebrate tree of life: what is known and unknown.

Vertebrates harbour microbes both internally and externally, and collectively, these microorganisms (the 'microbiome') contain genes that outnumber the host's genetic information 10-fold. The majority of the microorganisms associated with vertebrates are found within the gut, where they influence host physiology, immunity and development. The development of next-generation sequencing has led to a surge in effort to characterize the microbiomes of various vertebrate hosts, a necessary first step to determine the functional role these communities play in host evolution or ecology. This shift away from a culture-based microbiological approach, limited in taxonomic breadth, has resulted in the emergence of patterns suggesting a core vertebrate microbiome dominated by members of the bacterial phyla Bacteroidetes, Proteobacteria and Firmicutes. Still, there is a substantial variation in the methodology used to characterize the microbiome, from differences in sample type to issues of sampling captive or wild hosts, and the majority (>90%) of studies have characterized the microbiome of mammals, which represent just 8% of described vertebrate species. Here, we review the state of microbiome studies of nonmammalian vertebrates and provide a synthesis of emerging patterns in the microbiome of those organisms. We highlight the importance of collection methods, and the need for greater taxonomic sampling of natural rather than captive hosts, a shift in approach that is needed to draw ecologically and evolutionarily relevant inferences. Finally, we recommend future directions for vertebrate microbiome research, so that attempts can be made to determine the role that microbial communities play in vertebrate biology and evolution.

Molecular phylogeny of Panaspis and Afroablepharus skinks (Squamata: Scincidae) in the savannas of sub-Saharan Africa.

African snake-eyed skinks are relatively small lizards of the genera Panaspis and Afroablepharus. Species allocation of these genera frequently changed during the 20th century based on morphology, ecology, and biogeography. Members of these genera occur primarily in savanna habitats throughout sub-Saharan Africa and include species whose highly conserved morphology poses challenges for taxonomic studies. We sequenced two mitochondrial (16S and cyt b) and two nuclear genes (PDC and RAG1) from 76 Panaspis and Afroablepharus samples from across eastern, central, and southern Africa. Concatenated gene-tree and divergence-dating analyses were conducted to infer phylogenies and biogeographic patterns. Molecular data sets revealed several cryptic lineages, with most radiations occurring during the mid-Miocene to Pliocene. We infer that rifting processes (including the formation of the East African Rift System) and climatic oscillations contributed to the expansion and contraction of savannas, and caused cladogenesis in snake-eyed skinks. Species in Panaspis and Afroablepharus used in this study, including type species for both genera, formed a monophyletic group. As a result, the latter genus should be synonymized with the former, which has priority. Conservatively, we continue to include the West African species P. breviceps and P. togoensis within an expanded Panaspis, but note that they occur in relatively divergent clades, and their taxonomic status may change with improved taxon sampling. Divergence estimates and cryptic speciation patterns of snake-eyed skinks were consistent with previous studies of other savanna vertebrate lineages from the same areas examined in this study.

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 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.

Molecular systematics and historical biogeography of tree boas (Corallus spp.).

Inferring the evolutionary and biogeographic history of taxa occurring in a particular region is one way to determine the processes by which the biodiversity of that region originated. Tree boas of the genus Corallus are an ancient clade and occur throughout Central and South America and the Lesser Antilles, making it an excellent group for investigating Neotropical biogeography. Using sequenced portions of two mitochondrial and three nuclear loci for individuals of all recognized species of Corallus, we infer phylogenetic relationships, present the first molecular analysis of the phylogenetic placement of the enigmatic C. cropanii, develop a time-calibrated phylogeny, and explore the biogeographic history of the genus. We found that Corallus diversified within mainland South America, via over-water dispersals to the Lesser Antilles and Central America, and via the traditionally recognized Panamanian land bridge. Divergence time estimates reject the South American Caribbean-Track as a general biogeographic model for Corallus and implicate a role for events during the Oligocene and Miocene in diversification such as marine incursions and the uplift of the Andes. Our findings also suggest that recognition of the island endemic species, C. grenadensis and C. cookii, is questionable as they are nested within the widely distributed species, C. hortulanus. Our results highlight the importance of using widespread taxa when forming and testing biogeographic hypotheses in complex regions and further illustrate the difficulty of forming broadly applicable hypotheses regarding patterns of diversification in the Neotropical region.

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.

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.