Novel trophic niches drive variable progress towards ecological speciation within an adaptive radiation of pupfishes

John M, Chau D, Clair C, Chan H, Holzman R, Martin CH. A multi-peak performance landscape for scale biting in an adaptive radiation of pupfishes

The physical interactions between organisms and their environment ultimately shape diversification rates, but the contributions of biomechanics to evolutionary divergence are frequently overlooked. Here, we estimated a performance landscape for biting in an adaptive radiation of Cyprinodon pupfishes, including scale-biting and molluscivore specialists, and compared performance peaks with previous estimates of the fitness landscape in this system. We used high-speed video to film feeding strikes on gelatin cubes by scale eater, molluscivore, generalist and hybrid pupfishes and measured bite dimensions. We then measured five kinematic variables from 227 strikes using the SLEAP machine-learning model. We found a complex performance landscape with two distinct peaks best predicted gel-biting performance, corresponding to a significant non-linear interaction between peak gape and peak jaw protrusion. Only scale eaters and their hybrids were able to perform strikes within the highest performance peak, characterized by larger peak gapes and greater jaw protrusion. A performance valley separated this peak from a lower performance peak accessible to all species, characterized by smaller peak gapes and less jaw protrusion. However, most individuals exhibited substantial variation in strike kinematics and species could not be reliably distinguished by their strikes, indicating many-to-many mapping of morphology to performance. The two performance peaks observed in the lab were partially consistent with estimates of a two-peak fitness landscape measured in the wild, with the exception of the new performance peak for scale eaters. We thus reveal a new bimodal non-linear biomechanical model that connects morphology to performance to fitness in a sympatric radiation of trophic niche specialists.

Phenotypic covariation predicts diversification in an adaptive radiation of pupfishes

Phenotypic covariation among suites of traits may constrain or promote diversification both within and between species, yet few studies have empirically tested this relationship. In this study, we investigate whether phenotypic covariation of craniofacial traits is associated with diversification in an adaptive radiation of pupfishes found only on San Salvador Island, Bahamas (SSI). The radiation includes generalist, durophagous, and lepidophagous species. We compared phenotypic variation and covariation (i.e., the P matrix) between (1) allopatric populations of generalist pupfish from neighboring islands and estuaries in the Caribbean, (2) SSI pupfish allopatric lake populations with only generalist pupfish, and (3) SSI lake populations containing the full radiation in sympatry. Additionally, we examine patterns observed in the P matrices of two independent lab-reared F2 hybrid crosses of the two most morphologically distinct members of the radiation to make inferences about the underlying mechanisms contributing to the variation in craniofacial traits in SSI pupfishes. We found that the P matrix of SSI allopatric generalist populations exhibited higher levels of mean trait correlation, constraints, and integration with simultaneously lower levels of flexibility compared to allopatric generalist populations on other Caribbean islands and sympatric populations of all three species on SSI. We also document that while many craniofacial traits appear to result from additive genetic effects, variation in key traits such as head depth, maxilla length, and lower jaw length may be produced via non-additive genetic mechanisms. Ultimately, this study suggests that differences in phenotypic covariation significantly contribute to producing and maintaining organismal diversity.

Parallel evolution of integrated craniofacial traits in trophic specialist pupfishes

Populations may adapt to similar environments via parallel or non-parallel genetic changes, but the frequency of these alternative mechanisms and underlying contributing factors are still poorly understood outside model systems. We used QTL mapping to investigate the genetic basis of highly divergent craniofacial traits between the scale-eater (Cyprinodon desquamator) and molluscivore (C. brontotheroides) pupfish adapting to two different hypersaline lake environments on San Salvador Island, Bahamas. We lab-reared F2 scale-eater x molluscivore intercrosses from two different lake populations, estimated linkage maps, scanned for significant QTL for 29 skeletal and craniofacial traits, female mate preference, and sex. We compared the location of QTL between lakes to quantify parallel and non-parallel genetic changes. We detected significant QTL for six craniofacial traits in at least one lake. However, nearly all shared QTL loci were associated with a different craniofacial trait within each lake. Therefore, our estimate of parallel evolution of craniofacial genetic architecture could range from one out of six identical trait QTL (low parallelism) to five out of six integrated trait QTL (high parallelism). We suggest that pleiotropy and trait integration can affect estimates of parallel evolution, particularly within rapid radiations. We also observed increased adaptive introgression in shared QTL regions, suggesting that gene flow contributed to parallel evolution. Overall, our results suggest that the same genomic regions may contribute to parallel adaptation across integrated suites of craniofacial traits, rather than specific traits, and highlight the need for a more expansive definition of parallel evolution.

