Lack of parallel genetic patterns underlying the repeated ecological divergence of beach and stream‐spawning kokanee salmon

Recent progress in methods for detecting adaptive population divergence in situ shows promise for elucidating the conditions under which selection acts to generate intraspecific diversity. Rapid ecological diversification is common in fishes; however, the role of phenotypic plasticity and adaptation to local environments is poorly understood. It is now possible to investigate genetic patterns to make inferences regarding phenotypic traits under selection and possible mechanisms underlying ecotype divergence, particularly where similar novel phenotypes have arisen in multiple independent populations. Here, we employed a bottom‐up approach to test for signatures of directional selection associated with divergence of beach‐ and stream‐spawning kokanee, the obligate freshwater form of sockeye salmon (Oncorhynchus nerka). Beach‐ and stream‐spawners co‐exist in many post‐glacial lakes and exhibit distinct reproductive behaviours, life‐history traits and spawning habitat preferences. Replicate ecotype pairs across five lakes in British Columbia, Canada were genotyped at 57 expressed sequence tag‐linked and anonymous microsatellite loci identified in a previous genome scan. Fifteen loci exhibited signatures of directional selection (high F_ST outliers), four of which were identified in multiple lakes. However, the lack of parallel genetic patterns across all lakes may be a result of: 1) an inability to detect loci truly under selection; 2) alternative genetic pathways underlying ecotype divergence in this system; and/or 3) phenotypic plasticity playing a formative role in driving kokanee spawning habitat differences. Gene annotations for detected outliers suggest pathogen resistance and energy metabolism as potential mechanisms contributing to the divergence of beach‐ and stream‐spawning kokanee, but further study is required.     

Selection on life‐history traits and genetic population divergence in rotifers

A combination of founder effects and local adaptation – the Monopolization hypothesis – has been proposed to reconcile the strong population differentiation of zooplankton dwelling in ponds and lakes and their high dispersal abilities. The role genetic drift plays in genetic differentiation of zooplankton is well documented, but the impact of natural selection has received less attention. Here, we compare differentiation in neutral genetic markers (F_ST) and in quantitative traits (Q_ST) in six natural populations of the rotifer Brachionus plicatilis to assess the importance of natural selection in explaining genetic differentiation of life‐history traits. Five life‐history traits were measured in four temperature × salinity combinations in common‐garden experiments. Population differentiation for neutral genetic markers – 11 microsatellite loci – was very high (F_ST = 0.482). Differentiation in life‐history traits was higher in traits related to sexual reproduction than in those related to asexual reproduction. Q_ST values for diapausing egg production (a trait related to sexual reproduction) were higher than their corresponding F_ST in some pairs of populations. Our results indicate the importance of divergent natural selection in these populations and suggest local adaptation to the unpredictability of B. plicatilis habitats.     

A three‐component system incorporating Ppd‐D1, copy number variation at Ppd‐B1, and numerous small‐effect quantitative trait loci facilitates adaptation of heading time in winter wheat cultivars of worldwide origin

The broad adaptability of heading time has contributed to the global success of wheat in a diverse array of climatic conditions. Here, we investigated the genetic architecture underlying heading time in a large panel of 1,110 winter wheat cultivars of worldwide origin. Genome‐wide association mapping, in combination with the analysis of major phenology loci, revealed a three‐component system that facilitates the adaptation of heading time in winter wheat. The photoperiod sensitivity locus Ppd‐D1 was found to account for almost half of the genotypic variance in this panel and can advance or delay heading by many days. In addition, copy number variation at Ppd‐B1 was the second most important source of variation in heading, explaining 8.3% of the genotypic variance. Results from association mapping and genomic prediction indicated that the remaining variation is attributed to numerous small‐effect quantitative trait loci that facilitate fine‐tuning of heading to the local climatic conditions. Collectively, our results underpin the importance of the two Ppd‐1 loci for the adaptation of heading time in winter wheat and illustrate how the three components have been exploited for wheat breeding globally.     

Population‐genomic approach reveals adaptive floral divergence in discrete populations of a hawk moth‐pollinated violet

