Ecological drivers of ultraviolet colour evolution in snakes.

Ultraviolet (UV) colour patterns invisible to humans are widespread in nature. However, research bias favouring species with conspicuous colours under sexual selection can limit our assessment of other ecological drivers of UV colour, like interactions between predators and prey. Here we demonstrate widespread UV colouration across Western Hemisphere snakes and find stronger support for a predator defence function than for reproduction. We find that UV colouration has evolved repeatedly in species with ecologies most sensitive to bird predation, with no sexual dichromatism at any life stage. By modelling visual systems of potential predators, we find that snake conspicuousness correlates with UV colouration and predator cone number, providing a plausible mechanism for selection. Our results suggest that UV reflectance should not be assumed absent in "cryptically coloured" animals, as signalling beyond human visual capacities may be a key outcome of species interactions in many taxa for which UV colour is likely underreported.

Specimen collection is essential for modern science.

Natural history museums are vital repositories of specimens, samples and data that inform about the natural world; this Formal Comment revisits a Perspective that advocated for the adoption of compassionate collection practices, querying whether it will ever be possible to completely do away with whole animal specimen collection.

Evolution of Three-Finger Toxin Genes in Neotropical Colubrine Snakes (Colubridae).

Snake venom research has historically focused on front-fanged species (Viperidae and Elapidae), limiting our knowledge of venom evolution in rear-fanged snakes across their ecologically diverse phylogeny. Three-finger toxins (3FTxs) are a known neurotoxic component in the venoms of some rear-fanged snakes (Colubridae: Colubrinae), but it is unclear how prevalent 3FTxs are both in expression within venom glands and more broadly among colubrine species. Here, we used a transcriptomic approach to characterize the venom expression profiles of four species of colubrine snakes from the Neotropics that were dominated by 3FTx expression (in the genera Chironius, Oxybelis, Rhinobothryum, and Spilotes). By reconstructing the gene trees of 3FTxs, we found evidence of putative novel heterodimers in the sequences of Chironius multiventris and Oxybelis aeneus, revealing an instance of parallel evolution of this structural change in 3FTxs among rear-fanged colubrine snakes. We also found positive selection at sites within structural loops or "fingers" of 3FTxs, indicating these areas may be key binding sites that interact with prey target molecules. Overall, our results highlight the importance of exploring the venoms of understudied species in reconstructing the full evolutionary history of toxins across the tree of life.

Parsing variance by marker type: Testing biogeographic hypotheses and differential contribution of historical processes to population structure in a desert lizard.

A fundamental goal of population genetic studies is to identify historical biogeographic patterns and understand the processes that generate them. However, localized demographic events can skew population genetic inference. Assessing populations with multiple types of genetic markers, each with unique mutation rates and responses to changes in population size, can help to identify potentially confounding population-specific demographic processes. Here, we compared population structure and connectivity inferred from microsatellites and restriction site-associated DNA loci among 17 populations of an arid-specialist lizard, the desert night lizard, Xantusia vigilis, in central California to test among historical processes structuring population genetic diversity. We found that both marker types yielded generally concordant insights into population genetic structure including a major phylogenetic break maintained between two populations separated by less than 10 km, suggesting that either marker type could be used to understand generalized demographic patterns across the region for management purposes. However, we also found that the effects of demography on marker discordance could be used to elucidate population histories and distinguish among competing biogeographic hypotheses. Our results suggest that comparisons of within-population diversity across marker types provide powerful opportunities for leveraging marker discordance, particularly for understanding the creation and maintenance of contact zones among clades.

Leveraging natural history collections to understand the impacts of global change.

This joint Special Feature focuses on the contributions and potential of natural history collections to address global change questions.

Assembling microbial communities: a genomic analysis of a natural experiment in neotropical bamboo internodes.

