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

Adaptive radiation is recognized by a rapid burst of phenotypic, ecological and species diversification. However, it is unknown whether different species within an adaptive radiation evolve reproductive isolation at different rates. We compared patterns of genetic differentiation between nascent species within an adaptive radiation of Cyprinodon pupfishes using genotyping by sequencing. Similar to classic adaptive radiations, this clade exhibits rapid morphological diversification rates and two species are novel trophic specialists, a scale‐eater and hard‐shelled prey specialist (durophage), yet the radiation is <10 000 years old. Both specialists and an abundant generalist species all coexist in the benthic zone of lakes on San Salvador Island, Bahamas. Based on 13 912 single‐nucleotide polymorphisms (SNPs), we found consistent differences in genetic differentiation between each specialist species and the generalist across seven lakes. The scale‐eater showed the greatest genetic differentiation and clustered by species across lakes, whereas durophage populations often clustered with sympatric generalist populations, consistent with parallel speciation across lakes. However, we found strong evidence of admixture between durophage populations in different lakes, supporting a single origin of this species and genome‐wide introgression with sympatric generalist populations. We conclude that the scale‐eater is further along the speciation‐with‐gene‐flow continuum than the durophage and suggest that different adaptive landscapes underlying these two niche environments drive variable progress towards speciation within the same habitat. Our previous measurements of fitness surfaces in these lakes support this conclusion: the scale‐eating fitness peak may be more distant than the durophage peak on the complex adaptive landscape driving adaptive radiation.     

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

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

Changes in growth rates of oral jaw elements produce evolutionary novelty in bahamian pupfish

To understand the origins of novelty and the evolution of biological diversity, it is important to investigate the processes that generate phenotypic variation from genotypic variation. A number of path‐breaking studies have revealed the genetic basis for phenotypic differences between distantly related taxa, but how qualitative change is produced during the early stages of divergence is largely unexplored. Here, we focus on striking differences in jaw morphology exhibited by three closely related sympatric pupfish species (genus Cyprinodon) from San Salvador Island, Bahamas as a basis for investigating the genetic sources of morphological variation in recently diverged species. San Salvador Island pupfish are trophically diverse and display derived jaw morphologies distinct from any other species in the genus. We illustrate these qualitative morphological differences between species with 3D‐reconstructed CT‐images and camera lucida drawings of the skulls of wild‐caught fish. Quantitative data representing the size of individual bony skull elements in wild fish show how qualitatively novel morphologies arise as a consequence of changes to the size and shape of individual skull elements, particularly the dentary, premaxilla, and maxilla bones associated with the oral jaws. Consistent with these comparative data is that the growth rate of individual bony skull elements, measured on a developmental time series of lab‐reared fish, differs between species. Our data provide a critical foundation for future studies developing San Salvador Cyprinodon pupfishes as a model system to understand the evolution and development of novel morphologies at the species level. J. Morphol. 277:935–947, 2016. © 2016 Wiley Periodicals, Inc.     

THE GENETIC BASIS OF A COMPLEX FUNCTIONAL SYSTEM

The relationship between form and function can have profound effects on evolutionary dynamics and such effects may differ for simple versus complex systems. In particular, functions produced by multiple structural configurations (many‐to‐one mapping, MTOM) may dampen constituent trade‐offs and promote diversification. Unfortunately, we lack information about the genetic architecture of MTOM functional systems. The skulls of teleost fishes contain both simple (lower jaw levers) as well as more complex (jaws modeled as 4‐bar linkages) functional systems within the same craniofacial unit. We examined the mapping of form to function and the genetic basis of these systems by identifying quantitative trait loci (QTL) in hybrids of two Lake Malawi cichlid species. Hybrid individuals exhibited novelty (transgressive segregation) in morphological components and function of the simple and complex jaw systems. Functional novelty was proportional to the prevalence of extreme morphologies in the simple levers; by contrast, recombination of parental morphologies produced transgression in the MTOM 4‐bar linkage. We found multiple loci of moderate effect and epistasis controlling jaw phenotypes in both the simple and complex systems, with less phenotypic variance explained by QTL for the 4‐bar. Genetic linkage between components of the simple and complex systems partly explains phenotypic correlations and may constrain functional evolution.     

On the Measurement of Ecological Novelty: Scale-Eating Pupfish Are Separated by 168 my from Other Scale-Eating Fishes