Foraging-induced craniofacial plasticity is associated with an early, robust and dynamic transcriptional response

Phenotypic plasticity is the ability of a single genotype to vary its phenotype in response to the environment. Plasticity of the skeletal system in response to mechanical input is widely studied, but the timing of its transcriptional regulation is not well understood. Here, we used the cichlid feeding apparatus to examine the transcriptional dynamics of skeletal plasticity over time. Using three closely related species that vary in their ability to remodel bone and a panel of 11 genes, including well-studied skeletal differentiation markers and newly characterized environmentally sensitive genes, we examined plasticity at one, two, four and eight weeks following the onset of alternate foraging challenges. We found that the plastic species exhibited environment-specific bursts in gene expression beginning at one week, followed by a sharp decline in levels, while the species with more limited plasticity exhibited consistently low levels of gene expression. This trend held across nearly all genes, suggesting that it is a hallmark of the larger plasticity regulatory network. We conclude that plasticity of the cichlid feeding apparatus is not the result of slowly accumulating gene expression difference over time, but rather is stimulated by early bursts of environment-specific gene expression followed by a return to homeostatic levels.

John M, Chau D, Clair C, Chan H, Holzman R, Martin CH. 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.

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.

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.

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.

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.

The Role of Photobionts as Drivers of Diversification in an Island Radiation of Lichen-Forming Fungi

Speciation in oceanic islands has attracted the interest of scientists since the 19th century. One of the most striking evolutionary phenomena that can be studied in islands is adaptive radiation, that is, when a lineage gives rise to different species by means of ecological speciation. Some of the best-known examples of adaptive radiation are charismatic organisms like the Darwin finches of the Galapagos and the cichlid fishes of the great African lakes. In these and many other examples, a segregation of the trophic niche has been shown to be an important diversification driver. Radiations are known in other groups of organisms, such as lichen-forming fungi. However, very few studies have investigated their adaptive nature, and none have focused on the trophic niche. In this study, we explore the role of the trophic niche in a putative radiation of endemic species from the Macaronesian Region, the Ramalina decipiens group. The photobiont diversity was studied by Illumina MiSeq sequencing of the ITS2 region of 197 specimens spanning the phylogenetic breadth and geographic range of the group. A total of 66 amplicon sequence variants belonging to the four main clades of the algal genus Trebouxia were found. Approximately half of the examined thalli showed algal coexistence, but in most of them, a single main photobiont amounted to more than 90% of the reads. However, there were no significant differences in photobiont identity and in the abundance of ITS2 reads across the species of the group. We conclude that a segregation of the trophic niche has not occurred in the R. decipiens radiation.

Trophic niche shifts and phenotypic trait evolution are largely decoupled in Australasian parrots

BACKGROUND

Trophic shifts from one dietary niche to another have played major roles in reshaping the evolutionary trajectories of a wide range of vertebrate groups, yet their consequences for morphological disparity and species diversity differ among groups.

METHODS

Here, we use phylogenetic comparative methods to examine whether the evolution of nectarivory and other trophic shifts have driven predictable evolutionary pathways in Australasian psittaculid parrots in terms of ecological traits such as body size, beak shape, and dispersal capacity.

RESULTS

We found no evidence for an 'early-burst' scenario of lineage or morphological diversification. The best-fitting models indicate that trait evolution in this group is characterized by abrupt phenotypic shifts (evolutionary jumps), with no sign of multiple phenotypic optima correlating with different trophic strategies. Thus, our results point to the existence of weak directional selection and suggest that lineages may be evolving randomly or slowly toward adaptive peaks they have not yet reached.

