Molecular footprints of the Holocene retreat of dwarf birch in Britain

Past reproductive interactions among incompletely isolated species may leave behind a trail of introgressed alleles, shedding light on historical range movements. Betula pubescens is a widespread native tetraploid tree species in Britain, occupying habitats intermediate to those of its native diploid relatives, B. pendula and B. nana. Genotyping 1134 trees from the three species at 12 microsatellite loci, we found evidence of introgression from both diploid species into B. pubescens, despite the ploidy difference. Surprisingly, introgression from B. nana, a dwarf species whose present range is highly restricted in northern, high‐altitude peat bogs, was greater than introgression from B. pendula, which is morphologically similar to B. pubescens and has a substantially overlapping range. A cline of introgression from B. nana was found extending into B. pubescens populations far to the south of the current B. nana range. We suggest that this genetic pattern is a footprint of a historical decline and/or northwards shift in the range of B. nana populations due to climate warming in the Holocene. This is consistent with pollen records that show a broader, more southerly distribution of B. nana in the past. Ecological niche modelling predicts that B. nana is adapted to a larger range than it currently occupies, suggesting additional factors such as grazing and hybridization may have exacerbated its decline. We found very little introgression between B. nana and B. pendula, despite both being diploid, perhaps because their distributions in the past have rarely overlapped. Future conservation of B. nana may partly depend on minimization of hybridization with B. pubescens, and avoidance of planting B. pendula near B. nana populations.     

Essential oil in Betula spp. leaves naturally growing in Estonia

The essential oils of the leaves of Betula pendula Roth., Betula pubescens Ehrh., Betula humilis Schrank and Betula nana L. (Betulaceae), and B. pendula buds were obtained using hydrodistillation and their chemical compositions were analyzed using gas chromatography and gas chromatography – mass spectrometry. The yield of essential oil was the highest (0.27%) in leaves of. B. humilis and the lowest in the leaves of B. pendula (0.11%), B. pubescens (0.05%) and B. nana (traces). Overall, 35–60 compounds were identified representing more than 92% of the total volume of the oils investigated. The bicyclical sesquiterpenoids found in the essential oils extracted from the leaves of B. pendula, B. pubescens and B. humilis and from the buds of B. pendula where α-betulenol (20.4–33.1%), α-betulenol acetate (5,6–28.9%), β-betulenal (5.5–6.2%) and β-betulenol (2.0–5.8%). The buds of B. pendula contained more birkenal, (E)-β-caryophyllene, δ-cadinene and caryophyllene oxide, but less α-betulenol acetate, β-betulenal and aliphatic compounds than the oils extracted from the leaves investigated. The oil composition of B. humilis and B. nana leaves was studied for the first time. The content of essential oil of B. nana leaves was much more varied: the oil contained less bicyclical sesquiterpenoids and more aliphatic compounds than the other Betula species investigated. The findings of our study could be used in chemosystematics of birch species, especially for B. nana.     

Sap‐feeding insects on forest trees along latitudinal gradients in northern Europe: a climate‐driven patterns

Knowledge of the latitudinal patterns in biotic interactions, and especially in herbivory, is crucial for understanding the mechanisms that govern ecosystem functioning and for predicting their responses to climate change. We used sap‐feeding insects as a model group to test the hypotheses that the strength of plant–herbivore interactions in boreal forests decreases with latitude and that this latitudinal pattern is driven primarily by midsummer temperatures. We used a replicated sampling design and quantitatively collected and identified all sap‐feeding insects from four species of forest trees along five latitudinal gradients (750–1300 km in length, ten sites in each gradient) in northern Europe (59 to 70°N and 10 to 60°E) during 2008–2011. Similar decreases in diversity of sap‐feeding insects with latitude were observed in all gradients during all study years. The sap‐feeder load (i.e. insect biomass per unit of foliar biomass) decreased with latitude in typical summers, but increased in an exceptionally hot summer and was independent of latitude during a warm summer. Analysis of combined data from all sites and years revealed dome‐shaped relationships between the loads of sap‐feeders and midsummer temperatures, peaking at 17 °C in Picea abies, at 19.5 °C in Pinus sylvestris and Betula pubescens and at 22 °C in B. pendula. From these relationships, we predict that the losses of forest trees to sap‐feeders will increase by 0–45% of the current level in southern boreal forests and by 65–210% in subarctic forests with a 1 °C increase in summer temperatures. The observed relationships between temperatures and the loads of sap‐feeders differ between the coniferous and deciduous tree species. We conclude that climate warming will not only increase plant losses to sap‐feeding insects, especially in subarctic forests, but can also alter plant‐plant interactions, thereby affecting both the productivity and the structure of future forest ecosystems.     

