Measuring natural selection on genotypes and phenotypes in the wild

Operationalizing cultural adaptation to climate change: contemporary examples from United States agriculture

It has been proposed that climate adaptation research can benefit from an evolutionary approach. But related empirical research is lacking. We advance the evolutionary study of climate adaptation with two case studies from contemporary United States agriculture. First, we define 'cultural adaptation to climate change' as a mechanistic process of population-level cultural change. We argue this definition enables rigorous comparisons, yields testable hypotheses from mathematical theory and distinguishes adaptive change, non-adaptive change and desirable policy outcomes. Next, we develop an operational approach to identify 'cultural adaptation to climate change' based on established empirical criteria. We apply this approach to data on crop choices and the use of cover crops between 2008 and 2021 from the United States. We find evidence that crop choices are adapting to local trends in two separate climate variables in some regions of the USA. But evidence suggests that cover cropping may be adapting more to the economic environment than climatic conditions. Further research is needed to characterize the process of cultural adaptation, particularly the routes and mechanisms of cultural transmission. Furthermore, climate adaptation policy could benefit from research on factors that differentiate regions exhibiting adaptive trends in crop choice from those that do not. This article is part of the theme issue 'Climate change adaptation needs a science of culture'.

Rural selection drives the evolution of an urban-rural cline in coat color in gray squirrels

Phenotypic differences between urban and rural populations are well-documented, but the evolutionary processes driving trait variation along urbanization gradients are often unclear. We combined spatial data on abundance, trait variation, and measurements of fitness to understand cline structure and test for natural selection on heritable coat color morphs (melanic, gray) of eastern gray squirrels (Sciurus carolinensis) along an urbanization gradient. Population surveys using remote cameras and visual counts at 76 sites along the urbanization gradient revealed a significant cline in melanism, decreasing from 48% in the city center to <5% in rural woodlands. Among 76 squirrels translocated to test for phenotypic selection, survival was lower for the melanic than gray morph in rural woodlands, whereas there was no difference in survival between color morphs in the city. These results suggest the urban-rural cline in melanism is explained by natural selection favoring the gray morph in rural woodlands combined with relaxed selection in the city. Our study illustrates how trait variation between urban and rural populations can emerge from selection primarily in rural populations rather than adaptation to novel features of the urban environment.

Allele frequencies and selection coefficients in locally adapted populations

Understanding the role of natural selection in driving evolutionary change requires accurate estimates of the strength of selection acting at the genetic level in the wild. This is challenging to achieve but may be easier in the case of populations in migration-selection balance. When two populations are at equilibrium under migration-selection balance, there exist loci whose alleles are selected different ways in the two populations. Such loci can be identified from genome sequencing by their high values of F_ST. This raises the question of what is the strength of selection on locally-adaptive alleles. To answer this question we analyse a 1-locus 2-allele model of a population distributed between two niches. We show by simulation of selected cases that the outputs from finite-population models are essentially the same as those from deterministic infinite-population models. We then derive theory for the infinite-population model showing the dependence of selection coefficients on equilibrium allele frequencies, migration rates, dominance and relative population sizes in the two niches. An Excel spreadsheet is provided for the calculation of selection coefficients and their approximate standard errors from observed values of population parameters. We illustrate our results with a worked example, with graphs showing the dependence of selection coefficients on equilibrium allele frequencies, and graphs showing how F_ST depends on the selection coefficients acting on the alleles at a locus. Given the extent of recent progress in ecological genomics, we hope our methods may help those studying migration-selection balance to quantify the advantages conferred by adaptive genes.

Turning biology to life: some reflections

This essay presents various reflections on living systems, what they are and how they evolve, prompted by editing Teleonomy in Living Systems (a special issue of the Biological Journal of the Linnean Society). Conclusions include the suggestion that the linked notions of teleonomy and agency represent fundamental properties of matter that become apparent only when organized in the way that we consider to be that of a living system. As such, they are factors that form part of the intrinsic ‘a priori’ of living systems, as they evolve in form through space and time. Biology, the science of life and living systems, needs to be ‘biological’ if it is to be anything at all. Understanding the role of teleonomy (internal, inherent goal-seeking) will always play a necessary part in this endeavour: teleonomy represents one of the fundamental properties of living systems.

Moving past the challenges and misconceptions in urban adaptation research

