High levels of species' extirpation in an urban environment-A case study from Berlin, Germany, covering 1700-2023

Species loss is highly scale-dependent, following the species-area relationship. We analysed spatio-temporal patterns of species' extirpation on a multitaxonomic level using Berlin, the capital city of Germany. Berlin is one of the largest cities in Europe and has experienced a strong urbanisation trend since the late nineteenth century. We expected species' extirpation to be exceptionally high due to the long history of urbanisation. Analysing 37 regional Red Lists of Threatened Plants, Animals and Fungi of Berlin (covering 9498 species), we found that 16% of species were extirpated, a rate 5.9 times higher than at the German scale and 47.1 times higher than at the European scale. Species' extirpation in Berlin is comparable to that of another German city with a similarly broad taxonomic coverage, but much higher than in regional areas with less human impact. The documentation of species' extirpation started in the eighteenth century and is well documented for the nineteenth and twentieth centuries. We found an average annual extirpation of 3.6 species in the nineteenth century, 9.6 species in the twentieth century and the same number of extirpated species as in the nineteenth century were documented in the twenty-first century, despite the much shorter time period. Our results showed that species' extirpation is higher at small than on large spatial scales, and might be negatively influenced by urbanisation, with different effects on different taxonomic groups and habitats. Over time, we found that species' extirpation is highest during periods of high human alterations and is negatively affected by the number of people living in the city. But, there is still a lack of data to decouple the size of the area and the human impact of urbanisation. However, cities might be suitable systems for studying species' extirpation processes due to their small scale and human impact.

Sweepstakes reproduction facilitates rapid adaptation in highly fecund populations

Adaptation enables natural populations to survive in a changing environment. Understanding the mechanics of adaptation is therefore crucial for learning about the evolution and ecology of natural populations. We focus on the impact of random sweepstakes on selection in highly fecund haploid and diploid populations partitioned into two genetic types, with one type conferring selective advantage. For the diploid populations, we incorporate various dominance mechanisms. We assume that the populations may experience recurrent bottlenecks. In random sweepstakes, the distribution of individual recruitment success is highly skewed, resulting in a huge variance in the number of offspring contributed by the individuals present in any given generation. Using computer simulations, we investigate the joint effects of random sweepstakes, recurrent bottlenecks and dominance mechanisms on selection. In our framework, bottlenecks allow random sweepstakes to have an effect on the time to fixation, and in diploid populations, the effect of random sweepstakes depends on the dominance mechanism. We describe selective sweepstakes that are approximated by recurrent sweeps of strongly beneficial allelic types arising by mutation. We demonstrate that both types of sweepstakes reproduction may facilitate rapid adaptation (as defined based on the average time to fixation of a type conferring selective advantage conditioned on fixation of the type). However, whether random sweepstakes cause rapid adaptation depends also on their interactions with bottlenecks and dominance mechanisms. Finally, we review a case study in which a model of recurrent sweeps is shown to essentially explain population genomic data from Atlantic cod.

The diversification of butterfly wing patterns: progress and prospects

Butterfly wings display rich phenotypic diversity and are associated with complex biological functions, thus serving as an important evolutionary system to address the genetic basis and evolution of phenotypic diversification. We review recent butterfly studies that revealed complex functions underlying diversified wing patterns and describe the genetic and environmental factors involved in wing pattern determinations. These factors lead to inter-specific divergence, genetic polymorphism, and phenotypic plasticity, which in many cases are decided by several key genes. We also summarize the research advances on gene co-option as an important origin of functional complexity and evolutionary novelty. These findings reveal a pattern of evolutionary innovation within a constrained developmental framework during butterfly wing morphogenesis, but further research is required to gain a systematic and comprehensive understanding.

Genomic evidence for domestication selection in three hatchery populations of Chinook salmon, Oncorhynchus tshawytscha

Fish hatcheries are widely used to enhance fisheries and supplement declining wild populations. However, substantial evidence suggests that hatchery fish are subject to differential selection pressures compared to their wild counterparts. Domestication selection, or adaptation to the hatchery environment, poses a risk to wild populations if traits specific to success in the hatchery environment have a genetic component and there is subsequent introgression between hatchery and wild fish. Few studies have investigated domestication selection in hatcheries on a genomic level, and even fewer have done so in parallel across multiple hatchery-wild population pairs. In this study, we used low-coverage whole-genome sequencing to investigate signals of domestication selection in three separate hatchery populations of Chinook salmon, Oncorhynchus tshawytscha, after approximately seven generations of divergence from their corresponding wild progenitor populations. We sequenced 192 individuals from populations across Southeast Alaska and estimated genotype likelihoods at over six million loci. We discovered a total of 14 outlier peaks displaying high genetic differentiation (F ST) between hatchery-wild pairs, although no peaks were shared across the three comparisons. Peaks were small (53 kb on average) and often displayed elevated absolute genetic divergence (D xy) and linkage disequilibrium, suggesting some level of domestication selection has occurred. Our study provides evidence that domestication selection can lead to genetic differences between hatchery and wild populations in only a few generations. Additionally, our data suggest that population-specific adaptation to hatchery environments likely occurs through different genetic pathways, even for populations with similar standing genetic variation. These results highlight the need to collect paired genotype-phenotype data to understand how domestication may be affecting fitness and to identify potential management practices that may mitigate genetic risks despite multiple pathways of domestication.

Fire-driven animal evolution in the Pyrocene

Fire regimes are a major agent of evolution in terrestrial animals. Changing fire regimes and the capacity for rapid evolution in wild animal populations suggests the potential for rapid, fire-driven adaptive animal evolution in the Pyrocene. Fire drives multiple modes of evolutionary change, including stabilizing, directional, disruptive, and fluctuating selection, and can strongly influence gene flow and genetic drift. Ongoing and future research in fire-driven animal evolution will benefit from further development of generalizable hypotheses, studies conducted in highly responsive taxa, and linking fire-adapted phenotypes to their underlying genetic basis. A better understanding of evolutionary responses to fire has the potential to positively influence conservation strategies that embrace evolutionary resilience to fire in the Pyrocene.

High-quality haploid genomes corroborate 29 chromosomes and highly conserved synteny of genes in Hyles hawkmoths (Lepidoptera: Sphingidae)

BACKGROUND

Morphological and traditional genetic studies of the young Pliocene genus Hyles have led to the understanding that despite its importance for taxonomy, phenotypic similarity of wing patterns does not correlate with phylogenetic relationship. To gain insights into various aspects of speciation in the Spurge Hawkmoth (Hyles euphorbiae), we assembled a chromosome-level genome and investigated some of its characteristics.