Microbes participate in ecological communities, much like multicellular organisms. However, microbial communities lack the centuries of observation and theory describing and predicting ecological processes available for multicellular organisms. Here, we examine early bacterial community assembly in the water-filled internodes of Amazonian bamboos from the genus Guadua. Bamboo stands form distinct habitat patches within the lowland Amazonian rainforest and provide habitat for a suite of vertebrate and invertebrate species. Guadua bamboos develop sealed, water-filled internodes as they grow. Internodes are presumed sterile or near sterile while closed, but most are eventually opened to the environment by animals, after which they are colonized by microbes. We find that microbial community diversity increases sharply over the first few days of environmental exposure, and taxonomic identity of the microbes changes through this time period as is predicted for early community assembly in macroscopic communities. Microbial community taxonomic turnover is consistent at the bacteria phylum level, but at the level of Operational Taxonomic Units (OTUs), internode communities become increasingly differentiated through time. We argue that these tropical bamboos form an ideal study system for microbial community ecology due to their near-sterile condition prior to opening, relatively consistent environment after opening, and functionally limitless possibilities for replicates. Given the possible importance of opened internode habitats as locations of transmission for both pathogenic and beneficial microbes among animals, understanding the microbial dynamics of the internode habitat is a key conservation concern for the insect and amphibian species that use this microhabitat.

Divergent Specialization of Simple Venom Gene Profiles among Rear-Fanged Snake Genera (Helicops and Leptodeira, Dipsadinae, Colubridae).

Many venomous animals express toxins that show extraordinary levels of variation both within and among species. In snakes, most studies of venom variation focus on front-fanged species in the families Viperidae and Elapidae, even though rear-fanged snakes in other families vary along the same ecological axes important to venom evolution. Here we characterized venom gland transcriptomes from 19 snakes across two dipsadine rear-fanged genera (Leptodeira and Helicops, Colubridae) and two front-fanged genera (Bothrops, Viperidae; Micrurus, Elapidae). We compared patterns of composition, variation, and diversity in venom transcripts within and among all four genera. Venom gland transcriptomes of rear-fanged Helicops and Leptodeira and front-fanged Micrurus are each dominated by expression of single toxin families (C-type lectins, snake venom metalloproteinase, and phospholipase A2, respectively), unlike highly diverse front-fanged Bothrops venoms. In addition, expression patterns of congeners are much more similar to each other than they are to species from other genera. These results illustrate the repeatability of simple venom profiles in rear-fanged snakes and the potential for relatively constrained venom composition within genera.

The Batrachian Barf Bowl: An authentic research experience using ecological data from frog diets.

Authentic research experiences (AREs) are a powerful strategy for inspiring and retaining students in science, technology, engineering, and math (STEM) fields. However, recent demand for virtual learning has emphasized the need for remote AREs that also foster a sense of community and interpersonal connections among participants. Here, we describe an ARE activity that leverages digitized diet data from natural history collections to provide students with collaborative research experience across any learning environment. Using magnified photographs of frog stomach contents collected in the Peruvian Amazon, we designed an open-source "bowl game" competition that challenges students to identify, measure, and compare diet items across vouchered frog specimens ("Batrachian Barf Bowl"). To demonstrate learning outcomes, we ran this activity with 39 herpetology class students from the University of Notre Dame and the University of Michigan. We used pre- and post-activity assessments to evaluate effectiveness, scientific accuracy of results, and impact on student well-being. With minimal preparation and training in invertebrate identification, students were successful in identifying hundreds of frog diet items to taxonomic order, although accuracy varied among clades (global accuracy ~70%). While we found no difference in science identity, community, or self-efficacy between the two institutions at either time point (pre- and post-activity), we found that well-being was significantly higher for both sets of students after the activity. Overall, this approach offers a model for combining active learning with museum collections to provide experiential research opportunities that highlight the power of scientific collaboration.

Batch-Mask: Automated Image Segmentation for Organisms with Limbless or Non-Standard Body Forms.

