An exuberant, undescribed colour polymorphism in Theridion californicum (Araneae, Theridiidae): implications for a theridiid pattern ground plan and the convergent evolution of visible morphs

An opisthosomal (abdominal) colour polymorphism is described in the North American spider, Theridion californicum , comprising a plain Yellow morph and (at least) ten patterned morphs, which exhibit areas of red or black pigments superimposed on the yellow background, or no pigment (white). The polymorphism appears to be present throughout the species' range. The Yellow morph is the most frequent in populations, with patterned morphs all, individually, being rather rare. Progeny from known mothers were reared and indicate that the polymorphism is genetic and that Yellow is probably recessive to patterned morphs. Similar to other theridiids with well-studied colour polymorphisms, T. californicum occupies an under-leaf habitat and the variation in all these cases might be maintained by sight-hunting predators exerting negative frequency-dependent (apostatic) selection. In T. californicum , blocks of guanine underlying the pigmented hypodermis indicate a segmental patterning, which is not usually apparent in adult spiders. These segments, plus dorso-lateral divisions, permit the dorsal surface of the opisthosoma to be divided up into two mirror-image halves, each comprising 12 compartments. Each compartment can either lack pigment (thus appearing white as a result of underlying guanine) or be yellow, red, or black. All patterns in T. californicum can be derived from this ground plan, as can the morphs of other colour-polymorphic theridiids. It is suggested that selection for polymorphism, combined with constraints imposed by this theridiid ground plan, may have led to the convergent evolution of colour patterns across the family.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 23–34.     

Population structure and dispersal in a patchy landscape: nuclear and mitochondrial markers reveal area effects in the spider Theridion californicum (Araneae: Theridiidae)

Changes in land use have a major effect on patterns of biodiversity. However, few studies have examined the demographic and genetic shifts associated with a return to semi‐natural habitat following extended periods of human disturbance. Here we examine patterns of population structure in a spider restricted to the Pacific coastal strip of North America that exhibits an exuberant colour polymorphism. We use mitochondrial DNA and AFLP markers to examine genetic structure and estimate gene flow. The results show contrasting, gender‐specific patterns between these markers that suggest limited dispersal, combined with area effects most likely caused by expansion from refugial habitat patches following land‐management changes in a region of the San Francisco East Bay. Colour‐morph frequencies are not correlated with this complex genetic structure. Thus, unlike the classical area effects that were based on colour morphs, we demonstrate in T. californicum signals of historical contingency at neutral loci but not at the Colour locus, where traces of past events have been obliterated by balancing selection. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 104 , 600–620.     

SELECTION ON THE COLOR POLYMORPHISM IN HAWAIIAN HAPPY-FACE SPIDERS: EVIDENCE FROM GENETIC STRUCTURE AND TEMPORAL FLUCTUATIONS

Throughout this century genetic polymorphisms for color have been widely used as a research tool to allow insights into key evolutionary processes. Although color variants can often be diverse within populations, frequencies of different morphs may be similar across populations, either as a result of balancing selection or gene flow. Under these circumstances selection can be extremely difficult to demonstrate. Here we test for balancing selection on the naturally occurring color forms of the Hawaiian happy-face spider, Theridion grallator with two approaches. First, allozyme loci are used to generate a null model against which to test selection. Frequencies of alleles involved in the color polymorphism of T. grallator are used to generate another estimate for comparison. The results suggest that statistically similar frequencies of color morphs among populations of T. grallator may be maintained by some form of balancing selection. Second, we make use of an unusual event in which the normally stable frequencies of unpatterned and patterned morphs within a population were found to have shifted toward an excess of unpatterned morphs. We scored offspring of all fertilized, unpatterned (bottom-recessive) females found during this period of skewed morph frequencies and also in a year when morph frequencies were normal to deduce paternal color phenotypes. Mating was found to be random in the normal year, but in the perturbed year females had mated with rare (patterned) males twice as frequently as expected on the basis of the frequency of this morph type in the population. Both of these results are consistent with selection operating on the color polymorphism, and we speculate that apostatic selection, perhaps mediated by bird predators, may provide the mechanism.     

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.     

Sexual selection and genetic colour polymorphisms in animals