CONCLUSIONS

This study adds to a growing body of evidence indicating that the relationship between avian morphology and feeding ecology may be more complex than usually assumed and highlights the importance of adding more flexible models to the macroevolutionary toolbox.

Linking population-level and microevolutionary processes to understand speciation dynamics at the macroevolutionary scale

Although speciation dynamics have been described for several taxonomic groups in distinct geographic regions, most macroevolutionary studies still lack a detailed mechanistic view on how or why speciation rates change. To help partially fill this gap, we suggest that the interaction between the time taken by a species to geographically expand and the time populations take to evolve reproductive isolation should be considered when we are trying to understand macroevolutionary patterns. We introduce a simple conceptual index to guide our discussion on how demographic and microevolutionary processes might produce speciation dynamics at macroevolutionary scales. Our framework is developed under different scenarios: when speciation is mediated by geographical or resource-partitioning opportunities, and when diversity is limited or not. We also discuss how organismal intrinsic properties and different overall geographical settings can influence the tempo and mode of speciation. We argue that specific conditions observed at the microscale might produce a pulse in speciation rates even without a pulse in either climate or physical barriers. We also propose a hypothesis to reconcile the apparent inconsistency between speciation measured at the microscale and macroscale, and emphasize that diversification rates are better seen as an emergent property. We hope to bring the reader's attention to interesting mechanisms to be further studied, to motivate the development of new theoretical models that connect microevolution and macroevolution, and to inspire new empirical and methodological approaches to more adequately investigate speciation dynamics either using neontological or paleontological data.

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.

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.

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.

Rapid adaptive evolution of scale-eating kinematics to a novel ecological niche

The origins of novel trophic specialization, in which organisms begin to exploit resources for the first time, may be explained by shifts in behavior such as foraging preferences or feeding kinematics. One way to investigate behavioral mechanisms underlying ecological novelty is by comparing prey capture kinematics among species. We investigated the contribution of kinematics to the origins of a novel ecological niche for scale-eating within a microendemic adaptive radiation of pupfishes on San Salvador Island, Bahamas. We compared prey capture kinematics across three species of pupfish while they consumed shrimp and scales in the lab, and found that scale-eating pupfish exhibited peak gape sizes twice as large as in other species, but also attacked prey with a more obtuse angle between their lower jaw and suspensorium. We then investigated how this variation in feeding kinematics could explain scale-biting performance by measuring bite size (surface area removed) from standardized gelatin cubes. We found that a combination of larger peak gape and more obtuse lower jaw and suspensorium angles resulted in approximately 40% more surface area removed per strike, indicating that scale-eaters may reside on a performance optimum for scale biting. To test whether feeding performance could contribute to reproductive isolation between species, we also measured F1 hybrids and found that their kinematics and performance more closely resembled generalists, suggesting that F1 hybrids may have low fitness in the scale-eating niche. Ultimately, our results suggest that the evolution of strike kinematics in this radiation is an adaptation to the novel niche of scale eating.

Historical allopatry and secondary contact or primary intergradation in the Puerto Rican crested anole, Anolis cristatellus, on Vieques Island in the Caribbean

Recent work has revealed surprisingly deep mitochondrial genetic divergence in the lizard Anolis cristatellus among samples obtained from the small Caribbean island of Vieques. Here we sought to determine whether this had resulted from natural or anthropogenic causes, and (if the former) whether divergence occurred in a biogeographical context of allopatry followed by secondary contact, or via isolation-by-distance across the species’ historical range. We first estimated a mitochondrial gene tree for 379 samples and then genotyped 3407 single nucleotide polymorphic sites from 48 individuals using a modified genotyping-by-sequencing approach. We found that A. cristatellus samples from Vieques belong to two highly divergent mitochondrial subclades, but the geographical distribution of these haplogroups indicates that this pattern is probably natural in origin. Analysis of our single nucleotide polymorphic dataset revealed differentiation that is consistent with isolation-by-distance between the western and eastern ends of Vieques, suggesting that the overall pattern of divergence probably reflects primary intergradation with a mitochondrial break on the historical Puerto Rico Bank palaeo-island that happened to coincide with the present-day location of Vieques. Our findings help to underline the growing consensus that results from a single genetic marker can prove highly misleading in studies of historical population genetic structure.