Efficient comparisons of biological color patterns are critical for understanding the mechanisms by which organisms evolve in nature, including sexual selection, predator-prey interactions, and thermoregulation. However, limbless, elongate, or spiral-shaped organisms do not conform to the standard orientation and photographic techniques required for many automated analyses. Currently, large-scale color analysis of elongate animals requires time-consuming manual landmarking, which reduces their representation in coloration research despite their ecological importance. We present Batch-Mask: an automated, customizable workflow to automatically analyze large photographic datasets to isolate non-standard biological organisms from the background. Batch-Mask is completely open-source and does not depend on any proprietary software. We also present a user guide for fine-tuning weights to a custom dataset and incorporating existing manual visual analysis tools (e.g., micaToolbox) into a single automated workflow for comparing color patterns across images. Batch-Mask was 60x faster than manual landmarking and produced masks that correctly identified 96% of all snake pixels. To validate our approach, we used micaToolbox to compare pattern energy in a sample set of snake photographs segmented by Batch-Mask and humans and found no significant difference in the output results. The fine-tuned weights, user guide, and automated workflow substantially decrease the amount of time and attention required to quantitatively analyze non-standard biological subjects. With these tools, biologists can compare color, pattern, and shape differences in large datasets that include significant morphological variation in elongate body forms. This advance is especially valuable for comparative analyses of natural history collections across a broad range of morphologies. Through landmark-free automation, Batch-Mask can greatly expand the scale of space, time, or taxonomic breadth across which color variation can be quantitatively examined.

A guide for optimal iodine staining and high-throughput diceCT scanning in snakes.

Diffusible iodine-based contrast-enhanced computed tomography (diceCT) visualizes soft tissue from micro-CT (µCT) scans of specimens to uncover internal features and natural history information without incurring physical damage via dissection. Unlike hard-tissue imaging, taxonomic sampling within diceCT datasets is currently limited. To initiate best practices for diceCT in a nonmodel group, we outline a guide for staining and high-throughput µCT scanning in snakes. We scanned the entire body and one region of interest (i.e., head) for 23 specimens representing 23 species from the clades Aniliidae, Dipsadinae, Colubrinae, Elapidae, Lamprophiidae, and Viperidae. We generated 82 scans that include 1.25% Lugol's iodine stained (soft tissue) and unstained (skeletal) data for each specimen. We found that duration of optimal staining time increased linearly with body size; head radius was the best indicator. Postreconstruction of scans, optimal staining was evident by evenly distributed grayscale values and clear differentiation among soft-tissue anatomy. Under and over stained specimens produced poor contrast among soft tissues, which was often exacerbated by user bias during "digital dissections" (i.e., segmentation). Regardless, all scans produced usable data from which we assessed a range of downstream analytical applications within ecology and evolution (e.g., predator-prey interactions, life history, and morphological evolution). Ethanol destaining reversed the known effects of iodine on the exterior appearance of physical specimens, but required substantially more time than reported for other destaining methods. We discuss the feasibility of implementing diceCT techniques for a new user, including approximate financial and temporal commitments, required facilities, and potential effects of staining on specimens. We present the first high-throughput workflow for full-body skeletal and diceCT scanning in snakes, which can be generalized to any elongate vertebrates, and increases publicly available diceCT scans for reptiles by an order of magnitude.

Genetic mechanisms and correlational selection structure trait variation in a coral snake mimic.

Covariation among traits shapes both phenotypic evolution and ecological interactions across space and time. However, rampant geographical variation in the strength and direction of such correlations can be particularly difficult to explain through generalized mechanisms. By integrating population genomics, surveys of natural history collections and spatially explicit analyses, we tested multiple drivers of trait correlations in a coral snake mimic that exhibits remarkable polymorphism in mimetic and non-mimetic colour traits. We found that although such traits co-occur extensively across space, correlations were best explained by a mixture of genetic architecture and correlational selection, rather than by any single mechanism. Our findings suggest that spatially complex trait distributions may be driven more by the simple interaction between multiple processes than by complex variation in one mechanism alone. These interactions are particularly important in mimicry systems, which frequently generate striking geographical variation and genetic correlations among colour pattern traits.

What makes a fang? Phylogenetic and ecological controls on tooth evolution in rear-fanged snakes.