Genetic colour polymorphisms are widespread across animals and often subjected to complex selection regimes. Traditionally, colour morphs were used as simple visual markers to measure allele frequency changes in nature, selection, population divergence and speciation. With advances in sequencing technology and analysis methods, several model systems are emerging where the molecular targets of selection are being described. Here, we discuss recent studies on the genetics of sexually selected colour polymorphisms, aiming at (i) reviewing the evidence of sexual selection on colour polymorphisms, (ii) highlighting the genetic architecture, molecular and developmental basis underlying phenotypic colour diversification and (iii) discuss how the maintenance of such polymorphisms might be facilitated or constrained by these. Studies of the genetic architecture of colour polymorphism point towards the importance of tight clustering of colour loci with other trait loci, such as in the case of inversions and supergene structures. Other interesting findings include linkage between colour loci and mate preferences or sex determination, and the role of introgression and regulatory variation in fuelling polymorphisms. We highlight that more studies are needed that explicitly integrate fitness consequences of sexual selection on colour with the underlying molecular targets of colour to gain insights into the evolutionary consequences of sexual selection on polymorphism maintenance.     

Continental-Scale Footprint of Balancing and Positive Selection in a Small Rodent (Microtus arvalis)

Genetic adaptation to different environmental conditions is expected to lead to large differences between populations at selected loci, thus providing a signature of positive selection. Whereas balancing selection can maintain polymorphisms over long evolutionary periods and even geographic scale, thus leads to low levels of divergence between populations at selected loci. However, little is known about the relative importance of these two selective forces in shaping genomic diversity, partly due to difficulties in recognizing balancing selection in species showing low levels of differentiation. Here we address this problem by studying genomic diversity in the European common vole (Microtus arvalis) presenting high levels of differentiation between populations (average F _ST = 0.31). We studied 3,839 Amplified Fragment Length Polymorphism (AFLP) markers genotyped in 444 individuals from 21 populations distributed across the European continent and hence over different environmental conditions. Our statistical approach to detect markers under selection is based on a Bayesian method specifically developed for AFLP markers, which treats AFLPs as a nearly codominant marker system, and therefore has increased power to detect selection. The high number of screened populations allowed us to detect the signature of balancing selection across a large geographic area. We detected 33 markers potentially under balancing selection, hence strong evidence of stabilizing selection in 21 populations across Europe. However, our analyses identified four-times more markers (138) being under positive selection, and geographical patterns suggest that some of these markers are probably associated with alpine regions, which seem to have environmental conditions that favour adaptation. We conclude that despite favourable conditions in this study for the detection of balancing selection, this evolutionary force seems to play a relatively minor role in shaping the genomic diversity of the common vole, which is more influenced by positive selection and neutral processes like drift and demographic history.     

Spatial variation in colour morph, spotting and allozyme frequencies in the candy-stripe spider, Enoplognatha ovata (Theridiidae) on two Swedish archipelagos

The selective significance, if any, of many invertebrate visible polymorphisms is still not fully understood. Here we examine colour- and black spotting-morph frequencies in the spider Enoplognatha ovata in populations on two Swedish archipelagos with respect to different spatial scales and, in one archipelago, against the background of variation at four putative neutral allozyme marker loci. Every population studied was polymorphic for colour and 28 out of 30 contained all three colour morphs – lineata, redimita and ovata. We found no evidence for a breakdown in the traditional colour morph designation previously suggested for other northern European populations of this species. For colour there is no significant heterogeneity at spatial scales greater than between local populations within islands. Black spotting frequencies show a similar lack of pattern over larger spatial scales except that there are significant differences between the Stockholm and Göteborg archipelagos. Measures of population differentiation (θ) within the Stockholm islands for the two visible systems show them to be significantly more differentiated than the neutral markers, suggesting local selection acting on them in a population-specific manner. On the basis of previous observations and the distribution of spotting phenotypes on a European scale, it is argued that thermal selection might operate on black spotting during the juvenile stages favouring more spots in continental climates. It is not clear what selective forces act on colour.     

Cast adrift on an island: introduced populations experience an altered balance between selection and drift