John ME. 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.

Hybrid gene misregulation in multiple developing tissues within a recent adaptive radiation of Cyprinodon pupfishes

Genetic incompatibilities constitute the final stages of reproductive isolation and speciation, but little is known about incompatibilities that occur within recent adaptive radiations among closely related diverging populations. Crossing divergent species to form hybrids can break up coadapted variation, resulting in genetic incompatibilities within developmental networks shaping divergent adaptive traits. We crossed two closely related sympatric Cyprinodon pupfish species–a dietary generalist and a specialized molluscivore–and measured expression levels in their F1 hybrids to identify regulatory variation underlying the novel craniofacial morphology found in this recent microendemic adaptive radiation. We extracted mRNA from eight day old whole-larvae tissue and from craniofacial tissues dissected from 17–20 day old larvae to compare gene expression between a total of seven F1 hybrids and 24 individuals from parental species populations. We found 3.9% of genes differentially expressed between generalists and molluscivores in whole-larvae tissues and 0.6% of genes differentially expressed in craniofacial tissue. We found that 2.1% of genes were misregulated in whole-larvae hybrids whereas 19.1% of genes were misregulated in hybrid craniofacial tissues, after correcting for sequencing biases. We also measured allele specific expression across 15,429 heterozygous sites to identify putative compensatory regulatory mechanisms underlying differential expression between generalists and molluscivores. Together, our results highlight the importance of considering misregulation as an early indicator of genetic incompatibilities in the context of rapidly diverging adaptive radiations and suggests that compensatory regulatory divergence drives hybrid gene misregulation in developing tissues that give rise to novel craniofacial traits.

The behavioral origins of novelty: did increased aggression lead to scale-eating in pupfishes?

Behavioral changes in a new environment are often assumed to precede the origins of evolutionary novelties. Here, we examined whether an increase in aggression is associated with a novel scale-eating trophic niche within a recent radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas. We measured aggression using multiple behavioral assays and used transcriptomic analyses to identify differentially expressed genes in aggression and other behavioral pathways across 3 sympatric species in the San Salvador radiation (generalist, snail-eating specialist, and scale-eating specialist) and 2 generalist outgroups. Surprisingly, we found increased behavioral aggression and differential expression of aggression-related pathways in both the scale-eating and snail-eating specialists, despite their independent evolutionary origins. Increased behavioral aggression varied across both sex and stimulus context in both species. Our results indicate that aggression is not unique to scale-eating specialists. Instead, selection may increase aggression in other contexts such as niche specialization in general or mate competition. Alternatively, increased aggression may result from indirect selection on craniofacial traits, pigmentation, or metabolism-all traits which are highly divergent, exhibit signs of selective sweeps, and are affected by aggression-related genetic pathways which are differentially expressed in this system. In conclusion, the evolution of a novel predatory trophic niche within a recent adaptive radiation does not have clear-cut behavioral origins as previously assumed, highlighting the multivariate nature of adaptation and the complex integration of behavior with other phenotypic traits.

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.

Genomic data reject the hypothesis of sympatric ecological speciation in a clade of Desmognathus salamanders

