Trophic novelty is linked to exceptional rates of morphological diversification in two adaptive radiations of Cyprinodon pupfish

Evolutionary lability of a key innovation spurs rapid diversification

Rates of lineage diversification vary considerably across the tree of life, often as a result of evolutionary innovations1-5. Although the ability to produce new traits can vary between clades and may drive ecological transitions6-9, the impact of differences in the pace at which innovations evolve at macroevolutionary scales has been overlooked. Complex teeth are one innovation that contributed to the evolutionary success of major vertebrate lineages10-12. Here we show that evolutionary lability of tooth complexity, but not complexity itself, spurs rapid diversification across ray-finned fishes. Speciation rates are five times higher when transitions between simple and complex teeth occur rapidly. We find that African cichlids are unique among all fishes; they are dominated by lineages that transition between simple and complex teeth at unparalleled rates. This innovation interacted with the ecological versatility of complex teeth to spur rapid adaptive radiations in lakes Malawi, Victoria and Barombi Mbo. The marked effect on diversification stems from the tight association of tooth complexity with microhabitat and diet. Our results show that phylogenetic variation in how innovations evolve can have a stronger effect on patterns of diversification than the innovation itself. Investigating the impact of innovations from this new perspective will probably implicate more traits in causing heterogeneous diversification rates across the tree of life.

Why Do Some Lineages Radiate While Others Do Not? Perspectives for Future Research on Adaptive Radiations

Understanding the processes that drive phenotypic diversification and underpin speciation is key to elucidating how biodiversity has evolved. Although these processes have been studied across a wide array of clades, adaptive radiations (ARs), which are systems with multiple closely related species and broad phenotypic diversity, have been particularly fruitful for teasing apart the factors that drive and constrain diversification. As such, ARs have become popular candidate study systems for determining the extent to which ecological features, including aspects of organisms and the environment, and inter- and intraspecific interactions, led to evolutionary diversification. Despite substantial past empirical and theoretical work, understanding mechanistically how ARs evolve remains a major challenge. Here, we highlight a number of understudied components of the environment and of lineages themselves, which may help further our understanding of speciation and AR. We also outline some substantial remaining challenges to achieving a detailed understanding of adaptation, speciation, and the role of ecology in these processes. These major challenges include identifying factors that have a causative impact in promoting or constraining ARs, gaining a more holistic understanding of features of organisms and their environment that interact resulting in adaptation and speciation, and understanding whether the role of these organismal and environmental features varies throughout the radiation process. We conclude by providing perspectives on how future investigations into the AR process can overcome these challenges, allowing us to glean mechanistic insights into adaptation and speciation.

Preying on cyprinid snout warts (pearl organs) as a novel and peculiar habit in the Lake Malawi cichlid Docimodus evelynae

Cichlid fishes in the African Great Lakes have undergone explosive speciation, acquiring markedly varying ecologies and diets. There are multiple lineages of scale-eating cichlids, and their natural history and evolutionary ecology is only partially understood. We examined the feeding habit of Docimodus evelynae, a known scale eater, in Lake Malawi. The stomach contents of young individuals mainly consisted of unknown 1 mm hard, white warts (> 30%). To clarify the origin of these warts, we conducted an X-ray fluorometer analysis, and found they were rich in sulphur but low in silicon and calcium, suggesting they were epidermal tissues. Histological and morphological analyses revealed they were multicellular and cup-shaped. These characteristics matched only those of the pearl organs of the coexisting cyprinid Labeo cylindricus. DNA was extracted from the warts found in the stomach of five D. evelynae individuals, followed by PCR using primers targeting the partial COI gene of L. cylindricus. The resulting sequences exhibited 98% similarity to those of L. cylindricus. Pearl organs, never reported as a primary food for fish, could offer a substantial nutritional source based on calorific calculations. Understanding how this peculiar diet is foraged is essential for full comprehension of the food-web structure in this lake.

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.

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

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

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.

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

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

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.

Functional Trade-offs Asymmetrically Promote Phenotypic Evolution