A long-standing question in evolutionary biology is what becomes of adaptive traits when a species expands its range into novel environments. Here, we report the results of a study on an adaptive colour pattern polymorphism (stripes) of the coqui frog, Eleutherodactylus coqui, following its introduction to Hawaii from Puerto Rico. We compared population differentiation ( ) for the stripes locus-which underlies this colour pattern polymorphism-with neutral microsatellite loci to test for a signature of selection among native and introduced populations. Among native populations, for stripes were lower than expected under the neutral model, suggesting uniform balancing selection. Alternatively, among introduced populations, for stripes did not differ from the neutral model. These results suggest that the evolutionary dynamics of this previously adaptive trait have become dominated by random genetic drift following the range expansion.     

Heterogeneous genomic differentiation between walking-stick ecotypes: "isolation by adaptation" and multiple roles for divergent selection

Genetic differentiation can be highly variable across the genome. For example, loci under divergent selection and those tightly linked to them may exhibit elevated differentiation compared to neutral regions. These represent "outlier loci" whose differentiation exceeds neutral expectations. Adaptive divergence can also increase genome-wide differentiation by promoting general barriers to neutral gene flow, thereby facilitating genomic divergence via genetic drift. This latter process can yield a positive correlation between adaptive phenotypic divergence and neutral genetic differentiation (described here as "isolation-by-adaptation"). Here, we examine both these processes by combining an AFLP genome scan of two host plant ecotypes of Timema cristinae walking-sticks with existing data on adaptive phenotypic divergence and ecological speciation in these insects. We found that about 8% of loci are outliers in multiple population comparisons. Replicated comparisons between population-pairs using the same versus different host species revealed that 1-2% of loci are subject to host-related selection specifically. Locus-specific analyses revealed that up to 10% of putatively neutral (nonoutlier) AFLP loci exhibit significant isolation-by-adaptation. Our results suggest that selection may affect differentiation directly, via linkage, or by facilitating genetic drift. They thus illustrate the varied and sometimes nonintuitive contributions of selection to heterogeneous genomic differentiation.     

Evidence for adaptation from standing genetic variation on an antimicrobial peptide gene in the mussel Mytilus edulis

Genome scans of population differentiation identify candidate loci for adaptation but provide little information on how selection has influenced the genetic structure of these loci. Following a genome scan, we investigated the nature of the selection responsible for the outlying differentiation observed between populations of the marine mussel Mytilus edulis at a leucine/arginine polymorphism (L31R) in the antimicrobial peptide MGD2. We analysed DNA sequence polymorphisms, allele frequencies and population differentiation of polymorphisms closely linked to L31R, and pairwise and third‐order linkage disequilibria. An outlying level of population differentiation was observed at L31R only, while no departure from panmixia was observed at linked loci surrounding L31R, as in most of the genome. Selection therefore seems to affect L31R directly. Three hypotheses can explain the lack of differentiation in the chromosomal region close to L31R: (i) hitchhiking has occurred but migration and recombination subsequently erased the signal, (ii) selection was weak enough and recombination strong enough to limit the hitchhiking effect to a very small chromosomal region or (iii) selection acted on a pre‐existing polymorphism (i.e. standing variation) at linkage equilibrium with its background. Linkage equilibrium was observed between L31R and linked polymorphisms in every population analysed, as expected under the three hypotheses. However, linkage disequilibrium was observed in some populations between pairs of loci located upstream and downstream to L31R, generating a complex pattern of third‐order linkage disequilibria which is best explained by the hypothesis of selection on a pre‐existing polymorphism. We hypothesise that selection could be either balanced, maintaining alleles at different frequencies depending on the pathogen community encountered locally by mussels, or intermittent, resulting in sporadic fluctuations in allele frequency.     

Toll‐like receptor variation in the bottlenecked population of the Seychelles warbler: computer simulations see the ‘ghost of selection past’ and quantify the ‘drift debt’

