The peppered moth and industrial melanism: evolution of a natural selection case study.

From the outset multiple causes have been suggested for changes in melanic gene frequency in the peppered moth Biston betularia and other industrial melanic moths. These have included higher intrinsic fitness of melanic forms and selective predation for camouflage. The possible existence and origin of heterozygote advantage has been debated. From the 1950s, as a result of experimental evidence, selective predation became the favoured explanation and is undoubtedly the major factor driving the frequency change. However, modelling and monitoring of declining melanic frequencies since the 1970s indicate either that migration rates are much higher than existing direct estimates suggested or else, or in addition, non-visual selection has a role. Recent molecular work on genetics has revealed that the melanic (carbonaria) allele had a single origin in Britain, and that the locus is orthologous to a major wing patterning locus in Heliconius butterflies. New methods of analysis should supply further information on the melanic system and on migration that will complete our understanding of this important example of rapid evolution.

Linkage map of the peppered moth, Biston betularia (Lepidoptera, Geometridae): a model of industrial melanism.

We have constructed a linkage map for the peppered moth (Biston betularia), the classical ecological genetics model of industrial melanism, aimed both at localizing the network of loci controlling melanism and making inferences about chromosome dynamics. The linkage map, which is based primarily on amplified fragment length polymorphisms (AFLPs) and genes, consists of 31 linkage groups (LGs; consistent with the karyotype). Comparison with the evolutionarily distant Bombyx mori suggests that the gene content of chromosomes is highly conserved. Gene order is conserved on the autosomes, but noticeably less so on the Z chromosome, as confirmed by physical mapping using bacterial artificial chromosome fluorescence in situ hybridization (BAC-FISH). Synteny mapping identified three pairs of B. betularia LGs (11/29, 23/30 and 24/31) as being orthologous to three B. mori chromosomes (11, 23 and 24, respectively). A similar finding in an outgroup moth (Plutella xylostella) indicates that the B. mori karyotype (n=28) is a phylogenetically derived state resulting from three chromosome fusions. As with other Lepidoptera, the B. betularia W chromosome consists largely of repetitive sequence, but exceptionally we found a W homolog of a Z-linked gene (laminin A), possibly resulting from ectopic recombination between the sex chromosomes. The B. betularia linkage map, featuring the network of known melanization genes, serves as a resource for melanism research in Lepidoptera. Moreover, its close resemblance to the ancestral lepidopteran karyotype (n=31) makes it a useful reference point for reconstructing chromosome dynamic events and ancestral genome architectures. Our study highlights the unusual evolutionary stability of lepidopteran autosomes; in contrast, higher rates of intrachromosomal rearrangements support a special role of the Z chromosome in adaptive evolution and speciation.

Selective bird predation on the peppered moth: the last experiment of Michael Majerus.

Colour variation in the peppered moth Biston betularia was long accepted to be under strong natural selection. Melanics were believed to be fitter than pale morphs because of lower predation at daytime resting sites on dark, sooty bark. Melanics became common during the industrial revolution, but since 1970 there has been a rapid reversal, assumed to have been caused by predators selecting against melanics resting on today's less sooty bark. Recently, these classical explanations of melanism were attacked, and there has been general scepticism about birds as selective agents. Experiments and observations were accordingly carried out by Michael Majerus to address perceived weaknesses of earlier work. Unfortunately, he did not live to publish the results, which are analysed and presented here by the authors. Majerus released 4864 moths in his six-year experiment, the largest ever attempted for any similar study. There was strong differential bird predation against melanic peppered moths. Daily selection against melanics (s ≈ 0.1) was sufficient in magnitude and direction to explain the recent rapid decline of melanism in post-industrial Britain. These data provide the most direct evidence yet to implicate camouflage and bird predation as the overriding explanation for the rise and fall of melanism in moths.

Evolutionary dynamics of multilocus microsatellite arrangements in the genome of the butterfly Bicyclus anynana, with implications for other Lepidoptera.