Although the field of urban evolutionary ecology has recently expanded, much progress has been made in identifying adaptations that arise as a result of selective pressures within these unique environments. However, as studies within urban environments have rapidly increased, researchers have recognized that there are challenges and opportunities in characterizing urban adaptation. Some of these challenges are a consequence of increased direct and indirect human influence, which compounds long-recognized issues with research on adaptive evolution more generally. In this perspective, we discuss several common research challenges to urban adaptation related to (1) methodological approaches, (2) trait-environment relationships and the natural history of organisms, (3) agents and targets of selection, and (4) habitat heterogeneity. Ignoring these challenges may lead to misconceptions and further impede our ability to draw conclusions regarding evolutionary and ecological processes in urban environments. Our goal is to first shed light on the conceptual challenges of conducting urban adaptation research to help avoid the propagation of these misconceptions. We further summarize potential strategies to move forward productively to construct a more comprehensive picture of urban adaptation, and discuss how urban environments also offer unique opportunities and applications for adaptation research.

Experimental evolution of a more restrained clutch size when filial cannibalism is prevented in burying beetles Nicrophorus vespilloides

The overproduction of offspring is commonly associated with high hatching failure and a mechanism for dispensing with surplus young. We used experimental evolution of burying beetle populations Nicrophorus vespilloides to determine causality in these correlations. We asked does eliminating the mechanism for killing “spare” offspring cause the evolution of a more restrained clutch size and consequently select for reduced hatching failure? N. vespilloides typically overproduces eggs but kills 1^st instar larvae through partial filial cannibalism during brood care. We established replicate evolving populations that either could practice filial cannibalism (Full Care) or could not, by removing parents before their young hatched (No Care). After 20+ generations of experimental evolution, we measured clutch size and hatching success. We found that No Care females produced fewer eggs than Full Care females when allowed to breed on a small corpse, a finding not explained by differences in female quality. On larger corpses, females from both populations laid similar numbers of eggs. Furthermore, hatching success was greater in the No Care populations on small corpses. Our results suggest that the adaptive overproduction of offspring depends on a mechanism for eliminating surplus young and that killing offspring, in turn, relaxes selection against hatching failure.

Avian introgression patterns are consistent with Haldane's Rule

According to Haldane’s Rule, the heterogametic sex will show the greatest fitness reduction in a hybrid cross. In birds, where sex is determined by a ZW system, female hybrids are expected to experience lower fitness compared to male hybrids. This pattern has indeed been observed in several bird groups, but it is unknown whether the generality of Haldane’s Rule also extends to the molecular level. First, given the lower fitness of female hybrids, we can expect maternally inherited loci (i.e., mitochondrial and W-linked loci) to show lower introgression rates than biparentally inherited loci (i.e., autosomal loci) in females. Second, the faster evolution of Z-linked loci compared to autosomal loci and the hemizygosity of the Z-chromosome in females might speed up the accumulation of incompatible alleles on this sex chromosome, resulting in lower introgression rates for Z-linked loci than for autosomal loci. I tested these expectations by conducting a literature review which focused on studies that directly quantified introgression rates for autosomal, sex-linked, and mitochondrial loci. Although most studies reported introgression rates in line with Haldane’s Rule, it remains important to validate these genetic patterns with estimates of hybrid fitness and supporting field observations to rule out alternative explanations. Genomic data provide exciting opportunities to obtain a more fine-grained picture of introgression rates across the genome, which can consequently be linked to ecological and behavioral observations, potentially leading to novel insights into the genetic mechanisms underpinning Haldane’s Rule.

Parallel evolution of urban-rural clines in melanism in a widespread mammal

Urbanization is the dominant trend of global land use change. The replicated nature of environmental change associated with urbanization should drive parallel evolution, yet insight into the repeatability of evolutionary processes in urban areas has been limited by a lack of multi-city studies. Here we leverage community science data on coat color in > 60,000 eastern gray squirrels (Sciurus carolinensis) across 43 North American cities to test for parallel clines in melanism, a genetically based trait associated with thermoregulation and crypsis. We show the prevalence of melanism was positively associated with urbanization as measured by impervious cover. Urban–rural clines in melanism were strongest in the largest cities with extensive forest cover and weakest or absent in cities with warmer winter temperatures, where thermal selection likely limits the prevalence of melanism. Our results suggest that novel traits can evolve in a highly repeatable manner among urban areas, modified by factors intrinsic to individual cities, including their size, land cover, and climate.

Patterns of genetic divergence and demographic history shed light on island-mainland population dynamics and melanic plumage evolution in the white-winged Fairywren