Balancing selection can maintain immunogenetic variation within host populations, but detecting its signal in a postbottlenecked population is challenging due to the potentially overriding effects of drift. Toll‐like receptor genes (TLRs) play a fundamental role in vertebrate immune defence and are predicted to be under balancing selection. We previously characterized variation at TLR loci in the Seychelles warbler (Acrocephalus sechellensis), an endemic passerine that has undergone a historical bottleneck. Five of seven TLR loci were polymorphic, which is in sharp contrast to the low genomewide variation observed. However, standard population genetic statistical methods failed to detect a contemporary signature of selection at any TLR locus. We examined whether the observed TLR polymorphism could be explained by neutral evolution, simulating the population's demography in the software DIYABC. This showed that the posterior distributions of mutation rates had to be unrealistically high to explain the observed genetic variation. We then conducted simulations with an agent‐based model using typical values for the mutation rate, which indicated that weak balancing selection has acted on the three TLR genes. The model was able to detect evidence of past selection elevating TLR polymorphism in the prebottleneck populations, but was unable to discern any effects of balancing selection in the contemporary population. Our results show drift is the overriding evolutionary force that has shaped TLR variation in the contemporary Seychelles warbler population, and the observed TLR polymorphisms might be merely the ‘ghost of selection past’. Forecast models predict immunogenetic variation in this species will continue to be eroded in the absence of contemporary balancing selection. Such ‘drift debt’ occurs when a gene pool has not yet reached its new equilibrium level of polymorphism, and this loss could be an important threat to many recently bottlenecked populations.     

Identifying adaptive genetic divergence among populations from genome scans

The identification of signatures of natural selection in genomic surveys has become an area of intense research, stimulated by the increasing ease with which genetic markers can be typed. Loci identified as subject to selection may be functionally important, and hence (weak) candidates for involvement in disease causation. They can also be useful in determining the adaptive differentiation of populations, and exploring hypotheses about speciation. Adaptive differentiation has traditionally been identified from differences in allele frequencies among different populations, summarised by an estimate of FST. Low outliers relative to an appropriate neutral population-genetics model indicate loci subject to balancing selection, whereas high outliers suggest adaptive (directional) selection. However, the problem of identifying statistically significant departures from neutrality is complicated by confounding effects on the distribution of FST estimates, and current methods have not yet been tested in large-scale simulation experiments. Here, we simulate data from a structured population at many unlinked, diallelic loci that are predominantly neutral but with some loci subject to adaptive or balancing selection. We develop a hierarchical-Bayesian method, implemented via Markov chain Monte Carlo (MCMC), and assess its performance in distinguishing the loci simulated under selection from the neutral loci. We also compare this performance with that of a frequentist method, based on moment-based estimates of FST. We find that both methods can identify loci subject to adaptive selection when the selection coefficient is at least five times the migration rate. Neither method could reliably distinguish loci under balancing selection in our simulations, even when the selection coefficient is twenty times the migration rate.     

Functional and geographical differentiation of candidate balanced polymorphisms in Arabidopsis thaliana

Molecular population genetic analysis of three chromosomal regions in Arabidopsis thaliana suggested that balancing selection might operate to maintain variation at three novel candidate adaptive trait genes, including SOLUBLE STARCH SYNTHASE I (SSI) , PLASTID TRANSCRIPTIONALLY ACTIVE 7(PTAC7) , and BELL-LIKE HOMEODOMAIN 10 (BLH10). If balanced polymorphisms are indeed maintained at these loci, then we would expect to observe functional variation underlying the previously detected signatures of selection. We observe multiple replacement polymorphisms within and in the 32 amino acids just upstream of the protein–protein interacting BELL domain at the BLH10 locus. While no clear protein sequence differences are found between allele types in SSI and PTAC7, these two genes show evidence for allele-specific variation in expression levels. Geographical patterns of allelic differentiation seem consistent with population stratification in this species and a significant longitudinal cline was observed at all three candidate loci. These data support a hypothesis of balancing selection at all three candidate loci and provide a basis for more detailed functional work by identifying possible functional differences that might be selectively maintained.     

Comparative genome scan detects host‐related divergent selection in the grasshopper Hesperotettix viridis