The colonization of new adaptive zones is widely recognized as one of the hallmarks of adaptive radiation. However, the adoption of novel resources during this process is rarely distinguished from phenotypic change because morphology is a common proxy for ecology. How can we quantify ecological novelty independent of phenotype? Our study is split into two parts: we first document a remarkable example of ecological novelty, scale-eating (lepidophagy), within a rapidly-evolving adaptive radiation of Cyprinodon pupfishes on San Salvador Island, Bahamas. This specialized predatory niche is known in several other fish groups, but is not found elsewhere among the 1,500 species of atherinomorphs. Second, we quantify this ecological novelty by measuring the time-calibrated phylogenetic distance in years to the most closely-related species with convergent ecology. We find that scale-eating pupfish are separated by 168 million years of evolution from the nearest scale-eating fish. We apply this approach to a variety of examples and highlight the frequent decoupling of ecological novelty from phenotypic divergence. We observe that novel ecology is not always tightly correlated with rates of phenotypic or species diversification, particularly within recent adaptive radiations, necessitating the use of additional measures of ecological novelty independent of phenotype.     

Fitness landscape of sympatric pupfishes a useful tool for visualizing speciation

Arguably the most useful model of evolution emerged from the mind of Sewall Wright when he invented the fitness landscape (Wright 1932 ). In a recent issue of Molecular Ecology, Martin & Feinstein ( 2014 ) investigate the genetics and demographic history of an adaptive radiation of pupfish on San Salvador Island. Since the founder species colonized the island 10 000 years ago, two descendent species have appeared and in several lakes all three species (a durophage, a scale‐eater, and the generalist ancestral form) coexist. The three species are thought to occupy three distinct fitness peaks. The durophage and generalists' peaks are close, whereas the scale‐eater's peak is predicted to be distant and separated from the other two by a deep valley. Consistent with this view, gene flow between the two species on close fitness peaks is greater than the gene flow between these two species and the third species on a more distant peak. Correspondingly, the inferred fitness landscape predicts progress towards speciation, with more limited separation of species on close peaks, and that speciation is more complete for the scale‐eater. The article provides an illustrative example of the power afforded by analysis of large numbers of SNPs for estimating key parameters underlying evolutionary divergence.     

Modular Skeletal Evolution in Sticklebacks Is Controlled by Additive and Clustered Quantitative Trait Loci

Understanding the genetic architecture of evolutionary change remains a long-standing goal in biology. In vertebrates, skeletal evolution has contributed greatly to adaptation in body form and function in response to changing ecological variables like diet and predation. Here we use genome-wide linkage mapping in threespine stickleback fish to investigate the genetic architecture of evolved changes in many armor and trophic traits. We identify >100 quantitative trait loci (QTL) controlling the pattern of serially repeating skeletal elements, including gill rakers, teeth, branchial bones, jaws, median fin spines, and vertebrae. We use this large collection of QTL to address long-standing questions about the anatomical specificity, genetic dominance, and genomic clustering of loci controlling skeletal differences in evolving populations. We find that most QTL (76%) that influence serially repeating skeletal elements have anatomically regional effects. In addition, most QTL (71%) have at least partially additive effects, regardless of whether the QTL controls evolved loss or gain of skeletal elements. Finally, many QTL with high LOD scores cluster on chromosomes 4, 20, and 21. These results identify a modular system that can control highly specific aspects of skeletal form. Because of the general additivity and genomic clustering of major QTL, concerted changes in both protective armor and trophic traits may occur when sticklebacks inherit either marine or freshwater alleles at linked or possible “supergene” regions of the stickleback genome. Further study of these regions will help identify the molecular basis of both modular and coordinated changes in the vertebrate skeleton.     

Selection on QTL and complex traits in complex environments

Understanding genetic variation for complex traits in heterogeneous environments is a fundamental problem in biology. In this issue of Molecular Ecology, Fournier‐Level et al. ( 2013 ) analyse quantitative trait loci (QTL) influencing ecologically important phenotypes in mapping populations of Arabidopsis thaliana grown in four habitats across its native European range. They used causal modelling to quantify the selective consequences of life history and morphological traits and QTL on components of fitness. They found phenology QTL colocalizing with known flowering time genes as well as novel loci. Most QTL influenced fitness via life history and size traits, rather than QTL having direct effects on fitness. Comparison of phenotypes among environments found no evidence for genetic trade‐offs for phenology or growth traits, but genetic trade‐offs for fitness resulted because flowering time had opposite fitness effects in different environments. These changes in QTL effects and selective consequences may maintain genetic variation among populations.     

QTL mapping reveals the genetic architecture of loci affecting pre- and post-zygotic isolating barriers in Louisiana Iris