Trade-offs are thought to bias evolution and are core features of many anatomical systems. Therefore, trade-offs may have far-reaching macroevolutionary consequences, including patterns of morphological, functional, and ecological diversity. Jaws, like many complex anatomical systems, are comprised of elements involved in biomechanical trade-offs. We test the impact of a core mechanical trade-off, transmission of velocity versus force (i.e., mechanical advantage), on rates of jaw evolution in Neotropical cichlids. Across 130 species representing a wide array of feeding ecologies, we find that the velocity-force trade-off impacts evolution of the surrounding jaw system. Specifically, rates of jaw evolution are faster at functional extremes than in more functionally intermediate or unspecialized jaws. Yet, surprisingly, the effect on jaw evolution is uneven across the extremes of the velocity-force continuum. Rates of jaw evolution are 4 to 10-fold faster in velocity-modified jaws, whereas force-modified jaws are 7 to 18-fold faster, compared to unspecialized jaws, depending on the extent of specialization. Further, we find that a more extreme mechanical trade-off resulted in faster rates of jaw evolution. The velocity-force trade-off reflects a gradient from specialization on capture-intensive (e.g., evasive or buried) to processing-intensive prey (e.g., attached or shelled), respectively. The velocity extreme of the trade-off is characterized by large magnitudes of trait change leading to functionally divergent specialists and ecological stasis. By contrast, the force extreme of the trade-off is characterized by enhanced ecological lability made possible by phenotypes more readily co-opted for different feeding ecologies. This asymmetry of macroevolutionary outcomes along each extreme is likely the result of an enhanced utility of the pharyngeal jaw system as force-modified oral jaws are adapted for prey that require intensive processing (e.g., algae, detritus, and molluscs). The velocity-force trade-off, a fundamental feature of many anatomical systems, promotes rapid phenotypic evolution of the surrounding jaw system in a canonical continental adaptive radiation. Considering that the velocity-force trade-off is an inherent feature of all jaw systems that involve a lower element that rotates at a joint, spanning the vast majority of vertebrates, our results may be widely applicable across the tree of life. [adaptive radiation; constraint; decoupling; jaws; macroevolution; specialization].

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.

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

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

The enrichment paradox in adaptive radiations: Emergence of predators hinders diversification in resource rich environments

Adaptive radiations are known for rapid niche diversification in response to ecological opportunity. While most resources usually exist prior to adaptive radiation, novel niches associated with novel resources can be created as a clade diversifies. For example, in African lake cichlid radiations some species prey upon other species of the clade (intraclade consumers). Using a trait-based eco-evolutionary model, we investigate the evolution of intraclade consumers in adaptive radiations and the effect of this novel trophic interaction on the diversification process of the radiating clade. We find that the evolutionary emergence of intraclade consumers halts the diversification processes of other ecomorphs as a result of increased top-down control of density. Because high productivity enables earlier evolution of intraclade consumers, highly productive environments come to harbour less species-rich radiations than comparable radiations in less productive environments. Our results reveal how macroevolutionary and community patterns can emerge from ecological and microevolutionary processes.

A phylogeny of the genus Limia (Teleostei: Poeciliidae) suggests a single-lake radiation nested in a Caribbean-wide allopatric speciation scenario

OBJECTIVE

The Caribbean is an important global biodiversity hotspot. Adaptive radiations there lead to many speciation events within a limited period and hence are particularly prominent biodiversity generators. A prime example are freshwater fish of the genus Limia, endemic to the Greater Antilles. Within Hispaniola, nine species have been described from a single isolated site, Lake Miragoâne, pointing towards extraordinary sympatric speciation. This study examines the evolutionary history of the Limia species in Lake Miragoâne, relative to their congeners throughout the Caribbean.

RESULTS

For 12 Limia species, we obtained almost complete sequences of the mitochondrial cytochrome b gene, a well-established marker for lower-level taxonomic relationships. We included sequences of six further Limia species from GenBank (total N  = 18 species). Our phylogenies are in concordance with other published phylogenies of Limia. There is strong support that the species found in Lake Miragoâne in Haiti are monophyletic, confirming a recent local radiation. Within Lake Miragoâne, speciation is likely extremely recent, leading to incomplete lineage sorting in the mtDNA. Future studies using multiple unlinked genetic markers are needed to disentangle the relationships within the Lake Miragoâne clade.

Pervasive admixture and the spread of a large-lipped form in a cichlid fish radiation

