Museum genomics reveals the hybrid origin of an extinct crater lake endemic.

Crater lake fishes are common evolutionary model systems, with recent studies suggesting a key role for gene flow in promoting rapid adaptation and speciation. However, the study of these young lakes can be complicated by human-mediated extinctions. Museum genomics approaches integrating genetic data from recently extinct species are therefore critical to understanding the complex evolutionary histories of these fragile systems. Here, we examine the evolutionary history of an extinct Southern Hemisphere crater lake endemic, the rainbowfish Melanotaenia eachamensis. We undertook comprehensive sampling of extant rainbowfish populations of the Atherton Tablelands of Australia alongside historical museum material to understand the evolutionary origins of the extinct crater lake population and the dynamics of gene flow across the ecoregion. The extinct crater lake species is genetically distinct from all other nearby populations due to historic introgression between two proximate riverine lineages, similar to other prominent crater lake speciation systems, but this historic gene flow has not been sufficient to induce a species flock. Our results suggest that museum genomics approaches can be successfully combined with extant sampling to unravel complex speciation dynamics involving recently extinct species.

Body size shapes song in honeyeaters.

Birdsongs are among the most distinctive animal signals. Their evolution is thought to be shaped simultaneously by habitat structure and by the constraints of morphology. Habitat structure affects song transmission and detectability, thus influencing song (the acoustic adaptation hypothesis), while body size and beak size and shape necessarily constrain song characteristics (the morphological constraint hypothesis). Yet, support for the acoustic adaptation and morphological constraint hypotheses remains equivocal, and their simultaneous examination is infrequent. Using a phenotypically diverse Australasian bird clade, the honeyeaters (Aves: Meliphagidae), we compile a dataset consisting of song, environmental, and morphological variables for 163 species and jointly examine predictions of these two hypotheses. Overall, we find that body size constrains song frequency and pace in honeyeaters. Although habitat type and environmental temperature influence aspects of song, that influence is indirect, likely via effects of environmental variation on body size, with some evidence that elevation constrains the evolution of song peak frequency. Our results demonstrate that morphology has an overwhelming influence on birdsong, in support of the morphological constraint hypothesis, with the environment playing a secondary role generally via body size rather than habitat structure. These results suggest that changing body size (a consequence of both global effects such as climate change and local effects such as habitat transformation) will substantially influence the nature of birdsong.

Convergent genomic signatures associated with vertebrate viviparity.

BACKGROUND: Viviparity-live birth-is a complex and innovative mode of reproduction that has evolved repeatedly across the vertebrate Tree of Life. Viviparous species exhibit remarkable levels of reproductive diversity, both in the amount of care provided by the parent during gestation, and the ways in which that care is delivered. The genetic basis of viviparity has garnered increasing interest over recent years; however, such studies are often undertaken on small evolutionary timelines, and thus are not able to address changes occurring on a broader scale. Using whole genome data, we investigated the molecular basis of this innovation across the diversity of vertebrates to answer a long held question in evolutionary biology: is the evolution of convergent traits driven by convergent genomic changes? RESULTS: We reveal convergent changes in protein family sizes, protein-coding regions, introns, and untranslated regions (UTRs) in a number of distantly related viviparous lineages. Specifically, we identify 15 protein families showing evidence of contraction or expansion associated with viviparity. We additionally identify elevated substitution rates in both coding and noncoding sequences in several viviparous lineages. However, we did not find any convergent changes-be it at the nucleotide or protein level-common to all viviparous lineages. CONCLUSIONS: Our results highlight the value of macroevolutionary comparative genomics in determining the genomic basis of complex evolutionary transitions. While we identify a number of convergent genomic changes that may be associated with the evolution of viviparity in vertebrates, there does not appear to be a convergent molecular signature shared by all viviparous vertebrates. Ultimately, our findings indicate that a complex trait such as viviparity likely evolves with changes occurring in multiple different pathways.

Cycles of fusion and fission enabled rapid parallel adaptive radiations in African cichlids.

Although some lineages of animals and plants have made impressive adaptive radiations when provided with ecological opportunity, the propensities to radiate vary profoundly among lineages for unknown reasons. In Africa's Lake Victoria region, one cichlid lineage radiated in every lake, with the largest radiation taking place in a lake less than 16,000 years old. We show that all of its ecological guilds evolved in situ. Cycles of lineage fusion through admixture and lineage fission through speciation characterize the history of the radiation. It was jump-started when several swamp-dwelling refugial populations, each of which were of older hybrid descent, met in the newly forming lake, where they fused into a single population, resuspending old admixture variation. Each population contributed a different set of ancient alleles from which a new adaptive radiation assembled in record time, involving additional fusion-fission cycles. We argue that repeated fusion-fission cycles in the history of a lineage make adaptive radiation fast and predictable.

