Diabolical survival in Death Valley: recent pupfish colonization, gene flow and genetic assimilation in the smallest species range on earth

Inferring the evolution of reproductive isolation in a lineage of fossil threespine stickleback, Gasterosteus doryssus

Darwin attributed the absence of species transitions in the fossil record to his hypothesis that speciation occurs within isolated habitat patches too geographically restricted to be captured by fossil sequences. Mayr's peripatric speciation model added that such speciation would be rapid, further explaining missing evidence of diversification. Indeed, Eldredge and Gould's original punctuated equilibrium model combined Darwin's conjecture, Mayr's model and 124 years of unsuccessfully sampling the fossil record for transitions. Observing such divergence, however, could illustrate the tempo and mode of evolution during early speciation. Here, we investigate peripatric divergence in a Miocene stickleback fish, Gasterosteus doryssus. This lineage appeared and, over approximately 8000 generations, evolved significant reduction of 12 of 16 traits related to armour, swimming and diet, relative to its ancestral population. This was greater morphological divergence than we observed between reproductively isolated, benthic-limnetic ecotypes of extant Gasterosteus aculeatus. Therefore, we infer that reproductive isolation was evolving. However, local extinction of G. doryssus lineages shows how young, isolated, speciating populations often disappear, supporting Darwin's explanation for missing evidence and revealing a mechanism behind morphological stasis. Extinction may also account for limited sustained divergence within the stickleback species complex and help reconcile speciation rate variation observed across time scales.

Strong species structure but weak geographical structure in demersal Lake Victoria cichlids

