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.

Jaw size variation is associated with a novel craniofacial function for galanin receptor 2 in an adaptive radiation of pupfishes.

Understanding the genetic basis of novel adaptations in new species is a fundamental question in biology. Here we demonstrate a new role for galr2 in vertebrate craniofacial development using an adaptive radiation of trophic specialist pupfishes endemic to San Salvador Island, Bahamas. We confirmed the loss of a putative Sry transcription factor binding site upstream of galr2 in scale-eating pupfish and found significant spatial differences in galr2 expression among pupfish species in Meckel's cartilage using in situ hybridization chain reaction (HCR). We then experimentally demonstrated a novel role for Galr2 in craniofacial development by exposing embryos to Garl2-inhibiting drugs. Galr2-inhibition reduced Meckel's cartilage length and increased chondrocyte density in both trophic specialists but not in the generalist genetic background. We propose a mechanism for jaw elongation in scale-eaters based on the reduced expression of galr2 due to the loss of a putative Sry binding site. Fewer Galr2 receptors in the scale-eater Meckel's cartilage may result in their enlarged jaw lengths as adults by limiting opportunities for a circulating Galr2 agonist to bind to these receptors during development. Our findings illustrate the growing utility of linking candidate adaptive SNPs in non-model systems with highly divergent phenotypes to novel vertebrate gene functions.

Jaw size variation is associated with a novel craniofacial function for galanin receptor 2 in an adaptive radiation of pupfishes.

Understanding the genetic basis of novel adaptations in new species is a fundamental question in biology that also provides an opportunity to uncover new genes and regulatory networks with potential clinical relevance. Here we demonstrate a new role for galr2 in vertebrate craniofacial development using an adaptive radiation of trophic specialist pupfishes endemic to San Salvador Island in the Bahamas. We confirmed the loss of a putative Sry transcription factor binding site in the upstream region of galr2 in scale-eating pupfish and found significant spatial differences in galr2 expression among pupfish species in Meckel's cartilage and premaxilla using in situ hybridization chain reaction (HCR). We then experimentally demonstrated a novel function for Galr2 in craniofacial development and jaw elongation by exposing embryos to drugs that inhibit Galr2 activity. Galr2-inhibition reduced Meckel's cartilage length and increased chondrocyte density in both trophic specialists but not in the generalist genetic background. We propose a mechanism for jaw elongation in scale-eaters based on the reduced expression of galr2 due to the loss of a putative Sry binding site. Fewer Galr2 receptors in the scale-eater Meckel's cartilage may result in their enlarged jaw lengths as adults by limiting opportunities for a postulated Galr2 agonist to bind to these receptors during development. Our findings illustrate the growing utility of linking candidate adaptive SNPs in non-model systems with highly divergent phenotypes to novel vertebrate gene functions.

Multiple performance peaks for scale-biting in an adaptive radiation of pupfishes.

The physical interactions between organisms and their environment ultimately shape their rate of speciation and adaptive radiation, but the contributions of biomechanics to evolutionary divergence are frequently overlooked. Here we investigated an adaptive radiation of Cyprinodon pupfishes to measure the relationship between feeding kinematics and performance during adaptation to a novel trophic niche, lepidophagy, in which a predator removes only the scales, mucus, and sometimes tissue from their prey using scraping and biting attacks. We used high-speed video to film scale-biting strikes on gelatin cubes by scale-eater, molluscivore, generalist, and hybrid pupfishes and subsequently measured the dimensions of each bite. We then trained the SLEAP machine-learning animal tracking model to measure kinematic landmarks and automatically scored over 100,000 frames from 227 recorded strikes. Scale-eaters exhibited increased peak gape and greater bite length; however, substantial within-individual kinematic variation resulted in poor discrimination of strikes by species or strike type. Nonetheless, a complex performance landscape with two distinct peaks best predicted gel-biting performance, corresponding to a significant nonlinear interaction between peak gape and peak jaw protrusion in which scale-eaters and their hybrids occupied a second performance peak requiring larger peak gape and greater jaw protrusion. A bite performance valley separating scale-eaters from other species may have contributed to their rapid evolution and is consistent with multiple estimates of a multi-peak fitness landscape in the wild. We thus present an efficient deep-learning automated pipeline for kinematic analyses of feeding strikes and a new biomechanical model for understanding the performance and rapid evolution of a rare trophic niche.