Local adaptation to contrasting biotic or abiotic environments is an important evolutionary step that presumably precedes floral diversification at the species level, yet few studies have demonstrated the adaptive nature of intraspecific floral divergence in wild plant populations. We combine a population‐genomic approach with phenotypic information on floral traits to examine whether the differentiation in metric floral traits exhibited by 14 populations of the southern Spanish hawk moth‐pollinated violet Viola cazorlensis reflects adaptive divergence. Screening of many amplified fragment length polymorphism (AFLP) loci using a multiple‐marker‐based neutrality test identified nine outlier loci (2.6% of the total) that departed from neutral expectations and were potentially under selection. Generalized analysis of molecular variance revealed significant relationships between genetic distance and population divergence in three floral traits when genetic distance was based on outlier loci, but not when it was based on neutral ones. Population means of floral traits were closely correlated with population scores on the first principal coordinate axis of the genetic distance matrix using outlier loci, and with the allelic frequencies of four of the outlier loci. Results strongly support the adaptive nature of intraspecific floral divergence exhibited by V. cazorlensis and illustrate the potential of genome scans to identify instances of adaptive divergence when used in combination with phenotypic information.     

A highly pleiotropic amino acid polymorphism in the Drosophila insulin receptor contributes to life‐history adaptation

Finding the specific nucleotides that underlie adaptive variation is a major goal in evolutionary biology, but polygenic traits pose a challenge because the complex genotype–phenotype relationship can obscure the effects of individual alleles. However, natural selection working in large wild populations can shift allele frequencies and indicate functional regions of the genome. Previously, we showed that the two most common alleles of a complex amino acid insertion–deletion polymorphism in the Drosophila insulin receptor show independent, parallel clines in frequency across the North American and Australian continents. Here, we report that the cline is stable over at least a five‐year period and that the polymorphism also demonstrates temporal shifts in allele frequency concurrent with seasonal change. We tested the alleles for effects on levels of insulin signaling, fecundity, development time, body size, stress tolerance, and life span. We find that the alleles are associated with predictable differences in these traits, consistent with patterns of Drosophila life‐history variation across geography that likely reflect adaptation to the heterogeneous climatic environment. These results implicate insulin signaling as a major mediator of life‐history adaptation in Drosophila, and suggest that life‐history trade‐offs can be explained by extensive pleiotropy at a single locus.     

Candidate loci reveal genetic differentiation between temporally divergent migratory runs of Chinook salmon (Oncorhynchus tshawytscha)

Local adaptation is a dynamic process driven by selection that can vary both in space and time. One important temporal adaptation for migratory animals is the time at which individuals return to breeding sites. Chinook salmon (Oncorhynchus tshawytscha) are excellent subjects for studying the genetic basis of temporal adaptation because their high seasonal homing fidelity promotes reproductive isolation leading to the formation of local populations across diverse environments. We tested for adaptive genetic differentiation between seasonal runs of Chinook salmon using two candidate loci; the circadian rhythm gene, OtsClock1b, and Ots515NWFSC, a microsatellite locus showing sequence identity to three salmonid genes central to reproductive development. We found significant evidence for two genetically distinct migratory runs in the Feather River, California (OtsClock1b: F(ST)=0.042, P=0.02; Ots515NWFSC: F(ST)=0.058, P=0.003). In contrast, the fall and threatened spring runs are genetically homogenous based on neutral microsatellite data (F(ST)=-0.0002). Similarly, two temporally divergent migratory runs of Chinook salmon from New Zealand are genetically differentiated based on polymorphisms in the candidate loci (OtsClock1b: F(ST)=0.083, P-value=0.001; Ots515NWFSC: F(ST)=0.095, P-value=0.000). We used an individual-based assignment method to confirm that these recently diverged populations originated from a single source in California. Tests for selective neutrality indicate that OtsClock1b and Ots515NWFSC exhibit substantial departures from neutral expectations in both systems. The large F(ST )estimates could therefore be the result of directional selection. Evidence presented here suggests that OtsClock1b and Ots515NWFSC may influence migration and spawning timing of Chinook salmon in these river systems.     

Explaining the sawtooth: latitudinal periodicity in a circadian gene correlates with shifts in generation number

Many temperate insects take advantage of longer growing seasons at lower latitudes by increasing their generation number or voltinism. In some insects, development time abruptly decreases when additional generations are fit into the season. Consequently, latitudinal ‘sawtooth’ clines associated with shifts in voltinism are seen for phenotypes correlated with development time, like body size. However, latitudinal variation in voltinism has not been linked to genetic variation at specific loci. Here, we show a pattern in allele frequency among voltinism ecotypes of the European corn borer moth (Ostrinia nubilalis) that is reminiscent of a sawtooth cline. We characterized 145 autosomal and sex‐linked SNPs and found that period, a circadian gene that is genetically linked to a major QTL determining variation in post‐diapause development time, shows cyclical variation between voltinism ecotypes. Allele frequencies at an unlinked circadian clock gene cryptochrome1 were correlated with period. These results suggest that selection on development time to ‘fit’ complete life cycles into a latitudinally varying growing season produces oscillations in alleles associated with voltinism, primarily through changes at loci underlying the duration of transitions between diapause and other life history phases. Correlations among clock loci suggest possible coupling between the circadian clock and the circannual rhythms for synchronizing seasonal life history. We anticipate that latitudinal oscillations in allele frequency will represent signatures of adaptation to seasonal environments in other insects and may be critical to understanding the ecological and evolutionary consequences of variable environments, including response to global climate change.     