Adaptive radiations have proven important for understanding the mechanisms and processes underlying biological diversity. The convergence of form and function, as well as admixture and adaptive introgression, are common in adaptive radiations. However, distinguishing between these two scenarios remains a challenge for evolutionary research. The Midas cichlid species complex (Amphilophus spp.) is a prime example of adaptive radiation, with phenotypic diversification occurring at various stages of genetic differentiation. One species, A. labiatus, has large fleshy lips, is associated with rocky lake substrates, and occurs patchily within Lakes Nicaragua and Managua. By contrast, the similar, but thin-lipped, congener, A. citrinellus, is more common and widespread. We investigated the evolutionary history of the large-lipped form, specifically regarding whether the trait has evolved independently in both lakes from ancestral thin-lipped populations, or via dispersal and/or admixture events. We collected samples from distinct locations in both lakes, and assessed differences in morphology and ecology. Using RAD-seq, we genotyped thousands of SNPs to measure population structure and divergence, demographic history, and admixture. We found significant between-species differences in ecology and morphology, local intraspecific differences in body shape and trophic traits, but only limited intraspecific variation in lip shape. Despite clear ecological differences, our genomic approach uncovered pervasive admixture between the species and low genomic differentiation, with species within lakes being genetically more similar than species between lakes. Taken together, our results suggest a single origin of large-lips, followed by pervasive admixture and adaptive introgression, with morphology being driven by local ecological opportunities, despite ongoing gene-flow.

Variability and the primacy of the genotype

Cave animals and species flocks exhibit common evolutionary principles. In caves, all traits dependent on the information derived from light lose their biological function. Mutations destructive for such traits, but neutral for the organism as a whole, can persist and accumulate until a trait has vanished. Adaptive radiations start in ecosystems containing open niches. Here, selection on niche-specific traits, such as the viscerocranium in fish, is relaxed owing to the absence of competing species, and viscerocranial variability arises. It is transitorily high in recent and phylogenetically younger flocks, providing new phenotypes. It lessens and is completely lost after directional selection promotes the fixation of phenotypes that are best adapted. In cave animals and species flocks, single traits manifest phenotypic variability owing to relaxed selection. Like the eye in cave species, the viscerocranium can be classified a module, the development of which is encoded in gene regulatory networks. Mutations in these genes can result in new phenotypes. Regarding functionality, these mutations might be destructive and eliminated by selection, neutral and thus persisting, or beneficial and promoted to fixation by directional selection. Given the ancient heritage of teleostean fish, these gene regulatory networks might be prone to mutations at the same loci or to developmental reactions resulting in similar phenotypes in closely related or taxonomically and geographically distant species.

Differences in Cell Proliferation and Craniofacial Phenotype of Closely Related Species in the Pupfish Genus Cyprinodon

Understanding the genetic basis for phenotypic differences is fundamental to the study of macroevolutionary patterns of biological diversity. While technological advances in DNA sequencing have made researching genetic variation in wild taxa routine, fully understanding how these variants affect phenotype requires taking the next step to investigate how genetic changes alter cell and tissue interactions that ultimately produce phenotypes. In this article, we investigate a role for cell proliferation as a developmental source of craniofacial diversity in a radiation of 3 species of Cyprinodon from San Salvador Island, Bahamas. Patterns of cell proliferation in the heads of hatching-age fish differ among species of Cyprinodon, and correlate with differences in allometric growth rate among the jaws of 3 distinct species. Regional patterns of cell proliferation in the head are complex, resulting in an unintuitive result in which lower levels of cell proliferation in the posterior head region are associated with the development of relatively larger jaws in one species. We combine these data with previously published morphological and genomic data to show how studying the mechanisms generating phenotype at the cellular and tissue levels of biological organization can help mechanistically link genomic studies with classic morphological studies.

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

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

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

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

Ecomorphological divergence and trophic resource partitioning in 15 syntopic Indo-Pacific parrotfishes (Labridae: Scarini)

Adaptive diversification is a product of both phylogenetic constraint and ecological opportunity. The species-rich parrotfish genera Scarus and Chlorurus display considerable variation in trophic cranial morphology, but these parrotfishes are often described as generalist herbivores. Recent work has suggested that parrotfish partition trophic resources at very fine spatial scales, raising the question of whether interspecific differences in cranial morphology reflect trophic partitioning. We tested this hypothesis by comparing targeted feeding substrata with a previously published dataset of nine cranial morphological traits. We sampled feeding substrata of 15 parrotfish species at Lizard Island, Great Barrier Reef, Australia, by following individuals until focused biting was observed, then extracting a bite core 22 mm in diameter. Three indices were parameterized for each bite core: substratum taphonomy, maximum turf height and cover of crustose coralline algae. Parrotfish species were spread along a single axis of variation in feeding substrata: successional status of the substratum taphonomy and epilithic and endolithic biota. This axis of trophic variation was significantly correlated with cranial morphology, indicating that morphological disparity within this clade is associated with interspecific partitioning of feeding substrata. Phylogenetic signal and phylomorphospace analyses revealed that the evolution of this clade involved a hitherto-unrecognized level of trophic diversification.