The influence of hybridization on epidermal properties of birch species and the consequences for palaeoclimatic interpretations

The Fennoscandian birch population primarily consists of Betula nana, B. pendula and B. pubescens ssp. czerepanovii, the Mountain birch. Frequent hybridization between the Mountain birch and B. nana generates a wide range of genotypic and phenotypic plasticity in the subarctic birch zone of Fennoscandia. Phases of subarctic conditions prevailed during the Late Glacial in large parts of NW Europe, and palynological as well as macrofossil analysis provide some evidence for the occurrence of birch hybrids during these intervals. Leaves from genetically controlled specimens of Betula pendula, B. pubescens ssp. czerepanovii, B. nana and the hybrids B. pubescens ssp. czerepanovii × nana and B. nana × pubescens ssp. czerepanovii are investigated for their specific characteristics of the epidermis morphology. Frequency and size of epidermal cells and stomata reveal a close affinity of both hybrids to B. nana and allow a differentiation of the intermediate forms between B. nana and the Mountain birch. With respect to palaeoatmospheric CO₂ reconstructions based on stomatal index, epidermal analysis shows that a possible occurrence of hybrids in fossil leaf assemblages has no profound consequences for combined species records. However, the significant differences observed in B. nana demand the separation of this species. A comparison of the cuticle properties of B. pendula and B. pubescens from Finnish Lapland and leaf material from The Netherlands reveals a divergence of the stomatal index that may be due to differences in day light length.     

Actual and potential tree-line in the North Atlantic region,especially in Greenland and the Faroes

Betula pubescens Ehrh. s.l. is the only forest-forming tree-species which has spread naturally to Iceland and southwesternmost Greenland. No trees have immigrated naturally to the Faroe Isles. Afforestation experiments and gardening have resulted in successful introduction of several tree species. Recent experiments comprise planting of Rocky Mountains tree-line species in SW-Greenland northward to Søndre Strømfjord (67°N), and the introduction to the Faroe Isles (62°N) of trees and shrubs from cool temperate forests of the southern hemisphere, e.g. Nothofagus spp. and Embothrium coccineum Forst. (Proteaceae).     

Size of the leaf as a marker of birch productivity at a distance from the climatic optimum

Morphometric parameters of mature leaves (area, venation density, and shape coefficient) and trunks (size and biomass) of two edificatory birch species, Betula pendula Roth and B. pubescens Ehrh., were characterized. These species inhabit a 1600-km transect of the Urals, from forest-tundra to forest-steppe. In both species, the highest trunk biomass was observed in the subzone of south taiga. Trunk parameters and leaf sizes of B. pendula varied in a narrower range than in B. pubescens. As distinct from leaf shape, leaf size was correlated with the average multiyear climate characteristics but not with the weather during a given growth season. In B. pubescens, which range extends wide from its climatic optimum, leaf area was positively correlated with trunk biomass. In B. pendula, growing mainly within the area of its climatic optimum, such correlation was not observed. We concluded that such parameter as mature leaf size could be used for prediction of tree productivity growing outside of the zone of their climatic optimum.     

Leaf surface traits: overlooked determinants of birch resistance to herbivores and foliar micro-fungi?

Previous studies on the mechanisms of birch resistance to herbivores and foliar micro-fungi (both pathogenic and endophytic) have focused mainly on the role of internal leaf chemistry. In the present study, we examined genetic correlations between leaf surface traits (glandular trichome density and total concentrations of surface flavonoid aglycones) and occurrence of three species of foliar micro-fungi, one pathogenic rust (Melampsoridium betulinum) and two endophytic fungi (Fusicladium sp. and Melanconium sp.), and performance of autumnal moth larvae (Epirrita autumnata) in two birch species, Betula pubescens ssp. czerepanovii and B. pendula. The performance of autumnal moth larvae on B. pubescens ssp. czerepanovii was negatively correlated with density of glandular trichomes (RGR: r=–0.855; pupal mass: r=–0.709). In addition, rust infection was negatively correlated with trichome density in B. pendula (r=–0.675) and with epicuticular flavonoid aglycones in B. pubescens ssp. czerepanovii (r=–0.855). The frequency of the endophytic fungus Fusicladium sp., was related to epicuticular flavonoid aglycones (r=–0.782), while another endophytic fungus, Melanconium sp., showed no associations with any of the studied variables in B. pubescens ssp. czerepanovii. Our results indicate that leaf surface traits may be at least as important determinants of herbivore performance and micro-fungi abundance in birch as leaf internal chemistry.     