The existence of distinct traits in island versus mainland populations offers opportunities to gain insights into how eco-evolutionary processes operate under natural conditions. We used two island colonization events in the white-winged fairywren (Malurus leucopterus) to investigate the genomic and demographic origin of melanic plumage. This avian species is distributed across most of Australia, and males of the mainland subspecies (M. l. leuconotus) exhibit a blue nuptial plumage in contrast to males of two island subspecies - M. l. leucopterus on Dirk Hartog Island and M. l. edouardi on Barrow Island - that exhibit a black nuptial plumage. We used reduced-representation sequencing to explore differentiation and demographic history in this species and found clear patterns of divergence between mainland and island populations, with additional substructuring on the mainland. Divergence between the mainland and Dirk Hartog was approximately 10 times more recent than the split between the mainland and Barrow Island, supporting two independent colonizations. In both cases, estimated gene flow between the mainland and the islands was low, contributing to signals of divergence among subspecies. Our results present demographic reconstructions of mainland-island dynamics and associated plumage variation in white-winged fairywrens, with broader implications regarding our understanding of convergent evolution in insular populations.

Fitness maps to a large-effect locus in introduced stickleback populations

Mutations of small effect underlie most adaptation to new environments, but beneficial variants with large fitness effects are expected to contribute under certain conditions. Genes and genomic regions having large effects on phenotypic differences between populations are known from numerous taxa, but fitness effect sizes have rarely been estimated. We mapped fitness over a generation in an F2 intercross between a marine and a lake stickleback population introduced to a freshwater pond. A quantitative trait locus map of the number of surviving offspring per F2 female detected a single, large-effect locus near Ectodysplasin (Eda), a gene having an ancient freshwater allele causing reduced bony armor and other changes. F2 females homozygous for the freshwater allele had twice the number of surviving offspring as homozygotes for the marine allele, producing a large selection coefficient, s = 0.50 ± 0.09 SE. Correspondingly, the frequency of the freshwater allele increased from 0.50 in F2 mothers to 0.58 in surviving offspring. We compare these results to allele frequency changes at the Eda gene in an Alaskan lake population colonized by marine stickleback in the 1980s. The frequency of the freshwater Eda allele rose steadily over multiple generations and reached 95% within 20 y, yielding a similar estimate of selection, s = 0.49 ± 0.05, but a different degree of dominance. These findings are consistent with other studies suggesting strong selection on this gene (and/or linked genes) in fresh water. Selection on ancient genetic variants carried by colonizing ancestors is likely to increase the prevalence of large-effect fitness variants in adaptive evolution.

Genomic insights into historical population dynamics, local adaptation, and climate change vulnerability of the East Asian Tertiary relict Euptelea (Eupteleaceae)

The warm-temperate and subtropical climate zones of East Asia are a hotspot of plant species richness and endemism, including a noticeable number of species-poor Tertiary relict tree genera. However, little is understood about when East Asian Tertiary relict plants diversified, how they responded demographically to past environmental change, and to what extent their current genomic composition (and adaptive capacity) might mitigate the effects of global warming. Here, we obtained genomic (RAD-SNP) data for 171 samples from two extant species of Euptelea in China (24 E. pleiosperma populations) and Japan (11 E. polyandra populations) to elucidate their divergence and demographic histories, genome-wide associations with current environmental variables, and genomic vulnerability to future climate change. Our results indicate that Late Miocene changes in climate and/or sea level promoted species divergence, whereas Late Pliocene uplifting in southwest China likely fostered lineage divergence within E. pleiosperma. Its subsequent range expansion into central/east (CE) China bears genomic signatures of climate-driven selection, yet extant CE populations are predicted to be most vulnerable to future climate change. For E. polyandra, geography was the only significant predictor of genomic variation. Our findings indicate a profound impact of Late Neogene geological and climate change on the evolutionary history of Euptelea, with much stronger signals of local adaptation left in China than in Japan. This study deepens our understanding of the complex evolutionary forces that influence the distribution of genetic variation of Tertiary relict trees, and provides insights into their susceptibility to global change and potential for adaptive responses. Our results lay the groundwork for future conservation and restoration programs for Euptelea.

Selective sweeps on novel and introgressed variation shape mimicry loci in a butterfly adaptive radiation

Natural selection leaves distinct signatures in the genome that can reveal the targets and history of adaptive evolution. By analysing high-coverage genome sequence data from 4 major colour pattern loci sampled from nearly 600 individuals in 53 populations, we show pervasive selection on wing patterns in the Heliconius adaptive radiation. The strongest signatures correspond to loci with the greatest phenotypic effects, consistent with visual selection by predators, and are found in colour patterns with geographically restricted distributions. These recent sweeps are similar between co-mimics and indicate colour pattern turn-over events despite strong stabilising selection. Using simulations, we compare sweep signatures expected under classic hard sweeps with those resulting from adaptive introgression, an important aspect of mimicry evolution in Heliconius butterflies. Simulated recipient populations show a distinct 'volcano' pattern with peaks of increased genetic diversity around the selected target, characteristic of sweeps of introgressed variation and consistent with diversity patterns found in some populations. Our genomic data reveal a surprisingly dynamic history of colour pattern selection and co-evolution in this adaptive radiation.