Closely related taxa with dissimilar morphologies are often considered to have diverged via natural selection favoring different phenotypes. However, some studies have found these scenarios to be paired with limited or no genetic differentiation. Desmognathus quadramaculatus and D. marmoratus are sympatric salamander species thought to represent a case of ecological speciation based on distinct morphologies, but the results of previous studies have not resolved corresponding patterns of lineage divergence. Here, we use genome-wide data to test this hypothesis of ecological speciation. Population structure analyses partitioned individuals geographically, but not morphologically, into two adjacent regions of western North Carolina: Pisgah and Nantahala. Phylogenetic analyses confirmed the nominal species are nonmonophyletic and resolved deep divergence between the two geographic clusters. Model-testing overwhelmingly supported the hypothesis that lineage divergence followed geography. Finally, ecological niche modeling showed that Pisgah and Nantahala individuals occupy different climatic niches, and geographic boundaries for the two lineages correspond to differences in precipitation regimes across southern Appalachia. Overall, we reject the previous hypothesis of ecological speciation based on microhabitat partitioning. Instead, our results suggest that there are two cryptic lineages, each containing the same pair of morphotypes.

Disentangling the drivers of diversification in an imperiled group of freshwater fishes (Cyprinodontiformes: Goodeidae)

BACKGROUND

One of the most perplexing questions in evolutionary biology is why some lineages diversify into many species, and others do not. In many cases, ecological opportunity has played an important role, leading to diversification along trophic or habitat-based axes. The Goodeidae (Teleostomi: Cyprinodontiformes) are a family of freshwater fishes with two subfamilies: Goodeinae (42 species, viviparous, heterogeneous habitats, Mesa Central of Mexico) and Empetrichthyinae (4 species, oviparous, homogeneous habitats, Great Basin of the United States). These discrepant sets of characteristics and their sister-group relationship make the goodeids amenable to a comparative study of diversification. We gathered lateral body images from more than 1600 specimens of all extant species in the family. Geometric morphometric, and phylogenetic comparative analyses were used to address whether higher species diversity correlates with higher rates of morphological shape evolution and whether there are differences in functional/habitat modules between the two subfamilies.

RESULTS

This study recovered a higher rate of overall body shape evolution in the Goodeinae that is nearly double in magnitude compared to the Empetrichthyinae. A modularity test indicated that the Goodeinae displayed elevated rates of morphological evolution in comparison to the Empetrichthyinae when only trunk (locomotor) regions were compared between subfamilies. No significant differences in evolutionary shape rates were recovered when the trophic (head) regions were compared between subfamilies.

DISCUSSION

These results support the hypothesis that Mexican goodeids radiated via an ecological opportunity scenario into a wide-array of novel habitats in the island-like Mesa Central as evidenced by their high rate of shape evolution, relative to the Empetrichthyinae. This study quantitatively unraveled the drivers of evolution and eliminated trophic specialization as a driving force within the Goodeidae.

CONCLUSIONS

A combination of phylogenetic and morphometric data, and phylogenetic comparative analyses were used to examine body shape rate evolution within the Goodeidae. Results support the hypothesis that species in the subfamily Goodeinae on the central Mexican plateau had a higher rate of body shape evolution relative to its sister subfamily Empetrichthyinae in the Great Basin suggesting that the Goodeinae diversified via an ecological opportunity scenario along habitat, rather than trophic axes.

Speciation through the lens of biomechanics: locomotion, prey capture and reproductive isolation

Speciation is a multifaceted process that involves numerous aspects of the biological sciences and occurs for multiple reasons. Ecology plays a major role, including both abiotic and biotic factors. Whether populations experience similar or divergent ecological environments, they often adapt to local conditions through divergence in biomechanical traits. We investigate the role of biomechanics in speciation using fish predator-prey interactions, a primary driver of fitness for both predators and prey. We highlight specific groups of fishes, or specific species, that have been particularly valuable for understanding these dynamic interactions and offer the best opportunities for future studies that link genetic architecture to biomechanics and reproductive isolation (RI). In addition to emphasizing the key biomechanical techniques that will be instrumental, we also propose that the movement towards linking biomechanics and speciation will include (i) establishing the genetic basis of biomechanical traits, (ii) testing whether similar and divergent selection lead to biomechanical divergence, and (iii) testing whether/how biomechanical traits affect RI. Future investigations that examine speciation through the lens of biomechanics will propel our understanding of this key process.

Building trophic specializations that result in substantial niche partitioning within a young adaptive radiation