Comparative analyses of plastid and AFLP data suggest different colonization history and asymmetric hybridization between Betula pubescens and B. nana

Birches (Betula spp.) hybridize readily, confounding genetic signatures of refugial isolation and postglacial migration. We aimed to distinguish hybridization from range‐shift processes in the two widespread and cold‐adapted species B. nana and B. pubescens, previously shown to share a similarly east–west‐structured variation in plastid DNA (pDNA). We sampled the two species throughout their ranges and included reference samples of five other Betula species and putative hybrids. We analysed 901 individual plants using mainly nuclear high‐resolution markers (amplified fragment length polymorphisms; AFLPs); a subset of 64 plants was also sequenced for two pDNA regions. Whereas the pDNA variation as expected was largely shared between B. nana and B. pubescens, the two species were distinctly differentiated at AFLP loci. In B. nana, both the AFLP and pDNA results corroborated the former pDNA‐based hypothesis that it expanded from at least two major refugia in Eurasia, one south of and one east of the North European ice sheets. In contrast, B. pubescens showed a striking lack of geographic structuring of its AFLP variation. We identified a weak but significant increase in nuclear (AFLP) gene flow from B. nana into B. pubescens with increasing latitude, suggesting hybridization has been most frequent at the postglacial expansion front of B. pubescens and that hybrids mainly backcrossed to B. pubescens. Incongruence between pDNA and AFLP variation in B. pubescens can be explained by efficient expansion from a single large refugium combined with leading‐edge hybridization and plastid capture from B. nana during colonization of new territory already occupied by this more cold‐tolerant species.     

A Lichenological Excursion to the West of Scotland

ABSTRACT

Background: Hybridisation can be a threat for the survival of a rare species because, in the case of insufficient numbers of appropriate mates, a rare form is much more likely to cross with a widespread taxon.

Aims: In the present study, we tested hypotheses concerning the level of hybridisation between endangered Betula humilis and its widespread congeners: B. pendula and B. pubescens as a function of habitat conditions.

Methods: We genotyped 312 individuals of three species using AFLP markers. B. humilis specimens were sampled in populations with low and high groundwater levels. Morphological identification of B. pubescens and B. pendula was verified using the Atkinson discriminant function.

Results: Altogether, 15 individuals (4.8%) were indicated as putative hybrids. The B. humilis hybrids were found in dry habitats and they could be classified as F1 or F2 generation. Tree hybrids could represent backcrosses to either B. pendula or B. pubescens.

Conclusions: Genetic analyses contradicted the idea that hybridisation between B. humilis and its close relatives was extensive. On the other hand, the presence of introgressed individuals in the populations in areas with low groundwater levels implied that pollen swamping might be a threat for declining B. humilis stands.     

A study of a population of birches in Glen Gairn

Summary

A small population of birches was studied with a view to measuring the rangeof morphological variation in Betula pendula Roth, and B. pubescens Ehrh. The general morphological features of crown and bark of trees in the area indicated that both species were present. A more detailed study was made of leaf and fruit characters, accepted by other investigators as being useful in discriminating between the species. It was shown that such characters varied continuously throughout the population and no tree was wholly pubescens- like or pendula-like in all characters.

The implications of these results are discussed in relation to the situation in Glen Gairn and to the general taxonomic situation of the British birches.     

Long‐distance dispersal and genetic structure of natural populations: an assessment of the inverse isolation hypothesis in peat mosses