RESULTS

The genome of a male H. euphorbiae was sequenced using PacBio and Hi-C data, yielding a 504 Mb assembly (scaffold N50 of 18.2 Mb) with 99.9% of data represented by the 29 largest scaffolds forming the haploid chromosome set. Consistent with this, FISH analysis of the karyotype revealed n = 29 chromosomes and a WZ/ZZ (female/male) sex chromosome system. Estimates of chromosome length based on the karyotype image provided an additional quality metric of assembled chromosome size. Rescaffolding the published male H. vespertilio genome resulted in a high-quality assembly (651 Mb, scaffold N50 of 22 Mb) with 98% of sequence data in the 29 chromosomes. The larger genome size of H. vespertilio (average 1C DNA value of 562 Mb) was accompanied by a proportional increase in repeats from 45% in H. euphorbiae (measured as 472 Mb) to almost 55% in H. vespertilio. Several wing pattern genes were found on the same chromosomes in the two species, with varying amounts and positions of repetitive elements and inversions possibly corrupting their function.

CONCLUSIONS

Our two-fold comparative genomics approach revealed high gene synteny of the Hyles genomes to other Sphingidae and high correspondence to intact Merian elements, the ancestral linkage groups of Lepidoptera, with the exception of three simple fusion events. We propose a standardized approach for genome taxonomy using nucleotide homology via scaffold chaining as the primary tool combined with Oxford plots based on Merian elements to infer and visualize directionality of chromosomal rearrangements. The identification of wing pattern genes promises future understanding of the evolution of forewing patterns in the genus Hyles, although further sequencing data from more individuals are needed. The genomic data obtained provide additional reliable references for further comparative studies in hawkmoths (Sphingidae).

Genetic Architecture of Flowering Time Differs Between Populations With Contrasting Demographic and Selective Histories

Understanding the evolutionary factors that impact the genetic architecture of traits is a central goal of evolutionary genetics. Here, we investigate how quantitative trait variation accumulated over time in populations that colonized a novel environment. We compare the genetic architecture of flowering time in Arabidopsis populations from the drought-prone Cape Verde Islands and their closest outgroup population from North Africa. We find that trait polygenicity is severely reduced in the island populations compared to the continental North African population. Further, trait architectures and reconstructed allelic histories best fit a model of strong directional selection in the islands in accord with a Fisher-Orr adaptive walk. Consistent with this, we find that large-effect variants that disrupt major flowering time genes (FRI and FLC) arose first, followed by smaller effect variants, including ATX2 L125F, which is associated with a 4-day reduction in flowering time. The most recently arising flowering time-associated loci are not known to be directly involved in flowering time, consistent with an omnigenic signature developing as the population approaches its trait optimum. Surprisingly, we find no effect in the natural population of EDI-Cvi-0 (CRY2 V367M), an allele for which an effect was previously validated by introgression into a Eurasian line. Instead, our results suggest the previously observed effect of the EDI-Cvi-0 allele on flowering time likely depends on genetic background, due to an epistatic interaction. Altogether, our results provide an empirical example of the effects demographic history and selection has on trait architecture.

Species- and sex-dependent changes in body size between 1892 and 2017, and recent biochemical signatures in rural and urban populations of two ground beetle species

Increasing urbanisation and intensified agriculture lead to rapid transitions of ecosystems. Species that persist throughout rapid transitions may respond to environmental changes across space and/or time, for instance by altering morphological and/or biochemical traits. We used natural history museum specimens, covering the Anthropocene epoch, to obtain long-term data combined with recent samples. We tested whether rural and urban populations of two ground beetle species, Harpalus affinis and H. rufipes, exhibit spatio-temporal intraspecific differences in body size. On a spatial scale, we tested signatures of nitrogen and carbon stable isotopes enrichments in different tissues and body components in recent populations of both species from urban and agricultural habitats. For body size examinations, we used beetles, collected from the early 20th century until 2017 in the Berlin-Brandenburg region, Germany, where urbanisation and agriculture have intensified throughout the last century. For stable isotope examinations, we used recent beetles from urban and agricultural habitats. Our results revealed no spatio-temporal changes in body size in both species' females. Body size of H. rufipes males decreased in the city but remained constant in rural areas over time. We discuss our findings with respect to habitat quality, urban heat and interspecific differences in activity pattern. Although nitrogen isotope ratios were mostly higher in specimens from agricultural habitats, some urban beetles reached equal enrichments. Carbon signatures of both species did not differ between habitats, detecting no differences in energy sources. Our results indicate that increasing urbanisation and intensified agriculture are influencing species' morphology and/or biochemistry. However, changes may be species- and sex-specific.

Development and Validation of the Questionnaire for Adaptive Hyperactivity and Goal Achievement (AHGA)

Objective

This paper illustrates the preliminary psychometric properties of the Questionnaire for Adaptive Hyperactivity and Goal Achievement (AHGA), aimed at measuring adaptive characteristics of hyperactivity and goal pursuit in older adults.

Methods

The 12-item scale was administered to a sample of 120 subjects (older adults) between February 2022 and June 2022. The reliability of AHGA was measured using Cronbach's alpha, and factor structure was established using parallel analysis (PA) and principal component analysis (PCA). Convergent validity was tested against the Biological Rhythms Interview of Assessment in Neuropsychiatry (BRIAN).

Results

