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.

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.

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.

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.

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.

Convergence and divergence in anti-predator displays: a novel approach to quantitative behavioural comparison in snakes

Animals in nature use many strategies to evade or deter their predators. Within venomous snake mimicry, stereotyped anti-predator behaviours are hypothesized to be effective warning signals under strong selection for independent convergence across species. However, no studies have systematically quantified snake anti-predator displays across taxonomically broad clades to examine how these behaviours evolve within a comparative methods framework. Here we describe a new high-throughput approach for collecting and quantifying anti-predator displays in snakes that demonstrates both low observer bias and infinite extension. Then, we show this method’s utility by comparing 20 species spanning six taxonomic families from Peru. We found that a simple experimental set-up varying simulated predator cues was successful in eliciting displays across species and that high-speed videography captured a great diversity of anti-predator responses. Although display components show complicated patterns of covariance, we found support for behavioural convergence in anti-predator displays among elapid coral snakes and their distantly related mimics. Our approach provides new analytical opportunities for both behaviour and kinematics, especially macroevolutionary analyses across clades with similar difficulty in observing or comparing trait diversity.

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.

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.

Unlinked Mendelian inheritance of red and black pigmentation in snakes: Implications for Batesian mimicry

Identifying the genetic basis of mimetic signals is critical to understanding both the origin and dynamics of mimicry over time. For species not amenable to large laboratory breeding studies, widespread color polymorphism across natural populations offers a powerful way to assess the relative likelihood of different genetic systems given observed phenotypic frequencies. We classified color phenotype for 2175 ground snakes (Sonora semiannulata) across the continental United States to analyze morph ratios and test among competing hypotheses about the genetic architecture underlying red and black coloration in coral snake mimics. We found strong support for a two-locus model under simple Mendelian inheritance, with red and black pigmentation being controlled by separate loci. We found no evidence of either linkage disequilibrium between loci or sex linkage. In contrast to Batesian mimicry systems such as butterflies in which all color signal components are linked into a single "supergene," our results suggest that the mimetic signal in colubrid snakes can be disrupted through simple recombination and that color evolution is likely to involve discrete gains and losses of each signal component. Both outcomes are likely to contribute to the exponential increase in rates of color evolution seen in snake mimicry systems over insect systems.

Report from the first snake genomics and integrative biology meeting

This report summarizes the proceedings of the 1st Snake Genomics and Integrative Biology Meeting held in Vail, CO USA, 5-8 October 2011. The meeting had over twenty registered participants, and was conducted as a single session of presentations. Goals of the meeting included coordination of genomic data collection and fostering collaborative interactions among researchers using snakes as model systems.

Direct fitness correlates and thermal consequences of facultative aggregation in a desert lizard

Social aggregation is a common behavioral phenomenon thought to evolve through adaptive benefits to group living. Comparing fitness differences between aggregated and solitary individuals in nature--necessary to infer an evolutionary benefit to living in groups--has proven difficult because communally-living species tend to be obligately social and behaviorally complex. However, these differences and the mechanisms driving them are critical to understanding how solitary individuals transition to group living, as well as how and why nascent social systems change over time. Here we demonstrate that facultative aggregation in a reptile (the Desert Night Lizard, Xantusia vigilis) confers direct reproductive success and survival advantages and that thermal benefits of winter huddling disproportionately benefit small juveniles, which can favor delayed dispersal of offspring and the formation of kin groups. Using climate projection models, however, we estimate that future aggregation in night lizards could decline more than 50% due to warmer temperatures. Our results support the theory that transitions to group living arise from direct benefits to social individuals and offer a clear mechanism for the origin of kin groups through juvenile philopatry. The temperature dependence of aggregation in this and other taxa suggests that environmental variation may be a powerful but underappreciated force in the rapid transition between social and solitary behavior.