Dietary partitioning often accompanies the increased morphological diversity seen during adaptive radiations within aquatic systems. While such niche partitioning would be expected in older radiations, it is unclear how significant morphological divergence occurs within a shorter time period. Here we show how differential growth in key elements of the feeding mechanism can bring about pronounced functional differences among closely related species. An incredibly young adaptive radiation of three Cyprinodon species residing within hypersaline lakes in San Salvador Island, Bahamas, has recently been described. Characterized by distinct head shapes, gut content analyses revealed three discrete feeding modes in these species: basal detritivory as well as derived durophagy and lepidophagy (scale-feeding). We dissected, cleared and stained, and micro-CT scanned species to assess functionally relevant differences in craniofacial musculoskeletal elements. The widespread feeding mode previously described for cyprinodontiforms, in which the force of the bite may be secondary to the requisite dexterity needed to pick at food items, is modified within both the scale specialist and the durophagous species. While the scale specialist has greatly emphasized maxillary retraction, using it to overcome the poor mechanical advantage associated with scale-eating, the durophage has instead stabilized the maxilla. In all species the bulk of the adductor musculature is composed of AM A1. However, the combined masses of both adductor mandibulae (AM) A1 and A3 in the scale specialist were five times that of the other species, showing the importance of growth in functional divergence. The scale specialist combines plesiomorphic jaw mechanisms with both a hypertrophied AM A1 and a slightly modified maxillary anatomy (with substantial functional implications) to generate a bite that is both strong and allows a wide range of motion in the upper jaw, two attributes that normally tradeoff mechanically. Thus, a significant feeding innovation (scale-eating, rarely seen in fishes) may evolve based largely on allometric changes in ancestral structures. Alternatively, the durophage shows reduced growth with foreshortened jaws that are stabilized by an immobile maxilla. Overall, scale specialists showed the most divergent morphology, suggesting that selection for scale-biting might be stronger or act on a greater number of traits than selection for either detritivory or durophagy. The scale specialist has colonized an adaptive peak that few lineages have climbed. Thus, heterochronic changes in growth can quickly produce functionally relevant change among closely related species.

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.

The Cyprinodon variegatus genome reveals gene expression changes underlying differences in skull morphology among closely related species

BACKGROUND

Understanding the genetic and developmental origins of phenotypic novelty is central to the study of biological diversity. In this study we identify modifications to the expression of genes at four developmental stages that may underlie jaw morphological differences among three closely related species of pupfish (genus Cyprinodon) from San Salvador Island, Bahamas. Pupfishes on San Salvador Island are trophically differentiated and include two endemic species that have evolved jaw morphologies unlike that of any other species in the genus Cyprinodon.

RESULTS

We find that gene expression differs significantly across recently diverged species of pupfish. Genes such as Bmp4 and calmodulin, previously implicated in jaw diversification in African cichlid fishes and Galapagos finches, were not found to be differentially expressed among species of pupfish. Instead we find multiple growth factors and cytokine/chemokine genes to be differentially expressed among these pupfish taxa. These include both genes and pathways known to affect craniofacial development, such as Wnt signaling, as well as novel genes and pathways not previously implicated in craniofacial development. These data highlight both shared and potentially unique sources of jaw diversity in pupfish and those identified in other evolutionary model systems such as Galapagos finches and African cichlids.

CONCLUSIONS

We identify modifications to the expression of genes involved in Wnt signaling, Igf signaling, and the inflammation response as promising avenues for future research. Our project provides insight into the magnitude of gene expression changes contributing to the evolution of morphological novelties, such as jaw structure, in recently diverged pupfish species.

Molecular ecology studies of species radiations: current research gaps, opportunities and challenges