Inbreeding alters context‐dependent reproductive effort and immunity in male crickets

Infection can cause hosts to drastically alter their investment in key life‐history traits of reproduction and defence. Infected individuals are expected to increase investment in defence (e.g., by increasing immune function) and, due to trade‐offs, investment in other traits (e.g., current reproduction) should decrease. However, the terminal investment hypothesis postulates that decreased lifespan due to infection and the associated reduction in the expectation for future offspring will favour increased investment towards current reproduction. Variation in intrinsic condition will likely influence shifts in reproductive investment post‐infection, but this is often not considered in such assessments. For example, the extent of inbreeding can significantly impact an individual's lifetime fitness and may influence its reproductive behaviour following a threat of infection. Here, we investigated the effects of inbreeding status on an individual's reproductive investment upon infection, including the propensity to terminally invest. Male crickets (Gryllodes sigillatus) from four genetically distinct inbred lines and one outbred line were subjected to a treatment from an increasing spectrum of simulated infection cue intensities, using heat‐killed bacteria. We then measured reproductive effort (calling effort), survival and immune function (antibacterial activity, circulating haemocytes and haemocyte microaggregations). Inbred and outbred males diverged in how they responded to a low‐dose infection cue: relative to unmanipulated males, outbred males decreased calling effort, whereas inbred males increased calling effort. Moreover, we found that inbred males exhibited higher antibacterial activity and numbers of circulating haemocytes compared with outbred males. These results suggest that an individual's inbreeding status may have consequences for context‐dependent shifts in reproductive strategies, such as those triggered by infection.     

Heritability of and strong single gene (Pgi) effects on life‐history traits in the Glanville fritillary butterfly

We estimated broad‐sense heritabilities (H²) of 13 female and seven male life‐history traits of the Glanville fritillary butterfly (Melitaea cinxia) under semi‐natural conditions in a large outdoor population cage. The analysis was based on full‐sib families collected as young larvae in the field and reared under common garden conditions. We found significant genetic variance in female lifespan, fecundity, number of matings and host‐plant preference as well as in male body mass and mobility. Apart from host‐plant preference, female traits that were more strongly correlated with lifetime reproductive success (LRS; measured as total number of eggs laid) had higher H². LRS itself exhibited significant heritability. Host‐plant preference had very high H², consistent with a previously reported genetically determined geographical cline in host‐plant preference in the study area. Lifespan and egg hatching rate were significantly associated with a SNP in the coding region of the Pgi gene, for which there is previous evidence for balancing selection. Selection on Pgi, which furthermore shows spatial and temporal variation, may maintain genetic variance in fitness‐related life‐history traits. In contrast, we found no strong evidence for life‐history trade‐offs.     

The influence of population structure on gene expression and flowering time variation in the ubiquitous weed Capsella bursa‐pastoris (Brassicaceae)

Population structure is a potential problem when testing for adaptive phenotypic differences among populations. The observed phenotypic differences among populations can simply be due to genetic drift, and if the genetic distance between them is not considered, the differentiation may be falsely interpreted as adaptive. Conversely, adaptive and demographic processes might have been tightly associated and correcting for the population structure may lead to false negatives. Here, we evaluated this problem in the cosmopolitan weed Capsella bursa‐pastoris. We used RNA‐Seq to analyse gene expression differences among 24 accessions, which belonged to a much larger group that had been previously characterized for flowering time and circadian rhythm and were genotyped using genotyping‐by‐sequencing (GBS) technique. We found that clustering of accessions for gene expression retrieved the same three clusters that were obtained with GBS data previously, namely Europe, the Middle East and Asia. Moreover, the three groups were also differentiated for both flowering time and circadian rhythm variation. Correction for population genetic structure when analysing differential gene expression analysis removed all differences among the three groups. This may suggest that most differences are neutral and simply reflect population history. However, geographical variation in flowering time and circadian rhythm indicated that the distribution of adaptive traits might be confounded by population structure. To bypass this confounding effect, we compared gene expression differentiation between flowering ecotypes within the genetic groups. Among the differentially expressed genes, FLOWERING LOCUS C was the strongest candidate for local adaptation in regulation of flowering time.     