Studying phenotypic and genetic differentiation between very young species can be very informative with regard to learning about processes of speciation. Identifying and characterizing genetic species structure and distinguishing it from spatial genetic structure within a species is a prerequisite for this and is often not given sufficient attention. Young radiations of cichlid fish are classical speciation study systems. However, it is only during the past decade that population genomics based on next-generation sequencing has begun to provide the power to resolve species and distinguish speciation from spatial population structure for the youngest of these radiations. The Lake Victoria haplochromine cichlids constitute the youngest large cichlid fish radiation, probably <20,000 years old. Earlier work showed that communities of rocky reef cichlids are composed of many reciprocally monophyletic species despite their very recent origins. Here, we build on this work by studying assemblages of offshore demersal cichlids, adding analyses of within-species spatial structure to the sympatric species structure. We sampled seven multispecies communities along a 6-km-long transect from one side of the Mwanza Gulf to the other side. We investigated whether phenotypically diagnosed putative species are reciprocally monophyletic and whether such monophyly is stable across species geographic ranges. We show that all species are genetically strongly differentiated in sympatry, that they are reciprocally monophyletic, and that monophyly is stable across distribution ranges. We found significant differentiation between geographically distinct populations in two species, but no or weak isolation by distance. We further found subtle but significant morphological differences between all species and a linear relationship between genomic and morphological distance which suggests that differences in morphology begin to accumulate after speciation has already affected genome-wide restrictions of gene flow.

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 (F_ROH = 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.

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.

Testing for genetic assimilation with phylogenetic comparative analysis: Conceptual, methodological, and statistical considerations

Genetic assimilation is a process that leads to reduced phenotypic plasticity during adaptation to novel conditions, a potentially important phenomenon under global environmental change. Null expectations when testing for genetic assimilation, however, are not always clear. For instance, the statistical artifact of regression to the mean could bias us toward detecting genetic assimilation when it has not occurred. Likewise, the specific mechanism underlying plasticity expression may affect null expectations under neutral evolution. We used macroevolutionary numerical simulations to examine both of these important issues and their interaction, varying whether plasticity evolves, the evolutionary mechanism, trait measurement error, and experimental design. We also modified an existing reaction norm correction method to account for phylogenetic nonindependence. We found (1) regression to the mean is pervasive and can generate spurious support for genetic assimilation; (2) experimental design and post hoc correction can minimize this spurious effect; and (3) neutral evolution can produce patterns consistent with genetic assimilation without constraint or selection, depending on the mechanism of plasticity expression. Additionally, we reanalyzed published macroevolutionary data supporting genetic assimilation, and found that support was reduced after proper correction. Considerable caution is thus required whenever investigating genetic assimilation and reaction norm evolution at macroevolutionary scales.

Fish out of water: Genomic insights into persistence of rainbowfish populations in the desert

How populations of aquatic fauna persist in extreme desert environments is an enigma. Individuals often breed and disperse during favorable conditions. Theory predicts that adaptive capacity should be low in small populations, such as in desert fishes. We integrated satellite-derived surface water data and population genomic diversity from 20,294 single-nucleotide polymorphisms across 344 individuals to understand metapopulation persistence of the desert rainbowfish (Melanotaenia splendida tatei) in central Australia. Desert rainbowfish showed very small effective population sizes, especially at peripheral populations, and low connectivity between river catchments. Yet, there was no evidence of population-level inbreeding and a signal of possible adaptive divergence associated with aridity was detected. Candidate genes for local adaptation included functions related to environmental cues and stressful conditions. Eco-evolutionary modeling showed that positive selection in refugial subpopulations combined with connectivity during flood periods can enable retention of adaptive diversity. Our study suggests that adaptive variation can be maintained in small populations and integrate with neutral metapopulation processes to allow persistence in the desert.

The time course of molecular acclimation to seawater in a euryhaline fish

The Arabian pupfish, Aphanius dispar, is a euryhaline fish inhabiting both inland nearly-freshwater desert ponds and highly saline Red Sea coastal lagoons of the Arabian Peninsula. Desert ponds and coastal lagoons, located respectively upstream and at the mouths of dry riverbeds (“wadies”), have been found to potentially become connected during periods of intense rainfall, which could allow the fish to migrate between these different habitats. Flash floods would therefore flush Arabian pupfish out to sea, requiring a rapid acclimation to a greater than 40 ppt change in salinity. To investigate the molecular pathways of salinity acclimation during such events, a Red Sea coastal lagoon and a desert pond population were sampled, with the latter exposed to a rapid increase in water salinity. Changes in branchial gene expression were investigated via genome-wide transcriptome measurements over time from 6 h to 21 days. The two natural populations displayed basal differences in genes related to ion transport, osmoregulation and immune system functions. These mechanisms were also differentially regulated in seawater transferred fish, revealing their crucial role in long-term adaptation. Other processes were only transiently activated shortly after the salinity exposure, including cellular stress response mechanisms, such as molecular chaperone synthesis and apoptosis. Tissue remodelling processes were also identified as transient, but took place later in the timeline, suggesting their importance to long-term acclimation as they likely equip the fish with lasting adaptations to their new environment. The alterations in branchial functional pathways displayed by Arabian pupfish in response to salinity increases are diverse. These reveal a large toolkit of molecular processes important for adaptation to hyperosmolarity that allow for successful colonization to a wide variety of different habitats.

The endangered White Sands pupfish (Cyprinodon tularosa) genome reveals low diversity and heterogenous patterns of differentiation

The White Sands pupfish (Cyprinodon tularosa), endemic to New Mexico in Southwestern North America, is of conservation concern due in part to invasive species, chemical pollution, and groundwater withdrawal. Herein, we developed a draft reference genome and use it to provide biological insights into the evolution and conservation of C. tularosa. We used our assembly to localize microsatellite markers previously used to demarcate evolutionary significant units (ESU), quantified genomic divergence and transposable element profiles between species, and compared C. tularosa genomic diversity related species. Our de novo assembly of PacBio Sequel II error-corrected reads resulted in a 1.08 Gb draft genome with a contig N50 of 1.4 Mb and 25,260 annotated protein coding genes, including 95% of the expected Actinopterygii conserved complete single-copy orthologues. Many of the C. tularosa microsatellite markers used for conservation assessments fell within, or near, genes and exhibited a pattern of increased heterozygosity near genic areas compared to those in intergenic regions. Nuclear alignments between these two species revealed 193 genes contained in rapidly diverging tracts; transposable element profiles were largely concordant and suggest a shared, rapid expansion of LINE and Gypsy elements. Genome-wide heterozygosity was markedly lower in C. tularosa compared to estimates from other related species, probably because of smaller long-term effective population sizes constrained by their isolated and limited habitat. Overall, these inferences provide new insights into C. tularosa that should help inform future management efforts.

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.

Optimizing Sample Size for Population Genomic Study in a Global Invasive Lady Beetle, Harmonia Axyridis

Finding optimal sample sizes is critical for the accurate estimation of genetic diversity of large invasive populations. Based on previous studies, we hypothesized that a minimal sample size of 3-8 individuals is sufficient to dissect the population architecture of the harlequin lady beetle, Harmonia axyridis, a biological control agent and an invasive alien species. Here, equipped with a type IIB endonuclease restriction site-associated (2b-RAD) DNA sequencing approach, we identified 13,766 and 13,929 single nucleotide polymorphisms (SNPs), respectively, among native and invasive H. axyridis populations. With this information we simulated populations using a randomly selected 3000 SNPs and a subset of individuals. From this simulation we finally determined that six individuals is the minimum sample size required for the accurate estimation of intra- and inter-population genetic diversity within and across H. axyridis populations. Our findings provide an empirical advantage for population genomic studies of H. axyridis in particular and suggest useful tactics for similar studies on multicellular organisms in general.

Hormones, developmental plasticity, and adaptive evolution: Endocrine flexibility as a catalyst for 'plasticity-first' phenotypic divergence

Explaining how populations adapt to environments is among the foremost objectives of evolutionary theory. Over generations, natural selection impels the phenotypic distribution of a population based on individual variation in phenotype and fitness. However, environmental conditions can also shape how individuals develop within their lifetime to influence which phenotypes are expressed in a population. It has been proposed that such environmentally-initiated phenotypic variation - also termed developmental plasticity - may enable adaptive evolution under some scenarios. As dynamic regulators of development and phenotypic expression, hormones are important physiological mediators of developmental plasticity. Patterns of hormone secretion, hormone transport, and the sensitivity of tissues to hormones can each be altered by environmental conditions, and understanding how endocrine regulation shapes phenotypic development in an ecologically-relevant context has much to contribute toward clarifying the role of plasticity in evolutionary adaptation. This article explores how the environmental sensitivity of endocrine regulation may facilitate 'plasticity-first' evolution by generating phenotypic variants that precede adaptation to altered or novel environments. Predictions arising from 'plasticity-first' evolution are examined in the context of thyroid hormone mediation of morphological plasticity in Cyprinodon pupfishes from the Death Valley region of California and Nevada, USA. This clade of extremophile fishes diversified morphologically over the last ~20,000 years, and observations that some populations experienced contemporary phenotypic differentiation under recent habitat change provide evidence that hormone-mediate plasticity preceded genetic assimilation of morphology in one of the region's species: the Devils Hole pupfish, Cyprinodon diabolis. This example illustrates how conceptualizing hormones not only as regulators of homeostasis, but also as developmental intermediaries between environment conditions and phenotypic variation at the individual-, population-, and species-levels can enrich our understanding of endocrine regulation both as a facilitator of phenotypic change under shifting environments, and as important proximate mechanisms that may initiate 'plasticity-first' evolutionary adaptation.

Phenotypic diversity in an endangered freshwater fish Squalius microlepis (Actinopterygii, Leuciscidae)

Squalius microlepis was examined from recent and historical collections within the known range of the species with special emphasis on intraspecific variability and variations, and compared to its closest relative species S. tenellus (in total, 193 specimens; 33 absolute and 52 proportional measurements and ratios, and 12 counts including vertebrae). Squalius tenellus was perfectly differentiated in all statistical analyses and can be diagnosed by 76-95 (vs. 64-80) scales in lateral series, 68-83 (vs. 58-77) lateral-line scales, (17)18-20 (vs. 13-16(17)) scales above lateral line, and (7)8-10 (vs. 4-7) scales below lateral line. Squalius microlepis was morphologically heterogeneous, with two phenotypes readily distinguishable (phenotype 1 corresponding to S. microlepis s. str. as defined by its lectotype) by a combination of many characters; those contributing most to the discrimination were number of gill rakers, length of lower jaw (% interorbital width), and head length (% SL). Only phenotype 1 was found in the Ričina-Prološko Blato-Vrljika karst system; most of the specimens from the lower Matica and the Tihaljina-Trebižat karst system were identified as phenotype 2; the sample from karstic poljes near Vrgorac contained both phenotype 1 and 2, and individuals of intermediate morphology. As very limited molecular data exist on the two phenotypes of S. microlepis, we refrain from any taxonomic conclusions until new molecular approaches (and new markers) are used. We also report on a dramatic reduction of the area of distribution and abundance of S. microlepis in recent years.

The paradox behind the pattern of rapid adaptive radiation: how can the speciation process sustain itself through an early burst?

Rapid adaptive radiation poses a distinct question apart from speciation and adaptation: what happens after one speciation event? That is, how are some lineages able to continue speciating through a rapid burst? This question connects global macroevolutionary patterns to microevolutionary processes. Here we review major features of rapid radiations in nature and their mismatch with theoretical models and what is currently known about speciation mechanisms. Rapid radiations occur on three major diversification axes - species richness, phenotypic disparity, and ecological diversity - with exceptional outliers on each axis. The paradox is that the hallmark early stage of adaptive radiation, a rapid burst of speciation and niche diversification, is contradicted by most existing speciation models which instead predict continuously decelerating speciation rates and niche subdivision through time. Furthermore, while speciation mechanisms such as magic traits, phenotype matching, and physical linkage of co-adapted alleles promote speciation, it is often not discussed how these mechanisms could promote multiple speciation events in rapid succession. Additional mechanisms beyond ecological opportunity are needed to understand how rapid radiations occur. We review the evidence for five emerging theories: 1) the 'transporter' hypothesis: introgression and the ancient origins of adaptive alleles, 2) the 'signal complexity' hypothesis: the dimensionality of sexual traits, 3) the connectivity of fitness landscapes, 4) 'diversity begets diversity', and 5) flexible stem/'plasticity first'. We propose new questions and predictions to guide future work on the mechanisms underlying the rare origins of rapid radiation.