Understanding the drivers and limits of species radiations is a crucial goal of evolutionary genetics and molecular ecology, yet research on this topic has been hampered by the notorious difficulty of connecting micro- and macroevolutionary approaches to studying the drivers of diversification. To chart the current research gaps, opportunities and challenges of molecular ecology approaches to studying radiations, we examine the literature in the journal Molecular Ecology and revisit recent high-profile examples of evolutionary genomic research on radiations. We find that available studies of radiations are highly unevenly distributed among taxa, with many ecologically important and species-rich organismal groups remaining severely understudied, including arthropods, plants and fungi. Most studies employed molecular methods suitable over either short or long evolutionary time scales, such as microsatellites or restriction site-associated DNA sequencing (RAD-seq) in the former case and conventional amplicon sequencing of organellar DNA in the latter. The potential of molecular ecology studies to address and resolve patterns and processes around the species level in radiating groups of taxa is currently limited primarily by sample size and a dearth of information on radiating nuclear genomes as opposed to organellar ones. Based on our literature survey and personal experience, we suggest possible ways forward in the coming years. We touch on the potential and current limitations of whole-genome sequencing (WGS) in studies of radiations. We suggest that WGS and targeted ('capture') resequencing emerge as the methods of choice for scaling up the sampling of populations, species and genomes, including currently understudied organismal groups and the genes or regulatory elements expected to matter most to species radiations.

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.

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.

Niche divergence facilitated by fine-scale ecological partitioning in a recent cichlid fish adaptive radiation

Ecomorphological differentiation is a key feature of adaptive radiations, with a general trend for specialization and niche expansion following divergence. Ecological opportunity afforded by invasion of a new habitat is thought to act as an ecological release, facilitating divergence, and speciation. Here, we investigate trophic adaptive morphology and ecology of an endemic clade of oreochromine cichlid fishes (Alcolapia) that radiated along a herbivorous trophic axis following colonization of an isolated lacustrine environment, and demonstrate phenotype‐environment correlation. Ecological and morphological divergence of the Alcolapia species flock are examined in a phylogenomic context, to infer ecological niche occupation within the radiation. Species divergence is observed in both ecology and morphology, supporting the importance of ecological speciation within the radiation. Comparison with an outgroup taxon reveals large‐scale ecomorphological divergence but shallow genomic differentiation within the Alcolapia adaptive radiation. Ancestral morphological reconstruction suggests lake colonization by a generalist oreochromine phenotype that diverged in Lake Natron to varied herbivorous morphologies akin to specialist herbivores in Lakes Tanganyika and Malawi.

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.

Multispecies Outcomes of Sympatric Speciation after Admixture with the Source Population in Two Radiations of Nicaraguan Crater Lake Cichlids

The formation of species in the absence of geographic barriers (i.e. sympatric speciation) remains one of the most controversial topics in evolutionary biology. While theoretical models have shown that this most extreme case of primary divergence-with-gene-flow is possible, only a handful of accepted empirical examples exist. And even for the most convincing examples uncertainties remain; complex histories of isolation and secondary contact can make species falsely appear to have originated by sympatric speciation. This alternative scenario is notoriously difficult to rule out. Midas cichlids inhabiting small and remote crater lakes in Nicaragua are traditionally considered to be one of the best examples of sympatric speciation and lend themselves to test the different evolutionary scenarios that could lead to apparent sympatric speciation since the system is relatively small and the source populations known. Here we reconstruct the evolutionary history of two small-scale radiations of Midas cichlids inhabiting crater lakes Apoyo and Xiloá through a comprehensive genomic data set. We find no signs of differential admixture of any of the sympatric species in the respective radiations. Together with coalescent simulations of different demographic models our results support a scenario of speciation that was initiated in sympatry and does not result from secondary contact of already partly diverged populations. Furthermore, several species seem to have diverged simultaneously, making Midas cichlids an empirical example of multispecies outcomes of sympatric speciation. Importantly, however, the demographic models strongly support an admixture event from the source population into both crater lakes shortly before the onset of the radiations within the lakes. This opens the possibility that the formation of reproductive barriers involved in sympatric speciation was facilitated by genetic variants that evolved in a period of isolation between the initial founding population and the secondary migrants that came from the same source population. Thus, the exact mechanisms by which these species arose might be different from what had been thought before.

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.

High levels of interspecific gene flow in an endemic cichlid fish adaptive radiation from an extreme lake environment