All included subjects have an average age of 74.1±5.1 years. AHGA reliability was good (Cronbach's alpha: 0.713 [95%CI: 0.630 to 0.783]). Factor analysis suggested two main components: goal achievement and hyperactivity, which explained 41% of the variance in the data. The results support the convergent validity of the scale: AHGA measures adaptive characteristics of hyperactivity and goal pursuit, in contrast to BRIAN, which measures pathological characteristics.

Conclusion

The reported findings represent an innovative approach to hyperthymic features by embracing a broader spectrum concept that conceptualizes the potential transition between pathological and adaptive aspects as a continuum.

Increasing ecological heterogeneity can constrain biopesticide resistance evolution

Microbial biopesticides containing living parasites are valuable emerging crop protection technologies against insect pests, but they are vulnerable to resistance evolution. Fortunately, the fitness of alleles that provide resistance, including to parasites used in biopesticides, frequently depends on parasite identity and environmental conditions. This context-specificity suggests a sustainable approach to biopesticide resistance management through landscape diversification. To mitigate resistance risks, we advocate increasing the range of biopesticides available to farmers, whilst simultaneously encouraging other aspects of landscape-wide crop heterogeneity that can generate variable selection on resistance alleles. This approach requires agricultural stakeholders to prioritize diversity as well as efficiency, both within agricultural landscapes and the biocontrol marketplace.

Polygenic adaptation dynamics in large, finite populations

Although many phenotypic traits are determined by a large number of genetic variants, how a polygenic trait adapts in response to a change in the environment is not completely understood. In the framework of diffusion theory, we study the steady state and the adaptation dynamics of a large but finite population evolving under stabilizing selection and symmetric mutations when selection and mutation are moderately large. We find that in the stationary state, the allele frequency distribution at a locus is unimodal if its effect size is below a threshold effect and bimodal otherwise; these results are the stochastic analog of the deterministic ones where the stable allele frequency becomes bistable when the effect size exceeds a threshold. It is known that following a sudden shift in the phenotypic optimum, in an infinitely large population, selective sweeps at a large-effect locus are prevented and adaptation proceeds exclusively via subtle changes in the allele frequency; in contrast, we find that the chance of sweep is substantially enhanced in large, finite populations and the allele frequency at a large-effect locus can reach a high frequency at short times even for small shifts in the phenotypic optimum.

The role of spatial structure in multi-deme models of evolutionary rescue

Genetic variation and population sizes are critical factors for successful adaptation to novel environmental conditions. Gene flow between sub-populations is a potent mechanism to provide such variation and can hence facilitate adaptation, for instance by increasing genetic variation or via the introduction of beneficial variants. On the other hand, if gene flow between different habitats is too strong, locally beneficial alleles may not be able to establish permanently. In the context of evolutionary rescue, intermediate levels of gene flow are therefore often optimal for maximizing a species chance for survival in metapopulations without spatial structure. To which extent and under which conditions gene flow facilitates or hinders evolutionary rescue in spatially structured populations remains unresolved. We address this question by studying the differences between evolutionary rescue in the island model and in the stepping stone model in a gradually deteriorating habitat. We show that evolutionary rescue is modulated by the rate of gene flow between different habitats, which in turn depends strongly on the spatial structure and the pattern of environmental deterioration. We use these insights to show that in many cases spatially structured models can be translated into a simpler island model using an appropriately scaled effective migration rate.

Adaptation to urban environments

Despite widespread evidence of urban evolution, the adaptive nature of these changes is often unclear. We review different phenotypic and molecular lines of evidence used for assessing urban adaptation, discussing the benefits and limitations of each approach, and rare examples of their integration. We then provide a synthesis of local adaptation to urban and rural environments. These data were drawn from phenotypic reciprocal transplant studies, the majority of which focus on insects and other arthropods. Broadly, we found support for local adaptation to urban and rural environments. However, there was asymmetry in the evidence for local adaptation depending on population of origin, with urban adaptation being less prevalent than rural adaptation, suggesting many urban populations are still adapting to urban environments. Further, the general patterns were underlain by considerable variation among study systems; we discuss how environmental heterogeneity and costs of adaptation might explain system-specific variation in urban-rural local adaptation. We then look to the future of urban adaptation research, considering the magnitude and direction of adaptation in context of different agents of selection including urban heat islands, chemical pollutants, and biotic interactions.

The evolution of larger size in high-altitude Drosophila melanogaster has a variable genetic architecture

Important uncertainties persist regarding the genetic architecture of adaptive trait evolution in natural populations, including the number of genetic variants involved, whether they are drawn from standing genetic variation, and whether directional selection drives them to complete fixation. Here, we take advantage of a unique natural population of Drosophila melanogaster from the Ethiopian highlands, which has evolved larger body size than any other known population of this species. We apply a bulk segregant quantitative trait locus mapping approach to 4 unique crosses between highland Ethiopian and lowland Zambian populations for both thorax length and wing length. Results indicated a persistently variable genetic basis for these evolved traits (with largely distinct sets of quantitative trait loci for each cross), and at least a moderately polygenic architecture with relatively strong effects present. We complemented these mapping experiments with population genetic analyses of quantitative trait locus regions and gene ontology enrichment analysis, generating strong hypotheses for specific genes and functional processes that may have contributed to these adaptive trait changes. Finally, we find that the genetic architectures indicated by our quantitative trait locus mapping results for size traits mirror those from similar experiments on other recently evolved traits in this species. Collectively, these studies suggest a recurring pattern of polygenic adaptation in this species, in which causative variants do not approach fixation and moderately strong effect loci are present.

Silencing of Adc and Ebony Causes Abnormal Darkening of Cuticle in Henosepilachna vigintioctopunctata

N-β-alanyldopamine (NBAD) is a precursor of N-acylquinone sclerotin utilized for cross-linking between cuticular proteins for cuticle during insect molting. The importance of NBAD in cuticle tanning has not been well compared among different developing stages of insects. Henosepilachna vigintioctopunctata, a typical polyphagous pest feeding on a large number of Solanaceae and Cucurbitaceae plants in Asian countries, displays diverse cuticle pigmentation patterns among developing stages and body regions. Here, we found that the expression of three genes (Hvadc, Hvebony, and Hvtan) involved in NBAD biosynthesis peaked in the 4-day-old pupae or 0-day-old adults of H. vigintioctopunctata. At the first, second, third, and fourth larval instar and pupal stage, their transcript levels were high just before and/or right after the molting. Moreover, they were more abundantly transcribed at the larval heads than in the bodies. RNA interference (RNAi) of either Hvadc or Hvebony at the third instar larvae selectively deepened the color of the larval head capsules, antennae, mouthpart, scoli, strumae, and legs; and depletion of the two genes blackened the pupal head capsules, antennae, mouthpart, and legs. However, the knockdown of either Hvadc or Hvebony darkened the whole bodies of the adults. Conversely, RNAi of Hvtan at the third instar stage had little influence on the pigmentation in the larvae, pupae, and adults. These findings demonstrated that Adc and Ebony are important in cuticle pigmentation of H. vigintioctopunctata and suggested that larger quantities of NBAD were present in adults and play more important roles in pigmentation than larvae/pupae.

Complex adaptive architecture underlies adaptation to quantitative host resistance in a fungal plant pathogen

Plant pathogens often adapt to plant genetic resistance so characterization of the architecture underlying such an adaptation is required to understand the adaptive potential of pathogen populations. Erosion of banana quantitative resistance to a major leaf disease caused by polygenic adaptation of the causal agent, the fungus Pseudocercospora fijiensis, was recently identified in the northern Caribbean region. Genome scan and quantitative genetics approaches were combined to investigate the adaptive architecture underlying this adaptation. Thirty-two genomic regions showing host selection footprints were identified by pool sequencing of isolates collected from seven plantation pairs of two cultivars with different levels of quantitative resistance. Individual sequencing and phenotyping of isolates from one pair revealed significant and variable levels of correlation between haplotypes in 17 of these regions with a quantitative trait of pathogenicity (the diseased leaf area). The multilocus pattern of haplotypes detected in the 17 regions was found to be highly variable across all the population pairs studied. These results suggest complex adaptive architecture underlying plant pathogen adaptation to quantitative resistance with a polygenic basis, redundancy, and a low level of parallel evolution between pathogen populations. Candidate genes involved in quantitative pathogenicity and host adaptation of P. fijiensis were identified in genomic regions by combining annotation analysis with available biological data.

Ethanol resistance in Drosophila melanogaster has increased in parallel cold-adapted populations and shows a variable genetic architecture within and between populations

Understanding the genetic properties of adaptive trait evolution is a fundamental crux of biological inquiry that links molecular processes to biological diversity. Important uncertainties persist regarding the genetic predictability of adaptive trait change, the role of standing variation, and whether adaptation tends to result in the fixation of favored variants. Here, we use the recurrent evolution of enhanced ethanol resistance in Drosophila melanogaster during this species' worldwide expansion as a promising system to add to our understanding of the genetics of adaptation. We find that elevated ethanol resistance has evolved at least three times in different cooler regions of the species' modern range-not only at high latitude but also in two African high-altitude regions. Applying a bulk segregant mapping framework, we find that the genetic architecture of ethanol resistance evolution differs substantially not only between our three resistant populations, but also between two crosses involving the same European population. We then apply population genetic scans for local adaptation within our quantitative trait locus regions, and we find potential contributions of genes with annotated roles in spindle localization, membrane composition, sterol and alcohol metabolism, and other processes. We also apply simulation-based analyses that confirm the variable genetic basis of ethanol resistance and hint at a moderately polygenic architecture. However, these simulations indicate that larger-scale studies will be needed to more clearly quantify the genetic architecture of adaptive evolution and to firmly connect trait evolution to specific causative loci.

Rapid Evolutionary Adaptation in Response to Selection on Quantitative Traits

Evolutionary adaptation after sudden environmental changes can occur very rapidly. The mechanisms facilitating rapid adaptation range from strong positive directional selection leading to large shifts in the allele frequencies at a few loci (selective sweeps) to polygenic selection causing small changes in allele frequencies at many loci. In addition, combinations of these two extreme mechanisms may also result in fast evolution. In recent years, following reports of new case studies of rapid adaptation, population genetic models have been proposed to explain these observations. In these models, the role of the major selective forces (positive directional and stabilizing selection) is highlighted as well as the genetic architecture of quantitative traits. Furthermore, the factors limiting the speed of adaptation are analyzed, in particular, the effects of random genetic drift and demography due to finite population size.

Many ways to make darker flies: Intra- and interspecific variation in Drosophila body pigmentation components

Body pigmentation is an evolutionarily diversified and ecologically relevant trait with substantial variation within and between species, and important roles in animal survival and reproduction. Insect pigmentation, in particular, provides some of the most compelling examples of adaptive evolution, including its ecological significance and genetic bases. Pigmentation includes multiple aspects of color and color pattern that may vary more or less independently, and can be under different selective pressures. We decompose Drosophila thorax and abdominal pigmentation, a valuable eco-evo-devo model, into distinct measurable traits related to color and color pattern. We investigate intra- and interspecific variation for those traits and assess its different sources. For each body part, we measured overall darkness, as well as four other pigmentation properties distinguishing between background color and color of the darker pattern elements that decorate each body part. By focusing on two standard D. melanogaster laboratory populations, we show that pigmentation components vary and covary in distinct manners depending on sex, genetic background, and temperature during development. Studying three natural populations of D. melanogaster along a latitudinal cline and five other Drosophila species, we then show that evolution of lighter or darker bodies can be achieved by changing distinct component traits. Our results paint a much more complex picture of body pigmentation variation than previous studies could uncover, including patterns of sexual dimorphism, thermal plasticity, and interspecific diversity. These findings underscore the value of detailed quantitative phenotyping and analysis of different sources of variation for a better understanding of phenotypic variation and diversification, and the ecological pressures and genetic mechanisms underlying them.

Chromosome-level genome reference and genome editing of the tea geometrid

The tea geometrid is a destructive insect pest on tea plants, which seriously affects tea production in terms of both yield and quality and causes severe economic losses. The tea geometrid also provides an important study system to address the ecological adaptive mechanisms underlying its unique host plant adaptation and protective resemblance. In this study, we fully sequenced and de novo assembled the reference genome of the tea geometrid, Ectropis grisescens, using long sequencing reads. We presented a highly continuous, near-complete genome reference (787.4 Mb; scaffold N50: 26.9 Mb), along with the annotation of 18,746 protein-coding genes and 53.3% repeat contents. Importantly, we successfully placed 97.8% of the assembly in 31 chromosomes based on Hi-C interactions and characterized the sex chromosome based on sex-biased sequencing coverage. Multiple quality-control assays and chromosome-scale synteny with the model species all supported the high quality of the presented genome reference. We focused biological annotations on gene families related to the host plant adaptation and camouflage in the tea geometrid and performed comparisons with other representative lepidopteran species. Important findings include the E. grisescens-specific expansion of CYP6 P450 genes that might be involved in metabolism of tea defensive chemicals and unexpected massive expansion of gustatory receptor gene families that suggests potential polyphagy for this tea pest. Furthermore, we developed an efficient genome editing system based on CRISPR/Cas9 technology and successfully implement mutagenesis of a Hox gene in the tea geometrid. Our study provides key genomic resources both for exploring unique mechanisms underlying the ecological adaptation of tea geometrids and for developing environment-friendly strategies for tea pest management.

Evolutionary developmental biology and sustainability: A biology of resilience

Evolutionary developmental biology, and especially ecological developmental biology, is essential for discussions of sustainability and the responses to global climate change. First, this paper explores examples of animals that have successfully altered their development to accommodate human-made changes to their environments. We next document the ability of global warming to disrupt the development of those organisms with temperature-dependent sex-determination or with phenologies coordinating that organism's development with those of other species. The thermotolerance of Homo sapiens is also related to key developmental factors concerning brain development and maintenance, and the development of corals, the keystone organisms of tropical reefs, is discussed in relation to global warming as well as to other anthropogenic changes. While teratogenic and endocrine-disrupting compounds are not discussed in this essay, the ability of glyphosate herbicides to block insect development is highlighted. Last, the paper discusses the need to creatively integrate developmental biology with ecological, political, religious, and economic perspectives, as the flourishing of contemporary species may require altering the ways that Western science has considered the categories of nature, culture, and self.

Dissecting the eQTL Micro-Architecture in Caenorhabditis elegans

The study of expression quantitative trait loci (eQTL) using natural variation in inbred populations has yielded detailed information about the transcriptional regulation of complex traits. Studies on eQTL using recombinant inbred lines (RILs) led to insights on cis and trans regulatory loci of transcript abundance. However, determining the underlying causal polymorphic genes or variants is difficult, but ultimately essential for the understanding of regulatory networks of complex traits. This requires insight into whether associated loci are single eQTL or a combination of closely linked eQTL, and how this QTL micro-architecture depends on the environment. We addressed these questions by testing for independent replication of previously mapped eQTL in Caenorhabditis elegans using new data from introgression lines (ILs). Both populations indicate that the overall heritability of gene expression, number, and position of eQTL differed among environments. Across environments we were able to replicate 70% of the cis- and 40% of the trans-eQTL using the ILs. Testing eight different simulation models, we suggest that additive effects explain up to 60-93% of RIL/IL heritability for all three environments. Closely linked eQTL explained up to 40% of RIL/IL heritability in the control environment whereas only 7% in the heat-stress and recovery environments. In conclusion, we show that reproducibility of eQTL was higher for cis vs. trans eQTL and that the environment affects the eQTL micro-architecture.

High-resolution characterization of male ornamentation and re-evaluation of sex linkage in guppies

Coloration plays a key role in the ecology of many species, influencing how an organism interacts with its environment, other species and conspecifics. Guppies are sexually dimorphic, with males displaying sexually selected coloration resulting from female preference. Previous work has suggested that much of guppy colour pattern variation is Y-linked. However, it remains unclear how many individual colour patterns are Y-linked in natural populations as much of the previous work has focused on phenotypes either not found in the wild, or aggregate measures such as total colour area. Moreover, ornaments have traditionally been identified and delineated by hand, and computational methods now make it possible to extract pixels and identify ornaments with automated methods, reducing the potential for human bias. Here we developed a pipeline for semi-automated ornament identification and high-resolution image analysis of male guppy colour patterns and applied it to a multigenerational pedigree. Our results show that loci controlling the presence or the absence of individual male ornaments in our population are not predominantly Y-linked. However, we find that ornaments of similar colour are not independent of each other, and modifier loci that affect whole animal coloration appear to be at least partially Y-linked. Considering these results, Y-linkage of individual ornaments may not be important in driving colour changes in natural populations of guppies, or in expansions of the non-recombining Y region, while Y-linked modifier loci that affect aggregate traits may well play an important role.

Dental Calculus as a Tool to Study the Evolution of the Mammalian Oral Microbiome

Dental calculus, the calcified form of the mammalian oral microbial plaque biofilm, is a rich source of oral microbiome, host, and dietary biomolecules and is well preserved in museum and archaeological specimens. Despite its wide presence in mammals, to date, dental calculus has primarily been used to study primate microbiome evolution. We establish dental calculus as a valuable tool for the study of nonhuman host microbiome evolution, by using shotgun metagenomics to characterize the taxonomic and functional composition of the oral microbiome in species as diverse as gorillas, bears, and reindeer. We detect oral pathogens in individuals with evidence of oral disease, assemble near-complete bacterial genomes from historical specimens, characterize antibiotic resistance genes, reconstruct components of the host diet, and recover host genetic profiles. Our work demonstrates that metagenomic analyses of dental calculus can be performed on a diverse range of mammalian species, which will allow the study of oral microbiome and pathogen evolution from a comparative perspective. As dental calculus is readily preserved through time, it can also facilitate the quantification of the impact of anthropogenic changes on wildlife and the environment.

Polygenic Adaptation in a Population of Finite Size

Polygenic adaptation in response to selection on quantitative traits has become an important topic in evolutionary biology. Here we review the recent literature on models of polygenic adaptation. In particular, we focus on a model that includes mutation and both directional and stabilizing selection on a highly polygenic trait in a population of finite size (thus experiencing random genetic drift). Assuming that a sudden environmental shift of the fitness optimum occurs while the population is in a stochastic equilibrium, we analyze the adaptation of the trait to the new optimum. When the shift is not too large relative to the equilibrium genetic variance and this variance is determined by loci with mostly small effects, the approach of the mean phenotype to the optimum can be approximated by a rapid exponential process (whose rate is proportional to the genetic variance). During this rapid phase the underlying changes to allele frequencies, however, may depend strongly on genetic drift. While trait-increasing alleles with intermediate equilibrium frequencies are dominated by selection and contribute positively to changes of the trait mean (i.e., are aligned with the direction of the optimum shift), alleles with low or high equilibrium frequencies show more of a random dynamics, which is expected when drift is dominating. A strong effect of drift is also predicted for population size bottlenecks. Our simulations show that the presence of a bottleneck results in a larger deviation of the population mean of the trait from the fitness optimum, which suggests that more loci experience the influence of drift.

Important role of genetic drift in rapid polygenic adaptation

We analyzed a model to determine the factors that facilitate or limit rapid polygenic adaptation. This model includes population genetic terms of mutation and both directional and stabilizing selection on a highly polygenic trait in a diploid population of finite size. First, we derived the equilibrium distribution of the allele frequencies of the multilocus model by diffusion approximation. This formula describing the equilibrium allele frequencies as a mutation-selection-drift balance was examined by computer simulation using parameter values inferred for human height, a well-studied polygenic trait. Second, assuming that a sudden environmental shift of the fitness optimum occurs while the population is in equilibrium, we analyzed the adaptation of the trait to the new optimum. The speed at which the trait mean approaches the new optimum increases with the equilibrium genetic variance. Thus, large population size and/or large mutation rate may facilitate rapid adaptation. Third, the contribution of an individual locus i to polygenic adaptation depends on the compound parameterγ i p i 0 q i 0 , where γ i is the effect size,p i 0 the equilibrium frequency of the trait-increasing allele of this locus, andq i 0 = 1 - p i 0 . Thus, only loci with large values of this parameter contribute coherently to polygenic adaptation. Given that mutation rates are relatively small, this is more likely in large populations, in which the effects of drift are limited.

Field cricket genome reveals the footprint of recent, abrupt adaptation in the wild

Evolutionary adaptation is generally thought to occur through incremental mutational steps, but large mutational leaps can occur during its early stages. These are challenging to study in nature due to the difficulty of observing new genetic variants as they arise and spread, but characterizing their genomic dynamics is important for understanding factors favoring rapid adaptation. Here, we report genomic consequences of recent, adaptive song loss in a Hawaiian population of field crickets (Teleogryllus oceanicus). A discrete genetic variant, flatwing, appeared and spread approximately 15 years ago. Flatwing erases sound‐producing veins on male wings. These silent flatwing males are protected from a lethal, eavesdropping parasitoid fly. We sequenced, assembled and annotated the cricket genome, produced a linkage map, and identified a flatwing quantitative trait locus covering a large region of the X chromosome. Gene expression profiling showed that flatwing is associated with extensive genome‐wide effects on embryonic gene expression. We found that flatwing male crickets express feminized chemical pheromones. This male feminizing effect, on a different sexual signaling modality, is genetically associated with the flatwing genotype. Our findings suggest that the early stages of evolutionary adaptation to extreme pressures can be accompanied by greater genomic and phenotypic disruption than previously appreciated, and highlight how abrupt adaptation might involve suites of traits that arise through pleiotropy or genomic hitchhiking.

Genetic convergence of industrial melanism in three geometrid moths

The rise of dark (melanic) forms of many species of moth in heavily coal-polluted areas of nineteenth- and twentieth-century Britain, and their post-1970s fall, point to a common selective pressure (camouflage against bird predators) acting at the community level. The extent to which this convergent phenotypic response relied on similar genetic and developmental mechanisms is unknown. We examine this problem by testing the hypothesis that the locus controlling melanism in Phigalia pilosaria and Odontopera bidentata, two species of geometrid moth that showed strong associations between melanism and coal pollution, is the same as that controlling melanism in Biston betularia, previously identified as the gene cortex. Comparative linkage mapping using family material supports the hypothesis for both species, indicating a deeply conserved developmental mechanism for melanism involving cortex. However, in contrast to the strong selective sweep signature seen in British B. betularia, no significant association was detected between cortex-region markers and melanic morphs in wild-caught samples of P. pilosaria and O. bidentata, implying much older, or diverse, origins of melanic morph alleles in these latter species.

The genetic architecture of adaptation: convergence and pleiotropy in Heliconius wing pattern evolution

Unravelling the genetic basis of adaptive traits is a major challenge in evolutionary biology. Doing so informs our understanding of evolution towards an adaptive optimum, the distribution of locus effect sizes, and the influence of genetic architecture on the evolvability of a trait. In the Müllerian co-mimics Heliconius melpomene and Heliconius erato some Mendelian loci affecting mimicry shifts are well known. However, several phenotypes in H. melpomene remain to be mapped, and the quantitative genetics of colour pattern variation has rarely been analysed. Here we use quantitative trait loci (QTL) analyses of crosses between H. melpomene races from Peru and Suriname to map, for the first time, the control of the broken band phenotype to WntA and identify a ~100 kb region controlling this variation. Additionally, we map variation in basal forewing red-orange pigmentation to a locus centred around the gene ventral veins lacking (vvl). The locus also appears to affect medial band shape variation as it was previously known to do in H. erato. This adds to the list of homologous regions controlling convergent phenotypes between these two species. Finally we show that Heliconius wing-patterning genes are strikingly pleiotropic among wing pattern traits. Our results demonstrate how genetic architecture can shape, aid and constrain adaptive evolution.

Local adaptation does not lead to genome-wide differentiation in lava flow lizards

Adaptation can occur with or without genome-wide differentiation. If adaptive loci are linked to traits involved in reproductive isolation, genome-wide divergence is likely, and speciation is possible. However, adaptation can also lead to phenotypic differentiation without genome-wide divergence if levels of ongoing gene flow are high. Here, we use the replicated occurrence of melanism in lava flow lizards to assess the relationship between local adaptation and genome-wide differentiation. We compare patterns of phenotypic and genomic divergence among lava flow and nonlava populations for three lizard species and three lava flows in the Chihuahuan Desert. We find that local phenotypic adaptation (melanism) is not typically accompanied by genome-wide differentiation. Specifically, lava populations do not generally exhibit greater divergence from nonlava populations than expected by geography alone, regardless of whether the lava formation is 5,000 or 760,000 years old. We also infer that gene flow between lava and nonlava populations is ongoing in all lava populations surveyed. Recent work in the isolation by environment and ecological speciation literature suggests that environmentally driven genome-wide differentiation is common in nature. However, local adaptation may often simply be local adaptation rather than an early stage of ecological speciation.

Evolutionary genomics of gypsy moth populations sampled along a latitudinal gradient

The European gypsy moth (Lymantria dispar L.) was first introduced to Massachusetts in 1869 and within 150 years has spread throughout eastern North America. This large-scale invasion across a heterogeneous landscape allows examination of the genetic signatures of adaptation potentially associated with rapid geographical spread. We tested the hypothesis that spatially divergent natural selection has driven observed changes in three developmental traits that were measured in a common garden for 165 adult moths sampled from six populations across a latitudinal gradient covering the entirety of the range. We generated genotype data for 91,468 single nucleotide polymorphisms based on double digest restriction-site associated DNA sequencing and used these data to discover genome-wide associations for each trait, as well as to test for signatures of selection on the discovered architectures. Genetic structure across the introduced range of gypsy moth was low in magnitude (F_ST  = 0.069), with signatures of bottlenecks and spatial expansion apparent in the rare portion of the allele frequency spectrum. Results from applications of Bayesian sparse linear mixed models were consistent with the presumed polygenic architectures of each trait. Further analyses indicated spatially divergent natural selection acting on larval development time and pupal mass, with the linkage disequilibrium component of this test acting as the main driver of observed patterns. The populations most important for these signals were two range-edge populations established less than 30 generations ago. We discuss the importance of rapid polygenic adaptation to the ability of non-native species to invade novel environments.

Testing the potential mechanisms for the maintenance of a genetic color polymorphism in bluefin killifish populations

The maintenance of genetic variation in the face of natural selection is a long-standing question in evolutionary biology. In the bluefin killifish Lucania goodei, male coloration is polymorphic. Males can produce either red or yellow coloration in their anal fins, and both color morphs are present in all springs. These 2 morphs are heritable and how they are maintained in nature is unknown. Here, we tested 2 mechanisms for the maintenance of the red/yellow color morphs. Negative frequency-dependent mating success predicts that rare males have a mating advantage over common males. Spatial variation in fitness predicts that different color morphs have an advantage in different microhabitat types. Using a breeding experiment, we tested these hypotheses by creating populations with different ratios of red to yellow males (5 red:1 yellow; 1 red:5 yellow) and determining male mating success on shallow and deep spawning substrates. We found no evidence of negative frequency-dependent mating success. Common morphs tended to have higher mating success, and this was particularly so on shallow spawning substrates. However, on deep substrates, red males enjoyed higher mating success than yellow males, particularly so when red males were rare. However, yellow males did not have an advantage at either depth nor when rare. We suggest that preference for red males is expressed in deeper water, possibly due to alterations in the lighting environment. Finally, male pigment levels were correlated with one another and predicted male mating success. Hence, pigmentation plays an important role in male mating success.

Thermal Physiology and Developmental Plasticity of Pigmentation in the Harlequin Bug (Hemiptera: Pentatomidae)

Traits that promote the maintenance of body temperatures within an optimal range provide advantages to ectothermic species. Pigmentation plasticity is found in many insects and enhances thermoregulatory potential as increased melanization can result in greater heat retention. The thermal melanism hypothesis predicts that species with developmental plasticity will have darker pigmentation in colder environments, which can be an important adaptation for temperate species experiencing seasonal variation in climate. The harlequin bug (Murgantia histrionica, Hemiptera: Pentatomidae, Hahn 1834) is a widespread invasive crop pest with variable patterning where developmental plasticity in melanization could affect performance. To investigate the impact of temperature and photoperiod on melanization and size, nymphs were reared under two temperatures and two photoperiods simulating summer and fall seasons. The size and degree of melanization of adults were quantified using digital imagery. To assess the effect of coloration on the amount of heat absorption, we monitored the temperature of adults in a heating experiment. Overall, our results supported the thermal melanism hypothesis and temperature had a comparatively larger effect on coloration and size than photoperiod. When heated, the body temperature of individuals with darker pigmentation increased more relative to the ambient air temperature than individuals with lighter pigmentation. These results suggest that colder temperatures experienced late in the season can induce developmental plasticity for a phenotype that improves thermoregulation in this species. Our work highlights environmental signals and consequences for individual performance due to thermal melanism in a common invasive species, where capacity to respond to changing environments is likely contributing to its spread.

Experimental evidence for rapid genomic adaptation to a new niche in an adaptive radiation

A substantial part of biodiversity is thought to have arisen from adaptive radiations in which one lineage rapidly diversified into multiple lineages adapted to many different niches. However, selection and drift reduce genetic variation during adaptation to new niches and may thus prevent or slow down further niche shifts. We tested whether rapid adaptation is still possible from a highly derived ecotype in the adaptive radiation of threespine stickleback on the Haida Gwaii archipelago, Western Canada. In a 19-years selection experiment, we let giant stickleback from a large blackwater lake evolve in a small clearwater pond without vertebrate predators. 56 whole genomes from the experiment and 26 natural populations revealed that adaptive genomic change was rapid in many small genomic regions and encompassed 75% of the adaptive genomic change between 12,000 years old ecotypes. Adaptive genomic change was as fast as phenotypic change in defence and trophic morphology and both were largely parallel between the short-term selection experiment and long-term natural adaptive radiation. Our results show that functionally relevant standing genetic variation can persist in derived adaptive radiation members, allowing adaptive radiations to unfold very rapidly.

Mate choice, sexual selection, and endocrine-disrupting chemicals

Humans have disproportionately affected the habitat and survival of species through environmental contamination. Important among these anthropogenic influences is the proliferation of organic chemicals, some of which perturb hormone systems, the latter referred to as endocrine-disrupting chemicals (EDCs). EDCs are widespread in the environment and affect all levels of reproduction, including development of reproductive organs, hormone release and regulation through the life cycle, the development of secondary sexual characteristics, and the maturation and maintenance of adult physiology and behavior. However, what is not well-known is how the confluence of EDC actions on the manifestation of morphological and behavioral sexual traits influences mate choice, a process that requires the reciprocal evaluation of and/or acceptance of a sexual partner. Moreover, the outcomes of EDC-induced perturbations are likely to influence sexual selection; yet this has rarely been directly tested. Here, we provide background on the development and manifestation of sexual traits, reproductive competence, and the neurobiology of sexual behavior, and evidence for their perturbation by EDCs. Selection acts on individuals, with the consequences manifest in populations, and we discuss the implications for EDC contamination of these processes, and the future of species.

Evolution of life in urban environments

Our planet is an increasingly urbanized landscape, with over half of the human population residing in cities. Despite advances in urban ecology, we do not adequately understand how urbanization affects the evolution of organisms, nor how this evolution may affect ecosystems and human health. Here, we review evidence for the effects of urbanization on the evolution of microbes, plants, and animals that inhabit cities. Urbanization affects adaptive and nonadaptive evolutionary processes that shape the genetic diversity within and between populations. Rapid adaptation has facilitated the success of some native species in urban areas, but it has also allowed human pests and disease to spread more rapidly. The nascent field of urban evolution brings together efforts to understand evolution in response to environmental change while developing new hypotheses concerning adaptation to urban infrastructure and human socioeconomic activity. The next generation of research on urban evolution will provide critical insight into the importance of evolution for sustainable interactions between humans and our city environments.

Robustness and Radiation Resistance of the Pale Grass Blue Butterfly from Radioactively Contaminated Areas: A Possible Case of Adaptive Evolution

The pale grass blue butterfly, Zizeeria maha, has been used to evaluate biological impacts of the Fukushima nuclear accident in March 2011. Here, we examined the possibility that butterflies have adapted to be robust in the contaminated environment. Larvae (n = 2432) were obtained from adult butterflies (n = 20) collected from 7 localities with various contamination levels in May 2012, corresponding to the 7th generation after the accident. When the larvae were reared on non-contaminated host plant leaves from Okinawa, the normality rates of natural exposure without artificial irradiation (as an indication of robustness) were high not only in the least contaminated locality but also in the most contaminated localities. The normality rates were similarly obtained when the larvae were reared on non-contaminated leaves with external irradiation or on contaminated leaves from Fukushima to deliver internal irradiation. The normality rate of natural exposure and that of external or internal exposure were correlated, suggesting that radiation resistance (or susceptibility) likely reflects general state of health. The normality rate of external or internal exposure was divided by the relative normality rate of natural exposure, being defined as the resistance value. The resistance value was the highest in the populations of heavily contaminated localities and was inversely correlated with the distance from the Fukushima Dai-ichi nuclear power plant. These results suggest that the butterfly population might have adapted to the contaminated environment within approximately 1 year after the accident. The present study may partly explain the decrease in mortality and abnormality rates later observed in the contaminated areas.

Natural Variation in the Distribution and Abundance of Transposable Elements Across the Caenorhabditis elegans Species

Transposons are mobile DNA elements that generate both adaptive and deleterious phenotypic variation thereby driving genome evolution. For these reasons, genomes have mechanisms to regulate transposable element (TE) activity. Approximately 12–16% of the Caenorhabditis elegans genome is composed of TEs, of which the majority are likely inactive. However, most studies of TE activity have been conducted in the laboratory strain N2, which limits our knowledge of the effects of these mobile elements across natural populations. We analyzed the distribution and abundance of TEs in 208 wild C. elegans strains to better understand how transposons contribute to variation in natural populations. We identified 3,397 TEs as compared with the reference strain, of which 2,771 are novel insertions and 241 are TEs that have been excised in at least one wild strain. Likely because of their hypothesized deleterious effects, we find that TEs are found at low allele frequencies throughout the population, and we predict functional effects of TE insertions. The abundances of TEs reflect their activities, and these data allowed us to perform both genome-wide association mappings and rare variant correlations to reveal several candidate genes that impact TE regulation, including small regulatory piwi-interacting RNAs and chromatin factors. Because TE variation in natural populations could underlie phenotypic variation for organismal and behavioral traits, the transposons that we identified and their regulatory mechanisms can be used in future studies to explore the genomics of complex traits and evolutionary changes.

Two extended haplotype blocks are associated with adaptation to high altitude habitats in East African honey bees

Understanding the genetic basis of adaption is a central task in biology. Populations of the honey bee Apis mellifera that inhabit the mountain forests of East Africa differ in behavior and morphology from those inhabiting the surrounding lowland savannahs, which likely reflects adaptation to these habitats. We performed whole genome sequencing on 39 samples of highland and lowland bees from two pairs of populations to determine their evolutionary affinities and identify the genetic basis of these putative adaptations. We find that in general, levels of genetic differentiation between highland and lowland populations are very low, consistent with them being a single panmictic population. However, we identify two loci on chromosomes 7 and 9, each several hundred kilobases in length, which exhibit near fixation for different haplotypes between highland and lowland populations. The highland haplotypes at these loci are extremely rare in samples from the rest of the world. Patterns of segregation of genetic variants suggest that recombination between haplotypes at each locus is suppressed, indicating that they comprise independent structural variants. The haplotype on chromosome 7 harbors nearly all octopamine receptor genes in the honey bee genome. These have a role in learning and foraging behavior in honey bees and are strong candidates for adaptation to highland habitats. Molecular analysis of a putative breakpoint indicates that it may disrupt the coding sequence of one of these genes. Divergence between the highland and lowland haplotypes at both loci is extremely high suggesting that they are ancient balanced polymorphisms that greatly predate divergence between the extant honey bee subspecies.

Identifying the genes and genetic changes responsible for environmental adaptation is an important step towards understanding how species evolve. The honey bee Apis mellifera has adapted to a variety of habitats across its worldwide geographical distribution. Here we aim to identify the genetic basis of adaptation in honey bees living at high altitudes in the mountains of East Africa, which differ in appearance and behavior from their lowland relatives. We compare whole genome sequences from highland and lowland populations and find that, although in general they are extremely similar, there are two specific chromosomal regions (representing 1.4% of the genome) where they are strongly differentiated. These regions appear to represent structural rearrangements that are strongly correlated with altitude and contain many genes. One of these genomic regions harbors a set of octopamine receptor genes, which we hypothesize regulate differences in learning and foraging behavior between highland and lowland bees. The extremely high divergence between highland and lowland genetic variants in these regions indicates that they have an ancient origin and were likely to have been involved in environmental adaptation even before honey bees came to inhabit their current range.

Rapid Adaptation of a Polygenic Trait After a Sudden Environmental Shift

Although a number of studies have shown that natural and laboratory populations initially well adapted to their environment can evolve rapidly when conditions suddenly change, the dynamics of rapid adaptation are not well understood. Here a population genetic model of polygenic selection is analyzed to describe the short-term response of a quantitative trait after a sudden shift of the phenotypic optimum. We provide explicit analytical expressions for the timescales over which the trait mean approaches the new optimum. We find that when the effect sizes are small relative to a scaled mutation rate, small to moderate allele frequency changes occur in the short-term phase in a synergistic fashion. In contrast, selective sweeps, i.e., dramatic changes in the allele frequency, may occur provided the size of the effect is sufficiently large. Applications of our theoretical results to the relationship between QTL and selective sweep mapping and to tests of fast polygenic adaptation are discussed.

Genetic dissection of adaptive form and function in rapidly speciating cichlid fishes

Genes of major phenotypic effects and strong genetic correlations can facilitate adaptation, direct selective responses, and potentially lead to phenotypic convergence. However, the preponderance of this type of genetic architecture in repeatedly evolved adaptations remains unknown. Using hybrids between Haplochromis chilotes (thick-lipped) and Pundamilia nyererei (thin-lipped) we investigated the genetics underlying hypertrophied lips and elongated heads, traits that evolved repeatedly in cichlids. At least 25 loci of small-to-moderate and mainly additive effects were detected. Phenotypic variation in lip and head morphology was largely independent. Although several QTL overlapped for lip and head morphology traits, they were often of opposite effects. The distribution of effect signs suggests strong selection on lips. The fitness implications of several detected loci were demonstrated using a laboratory assay testing for the association between genotype and variation in foraging performance. The persistence of low fitness alleles in head morphology appears to be maintained through antagonistic pleiotropy/close linkage with positive-effect lip morphology alleles. Rather than being based on few major loci with strong positive genetic correlations, our results indicate that the evolution of the Lake Victoria thick-lipped ecomorph is the result of selection on numerous loci distributed throughout the genome.