Comparative tests of the species‐genetic diversity correlation at neutral and nonneutral loci in four species of stream insect

A fundamental question linking population genetics and community ecology is how adaptive processes (e.g., natural selection) and neutral processes (e.g., drift‐migration equilibrium) underpin the species‐genetic diversity correlation (SGDC). Here, we combine genome scans and outlier loci detection with community analysis to separately test for neutral and nonneutral SGDCs in four species of stream insect. We sampled 60 localities in Japan and examined the relationships among population AFLP band richness (Br), taxon richness of the total community (S) and of the trophic guild (S_tr), and 15 habitat parameters that could potentially drive adaptation and influence richness. Neutral Br was positively correlated with S only in the dominant species of these communities, suggesting Br may be constrained when intraspecific competition is pronounced. Nonneutral Br was correlated with S_tr in a species restricted to high elevations where habitat heterogeneity was highest. Community distance and genetic distance (β‐SGDC) was correlated in two of the four species at both neutral and nonneutral loci. Distance‐based redundancy analysis found geographic isolation and elevation to drive divergence of both communities and populations. This suggests that both neutral and adaptive divergence occurred through the shared influences of geographic isolation and local adaptation at the two levels of diversity.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION DURING SPECIATION‐WITH‐GENE‐FLOW IN LAKE WHITEFISH SPECIES PAIRS ASSESSED BY RAD SEQUENCING

During speciation‐with‐gene‐flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome‐wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome‐wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation‐with‐gene‐flow mechanism.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.     

Life‐history tradeoffs revealed by seasonal declines in reproductive traits of Arctic‐breeding shorebirds

Seasonal declines in breeding performance are widespread in wild animals, resulting from temporal changes in environmental conditions or from individual variation. Seasonal declines might drive selection for early breeding, with implications for other stages of the annual cycle. Alternatively, selection on the phenology of nonbreeding stages could constrain timing of the breeding season and lead to seasonal changes in reproductive performance. We studied 25 taxa of migratory shorebirds (including five subspecies) at 16 arctic sites in Russia, Alaska, and Canada. We investigated seasonal changes in four reproductive traits, and developed a novel Bayesian risk‐partitioning model of daily nest survival to examine seasonal trends in two causes of nest failure. We found strong seasonal declines in reproductive traits for a subset of species. The probability of laying a full four‐egg clutch declined by 8–78% in 12 of 25 taxa tested, daily nest survival rates declined by 1–12% in eight of 22 taxa, incubation duration declined by 2.0–2.5% in two of seven taxa, and mean egg volume declined by 5% in one of 15 taxa. Temporal changes were not fully explained by individual variation. Across all species, the proportion of failed nests that were depredated declined over the season from 0.98 to 0.60, while the proportion abandoned increased from 0.01 to 0.35 and drove the seasonal declines in nest survival. An increase in abandonment of late nests is consistent with a life‐history tradeoff whereby either adult mortality increased or adults deserted the breeding attempt to maximize adult survival. In turn, seasonal declines in clutch size and incubation duration might be adaptive to hasten hatching of later nests. In other species of shorebirds, we found no seasonal patterns in breeding performance, suggesting that some species are not subject to selective pressure for early breeding.     

Flower power: Floral and resource manipulations reveal how and why reproductive trade‐offs occur for lowbush blueberry (Vaccinium angustifolium)

Plant reproductive trade‐offs are thought to be caused by resource limitations or other constraints, but more empirical support for these hypotheses would be welcome. Additionally, quantitative characterization of these trade‐offs, as well as consideration of whether they are linear, could yield additional insights. We expanded our flower removal research on lowbush blueberry (Vaccinium angustifolium) to explore the nature of and causes of its reproductive trade‐offs. We used fertilization, defoliation, positionally biased flower removal, and multiple flower removal levels to discern why reproductive trade‐offs occur in this taxon and to plot these trade‐offs along two continuous axes. We found evidence through defoliation that vegetative mass per stem may trade off with reproductive effort in lowbush blueberry because the two traits compete for limited carbon. Also, several traits including ripe fruit production per reproductive node and fruit titratable acidity may be “sink‐limited”—they decline with increasing reproductive effort because average reproductive structure quality declines. We found no evidence that reproductive trade‐offs were caused by nitrogen limitation. Use of reproductive nodes remaining per stem as a measure of reproductive effort indicated steeper trade‐offs than use of the proportion of nodes remaining. For five of six traits, we found evidence that the trade‐off could be concave down or up instead of strictly linear. Synthesis. To date, studies have aimed primarily at identifying plant reproductive trade‐offs. However, understanding how and why these trade‐offs occur represent the exciting and necessary next steps for this line of inquiry.     

