The target is genetic variation in a reaction norm: how pigmentation produced by a genotype changes with developmental temperature. A gene affecting pigmentation at every temperature is not necessarily a gene affecting plasticity. The design must estimate inherited differences in responsiveness separately from differences in baseline colour.
Raise replicated genotypes or families across several controlled developmental temperatures, distributing members of the same genetic background among environments. Standardize sex, age, nutrition, density and the body region scored; female posterior abdominal pigmentation is a useful defined trait. Randomize position and scoring order and quantify colour from calibrated images. Record viability and development time, since selective mortality can change the composition of survivors and imitate a plastic response. A parental common-garden generation and reciprocal crosses can help separate inherited nuclear effects from maternal or carry-over effects.
A simple mixed-effects model is
where genotype and environment effects describe baseline differences and the shared temperature response, their interaction describes genotype-specific responses, and represents an appropriate block effect. Treating temperature as a continuous predictor gives genotype-specific slopes, but a curved response should not be forced into one line. Fit reaction norms with uncertainty and test whether their shapes differ.
Figure 1.
Illustrative pigmentation reaction norms separating baseline genetic differences from genotype-by-temperature interaction
. Details.
Quantitative trait locus mapping in segregating crosses can locate loci affecting the temperature interaction. A genome-wide association study in a diverse panel can complement it, with controls for relatedness and population structure. Analyse the joint phenotype across temperatures, rather than comparing a significant association at one temperature with a nonsignificant association at another and calling the difference an interaction. Mapping resolution, power and multiple-testing correction determine how many candidate regions can be distinguished.
Follow candidates with staged expression measurements, chromatin assays and functional perturbations. A coding change might alter enzyme temperature sensitivity; a cis-regulatory change might alter where or when the enzyme is produced; a trans-acting chromatin or hormonal regulator could affect several pathways. Allele-specific expression in a shared trans environment can help distinguish cis from trans effects. Replacing an allele or rescuing a perturbation and then measuring the entire reaction norm tests causality more strongly than expression correlation alone.
The expectation is a context-dependent genetic network, potentially including pigmentation enzymes, developmental regulators and chromatin regulation. Jean-Michel Gibert and colleagues found temperature-sensitive regulation involving Abdominal-B, bric-à-brac and chromatin regulators in female abdominal pigmentation. That supports a distributed mechanism rather than a universal dedicated plasticity gene. It does not imply that the same causal variants explain every natural population or all body regions.[ref]
Thus one may find variants affecting baseline colour, slope or curvature, as well as environment-dependent dominance and epistasis. Some apparent interactions may disappear after controlling development stage or survival. Demonstrating heritable plasticity does not establish that it is adaptive: that further claim requires fitness comparisons in the relevant environments.
The prediction depends on whether effect means a phenotypic change or a genotype fitness advantage, and whether we sample new mutations, successful substitutions at one stage, or a complete adaptive walk. These distributions differ even in simple models. An adaptive substitution is a change that spreads because it increases genotype fitness in the relevant genetic background and environment; its effect on a measured phenotypic trait is a separate quantity.
Fisher's geometric model represents a phenotype as a point displaced from an optimum. As derived in question 3, a mutation is favorable when . If its length is and the current distance is , favorable directions require . Thus cannot improve an isotropic distance-based genotype fitness, whereas small steps have some favorable directions. With random directions in a high-dimensional phenotypic trait space, large steps are particularly unlikely to improve all the selected components together. This is a probability-of-benefit argument, not yet a prediction for what fixes: extremely small advantageous effects are easily lost through genetic drift.
During an adaptive walk toward a fixed optimum, the remaining distance shrinks. Substantial changes can be favorable early, whereas fine adjustments dominate late. In a classical treatment of this process, the distribution of magnitudes of the phenotypic factors fixed over the walk is approximately an exponential distribution, with many small changes and a diminishing tail of large ones. This can reconcile a few prominent changes with numerous minor contributors; it is not the assertion that every step must be tiny.[ref] The shape concerns the stated model and effect measure. A moving optimum, structured developmental effects, standing genetic variation or strong epistasis can change it.
A related but distinct argument uses extreme value theory. If the current genotype lies high in a distribution of possible genotype fitnesses, beneficial neighbors sample its upper tail. For a broad class of rapidly declining tails, excesses above a high threshold approximately follow an exponential distribution. But bounded tails and heavy tails need not give the same answer; they correspond to other cases of a generalized Pareto distribution. Neither an exponential tail nor the assumption of a high initial threshold is a law of nature. Also, changing the genetic background changes the set of available beneficial effects.
Even if the distribution of beneficial fitness effects among newly arising mutations is exponential at a particular stage, successful substitutions at that stage are biased toward larger advantages. In an ideal haploid population with independent reproduction and offspring-number variance near one, let a rare new allele have relative genotype fitness , with , and assume that it completes its fate before the next successful beneficial lineage competes with it. The usual branching approximation gives fixation probability . More generally, if includes the relative arrival rates of beneficial mutations,
For the illustrative small-effect approximation , this becomes
It is a gamma distribution of shape 2, not an exponential distribution. The mathematical tail extends beyond the small- approximation, so this is biologically appropriate when is small enough that that tail carries negligible mass. With , the fraction exceeding is among arrivals but among substitutions in this approximation. This calculation and the pooled phenotypic-walk result answer different questions and are not contradictory.
Tests can begin with experimental evolution. Start replicate biological populations from a known ancestor, archive samples, sequence through time, and reconstruct individual mutations in the ancestral and contemporaneous genetic backgrounds. Competitive assays against a marked reference estimate selection coefficients; swap reference labels and include unmodified reconstruction controls to detect marker and manipulation costs. A short assay that directly samples one-step mutants tests available effects, while time-series substitutions test the filtering caused by genetic drift and competing lineages. In particular, experiments on two viruses rejected an exponential model for the beneficial effects they sampled, showing why tail shape must be tested rather than assumed.[ref]
Natural adaptation can be studied with quantitative trait locus mapping, association mapping and functional reconstruction. However, a mapped interval can contain several causal changes, and detection favors large effects and common alleles. Effect estimates from selected discoveries suffer winner's curse; independent samples and confidence intervals are essential. The appropriate comparison fits the predicted distribution after accounting for detection thresholds, mutation rate differences and uncertainty, rather than drawing a histogram of significant hits and comparing its outline with an exponential curve. The biggest known genetic locus need not be the biggest actual causal substitution.
Further obstacles include distinguishing adaptive changes from neutral passengers, measuring lifetime genotype fitness in the historical environment, and recovering substitutions whose effects changed through epistasis. clonal interference favors large advantages when several beneficial lineages compete; genetic recombination changes that competition. diminishing-returns epistasis can make the same mutation less beneficial in fitter descendants. Sampling only surviving, detectable lineages misses small effects; observing a rise in allele frequency does not prove eventual allele fixation. Theory offers conditional distributional predictions that are testable only when the effect measure, sampling stage, population regime and observation process are matched.
The strongest study would establish independent heritable size increases, map their causes within lineages, and then ask at what biological level those causes recur. Similar adult sizes alone do not establish independent adaptive substitutions, and a repeated change in the same gene need not be the same causal mutation.
Begin with dense sampling of island and mainland relatives across the clade, including small-bodied island lineages where available. Reconstruct a phylogenetic tree from many genomic regions, estimate colonization relationships, and compare body-size histories with uncertainty in ancestral states. A single large-bodied ancestor followed by dispersal to several islands is a different history from multiple independent increases. Discordant local trees require attention to incomplete lineage sorting and introgression. Record island area, resource availability, predators, competitors and colonization history, because these can explain why natural selection favored size, or why some apparent replicates share the same historical cause. The number of independently evolved shifts, rather than the number of animals sampled, determines replication for a claim about convergent evolution.
Define the phenotype carefully. Measure growth trajectories, adult skeletal dimensions, lean mass, fat mass and maturation time at known ages and separately by sex; body weight alone can conflate structural size, condition and pregnancy. Measure individuals blind to ancestry. Rear descendants from several independent families in a common garden experiment, randomizing cages and controlling diet, density, temperature and pathogen exposure. Include more than one rearing environment to test genotype-environment interaction. Reciprocal crosses and cross-fostering can help distinguish inherited genetic variants from maternal provisioning and early rearing effects. Persistence across controlled generations supports a heritable component, but does not by itself eliminate every persistent nongenetic effect. Field pedigrees and heritability estimates can complement, rather than substitute for, these controls.
Generate comparable genome data. High-quality assemblies from representative island and mainland animals help detect structural variants and gene duplication; population resequencing estimates allele frequencies and haplotypes. Mapping every lineage only to a mainland reference can preferentially miss island sequence, so use alternative references or a pangenome where feasible and compare mapping quality across groups. Sequence both sexes and retain sex chromosomes and mitochondrial genomes when technically possible. Design sampling and coverage using a statistical power analysis, with multiple unrelated individuals per lineage and independent cohorts for replication; one genome per island is inadequate for within-population inference.
Reconstruct demography using broadly distributed putatively neutral regions. Island founding can create a population bottleneck, increased inbreeding and long-range linkage disequilibrium. These processes can imitate aspects of a selective sweep. Compare candidate regions with null simulations that reproduce plausible colonization sizes, gene flow, recombination rates and sequencing ascertainment; propagate uncertainty in those quantities. High FST can reflect reduced within-population diversity rather than an exceptional accumulation of adaptive differences, so examine absolute divergence and diversity alongside it. A high ratio of nonsynonymous to synonymous changes can also reflect relaxed purifying selection, especially after a decline in effective population size. No single outlier statistic identifies a body-size mechanism.
Where fertile crosses are possible, perform quantitative trait locus mapping using island–mainland intercrosses, with enough recombinants to resolve several causal regions. Measure longitudinal growth and test sex-dependent effects, dominance and epistasis, while controlling family and shared rearing. Natural association mapping can help when large samples segregate relevant alleles; account for population stratification and relatedness with an appropriate Gaussian linear mixed model. Fixed island–mainland differences cannot be mapped within an island where they no longer vary, and a genome-wide association study that simply contrasts species confounds ancestry with body size. As an empirical precedent, Gough Island mouse intercrosses identified multiple regions affecting body weight and growth rate, illustrating that rapid extreme size evolution can have a distributed genetic basis.[ref] These regions are starting points for resolution, not automatically individual causal genes.
Integrate mapping with measurements in relevant tissues and developmental windows. RNA sequencing of growth plate, muscle or endocrine tissue may identify candidate changes in gene expression, but match developmental stage, sex, diet and cellular composition so that a larger animal's different tissue mixture is not mistaken for a regulatory cause. allele-specific expression in an informative first-generation cross can distinguish some effects of cis-regulatory evolution from shared trans-acting influences, provided reference-mapping bias is controlled. Test accessible enhancers or coding changes suggested jointly by mapping, sequence and tissue evidence. Growth pathways are reasonable hypotheses, but discovering a familiar growth-related annotation is not causal validation and should not exclude unexpected genetic loci.
Resolve the level of convergence explicitly. Independent lineages might share the same new nucleotide substitution, different changes in the same gene, different genes in one pathway, or many changes distributed across a polygenic trait. Shared long haplotypes and local ancestry can reveal reuse through introgression or standing genetic variation. Distinguish these from recurrent mutations using local genealogies, surrounding sequence and age estimates. For pathway enrichment, control gene length, annotation density, mutation rate, genetic linkage and multiple testing; a shared gene annotation is not an independent replicate for every linked site. Failure to find repeated individual hits does not exclude convergence at a broader biological level.
Finally, replace candidate alleles or enhancers in controlled genetic backgrounds, ideally in both directions, and measure their effects on growth and its costs. Use unedited siblings, sham-edited lines, independent edited founders and a rescue or precise reversion to address manipulation and off-target effects. Test combinations if mapping indicates epistasis. Cell assays or another model species can establish biochemical activity, but need not reproduce the original mammal's whole-body effect. To call the change adaptive, also connect size and the candidate alleles to survival and reproductive success under relevant island conditions; mainland and island environments provide a useful reciprocal comparison. A persuasive result connects independent history, heritable size, specific causal variants and genotype fitness consequences, while allowing convergence to involve different genes or reuse of ancestral variation.
Apparent resemblance generates competing hypotheses, not a diagnosis of mimicry. Establishing receiver-relevant similarity, evolutionary history and the biological fitness effect of resemblance are separate tasks. Shared ancestry, shared developmental responses and similar ecological demands may all contribute, and the origin of a pattern need not be the same as the process currently maintaining it.
First establish that the sampled individuals represent distinct species and comparable life stages and sexes. Voucher records, repeated observations, morphology and independent genetic evidence help distinguish species from geographical forms or sexually different members of one species. Sampling should span localities, seasons and within-species variation, since choosing only strikingly similar specimens exaggerates community-wide resemblance. For each individual, record dorsal and ventral patterns, size, resting posture, flight behaviour, microhabitat and activity time. An identical-looking patch in a photograph may rarely be exposed to the same receiver in nature.
Similarity must then be measured for plausible receivers. Predatory birds and butterfly mates need not share human colour vision. Calibrated reflectance measurements across visible and ultraviolet wavelengths, natural illumination and relevant backgrounds can support a colour vision model. A simplified relative receptor quantum catch is:
where is illumination, surface reflectance and the effective spectral sensitivity of receptor class . The comparison should include receptor noise and spatial resolution, because two spectrally distinguishable patches may be unresolved at attack distance. An apparent human match may conceal an ultraviolet or spatial difference. Conversely, experimentally demonstrated predator generalization can establish functional similarity despite a difference obvious to a human. A model predicts discrimination; behavioural tests verify whether the animal actually uses it.
One hypothesis is retention of an ancestral pattern through biological homology. Infer relationships from many independent genomic regions, then reconstruct pattern changes with uncertainty on the resulting phylogenetic tree. Resemblance concentrated among close relatives is compatible with retention. Similarity among distant relatives, especially matching different local patterns in different places, supports repeated change, although uncertainty in ancestral states must be retained. A single wing-pattern gene tree cannot serve as the independent species history when that very region might have crossed species boundaries.
A second hypothesis is Müllerian mimicry. If several species are defended and share predators, a common warning pattern can reduce the number of individuals attacked while predators learn avoidance. A local mimicry ring is predicted, with newcomers to an unfamiliar pattern suffering more attacks. Measure defence rather than infer it from colour: chemical assays, controlled feeding responses and subsequent predator behaviour address different components of palatability. predators from different localities may have different experience, so both naive and experienced receivers are informative. Several defended species can differ in defence strength; their interaction need not fit a perfectly symmetric mutual benefit.
A third hypothesis is Batesian mimicry, in which a relatively undefended species gains protection by resembling a defended model. This predicts protection dependent on the model's presence, abundance and the receiver's experience. A resemblance experiment should therefore vary both visual matching and the frequency or density of defended and undefended prey. A palatable mimic that becomes too common can weaken learned avoidance, whereas sharing a well-defended signal can reinforce it. These are ecological predictions to test, not a rule that each natural assemblage falls into one simple category. Defence varies with diet, age and the predator used in the assay.
The critical causal experiment changes the signal while controlling the rest of the prey. Artificial prey or wing models can independently vary colour and geometry against a standardized food reward. Randomized placement and position changes control background and location; handling and paint controls are needed with live butterflies. Record detection, attack and rejection separately: fewer missing prey can reflect poor detection rather than avoidance after recognition. Camera records help identify actual predators and distinguish attack from weather damage or removal by other animals. Artificial models lack movement and some chemical information, so their results should be checked against observations of live prey where feasible.
A field precedent is the release of different Heliconius cydno colour morphs into communities with different abundant co-models. Survival patterns supported an advantage to local matching, and the experiment also examined the prediction that this disadvantage should decrease with increased release density. [ref] This is more informative than simply mapping similar butterflies in the same forest. However, disappearance after release can include dispersal and unequal detectability; recapture effort, movement observations and attack evidence are needed when adapting such a design. Replication belongs at the site or independently manipulated community level, rather than treating every observation within one site as an independent community comparison.
Other ecological explanations deserve distinct predictions. Shared crypsis could make species match the same leaves, bark or light environment without resembling one another as defended models. Test detection against matched and mismatched backgrounds, separating the visibility of the prey from a predator's willingness to consume it. Shared thermoregulatory demands could favour similar absorption properties; temperature and performance measurements would then be relevant. Similar mate-detection environments could favour similar signals, but mate preference experiments must establish the relevant receiver and show how species avoid reproductive confusion. These mechanisms can coexist with defensive mimicry.
A common developmental or environmental cause is another possibility. Rear offspring from different populations and species in standardized and experimentally varied conditions, controlling larval diet, temperature and age. Persistence of differences under common rearing supports inherited contributions; condition-dependent convergence supports phenotypic plasticity. Crosses and quantitative trait locus mapping can then test the genetic architecture of patterning. Work on Optix in Heliconius combined mapping, expression and population comparisons to identify a shared regulator of red wing patterning. [ref] Use of the same developmental gene does not by itself prove a single recent origin of the pattern, nor prove that resemblance is adaptive.
Finally, adaptive introgression could transfer a pattern-associated region between species. Genome-wide ancestry should be compared with local ancestry around pattern regions, testing alternatives such as incomplete lineage sorting. The Heliconius genome study found evidence of exchange associated with mimetic adaptations among lineages. [ref] This concerns the route by which a useful pattern is acquired; selection by predators can still explain its maintenance. The strongest explanation combines receiver-specific similarity, a tested ecological benefit and an independently supported developmental and evolutionary history; neither human resemblance nor a shared pattern-associated gene supplies all three.