Geographic range size and speciation in honeyeaters.

BACKGROUND: Darwin and others proposed that a species' geographic range size positively influences speciation likelihood, with the relationship potentially dependent on the mode of speciation and other contributing factors, including geographic setting and species traits. Several alternative proposals for the influence of range size on speciation rate have also been made (e.g. negative or a unimodal relationship with speciation). To examine Darwin's proposal, we use a range of phylogenetic comparative methods, focusing on a large Australasian bird clade, the honeyeaters (Aves: Meliphagidae). RESULTS: We consider the influence of range size, shape, and position (latitudinal and longitudinal midpoints, island or continental species), and consider two traits known to influence range size: dispersal ability and body size. Applying several analytical approaches, including phylogenetic Bayesian path analysis, spatiophylogenetic models, and state-dependent speciation and extinction models, we find support for both the positive relationship between range size and speciation rate and the influence of mode of speciation. CONCLUSIONS: Honeyeater speciation rate differs considerably between islands and the continental setting across the clade's distribution, with range size contributing positively in the continental setting, while dispersal ability influences speciation regardless of setting. These outcomes support Darwin's original proposal for a positive relationship between range size and speciation likelihood, while extending the evidence for the contribution of dispersal ability to speciation.

Analysis of ancestry heterozygosity suggests that hybrid incompatibilities in threespine stickleback are environment dependent.

Hybrid incompatibilities occur when interactions between opposite ancestry alleles at different loci reduce the fitness of hybrids. Most work on incompatibilities has focused on those that are "intrinsic," meaning they affect viability and sterility in the laboratory. Theory predicts that ecological selection can also underlie hybrid incompatibilities, but tests of this hypothesis using sequence data are scarce. In this article, we compiled genetic data for F2 hybrid crosses between divergent populations of threespine stickleback fish (Gasterosteus aculeatus L.) that were born and raised in either the field (seminatural experimental ponds) or the laboratory (aquaria). Because selection against incompatibilities results in elevated ancestry heterozygosity, we tested the prediction that ancestry heterozygosity will be higher in pond-raised fish compared to those raised in aquaria. We found that ancestry heterozygosity was elevated by approximately 3% in crosses raised in ponds compared to those raised in aquaria. Additional analyses support a phenotypic basis for incompatibility and suggest that environment-specific single-locus heterozygote advantage is not the cause of selection on ancestry heterozygosity. Our study provides evidence that, in stickleback, a coarse-albeit indirect-signal of environment-dependent hybrid incompatibility is reliably detectable and suggests that extrinsic incompatibilities can evolve before intrinsic incompatibilities.

Delayed Adaptive Radiation among New Zealand Stream Fishes: Joint Estimation of Divergence Time and Trait Evolution in a Newly Delineated Island Species Flock.

Adaptive radiations are generally thought to occur soon after a lineage invades a region offering high levels of ecological opportunity. However, few adaptive radiations beyond a handful of exceptional examples are known, so a comprehensive understanding of their dynamics is still lacking. Here, we present a novel case of an island species flock of freshwater fishes with a radically different tempo of adaptive history than that found in many popular evolutionary model systems. Using a phylogenomic data set combined with simultaneous Bayesian estimation of divergence times and trait-based speciation and extinction models, we show that the New Zealand Gobiomorphus gudgeons comprise a monophyletic assemblage, but surprisingly, the radiation did not fully occupy freshwater habitats and explosively speciate until more than 10 myr after the lineage invaded the islands. This shift in speciation rate was not accompanied by an acceleration in the rate of morphological evolution in the freshwater crown clade relative to the other species, but is correlated with a reduction in head pores and scales as well as an increase in egg size. Our results challenge the notion that clades always rapidly exploit ecological opportunities in the absence of competing lineages. Instead, we demonstrate that adaptive radiation can experience a slow start before undergoing accelerated diversification and that lineage and phenotypic diversification may be uncoupled in young radiations. [Adaptive radiation; Eleotridae; freshwater; Gobiomorphus; New Zealand.].

The ecological and genomic basis of explosive adaptive radiation.