Dynamic Evolution of Antimicrobial Peptides Underscores Trade-Offs Between Immunity and Ecological Fitness

There is a developing interest in how immune genes may function in other physiological roles, and how traditionally non-immune peptides may, in fact, be active in immune contexts. In the absence of infection, the induction of the immune response is costly, and there are well-characterized trade-offs between immune defense and fitness. The agents behind these fitness costs are less understood. Here we implicate antimicrobial peptides (AMPs) as particularly costly effectors of immunity using an evolutionary framework. We describe the independent loss of AMPs in multiple lineages of Diptera (true flies), tying these observations back to life history. We then focus on the intriguing case of the glycine-rich AMP, Diptericin, and find several instances of loss, pseudogenization, and segregating null alleles. We suggest that Diptericin may be a particularly toxic component of the Dipteran immune response lost in flies either with reduced pathogen pressure or other environmental factors. As Diptericins have recently been described to have neurological roles, these findings parallel a developing interest in AMPs as potentially harmful neuropeptides, and AMPs in other roles beyond immunity.

Unravelling the evolution of complex reproductive traits with phenotypic engineering

The social dynamics surrounding courtship, mating and parental care are complex enough when just a single male and female are involved, but for species that employ multiple strategies for achieving fertilization success, the network of interactions among rivals, allies and suitors can be utterly complicated. Such is the case in the ocellated wrasse, Symphodus ocellatus, in which males adopt one of three mating strategies. The large, colourful "nesting males" court females, defend territories and care for fertilized eggs until they hatch. The smaller "satellite males" help the nesting males court females and guard against the third morph, the "sneaker males", which sneak in when a nesting male is spawning with a female and surreptitiously release sperm. Sneaker males perform no courtship displays nor defend territories, so their reproductive investment is devoted entirely to sperm production. And these alternative male strategies work: 100% of nests contain some eggs fertilized by sneaker and satellite males, despite the fact that parental care is solely the responsibility of nesting males In this issue of Molecular Ecology, work to untangle the proximate mechanisms regulating the reproductive physiology of nesting males and their behaviour towards other males, which impacts the entire social network. Moreover, they describe how variation in neuroendocrine regulation can give rise to variation in reproductive traits, upon which sexual selection can act.

The Evolutionary History of Nebraska Deer Mice: Local Adaptation in the Face of Strong Gene Flow

The interplay of gene flow, genetic drift, and local selective pressure is a dynamic process that has been well studied from a theoretical perspective over the last century. Wright and Haldane laid the foundation for expectations under an island-continent model, demonstrating that an island-specific beneficial allele may be maintained locally if the selection coefficient is larger than the rate of migration of the ancestral allele from the continent. Subsequent extensions of this model have provided considerably more insight. Yet, connecting theoretical results with empirical data has proven challenging, owing to a lack of information on the relationship between genotype, phenotype, and fitness. Here, we examine the demographic and selective history of deer mice in and around the Nebraska Sand Hills, a system in which variation at the Agouti locus affects cryptic coloration that in turn affects the survival of mice in their local habitat. We first genotyped 250 individuals from 11 sites along a transect spanning the Sand Hills at 660,000 single nucleotide polymorphisms across the genome. Using these genomic data, we found that deer mice first colonized the Sand Hills following the last glacial period. Subsequent high rates of gene flow have served to homogenize the majority of the genome between populations on and off the Sand Hills, with the exception of the Agouti pigmentation locus. Furthermore, mutations at this locus are strongly associated with the pigment traits that are strongly correlated with local soil coloration and thus responsible for cryptic coloration.

Experimentally induced variation in neuroendocrine processes affects male reproductive behaviour, sperm characteristics and social interactions

While extensive research has focused on how social interactions evolve, the fitness consequences of the neuroendocrine mechanisms underlying these interactions have rarely been documented, especially in the wild. Here, we measure how the neuroendocrine mechanisms underlying male behaviour affect mating success and sperm competition in the ocellated wrasse (Symphodus ocellatus). In this species, males exhibit three alternative reproductive types. "Nesting males" provide parental care, defend territories and form cooperative associations with unrelated "satellites," who cheat by sneaking fertilizations but help by reducing sperm competition from "sneakers" who do not cooperate or provide care. To measure the fitness consequences of the mechanisms underlying these social interactions, we used "phenotypic engineering" that involved administering an androgen receptor antagonist (flutamide) to wild, free-living fish. Nesting males treated with flutamide shifted their aggression from sneakers to satellite males and experienced decreased submissiveness by sneaker males (which correlated with decreased nesting male mating success). The preoptic area (POA), a region controlling male reproductive behaviours, exhibited dramatic down-regulation of androgen receptor (AR) and vasotocin 1a receptor (V1aR) mRNA following experimental manipulation of androgen signalling. We did not find a direct effect of the manipulation on male mating success, paternity or larval production. However, variation in neuroendocrine mechanisms generated by the experimental manipulation was significantly correlated with changes in behaviour and mating success: V1aR expression was negatively correlated with satellite-directed aggression, and expression of its ligand arginine vasotocin (AVT) was positively correlated with courtship and mating success, thus revealing the potential for sexual selection on these mechanisms.