Severe inbreeding, increased mutation load and gene loss-of-function in the critically endangered Devils Hole pupfish.

Small populations with limited range are often threatened by inbreeding and reduced genetic diversity, which can reduce fitness and exacerbate population decline. One of the most extreme natural examples is the Devils Hole pupfish (Cyprinodon diabolis), an iconic and critically endangered species with the smallest known range of any vertebrate. This species has experienced severe declines in population size over the last 30 years and suffered major bottlenecks in 2007 and 2013, when the population shrunk to 38 and 35 individuals, respectively. Here, we analysed 30 resequenced genomes of desert pupfishes from Death Valley, Ash Meadows and surrounding areas to examine the genomic consequences of small population size. We found extremely high levels of inbreeding (FROH = 0.34-0.81) and an increased amount of potentially deleterious genetic variation in the Devils Hole pupfish as compared to other species, including unique, fixed loss-of-function alleles and deletions in genes associated with sperm motility and hypoxia. Additionally, we successfully resequenced a formalin-fixed museum specimen from 1980 and found that the population was already highly inbred prior to recent known bottlenecks. We thus document severe inbreeding and increased mutation load in the Devils Hole pupfish and identify candidate deleterious variants to inform management of this conservation icon.

Trophic specialization on unique resources despite limited niche divergence in a celebrated example of sympatric speciation.

Trophic niche partitioning is observed in many adaptive radiations and is hypothesized to be a central process underlying species divergence. However, patterns of dietary niche partitioning are inconsistent across radiations and there are few studies of niche partitioning in putative examples of sympatric speciation. Here, we conducted the first quantitative study of dietary niche partitioning using stomach contents and stable isotope analyses in one of the most celebrated examples of sympatric speciation: the cichlid radiation from crater lake Barombi Mbo, Cameroon. We found little evidence for trophic niche partitioning among cichlids, including the nine species coexisting in the narrow littoral zone. Stable isotope analyses supported these conclusions of substantial dietary overlap. Our data, however, did reveal that five of eleven species consume rare dietary items, including freshwater sponge, terrestrial ants, and nocturnal foraging on shrimp. Stomach contents of the spongivore (Pungu maclareni) were 20% freshwater sponge, notable considering that only 0.04% of all fishes consume sponges. Overall, we conclude that cichlid species in lake Barombi Mbo overlap considerably in broad dietary niches-in part due to the large proportion of detritus in the stomach contents of all species-but there is evidence for divergence among species in their diet specializations on unique resources. We speculate that these species may utilize these additional specialized resources during periods of low resource abundance in support of Liem's paradox.

Minimal overall divergence of the gut microbiome in an adaptive radiation of Cyprinodon pupfishes despite potential adaptive enrichment for scale-eating.

Adaptive radiations offer an excellent opportunity to understand the eco-evolutionary dynamics of gut microbiota and host niche specialization. In a laboratory common garden, we compared the gut microbiota of two novel derived trophic specialist pupfishes, a scale-eater and a molluscivore, to closely related and distant outgroup generalist populations, spanning both rapid trophic evolution within 10 kya and stable generalist diets persisting over 11 Mya. We predicted an adaptive and highly divergent microbiome composition in the trophic specialists reflecting their rapid rates of craniofacial and behavioral diversification. We sequenced 16S rRNA amplicons of gut microbiomes from lab-reared adult pupfishes raised under identical conditions and fed the same high protein diet. In contrast to our predictions, gut microbiota largely reflected phylogenetic distance among species, rather than generalist or specialist life history, in support of phylosymbiosis. However, we did find significant enrichment of Burkholderiaceae bacteria in replicated lab-reared scale-eater populations. These bacteria sometimes digest collagen, the major component of fish scales, supporting an adaptive shift. We also found some enrichment of Rhodobacteraceae and Planctomycetia in lab-reared molluscivore populations, but these bacteria target cellulose. Overall phylogenetic conservation of microbiome composition contrasts with predictions of adaptive radiation theory and observations of rapid diversification in all other trophic traits in these hosts, including craniofacial morphology, foraging behavior, aggression, and gene expression, suggesting that the functional role of these minor shifts in microbiota will be important for understanding the role of the microbiome in trophic diversification.