ABSTRACT: BACKGROUND: Hybridization among Louisiana Irises has been well established and the genetic architecture of reproductive isolation is known to affect the potential for and the directionality of introgression between taxa. Here we use co-dominant markers to identify regions where QTL are located both within and between backcross maps to compare the genetic architecture of reproductive isolation and fitness traits across treatments and years. RESULTS: QTL mapping was used to elucidate the genetic architecture of reproductive isolation between Iris fulva and Iris brevicaulis. Homologous co-dominant EST-SSR markers scored in two backcross populations between I. fulva and I. brevicaulis were used to generate genetic linkage maps. These were used as the framework for mapping QTL associated with variation in 11 phenotypic traits likely responsible for reproductive isolation and fitness. QTL were dispersed throughout the genome, with the exception of one region of a single linkage group (LG) where QTL for flowering time, sterility, and fruit production clustered. In most cases, homologous QTL were not identified in both backcross populations, however, homologous QTL for flowering time, number of growth points per rhizome, number of nodes per inflorescence, and number of flowers per node were identified on several linkage groups. CONCLUSION: Two different traits affecting reproductive isolation, flowering time and sterility, exhibit different genetic architectures, with numerous QTL across the Iris genome controlling flowering time and fewer, less distributed QTL affecting sterility. QTL for traits affecting fitness are largely distributed across the genome with occasional overlap, especially on LG 4, where several QTL increasing fitness and decreasing sterility cluster. Given the distribution and effect direction of QTL affecting reproductive isolation and fitness, we have predicted genomic regions where introgression may be more likely to occur (those regions associated with an increase in fitness and unlinked to loci controlling reproductive isolation) and those that are less likely to exhibit introgression (those regions linked to traits decreasing fitness and reproductive isolation).     

Reproductive isolation among endemic pupfishes (Cyprinodon) on San Salvador Island,Bahamas: microsatellite evidence

Three pupfish (Cyprinodon) morphotypes (two endemic) occur in some of the young (6000 ypb) saline lakes on the Bahamian island of San Salvador. The ‘normal’ morph, a detritivore/omnivore, is not different in its general features from Cyprinodon variegatus from other Bahamian islands. ‘Bulldog’ is a scale‐eater/piscivore that preys upon normal pupfish, and ‘bozo’ is a specialized molluskivore. Reproductive isolation among these morphs is not predicted by the evolutionary biology of congeneric species because sympatry of even morphogically and ecologically quite divergent pupfishes has usually resulted in hybridization/introgression. Survey of variation at eight microsatellite loci reveals that sympatric normal and bulldog populations are genetically distinctive by several criteria, and are therefore likely reproductively isolated. The bulldog morph in Crescent Pond is markedly divergent from those in Little Lake and Osprey Lake, a finding consistent with, although it does not prove, separate parallel origins of this morphotype. The data also suggest that the bulldogs in the latter two lakes did not evolve by intralacustrine speciation from the current sympatric normal populations. Some of the genetic data suggest that the bozo morph may also be reproductively isolated from the other two pupfishes, but only a small, pooled sample of this rare morphotype was available, and the issue is not resolved. Isolating mechanisms between bulldog and normal morphs are of special interest because of the possibility that they arose as a consequence of a predator–prey relationship. A strong correlation between reproductive isolation and predator–prey interactions could provide an important example of ecological speciation via direct selection against heterotypic interactions. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 566–582.     

Deciphering the genomic architecture of the stickleback brain with a novel multilocus gene‐mapping approach

Quantitative traits important to organismal function and fitness, such as brain size, are presumably controlled by many small‐effect loci. Deciphering the genetic architecture of such traits with traditional quantitative trait locus (QTL) mapping methods is challenging. Here, we investigated the genetic architecture of brain size (and the size of five different brain parts) in nine‐spined sticklebacks (Pungitius pungitius) with the aid of novel multilocus QTL‐mapping approaches based on a de‐biased LASSO method. Apart from having more statistical power to detect QTL and reduced rate of false positives than conventional QTL‐mapping approaches, the developed methods can handle large marker panels and provide estimates of genomic heritability. Single‐locus analyses of an F₂ interpopulation cross with 239 individuals and 15 198, fully informative single nucleotide polymorphisms (SNPs) uncovered 79 QTL associated with variation in stickleback brain size traits. Many of these loci were in strong linkage disequilibrium (LD) with each other, and consequently, a multilocus mapping of individual SNPs, accounting for LD structure in the data, recovered only four significant QTL. However, a multilocus mapping of SNPs grouped by linkage group (LG) identified 14 LGs (1–6 depending on the trait) that influence variation in brain traits. For instance, 17.6% of the variation in relative brain size was explainable by cumulative effects of SNPs distributed over six LGs, whereas 42% of the variation was accounted for by all 21 LGs. Hence, the results suggest that variation in stickleback brain traits is influenced by many small‐effect loci. Apart from suggesting moderately heritable (h² ≈ 0.15–0.42) multifactorial genetic architecture of brain traits, the results highlight the challenges in identifying the loci contributing to variation in quantitative traits. Nevertheless, the results demonstrate that the novel QTL‐mapping approach developed here has distinctive advantages over the traditional QTL‐mapping methods in analyses of dense marker panels.     