Speciation rates vary considerably among lineages, and our understanding of what drives the rapid succession of speciation events within young adaptive radiations remains incomplete1-11. The cichlid fish family provides a notable example of such variation, with many slowly speciating lineages as well as several exceptionally large and rapid radiations12. Here, by reconstructing a large phylogeny of all currently described cichlid species, we show that explosive speciation is solely concentrated in species flocks of several large young lakes. Increases in the speciation rate are associated with the absence of top predators; however, this does not sufficiently explain explosive speciation. Across lake radiations, we observe a positive relationship between the speciation rate and enrichment of large insertion or deletion polymorphisms. Assembly of 100 cichlid genomes within the most rapidly speciating cichlid radiation, which is found in Lake Victoria, reveals exceptional 'genomic potential'-hundreds of ancient haplotypes bear insertion or deletion polymorphisms, many of which are associated with specific ecologies and shared with ecologically similar species from other older radiations elsewhere in Africa. Network analysis reveals fundamentally non-treelike evolution through recombining old haplotypes, and the origins of ecological guilds are concentrated early in the radiation. Our results suggest that the combination of ecological opportunity, sexual selection and exceptional genomic potential is the key to understanding explosive adaptive radiation.

Rapid generation of ecologically relevant behavioral novelty in experimental cichlid hybrids.

The East African cichlid radiations are characterized by repeated and rapid diversification into many distinct species with different ecological specializations and by a history of hybridization events between nonsister species. Such hybridization might provide important fuel for adaptive radiation. Interspecific hybrids can have extreme trait values or novel trait combinations and such transgressive phenotypes may allow some hybrids to explore ecological niches neither of the parental species could tap into. Here, we investigate the potential of second-generation (F2) hybrids between two generalist cichlid species from Lake Malawi to exploit a resource neither parental species is specialized on: feeding by sifting sand. Some of the F2 hybrids phenotypically resembled fish of species that are specialized on sand sifting. We combined experimental behavioral and morphometric approaches to test whether the F2 hybrids are transgressive in both morphology and behavior related to sand sifting. We then performed a quantitative trait loci (QTL) analysis using RADseq markers to investigate the genetic architecture of morphological and behavioral traits. We show that transgression is present in several morphological traits, that novel trait combinations occur, and we observe transgressive trait values in sand sifting behavior in some of the F2 hybrids. Moreover, we find QTLs for morphology and for sand sifting behavior, suggesting the existence of some loci with moderate to large effects. We demonstrate that hybridization has the potential to rapidly generate novel and ecologically relevant phenotypes that may be suited to a niche neither of the parental species occupies. Interspecific hybridization may thereby contribute to the rapid generation of ecological diversity in cichlid radiations.

Rapid Divergence of Predator Functional Traits Affects Prey Composition in Aquatic Communities.

Identifying traits that underlie variation in individual performance of consumers (i.e., trait utility) can help reveal the ecological causes of population divergence and the subsequent consequences for species interactions and community structure. Here, we document a case of rapid divergence (over the past 100 generations, or ∼150 years) in foraging traits and feeding efficiency between a lake and stream population pair of threespine stickleback. Building on predictions from functional trait models of fish feeding, we analyzed foraging experiments with a Bayesian path analysis and elucidated the traits explaining variation in foraging performance and the species composition of ingested prey. Despite extensive previous research on the divergence of foraging traits among populations and ecotypes of stickleback, our results provide novel experimental evidence of trait utility for jaw protrusion, gill raker length, and gill raker spacing when foraging on a natural zooplankton assemblage. Furthermore, we discuss how these traits might contribute to the differential effects of lake and stream stickleback on their prey communities, observed in both laboratory and mesocosm conditions. More generally, our results illustrate how the rapid divergence of functional foraging traits of consumers can impact the biomass, species composition, and trophic structure of prey communities.

Reef fish functional traits evolve fastest at trophic extremes.

Trophic ecology is thought to exert a profound influence on biodiversity, but the specifics of the process are rarely examined at large spatial and evolutionary scales. We investigate how trophic position and diet breadth influence functional trait evolution in one of the most species-rich and complex vertebrate assemblages, coral reef fishes, within a large-scale phylogenetic framework. We show that, in contrast with established theory, functional traits evolve fastest in trophic specialists with narrow diet breadths at both very low and high trophic positions. Top trophic level specialists exhibit the most functional diversity, while omnivorous taxa with intermediate trophic positions and wide diet breadth have the least functional diversity. Our results reveal the importance of trophic position in shaping evolutionary dynamics while simultaneously highlighting the incredible trophic and functional diversity present in coral reef fish assemblages.

Structural Control of Nonnative Ligand Binding in Engineered Mutants of Phosphoenolpyruvate Carboxykinase.