A new theoretical performance landscape for suction feeding reveals adaptive kinematics in a natural population of reef damselfish.

Understanding how organismal traits determine performance and, ultimately, fitness is a fundamental goal of evolutionary eco-morphology. However, multiple traits can interact in non-linear and context-dependent ways to affect performance, hindering efforts to place natural populations with respect to performance peaks or valleys. Here, we used an established mechanistic model of suction-feeding performance (SIFF) derived from hydrodynamic principles to estimate a theoretical performance landscape for zooplankton prey capture. This performance space can be used to predict prey capture performance for any combination of six morphological and kinematic trait values. We then mapped in situ high-speed video observations of suction feeding in a natural population of a coral reef zooplanktivore, Chromis viridis, onto the performance space to estimate the population's location with respect to the topography of the performance landscape. Although the kinematics of the natural population closely matched regions of high performance in the landscape, the population was not located on a performance peak. Individuals were furthest from performance peaks on the peak gape, ram speed and mouth opening speed trait axes. Moreover, we found that the trait combinations in the observed population were associated with higher performance than expected by chance, suggesting that these combinations are under selection. Our results provide a framework for assessing whether natural populations occupy performance optima.

We get by with a little help from our friends: shared adaptive variation provides a bridge to novel ecological specialists during adaptive radiation.

Adaptive radiations involve astounding bursts of phenotypic, ecological and species diversity. However, the microevolutionary processes that underlie the origins of these bursts are still poorly understood. We report the discovery of an intermediate C. sp. 'wide-mouth' scale-eating ecomorph in a sympatric radiation of Cyprinodon pupfishes, illuminating the transition from a widespread algae-eating generalist to a novel microendemic scale-eating specialist. We first show that this ecomorph occurs in sympatry with generalist C. variegatus and scale-eating specialist C. desquamator on San Salvador Island, Bahamas, but is genetically differentiated, morphologically distinct and often consumes scales. We then compared the timing of selective sweeps on shared and unique adaptive variants in trophic specialists to characterize their adaptive walk. Shared adaptive regions swept first in both the specialist desquamator and the intermediate 'wide-mouth' ecomorph, followed by unique sweeps of introgressed variation in 'wide-mouth' and de novo variation in desquamator. The two scale-eating populations additionally shared 9% of their hard selective sweeps with the molluscivore C. brontotheroides, despite no single common ancestor among specialists. Our work provides a new microevolutionary framework for investigating how major ecological transitions occur and illustrates how both shared and unique genetic variation can provide a bridge for multiple species to access novel ecological niches.

Hybridization alters the shape of the genotypic fitness landscape, increasing access to novel fitness peaks during adaptive radiation.

Estimating the complex relationship between fitness and genotype or phenotype (i.e. the adaptive landscape) is one of the central goals of evolutionary biology. However, adaptive walks connecting genotypes to organismal fitness, speciation, and novel ecological niches are still poorly understood and processes for surmounting fitness valleys remain controversial. One outstanding system for addressing these connections is a recent adaptive radiation of ecologically and morphologically novel pupfishes (a generalist, molluscivore, and scale-eater) endemic to San Salvador Island, Bahamas. We leveraged whole-genome sequencing of 139 hybrids from two independent field fitness experiments to identify the genomic basis of fitness, estimate genotypic fitness networks, and measure the accessibility of adaptive walks on the fitness landscape. We identified 132 single nucleotide polymorphisms (SNPs) that were significantly associated with fitness in field enclosures. Six out of the 13 regions most strongly associated with fitness contained differentially expressed genes and fixed SNPs between trophic specialists; one gene (mettl21e) was also misexpressed in lab-reared hybrids, suggesting a potential intrinsic genetic incompatibility. We then constructed genotypic fitness networks from adaptive alleles and show that scale-eating specialists are the most isolated of the three species on these networks. Intriguingly, introgressed and de novo variants reduced fitness landscape ruggedness as compared to standing variation, increasing the accessibility of genotypic fitness paths from generalist to specialists. Our results suggest that adaptive introgression and de novo mutations alter the shape of the fitness landscape, providing key connections in adaptive walks circumventing fitness valleys and triggering the evolution of novelty during adaptive radiation.

