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.

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.

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.

Searching for Sympatric Speciation in the Genomic Era.

Sympatric speciation illustrates how natural and sexual selection may create new species in isolation without geographic barriers. However, recent genomic reanalyses of classic examples of sympatric speciation reveal complex histories of secondary gene flow from outgroups into the radiation. In contrast, the rich theoretical literature on this process distinguishes among a diverse range of models based on simple genetic histories and different types of reproductive isolating barriers. Thus, there is a need to revisit how to connect theoretical models of sympatric speciation and their predictions to empirical case studies in the face of widespread gene flow. Here, theoretical differences among different types of sympatric speciation and speciation-with-gene-flow models are reviewed and summarized, and genomic analyses are proposed for distinguishing which models apply to case studies based on the timing and function of adaptive introgression. Investigating whether secondary gene flow contributed to reproductive isolation is necessary to test whether predictions of theory are ultimately borne out in nature.

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.

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.

Variation Across Mitochondrial Gene Trees Provides Evidence for Systematic Error: How Much Gene Tree Variation Is Biological?

The use of large genomic data sets in phylogenetics has highlighted extensive topological variation across genes. Much of this discordance is assumed to result from biological processes. However, variation among gene trees can also be a consequence of systematic error driven by poor model fit, and the relative importance of biological vs. methodological factors in explaining gene tree variation is a major unresolved question. Using mitochondrial genomes to control for biological causes of gene tree variation, we estimate the extent of gene tree discordance driven by systematic error and employ posterior prediction to highlight the role of model fit in producing this discordance. We find that the amount of discordance among mitochondrial gene trees is similar to the amount of discordance found in other studies that assume only biological causes of variation. This similarity suggests that the role of systematic error in generating gene tree variation is underappreciated and critical evaluation of fit between assumed models and the data used for inference is important for the resolution of unresolved phylogenetic questions.

Speciation in sympatry with ongoing secondary gene flow and a potential olfactory trigger in a radiation of Cameroon cichlids.

The process of sympatric speciation in nature remains a fundamental unsolved problem. Cameroon crater lake cichlid radiations were long regarded as one of the most compelling examples; however, recent work showed that their origins were more complex than a single colonization event followed by isolation. Here, we performed a detailed investigation of the speciation history of a radiation of Coptodon cichlids from Lake Ejagham, Cameroon, using whole-genome sequencing data. The existence of the Lake Ejagham Coptodon radiation is remarkable as this 0.5 km2 lake offers limited scope for divergence across a shallow depth gradient, disruptive selection is currently weak, and the species are sexually monochromatic. We infer that Lake Ejagham was colonized by Coptodon cichlids soon after its formation 9,000 years ago, yet speciation occurred only in the last 1,000-2,000 years. We show that secondary gene flow from a nearby riverine species has been ongoing, into ancestral as well as extant lineages, and we identify and date river-to-lake admixture blocks. One block contains a cluster of olfactory receptor genes that introgressed near the time of the first speciation event and coincides with a higher overall rate of admixture. Olfactory signalling is a key component of mate choice and species recognition in cichlids. A functional role for this introgression event is consistent with previous findings that sexual isolation appears much stronger than ecological isolation in Ejagham Coptodon. We conclude that speciation in this radiation took place in sympatry, yet may have benefited from ongoing riverine gene flow.

Loci associated with cave-derived traits concentrate in specific regions of the Mexican cavefish genome.

A major goal of modern evolutionary biology is connecting phenotypic evolution with its underlying genetic basis. The Mexican cavefish (Astyanax mexicanus), a characin fish species comprised of a surface ecotype and a cave-derived ecotype, is well suited as a model to study the genetic mechanisms underlying adaptation to extreme environments. Here we map 206 previously published quantitative trait loci (QTL) for cave-derived traits in A. mexicanus to the newest version of the surface fish genome assembly, AstMex3. This analysis revealed that QTL cluster in the genome more than expected by chance, and this clustering is not explained by the distribution of genes in the genome. To investigate whether certain characteristics of the genome facilitate phenotypic evolution, we tested whether genomic characteristics, such as highly mutagenic CpG sites, are reliable predictors of the sites of trait evolution but did not find any significant trends. Finally, we combined the QTL map with previously collected expression and selection data to identify a list of 36 candidate genes that may underlie the repeated evolution of cave phenotypes, including rgrb which is predicted to be involved in phototransduction. We found this gene has disrupted exons in all non-hybrid cave populations but intact reading frames in surface fish. Overall, our results suggest specific "evolutionary hotspots" in the genome may play significant roles in driving adaptation to the cave environment in Astyanax mexicanus and demonstrate how this compiled dataset can facilitate our understanding of the genetic basis of repeated evolution in the Mexican cavefish.

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.