It is well accepted that the shape of the dispersal kernel, especially its tail, has a substantial effect on the genetic structure of species. Theory predicts that dispersal by fat‐tailed kernels reshuffles genetic material, and thus, preserves genetic diversity during colonization. Moreover, if efficient long‐distance dispersal is coupled with random colonization, an inverse isolation effect is predicted to develop in which increasing genetic diversity per colonizer is expected with increasing distance from a genetically variable source. By contrast, increasing isolation leads to decreasing genetic diversity when dispersal is via thin‐tailed kernels. Here, we use a well‐established model group for dispersal biology (peat mosses: genus Sphagnum) with a fat‐tailed dispersal kernel, and the natural laboratory of the Stockholm archipelago to study the validity of the inverse isolation hypothesis in spore‐dispersed plants in island colonization. Population genetic structure of three species (Sphagnum fallax, Sphagnum fimbriatum and Sphagnum palustre) with contrasting life histories and ploidy levels were investigated on a set of islands using microsatellites. Our data show (, amova , IBD) that dispersal of the two most abundant species can be well approximated by a random colonization model. We find that genetic diversity per colonizer on islands increases with distance from the mainland for S. fallax and S. fimbriatum. By contrast, S. palustre deviates from this pattern, owing to its restricted distribution in the region, affecting its source pool strength. Therefore, the inverse isolation effect appears to hold in natural populations of peat mosses and, likely, in other organisms with small diaspores.     

Species turnover (β‐diversity) in ectomycorrhizal fungi linked to uptake capacity

Ectomycorrhizal (EcM) fungal communities may be shaped by both deterministic and stochastic processes, potentially influencing ecosystem development and function. We evaluated community assembly processes for EcM fungi of Pseudotsuga menziesii among 12 sites up to 400 km apart in southwest British Columbia (Canada) by investigating species turnover (β‐diversity) in relation to soil nitrogen (N) availability and physical distance. We then examined functional traits for an N‐related niche by quantifying net fluxes of , and protons on excised root tips from three contrasting sites using a microelectrode ion flux measurement system. EcM fungal communities were well aligned with soil N availability and pH, with no effect of site proximity (distance–decay curve) on species assemblages. Species turnover was significant (β_1/2 = 1.48) along soil N gradients, with many more Tomentella species on high N than low N soils, in contrast to Cortinarius species. Ammonium uptake was greatest in the spring on the medium and rich sites and averaged over 190 nmol/m²/s for Tomentella species. The lowest uptake rates of were by nonmycorrhizal roots of axenically grown seedlings (10 nmol/m²/s), followed by Cortinarius species (60 nmol/m²/s). EcM roots from all sites displayed only marginal uptake of nitrate (8.3 nmol/m²/s). These results suggest uptake capacity is an important functional trait influencing the assembly of EcM fungal communities. The diversity of EcM fungal species across the region arguably provides critical belowground adaptations to organic and inorganic N supply that are integral to temperate rainforest ecology.     

Potassium nutrition and free polyamines of Betula pendula Roth and Betula pubescens Ehrh.

Potassium and free polyamine concentrations in the leaves of silver birch (Betula pendula Roth) and downy birch (Betula pubescens Ehrh) were followed during three successive growing seasons 1996, 1997 and 1998 in order to define K deficiency levels.The highest foliar K concentrations were found in June. In August, the K concentrations were lower and remained quite stable in Betula pendula but varied in Betula pubescens. In addition to a common diamine, putrescine, and the polyamines, spermidine and spermine, also a less common diamine and polyamine, 1,3-diaminopropane and norspermidine were found in the birch leaves. The accumulation of both diamines, putrescine and 1,3-diaminopropane, was used to define the critical levels of K nutrition in birch leaves. Foliar K concentrations below 7–8 mg g ^–1 DW were found to correlate with the accumulation of diamines at most sampling dates.     

Balancing selection and heterozygote advantage in major histocompatibility complex loci of the bottlenecked Finnish wolf population

Maintaining effective immune response is an essential factor in the survival of small populations. One of the most important immune gene regions is the highly polymorphic major histocompatibility complex (MHC). We investigated how a population bottleneck and recovery have influenced the diversity and selection in three MHC class II loci, DLA‐DRB1, DLA‐DQA1 and DLA‐DQB1, in the Finnish wolf population. We studied the larger Russian Karelian wolf population for comparison and used 17 microsatellite markers as reference loci. The Finnish and Karelian wolf populations did not differ substantially in their MHC diversities ( = 0.047, P = 0.377), but differed in neutral microsatellite diversities ( = 0.148, P = 0.008). MHC allele frequency distributions in the Finnish population were more even than expected under neutrality, implying balancing selection. In addition, an excess of nonsynonymous compared to synonymous polymorphisms indicated historical balancing selection. We also studied association between helminth (Trichinella spp. and Echinococcus canadensis) prevalence and MHC diversity at allele and SNP level. MHC‐heterozygous wolves were less often infected by Trichinella spp. and carriers of specific MHC alleles, SNP haplotypes and SNP alleles had less helminth infections. The associated SNP haplotypes and alleles were shared by different MHC alleles, which emphasizes the necessity of single‐nucleotide‐level association studies also in MHC. Here, we show that strong balancing selection has had similar effect on MHC diversities in the Finnish and Russian Karelian wolf populations despite significant genetic differentiation at neutral markers and small population size in the Finnish population.     