How to Investigate the Origins of Novelty: Insights Gained from Genetic, Behavioral, and Fitness Perspectives

Biologists are drawn to the most extraordinary adaptations in the natural world, often referred to as evolutionary novelties, yet rarely do we understand the microevolutionary context underlying the origins of novel traits, behaviors, or ecological niches. Here we discuss insights gained into the origins of novelty from a research program spanning biological levels of organization from genotype to fitness in Caribbean pupfishes. We focus on a case study of the origins of novel trophic specialists on San Salvador Island, Bahamas and place this radiation in the context of other rapid radiations. We highlight questions that can be addressed about the origins of novelty at different biological levels, such as measuring the isolation of novel phenotypes on the fitness landscape, locating the spatial and temporal origins of adaptive variation contributing to novelty, detecting dysfunctional gene regulation due to adaptive divergence, and connecting behaviors with novel traits. Evolutionary novelties are rare, almost by definition, and we conclude that integrative case studies can provide insights into this rarity relative to the dynamics of adaptation to more common ecological niches and repeated parallel speciation, such as the relative isolation of novel phenotypes on fitness landscapes and the transient availability of ecological, genetic, and behavioral opportunities.

Local ancestry analysis reveals genomic convergence in extremophile fishes

The molecular basis of convergent phenotypes is often unknown. However, convergence at a genomic level is predicted when there are large population sizes, gene flow among diverging lineages or strong genetic constraints. We used whole-genome resequencing to investigate genomic convergence in fishes ( Poecilia spp.) that have repeatedly colonized hydrogen sulfide (H₂S)-rich environments in Mexico. We identified genomic similarities in both single nucleotide polymorphisms (SNPs) and structural variants (SVs) among independently derived sulfide spring populations, with approximately 1.2% of the genome being shared among sulfidic ecotypes. We compared these convergent genomic regions to candidate genes for H₂S adaptation identified from transcriptomic analyses and found that a significant proportion of these candidate genes (8%) were also in regions where sulfidic individuals had similar SNPs, while only 1.7% were in regions where sulfidic individuals had similar SVs. Those candidate genes included genes involved in sulfide detoxification, the electron transport chain (the main toxicity target of H₂S) and other processes putatively important for adaptation to sulfidic environments. Regional genomic similarity across independent populations exposed to the same source of selection is consistent with selection on standing variation or introgression of adaptive alleles across divergent lineages. However, combined with previous analyses, our data also support that adaptive changes in mitochondrially encoded subunits arose independently via selection on de novo mutations. Pressing questions remain on what conditions ultimately facilitate the independent rise of adaptive alleles at the same loci in separate populations, and thus, the degree to which evolution is repeatable or predictable. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Plasticity-led evolution: evaluating the key prediction of frequency-dependent adaptation