The sequences flanking microsatellites isolated from the butterfly Bicyclus anynana display high levels of similarity among different loci. We examined sequence data for evidence of the two mechanisms most likely to generate these similarities, namely recombination mediated events, such as unequal crossing over or gene conversion and through transposition of mobile elements (MEs). Many sequences contained tandemly arranged microsatellites, lending support to recombination as the multiplication mechanism. There is, however, also support for ME-mediated multiplication of microsatellites and their flanking sequences. Homology with a known Lepidopteran ME was found in B. anynana microsatellite regions, and polymorphic microsatellite markers with partial similarities in their flanking sequences were passed on to the next generation independently, indicating that they are not linked. Therefore, the rise of these similarities appears to be mediated through both processes, either as an interaction between the two, or by each being responsible for part of the observations. A large proportion of microsatellites embedded in repetitive DNA is representative for most studied butterflies and moths, and a BLAST survey of the B. anynana sequences revealed four short microsatellite-associated sequences that were present in many species of Lepidoptera. The similarities usually start to deviate beyond these sequences, which suggests that they define the extremes of a repeated unit. Further study of these conserved sequences may help to understand the mechanism underlying the multiplication events, and answer the question of why these redundancies are predominantly found in this insect group.

Molecular and ecological evidence for small-scale isolation by distance in an endangered damselfly, Coenagrion mercuriale.

Coenagrion mercuriale (Charpentier) (Odonata: Zygoptera) is one of Europe's most threatened damselflies and is listed in the European Habitats directive. We combined an intensive mark-release-recapture (MRR) study with a microsatellite-based genetic analysis for C. mercuriale from the Itchen Valley, UK, as part of an effort to understand the dispersal characteristics of this protected species. MRR data indicate that adult damselflies are highly sedentary, with only a low frequency of interpatch movement that is predominantly to neighbouring sites. This restricted dispersal leads to significant genetic differentiation throughout most of the Itchen Valley, except between areas of continuous habitat, and isolation by distance (IBD), even though the core populations are separated by less than 10 km. An urban area separating some sites had a strong effect on the spatial genetic structure. Average pairwise relatedness between individual damselflies is positive at short distances, reflecting fine-scale genetic clustering and IBD both within- and between-habitat patches. Damselflies from a fragmented habitat have higher average kinship than those from a large continuous population, probably because of poorer dispersal and localized breeding in the former. Although indirect estimates of gene flow must be interpreted with caution, it is encouraging that our results indicate that the spatial pattern of genetic variation matches closely with that expected from direct observations of movement. These data are further discussed with respect to possible barriers to dispersal within the study site and the ecology and conservation of C. mercuriale. To our knowledge, this is the first report of fine-scale genetic structuring in any zygopteran species.

Phylogeography of Aglais urticae (Lepidoptera) based on DNA sequences of the mitochondrial COI gene and control region.

Mitochondrial DNA (mtDNA) sequences of the COI gene and the control region were used to examine the genetic population structure of Aglais urticae L. (Lepidoptera) over its entire geographic range, i.e., the Palaearctic. The phylogenetic relationships within and between A. urticae subspecies were determined and patterns of mtDNA divergence and ecological differentiation were compared. High gene flow together with a recent and sudden population expansion characterise the genetic population structure of this species. No geographically induced differentiation was observed, nor were subspecies identified as separate evolutionary units. The discrepancy between the genetic and ecological variation is most likely due to the slower rate of mtDNA evolution compared to ecological differentiation. The control region proved to be a less useful molecular marker for the population genetics and the phylogenetic reconstruction of closely related taxa in A. urticae than it has for other species. The extreme bias in adenine and thymine content (A+T=90.91%) probably renders this region highly susceptible to homoplasy, resulting in a less informative molecular marker.

Effects of bottlenecks on quantitative genetic variation in the butterfly Bicyclus anynana.