One of the main drivers of evolution is natural selection, which is when organisms better adapted to their environment are more likely to survive and reproduce. A common metaphor to explain this process is a landscape covered in peaks and valleys: the peaks represent genetic combinations or traits with high evolutionary fitness, while the valleys represent those with low fitness. As a population evolves and its environment changes, it moves among these peaks taking small steps across the landscape. However, there is a limit to how far an organism can travel in one leap. So, what happens when they need to cross a valley of low fitness to get to the next peak? To address this question, Patton et al. studied three young species of pupfish that recently evolved from a common ancestor and co-habit the same environment in the Caribbean. Patton et al. sequenced whole genomes of each new species and used this to build a genotypic fitness landscape, a network linking neighboring genotypes which each have a unique fitness value that was measured during field experiments. This revealed that most of the paths connecting the different species passed through valleys of low fitness. But there were rare, narrow ridges connecting each species. Next, Patton et al. found that new mutations as well as genetic variations that arose from mating with pupfish on other Caribbean islands altered genetic interactions and changed the shape of the fitness landscape. Ultimately, this significantly increased the accessibility of fitness peaks by both adding more ridges and decreasing the lengths of paths, expanding the realm of possible evolutionary outcomes. Understanding how fitness landscapes change during evolution could help to explain where new species come from. Other researchers could apply the same approach to estimate the genotypic fitness landscapes of other species, from bacteria to vertebrates. These networks could be used to visualize the complex fitness landscape that connects all lifeforms on Earth.

Linking ecomechanical models and functional traits to understand phenotypic diversity.

Physical principles and laws determine the set of possible organismal phenotypes. Constraints arising from development, the environment, and evolutionary history then yield workable, integrated phenotypes. We propose a theoretical and practical framework that considers the role of changing environments. This 'ecomechanical approach' integrates functional organismal traits with the ecological variables. This approach informs our ability to predict species shifts in survival and distribution and provides critical insights into phenotypic diversity. We outline how to use the ecomechanical paradigm using drag-induced bending in trees as an example. Our approach can be incorporated into existing research and help build interdisciplinary bridges. Finally, we identify key factors needed for mass data collection, analysis, and the dissemination of models relevant to this framework.

A vertebrate adaptive radiation is assembled from an ancient and disjunct spatiotemporal landscape.

To investigate the origins and stages of vertebrate adaptive radiation, we reconstructed the spatial and temporal histories of adaptive alleles underlying major phenotypic axes of diversification from the genomes of 202 Caribbean pupfishes. On a single Bahamian island, ancient standing variation from disjunct geographic sources was reassembled into new combinations under strong directional selection for adaptation to the novel trophic niches of scale-eating and molluscivory. We found evidence for two longstanding hypotheses of adaptive radiation: hybrid swarm origins and temporal stages of adaptation. Using a combination of population genomics, transcriptomics, and genome-wide association mapping, we demonstrate that this microendemic adaptive radiation of novel trophic specialists on San Salvador Island, Bahamas experienced twice as much adaptive introgression as generalist populations on neighboring islands and that adaptive divergence occurred in stages. First, standing regulatory variation in genes associated with feeding behavior (prlh, cfap20, and rmi1) were swept to fixation by selection, then standing regulatory variation in genes associated with craniofacial and muscular development (itga5, ext1, cyp26b1, and galr2) and finally the only de novo nonsynonymous substitution in an osteogenic transcription factor and oncogene (twist1) swept to fixation most recently. Our results demonstrate how ancient alleles maintained in distinct environmental refugia can be assembled into new adaptive combinations and provide a framework for reconstructing the spatiotemporal landscape of adaptation and speciation.

Few Fixed Variants between Trophic Specialist Pupfish Species Reveal Candidate Cis-Regulatory Alleles Underlying Rapid Craniofacial Divergence.

