Limitations of climatic data for inferring species boundaries: insights from speckled rattlesnakes

Phylogenomic discordance is driven by wide-spread introgression and incomplete lineage sorting during rapid species diversification within rattlesnakes (Viperidae: Crotalus and Sistrurus)

Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Venom comparisons of endemic and micro-endemic speckled rattlesnakes Crotalus mitchellii, C. polisi and C. thalassoporus from Baja California Peninsula

A high diversity of rattlesnake species can be found in the Baja California peninsula and the island of the Gulf of California, nevertheless, their venom has been poorly evaluated. The aim of this work was to present the first characterization of endemic Crotalus mitchellii, micro endemic C. polisi and C. thalassoporus venoms. All samples provoke human plasma coagulation showing doses in the rank of 2.3-41.0 μg and also produce rapid hydrolysis of the alpha chain of bovine fibrinogen while the beta chain is attacked at larger incubation periods by C. polisi and especially by C. thalassoporus. Phospholipase activity ranging from 23.2 to 173.8 U/mg. The venoms of C. thalassoporus and C. polisi show very high hemorrhagic activity (from 0.03 to 0.31 μg). A total of 130 toxin-related proteins were identified and classified into ten families. Crotalus mitchellii venom was characterized by high abundance of crotoxin-like and other phospholipase proteins (34.5%) and serine proteinases (29.8%). Crotalus polisi showed a similar proportion of metalloproteinases (34%) and serine proteinases (22.8%) components with important contribution of C-type lectins (14.3%) and CRiSP (14.0%) proteins. Venom of C. thalassoporus is dominated by metalloproteases that amount to more than 66% of total toxin proteins. These results provide a foundation for comprehending the biological, ecological and evolutionary significance of venom composition of speckled rattlesnake from the Baja California peninsula.

Assessing genomic and ecological differentiation among subspecies of the Rough-footed Mud Turtle, Kinosternon hirtipes

Combining genetic and ecological measures of differentiation can provide compelling evidence for ecological and genetic divergence among lineages. The Rough-footed Mud Turtle, Kinosternon hirtipes, is distributed from the Trans-Pecos region of Texas to the highlands of Central Mexico and contains six described subspecies, five of which are extant. We use ddRAD sequencing and species distribution models to assess levels of ecological and genetic differentiation among these subspecies. We also predict changes in climatically suitable habitat under different climate change scenarios and assess levels of genetic diversity and inbreeding within each lineage. Our results show that there is strong genetic and ecological differentiation among multiple lineages within K. hirtipes, and that this differentiation appears to be the result of vicariance associated with the Trans-Mexican Volcanic Belt. We propose changes to subspecies designations to more accurately reflect the evolutionary relationships among populations and assess threats to each subspecies.

VenomMaps: Updated species distribution maps and models for New World pitvipers (Viperidae: Crotalinae)

Beyond providing critical information to biologists, species distributions are useful for naturalists, curious citizens, and applied disciplines including conservation planning and medical intervention. Venomous snakes are one group that highlight the importance of having accurate information given their cosmopolitan distribution and medical significance. Envenomation by snakebite is considered a neglected tropical disease by the World Health Organization and venomous snake distributions are used to assess vulnerability to snakebite based on species occurrence and antivenom/healthcare accessibility. However, recent studies highlighted the need for updated fine-scale distributions of venomous snakes. Pitvipers (Viperidae: Crotalinae) are responsible for >98% of snakebites in the New World. Therefore, to begin to address the need for updated fine-scale distributions, we created VenomMaps, a database and web application containing updated distribution maps and species distribution models for all species of New World pitvipers. With these distributions, biologists can better understand the biogeography and conservation status of this group, researchers can better assess vulnerability to snakebite, and medical professionals can easily discern species found in their area.

Measurement(s)	Species Distributions
Technology Type(s)	Geographic Information System • Species Distribution Model (MaxEnt/kuenm)
Factor Type(s)	Occurrence Records • Environmental Data
Sample Characteristic - Organism	Crotalinae
Sample Characteristic - Location	North America • South America

Genome-wide data reveal extensive gene flow during the diversification of the western rattlesnakes (Viperidae: Crotalinae: Crotalus)

Hybridization and introgression are important, but often overlooked processes when inferring phylogenies. When these processes are not accounted for and a strictly diverging phylogenetic model is applied to groups with a history of hybridization, phylogenetic inference and parameter estimation can be inaccurate. Recent developments in phylogenetic network approaches coupled with the increasing availability of genomic data allow inferences of reticulate evolutionary histories across the tree of life. The western rattlesnake species group (C. viridis species complex, C. mitchellii species complex, C. scutulutas, and C. tigris) is an iconic snake lineage that is widespread across western North America. This group is composed of several species complexes with unclear species limits, likely the result of ongoing gene flow among nascent lineages. Here I generate reduced representation genomic data and test for a history of reticulation within this group. I demonstrate that all species have undergone hybridization with at least one other lineage, suggesting introgression is widespread in this group. Topologies differ between phylogenies estimated under the multispecies coalescent and multispecies network coalescent methods, indicating that gene flow has obscured phylogenetic relationships within this group. These past introgression events are predominantly restricted to species that co-occur geographically. However, within species that have a history of introgression, this signature is detected regardless of specimen sampling across geography. Overall, my results suggest the accumulation of reproductive isolating barriers occurs slowly in rattlesnakes which likely leads to the difficulty in delimiting species, furthermore, the results of this study have implications for trait evolution in this group.