Protein engineering to alter recognition underlying ligand binding and activity has enormous potential. Here, ligand binding for Escherichia coli phosphoenolpyruvate carboxykinase (PEPCK), which converts oxaloacetate into CO2 and phosphoenolpyruvate as the first committed step in gluconeogenesis, was engineered to accommodate alternative ligands as an exemplary system with structural information. From our identification of bicarbonate binding in the PEPCK active site at the supposed CO2 binding site, we probed binding of nonnative ligands with three oxygen atoms arranged to resemble the bicarbonate geometry. Crystal structures of PEPCK and point mutants with bound nonnative ligands thiosulfate and methanesulfonate along with strained ATP and reoriented oxaloacetate intermediates and unexpected bicarbonate were determined and analyzed. The mutations successfully altered the bound ligand position and orientation and its specificity: mutated PEPCKs bound either thiosulfate or methanesulfonate but never both. Computational calculations predicted a methanesulfonate binding mutant and revealed that release of the active site ordered solvent exerts a strong influence on ligand binding. Besides nonnative ligand binding, one mutant altered the Mn2+ coordination sphere: instead of the canonical octahedral ligand arrangement, the mutant in question had an only five-coordinate arrangement. From this work, critical features of ligand binding, position, and metal ion cofactor geometry required for all downstream events can be engineered with small numbers of mutations to provide insights into fundamental underpinnings of protein-ligand recognition. Through structural and computational knowledge, the combination of designed and random mutations aids in the robust design of predetermined changes to ligand binding and activity to engineer protein function.

Feeding ecology underlies the evolution of cichlid jaw mobility.

The fish feeding apparatus is among the most diverse functional systems in vertebrates. While morphological and mechanical variations of feeding systems are well studied, we know far less about the diversity of the motions that they produce. We explored patterns of feeding movements in African cichlids from Lakes Malawi and Tanganyika, asking whether the degree of kinesis is associated with dietary habits of species. We used geometric morphometrics to measure feeding kinesis as trajectories of shape change, based on 326 high-speed videos in 56 species. Cranial morphology was significantly related to feeding movements, both of which were distributed along a dietary axis associated with prey evasiveness. Small-mouthed cichlids that feed by scraping algae and detritus from rocks had low kinesis strikes, while large-mouthed species that eat large, evasive prey (fishes and shrimps) generated the greatest kinesis. Despite having higher overall kinesis, comparisons of trajectory shape (linearity) revealed that cichlids that eat mobile prey also displayed more kinematically conserved, or efficient, feeding motions. Our work indicates that prey evasiveness is strongly related to the evolution of cichlid jaw mobility, suggesting that this same relationship may explain the origins and diversity of highly kinetic jaws that characterize the super-radiation of spiny-rayed fishes.

Replicated divergence in cichlid radiations mirrors a major vertebrate innovation.

Decoupling of the upper jaw bones--jaw kinesis--is a distinctive feature of the ray-finned fishes, but it is not clear how the innovation is related to the extraordinary diversity of feeding behaviours and feeding ecology in this group. We address this issue in a lineage of ray-finned fishes that is well known for its ecological and functional diversity--African rift lake cichlids. We sequenced ultraconserved elements to generate a phylogenomic tree of the Lake Tanganyika and Lake Malawi cichlid radiations. We filmed a diverse array of over 50 cichlid species capturing live prey and quantified the extent of jaw kinesis in the premaxillary and maxillary bones. Our combination of phylogenomic and kinematic data reveals a strong association between biting modes of feeding and reduced jaw kinesis, suggesting that the contrasting demands of biting and suction feeding have strongly influenced cranial evolution in both cichlid radiations.

Body ram, not suction, is the primary axis of suction-feeding diversity in spiny-rayed fishes.

Suction-feeding fishes exhibit diverse prey-capture strategies that vary in their relative use of suction and predator approach (ram), which is often referred to as the ram-suction continuum. Previous research has found that ram varies more than suction distance among species, such that ram accounts for most differences in prey-capture behaviors. To determine whether these findings hold at broad evolutionary scales, we collected high-speed videos of 40 species of spiny-rayed fishes (Acanthomorpha) feeding on live prey. For each strike, we calculated the contributions of suction, body ram (swimming) and jaw ram (mouth movement relative to the body) to closing the distance between predator and prey. We confirm that the contribution of suction distance is limited even in this phylogenetically and ecologically broad sample of species, with the extreme suction area of prey-capture space conspicuously unoccupied. Instead of a continuum from suction to ram, we find that variation in body ram is the major factor underlying the diversity of prey-capture strategies among suction-feeding fishes. Independent measurement of the contribution of jaw ram revealed that it is an important component of diversity among spiny-rayed fishes, with a number of ecomorphologies relying heavily on jaw ram, including pivot feeding in syngnathiforms, extreme jaw protruders and benthic sit-and-wait ambush predators. A combination of morphological and behavioral innovations has allowed fish to invade the extreme jaw ram area of prey-capture space. We caution that while two-species comparisons may support a ram-suction trade-off, these patterns do not speak to broader patterns across spiny-rayed fishes.