Insectivory leads to functional convergence in a group of Neotropical rodents

The mandible of vertebrates serves as insertion area for masticatory muscles that originate on the skull, and its functional properties are subject to selective forces related to trophic ecology. The efficiency of masticatory muscles can be measured as mechanical advantage on the mandible, which, in turn, has the property of correlating with bite force and shape. In the present work, we quantify the mechanical advantage of the mandible of akodontine rodents, which present a diverse radiation of insectivorous specialists, to assess their relationship to the estimated bite force and diet. We also tested the degree of morphofunctional convergence in response to insectivory on the group. We found the mechanical advantages to be convergent on insectivorous species, and associated with the estimated bite force, with higher mechanical advantages in species with a stronger bite and short, robust mandibles and lower mechanical advantages in insectivorous species with weaker bites and more elongated, dorso-ventrally compressed mandibles. Insectivorous species of Akodontini are functional specialists for the consumption of live prey and may exploit the resources that shrews, moles and hedgehogs consume elsewhere.

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.

Macroevolutionary integration of phenotypes within and across ant worker castes

Phenotypic traits are often integrated into evolutionary modules: sets of organismal parts that evolve together. In social insect colonies, the concepts of integration and modularity apply to sets of traits both within and among functionally and phenotypically differentiated castes. On macroevolutionary timescales, patterns of integration and modularity within and across castes can be clues to the selective and ecological factors shaping their evolution and diversification. We develop a set of hypotheses describing contrasting patterns of worker integration and apply this framework in a broad (246 species) comparative analysis of major and minor worker evolution in the hyperdiverse ant genus Pheidole. Using geometric morphometrics in a phylogenetic framework, we inferred fast and tightly integrated evolution of mesosoma shape between major and minor workers, but slower and more independent evolution of head shape between the two worker castes. Thus, Pheidole workers are evolving as a mixture of intracaste and intercaste integration and rate heterogeneity. The decoupling of homologous traits across worker castes may represent an important process facilitating the rise of social complexity.

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.

Complex Ecological Phenotypes on Phylogenetic Trees: A Markov Process Model for Comparative Analysis of Multivariate Count Data

The evolutionary dynamics of complex ecological traits—including multistate representations of diet, habitat, and behavior—remain poorly understood. Reconstructing the tempo, mode, and historical sequence of transitions involving such traits poses many challenges for comparative biologists, owing to their multidimensional nature. Continuous-time Markov chains are commonly used to model ecological niche evolution on phylogenetic trees but are limited by the assumption that taxa are monomorphic and that states are univariate categorical variables. A necessary first step in the analysis of many complex traits is therefore to categorize species into a predetermined number of univariate ecological states, but this procedure can lead to distortion and loss of information. This approach also confounds interpretation of state assignments with effects of sampling variation because it does not directly incorporate empirical observations for individual species into the statistical inference model. In this study, we develop a Dirichlet-multinomial framework to model resource use evolution on phylogenetic trees. Our approach is expressly designed to model ecological traits that are multidimensional and to account for uncertainty in state assignments of terminal taxa arising from effects of sampling variation. The method uses multivariate count data across a set of discrete resource categories sampled for individual species to simultaneously infer the number of ecological states, the proportional utilization of different resources by different states, and the phylogenetic distribution of ecological states among living species and their ancestors. The method is general and may be applied to any data expressible as a set of observational counts from different categories. [Comparative methods; Dirichlet multinomial; ecological niche evolution; macroevolution; Markov model.]

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.

Macroevolutionary analyses indicate that repeated adaptive shifts towards predatory diets affect functional diversity in Neotropical cichlids