Allopatric divergence and secondary contact with gene flow: a recurring theme in rattlesnake speciation

The study of recently diverged lineages whose geographical ranges come into contact can provide insight into the early stages of speciation and the potential roles of reproductive isolation in generating and maintaining species. Such insight can also be important for understanding the strategies and challenges for delimiting species within recently diverged species complexes. Here, we use mitochondrial and nuclear genetic data to study population structure, gene flow and demographic history across a geographically widespread rattlesnake clade, the western rattlesnake species complex (Crotalus cerberus, Crotalus viridis, Crotalus oreganus and relatives), which contains multiple lineages with ranges that overlap geographically or contact one another. We find evidence that the evolutionary history of this group does not conform to a bifurcating tree model and that pervasive gene flow has broadly influenced patterns of present-day genetic diversity. Our results suggest that lineage diversity has been shaped largely by drift and divergent selection in isolation, followed by secondary contact, in which reproductive isolating mechanisms appear weak and insufficient to prevent introgression, even between anciently diverged lineages. The complexity of divergence and secondary contact with gene flow among lineages also provides new context for why delimiting species within this complex has been difficult and contentious historically.

A conservation checklist of the amphibians and reptiles of Sonora, Mexico, with updated species lists

Sonora has a rich natural diversity, including reptiles and amphibians. Sonora's location on the United States-Mexico border creates some unique conservation challenges for its wildlife. We compiled a list of the amphibian and reptile species currently known for Sonora, summarized the conservation status of these species, and compared our list of species with known species lists for adjacent states. The herpetofauna of Sonora comprises 200 species of amphibians and reptiles (38 amphibians and 162 reptiles). Overall, Sonora shares the most species with Chihuahua, Sinaloa, and Arizona. Approximately 11% of the amphibian and reptile species are IUCN listed, but 35.5% are placed in a protected category by SEMARNAT, and 32.6% are categorized as high risk by the Environmental Vulnerability Score.

DNA barcodes from snake venom: a broadly applicable method for extraction of DNA from snake venoms

DNA barcoding is a simple technique used to develop a large-scale system of classification that is broadly applicable across a wide variety of taxa. DNA-based analysis of snake venoms can provide a system of classification independent of currently accepted taxonomic relationships by generating DNA barcodes specific to each venom sample. DNA purification from dried snake venoms has previously required large amounts of starting material, has resulted in low yields and inconsistent amplification, and was possible with front-fanged snakes only. Here, we present a modified DNA extraction protocol applied to venoms of both front- and rear-fanged snakes that requires significantly less starting material (1 mg) and yields sufficient amounts of DNA for successful PCR amplification of regions commonly used for DNA barcoding. [Formula: see text].

Inferring Passenger Denial Behavior of Taxi Drivers from Large-Scale Taxi Traces

How to understand individual human actions is a fundamental question to modern science, which drives and incurs many social, technological, racial, religious and economic phenomena. Human dynamics tries to reveal the temporal pattern and internal mechanism of human actions in letter or electronic communications, from the perspective of continuous interactions among friends or acquaintances. For interactions between stranger to stranger, taxi industry provide fruitful phenomina and evidence to investigate the action decisions. In fact, one striking disturbing events commonly reported in taxi industry is passenger refusing or denial, whose reasons vary, including skin color, blind passenger, being a foreigner or too close destination, religion reasons and anti specific nationality, so that complaints about taxi passenger refusing have to be concerned and processed carefully by local governments. But more universal factors for this phenomena are of great significance, which might be fulfilled by big data research to obtain novel insights in this question. In this paper, we demonstrate the big data analytics application in revealing novel insights from massive taxi trace data, which, for the first time, validates the passengers denial in taxi industry and estimates the denial ratio in Beijing city. We first quantify the income differentiation facts among taxi drivers. Then we find out that choosing the drop-off places also contributes to the high income for taxi drivers, compared to the previous explanation of mobility intelligence. Moreover, we propose the pick-up, drop-off and grid diversity concepts and related diversity analysis suggest that, high income taxi drivers will deny passengers in some situations, so as to choose the passengers' destination they prefer. Finally we design an estimation method for denial ratio and infer that high income taxi drivers will deny passengers with 8.52% likelihood in Beijing. Our work exhibits the power of big data analysis in revealing some dark side investigation.

Including Fossils in Phylogenetic Climate Reconstructions: A Deep Time Perspective on the Climatic Niche Evolution and Diversification of Spiny Lizards (Sceloporus)

Fossils and other paleontological information can improve phylogenetic comparative method estimates of phenotypic evolution and generate hypotheses related to species diversification. Here, we use fossil information to calibrate ancestral reconstructions of suitable climate for Sceloporus lizards in North America. Integrating data from the fossil record, general circulation models of paleoclimate during the Miocene, climate envelope modeling, and phylogenetic comparative methods provides a geographically and temporally explicit species distribution model of Sceloporus-suitable habitat through time. We provide evidence to support the historic biogeographic hypothesis of Sceloporus diversification in warm North American deserts and suggest a relatively recent Sceloporus invasion into Mexico around 6 Ma. We use a physiological model to map extinction risk. We suggest that the number of hours of restriction to a thermal refuge limited Sceloporus from inhabiting Mexico until the climate cooled enough to provide suitable habitat at approximately 6 Ma. If the future climate returns to the hotter climates of the past, Mexico, the place of highest modern Sceloporus richness, will no longer provide suitable habitats for Sceloporus to survive and reproduce.