Evaluating genomic divergence and parallelism in replicate ecomorphs from young and old cichlid adaptive radiations.

Comparative genomic studies of closely related species typically focus on single species pairs at one given stage of divergence. That makes it difficult to infer the continuum of evolutionary process during speciation and beyond. Here, we use whole-genome resequencing to examine genomic patterns of divergence in three sympatric cichlid species pairs with very similar functional and ecological differentiation, but different ages. We find a strong signature of increasing genomic divergence with time in both the mitochondrial genome and the nuclear genome. In contrast to many other systems, we find that in these cichlids, regions of elevated relative differentiation also exhibit increased absolute differentiation. We detect a signature of convergent evolution in a comparison of outlier regions across all three species pair comparisons, but the extent of it is modest, and regions that are strongly divergent in any one pair tend to be only slightly elevated in the other pairs, consistent with a repeatable but polygenic basis of traits that characterize the ecomorphs. Our results suggest that strong functional phenotypic differentiation, as seen in all three species pairs, is generally associated with a clear signature of genomic divergence, even in the youngest species pair.

Intermediate Kinematics Produce Inferior Feeding Performance in a Classic Case of Natural Hybridization.

Selection on naturally occurring hybrid individuals is a key component of speciation theory, but few studies examine the functional basis of hybrid performance. We examine the functional consequences of hybridization in nature, using the freshwater sunfishes (Centrarchidae), where natural hybrids have been studied for more than a century and a half. We examined bluegill (Lepomis macrochirus), green sunfish (Lepomis cyanellus), and their naturally occurring hybrid, using prey-capture kinematics and morphology to parameterize suction-feeding simulations on divergent parental resources. Hybrid individuals exhibited kinematics intermediate between those of the two parental species. However, performance assays indicated that hybrids display performance most similar to the worse-performing species for a given parental resource. Our results show that intermediate hybrid phenotypes can be impaired by a less-than-intermediate performance and hence suffer a larger loss in fitness than could be inferred from morphology alone.

A pharyngeal jaw evolutionary innovation facilitated extinction in Lake Victoria cichlids.

Evolutionary innovations, traits that give species access to previously unoccupied niches, may promote speciation and adaptive radiation. Here, we show that such innovations can also result in competitive inferiority and extinction. We present evidence that the modified pharyngeal jaws of cichlid fishes and several marine fish lineages, a classic example of evolutionary innovation, are not universally beneficial. A large-scale analysis of dietary evolution across marine fish lineages reveals that the innovation compromises access to energy-rich predator niches. We show that this competitive inferiority shaped the adaptive radiation of cichlids in Lake Tanganyika and played a pivotal and previously unrecognized role in the mass extinction of cichlid fishes in Lake Victoria after Nile perch invasion.

Origins, Innovations, and Diversification of Suction Feeding in Vertebrates.

We review the origins, prominent innovations, and major patterns of diversification in suction feeding by vertebrates. Non-vertebrate chordates and larval lamprey suspension-feed by capturing small particles in pharyngeal mucous. In most of these lineages the gentle flows that transport particles are generated by buccal cilia, although larval lamprey and thaliacean urochordates have independently evolved a weak buccal pump to generate an oscillating flow of water that is powered by elastic recovery of the pharynx following compression by buccal muscles. The evolution of jaws and the hyoid facilitated powerful buccal expansion and high-performance suction feeding as found today throughout aquatic vertebrates. We highlight three major innovations in suction feeding. Most vertebrate suction feeders have mechanisms that occlude the corners of the open mouth during feeding. This produces a planar opening that is often nearly circular in shape. Both features contribute to efficient flow of water into the mouth and help direct the flow to the area directly in front of the mouth's aperture. Among several functions that have been identified for protrusion of the upper jaw, is an increase in the hydrodynamic forces that suction feeders exert on their prey. Protrusion of the upper jaw has evolved five times in ray-finned fishes, including in two of the most successful teleost radiations, cypriniforms and acanthomorphs, and is found in about 60% of living teleost species. Diversification of the mechanisms of suction feeding and of feeding behavior reveals that suction feeders with high capacity for suction rarely approach their prey rapidly, while slender-bodied predators with low capacity for suction show the full range of attack speeds. We hypothesize that a dominant axis of diversification among suction feeders involves a trade-off between the forces that are exerted on prey and the volume of water that is ingested.