During adaptive radiation, diversification within clades is limited by adaptation to the available ecological niches, and this may drive patterns of both trait and species diversity. However, adaptation to disparate niches may result in varied impacts on the timing, pattern and rate of morphological evolution. In this study, we examined the relationship between feeding ecology and functional diversification across a diverse clade of freshwater fishes, the Neotropical cichlids. Species dietary niches were ordinated via multivariate analysis of stomach content data. We investigated changes in the rate and pattern of morphological diversification associated with feeding, including dietary niche and degree of dietary specialization. A major division in dietary niche space was observed between predators that consume fish and macroinvertebrates vs. other groups with diets dominated by small invertebrates, detritus or vegetation. These trophic niches were strongly associated with groupings defined by functional morphospace. Clades within the piscivore/macroinvertivore group rarely transitioned to other dietary niches. Comparatively, high dietary specialization enhanced functional diversification, driving the evolution of more extreme morphologies. Divergent patterns of trophic diversification among Neotropical cichlids appear to derive from different performance demands in regional abiotic and biotic environments associated with biogeographical history.

The Evolutionary Dynamics of Mechanically Complex Systems

Animals use a diverse array of motion to feed, escape predators, and reproduce. Linking morphology, performance, and fitness is a foundational paradigm in organismal biology and evolution. Yet, the influence of mechanical relationships on evolutionary diversity remains unresolved. Here, I focus on the many-to-one mapping of form to function, a widespread, emergent property of many mechanical systems in nature, and discuss how mechanical redundancy influences the tempo and mode of phenotypic evolution. By supplying many possible morphological pathways for functional adaptation, many-to-one mapping can release morphology from selection on performance. Consequently, many-to-one mapping decouples morphological and functional diversification. In fish, for example, parallel morphological evolution is weaker for traits that contribute to mechanically redundant motions, like suction feeding performance, than for systems with one-to-one form-function relationships, like lower jaw lever ratios. As mechanical complexity increases, historical factors play a stronger role in shaping evolutionary trajectories. Many-to-one mapping, however, does not always result in equal freedom of morphological evolution. The kinematics of complex systems can often be reduced to variation in a few traits of high mechanical effect. In various different four-bar linkage systems, for example, mechanical output (kinematic transmission) is highly sensitive to size variation in one or two links, and insensitive to variation in the others. In four-bar linkage systems, faster rates of evolution are biased to traits of high mechanical effect. Mechanical sensitivity also results in stronger parallel evolution-evolutionary transitions in mechanical output are coupled with transition in linkages of high mechanical effect. In other words, the evolutionary dynamics of complex systems can actually approximate that of simpler, one-to-one systems when mechanical sensitivity is strong. When examined in a macroevolutionary framework, the same mechanical system may experience distinct selective pressures in different groups of organisms. For example, performance tradeoffs are stronger for organisms that use the same mechanical structure for more functions. In general, stronger performance tradeoffs result in less phenotypic diversity in the system and, sometimes, a slower rate of evolution. These macroevolutionary trends can contribute to unevenness in functional and lineage diversity across the tree of life. Finally, I discuss how the evolution of mechanical systems informs our understanding of the relative roles of determinism and contingency in evolution.

A Morphological Integration Perspective on the Evolution of Dimorphism among Sexes and Social Insect Castes

Many species have evolved alternate phenotypes, thus enabling individuals to conditionally produce phenotypes that are favorable for reproductive success. Examples of this phenomenon include sexual dimorphism, alternative reproductive strategies, and social insect castes. While the evolutionary functions and developmental mechanisms of dimorphic phenotypes have been studied extensively, little attention has focused on the evolutionary covariance between each phenotype. We extend the conceptual framework and methods of morphological integration to hypothesize that dimorphic traits tend to be less integrated between sexes or social castes. In the case of social insects, we describe results from our recent study of an ant genus in which workers have major and minor worker castes that perform different behavioral repertoires in and around the nest. In the case of birds, we describe a new analysis of a family of songbirds that exhibits plumage coloration that can differ greatly between males and females, with apparently independent changes in each sex. Ant head shape, which is highly specialized in each worker caste, was weakly integrated between worker castes, whereas thorax shape, which is more monomorphic, was tightly integrated. Similarly, in birds, we found a negative association between dimorphism and the degree of integration between sexes. We also found that integration decreased in fairy wrens (Malurus) for many feather patches that evolved greater dichromatism. Together, this suggests that the process of evolving increased dimorphism results in a decrease in integration between sexes and social castes. We speculate that once a mechanism for dimorphism evolves, that mechanism can create independent variation in one sex or caste upon which selection may act.