Fine-scale geographic patterns of gene flow and reproductive character displacement in Drosophila subquinaria and Drosophila recens

When two species are incompletely isolated, strengthening premating isolation barriers in response to the production of low fitness hybrids may complete the speciation process. Here, we use the sister species Drosophila subquinaria and Drosophila recens to study the conditions under which this reinforcement of species boundaries occurs in natural populations. We first extend the region of known sympatry between these species, and then we conduct a fine-scale geographic survey of mate discrimination coupled with estimates of gene flow within and admixture between species. Within D. subquinaria, reinforcement is extremely effective: we find variation in mate discrimination both against D. recens males and against conspecific allopatric males on the scale of a few kilometres and in the face of gene flow both from conspecific populations and introgression from D. recens. In D. recens, we do not find evidence for increased mate discrimination in sympatry, even where D. recens is rare, consistent with substantial gene flow throughout the species' range. Finally, we find that introgression between species is asymmetric, with more from D. recens into D. subquinaria than vice versa. Within each species, admixture is highest in the geographic region where it is rare relative to the other species, suggesting that when hybrids are produced they are of low fitness. In sum, reinforcement within D. subquinaria is effective at maintaining species boundaries, but even when reinforcing selection is strong it may not always result in a pattern of strong reproductive character displacement due to variation in the frequency of hybridization and gene flow from neighbouring populations.

New diagnostic SNP molecular markers for the Mytilus species complex

The development of diagnostic markers has been a long-standing interest of population geneticists as it allows clarification of taxonomic uncertainties. Historically, there has been much debate on the taxonomic status of species belonging to the Mytilus species complex (M. edulis, M. galloprovincialis and M. trossulus), and whether they are discrete species. We analysed reference pure specimens of M. edulis, M. galloprovincialis and M. trossulus, using Restriction site associated DNA (RAD) sequencing and identified over 6,000 SNP markers separating the three species unambiguously. We developed a panel of diagnostic SNP markers for the genotyping of Mytilus species complex as well as the identification of hybrids and interspecies introgression events in Mytilus species. We validated a panel of twelve diagnostic SNP markers which can be used for species genotyping. Being able to accurately identify species and hybrids within the Mytilus species complex is important for the selective mussel stock management, the exclusion of invasive species, basic physiology and bio-diversity studies.

Viral Fitness Correlates with the Magnitude and Direction of the Perturbation Induced in the Host's Transcriptome: The Tobacco Etch Potyvirus-Tobacco Case Study

Determining the fitness of viral genotypes has become a standard practice in virology as it is essential to evaluate their evolutionary potential. Darwinian fitness, defined as the advantage of a given genotype with respect to a reference one, is a complex property that captures, in a single figure, differences in performance at every stage of viral infection. To what extent does viral fitness result from specific molecular interactions with host factors and regulatory networks during infection? Can we identify host genes in functional classes whose expression depends on viral fitness? Here, we compared the transcriptomes of tobacco plants infected with seven genotypes of tobacco etch potyvirus that differ in fitness. We found that the larger the fitness differences among genotypes, the more dissimilar the transcriptomic profiles are. Consistently, two different mutations, one in the viral RNA polymerase and another in the viral suppressor of RNA silencing, resulted in significantly similar gene expression profiles. Moreover, we identified host genes whose expression showed a significant correlation, positive or negative, with the virus' fitness. Differentially expressed genes which were positively correlated with viral fitness activate hormone- and RNA silencing-mediated pathways of plant defense. In contrast, those that were negatively correlated with fitness affect metabolism, reducing growth, and development. Overall, these results reveal the high information content of viral fitness and suggest its potential use to predict differences in genomic profiles of infected hosts.

Contemporary evolution of a Lepidopteran species, Heliothis virescens, in response to modern agricultural practices

