Host-plant-derived variation in ultraviolet wing patterns influences mate selection by male butterflies
Lehrstuhl Tierökologie I, Universität Bayreuth, D-95440 Bayreuth, Germany
*Author for correspondence (e-mail: Helge.Knuettel{at}uni-bayreuth.de)
Accepted April 19, 2001
![]() |
Summary |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
We discuss the use of principal component analysis in analysing spectral data in the context of visual communication. We propose the alternative application of confidence intervals of averaged spectra as a novel straightforward statistical method for comparing groups of spectra in a manner that is independent of assumptions about the visual system of the receiver. In addition, they can be used to give confidence intervals to derived measures of colour such as quantum catch by photoreceptors.
Key words: insect, butterfly, animalplant interaction, vision, signal, colour, pigment, ultraviolet light, flavonoid, behaviour, mate choice, Lycaenidae, Polyommatus icarus, principal component analysis, spectral data, confidence intervals.
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
How are wing colours and wing patterns produced in the Lepidoptera? The wing coloration is caused by the wings cover of coloured scales, an apomorphy that also gave the order its scientific name. The colours of the scales have several bases. Physical colours arise from the very elaborate scale structure. Interference colours result from multiple layers, three-dimensional diffraction gratings, Tyndall-blues and plain back scattering of light (Ghiradella, 1991). The pigments in the scales make an even greater contribution to wing colour. Most of these pigments are synthesised de novo during scale development in the pupa. The processes of pattern formation and pigment synthesis are highly regulated by switching on and off genes in ontogenesis (Nijhout, 1991; Brakefield et al., 1998).
A second source of wing pigments is secondary plant compounds taken up from the larval diet (Wilson, 1987; Nijhout, 1991). These pigments need to be stored by the larvae, retained through all instars and then transferred to the developing wings in the pupa. Little information has been published on these processes and on their relevance for the animals.
During the last few years, data on the sequestration of flavonoids by blue butterflies (Lycaenidae) from their larval food plants have become available (Geuder et al., 1997; Schittko et al., 1999; Burghardt, 2000). This work was done mainly on the common blue butterfly Polyommatus icarus, which is also the focal species of this study. It was shown that some of the flavonoids in the larval diet are taken up, partially metabolised by glycosylation (a process that probably serves as a means of detoxification), and stored. Sequestered flavonoids are retained up to the pupal stage, probably in the fat body. During the late pupal stage, most of the flavonoids are transferred to the wings, more or less simultaneously with the onset of pigment synthesis in the wing scales (Kornmaier, 1999). The sequestration of flavonoids appears to be a very specific process, as is also suggested by flavonoids being kept through all instars and transferred specifically to the wings.
The flavonoids are a major class of secondary plant compounds. All flavonoids are ultraviolet-absorbing, and depending on the type and placement of functional groups and side chains bonded to the central flavan, may possess human-visible colours (Harborne, 1991). The anthocyanins, a large subclass of the flavonoids, are the most important flower pigments (Brouillard and Dangles, 1993). P. icarus seems to sequester only flavonols, a subclass of flavonoids, which appear colourless to yellow to the human eye (Burghardt, 2000).
The ecological function of flavonoid sequestration in P. icarus is unknown. The specific allocation of these pigments to the wings called for an investigation of their effects on wing coloration. After it became clear that flavonoids do alter wing colours (see below), questions arose as to whether this serves as a colour signal to conspecifics and what information may be transferred by such a signal. We therefore set out to investigate the effects of flavonoids on butterfly wing coloration in detail.
The appraisal of colour in the context of animal communication must be carried out in a way that is independent of the human visual system (Endler, 1990; Cuthill and Bennett, 1993; Bennett et al., 1994). Colours have to be measured physically for the whole wavelength range visible to the species under study. This includes the ultraviolet range from 300nm to 400nm, which is usually not perceivable by humans without technical help (Endler, 1990). However, ultraviolet vision is well known in many vertebrates and invertebrates (Menzel, 1979; Jacobs et al., 1991; Jacobs, 1992; Tovée, 1995), and in other species ultraviolet light has been demonstrated to have at least a behavioural influence (Sharma et al., 1999). We therefore determined wing coloration by measuring wing spectral reflectance in a manner that was independent of human vision. Reflectance spectra covered the human-visible and the ultraviolet range (300nm to 690nm). This is the range in which vision in butterflies is made possible by a combination of the light available in nature (lower limit approximately 300nm) (Henderson and Hodgkiss, 1963; Henderson, 1970) and the spectral sensitivities of butterfly photoreceptors (up to approximately 700nm) (e.g. Eguchi et al., 1982; Bernard and Remington, 1991).
Subsequent analyses of spectral data must also be performed independently of the human visual system. For most species, little or no data on visual physiology are available, so a statistical approach that is independent of any assumptions on a receivers psychology must be chosen (Endler, 1990; Cuthill et al., 1999). Principal component analysis (PCA) with subsequent multivariate analysis of variance (MANOVA) was introduced by Cuthill et al. (Cuthill et al., 1999) as a very sensitive tool for comparing spectra. Here, we will evaluate this method and compare it with an alternative, newly proposed procedure using averaged spectra and their confidence limits.
![]() |
Materials and methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
P. icarus butterflies exhibit a strong sexual dimorphism, with males having an iridescent bright blue to violet upperside and females having dark brown uppersides with a submarginal row of more or less pronounced orange spots (Fig.1). The undersides of both sexes are much less dimorphic, with submarginal orange spots and a number of black spots with white surrounds set on a background that is greyish in males and more brownish in females (Fig.2). In this study, we looked only at the wing colours of female butterflies.
|
|
Caterpillars were reared on an artificial diet based on peas and wheatgerm, with some dried alfalfa (Medicago sativa) leaves added as a feeding stimulus. This diet contains only traces of flavonoids, and none of these flavonoids is sequestered by P. icarus larvae (Burghardt, 2000). To this mixture, the flavonoid quercetin (quercetin dihydrate, purum; Roth, Karlsruhe, Germany) (2.5% dry mass) was added to yield a second diet that differed only in flavonoid content. Larvae fed this quercetin-spiked diet sequester large amounts of this flavonol, most of which is stored as quercetin-3-O-galactoside (Burghardt, 2000).
Reflectance spectroradiometry
The spectral reflectance of wing colours was measured using an InstaspecII diode-array photometer equipped with an MS125 spectrograph (400 lines mm-1 grating) (Oriel, Stratford). The measurement range was 300nm to 690nm, with a resolution of approximately 0.5nm. Reflectance was measured relative to a Spectralon99% white reflectance standard (Labsphere, Congleton). An Oriel 77184 sighting optics with fibre optics, extension tubes and an Ultra-Fluar 10/0.20 quartz microscope objective (Zeiss, Jena) allowed simultaneous viewing and measurement of wing patterns. All measurements were taken from a direction normal to the wing plane. Measured spot diameters ranged from 0.11mm to 0.21mm, and the numerical apertures (N.A.) of the measuring beam lay between 0.14 and 0.15. Illuminating light was delivered from an XBO75W/2OFR xenon arc lamp (Osram, Munich) with an Oriel 68806 stabilised power supply, Oriel fibre optics and focusing beam probe. Illuminating light was focused onto the wings (N.A.=0.070.09) at 45° to the wing plane. Illumination was always from the distal direction. To ensure exactly the same wing orientation for all measurements, the wings were positioned with their hind margins parallel to the projection of the illuminating beam into the wing plane. The wings were mounted on black velvet to eliminate stray light. The wings were handled very carefully with pointed forceps at the most proximal end only, where no measurements were taken.
In all subsequent analyses, only the spectral reflectance data for the hindwings are presented. A clear gradient of ultraviolet reflectance was present in the forewings, while no such spatial variation was apparent in the hindwings within a given food treatment (see Results). For each individual, a randomly chosen hindwing was used in the measurements. For each of the human-visible colour categories (black, white, orange and background for undersides; brown and orange for uppersides), a number of measurements were taken from each wing (usually 5 for black, orange and brown and 10 for white and background). Individual measurements of a given colour category were sampled from all areas of a wing. For every measurement, the exact position (i.e. pattern element or position relative to veins) on the wing was recorded.
Comparing spectral data
Principal component analysis (PCA) (for an introduction to the method see) (Kim and Mueller, 1978a; Kim and Mueller, 1978b; Chatfield and Collins, 1980) has become increasingly popular for comparing spectra (Young, 1986; Bakker and Arnold, 1993; Cronin et al., 1997; Cuthill et al., 1999; Grill and Rush, 2000). PCA reduces the amount of data to a small number of variables for each spectrum. These new variables can then be used for statistical analyses, such as MANOVA (Kim and Mueller, 1978a; Kim and Mueller, 1978b; Chatfield and Collins, 1980; Cuthill et al., 1999). This method has the advantage of being independent of any assumptions about visual physiology. It simply makes a statistical comparison of spectra. We will explain this method briefly. (For more details on the usage of principal component analysis with spectra, see) (Endler, 1990; Cuthill et al., 1999).
We start with spectra, which are measured values for a number of wavelength intervals. These wavelength intervals can be seen as variables that are highly intercorrelated because most biological objects exhibit broadband absorption or reflection of light. The idea of PCA is to transform these i dependent variables into i new, orthogonal variables, the principal components (PCs) using linear combinations. The i data points of a spectrum are transformed into i PC scores during this operation. The PCs are chosen in sequence so that each will extract the maximum variation in the data matrix. For reflectance spectra, the first three PCs (PC1PC3) will usually explain approximately 95% of the total variation (Endler, 1990; Cuthill et al., 1999; see Results below). So the first few PCs can then be used in statistical analyses, in our case a comparison between different food treatments. The meaning of a principal component is expressed in its PC coefficients, which are the values used to transform the original data points into PC scores.
Before PCA, we reduced the wavelength resolution of our spectra, which originally consisted of 795 points each, by calculating mean values for 10nm intervals. Means represented our data better than medians, and there were no outlying datapoints to pose a problem. The original and reduced data were almost identical, as judged by visual inspection of the spectra. We decided to separate the analysis of mean reflectance of spectra from that of spectral shape, as proposed by Endler (Endler, 1990) and applied for the first time by Cuthill et al. (Cuthill et al., 1999). To achieve this, we subtracted from every spectrum its mean reflectance and then used these standardised spectra for PCA (StatSoft, 1999a).
Mean reflectance, as a measure of brightness, was analysed by analysis of variance (ANOVA). The first three PCs of standardised spectra were analysed separately, as a measure of spectral shape, by MANOVA. ANOVA and MANOVA (StatSoft, 1999a) were both repeated-measures designs with larval diet as a between-group factor and spectra as replicate measurements from individual butterflies. PCA and statistical analyses were performed separately for every colour category.
An alternative method of comparing groups of spectra is to use the means of the spectra and the confidence intervals of these means. Calculation of these statistics is simple and straightforward. For a given group of n spectra, each made up of i data points (measured values at wavelength i), one calculates the arithmetic mean and standard error for every value of
i. Multiplication of the tabulated value of Students t for n-1 degrees of freedom at a chosen probability level, e.g. P=95%, by the standard errors yields the confidence intervals of the mean spectrum (Sokal and Rohlf, 1981; Sachs, 1997). If the confidence intervals of two groups of spectra do not overlap, one can reject the null hypothesis that the parametric means of both groups are the same at a probability level P (if n>10 for both groups). One can compare for overlap of confidence intervals in specific wavelength ranges, and at the same time, relate to absolute reflectance values in these ranges.
To demonstrate the use of confidence intervals of spectra in derived measures of colour, we calculated the quantum catch of hypothetical photoreceptors for mean spectra and their confidence limits (Bernard and Remington, 1991). Mean spectra and confidence limits transform into mean values of quantum catch with confidence intervals for the wing colours. We used the nomogram templates given by Stavenga et al. (Stavenga et al., 1993). Values for the maximum sensitivity of the photoreceptors were taken from the work of Bernard and Remington (Bernard and Remington, 1991) on two Lycaena species because these are the only such data published for lycaenid butterflies. The illuminating light was the generalised CIE daylight spectrum D65 (Schanda, 1997) transformed to quantum flux and normalised to unity.
Courtship and mate choice in Polyommatus icarus
At favourable ambient temperatures, male P. icarus patrol the habitat in search for females. After locating a female, a male approaches her, usually exhibits a distinct fluttering behaviour next to the female, alights, sometimes still fluttering, and tries to copulate (see Pellmyr, 1982). Female polyommatine blues usually mate only once (Häuser, 1990; Drummond, 1984). Males may mate several times, but require at least a day for the production of a new spermatophore (for related species, see Drummond, 1984). Male approaches to females occur quite frequently but are only very rarely successful. The courtship sequence is almost never completed. Usually the male decides to stop at some point in the sequence before trying to copulate and goes on searching for another female. Alternatively, the female refuses copulation when the male tries to grasp her with his genitalia. Rape is impossible because in the Polyommatus section the female has to actively evert a sclerotised gonoporus for copulation (Häuser, 1993). Therefore, both sexes exhibit mate choice.
We observed male reactions to female dummies made from dead dried females. Freshly emerged females were allowed to stretch and harden their wings, then killed by freezing, placed on insect pins and dried at room temperature. Dummies were prepared from females that had been raised either on an artificial diet from which no flavonoids are sequestered or on the same artificial diet but with quercetin added. Dummies were prepared with their wings in postures found in the field: the wings closed above the body, the wings opened at an angle of approximately 90° or the wings spread at almost 180°. Prior to the behavioural experiments, the dummies were dried for 2 days at 4050°C to evaporate any volatile substances.
Mate-choice experiments in the laboratory
Experiments were performed in a flight cage (2.1m x 1.1m x 1.0m high) in an environmental chamber at 30°C. The cage was lit by fluorescent lamps (at >10kHz) with visible and ultraviolet light (eight Osram Biolux L58W/72-965 and two Osram L36W/73). Lamp tubes were positioned 20cm above the cage ceiling made of a transparent plastic screen. In this cage the butterflies exhibited essentially all the behaviour patterns observed in the field, and they courted and mated readily.
In the centre of the cage, nine wooden poles (diameter 10mm, length 23cm) were placed upright in a square with one pole in the centre and three poles along each side of the square (or three rows of three sticks each; both poles within a row and rows 10cm apart). The dummies were pinned onto these poles during the experiments. Three replicates of each of the two dummy types (differing in either food treatment, size or alteration of the external appearance; see below) were used in each test. The remaining three positions were left empty. The dummies were positioned so that the nearest neighbours of every dummy were always either (i) one of the other dummy types or (ii) empty poles. Positions were changed every 30min to avoid positional bias. A complete experiment included nine such rounds of 30min each, resulting in a total observation time of 270min per experiment. During an experiment, a single dummy was placed in all nine positions and all positions were occupied for the same time by the different dummy types. The size of the dummies was controlled, with the mean wing lengths of dummy types within an experiment differing by no more than 0.5mm (except where the effect of dummy size was tested).
The experiments were performed with our breeding population from the flight cage. However, most of the females in the cage were removed just before an experiment and immediately replaced afterwards because they tended to sit and bask for prolonged periods on the poles and thus to block them (see below). 5080 butterflies were in the cage during the experiments, but not all of these were active. The numbers of active butterflies were more or less constant during an experiment but differed between experiments.
During the experiments, the reactions of the males to the dummies were counted and recorded on a dictaphone. The observer sat outside the flight cage and watched the animals through a white mosquito net from a darker part of the room. The presence of the observer had no noticeable effect on the behaviour of the butterflies. The reactions of the males fell into several clearly distinguishable categories. If a male changed his flight path and approached a female dummy or steered directly towards it, this was scored as approach. If a male exhibited the distinct fluttering behaviour next to the dummy, usually circling around it, this was counted as fluttering. Landings and copulation attempts were the other categories distinguished. When a butterfly landed on a pole and remained there, this position was considered as occupied because the reactions of other males to this pole might have been caused by the presence of the living butterfly rather than the dummy. Reactions to occupied poles were therefore not scored, and the duration of occupancy recorded. Similarly, reactions to a pole were not scored when another individual was flying in the vicinity.
In a series of six experiments (three experimental designs, each replicated twice), we tested male reactions to dummies made from flavonoid-free versus flavonoid-rich females. In each experiment, only dummies of the same posture type (wings opened at 90°, wings fully spread, or wings closed) were presented. With the wings opened at 90°, both ventral and dorsal wing surfaces were visible. In experiments with the wings fully spread or closed, each dummy was placed on a round cardboard disk (diameter 3.5cm) mounted on top of the pole. With the wings fully spread, the undersides of the wings were mostly hidden by the white, ultraviolet-absorbing cardboard. In experiments with the wings closed, only the undersides of the wings were visible.
A seventh experiment tested the influence of dummy size. A group of small (mean wing length 11.5±0.2mm) versus large (mean wing length 15.4±0.3mm; means ± S.D., N=6) flavonoid-rich dummies was presented (wings spread at 90°).
In an eighth experiment, we used dummies with altered external appearance. Dummies were made as usual from flavonoid-free females with the wings opened at 90°. The undersides of the dummies were then sprayed with a solution of the quercetin-3-O-diglycoside rutin (purum; Roth, Karlsruhe, Germany; 2mgml-1 in ethanol). Rutin is frequently found in adult P. icarus (Burghardt, 2000). The absorption spectrum (300690nm) of rutin is virtually identical to that of other quercetin-3-O-glycosides, particularly to quercetin-3-O-galactoside, which is the dominant flavonoid in butterflies reared on the quercetin-spiked diet (Burghardt, 2000; H. Knüttel, unpublished results). The externally applied rutin reduces ultraviolet reflectance of the wings, and therefore mimics the light-absorbing effect of dietary flavonoids on wing coloration. At the same time, any other potential physiological effects of dietary flavonoids which might alter the constitution of the dummies were precluded. The rutin solution was cautiously delivered from a vaporiser and dried off almost immediately, leaving a layer of tiny rutin crystals on the wing scales. This was repeated until reflectance spectroradiometry indicated colours very similar to those of naturally flavonoid-rich butterflies. The second group of dummies in this experiment was sprayed with ethanol only.
The use of Pearsons chi-squared tests in the statistical analysis of frequency-type behavioural data, although being chosen quite frequently by many researchers, is highly debated (e.g. Hurlbert, 1984; Lombardi and Hurlbert, 1996; Wilson and Dugatkin, 1996). Therefore, we decided to use an alternative approach (Kramer and Schmidhammer, 1992). Each of the experiments described above was divided into subsets. For each subset the proportions of reactions to the two dummy types (corrected for actual observation time; see below) were calculated, thereby collapsing all reactions of the males in a subset into a single data point for each dummy type. Subsequently, using a two-sided t-test for dependent data sets, we tested the null-hypothesis that these proportions did not differ between the two dummy types in an experiment, i.e. that the reactions of the male population in the flight cage were not biased with respect to dummy type. In the experiments that were run twice, the proportions of both experiments were collated and subjected to a single t-test.
In the experiments which tested the effect of dummy size, or of external application of rutin, each subset for statistical analysis comprised one round of 30min (resulting in nine estimates of relative proportions of reactions). In the other experiments (and their replicates), the frequency of reactions within a round of 30min was sometimes rather low, so we pooled three consecutive rounds for a subset (resulting in six estimates per test situation).
The numbers of reactions and duration of occupancy by butterflies were summed for each subset and for each of the two dummy types tested in an experiment. Actual observation times for every dummy type were calculated by subtracting the duration of occupancy from the total observation time during a subset (3 dummiesx30min=90min, or 3 dummiesx90min=270min). To account for variable actual observation times, the numbers of reactions scored during a subset were standardised to a total observation time of 90min, or 270min, to yield corrected numbers of reactions. From these values, we calculated the proportions of reactions to the two dummy types in each subset.
Mate-choice experiments in the field
We tested the effects of the flavonoid content of female dummies on their attractiveness for wild free-flying males. Experiments were performed as described by Burghardt et al. (Burghardt et al., 2000) with dummies made from females reared on the artificial diets detailed above.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
These findings are corroborated by the reflectance spectra (Fig.3). The spectra of the dark brown uppersides of the wings increased monotonically with wavelength (Fig.3F). The orange spots showed very low and flat reflectance spectra up to about 500nm with a sharp increase in reflectance at higher wavelengths and a short drop at approximately 660nm (Fig.3D,E). The mean spectra for both food treatments and their confidence intervals were almost the same except for the orange spots from the undersides of the wings (Fig.3DF). In these, reflectance was lower in the ultraviolet range (300400nm) in flavonoid-rich animals (Fig.3D). This difference was significant, as indicated by non-overlapping confidence intervals, but very small in absolute value (maximum approximately 1%). The reflectance spectra of white spots (Fig.3A) distinctly differed between food treatments. The spectral reflectance of the white spots of flavonoid-free animals was more or less flat in the visible range (400690nm) and decreased slowly at shorter wavelengths. In contrast, the reflectance of the white spots of flavonoid-rich animals was very low in the ultraviolet region but increased drastically from approximately 400nm upwards to reach the values of flavonoid-free animals at approximately 450nm. At even longer wavelengths, it increased slightly, to exceed the values of the flavonoid-free animals above approximately 620nm. The spectra of the background coloration (Fig.3B) monotonically increased, those of flavonoid-free animals being flatter. The background spectra of flavonoid-rich animals remained very low in the ultraviolet region but exceed those of flavonoid-free animals above 420nm. A similar pattern was found for the spectra of the black spots (Fig.3C), but with much lower absolute reflectance values.
|
|
Eigenvalues for the first few PCs were very high, i.e. they explained most of the between-spectra variation in spectral shape. All PC1s extracted more than 67% of the variation, and the first three PCs together always accounted for more than 95% of the variation (Table2).
|
|
We compared confidence limits of the calculated quantum catch for hypothetical photoreceptors (Fig.5). As expected from the analyses presented above, values for the ultraviolet receptor (P360) were significantly different for all colours of the undersides of the wings, as were values for the blue receptor (P437) for white spots. In the background coloration of the undersides of the wings, the blue receptor exhibited no significant difference but the green (P500) and red (P568) receptors did. The values for the red receptor were significantly different for the black spots, too.
|
|
Mate-choice experiments in the field
In the field, we presented female dummies to wild, free-flying males. Of the males passing a dummy within a radius of 15cm, 51.1% showed a noticeable response, while 49.9% ignored the dummies (both flavonoid-rich and flavonoid-free). The fraction of by-passing males did not differ between flavonoid-rich and flavonoid-free dummies (Fig.6). Hence, males were not more attracted to either dummy category. There was, however, a very clear difference between the two dummy types in the mode of reaction of responsive males. At flavonoid-rich dummies, 65.3% of the male responses ended in the more intense fluttering behaviour, but this happened at only 41.6% of the flavonoid-free dummies. Therefore, males were more interested in flavonoid-rich dummies, spending more time and energy there (highly significant after BonferroniHochberg correction for two dependent tests on the same data set).
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Although shown only for females in this study, flavonoids influence wing colours in both sexes of this species and for a broad range of natural diets (H. Knüttel, unpublished results). Hence, there is a large intraspecific variability in wing coloration due to individual ontogenesis in this butterfly species. Flowers and buds are the preferred diet of P. icarus and more flavonoids are generally sequestered from these parts of the plant (Burghardt et al., 1997; Burghardt et al., 2001). Larvae grow faster and larger and beneficial myrmecophilous interactions with ants are more strongly developed when they are fed protein-rich flowers rather than foliage (Burghardt and Fiedler, 1996). The amount of flavonoids sequestered may, under natural conditions, be correlated with the quality of the food, which may in turn relate to higher quality as a potential mate. In butterflies, both the body size and the protein content derived from larval nutrition typically indicate high female fecundity (Honek, 1993; Wheeler, 1996). In P. icarus, flavonoid content and concentration is generally positively related to body mass, further strengthening the idea of a link between flavonoid sequestration during ontogenesis and fitness-related characters (Burghardt et al., 2001).
In laboratory experiments, we found that males indeed preferred flavonoid-rich female dummies. It seems very probable that this preference is mediated visually and not by chemoreception. Flavonoids, and their glycosides in particular, are not volatile (Harborne, 1991), and we attempted to evaporate any other volatile substances from our dummies prior to the behavioural experiments. We found differences in wing colour only in the undersides of the wings. However, a male preference for flavonoid-rich dummies was found irrespective of which sides of the wings were visible. In the dummies with the wings fully spread, the undersides of the wings were not completely hidden, and body coloration, which is also influenced by flavonoids (Fig.2), may have served as an additional cue to the males. In these dummies, the proportion of male reactions to the flavonoid-rich dummies was lower than in the dummies of the other posture types, in which the undersides of the wings were readily visible. Moreover, visual signals from the undersides of the wings seem to play the major role in mate recognition from a distance in many lycaenid butterflies (Douwes, 1976; Wago et al., 1976; Pellmyr, 1982). Therefore, it seems likely that visual signals from the wing undersides are mainly responsible for the males preference for flavonoid-rich female dummies.
Still, it is possible that flavonoids in the diet may trigger the production of, or serve as precursors of other unknown substances that are sensed and preferred by the males. So, in another experiment, we used artificially altered dummies that were originally flavonoid-free but had the undersides of their wings sprayed with the flavonoid rutin frequently found in P. icarus. In these dummies, any potential physiological change in chemical composition due to the sequestration of flavonoids during ontogenesis was precluded. Rutin-sprayed, ultraviolet-absorbing dummies were again highly preferred, eliminating the possibility of olfactory cues being the stimuli mediating male preference.
We then tested male preference in the field to avoid possible limits imposed by indoor behavioural experiments. In spite of the presence of living females in the natural habitat, the males frequently responded to the dummies. Neither flavonoid-rich nor flavonoid-free dummies attracted more males, but males spent more time and energy closely inspecting the flavonoid-rich dummies. Exactly the same results were obtained when using flavonoid-rich or flavonoid-poor dummies reared on natural host-plants (Burghardt et al., 2000).
From these results, we conclude that, once having detected a female, males of P. icarus are more strongly attracted to flavonoid-rich females. This preference is mediated visually and is probably due mostly to the colours of the undersides of the wings. As flavonoid-induced colour differences are most pronounced in the ultraviolet range, this seems to be the spectral range most important to this behaviour. The ultraviolet coloration of undersides of the wings is to a large degree environmentally determined, and intraspecific variation is enormous. It is not therefore a species-specific characteristic, as has been assumed previously and used taxonomically (see Knüttel and Fiedler, 2001; for Gonepteryx spp. see also Brunton and Majerus, 1995).
Our study is the first case known to us in which secondary plant compounds have been shown to alter the external appearance of an insect and where this signal is used in visually mediated mate choice. In contrast, this phenomenon is well known in vertebrates, in which carotenoids are of special importance in this context (e.g. Grether et al., 1999; Blount et al., 2000).
One of the goals of this study was to evaluate the statistical analysis of spectra by PCA and MANOVA as introduced by Cuthill et al. (Cuthill et al., 1999) and to compare this approach with the use of mean spectra and their confidence intervals, used here for the first time. It soon became apparent that PCA of spectra in combination with MANOVA is a very sensitive tool in grouping spectra, as Cuthill et al. (Cuthill et al., 1999) found in an analysis of iridescent bird coloration. This method has the advantage of being independent of any assumptions about the visual system of the receiver of a signal. A serious drawback of PCA is the lack of any straightforward explanation of the meaning of the principal components and that this meaning is a different one for every set of data analysed. It is not possible to compare PCs resulting from different studies (Endler, 1990; Cuthill et al., 1999). The addition of a few spectra to a set of data may completely change the result of a PCA (Endler, 1990).
We do, however, consider it valuable not only to be able to show that there is a difference between groups of spectra, but also to be able to locate the spectral origin of this difference precisely. The PC coefficients are the values used to transform the original data into the PCs. From these, it is possible to attempt to identify the meaning of a PC. However, the spectra of PC coefficients often are quite complex in shape, even for very smooth and straight (reflectance) spectra (Fig.3, Fig.4), so this task is not trivial. We did not accomplish very much with our data. The situation becomes even more complicated when one tries to take into account the factor loadings, which are the correlations between the original variables and the PCs. Here, the picture becomes even worse. One not only really does not know what a PC actually means, but there is no mathematical justification stating how many of the PCs, and therefore how much of the variation in the data, should be included in the subsequent statistical analysis. However, the MANOVA results may depend on this choice.
There are only rules of thumb as to how many PCs to include in later analyses. Either PCs are considered with an eigenvalue larger than 1 or that account for more than 5% of total variation (Kaiser, 1960; Chatfield and Collins, 1980; StatSoft, 1999b). The decision is arbitrary and depends on the experimenter. Biologically relevant information may, however, be found in a PC that accounts for a very small amount of the variation and is therefore not considered for statistical analysis (Cuthill et al., 1999). We decided a priori to always include the first three PCs in our analyses. If we had chosen to include only PCs accounting for more than 5% of the variation, we would have had to leave out the PC3 in all analyses except for the black spots. Thus, we would have missed a significant difference between food treatments in the orange spots of the undersides of the wings. Although very small in absolute value, this difference was found and located in the spectrum simply by looking at the confidence limits of the mean spectra of these spots.
The wavelength resolution of the spectra and, therefore, the number of variables used for PCA need to be considered in this context. The eigenvalue of a PC is the fraction of total variation it explains multiplied by the number of variables (StatSoft, 1999b). The effect of the decision to include PCs with an eigenvalue larger than 1 then depends heavily on the wavelength resolution chosen.
As a result of all these complications involved in PCA and MANOVA, we found the use of mean spectra and their confidence intervals very valuable. The calculation is easily accomplished, and the results can be plotted together with the original data for purposes of analysis and illustration. The concept of confidence intervals is easy to comprehend and quite powerful. At the same time, the comparison of confidence intervals is a very sharp tool because we found significant differences in the same groups of spectra as we found using PCA and MANOVA. In contrast to using PCA in combination with MANOVA, it was possible directly to locate the spectral ranges that differed between treatments. Hence, it was possible to relate statistical differences to absolute variation in reflectance values within a specific range. For example, one can question whether the minute, but nevertheless statistically significant, differences found in the orange and black spots, which are surrounded by much brighter colours, are of much biological relevance. Another important advantage of mean spectra and their confidence intervals is that they can be used in measures of colour derived from these spectra, as was illustrated for the quantum catch of hypothetical photoreceptors. Using this latter approach, we do not suggest that P. icarus possesses photoreceptors with the properties assumed in our calculations, nor do we imply specific neuronal interactions between the outputs of the photoreceptors. Rather, we want to demonstrate the potential for further use of confidence intervals in species in which more is known about their visual physiology. There is evidence that true colour vision is involved in the feeding of P. icarus (A. Kelber and H. Knüttel, unpublished results). Finally, unlike the results of PCA, mean spectral curves with their confidence intervals can readily be compared between different studies.
In summary, the present study demonstrates that host-plant-derived allochthonous pigments can radically change colour patterns in butterflies, which in turn has significant consequences for mate selection. The effects of flavonoids on colours mostly occur in the ultraviolet range and their precise measurement and quantitative analysis still pose some challenges. However, tools are accumulating to deal with these complexities, and with further progress in visual physiology, the underlying mechanism will be revealed. Even in the absence of such knowledge, the profound consequences of host-plant-derived coloration for ecology and evolutionary biology are becoming evident.
![]() |
Acknowledgments |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Bakker, J. and Arnold, G. M. (1993). Analysis of sensory and chemical data for color evaluation of a range of red port wines. Am. J. Enol. Viticult. 44, 2734.
Bennett, A. T. D., Cuthill, I. C. and Norris, K. J. (1994). Sexual selection and the mismeasure of color. Am. Nat. 144, 848860.
Bernard, G. D. and Remington, C. L. (1991). Color vision in Lycaena butterflies: Spectral tuning of receptor arrays in relation to behavioral ecology. Proc. Natl. Acad. Sci. USA 88, 27832787.[Abstract]
Blount, J. D., Houston, D. C. and Moller, A. P. (2000). Why egg yolk is yellow. Trends Ecol. Evol. 15, 4749.[Medline]
Brakefield, P. M., Kesbeke, F. and Koch, P. B. (1998). The regulation of phenotypic plasticity of eyespots in the butterfly Bicyclus anynana. Am. Nat. 152, 854860.
Brouillard, R. and Dangles, O. (1993). Flavonoids and flower colour. In The Flavonoids: Advances in Research Since 1986 (ed. J. B. Harborne), pp. 565588. London: Chapman & Hall.
Brunton, C. F. A. and Majerus, M. E. N. (1995). Ultraviolet colours in butterflies: intra- or inter-specific communication? Proc. R. Soc. Lond. B 260, 199204.
Burghardt, F. (2000). Stoffwechsel und ökologische Funktion von Flavonoiden im Bläuling Polyommatus icarus (Lepidoptera: Lycaenidae). PhD Dissertation. Bayerische Ludwig-Maximilians-Universität, Würzburg.
Burghardt, F. and Fiedler, K. (1996). The influence of diet on growth and secretion behaviour of myrmecophilous Polyommatus icarus caterpillars (Lepidoptera: Lycaenidae). Ecol. Ent. 21, 18.
Burghardt, F., Fiedler, K. and Proksch, P. (1997). Flavonoidaufnahme im Bäuling Polyommatus icarus: Der Einfluß von Nahrungspflanze, Größe und Geschlecht. Mitt. Dt. Ges. Allg. Angew. Ent. 11, 533536.
Burghardt, F., Knüttel, H., Becker, M. and Fiedler, K. (2000). Flavonoid wing pigments increase attractiveness of female common blue (Polyommatus icarus) butterflies to mate-searching males. Naturwissenschaften 87, 304307.[Medline]
Burghardt, F., Proksch, P. and Fiedler, K. (2001). Flavonoid sequestration by the common blue butterfly Polyommatus icarus: quantitative intraspecific variation in relation to larval hostplant, sex and body size. Biochem. Syst. Ecol. 29 (in press).
Chatfield, C. and Collins, A. J. (1980). Introduction to Multivariate Analysis. London: Chapman & Hall.
Cronin, T. W., Ruderman, D. L., Chiao, C. C. and Robinson, P. R. (1997). Human color vision and natural images: chromatic and spatial signals are not correlated. Invest. Ophthal. Vis. Sci. 38, S254S254.
Cuthill, I. C. and Bennett, A. T. D. (1993). Mimicry and the eye of the beholder. Proc. R. Soc. Lond. B 253, 203204.
Cuthill, I. C., Bennett, A. T. D., Partridge, J. C. and Maier, E. J. (1999). Plumage reflectance and the objective assessment of avian sexual dichromatism. Am. Nat. 153, 183200.
Darwin, C. (1859). On the Origin of Species by Means of Natural Selection, or the preservation of favoured races in the struggle for life, first edition. London: John Murray. (Facsimile, Cambridge, Mass.: Harvard University Press. 1964).
Douwes, P. (1976). Mating behaviour in Heodes virgaureae with particular reference to the stimuli from the female (Lepidoptera: Lycaenidae). Ent. Germ. 2, 232241.
Drummond, B. A., III (1984). Multiple mating and sperm competition in the Lepidoptera. In Sperm Competition and the Evolution of Animal Mating Systems (ed. R. L. Smith), pp. 291370. San Diego, CA: Academic Press.
Eguchi, E., Watanabe, K., Hariyama, T. and Yamamoto, K. (1982). A comparison of electrophysiologically determined spectral responses in 35 species of Lepidoptera. J. Insect Physiol. 28, 675682.
Endler, J. A. (1990). On the measurement and classification of colour in studies of animal colour patterns. Biol. J. Linn. Soc. 41, 315352.
Geuder, M., Wray, V., Fiedler, K. and Proksch, P. (1997). Sequestration and metabolism of host-plant flavonoids by the lycaenid butterfly Polyommatus bellargus. J. Chem. Ecol. 23, 13611372.
Ghiradella, H. (1991). Light and colour on the wing: structural colours in butterflies and moths. Appl. Opt. 30, 34923500.
Grether, G. F., Hudon, J. and Millie, D. F. (1999). Carotenoid limitation of sexual coloration along an environmental gradient in guppies. Proc. R. Soc. Lond. B 266, 13171322.
Grill, C. P. and Rush, V. N. (2000). Analysing spectral data: comparison and application of two techniques. Biol. J. Linn. Soc. 69, 121138.
Harborne, J. B. (1991). Flavonoid pigments. In Herbivores, their Interactions with Secondary Plant Metabolites, vol. 1, The Chemical Participants (ed. G. A. Rosenthal and M. R. Berenbaum), pp. 389429. San Diego: Academic Press.
Häuser, C. (1990). Vergleichende Untersuchungen zur Morphologie, systematischen Bedeutung und Phylogenese der weiblichen Genitalorgane der Lepidoptera (Insecta). PhD Dissertation. Albert-Ludwigs-Universität, Freiburg im Breisgau.
Häuser, C. (1993). Die inneren weiblichen Organe der Tagfalter (Rhopalocera): Vergleichende Morphologie und phylogenetische Interpretation (Insecta, Lepidoptera). Zool. Jb. Syst. 120, 389439.
Henderson, S. T. (1970). Daylight and its Spectrum. London: Adam Hilger Ltd.
Henderson, S. T. and Hodgkiss, D. (1963). The spectral energy distribution of daylight. Br. J. Appl. Phys. 14, 125131.
Honek, A. (1993). Intraspecific variation in body size and fecundity in insects: a general relationship. Oikos 66, 483492.
Hurlbert, S. H. (1984). Pseudoreplication and the design of ecological field experiments. Ecol. Monogr. 54, 187211.
Jacobs, G. H. (1992). Ultraviolet vision in vertebrates. Am. Zool. 32, 544554.
Jacobs, G. H., Neitz, J. and Deegan, J. F. (1991). Retinal receptors in rodents maximally sensitive to ultraviolet light. Nature 353, 655656.[Medline]
Kaiser, H. F. (1960). The application of electronic computers to factor analysis. Educational and Psychological Measurement 20, 141151.
Kim, J.-S. and Mueller, C. W. (1978a). Introduction to Factor Analysis. What it is and how to do it. Sage University Paper Series on Quantitative Applications in the Social Sciences, vol. 07013. Beverly Hills & London: Sage Publications.
Kim, J.-S. and Mueller, C. W. (1978b). Factor Analysis. Statistical Methods and Practical Issues. Sage University Paper Series on Quantitative Applications in the Social Sciences, vol. 07014. Beverly Hills & London: Sage Publications.
Knüttel, H. and Fiedler, K. (2001). On the use of UV wing patterns in butterfly morphology and systematics. J. Lepid. Soc. (in press).
Kornmaier, B. (1999). Einlagerung von Flavonoiden in die Flügelanlagen während der Puppenphase bei Polyommatus icarus (Lepidoptera: Lycaenidae). Diploma thesis. Universität Bayreuth.
Kramer, M. and Schmidhammer, J. (1992). The chi-squared statistic in ethology: use and misuse. Anim. Behav. 44, 833841.
Lombardi, C. M. and Hurlbert, S. H. (1996). Sunfish cognition and pseudoreplication. Anim. Behav. 52, 419422.
Martín Cano, J. (1984). Biología comparada de Lampides boeticus (L.), Syntarucus pirithous (L.) y Polyommatus icarus (Rot.) (Lep., Lycaenidae). Graellsia 40, 163193.
Menzel, R. (1979). Spectral sensitivity and color vision in invertebrates. In Handbook of Sensory Physiology, vol. VII/6A (ed. H. Autrum), pp. 503580. Berlin, Heidelberg, New York: Springer-Verlag.
Nijhout, H. F. (1991). The Development and Evolution of Butterfly Wing Patterns. Washington: Smithsonian Institution Press.
Pellmyr, O. (1982). Plebeian courtship revisited: Studies on the female-produced male behavior-eliciting signals in Lycaeides idas courtship (Lycaenidae). J. Res. Lepid. 21, 147157.
Sachs, L. (1997). Angewandte Statistik, eighth edition. Berlin: Springer-Verlag.
Schanda, J. D. (1997). Colorimetry. In Handbook of Applied Photometry (ed. C. DeCusatis), pp. 327412. Woodbury, New York: AIP Press, American Institute of Physics.
Schittko, U., Burghardt, F., Fiedler, K., Wray, V. and Proksch, P. (1999). Sequestration and distribution of flavonoids in the common blue butterfly Polyommatus icarus reared on Trifolium repens. Phytochemistry 51, 609614.
Sharma, V. K., Singaravel, M. and Subbaraj, R. (1999). Ultraviolet light-induced phase response curve for the locomotor activity rhythm of the field mouse Mus booduga. Naturwissenschaften 86, 9697.[Medline]
Silberglied, R. E. (1984). Visual communication and sexual selection among butterflies. In The Biology of Butterflies (ed. R. I. Vane-Wright and P. R. Ackery), pp. 207223. London: Academic Press.
Sokal, R. R. and Rohlf, F. J. (1981). Biometry, second edition. New York: Freeman.
StatSoft (1999a). STATISTICA für Windows. Tulsa, Oklahoma: StatSoft.
StatSoft (1999b). STATISTICA für Windows [Computer-Programm-Handbuch]. Tulsa, Oklahoma: StatSoft.
Stavenga, D. G., Smits, R. P. and Hoenders, B. J. (1993). Simple exponential functions describing the absorbance bands of visual pigment spectra. Vision Res. 33, 10111017.[Medline]
Tovée, M. J. (1995). Ultra-violet photoreceptors in the animal kingdom: Their distribution and function. Trends Ecol. Evol. 10, 455460.
Vane-Wright, R. I. and Ackery, P. R. (1984). The Biology of Butterflies. London: Academic Press.
Vane-Wright, R. I. and Boppré, M. (1993). Visual and chemical signalling in butterflies: functional and phylogenetic perspectives. Phil. Trans. R. Soc. Lond. B 340, 197205.
Wago, H., Unno, K. and Suzuki, Y. (1976). Studies of the mating behavior of the Pale Grass Blue, Zizeeria maha argia (Lepidoptera: Lycaenidae). I. Recognition of conspecific individuals by flying males. Appl. Ent. Zool. 11, 302311.
Wheeler, D. (1996). The role of nourishment in oogenesis. Annu. Rev. Ent. 41, 407431.
Wilson, A. (1987). Flavonoid pigments in chalkhill blue (Lysandra coridon Poda) and other lycaenid butterflies. J. Chem. Ecol. 13, 473493.
Wilson, D. S. and Dugatkin, L. A. (1996). A reply to Lomabardi & Hurlbert. Anim. Behav. 52, 423425.
Young, R. D. (1986). Principal-component analysis of macaque lateral geniculate nucleus chromatic data. J. Opt. Soc. Am. A 3, 17351742.[Medline]