Phylogenetic Comparative Methods and the Evolution of Multivariate Phenotypes

Evolutionary biology is multivariate, and advances in phylogenetic comparative methods for multivariate phenotypes have surged to accommodate this fact. Evolutionary trends in multivariate phenotypes are derived from distances and directions between species in a multivariate phenotype space. For these patterns to be interpretable, phenotypes should be characterized by traits in commensurate units and scale. Visualizing such trends, as is achieved with phylomorphospaces, should continue to play a prominent role in macroevolutionary analyses. Evaluating phylogenetic generalized least squares (PGLS) models (e.g., phylogenetic analysis of variance and regression) is valuable, but using parametric procedures is limited to only a few phenotypic variables. In contrast, nonparametric, permutation-based PGLS methods provide a flexible alternative and are thus preferred for high-dimensional multivariate phenotypes. Permutation-based methods for evaluating covariation within multivariate phenotypes are also well established and can test evolutionary trends in phenotypic integration. However, comparing evolutionary rates and modes in multivariate phenotypes remains an important area of future development.

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

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

Evolution of a soft-tissue foraging adaptation in African cichlids: Roles for novelty, convergence, and constraint

Understanding the origins of biodiversity demands consideration of both extrinsic (e.g., ecological opportunity) and intrinsic (e.g., developmental constraint) factors. Here, we use a combination of phylogenetic and genetic tools to address the origin of novelty in African cichlids. In particular, we focus on an extreme hypertrophied snout that is structurally integrated with the upper jaw. We show that this bizarre trait has evolved independently in at least two distinct and ecologically successful cichlid clades. We find that snout dimensions are decoupled both phenotypically and genetically, which has enabled it to evolve independently in multiple directions. Further, patterns of variation among species and within a genetic mapping pedigree suggest that relative to snout length, depth is under greater genetic and/or developmental constraint. Models of evolution suggest that snout shape is under selection for feeding behavior, with snout depth being important for algae scraping and snout length for sand sifting. Indeed, the deep snout of some algivores is achieved via an expansion of the intermaxillary ligament, which is important for jaw stability and may increase feeding performance. Overall, our data imply that the evolution of exaggerated snout depth required overcoming a genetic/developmental constraint, which led to expanded ecological opportunity via foraging adaptation.

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.

Simultaneous detection of macroevolutionary patterns in phenotypic means and rate of change with and within phylogenetic trees including extinct species

Recognizing evolutionary trends in phenotypic means and rates requires the application of phylogenetic comparative methods (PCMs). Most PCMs are unsuited to make full use of fossil information, which is a drawback, given the inclusion of such data improves, and in some cases even corrects, the proper understanding of trait evolution. Here we present a new computer application, written in R, that allows the simultaneous computation of temporal trends in phenotypic mean and evolutionary rate along a phylogeny, and to contrast such patterns among different clades within the tree. By using simulation experiments, we show the new implementation, names search.trend is as powerful as existing PCM tools in discerning macroevolutionary patterns in phenotypic means and rates, but differently from any other PCM allows comparing individual clades to each other, and provides rich information about trait evolution for all lineages in the tree.

Extreme and rapid bursts of functional adaptations shape bite force in amniotes

Adaptation is the fundamental driver of functional and biomechanical evolution. Accordingly, the states of biomechanical traits (absolute or relative trait values) have long been used as proxies for adaptations in response to direct selection. However, ignoring evolutionary history, in particular ancestry, passage of time and the rate of evolution, can be misleading. Here, we apply a recently developed phylogenetic statistical approach using significant rate shifts to detect instances of exceptional rates of adaptive changes in bite force in a large group of terrestrial vertebrates, the amniotes. Our results show that bite force in amniotes evolved through multiple bursts of exceptional rates of adaptive changes, whereby whole groups-including Darwin's finches, maniraptoran dinosaurs (group of non-avian dinosaurs including birds), anthropoids and hominins (fossil and modern humans)-experienced significant rate increases compared to the background rate. However, in most parts of the amniote tree of life, we find no exceptional rate increases, indicating that coevolution with body size was primarily responsible for the patterns observed in bite force. Our approach represents a template for future studies in functional morphology and biomechanics, where exceptional rates of adaptive changes can be quantified and potentially linked to specific ecological factors underpinning major evolutionary radiations.

Correlates of rate heterogeneity in avian ecomorphological traits