Studying recent adaptive radiations in isolated insular systems avoids complicating causal events and thus may offer clearer insight into mechanisms generating biological diversity. Here, we investigate evolutionary relationships and genomic differentiation within the recent radiation of Alcolapia cichlid fish that exhibit extensive phenotypic diversification, and which are confined to the extreme soda lakes Magadi and Natron in East Africa. We generated an extensive RAD data set of 96 individuals from multiple sampling sites and found evidence for genetic admixture between species within Lake Natron, with the highest levels of admixture between sympatric populations of the most recently diverged species. Despite considerable environmental separation, populations within Lake Natron do not exhibit isolation by distance, indicating panmixia within the lake, although individuals within lineages clustered by population in phylogenomic analysis. Our results indicate exceptionally low genetic differentiation across the radiation despite considerable phenotypic trophic variation, supporting previous findings from smaller data sets; however, with the increased power of densely sampled SNPs, we identify genomic peaks of differentiation (FST outliers) between Alcolapia species. While evidence of ongoing gene flow and interspecies hybridization in certain populations suggests that Alcolapia species are incompletely reproductively isolated, the identification of outlier SNPs under diversifying selection indicates the radiation is undergoing adaptive divergence.

Complex histories of repeated gene flow in Cameroon crater lake cichlids cast doubt on one of the clearest examples of sympatric speciation

One of the most celebrated examples of sympatric speciation in nature are monophyletic radiations of cichlid fishes endemic to Cameroon crater lakes. However, phylogenetic inference of monophyly may not detect complex colonization histories involving some allopatric isolation, such as double invasions obscured by genome-wide gene flow. Population genomic approaches are better suited to test hypotheses of sympatric speciation in these cases. Here, we use comprehensive sampling from all four sympatric crater lake cichlid radiations in Cameroon and outgroups across Africa combined with next-generation sequencing to genotype tens of thousands of SNPs. We find considerable evidence of gene flow between all four radiations and neighboring riverine populations after initial colonization. In a few cases, some sympatric species are more closely related to outgroups than others, consistent with secondary gene flow facilitating their speciation. Our results do not rule out sympatric speciation in Cameroon cichlids, but rather reveal a complex history of speciation with gene flow, including allopatric and sympatric phases, resulting in both reproductively isolated species and incipient species complexes. The best remaining non-cichlid examples of sympatric speciation all involve assortative mating within microhabitats. We speculate that this feature may be necessary to complete the process of sympatric speciation in nature.

AftrRAD: a pipeline for accurate and efficient de novo assembly of RADseq data

An increase in studies using restriction site-associated DNA sequencing (RADseq) methods has led to a need for both the development and assessment of novel bioinformatic tools that aid in the generation and analysis of these data. Here, we report the availability of AftrRAD, a bioinformatic pipeline that efficiently assembles and genotypes RADseq data, and outputs these data in various formats for downstream analyses. We use simulated and experimental data sets to evaluate AftrRAD's ability to perform accurate de novo assembly of loci, and we compare its performance with two other commonly used programs, stacks and pyrad. We demonstrate that AftrRAD is able to accurately assemble loci, while accounting for indel variation among alleles, in a more computationally efficient manner than currently available programs. AftrRAD run times are not strongly affected by the number of samples in the data set, making this program a useful tool when multicore systems are not available for parallel processing, or when data sets include large numbers of samples.

Speciation in Freshwater Fishes

The extraordinary species richness of freshwater fishes has attracted much research on mechanisms and modes of speciation. We here review research on speciation in freshwater fishes in light of speciation theory, and place this in a context of broad-scale diversity patterns in freshwater fishes. We discuss several major repeated themes in freshwater fish speciation and the speciation mechanisms they are frequently associated with. These include transitions between marine and freshwater habitats, transitions between discrete freshwater habitats, and ecological transitions within habitats, as well as speciation without distinct niche shifts. Major research directions in the years to come include understanding the transition from extrinsic environment-dependent to intrinsic reproductive isolation and its influences on species persistence and understanding the extrinsic and intrinsic constraints to speciation and how these relate to broad-scale diversification patterns through time.