Relationship between pre‐ and post‐copulatory traits in Salvator rufescens (Squamata: Teiidae)

Understanding pre‐ and post‐copulatory mechanisms of sexual selection can provide insights into the evolution of male reproductive strategies. The phenotype‐linked fertility hypothesis postulates that male sperm quality and secondary sexual characteristics will positively co‐vary, whereas the sperm competition hypothesis predicts a negative association between those traits. Male reproductive traits often show variation throughout the reproductive period, suggesting that the relationship between pre‐ and post‐copulatory sexual selection may vary temporally. Here, we evaluated the relationship between secondary sexual character and sperm traits and its temporal variation in Salvator rufescens, a south American lizard. We observed a negative relationship between jaw muscle and principal piece length of sperm and a variation in the relationship between pre‐ and post‐copulatory traits throughout the reproductive period. Collectively, our results evidenced a trade‐off between pre‐ and post‐copulatory traits and a strong seasonal flexibility of male reproductive strategies in this lizard species.     

Genetic dissection of pre‐anthesis sub‐phase durations during the reproductive spike development of wheat

Flowering time is an important factor affecting grain yield in wheat. In this study, we divided reproductive spike development into eight sub‐phases. These sub‐phases have the potential to be delicately manipulated to increase grain yield. We measured 36 traits with regard to sub‐phase durations, determined three grain yield‐related traits in eight field environments and mapped 15 696 single nucleotide polymorphism (SNP, based on 90k Infinium chip and 35k Affymetrix chip) markers in 210 wheat genotypes. Phenotypic and genetic associations between grain yield traits and sub‐phase durations showed significant consistency (Mantel test; r = 0.5377, P < 0.001). The shared quantitative trait loci (QTLs) revealed by the genome‐wide association study suggested a close association between grain yield and sub‐phase duration, which may be attributed to effects on spikelet initiation/spikelet number (double ridge to terminal spikelet stage, DR‐TS) and assimilate accumulation (green anther to anthesis stage, GA‐AN). Moreover, we observed that the photoperiod‐sensitivity allele at the Ppd‐D1 locus on chromosome 2D markedly extended all sub‐phase durations, which may contribute to its positive effects on grain yield traits. The dwarfing allele at the Rht‐D1 (chromosome 4D) locus altered the sub‐phase duration and displayed positive effects on grain yield traits. Data for 30 selected genotypes (from among the original 210 genotypes) in the field displayed a close association with that from the greenhouse. Most importantly, this study demonstrated specific connections to grain yield in narrower time windows (i.e. the eight sub‐phases), rather than the entire stem elongation phase as a whole.     

Reproductive and post‐reproductive life history of wild‐caught Drosophila melanogaster under laboratory conditions

The life history of the fruit fly (Drosophila melanogaster) is well understood, but fitness components are rarely measured by following single individuals over their lifetime, thereby limiting insights into lifetime reproductive success, reproductive senescence and post‐reproductive lifespan. Moreover, most studies have examined long‐established laboratory strains rather than freshly caught individuals and may thus be confounded by adaptation to laboratory culture, inbreeding or mutation accumulation. Here, we have followed the life histories of individual females from three recently caught, non‐laboratory‐adapted wild populations of D. melanogaster. Populations varied in a number of life‐history traits, including ovariole number, fecundity, hatchability and lifespan. To describe individual patterns of age‐specific fecundity, we developed a new model that allowed us to distinguish four phases during a female's life: a phase of reproductive maturation, followed by a period of linear and then exponential decline in fecundity and, finally, a post‐ovipository period. Individual females exhibited clear‐cut fecundity peaks, which contrasts with previous analyses, and post‐peak levels of fecundity declined independently of how long females lived. Notably, females had a pronounced post‐reproductive lifespan, which on average made up 40% of total lifespan. Post‐reproductive lifespan did not differ among populations and was not correlated with reproductive fitness components, supporting the hypothesis that this period is a highly variable, random ‘add‐on’ at the end of reproductive life rather than a correlate of selection on reproductive fitness. Most life‐history traits were positively correlated, a pattern that might be due to genotype by environment interactions when wild flies are brought into a novel laboratory environment but that is unlikely explained by inbreeding or positive mutational covariance caused by mutation accumulation.     

Genome‐wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America

Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome‐wide next‐generation sequencing of DNA pools (‘pool‐seq’) from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFβ/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome‐wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well‐known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.