BACKGROUND: Fangs are a putative key innovation that revolutionized prey capture and feeding in snakes, and - along with their associated venom phenotypes - have made snakes perhaps the most medically-significant vertebrate animals. Three snake clades are known for their forward-positioned fangs, and these clades (Elapidae, Viperidae, and Atractaspidinae) contain the majority of snakes that are traditionally considered venomous. However, many other snakes are "rear-fanged": they possess potentially venom-delivering teeth situated at the rear end of the upper jaw. Quantification of fang phenotypes - and especially those of rear-fanged species - has proved challenging or impossible owing to the small size and relative rarity of many such snakes. Consequently, it has been difficult to understand the evolutionary history of both venom and prey-capture strategies across extant snakes. We quantified variation in the dentition of 145 colubriform ("advanced") snake species using microCT scanning and compared dental characters with ecological data on species' diet and prey capture method(s) to understand broader patterns in snake fang evolution. RESULTS: Dental traits such as maxilla length, tooth number, and fang size show strong phylogenetic signal across Colubriformes. We find extreme heterogeneity and evolutionary lability in the rear-fanged phenotype in colubrid (colubrine, dipsadine, and natricine lineages) and lamprophiid snakes, in contrast to relative uniformity in the front fanged phenotypes of other groups (vipers and, to a lesser extent, elapids). Fang size and position are correlated with venom-use in vipers, elapids, and colubrid snakes, with the latter group shifting fangs anteriorly by shortening the entire maxillary bone. We find that maxilla length and tooth number may also be correlated with the evolution of dietary specialization. Finally, an ancestral state reconstruction suggests that fang loss is a recurring phenomenon in colubrid snakes, likely accompanied by shifts in diet and prey capture mode. CONCLUSIONS: Our study provides a framework for quantifying the complex morphologies associated with venom use in snakes. Our results suggest that fang phenotypes, and particularly the rear-fanged phenotype, in snakes are both diverse and labile, facilitating a wide range of ecological strategies and contributing to spectacular radiations of these organisms in tropical and subtropical biomes worldwide.

Metabolically similar cohorts of bacteria exhibit strong cooccurrence patterns with diet items and eukaryotic microbes in lizard guts.

Gut microbiomes perform essential services for their hosts, including helping them to digest food and manage pathogens and parasites. Performing these services requires a diverse and constantly changing set of metabolic functions from the bacteria in the microbiome. The metabolic repertoire of the microbiome is ultimately dependent on the outcomes of the ecological interactions of its member microbes, as these interactions in part determine the taxonomic composition of the microbiome. The ecological processes that underpin the microbiome's ability to handle a variety of metabolic challenges might involve rapid turnover of the gut microbiome in response to new metabolic challenges, or it might entail maintaining sufficient diversity in the microbiome that any new metabolic demands can be met from an existing set of bacteria. To differentiate between these scenarios, we examine the gut bacteria and resident eukaryotes of two generalist-insectivore lizards, while simultaneously identifying the arthropod prey each lizard was digesting at the time of sampling. We find that the cohorts of bacteria that occur significantly more or less often than expected with arthropod diet items or eukaryotes include bacterial species that are highly similar to each other metabolically. This pattern in the bacterial microbiome could represent an early step in the taxonomic shifts in bacterial microbiome that occur when host lineages change their diet niche over evolutionary timescales.

Widespread chytrid infection across frogs in the Peruvian Amazon suggests critical role for low elevation in pathogen spread and persistence.