Plasticity-led evolution occurs when a change in the environment triggers a change in phenotype via phenotypic plasticity, and this pre-existing plasticity is subsequently refined by selection into an adaptive phenotype. A critical, but largely untested prediction of plasticity-led evolution (and evolution by natural selection generally) is that the rate and magnitude of evolutionary change should be positively associated with a phenotype's frequency of expression in a population. Essentially, the more often a phenotype is expressed and exposed to selection, the greater its opportunity for adaptive refinement. We tested this prediction by competing against each other spadefoot toad tadpoles from different natural populations that vary in how frequently they express a novel, environmentally induced carnivore ecomorph. As expected, laboratory-reared tadpoles whose parents were derived from populations that express the carnivore ecomorph more frequently were superior competitors for the resource for which this ecomorph is specialized-fairy shrimp. These tadpoles were better at using this resource both because they were more efficient at capturing and consuming shrimp and because they produced more exaggerated carnivore traits. Moreover, they exhibited these more carnivore-like features even without experiencing the inducing cue, suggesting that this ecomorph has undergone an extreme form of plasticity-led evolution-genetic assimilation. Thus, our findings provide evidence that the frequency of trait expression drives the magnitude of adaptive refinement, thereby validating a key prediction of plasticity-led evolution specifically and adaptive evolution generally.

An embryonic staging series up to hatching for Cyprinodon variegatus: An emerging fish model for developmental, evolutionary, and ecological research

Using multiple taxa to research development is necessary for making general conclusions about developmental patterns and mechanisms. We present a staging series for Cyprinodon variegatus as a basis for further study of the developmental biology of fishes in the genus Cyprinodon and for comparative work on teleost fishes beyond the standard models. Cyprinodon are small, euryhaline fishes, widely distributed in fresh, brackish, and hypersaline waters of southern and eastern North America. Cyprinodontids are closely related to fundulids, providing a comparative reference point to the embryological model, Fundulus heteroclitus. Ecologists and evolutionary biologists commonly study Cyprinodon, and we have been using Cyprinodon to study skull variation and its genetic basis among closely related species. We divided embryonic development of C. variegatus into 34 morphologically identifiable stages. We reference our staging series to that already defined for a related model species, Oryzias latipes (medaka) that is studied by a large community of researchers. We provide a description of the early chondrogenesis and ossification of skull and caudal fin bones during the latter stages of embryonic development. We show that Cyprinodon are tractable for studying development. Eggs can be obtained easily from breeding pairs and our study provides a staging system to facilitate future developmental studies.

The genomic landscape at a late stage of stickleback speciation: High genomic divergence interspersed by small localized regions of introgression

Speciation is a continuous process and analysis of species pairs at different stages of divergence provides insight into how it unfolds. Previous genomic studies on young species pairs have revealed peaks of divergence and heterogeneous genomic differentiation. Yet less known is how localised peaks of differentiation progress to genome-wide divergence during the later stages of speciation in the presence of persistent gene flow. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence builds up during speciation. However, attention has largely focused on young postglacial species pairs, with little knowledge of the genomic signatures of divergence and introgression in older stickleback systems. The Japanese stickleback species pair, composed of the Pacific Ocean three-spined stickleback (Gasterosteus aculeatus) and the Japan Sea stickleback (G. nipponicus), which co-occur in the Japanese islands, is at a late stage of speciation. Divergence likely started well before the end of the last glacial period and crosses between Japan Sea females and Pacific Ocean males result in hybrid male sterility. Here we use coalescent analyses and Approximate Bayesian Computation to show that the two species split approximately 0.68-1 million years ago but that they have continued to exchange genes at a low rate throughout divergence. Population genomic data revealed that, despite gene flow, a high level of genomic differentiation is maintained across the majority of the genome. However, we identified multiple, small regions of introgression, occurring mainly in areas of low recombination rate. Our results demonstrate that a high level of genome-wide divergence can establish in the face of persistent introgression and that gene flow can be localized to small genomic regions at the later stages of speciation with gene flow.