Role of basidiospores as inocula of mycorrhizal fungi of birch

Summary Fruitbodies of sheathing mycorrhizal fungi collected under birch (Betula pendula andB. pubescens) were suspended over pots of soil and the resulting spore-supplemented soils were planted with sedlings ofB. pendula. Inocybe lacera, I. lanuginella, Hebeloma sacchariolens andH. leucosarx formed mycorrhizas readily.Lactarius pubescens andLeccinum roseofracta did not form mycorrhizas from basidiospore inocula, even after prolonged periods of seedling growth.Paxillus involutus gave equivocal results, perhaps because the soil was unsuitable for this species. Storage of the basidiospore-supplemented soils for 6 months in outdoor conditions or in a growth room at 18°C did not materially alter the results.The results are discussed in relation to the concept of mycorrhizal succession.     

Species-specific effect of UV-B radiation on the temporal pattern of leaf growth

Recent molecular and physiological studies have demonstrated that ultraviolet‐B radiation (UV‐B) can affect some of the processes involved in leaf growth, but the phases of leaf growth affected have not been clearly delimited. We used functional growth analysis to assess the effects of UV‐B radiation on the time course of leaf growth in seedlings of two birch species (Betula pendula and Betula pubescens). Our aim was to identify the phase(s) of leaf development affected by UV‐B radiation. In a greenhouse study, 1‐year‐old birch seedlings were subjected to three daily doses of supplemental UV‐B radiation treatments (UV‐B⁺) and no UV‐B radiation controls (UV‐B^?). Leaf growth measurements every 2 days were complemented by assessment of other functional traits over a 4‐week period at the start of the growing season. Using fitted curves, we were able to determine that the rate of leaf expansion was slowed by the UV‐B⁺ treatment in leaves of B. pendula because of a slower maximum leaf growth rate compared with plants under the UV‐B^? controls, but that compensation toward the end of the period of expansion negated this difference when leaves reached their final size. UV‐B⁺ had little effect on the rate of B. pubescens leaf growth despite a larger reduction in leaf final size due to UV‐B⁺ than occurred in B. pendula leaves. In conclusion, effective regulation ameliorated the effects of UV‐B radiation on leaf and seedling growth in B. pendula, whereas in B. pubescens, reductions in leaf final size under UV‐B⁺ were consistent with a slightly reduced rate of height growth.     

The genomic consequences of adaptive divergence and reproductive isolation between species of manakins

The processes of adaptation and speciation are expected to shape genomic variation within and between diverging species. Here we analyze genomic heterogeneity of genetic differentiation and introgression in a hybrid zone between two bird species (Manacus candei and M. vitellinus) using 59 100 SNPs, a whole genome assembly, and Bayesian models. Measures of genetic differentiation () and introgression (genomic cline center [α] and rate [β]) were highly heterogeneous among loci. We identified thousands of loci with elevated parameter estimates, some of which are likely to be associated with variation in fitness in Manacus populations. To analyze the genomic organization of differentiation and introgression, we mapped SNPs onto a draft assembly of the M. vitellinus genome. Estimates of , α, and β were autocorrelated at very short physical distances (< 100 bp), but much less so beyond this. In addition, average statistical associations (linkage disequilibrium) between SNPs were generally low and were not higher in admixed populations than in populations of the parental species. Although they did not occur with a constant probability across the genome, loci with elevated , α, and β were not strongly co-localized in the genome. Contrary to verbal models that predict clustering of loci involved in adaptation and isolation in discrete genomic regions, these results are consistent with the hypothesis that genetic regions involved in adaptive divergence and reproductive isolation are scattered throughout the genome. We also found that many loci were characterized by both exceptional genetic differentiation and introgression, consistent with the hypothesis that loci involved in isolation are also often characterized by a history of divergent selection. However, the concordance between isolation and differentiation was only partial, indicating a complex architecture and history of loci involved in isolation.     