Heterogeneity in rates of trait evolution is widespread, but it remains unclear which processes drive fast and slow character divergence across global radiations. Here, we test multiple hypotheses for explaining rate variation in an ecomorphological trait (beak shape) across a globally distributed group (birds). We find low support that variation in evolutionary rates of species is correlated with life history, environmental mutagenic factors, range size, number of competitors, or living on islands. Indeed, after controlling for the negative effect of species' age, 80% of variation in species‐specific evolutionary rates remains unexplained. At the clade level, high evolutionary rates are associated with unusual phenotypes or high species richness. Taken together, these results imply that macroevolutionary rates of ecomorphological traits are governed by both ecological opportunity in distinct adaptive zones and niche differentiation among closely related species.

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.

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

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

Study of morphological variation of northern Neotropical Ariidae reveals conservatism despite macrohabitat transitions

BACKGROUND

Morphological convergence triggered by trophic adaptations is a common pattern in adaptive radiations. The study of shape variation in an evolutionary context is usually restricted to well-studied fish models. We take advantage of the recently revised systematics of New World Ariidae and investigate skull shape evolution in six genera of northern Neotropical Ariidae. They constitute a lineage that diversified in the marine habitat but repeatedly adapted to freshwater habitats. 3D geometric morphometrics was applied for the first time in catfish skulls and phylogenetically informed statistical analyses were performed to test for the impact of habitat on skull diversification after habitat transition in this lineage.

RESULTS

We found that skull shape is conserved throughout phylogeny. A morphospace analysis revealed that freshwater and marine species occupy extreme ends of the first principal component axis and that they exhibit similar Procrustes variances. Yet freshwater species occupy the smallest shape space compared to marine and brackish species (based on partial disparity), and marine and freshwater species have the largest Procrustes distance to each other. We observed a single case of shape convergence as derived from 'C-metrics', which cannot be explained by the occupation of the same habitat.

CONCLUSIONS

Although Ariidae occupy such a broad spectrum of different habitats from sea to freshwater, the morphospace analysis and analyses of shape and co-variation with habitat in a phylogenetic context shows that conservatism dominates skull shape evolution among ariid genera.

The roles of inter- and intra-sexual selection in behavioral isolation between native and invasive pupfishes

Male-male competition and female mate choice may both play important roles in driving and maintaining reproductive isolation between species. When previously allopatric species come into secondary contact with each other due to introductions, they provide an opportunity to evaluate the identity and strength of reproductive isolating mechanisms. If reproductive isolation is not maintained, hybridization may occur. We examined how reproductive isolating mechanisms mediate hybridization between endemic populations of the Red River pupfish Cyprinodon rubrofluviatilis and the recently introduced sheepshead minnow C. variegatus. In lab-based dominance trials, males of both species won the same number of competitions. However, male C. rubrofluviatilis that won competitions were more aggressive than C. variegatus winners, and more aggression was needed to win against competitor C. variagatus than allopatric C. rubrofluviatilis. Duration of fights also differed based on the relatedness of the competitor. In dichotomous mate choice trials, there were no conspecific or heterospecific preferences expressed by females of either species. Our findings that male-male aggression differs between closely and distantly related groups, but female choice does not suggest that male-male competition may be the more likely mechanism to impede gene flow in this system.

Trophic specialization drives morphological evolution in sea snakes

Viviparous sea snakes are the most rapidly speciating reptiles known, yet the ecological factors underlying this radiation are poorly understood. Here, we reconstructed dated trees for 75% of sea snake species and quantified body shape (forebody relative to hindbody girth), maximum body length and trophic diversity to examine how dietary specialization has influenced morphological diversification in this rapid radiation. We show that sea snake body shape and size are strongly correlated with the proportion of burrowing prey in the diet. Specialist predators of burrowing eels have convergently evolved a 'microcephalic' morphotype with dramatically reduced forebody relative to hindbody girth and intermediate body length. By comparison, snakes that predominantly feed on burrowing gobies are generally short-bodied and small-headed, but there is no evidence of convergent evolution. The eel specialists also exhibit faster rates of size and shape evolution compared to all other sea snakes, including those that feed on gobies. Our results suggest that trophic specialization to particular burrowing prey (eels) has invoked strong selective pressures that manifest as predictable and rapid morphological changes. Further studies are needed to examine the genetic and developmental mechanisms underlying these dramatic morphological changes and assess their role in sea snake speciation.