Trophic apparatus in cyprinodontiform fishes: Functional specializations for picking and scraping behaviors

Cyprinodontiforms are a diverse and speciose order that includes topminnows, pupfishes, swordtails, mosquitofishes, guppies, and mollies. Sister group to the Beloniformes and Atheriniformes, Cyprinodontiformes contains approximately twice the number of species of these other two orders combined. Recent studies suggest that this group is well suited to capturing prey by “picking” small items from the water surface, water column, and the substrate. Because picking places unusual performance demands on the feeding apparatus, this mode of prey capture may rely upon novel morphological modifications not found in more widespread ram‐ or suction‐based feeding mechanisms. To assess this evolutionary hypothesis, we describe the trophic anatomy of 16 cyprinodontiform species, selected to broadly represent the order as well as capture intrageneric variation. The group appears to have undergone gradual morphological changes to become increasingly specialized for picking and scraping behaviors. We also identify a suite of functional characters related to the acquisition of a novel and previously undescribed mechanism of premaxillary protrusion and retraction, including: modification of the “premaxillomandibular” ligament (which connects each side of the premaxilla to the ipsilateral mandible, or lower jaw), a novel architecture of the ligaments and bony elements that unite the premaxillae, maxillae and palatine bones, and novel insertions of the adductor muscles onto the jaws. These morphological changes to both the upper and lower jaws suggest an evolutionary trend within this group toward increased reliance on picking individual prey from the water column/substrate or for scraping encrusting material from the substrate. We propose that the suite of morphological characters described here enable a functional innovation, “picking,” which leads to novel trophic habits. J. Morphol., 2009. © 2008 Wiley‐Liss, Inc.     

Paths to selection on life history loci in different natural environments across the native range of Arabidopsis thaliana

Selection on quantitative trait loci (QTL) may vary among natural environments due to differences in the genetic architecture of traits, environment‐specific allelic effects or changes in the direction and magnitude of selection on specific traits. To dissect the environmental differences in selection on life history QTL across climatic regions, we grew a panel of interconnected recombinant inbred lines (RILs) of Arabidopsis thaliana in four field sites across its native European range. For each environment, we mapped QTL for growth, reproductive timing and development. Several QTL were pleiotropic across environments, three colocalizing with known functional polymorphisms in flowering time genes (CRY2, FRI and MAF2‐5), but major QTL differed across field sites, showing conditional neutrality. We used structural equation models to trace selection paths from QTL to lifetime fitness in each environment. Only three QTL directly affected fruit number, measuring fitness. Most QTL had an indirect effect on fitness through their effect on bolting time or leaf length. Influence of life history traits on fitness differed dramatically across sites, resulting in different patterns of selection on reproductive timing and underlying QTL. In two oceanic field sites with high prereproductive mortality, QTL alleles contributing to early reproduction resulted in greater fruit production, conferring selective advantage, whereas alleles contributing to later reproduction resulted in larger size and higher fitness in a continental site. This demonstrates how environmental variation leads to change in both QTL effect sizes and direction of selection on traits, justifying the persistence of allelic polymorphism at life history QTL across the species range.     

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

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

GENOME DIVERGENCE AND THE GENETIC ARCHITECTURE OF BARRIERS TO GENE FLOW BETWEEN LYCAEIDES IDAS AND L. MELISSA