In this study, we used a comparative genome scan to examine patterns of population differentiation with respect to host plant use in Hesperotettix viridis, a Nearctic oligophagous grasshopper locally specialized on various Asteraceae including Solidago, Gutierrezia, and Ericameria. We identified amplified fragment length polymorphism (AFLP) loci with significantly elevated F_ST (outlier loci) in multiple different‐host and same‐host comparisons of populations while controlling for geographic distance. By comparing the number and identities of outlier loci in different‐host vs. same‐host comparisons, we found evidence of host plant‐related divergent selection for some population comparisons (Solidago‐ vs. Gutierrezia‐feeders), while other comparisons (Ericameria‐ vs. Gutierrezia‐feeders) failed to demonstrate a strong role for host association in population differentiation. In comparisons of Solidago‐ vs. Gutierrezia‐feeding populations, a relatively high number of outlier loci observed repeatedly in different‐host comparisons (35% of all outliers and 2.7% of all 625 AFLP loci) indicated a significant role for host‐related selection in contributing to overall genomic differentiation in this grasshopper. Mitochondrial DNA sequence data revealed a star‐shaped phylogeny with no host‐ or geography‐related structure, low nucleotide diversity, and high haplotype diversity, suggesting a recent population expansion. mtDNA data do not suggest a long period of isolation in separate glacial refugia but are instead more compatible with a single glacial refugium and more recent divergence in host use. Our study adds to research documenting heterogeneity in differentiation across the genome as a consequence of divergent natural selection, a phenomenon that may occur as part of the process of ecological speciation.     

Isolation by adaptation in Neochlamisus leaf beetles: host-related selection promotes neutral genomic divergence

This study tests how divergent natural selection promotes genomic differentiation during ecological speciation. Specifically, we use adaptive ecological divergence (here, population divergence in host plant use and preference) as a proxy for selection strength and evaluate the correlation between levels of adaptive and genetic differentiation across pairwise population comparisons. Positive correlations would reveal the pattern predicted by our hypothesis, that of 'isolation by adaptation' (IBA). Notably, IBA is predicted not only for selected loci but also for neutral loci. This may reflect the effects of divergent selection on neutral loci that are 'loosely linked' to divergently selected loci or on geneflow restriction that facilitates genetic drift at all loci, including neutral loci that are completely unlinked to those evolving under divergent selection. Here, we evaluate IBA in maple- and willow-associated populations of Neochlamisus bebbianae leaf beetles. To do so, we collected host preference data to construct adaptive divergence indices and used AFLPs (amplified fragment length polymorphisms) and mitochondrial sequences to quantify genetic differentiation. Partial Mantel tests showed significant IBA in 'pooled' analyses of putatively neutral and of putatively selected ('outlier') AFLP loci. This pattern was also recovered in 12% of 'locus-specific' analyses that separately evaluated genetic differentiation at individual neutral loci. These results provided evidence for widespread effects of selection on neutral genomic divergence. Our collective findings indicate that host-related selection may play important roles in the population genomic differentiation of both neutral and selected gene regions in herbivorous insects.     

Differential gene exchange between parapatric morphs of Littorina saxatilis detected using AFLP markers

Speciation requires the acquisition of reproductive isolation, and the circumstances under which this could evolve are of great interest. Are new species formed after the acquisition of generalized incompatibility arising between physically separated populations, or may they arise as a result of the action of disruptive selection beginning with the divergence of a rather restricted set of gene loci? Here we apply the technique of amplified fragment length polymorphism (AFLP) analysis to an intertidal snail whose populations display a cline in shell shape across vertical gradients on rocky shores. We compare the F_ST values for 306 AFLP loci with the distribution of F_ST estimated from a simulation model using values of mutation and migration derived from the data. We find that about 5% of these loci show greater differentiation than expected, providing evidence of the effects of selection across the cline, either direct or indirect through linkage. This is consistent with expectations from nonallopatric speciation models that propose an initial divergence of a small part of the genome driven by strong disruptive selection while divergence at other loci is prevented by gene flow. However, the pattern could also be the result of differential introgression after secondary contact.     

A comparison of AFLPs and microsatellites to identify the population structure of brown trout (Salmo trutta L.) populations from Hardangervidda, Norway