The contribution of trees to ecosystem methane emissions in a temperate forested wetland

Wetland‐adapted trees are known to transport soil‐produced methane (CH₄), an important greenhouse gas to the atmosphere, yet seasonal variations and controls on the magnitude of tree‐mediated CH₄ emissions remain unknown for mature forests. We examined the spatial and temporal variability in stem CH₄ emissions in situ and their controls in two wetland‐adapted tree species (Alnus glutinosa and Betula pubescens) located in a temperate forested wetland. Soil and herbaceous plant‐mediated CH₄ emissions from hollows and hummocks also were measured, thus enabling an estimate of contributions from each pathway to total ecosystem flux. Stem CH₄ emissions varied significantly between the two tree species, with Alnus glutinosa displaying minimal seasonal variations, while substantial seasonal variations were observed in Betula pubescens. Trees from each species emitted similar quantities of CH₄ from their stems regardless of whether they were situated in hollows or hummocks. Soil temperature and pore‐water CH₄ concentrations best explained annual variability in stem emissions, while wood‐specific density and pore‐water CH₄ concentrations best accounted for between‐species variations in stem CH₄ emission. Our study demonstrates that tree‐mediated CH₄ emissions contribute up to 27% of seasonal ecosystem CH₄ flux in temperate forested wetland, with the largest relative contributions occurring in spring and winter. Tree‐mediated CH₄ emissions currently are not included in trace gas budgets of forested wetland. Further work is required to quantify and integrate this transport pathway into CH₄ inventories and process‐based models.     

Genetic linkage maps of white birches (<Emphasis Type="Italic">Betula platyphylla</Emphasis> Suk. and <Emphasis Type="Italic">B. pendula</Emphasis> Roth) based on RAPD and AFLP markers

A pseudo-testcross mapping strategy was used in combination with the random amplified polymorphism DNA (RAPD) and amplified fragment length polymorphism (AFLP) genotyping methods to develop two moderately dense genetic linkage maps for Betula platyphylla Suk. (Asian white birch) and B. pendula Roth (European white birch). Eighty F₁ progenies were screened with 291 RAPD markers and 451 AFLP markers. We selected 230 RAPD and 362 AFLP markers with 1:1 segregation and used them for constructing the parent-specific linkage maps. The resultant map for B. platyphylla was composed of 226 markers in 24 linkage groups (LGs), and spanned 2864.5 cM with an average of 14.3 cM between adjacent markers. The linkage map for B. pendula was composed of 226 markers in 23 LGs, covering 2489.7 cM. The average map distance between adjacent markers was 13.1 cM. Clustering of AFLP markers was observed on several LGs. The availability of these white birch linkage maps will contribute to the molecular genetics and the implementation of marker-assisted selection in these important forest species.     

Flowers for bees,autumn colours for whom? An experimental test with autumn colouration in mountain birches (Betula pubescens ssp. czerepanovii) as a signal against the moth Epirrita autumnata

In 2001, Hamilton and Brown proposed a controversial hypothesis of handicap signalling to potential insect parasites as an adaptive explanation for autumn leaf colouration. In subsequent studies there has been little attention to the costs and benefits of early autumnal colour change. Yet, in an observational study by Hagen et al. (2003) it was demonstrated that birch trees [Betula pubescens ssp. czerepanovii Ehrhart (Betulaceae)] turning yellow early in autumn had less damage from insects chewing on leaves the subsequent summer. Here, two experiments are presented which test the mechanisms in this model. The first addresses the proposed defence of leaves of B. pubescens ssp. czerepanovii by letting caterpillars of Epirrita autumnata Borkhausen (Lepidoptera: Geometridae), the birches’ most common insect parasites, choose between leaves from trees that either turned yellow late or early the foregoing autumn. The second experiment addresses whether adult female E. autumnata choose between early or late senescent (i.e., yellow or green) ‘twigs’ when ovipositing in autumn. We could not find evidence of preferences in either larvae or females, suggesting that timing of colour change in B. pubescens ssp. czerepanovii is not a warning signal that elicits a response in E. autumnata.