Genome divergence during speciation is a dynamic process that is affected by various factors, including the genetic architecture of barriers to gene flow. Herein we quantitatively describe aspects of the genetic architecture of two sets of traits, male genitalic morphology and oviposition preference, that putatively function as barriers to gene flow between the butterfly species Lycaeides idas and L. melissa. Our analyses are based on unmapped DNA sequence data and a recently developed Bayesian regression approach that includes variable selection and explicit parameters for the genetic architecture of traits. A modest number of nucleotide polymorphisms explained a small to large proportion of the variation in each trait, and average genetic variant effects were nonnegligible. Several genetic regions were associated with variation in multiple traits or with trait variation within‐ and among‐populations. In some instances, genetic regions associated with trait variation also exhibited exceptional genetic differentiation between species or exceptional introgression in hybrids. These results are consistent with the hypothesis that divergent selection on male genitalia has contributed to heterogeneous genetic differentiation, and that both sets of traits affect fitness in hybrids. Although these results are encouraging, we highlight several difficulties related to understanding the genetics of speciation.     

A quantitative trait locus analysis of personality in wild bighorn sheep

Personality, the presence of persistent behav105 ioral differences among individuals over time or contexts, potentially has important ecological and evolutionary consequences. However, a lack of knowledge about its genetic architecture limits our ability to understand its origin, evolution, and maintenance. Here, we report on a genome‐wide quantitative trait locus (QTL) analysis for two personality traits, docility and boldness, in free‐living female bighorn sheep from Ram Mountain, Alberta, Canada. Our variance component linkage analysis based on 238 microsatellite loci genotyped in 310 pedigreed individuals identified suggestive docility and boldness QTL on sheep chromosome 2 and 6, respectively. A lack of QTL overlap indicated that genetic covariance between traits was not modulated by pleiotropic effects at a major locus and may instead result from linkage disequilibrium or pleiotropic effects at QTL of small effects. To our knowledge, this study represents the first attempt to dissect the genetic architecture of personality in a free‐living wildlife population, an important step toward understanding the link between molecular genetic variation in personality and fitness and the evolutionary processes maintaining this variation.     

QTL mapping of freezing tolerance: links to fitness and adaptive trade‐offs

Local adaptation, defined as higher fitness of local vs. nonlocal genotypes, is commonly identified in reciprocal transplant experiments. Reciprocally adapted populations display fitness trade‐offs across environments, but little is known about the traits and genes underlying fitness trade‐offs in reciprocally adapted populations. We investigated the genetic basis and adaptive significance of freezing tolerance using locally adapted populations of Arabidopsis thaliana from Italy and Sweden. Previous reciprocal transplant studies of these populations indicated that subfreezing temperature is a major selective agent in Sweden. We used quantitative trait locus (QTL) mapping to identify the contribution of freezing tolerance to previously demonstrated local adaptation and genetic trade‐offs. First, we compared the genomic locations of freezing tolerance QTL to those for previously published QTL for survival in Sweden, and overall fitness in the field. Then, we estimated the contributions to survival and fitness across both field sites of genotypes at locally adaptive freezing tolerance QTL. In growth chamber studies, we found seven QTL for freezing tolerance, and the Swedish genotype increased freezing tolerance for five of these QTL. Three of these colocalized with locally adaptive survival QTL in Sweden and with trade‐off QTL for overall fitness. Two freezing tolerance QTL contribute to genetic trade‐offs across environments for both survival and overall fitness. A major regulator of freezing tolerance, CBF2, is implicated as a candidate gene for one of the trade‐off freezing tolerance QTL. Our study provides some of the first evidence of a trait and gene that mediate a fitness trade‐off in nature.     

Thyroid hormone modulation during zebrafish development recapitulates evolved diversity in danionin jaw protrusion mechanics

Protrusile jaws are a highly useful innovation that has been linked to extensive diversification in fish feeding ecology. Jaw protrusion can enhance the performance of multiple functions, such as suction production and capturing elusive prey. Identifying the developmental factors that alter protrusion ability will improve our understanding of fish diversification. In the zebrafish protrusion arises postmetamorphosis. Fish metamorphosis typically includes significant changes in trophic morphology, accompanies a shift in feeding niche and coincides with increased thyroid hormone production. We tested whether thyroid hormone affects the development of zebrafish feeding mechanics. We found that it affected all developmental stages examined, but that effects were most pronounced after metamorphosis. Thyroid hormone levels affected the development of jaw morphology, feeding mechanics, shape variation, and cranial ossification. Adult zebrafish utilize protrusile jaws, but an absence of thyroid hormone impaired development of the premaxillary bone, which is critical to jaw protrusion. Premaxillae from early juvenile zebrafish and hypothyroid adult zebrafish resemble those from adults in the genera Danionella, Devario, and Microdevario that show little to no jaw protrusion. Our findings suggest that evolutionary changes in how the developing skulls of danionin minnows respond to thyroid hormone may have promoted diversification into different feeding niches.