A total of 495 fish from 11 Hardangervidda lakes were genotyped in order to compare amplified fragment length polymorphisms (AFLP) and microsatellites in terms of their capacity to infer population genetic structure. The 11 microsatellites used in this study gave a greater polymorphism information content and greater gene diversity, with an average of 14.8 alleles per locus, than the six AFLP primer combinations used. However, the AFLPs resulted in 178 polymorphic loci and a 3.1 times larger marker index (effective multiplex ratio multiplied with the gene diversity). Comparable population structuring, for example in terms of distinguishing fish from the different river systems, was obtained with both marker systems. An AFLP and microsatellite multilocus Bayesian assignment test with the structure program divided the fish into six groups largely concurrent with main branches on a population neighbour-joining tree. Yet, the admixture status of individuals is mostly contradictory in the AFLP and the microsatellite analyses. The results are discussed concerning migration between lake populations.     

Signatures of balancing selection are maintained at disease resistance loci following mating system evolution and a population bottleneck in the genus Capsella

ABSTRACT: BACKGROUND: Population bottlenecks can lead to a loss of variation at disease resistance loci, which could have important consequences for the ability of populations to adapt to pathogen pressure. Alternatively, current or past balancing selection could maintain high diversity, creating a strong heterogeneity in the retention of polymorphism across the genome of bottlenecked populations. We sequenced part of the LRR region of 9 NBS-LRR disease resistance genes in the outcrossing Capsella grandiflora and the recently derived, bottlenecked selfing species Capsella rubella, and compared levels and patterns of nucleotide diversity and divergence with genome-wide reference loci. RESULTS: In strong contrast with reference loci, average diversity at resistance loci was comparable between C. rubella and C. grandiflora, primarily due to two loci with highly elevated diversity indicative of past or present balancing selection. Average between-species differentiation was also reduced at the set of R-genes compared with reference loci, which is consistent with the maintenance of ancestral polymorphism. CONCLUSIONS: Historical or ongoing balancing selection on plant disease resistance genes is a likely contributor to the retention of ancestral polymorphism in some regions of the bottlenecked Capella rubella genome.     

AFLP utility for population assignment studies: analytical investigation and empirical comparison with microsatellites

Individual-based population assignment tests have thus far mainly relied on the use of microsatellite loci. However, the logistic difficulty of screening large numbers of loci required to reach sufficient statistical power hampers the usefulness of microsatellites in situations of weak population structuring. Amplified fragment length polymorphisms (AFLP) represents an alternative for overcoming this logistical issue as the technique allows the user to characterize a much larger number of loci with a comparable analytical effort. In this study, an assignment test based on maximum likelihood for dominant markers was used to investigate the potential usefulness of AFLP for population assignment. We also compared assignment success achieved with AFLP with that obtained using microsatellites in a case study of low population differentiation involving whitefish (Coregonus clupeaformis) sympatric ecotypes. The analytical investigation showed that the minimum number of AFLP loci required to reach an assignment success of 95% stood within values that are easily achievable in many situations. This also showed how assignment success varied according to the number of AFLP loci used, their absolute frequency and their frequency differential and sampling errors, as well as the number of putative source populations. The case study showed that given a comparable analytical effort in the laboratory, AFLP were much more efficient than the microsatellite loci in discriminating the source of an individual among putative populations. AFLP resulted in higher assignment success at all levels of stringency and the log-likelihood differences between populations obtained with AFLP for each individual were much larger than those obtained with microsatellites. These results indicate that research involving individual-based population assignment methods should benefit importantly from the use of AFLP markers, especially in systems characterized by weak population structuring.     

Long‐term balancing selection on chromosomal variants associated with crypsis in a stick insect

How polymorphisms are maintained within populations over long periods of time remains debated, because genetic drift and various forms of selection are expected to reduce variation. Here, we study the genetic architecture and maintenance of phenotypic morphs that confer crypsis in Timema cristinae stick insects, combining phenotypic information and genotyping‐by‐sequencing data from 1,360 samples across 21 populations. We find two highly divergent chromosomal variants that span megabases of sequence and are associated with colour polymorphism. We show that these variants exhibit strongly reduced effective recombination, are geographically widespread and probably diverged millions of generations ago. We detect heterokaryotype excess and signs of balancing selection acting on these variants through the species’ history. A third chromosomal variant in the same genomic region likely evolved more recently from one of the two colour variants and is associated with dorsal pattern polymorphism. Our results suggest that large‐scale genetic variation associated with crypsis has been maintained for long periods of time by potentially complex processes of balancing selection.