Selection-driven trait loss in independently evolved cavefish populations.

Laboratory studies have demonstrated that a single phenotype can be produced by many different genotypes; however, in natural systems, it is frequently found that phenotypic convergence is due to parallel genetic changes. This suggests a substantial role for constraint and determinism in evolution and indicates that certain mutations are more likely to contribute to phenotypic evolution. Here we use whole genome resequencing in the Mexican tetra, Astyanax mexicanus, to investigate how selection has shaped the repeated evolution of both trait loss and enhancement across independent cavefish lineages. We show that selection on standing genetic variation and de novo mutations both contribute substantially to repeated adaptation. Our findings provide empirical support for the hypothesis that genes with larger mutational targets are more likely to be the substrate of repeated evolution and indicate that features of the cave environment may impact the rate at which mutations occur.

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.

Marketplace shrimp mislabeling in North Carolina.

Seafood mislabeling occurs in a wide range of seafood products worldwide, resulting in public distrust, economic fraud, and health risks for consumers. We quantified the extent of shrimp mislabeling in coastal and inland North Carolina. We used standard DNA barcoding procedures to determine the species identity of 106 shrimp sold as "local" by 60 vendors across North Carolina. Thirty-four percent of the purchased shrimp was mislabeled, and surprisingly the percentage did not differ significantly between coastal and inland counties. One third of product incorrectly marketed as "local" was in fact whiteleg shrimp: an imported and globally farmed species native to the eastern Pacific, not found in North Carolina waters. In addition to the negative ecosystem consequences of shrimp farming (e.g., the loss of mangrove forests and the coastal buffering they provide), North Carolina fishers-as with local fishers elsewhere-are negatively impacted when vendors label farmed, frozen, and imported shrimp as local, fresh, and wild-caught.

A high proportion of red snapper sold in North Carolina is mislabeled.

Seafood mislabeling occurs when a market label is inaccurate, primarily in terms of species identity, but also regarding weight, geographic origin, or other characteristics. This widespread problem allows cheaper or illegally-caught species to be marketed as species desirable to consumers. Previous studies have identified red snapper (Lutjanus campechanus) as one of the most frequently mislabeled seafood species in the United States. To quantify how common mislabeling of red snapper is across North Carolina, the Seafood Forensics class at the University of North Carolina at Chapel Hill used DNA barcoding to analyze samples sold as "red snapper" from restaurants, seafood markets, and grocery stores purchased in ten counties. Of 43 samples successfully sequenced and identified, 90.7% were mislabeled. Only one grocery store chain (of four chains tested) accurately labeled red snapper. The mislabeling rate for restaurants and seafood markets was 100%. Vermilion snapper (Rhomboplites aurorubens) and tilapia (Oreochromis aureus and O. niloticus) were the species most frequently substituted for red snapper (13 of 39 mislabeled samples for both taxa, or 26 of 39 mislabeled total). This study builds on previous mislabeling research by collecting samples of a specific species in a confined geographic region, allowing local vendors and policy makers to better understand the scope of red snapper mislabeling in North Carolina. This methodology is also a model for other academic institutions to engage undergraduate researchers in mislabeling data collection, sample processing, and analysis.

Don't throw out the sympatric speciation with the crater lake water: fine-scale investigation of introgression provides equivocal support for causal role of secondary gene flow in one of the clearest examples of sympatric speciation.

Genomic data has revealed complex histories of colonization and repeated gene flow previously unrecognized in some of the most celebrated examples of sympatric speciation and radiation. However, much of the evidence for secondary gene flow into these radiations comes from summary statistics calculated from sparse genomic sampling without knowledge of which specific genomic regions introgressed. This tells us little about how gene flow potentially influenced sympatric diversification. Here, we investigated whole genomes of Barombi Mbo crater lake cichlids for fine-scale patterns of introgression with neighboring riverine cichlid populations. We found evidence of secondary gene flow into the radiation scattered across <0.24% of the genome; however, from our analyses, it is not clear if the functional diversity in these regions contributed to the ecological, sexual, and morphological diversity found in the lake. Unlike similar studies, we found no obvious candidate genes for adaptive introgression and we cannot rule out that secondary gene flow was predominantly neutral with respect to the diversification process. We also found evidence for differential assortment of ancestral polymorphisms found in riverine populations between sympatric sister species, suggesting the presence of an ancestral hybrid swarm. Although the history of gene flow and colonization is more complicated than previously assumed, the lack of compelling evidence for secondary gene flow's role in species diversification suggests that we should not yet rule out one of the most celebrated examples of sympatric speciation in nature without a more thorough investigation of the timing and functional role of each introgressed region.