The effects of a single population bottleneck of differing severity on heritability and additive genetic variance was investigated experimentally using a butterfly. An outbred laboratory stock was used to found replicate lines with one pair, three pairs and 10 pairs of adults, as well as control lines with approximately 75 effective pairs. Heritability and additive genetic variance of eight wing pattern characters and wing size were estimated using parent-offspring covariances in the base population and in all daughter lines. Individual morphological characters and principal components of the nine characters showed a consistent pattern of treatment effects in which average heritability and additive genetic variance was lower in one pair and three pair lines than in 10 pair and control lines. Observed losses in heritability and additive genetic variance were significantly greater than predicted by the neutral additive model when calculated with coefficients of inbreeding estimated from demographic parameters alone. However, use of molecular markers revealed substantially more inbreeding, generated by increased variance in family size and background selection. Conservative interpretation of a statistical analysis incorporating this previously undetected inbreeding led to the conclusion that the response to inbreeding of the morphological traits studied showed no significant departure from the neutral additive model. This result is consistent with the evidence for minimal directional dominance for these traits. In contrast, egg hatching rate in the same experimental lines showed strong inbreeding depression, increased phenotypic variance and rapid response to selection, highly indicative of an increase in additive genetic variance due to dominance variance conversion.

Effective population size, reproductive success and sperm precedence in the butterfly, Bicyclus anynana, in captivity.

A pedigree approach is used to estimate the effective population size yn two population cages of the butterfly, Bicyclus anynana. Each cage was founded with 54 individually marked adults of each sex. Matings were recorded over a 3-day period. Eggs were then collected from each female over a similar period before the numbers of hatching larvae were counted to assess progeny number. The males showed a higher variance in reproductive success than the females. Since about one-quarter of all females mated more than once, we also examined the pattern of sperm precedence using molecular markers or, in separate crossing experiments, wing pattern mutants. Both instances of complete first and last male sperm precedence, as well as of sperm mixing, were found. In some crosses a 'leakiness' was found in which some of the early eggs laid by a female were fertilized by a male partner which was subsequently completely unsuccessful. However, the estimates of effective population size were largely unaffected by the pattern of sperm precedence. Estimates for Ne : N in each cage were close to 0.60. The possibility of obtaining comparable estimates in selected natural populations of butterflies is discussed.

Inbreeding of bottlenecked butterfly populations. Estimation using the likelihood of changes in marker allele frequencies.

Polymorphic enzyme and minisatellite loci were used to estimate the degree of inbreeding in experimentally bottlenecked populations of the butterfly, Bicyclus anynana (Satyridae), three generations after founding events of 2, 6, 20, or 300 individuals, each bottleneck size being replicated at least four times. Heterozygosity fell more than expected, though not significantly so, but this traditional measure of the degree of inbreeding did not make full use of the information from genetic markers. It proved more informative to estimate directly the probability distribution of a measure of inbreeding, sigma2, the variance in the number of descendants left per gene. In all bottlenecked lines, sigma2 was significantly larger than in control lines (300 founders). We demonstrate that this excess inbreeding was brought about both by an increase in the variance of reproductive success of individuals, but also by another process. We argue that in bottlenecked lines linkage disequilibrium generated by the small number of haplotypes passing through the bottleneck resulted in hitchhiking of particular marker alleles with those haplotypes favored by selection. In control lines, linkage disequilibrium was minimal. Our result, indicating more inbreeding than expected from demographic parameters, contrasts with the findings of previous (Drosophila) experiments in which the decline in observed heterozygosity was slower than expected and attributed to associative overdominance. The different outcomes may both be explained as a consequence of linkage disequilibrium under different regimes of inbreeding. The likelihood-based method to estimate inbreeding should be of wide applicability. It was, for example, able to resolve small differences in sigma2 among replicate lines within bottleneck-size treatments, which could be related to the observed variation in reproductive viability.