Investigating closely related species that rapidly evolved divergent feeding morphology is a powerful approach to identify genetic variation underlying variation in complex traits. This can also lead to the discovery of novel candidate genes influencing natural and clinical variation in human craniofacial phenotypes. We combined whole-genome resequencing of 258 individuals with 50 transcriptomes to identify candidate cis-acting genetic variation underlying rapidly evolving craniofacial phenotypes within an adaptive radiation of Cyprinodon pupfishes. This radiation consists of a dietary generalist species and two derived trophic niche specialists-a molluscivore and a scale-eating species. Despite extensive morphological divergence, these species only diverged 10 kya and produce fertile hybrids in the laboratory. Out of 9.3 million genome-wide SNPs and 80,012 structural variants, we found very few alleles fixed between species-only 157 SNPs and 87 deletions. Comparing gene expression across 38 purebred F1 offspring sampled at three early developmental stages, we identified 17 fixed variants within 10 kb of 12 genes that were highly differentially expressed between species. By measuring allele-specific expression in F1 hybrids from multiple crosses, we found that the majority of expression divergence between species was explained by trans-regulatory mechanisms. We also found strong evidence for two cis-regulatory alleles affecting expression divergence of two genes with putative effects on skeletal development (dync2li1 and pycr3). These results suggest that SNPs and structural variants contribute to the evolution of novel traits and highlight the utility of the San Salvador Island pupfish system as an evolutionary model for craniofacial development.

Surprising spatiotemporal stability of a multi-peak fitness landscape revealed by independent field experiments measuring hybrid fitness.

The effect of the environment on fitness in natural populations is a fundamental question in evolutionary biology. However, experimental manipulations of both environment and phenotype at the same time are rare. Thus, the relative importance of the competitive environment versus intrinsic organismal performance in shaping the location, height, and fluidity of fitness peaks and valleys remains largely unknown. Here, we experimentally tested the effect of competitor frequency on the complex fitness landscape driving adaptive radiation of a generalist and two trophic specialist pupfishes, a scale-eater and molluscivore, endemic to hypersaline lakes on San Salvador Island (SSI), Bahamas. We manipulated phenotypes, by generating 3407 F4/F5 lab-reared hybrids, and competitive environment, by altering the frequency of rare transgressive hybrids between field enclosures in two independent lake populations. We then tracked hybrid survival and growth rates across these four field enclosures for 3-11 months. In contrast to competitive speciation theory, we found no evidence that the frequency of hybrid phenotypes affected their survival. Instead, we observed a strikingly similar fitness landscape to a previous independent field experiment, each supporting multiple fitness peaks for generalist and molluscivore phenotypes and a large fitness valley isolating the divergent scale-eater phenotype. These features of the fitness landscape were stable across manipulated competitive environments, multivariate trait axes, and spatiotemporal heterogeneity. We suggest that absolute performance constraints and divergent gene regulatory networks shape macroevolutionary (interspecific) fitness landscapes in addition to microevolutionary (intraspecific) competitive dynamics. This interplay between organism and environment underlies static and dynamic features of the adaptive landscape.

Work that body: fin and body movements determine herbivore feeding performance within the natural reef environment.

Herbivorous fishes form a keystone component of reef ecosystems, yet the functional mechanisms underlying their feeding performance are poorly understood. In water, gravity is counter-balanced by buoyancy, hence fish are recoiled backwards after every bite they take from the substrate. To overcome this recoil and maintain contact with the algae covered substrate, fish need to generate thrust while feeding. However, the locomotory performance of reef herbivores in the context of feeding has hitherto been ignored. We used a three-dimensional high-speed video system to track mouth and body kinematics during in situ feeding strikes of fishes in the genus Zebrasoma, while synchronously recording the forces exerted on the substrate. These herbivores committed stereotypic and coordinated body and fin movements when feeding off the substrate and these movements determined algal biomass removed. Specifically, the speed of rapidly backing away from the substrate was associated with the magnitude of the pull force and the biomass of algae removed from the substrate per feeding bout. Our new framework for measuring biting performance in situ demonstrates that coordinated movements of the body and fins play a crucial role in herbivore foraging performance and may explain major axes of body and fin shape diversification across reef herbivore guilds.

Ecological divergence in sympatry causes gene misexpression in hybrids.

