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 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.