Outbreaks of emerging infectious diseases are becoming more frequent as climate changes wildlife communities at unprecedented rates, driving population declines and raising concerns for species conservation. One critical disease is the global pandemic of chytridiomycosis in frogs, which can be caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Although there is clear evidence for Bd-induced mortality across high-elevation frog communities, little attention is given to the role of lowlands in Bd's persistence and spread because low elevations are assumed to be too warm to harbor significant levels of Bd. Here, we report widespread Bd infection across 80 frog species from three sites in the lowland Peruvian Amazon, an area with no documented Bd-related amphibian declines. Despite observing no clinical signs of infection in the field, we found that 24-46% of individuals were infected per site (up to ≈105,000 zoospore equivalents per frog) by three Bd strains from the global pandemic lineage (Bd-GPL). We also found collection site and seasonal effects to be only weak predictors of Bd prevalence and load, with lower elevation and drier habitats marginally decreasing both prevalence and load. We found no further effect of host phylogeny, ecotype, or body size. Our results showing high and widespread prevalence across a lowland tropical ecosystem contradict the expectations based on the global pattern of pathogenicity of Bd that is largely restricted to higher elevations and colder temperatures. These findings imply that the lowlands may play a critical role in the spread and persistence of Bd over time and space.

Environmental heterogeneity and not vicariant biogeographic barriers generate community-wide population structure in desert-adapted snakes.

Genetic structure can be influenced by local adaptation to environmental heterogeneity and biogeographic barriers, resulting in discrete population clusters. Geographic distance among populations, however, can result in continuous clines of genetic divergence that appear as structured populations. Here, we evaluate the relevant importance of these three factors over a landscape characterized by environmental heterogeneity and the presence of a hypothesized biogeographic barrier in producing population genetic structure within 13 codistributed snake species using a genomic data set. We demonstrate that geographic distance and environmental heterogeneity across western North America contribute to population genomic divergence. Surprisingly, landscape features long thought to contribute to biogeographic barriers play little role in divergence community wide. Our results suggest that isolation by environment is the most important contributor to genomic divergence. Furthermore, we show that models of population clustering that incorporate spatial information consistently outperform nonspatial models, demonstrating the importance of considering geographic distances in population clustering. We argue that environmental and geographic distances as drivers of community-wide divergence should be explored before assuming the role of biogeographic barriers.

Natural history bycatch: a pipeline for identifying metagenomic sequences in RADseq data.

BACKGROUND: Reduced representation genomic datasets are increasingly becoming available from a variety of organisms. These datasets do not target specific genes, and so may contain sequences from parasites and other organisms present in the target tissue sample. In this paper, we demonstrate that (1) RADseq datasets can be used for exploratory analysis of tissue-specific metagenomes, and (2) tissue collections house complete metagenomic communities, which can be investigated and quantified by a variety of techniques. METHODS: We present an exploratory method for mining metagenomic "bycatch" sequences from a range of host tissue types. We use a combination of the pyRAD assembly pipeline, NCBI's blastn software, and custom R scripts to isolate metagenomic sequences from RADseq type datasets. RESULTS: When we focus on sequences that align with existing references in NCBI's GenBank, we find that between three and five percent of identifiable double-digest restriction site associated DNA (ddRAD) sequences from host tissue samples are from phyla to contain known blood parasites. In addition to tissue samples, we examine ddRAD sequences from metagenomic DNA extracted snake and lizard hind-gut samples. We find that the sequences recovered from these samples match with expected bacterial and eukaryotic gut microbiome phyla. DISCUSSION: Our results suggest that (1) museum tissue banks originally collected for host DNA archiving are also preserving valuable parasite and microbiome communities, (2) that publicly available RADseq datasets may include metagenomic sequences that could be explored, and (3) that restriction site approaches are a useful exploratory technique to identify microbiome lineages that could be missed by primer-based approaches.

Predator Perspective Drives Geographic Variation in Frequency-Dependent Polymorphism.