Adaptation to human-induced environmental change has the potential to profoundly influence the genomic architecture of affected species. This is particularly true in agricultural ecosystems, where anthropogenic selection pressure is strong. Heliothis virescens primarily feeds on cotton in its larval stages, and US populations have been declining since the widespread planting of transgenic cotton, which endogenously expresses proteins derived from Bacillus thuringiensis (Bt). No physiological adaptation to Bt toxin has been found in the field, so adaptation in this altered environment could involve (i) shifts in host plant selection mechanisms to avoid cotton, (ii) changes in detoxification mechanisms required for cotton-feeding vs. feeding on other hosts or (iii) loss of resistance to previously used management practices including insecticides. Here, we begin to address whether such changes occurred in H. virescens populations between 1997 and 2012, as Bt-cotton cultivation spread through the agricultural landscape. For our study, we produced an H. virescens genome assembly and used this in concert with a ddRAD-seq-enabled genome scan to identify loci with significant allele frequency changes over the 15-year period. Genetic changes at a previously described H. virescens insecticide target of selection were detectable in our genome scan and increased our confidence in this methodology. Additional loci were also detected as being under selection, and we quantified the selection strength required to elicit observed allele frequency changes at each locus. Potential contributions of genes near loci under selection to adaptive phenotypes in the H. virescens cotton system are discussed.

Niche construction, sources of selection and trait coevolution

Organisms modify and choose components of their local environments. This 'niche construction' can alter ecological processes, modify natural selection and contribute to inheritance through ecological legacies. Here, we propose that niche construction initiates and modifies the selection directly affecting the constructor, and on other species, in an orderly, directed and sustained manner. By dependably generating specific environmental states, niche construction co-directs adaptive evolution by imposing a consistent statistical bias on selection. We illustrate how niche construction can generate this evolutionary bias by comparing it with artificial selection. We suggest that it occupies the middle ground between artificial and natural selection. We show how the perspective leads to testable predictions related to: (i) reduced variance in measures of responses to natural selection in the wild; (ii) multiple trait coevolution, including the evolution of sequences of traits and patterns of parallel evolution; and (iii) a positive association between niche construction and biodiversity. More generally, we submit that evolutionary biology would benefit from greater attention to the diverse properties of all sources of selection.

Living in an adaptive world: Genomic dissection of the genus Homo and its immune response

More than a decade after the sequencing of the human genome, a deluge of genome-wide population data are generating a portrait of human genetic diversity at an unprecedented level of resolution. Genomic studies have provided new insight into the demographic and adaptive history of our species, Homo sapiens, including its interbreeding with other hominins, such as Neanderthals, and the ways in which natural selection, in its various guises, has shaped genome diversity. These studies, combined with functional genomic approaches, such as the mapping of expression quantitative trait loci, have helped to identify genes, functions, and mechanisms of prime importance for host survival and involved in phenotypic variation and differences in disease risk. This review summarizes new findings in this rapidly developing field, focusing on the human immune response. We discuss the importance of defining the genetic and evolutionary determinants driving immune response variation, and highlight the added value of population genomic approaches in settings relevant to immunity and infection.

Survival by genotype: patterns at Mc1r are not black and white at the White Sands ecotone

Measuring links among genotype, phenotype and survival in the wild has long been a focus of studies of adaptation. We conducted a 4-year capture-recapture study to measure survival by genotype and phenotype in the Southwestern Fence Lizard (Sceloporus cowlesi) at the White Sands ecotone (transition area between white sands and dark soil habitats). We report several unanticipated findings. First, in contrast with previous work showing that cryptic blanched coloration in S. cowlesi from the heart of the dunes is associated with mutations in the melanocortin-1 receptor gene (Mc1r), ecotonal S. cowlesi showed minimal association between colour phenotype and Mc1r genotype. Second, the frequency of the derived Mc1r allele in ecotonal S. cowlesi appeared to decrease over time. Third, our capture-recapture data revealed a lower survival rate for S. cowlesi individuals with the derived Mc1r allele. Thus, our results suggest that selection at the ecotone may have favoured the wild-type allele in recent years. Even in a system where a genotype-phenotype association appeared to be black and white, our study suggests that additional factors - including phenotypic plasticity, epistasis, pleiotropy and gene flow - may play important roles at the White Sands ecotone. Our study highlights the importance of linking molecular, genomic and organismal approaches for understanding adaptation in the wild. Furthermore, our findings indicate that dynamics of natural selection can be particularly complex in transitional habitats like ecotones and emphasize the need for future research that examines the patterns of ongoing selection in other ecological 'grey' zones.

Bud phenology and growth are subject to divergent selection across a latitudinal gradient in Populus angustifolia and impact adaptation across the distributional range and associated arthropods