Parallel evolution of gene expression between trophic specialists despite divergent genotypes and morphologies

Parallel evolution of gene expression commonly underlies convergent niche specialization, but parallel changes in expression could also underlie divergent specialization. We investigated divergence in gene expression and whole-genome genetic variation across three sympatric Cyprinodon pupfishes endemic to San Salvador Island, Bahamas. This recent radiation consists of a generalist and two derived specialists adapted to novel niches: a scale-eating and a snail-eating pupfish. We sampled total mRNA from all three species at two early developmental stages and compared gene expression with whole-genome genetic differentiation among all three species in 42 resequenced genomes. Eighty percent of genes that were differentially expressed between snail-eaters and generalists were up or down regulated in the same direction between scale-eaters and generalists; however, there were no fixed variants shared between species underlying these parallel changes in expression. Genes showing parallel evolution of expression were enriched for effects on metabolic processes, whereas genes showing divergent expression were enriched for effects on cranial skeleton development and pigment biosynthesis, reflecting the most divergent phenotypes observed between specialist species. Our findings reveal that even divergent niche specialists may exhibit convergent adaptation to higher trophic levels through shared genetic pathways. This counterintuitive result suggests that parallel evolution in gene expression can accompany divergent ecological speciation during adaptive radiation.

Best available science still supports an ancient common origin of Devils Hole and Devils Hole pupfish

The age of DHP and how pupfish colonized Devils Hole have always been a topic of interest. Recently, two different publications (Martin, Crawford, Turner, & Simons, & Sağlam et al., ) tackled this issue using genomic data sets and demographic models but came to widely different conclusions. In their comment, Martin and Höhne () argue that our results (Sağlam et al., ) were misleading because we used inappropriate calibration information and biased a priori assumptions. They then re-analysed our data using a "biologically informed" mutation rate prior and concluded that our data support a much younger age of DHP (12.6 kya) as opposed to 60 kya reported in our study. Below we will summarize why their arguments do not hold up and explore some of the inconsistencies between their claims and what was actually presented in our study. Furthermore, we will demonstrate their re-analyses provide no new information compared to what was presented in our original manuscript and reinforce our estimate of a 60 kya divergence of DHP as outweighing competing hypotheses.

Fantastic Beasts and How To Sequence Them: Ecological Genomics for Obscure Model Organisms

The application of genomic approaches to 'obscure model organisms' (OMOs), meaning species with no prior genomic resources, enables increasingly sophisticated studies of the genomic basis of evolution, acclimatization, and adaptation in real ecological contexts. I consider here ecological questions that can be addressed using OMOs, and indicate optimal sequencing and data-handling solutions for each case. With this I hope to promote the diversity of OMO-based projects that would capitalize on the peculiarities of the natural history of OMOs and could feasibly be completed within the scope of a single PhD thesis.

Genome-scale data reveal that endemic Poecilia populations from small sulphidic springs display no evidence of inbreeding

Populations with limited ranges can be highly vulnerable to changes in their environment and are, thus, of high conservation concern. Populations that experience human-induced range reductions are often highly inbred and lack genetic diversity, but it is unknown whether this is also the case for populations with naturally small ranges. The fishes Poecilia sulphuraria (listed as critically endangered) and Poecilia thermalis, which are endemic to small hydrogen sulphide-rich springs in southern Mexico, are examples of such populations with inherently small habitats. We used geometric morphometrics and population genetics to quantify phenotypic and genetic variation within and among two populations of P. sulphuraria and one population of P. thermalis. Principal component analyses revealed phenotypic and genetic differences among the populations. Evidence for inbreeding was low compared to populations that have undergone habitat reduction. The genetic data were also used to infer the demographic history of these populations to obtain estimates for effective population sizes and migration rates. Effective population sizes were large given the small habitats of these populations. Our results imply that these three endemic extremophile populations should each be considered separately for conservation purposes. Additionally, this study suggests that populations in naturally small habitats may have lower rates of inbreeding and higher genetic diversity than expected, and therefore may be better equipped to handle environmental perturbations than anticipated. We caution, however, that the inferred lack of inbreeding and the large effective population sizes could potentially be a result of colonization by genetically diverse ancestors.

Adaptive introgression from distant Caribbean islands contributed to the diversification of a microendemic adaptive radiation of trophic specialist pupfishes