Color polymorphism in natural populations can manifest as a striking patchwork of phenotypes in space, with neighboring populations characterized by dramatic differences in morph composition. These geographic mosaics can be challenging to explain in the absence of localized selection because they are unlikely to result from simple isolation-by-distance or clinal variation in selective regimes. To identify processes that can lead to the formation of geographic mosaics, we developed a simulation-based model to explore the influence of predator perspective, selection, migration, and genetic linkage of color loci on allele frequencies in polymorphic populations over space and time. Using simulated populations inspired by the biology of Heliconius longwing butterflies, Cepaea land snails, Oophaga poison frogs, and Sonora ground snakes, we found that the relative sizes of predator and prey home ranges can produce large differences in morph composition between neighboring populations under both positive and negative frequency-dependent selection. We also demonstrated the importance of the interaction of predator perspective with the type of frequency dependence and localized directional selection across migration and selection intensities. Our results show that regional-scale predation can promote the formation of phenotypic mosaics in prey species, without the need to invoke spatial variation in selective regimes. We suggest that predator behavior can play an important and underappreciated role in the formation and maintenance of geographic mosaics in polymorphic species.

Anolis lizards as biocontrol agents in mainland and island agroecosystems.

Our knowledge of ecological interactions that bolster ecosystem function and productivity has broad applications to the management of agricultural systems. Studies suggest that the presence of generalist predators in agricultural landscapes leads to a decrease in the abundance of herbivorous pests, but our understanding of how these interactions vary across taxa and along gradients of management intensity and eco-geographic space remains incomplete. In this study, we assessed the functional response and biocontrol potential of a highly ubiquitous insectivore (lizards in the genus Anolis) on the world's most important coffee pest, the coffee berry borer (Hypothalemus hampei). We conducted field surveys and laboratory experiments to examine the impact of land-use intensification on species richness and abundance of anoles and the capacity of anoles to reduce berry borer infestations in mainland and island coffee systems. Our results show that anoles significantly reduce coffee infestation rates in laboratory settings (Mexico, p = .03, F = 5.13 df = 1, 35; Puerto Rico, p = .014, F = 8.82, df = 1, 10) and are capable of consuming coffee berry borers in high abundance. Additionally, diversified agroecosystems bolster anole abundance, while high-intensity practices, including the reduction of vegetation complexity and the application of agrochemicals were associated with reduced anole abundance. The results of this study provide supporting evidence of the positive impact of generalist predators on the control of crop pests in agricultural landscapes, and the role of diversified agroecosystems in sustaining both functionally diverse communities and crop production in tropical agroecosystems.

Historical Environment Is Reflected in Modern Population Genetics and Biogeography of an Island Endemic Lizard (Xantusia riversiana reticulata).

The restricted distribution and isolation of island endemics often produces unique genetic and phenotypic diversity of conservation interest to management agencies. However, these isolated species, especially those with sensitive life history traits, are at high risk for the adverse effects of genetic drift and habitat degradation by non-native wildlife. Here, we study the population genetic diversity, structure, and stability of a classic "island giant" (Xantusia riversiana, the Island Night Lizard) on San Clemente Island, California following the removal of feral goats. Using DNA microsatellites, we found that this population is reasonably genetically robust despite historical grazing, with similar effective population sizes and genetic diversity metrics across all sampling locations irrespective of habitat type and degree of degradation. However, we also found strong site-specific patterns of genetic variation and low genetic diversity compared to mainland congeners, warranting continued special management as an island endemic. We identify both high and low elevation areas that remain valuable repositories of genetic diversity and provide a case study for other low-dispersal coastal organisms in the face of future climate change.

Coral snakes predict the evolution of mimicry across New World snakes.

Batesian mimicry, in which harmless species (mimics) deter predators by deceitfully imitating the warning signals of noxious species (models), generates striking cases of phenotypic convergence that are classic examples of evolution by natural selection. However, mimicry of venomous coral snakes has remained controversial because of unresolved conflict between the predictions of mimicry theory and empirical patterns in the distribution and abundance of snakes. Here we integrate distributional, phenotypic and phylogenetic data across all New World snake species to demonstrate that shifts to mimetic coloration in nonvenomous snakes are highly correlated with coral snakes in both space and time, providing overwhelming support for Batesian mimicry. We also find that bidirectional transitions between mimetic and cryptic coloration are unexpectedly frequent over both long- and short-time scales, challenging traditional views of mimicry as a stable evolutionary 'end point' and suggesting that insect and snake mimicry may have different evolutionary dynamics.