Temperate forest tree species that span large geographical areas and climatic gradients often have high levels of genetic variation. Such species are ideal for testing how neutral demographic factors and climate‐driven selection structure genetic variation within species, and how this genetic variation can affect ecological communities. Here, we quantified genetic variation in vegetative phenology and growth traits in narrowleaf cottonwood, Populus angustifolia, using three common gardens planted with genotypes originating from source populations spanning the species' range along the Rocky Mountains of North America (ca. 1700 km). We present three main findings. First, we found strong evidence of divergent selection (Q_ST > F_ST) on fall phenology (bud set) with adaptive consequences for frost avoidance. We also found evidence for selection on bud flush duration, tree height, and basal diameter, resulting in population differentiation. Second, we found strong associations with climate variables that were strongly correlated with latitude of origin. More strongly differentiated traits also showed stronger climate correlations, which emphasizes the role that climate has played in divergent selection throughout the range. We found population × garden interaction effects; for some traits, this accounted for more of the variance than either factor alone. Tree height was influenced by the difference in climate of the source and garden locations and declined with increasing transfer distance. Third, growth traits were correlated with dependent arthropod community diversity metrics. Synthesis. Overall, we conclude that climate has influenced genetic variation and structure in phenology and growth traits and leads to local adaptation in P. angustifolia, which can then impact dependent arthropod species. Importantly, relocation of genotypes far northward or southward often resulted in poor growth, likely due to a phenological mismatch with photoperiod, the proximate cue for fall growth cessation. Genotypes moved too far southward suffer from early growth cessation, whereas those moved too far northward are prone to fall frost and winter dieback. In the face of current and forecasted climate change, habitat restoration, forestry, and tree breeding efforts should utilize these findings to better match latitudinal and climatic source environments with management locations for optimal future outcomes.

Culturomics as a data playground for tests of selection: Mathematical approaches to detecting selection in word use

In biological evolution traits may rise and fall in frequency due to genetic drift, where variant frequencies change by chance, or by selection where advantageous variants will rise in frequency. The neutral model of evolution, first developed by Kimura in the 1960s, has become the standard against which selection is detected. While the balance between these two important forces - drift and selection - has been well established in biology there are other domains where the contribution of these processes is still coming together. Although the idea of natural selection has been applied to the cultural domain since the time of Darwin, it has proven more challenging to positively identify cultural traits under selection both because of a lack of established tests for selection and a lack of large cultural data sets. However, in recent years with the accumulation of large cultural data sets many cultural features from pre-history pottery to modern baby names have been shown to evolve according to the neutral theory. But there is accumulating empirical evidence from cultural processes suggesting that the neutral theory alone cannot account for all features of the data. As such, there has been a renewed interest in determining whether there is selection amidst drift. Here we analyze a subset English word frequencies, and determine whether frequency change reveals processes of selection. Inspired by the Moran and Wright-Fisher models in population genetics, we developed a neutral model of word frequency variation to assess when linguistic data appears to depart from neutral evolution. As such, our model represents a possible "test for selection" in the linguistic domain. We explore how the distribution of word use has changed for sets of words in English for more than 100 years (1901-2008) as expressed in vocabulary usage in published books, made available by Google Ngram. When comparing empirical word frequency changes to our neutral model we find pervasive and systematic departures from neutrality.

Natural selection in a postglacial range expansion: the case of the colour cline in the European barn owl

Gradients of variation--or clines--have always intrigued biologists. Classically, they have been interpreted as the outcomes of antagonistic interactions between selection and gene flow. Alternatively, clines may also establish neutrally with isolation by distance (IBD) or secondary contact between previously isolated populations. The relative importance of natural selection and these two neutral processes in the establishment of clinal variation can be tested by comparing genetic differentiation at neutral genetic markers and at the studied trait. A third neutral process, surfing of a newly arisen mutation during the colonization of a new habitat, is more difficult to test. Here, we designed a spatially explicit approximate Bayesian computation (ABC) simulation framework to evaluate whether the strong cline in the genetically based reddish coloration observed in the European barn owl (Tyto alba) arose as a by-product of a range expansion or whether selection has to be invoked to explain this colour cline, for which we have previously ruled out the actions of IBD or secondary contact. Using ABC simulations and genetic data on 390 individuals from 20 locations genotyped at 22 microsatellites loci, we first determined how barn owls colonized Europe after the last glaciation. Using these results in new simulations on the evolution of the colour phenotype, and assuming various genetic architectures for the colour trait, we demonstrate that the observed colour cline cannot be due to the surfing of a neutral mutation. Taking advantage of spatially explicit ABC, which proved to be a powerful method to disentangle the respective roles of selection and drift in range expansions, we conclude that the formation of the colour cline observed in the barn owl must be due to natural selection.

Adaptations between ecotypes and along environmental gradients in Panicum virgatum