Rapid diversification often involves complex histories of gene flow that leave variable and conflicting signatures of evolutionary relatedness across the genome. Identifying the extent and source of variation in these evolutionary relationships can provide insight into the evolutionary mechanisms involved in rapid radiations. Here we compare the discordant evolutionary relationships associated with species phenotypes across 42 whole genomes from a sympatric adaptive radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas and several outgroup pupfish species in order to understand the rarity of these trophic specialists within the larger radiation of Cyprinodon. 82% of the genome depicts close evolutionary relationships among the San Salvador Island species reflecting their geographic proximity, but the vast majority of variants fixed between specialist species lie in regions with discordant topologies. Top candidate adaptive introgression regions include signatures of selective sweeps and adaptive introgression of genetic variation from a single population in the northwestern Bahamas into each of the specialist species. Hard selective sweeps of genetic variation on San Salvador Island contributed 5 times more to speciation of trophic specialists than adaptive introgression of Caribbean genetic variation; however, four of the 11 introgressed regions came from a single distant island and were associated with the primary axis of oral jaw divergence within the radiation. For example, standing variation in a proto-oncogene (ski) known to have effects on jaw size introgressed into one San Salvador Island specialist from an island 300 km away approximately 10 kya. The complex emerging picture of the origins of adaptive radiation on San Salvador Island indicates that multiple sources of genetic variation contributed to the adaptive phenotypes of novel trophic specialists on the island. Our findings suggest that a suite of factors, including rare adaptive introgression, may be necessary for adaptive radiation in addition to ecological opportunity.

Characterization and phylogenetic analysis of complete mitochondrial genomes for two desert cyprinodontoid fishes, Empetrichthys latos and Crenichthys baileyi

The Pahrump poolfish (Empetrichthys latos) and White River springfish (Crenichthys baileyi) are small-bodied teleost fishes (order Cyprinodontiformes) endemic to the arid Great Basin and Mojave Desert regions of western North America. These taxa survive as small, isolated populations in remote streams and springs and evolved to tolerate extreme conditions of high temperature and low dissolved oxygen. Both species have experienced severe population declines over the last 50-60years that led to some subspecies being categorized with protected status under the U.S. Endangered Species Act. Here we report the first sequencing of the complete mitochondrial DNA genomes for both E. l. latos and the moapae subspecies of C. baileyi. Complete mitogenomes of 16,546bp nucleotides were obtained from two E. l. latos individuals collected from introduced populations at Spring Mountain Ranch State Park and Shoshone Ponds Natural Area, Nevada, USA, while a single mitogenome of 16,537bp was sequenced for C. b. moapae. The mitogenomes of both species contain 13 protein-encoding genes, twenty-two tRNAs, and two rRNAs (12S and 18S) following the syntenic arrangement typical of Actinopterygiian fish mitogenomes, as well as D-loop control regions of 858bp for E. latos and 842bp for C. baileyi moapae. The two E. latos individuals exhibited only 0.0181% nucleotide sequence divergence across the entire mitogenome, implying little intraspecific mtDNA genetic variation. Comparative phylogenetic analysis of the poolfish and springfish mitochondrial genomes to available mitogenomes of other Cyprinodontoid fishes confirmed the close relationship of these oviparous Empetrichthys and Crenichthys genera to the viviparous goodeid fishes of central Mexico, and showed the combined clade of these fishes to be a sister group to the Profundulidae killifishes. Despite several significant life history and morphological differences between the Empetrichthyinae and Goodienae, estimates of evolutionary genetic distances using two partial regions of mtDNA point to inclusion of the Empetrichthys and Crenichthys genera within the family Goodeidae along with the goodeid fishes of central Mexico.

Novel Candidate Genes Underlying Extreme Trophic Specialization in Caribbean Pupfishes

The genetic changes responsible for evolutionary transitions from generalist to specialist phenotypes are poorly understood. Here we examine the genetic basis of craniofacial traits enabling novel trophic specialization in a sympatric radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas. This recent radiation consists of a generalist species and two novel specialists: a small-jawed "snail-eater" and a large-jawed "scale-eater." We genotyped 12 million single nucleotide polymorphisms (SNPs) by whole-genome resequencing of 37 individuals of all three species from nine populations and integrated genome-wide divergence scans with association mapping to identify divergent regions containing putatively causal SNPs affecting jaw size-the most rapidly diversifying trait in this radiation. A mere 22 fixed variants accompanied extreme ecological divergence between generalist and scale-eater species. We identified 31 regions (20 kb) containing variants fixed between specialists that were significantly associated with variation in jaw size which contained 11 genes annotated for skeletal system effects and 18 novel candidate genes never previously associated with craniofacial phenotypes. Six of these 31 regions showed robust signs of hard selective sweeps after accounting for demographic history. Our data are consistent with predictions based on quantitative genetic models of adaptation, suggesting that the effect sizes of regions influencing jaw phenotypes are positively correlated with distance between fitness peaks on a complex adaptive landscape.

Rapid genetic and morphologic divergence between captive and wild populations of the endangered Leon Springs pupfish, Cyprinodon bovinus

The Leon Springs pupfish (Cyprinodon bovinus) is an endangered species currently restricted to a single desert spring and a separate captive habitat in southwestern North America. Following establishment of the captive population from wild stock in 1976, the wild population has undergone natural population size fluctuations, intentional culling to purge genetic contamination from an invasive congener (Cyprinodon variegatus) and augmentation/replacement of wild fish from the captive stock. A severe population decline following the most recent introduction of captive fish prompted us to examine whether the captive and wild populations have differentiated during the short time they have been isolated from one another. If so, the development of divergent genetic and/or morphologic traits between populations could contribute to a diminished ability of fish from one location to thrive in the other. Examination of genomewide single nucleotide polymorphisms and morphologic variation revealed no evidence of residual C. variegatus characteristics in contemporary C. bovinus samples. However, significant genetic and morphologic differentiation was detected between the wild and captive populations, some of which might reflect local adaptation. Our results indicate that genetic and physical characteristics can diverge rapidly between isolated subdivisions of managed populations, potentially compromising the value of captive stock for future supplementation efforts. In the case of C. bovinus, our findings underscore the need to periodically inoculate the captive population with wild genetic material to help mitigate genetic, and potentially morphologic, divergence between them and also highlight the utility of parallel morphologic and genomic evaluation to inform conservation management planning.