Ecological speciation occurs when reproductive isolation evolves as a byproduct of adaptive divergence between populations. Selection favouring gene regulatory divergence between species could result in transgressive levels of gene expression in F1 hybrids that may lower hybrid fitness. We combined 58 resequenced genomes with 124 transcriptomes to identify patterns of hybrid gene misexpression that may be driven by adaptive regulatory divergence within a young radiation of Cyprinodon pupfishes, which consists of a dietary generalist and two trophic specialists-a molluscivore and a scale-eater. We found more differential gene expression between closely related sympatric specialists than between allopatric generalist populations separated by 1,000 km. Intriguingly, 9.6% of genes that were differentially expressed between sympatric species were also misexpressed in F1 hybrids. A subset of these genes were in highly differentiated genomic regions and enriched for functions important for trophic specialization, including head, muscle and brain development. These regions also included genes that showed evidence of hard selective sweeps and were significantly associated with oral jaw length-the most rapidly diversifying skeletal trait in this radiation. Our results indicate that divergent ecological selection in sympatry can contribute to hybrid gene misexpression which may act as a reproductive barrier between nascent species.

Hydrodynamic Simulations of the Performance Landscape for Suction-Feeding Fishes Reveal Multiple Peaks for Different Prey Types.

The complex interplay between form and function forms the basis for generating and maintaining organismal diversity. Fishes that rely on suction-feeding for prey capture exhibit remarkable phenotypic and trophic diversity. Yet the relationships between fish phenotypes and feeding performance on different prey types are unclear, partly because the morphological, biomechanical, and hydrodynamic mechanisms that underlie suction-feeding are complex. Here we demonstrate a general framework to investigate the mapping of multiple phenotypic traits to performance by mapping kinematic variables to suction-feeding capacity. Using a mechanistic model of suction-feeding that is based on core physical principles, we predict prey capture performance across a broad range of phenotypic trait values, for three general prey types: mollusk-like prey, copepod-like prey, and fish-like prey. Mollusk-like prey attach to surfaces, copepod-like prey attempt to escape upon detecting the hydrodynamic disturbance produced by the predator, and fish-like prey attempt to escape when the predator comes within a threshold distance. This approach allowed us to evaluate suction-feeding performance for any combination of six key kinematic traits, irrespective of whether these trait combinations were observed in an extant species, and to generate a multivariate mapping of phenotype to performance. We used gradient ascent methods to explore the complex topography of the performance landscape for each prey type, and found evidence for multiple peaks. Characterization of phenotypes associated with performance peaks indicates that the optimal kinematic parameter range for suction-feeding on different prey types are narrow and distinct from each other, suggesting different functional constraints for the three prey types. These performance landscapes can be used to generate hypotheses regarding the distribution of extant species in trait space and their evolutionary trajectories toward adaptive peaks on macroevolutionary fitness landscapes.

Oral shelling within an adaptive radiation of pupfishes: Testing the adaptive function of a novel nasal protrusion and behavioural preference.

Dietary specialization on hard prey items, such as mollusks and crustaceans, is commonly observed in a diverse array of fish species. Many fish consume these types of prey by crushing the shell to consume the soft tissue within, but a few fishes extricate the soft tissue without breaking the shell using a method known as oral shelling. Oral shelling involves pulling a mollusc from its shell and it may be a way to subvert an otherwise insurmountable shell defence. However, the biomechanical requirements and potential adaptations for oral shelling are unknown. Here, we test the hypothesis that a novel nasal protrusion is an adaptation for oral shelling in the durophagous pupfish (Cyprinodon brontotheroides). We first demonstrate oral shelling in this species and then predict that a larger nasal protrusion would allow pupfish to consume larger snails. Durophagous pupfish are found within an endemic radiation of pupfish on San Salvador Island, Bahamas. We took advantage of closely related sympatric species and outgroups to test: (a) whether durophagous pupfish shell and consume more snails than other species, (b) if F1 and F2 durophagous hybrids consume similar amounts of snails as purebred durophagous pupfish, and (c) if nasal protrusion size in parental and hybrid populations increases the maximum size of consumed snails. We found that durophagous pupfish and their hybrids consumed the most snails, but did not find a strong association between nasal protrusion size and maximum snail size consumed within the parental or F2 hybrid population, suggesting that the size of their novel nasal protrusion does not provide a major benefit in oral shelling. Instead, we suggest that the nasal protrusion may increase feeding efficiency, act as a sensory organ, or is a sexually selected trait, and that a strong feeding preference may be most important for oral shelling.