Determining the patterns and mechanisms of natural selection in the wild is of fundamental importance to understanding the differentiation of populations and the evolution of new species. However, it is often unknown the extent to which adaptive genetic variation is distributed among ecotypes between distinct habitats versus along large-scale geographic environmental gradients, such as those that track latitude. Classic studies of selection in the wild in switchgrass, Panicum virgatum, tested for adaptation at both of these levels of natural variation. Here we review what these field experiments and modern agronomic field trials have taught us about natural variation and selection at both the ecotype and environmental gradient levels in P. virgatum. With recent genome sequencing efforts in P. virgatum, it is poised to become an excellent system for understanding the adaptation of grassland species across the eastern half of North America. The identification of genetic loci involved in different types of adaptations will help to understand the evolutionary mechanisms of diversification within P. virgatum and provide useful information for the breeding of high-yielding cultivars for different ecoregions.

Signatures of rapid evolution in urban and rural transcriptomes of white-footed mice (Peromyscus leucopus) in the New York metropolitan area

Urbanization is a major cause of ecological degradation around the world, and human settlement in large cities is accelerating. New York City (NYC) is one of the oldest and most urbanized cities in North America, but still maintains 20% vegetation cover and substantial populations of some native wildlife. The white-footed mouse, Peromyscusleucopus, is a common resident of NYC's forest fragments and an emerging model system for examining the evolutionary consequences of urbanization. In this study, we developed transcriptomic resources for urban P. leucopus to examine evolutionary changes in protein-coding regions for an exemplar "urban adapter." We used Roche 454 GS FLX+ high throughput sequencing to derive transcriptomes from multiple tissues from individuals across both urban and rural populations. From these data, we identified 31,015 SNPs and several candidate genes potentially experiencing positive selection in urban populations of P. leucopus. These candidate genes are involved in xenobiotic metabolism, innate immune response, demethylation activity, and other important biological phenomena in novel urban environments. This study is one of the first to report candidate genes exhibiting signatures of directional selection in divergent urban ecosystems.

Genetic signature of adaptive peak shift in threespine stickleback

Transition of an evolving population to a new adaptive optimum is predicted to leave a signature in the distribution of effect sizes of fixed mutations. If they affect many traits (are pleiotropic), large effect mutations should contribute more when a population evolves to a farther adaptive peak than to a nearer peak. We tested this prediction in wild threespine stickleback fish (Gasterosteus aculeatus) by comparing the estimated frequency of large effect genetic changes underlying evolution as the same ancestor adapted to two lake types since the end of the ice age. A higher frequency of large effect genetic changes (quantitative trait loci) contributed to adaptive evolution in populations that adapted to lakes representing a more distant optimum than to lakes in which the optimum phenotype was nearer to the ancestral state. Our results also indicate that pleiotropy, not just optimum overshoot, contributes to this difference. These results suggest that a series of adaptive improvements to a new environment leaves a detectable mark in the genome of wild populations. Although not all assumptions of the theory are likely met in natural systems, the prediction may be robust enough to the complexities of natural environments to be useful when forecasting adaptive responses to large environmental changes.

Challenges in estimating insecticide selection pressures from mosquito field data

Insecticide resistance has the potential to compromise the enormous effort put into the control of dengue and malaria vector populations. It is therefore important to quantify the amount of selection acting on resistance alleles, their contributions to fitness in heterozygotes (dominance) and their initial frequencies, as a means to predict the rate of spread of resistance in natural populations. We investigate practical problems of obtaining such estimates, with particular emphasis on Mexican populations of the dengue vector Aedes aegypti. Selection and dominance coefficients can be estimated by fitting genetic models to field data using maximum likelihood (ML) methodology. This methodology, although widely used, makes many assumptions so we investigated how well such models perform when data are sparse or when spatial and temporal heterogeneity occur. As expected, ML methodologies reliably estimated selection and dominance coefficients under idealised conditions but it was difficult to recover the true values when datasets were sparse during the time that resistance alleles increased in frequency, or when spatial and temporal heterogeneity occurred. We analysed published data on pyrethroid resistance in Mexico that consists of the frequency of a Ile1,016 mutation. The estimates for selection coefficient and initial allele frequency on the field dataset were in the expected range, dominance coefficient points to incomplete dominance as observed in the laboratory, although these estimates are accompanied by strong caveats about possible impact of spatial and temporal heterogeneity in selection.

Repeated modification of early limb morphogenesis programmes underlies the convergence of relative limb length in Anolis lizards

The independent evolution of similar morphologies has long been a subject of considerable interest to biologists. Does phenotypic convergence reflect the primacy of natural selection, or does development set the course of evolution by channelling variation in certain directions? Here, we examine the ontogenetic origins of relative limb length variation among Anolis lizard habitat specialists to address whether convergent phenotypes have arisen through convergent developmental trajectories. Despite the numerous developmental processes that could potentially contribute to variation in adult limb length, our analyses reveal that, in Anolis lizards, such variation is repeatedly the result of changes occurring very early in development, prior to formation of the cartilaginous long bone anlagen.