The genetic architecture of novel trophic specialists: larger effect sizes are associated with exceptional oral jaw diversification in a pupfish adaptive radiation

The genetic architecture of adaptation is fundamental to understanding the mechanisms and constraints governing diversification. However, most case studies focus on loss of complex traits or parallel speciation in similar environments. It is still unclear how the genetic architecture of these local adaptive processes compares to the architecture of evolutionary transitions contributing to morphological and ecological novelty. Here, we identify quantitative trait loci (QTL) between two trophic specialists in an excellent case study for examining the origins of ecological novelty: a sympatric radiation of pupfishes endemic to San Salvador Island, Bahamas, containing a large-jawed scale-eater and a short-jawed molluscivore with a skeletal nasal protrusion. These specialized niches and trophic traits are unique among over 2000 related species. Measurements of the fitness landscape on San Salvador demonstrate multiple fitness peaks and a larger fitness valley isolating the scale-eater from the putative ancestral intermediate phenotype of the generalist, suggesting that more large-effect QTL should contribute to its unique phenotype. We evaluated this prediction using an F2 intercross between these specialists. We present the first linkage map for pupfishes and detect significant QTL for sex and eight skeletal traits. Large-effect QTL contributed more to enlarged scale-eater jaws than the molluscivore nasal protrusion, consistent with predictions from the adaptive landscape. The microevolutionary genetic architecture of large-effect QTL for oral jaws parallels the exceptional diversification rates of oral jaws within the San Salvador radiation observed over macroevolutionary timescales and may have facilitated exceptional trophic novelty in this system.

Endocrine and metabolic impacts of warming aquatic habitats: differential responses between recently isolated populations of a eurythermal desert pupfish

Temperatures of inland aquatic habitats are increasing with climate change, and understanding how fishes respond physiologically to thermal stress will be crucial for identifying species most susceptible to these changes. Desert fishes may be particularly vulnerable to rising temperatures because many species occupy only a fraction of their historical range and occur in habitats with already high temperatures. Here, we examined endocrine and metabolic responses to elevated temperature in Amargosa pupfish, Cyprinodon nevadensis amargosae. We studied C. n. amargosae from two habitats with distinct thermal conditions: the Amargosa River, which experiences diurnally and seasonally variable temperatures (0.2-40°C); and Tecopa Bore, a spring and marsh fed by hot groundwater (47.5°C) from an artesian borehole. These allopatric populations differ in morphology, and prior evidence suggests that temperature might contribute to these differences via altered thyroid hormone (TH) regulation of morphological development. Here, we document variation in hepatic iodothyronine deiodinase type 2 (dio2) and type 3 (dio3) and TH receptor β (trβ) gene transcript abundance between the Amargosa River and Tecopa Bore wild populations. Fish from these populations acclimated to 24 or 34°C retained differences in hepatic dio2, dio3 and trβ mRNAs and also varied in transcripts encoding the TH membrane transporters monocarboxylate transporter 8 (mct8) and organic anion-transporting protein 1c1 (oatp1c1). Tecopa Bore pupfish also exhibited higher dio2 and trβ mRNA levels in skeletal muscle relative to Amargosa River fish. Muscle citrate synthase activity was lower at 34°C for both populations, whereas lactate dehydrogenase activity and lactate dehydrogenase A-chain (ldhA) transcripts were both higher and 3,5,3'-triiodothryonine responsive in Tecopa Bore pupfish only. These findings reveal that local population variation and thermal experience interact to shape how pupfish respond to elevated temperatures, and point to the need to consider such interactions in management actions for desert fishes under a changing climate.

Successful population establishment from small introductions appears to be less common than believed

Although small populations are at high risk of extinction, there are regular reports in the scientific literature of purported small, isolated, persistent populations. One source of evidence of the viability of small populations comes from the alleged successful introduction of species to areas outside their original range from introductions of few individuals. We reviewed the examples from introduction compendia on deliberate translocations of birds, and the original sources, to identify and evaluate purported examples of successful establishments from small introductions. We found 23 purportedly successful introductions from few (<30) individuals. After assessing original sources, we found that two of the claims were substantiated; the rest were ambiguous or could be rejected as examples, primarily due to a lack of evidence in original sources of the number of birds released and because of supplemental individuals from other releases, releases in nearby regions, and the possibility of natural invasion. Our results suggest that reports of successful establishment of birds from introductions of few individuals have been overstated. These results strengthen the relationship previously reported between propagule pressure and likelihood of establishment, and support the lack of viability of small populations presumed by population theory. We suggest that analyses of introduction failure and success would benefit from excluding studies where introduction effort is unknown or unreliably documented.

The cryptic origins of evolutionary novelty: 1000-fold faster trophic diversification rates without increased ecological opportunity or hybrid swarm

Ecological opportunity is frequently proposed as the sole ingredient for adaptive radiation into novel niches. An additional trigger may be genome-wide hybridization resulting from "hybrid swarm." However, these hypotheses have been difficult to test due to the rarity of comparable control environments lacking adaptive radiations. Here I exploit such a pattern in microendemic radiations of Caribbean pupfishes. I show that a sympatric three species radiation on San Salvador Island, Bahamas diversified 1445 times faster than neighboring islands in jaw length due to the evolution of a novel scale-eating adaptive zone from a generalist ancestral niche. I then sampled 22 generalist populations on seven neighboring islands and measured morphological diversity, stomach content diversity, dietary isotopic diversity, genetic diversity, lake/island areas, macroalgae richness, and Caribbean-wide patterns of gene flow. None of these standard metrics of ecological opportunity or gene flow were associated with adaptive radiation, except for slight increases in macroalgae richness. Thus, exceptional trophic diversification is highly localized despite myriad generalist populations in comparable environmental and genetic backgrounds. This study provides a strong counterexample to the ecological and hybrid swarm theories of adaptive radiation and suggests that diversification of novel specialists on a sparse fitness landscape is constrained by more than ecological opportunity and gene flow.

Sequencing and characterization of the complete mitochondrial genome of the endangered Devils Hole pupfish Cyprinodon diabolis (Cyprinodontiformes: Cyprinodontidae)

The Devils Hole pupfish (Cyprinodon diabolis) is an imperiled cypinodontid fish endemic to Devils Hole, a single groundwater-fed pool in western North America characterized by temperature and oxygen conditions lethal to most fishes. In this study, we obtained the 16,499 base pair (bp) mitochondrial DNA sequence of C. diabolis. The mitogenome encodes 13 polypeptides, 22 tRNAs, 12S and 18S rRNAs, and an 832 bp D-loop region, and has a nucleotide composition of A, 25.69%; T, 27.32%; G, 17.49%; and C, 29.50% (GC content: 46.99%). The availability of this mitogenome will facilitate evaluations of C. diabolis genetic structure for management and conservation of this iconic endangered species.

A New Species of the Bay Goby Genus Eucyclogobius, Endemic to Southern California: Evolution, Conservation, and Decline

A geographically isolated set of southern localities of the formerly monotypic goby genus Eucyclogobius is known to be reciprocally monophyletic and substantially divergent in mitochondrial sequence and nuclear microsatellite-based phylogenies relative to populations to the north along the California coast. To clarify taxonomic and conservation status, we conducted a suite of analyses on a comprehensive set of morphological counts and measures from across the range of Eucyclogobius and describe the southern populations as a new species, the Southern Tidewater Goby, Eucyclogobius kristinae, now separate from the Northern Tidewater Goby Eucyclogobius newberryi (Girard 1856). In addition to molecular distinction, adults of E. kristinae are diagnosed by: 1) loss of the anterior supratemporal lateral-line canals resulting in higher neuromast counts, 2) lower pectoral and branched caudal ray counts, and 3) sets of measurements identified via discriminant analysis. These differences suggest ecological distinction of the two species. Previous studies estimated lineage separation at 2-4 million years ago, and mitochondrial sequence divergence exceeds that of other recognized fish species. Fish from Santa Monica Artesian Springs (Los Angeles County) northward belong to E. newberryi; those from Aliso Creek (Orange County) southward constitute E. kristinae. The lagoonal habitat of Eucyclogobius has been diminished or degraded, leading to special conservation status at state and federal levels beginning in 1980. Habitat of the newly described species has been impacted by a range of anthropogenic activities, including the conversion of closing lagoons to open tidal systems in the name of restoration. In the last 30 years, E. kristinae has only been observed in nine intermittently occupied lagoonal systems in northern San Diego County; it currently persists in only three sites. Thus, the new species is in imminent danger of extinction and will require ongoing active management.

Context dependence in complex adaptive landscapes: frequency and trait-dependent selection surfaces within an adaptive radiation of Caribbean pupfishes

The adaptive landscape provides the foundational bridge between micro- and macroevolution. One well-known caveat to this perspective is that fitness surfaces depend on ecological context, including competitor frequency, traits measured, and resource abundance. However, this view is based largely on intraspecific studies. It is still unknown how context-dependence affects the larger features of peaks and valleys on the landscape which ultimately drive speciation and adaptive radiation. Here, I explore this question using one of the most complex fitness landscapes measured in the wild in a sympatric pupfish radiation endemic to San Salvador Island, Bahamas by tracking survival and growth of laboratory-reared F2 hybrids. I present new analyses of the effects of competitor frequency, dietary isotopes, and trait subsets on this fitness landscape. Contrary to expectations, decreasing competitor frequency increased survival only among very common phenotypes, whereas less common phenotypes rarely survived despite few competitors, suggesting that performance, not competitor frequency, shapes large-scale features of the fitness landscape. Dietary isotopes were weakly correlated with phenotype and growth, but did not explain additional survival variation. Nonlinear fitness surfaces varied substantially among trait subsets, revealing one-, two-, and three-peak landscapes, demonstrating the complexity of selection in the wild, even among similar functional traits.