When night never falls: female sexual signalling in a nocturnal insect along a latitudinal gradient

The environment can play an important role in animal communication by affecting signal transmission and detection. Variation in the signalling environment is expected to be especially pronounced in widely distributed species, potentially affecting how their signals are detected. Such environmental variability is presumably relevant for sedentary females of a nocturnal capital breeder, the European common glow-worm (Lampyris noctiluca), which produce green light during the night to attract flying males to mate. Being widely distributed in Europe, glow-worm populations are exposed to both rapidly descending, darker summer nights in the south, and slowly dimming, brighter summer nights further north, with the latter potentially posing challenges to the visibility of the female glow. To test how female signalling is affected by latitude, we sampled glowing females during summer nights along a latitudinal gradient in Finland, Northern Europe, and used a novel apparatus to measure the intensity and peak wavelength (hue/colour) of their glow. Surprisingly, females at higher latitudes, similar to those at lower latitudes, were commonly glowing during the brightest (and hence the shortest) nights of the year. Females also glowed brighter in more northern areas, partly due to their larger body size, whereas the colour of their glow was not associated with latitude. Since females glow even during midsummer, independent of latitude, the increase in glow intensity at higher latitudes presumably serves to maintain signal visibility in brighter signalling conditions. Overall, these findings highlight the influence of environmental conditions on the evolution of sexual signals, especially in the context of species distribution range. When environmental conditions impact signal transmission and perception, local conditions can have a crucial role in shaping animal communication and signal evolution. To analyse how dark-dependant common glow-worm females cope with variable nocturnal light environments, we used a novel apparatus, presumably not applied to living animals before, to measure female glow intensity at various latitudes along a latitudinal gradient. Interestingly, females did not avoid signalling during the brightest summer nights, but instead, their glow intensity and body size both increased with latitude. These findings suggest that females can ensure visibility to mate-searching males over a range of local conditions. Our study therefore shows how females can adapt to environmental constraints on signal visibility, and how the expression of sexual signals is shaped not only by social interactions but also by the signalling environment.


Introduction
Signals lie at the core of animal communication. They include morphological, physiological, and behavioural traits that are used to pass information between individuals (Otte 1974). In the context of sexual selection, signals may confer information about the signaller's quality and sexual availability to members of the opposite sex (Andersson 1994). Sexual selection favours signals that can most effectively deliver the message in the signalling environment and arouse the wanted response from potential mating partners, increasing the fitness Communicated by D. Kemp of the signaller (Endler 1992(Endler , 2000. Larger, brighter, and louder signals are the easiest to detect and often the most efficient at both active mate choice and passive mate attraction (Bateson 1983;Andersson 1994). The receiver's sensory and cognitive abilities, however, affect its mate choice processes, limiting the potential functions of the transmitted signal and, therefore, signal evolution (Ryan 1990). As an example of such limitations, human eyes cannot detect UV light, and hence, we fail to respond to any signals solely from the UV range. Moreover, signal processing may impact aspects of the receiver's sensory system, including its photoreceptors' spectral sensitivity (Endler 1992;Rowe 1999;Cronin et al. 2014). For instance, a signal consisting of a particular range of wavelengths may select for the receiver's spectral sensitivity to match the signal, allowing more efficient signal processing (a higher signal to noise ratio), as seen in some species of fish (Lythgoe and Partridge 1989;Bowmaker 2008), terrestrial vertebrates (Fleishman 1986;Ryan 1990), and insects (Alexander 1962;Seliger et al. 1982b;Bernard and Remington 1991;Cronin et al. 2000).
Besides the interactions between signallers and receivers, the signalling environment may also shape signal evolution by affecting signal detection (Endler 1993a(Endler , 2000Endler and Basolo 1998;Brumm 2004;Cole 2013). Signals travel from a sender to a receiver through a medium, such as air or water. The properties of the medium, together with the local conditions, are likely to affect whether the signal reaches the receiver, and if it does, how it is perceived (Endler 1990(Endler , 1991Gamble et al. 2003;Smit et al. 2019). Therefore, a certain signal property (for example, an acoustic frequency) may be conspicuous in one environment, but appears different to, or disappears altogether from, a receiver in another. In other words, signal conspicuousness and efficiency are relative measures (Endler 1992(Endler , 1993aJennions and Petrie 1997). When the signaller cannot efficiently control its signalling environment, it may increase its signal conspicuousness by adjusting the contrast between the signal and the signalling environment. This is possible by adjusting the signalling behaviour, for example, by choosing a signalling time or environment that maximises the signal to noise ratio, or by modifying the signal itself (Endler 2000). For instance, the effects of background noise can be minimised by signals that greatly differ from the signalling environment, as seen in Anolis lizards, which bob their head at a frequency different from the movements of the surrounding vegetation (Endler 1992). In the context of colours used in communication, highly represented wavelengths of ambient light can also be relevant to improve signal to noise ratio. For instance, a colour patch of an animal is perceived brighter when it matches the ambient light spectrum, increasing the overall contrast between the signaller and the environment. Similarly, a mismatch between the reflectance and ambient light spectra can result in a duller signal (Endler 1990(Endler , 1991. From a receiver's point of view, its spectral sensitivity may, in a relatively stable environment, evolve to match light conditions of the environment (Endler 1992;Leal and Fleishman 2002;McComb et al. 2010), which demonstrates the intimate links between the signaller, the receiver, and the signalling environment (Endler 2000;Cole 2013).
Light is a key environmental factor affecting living organisms. Depending on latitude, local populations are affected differently by the temporal distribution of light, for instance, in terms of the timing of reproduction, migratory behaviour, and mate attraction signals (Leal and Fleishman 2004;Hansson et al. 2007;Bradshaw and Holzapfel 2010;Shieh et al. 2017). During summers at high latitudes, nights are short and relatively bright, dusk and dawn are long, and ambient light consists of longer wavelengths than closer to the equator (Endler 1993b;Mills 2008). An individual close to London (N51.5°) at the beginning of July experiences about 2 h less of night-time than one close to Madrid (N40.4°) (as measured from the apparent sunset to the apparent sunrise; https://www. esrl.noaa.gov/gmd/grad/solcalc/). An individual to the north of N66.5°would, in turn, not experience any dark period at the end of June. Such variation in environmental light conditions may have a significant impact on visual signalling in nocturnal and twilight-active species (Biggley et al. 1967;Lall et al. 1980;Hall et al. 2016).
Most visual signals function by reflecting ambient light and thus are transmitted during daytime. Because visual signals have limited use for nocturnal animals, olfactory and acoustic signals are often their dominating signal types. However, the relationship between signals and the signalling environment is very relevant for animals that signal with bioluminescent light, such as fireflies and glow-worms (Coleoptera: Lampyridae) (Endler 1992). Most Lampyrid beetles are nocturnal or duskactive and exhibit large among-species variation in glowing properties (Lewis and Cratsley 2008). Notably, the peak wavelength (λ max ) of the glow ranges from 546 nm (green) to 635 nm (red) (Biggley et al. 1967;Hastings 1996), and to a large extent depends on the timing of the signaller's peak activity period (Seliger et al. 1964;Biggley et al. 1967;Lall et al. 1980;Hall et al. 2016). Nocturnal lampyrid species typically emit green light that matches the colour of the surrounding vegetation, resulting in a greater amount of the light being reflected from the vegetation, and therefore a greater likelihood of attracting a mate (Seliger et al. 1982a, b). The yellower (higher wavelength) glow of dusk-active species, in turn, is more advantageous during sunsets when the signaller is surrounded by the yellower light of the setting sun. Under such conditions, the yellower glow contrasts with the wavelengths reflected from the surrounding vegetation, resulting in a pronounced signal (Lall et al. 1980;Hall et al. 2016). In this respect, the wavelength variation of bioluminescence in fireflies has presumably developed to enhance the visibility and efficiency of their sexual signals in different signalling environments (Lall et al. 1980;Hall et al. 2016). Thus far, it has remained unknown whether glow intensity within a species, as well, is adjustable to geographical variation in light conditions.
As light conditions change with latitude, populations of widely distributed species may be expected to adjust their signals to local light environments to ensure signal visibility. We addressed this question using a nocturnal beetle, the European common glow-worm (Lampyris noctiluca L., Coleoptera: Lampyridae), in which flightless females glow during summer nights to attract flying males (Schwalb 1961;Tyler 2002). Their distribution in Europe extends as far north as central Finland (N64°) (information from the Finnish biodiversity info facility, https://laji.fi/). Light conditions may affect detection of female glow, with females needing to attain an adequate contrast with their signalling environment to be visible to mate-searching males. In particular, summer nights at higher latitudes are notably brighter and shorter (for instance, ambient light in the middle of the night can be as high as 5 lx, bright enough to allow reading during the~4-h daily dark period in June at N62°), and the ambient light also consists of longer wavelengths than in more southern regions. It is expected that such conditions at higher latitudes require brighter glow and/or glow with longer wavelengths (more yellow, as during dusk) for the signal to reach similar detectability as glow signals emitted under light conditions of lower latitudes (Lall et al. 1980;Seliger et al. 1982a, b;Hall et al. 2016).
The aim of this study was to investigate how adult female glow-worms cope with signalling at different latitudes and their light conditions. We expected that brighter summer nights at higher latitudes select for (1) a later emergence of females in summer to avoid glowing when the nights are the brightest, and/or (2) an increased female glow intensity and peak wavelength (λ max , a measure of glow hue, or more generally speaking, colour). Each of these traits can be expected to increase the visibility of females to mate-searching males. We tested these hypotheses by locating glowing females and by measuring their glow intensity, λ max , and ambient light intensity, along a latitudinal gradient in Finland.

Study species
The common glow-worm is widely distributed in the northern hemisphere, ranging between latitudes N31°(central Iran) and N64°(central Finland) (Samin et al. 2018; the Finnish biodiversity info facility, https://laji.fi/). In many areas within this range, including England (Hickmott and Tyler 2011) and southern Finland (personal observations), adults are most active between June and July when summer nights are at their shortest and brightest. During the narrow timeframe when the sun is set, flying males search for females, with the brightest, and thus largest (Hopkins et al. 2015) females having the highest probability of attracting a mate (Tyler 2002;Hopkins et al. 2015;Lehtonen and Kaitala 2020). Females produce their mate attraction glow with chemical chain reactions within a light organ in the tip of their abdomen (Tyler 2002). Once a female has reached its maximum brightness, it glows constantly throughout the dark period (Tyler 2002). The glow has a yellowish green colour, peaking at 546-570 nm in wavelength (Booth et al. 2004;De Cock 2004;Bird and Parker 2014). Delayed mating is a considerable cost for females because they do not eat as adults (Tyler 2002) and lose eggs until succeeding to mate (Wing 1989;Hopkins 2018). Therefore, good visibility of their mate attraction signal is crucial.

Measurements
We measured female glow and ambient light at females' chosen glowing spots in June and July of 2017 and 2018. In total, we sampled 11 sites within four separate study areas, along a latitudinal gradient in Finland ( Fig. 1). We started searching for glowing females at each site late in the evening, but before the first ones began to glow. Independent of latitude, females typically glowed in open areas, either on a rocky surface (often partly surrounded by lichen, moss, or other vegetation) or hanging from a blade of grass. In 2017, we measured females in the southern sites (study areas A and B in Fig. 1) between the 16th of June and the 9th of July, and in northern sites (study areas C and D in Fig. 1), where the growing season starts later, between the 10th and the 31st of July. In 2018, we found and measured females starting from the southernmost study sites on the 3rd of June, and then continued north, keeping up with the progression of summer at different latitudes. We arrived the northernmost study sites on the 24th of June to search for and measure females until the 7th of July, after which we returned to the southern sites and continued measurements there (Table 1). The dark period was shorter, and dawn and dusk lasted longer further north (personal observations, https://www.esrl.noaa. gov/gmd/grad/solcalc/). Thus far, it has been challenging to collect glow intensity and colour data without handling the insect during the measurements. Moreover, previous measurement procedures have most likely included a considerable risk of measurement biases due to extraneous light (Hall et al. 2016). We successfully introduced a novel method to measure the absolute irradiance of female glow in field conditions (from here on called "glow intensity"). This method allowed us to measure glow intensities of females in their natural mate attraction pose with only a little or no handling of the females, minimising environmental light biases. In particular, we recorded the glow spectra of 85 females in the field over the 2 years (2017: N = 34, 2018: N = 51) using an integrating sphere 50 mm in diameter (IC-2, StellarNet Inc., USA) fibre-optically connected to a spectrometer (Flame-S, Ocean Optics Inc., USA), controlled by the operating software OceanView (version 1.5.2, Ocean Optics Inc., USA) on a laptop. The integrating sphere with its barium sulphate coating ensures total integration and mixing of the light entering the sphere, producing a highly diffuse reflectance (Fig. 6 in Endler 1990). Its small opening for the object to be measured also blocks direct ambient light from entering the device. These measurements were done on a 24 cm × 32 cm cardboard platform that was covered with a sheet of white paper and had a 3.0-cm-long and 0.5-cm-wide white hollow plastic stick protruding through the paper and pointing upwards (Fig. 2).
When we located a glowing female, we placed it gently on the platform next to the stick. If the female did not climb onto the stick, we gently placed the female directly on it. When we found a female that was glowing on grass, moss, or leaf, we took some of the plant material with the female to avoid touching and disturbing it. If it was not practical to remove some of the material the female was glowing on, we moved the female to the platform gently with a plastic spoon. When attracting males, a female raises its abdomen upwards and stays rather still. Therefore, to minimise the possibility of stress from handling, we proceeded with the measurement only after the female had settled down and raised its abdomen. We then lowered the integrating sphere ( Fig. 2) so that the female only just disappeared inside (its glow not directed to the photodetector of the collecting optic fibre to avoid measurement error) and initiated the measurement. We discarded the recordings in which the female lowered its abdomen during (or just before) the measurement, as detected with the OceanView software visualisation on the laptop screen, and took a new one.
Each measurement included six scans with an integration period of 0.5 s (the time over which the spectrometer gathers light during one scan), and the average over the six scans was used in the analyses. When the female glow was too dim to produce a detectable emission spectrum, we increased the integration time to 0.8 (N = 1), 1.0 (N = 5), or 2.0 s (N = 1). If the glowing intensity was at the border of the equipment's accuracy, we were not able to reliably measure the glow of the individual. This happened mostly at the start of the study period in June in the southernmost study area, likely because night temperatures Fig. 1 Glow-worm collection and study areas in Finland. In 2017, we collected the data from study areas A and D, and from all the four study areas in 2018. The distance between the farthest areas is circa 350 km. The data for creating the map was extracted from GADM database (www. gadm.org), version 3.4, in April 2018 were low, and the female glowing capacity may have been challenged. After the glow measurement, we measured the width of the back shield (the pronotum) with a Vernier calliper, as it is a good proxy of female body size (Hopkins et al. 2015;Borshagovski et al. 2019). We measured each female only once and then returned it to its original spot of collection. The measurement of the ambient light where the female was found was conducted immediately after the measurement of the pronotum. We used the same devices as with female glow measurements but replaced the integrating sphere with a cosine-corrected irradiance probe (CC-3-UV-S, Ocean Optics Inc., USA), which gathers light from a 180°angle (Fig. 5 in Endler 1990; Online Resource 1, Fig. S1). The spectrometer was oriented so that the cosine-corrected irradiance probe was pointing upwards at the spot where the female had glowed. We positioned the laptop and measurers as far away from the measuring spot as allowed by the 1.5-m cable connecting the laptop and spectrometer, to avoid light reflecting from the measurers and affecting the measurements. Each ambient light measurement included three scans with an integration time of 3 s, from which we set the programme to calculate the average values. If the ambient light was too dim to produce a detectable emission spectrum, for instance when a female was glowing in a dense bush, we increased the integration time to 4.0 (N = 12), 5.0 (N = 4), 6.0 (N = 2), or 10.0 s (N = 1). Such sensitivity adjustments are not expected to impact the result of the measurement. We successfully measured the ambient light level N = 75 times. We did not measure it when it was too dark for our measurement device to produce a reliable measurement. In addition, when the female glowed close to a streetlamp (N = 6), ambient light emission spectra peaked at significantly higher wavelengths than under natural light, and we therefore excluded these ambient light measurements from the analyses.
For the analyses of light data, we produced a script for MATLAB R2016b (the MathWorks Inc., USA). In brief, each data pair of female glow and ambient light was processed as follows. Each measurement was assessed between 300 and 700 nm, averaged and low-pass filtered (Savitzky-Golay, S-G), and depicted in emission spectrum figures (Fig. 3 for the female spectrum, and Fig. S2 in Online Resource 1 for the ambient light spectrum). If some of the six scans of a female resulted in notably lower glow intensity values than the rest (potentially biasing the intensity measurement result), we discarded them, one at a time, until the average of the remaining intensity values equated ≥ 85% of the highest intensity value for that female (N = 6). To determine an accurate peak location of the female glow spectrum (λ max , measure of hue/ colour), a skewed Gaussian fitting was applied. Absolute irradiance of the glow in photons/s, and the ambient light in photons/s (also in lux, to enable comparisons with previous findings) were calculated from the surface area of the Gaussian fit and over the 300-700-nm S-G-filtered average spectrum, respectively (Online Resource 2).
We calibrated the integrating sphere before field measurements with a light source, LS-1-CAL (Ocean Optics Inc., USA). The cosine-corrected irradiance probe was calibrated by Ocean Optics with a standard light source. Ideally, the equipment should be calibrated before every measurement. However, due to the logistic challenges it would have caused under field conditions, and its time-consuming nature, we first conducted several test measurements by reassembling the measurement equipment with a calibration against a stable light source to see if any significant differences would appear in the results. Because we did not find any significant shifts in the results, we produced all the measurements for both study years with the same calibration, conducted just before the 2017 field season.
Overall, the dataset for our analyses included 85 female glow intensities and λ max values, 75 ambient light intensities, 74 ambient λ max values, and 81 female size measurements. In addition, for each of these measurements, we also retrospectively assigned a temperature value from the data of the Finnish Meteorological Institute (available for each hour at https://en.ilmatieteenlaitos. fi/). Because temperature values are not additive quantities (in relation to female glow intensity and λ max ), we used temperature as a categorical variable in the data analysis by classifying the temperature values at the time of glow measurements into four categories, (1) < 10°C, (2) 10-13.9°C, (3) 14-17.9°C, and (4) ≥ 18°C (the temperatures for each site are given in Table 1).

Statistical analyses
We conducted all analyses with R software version 3.6.1 (The R Foundation for Statistical Computing, http://www.r-project. org). To analyse the factors affecting female glow (intensity and peak wavelength), we performed two linear models (LM).
For this purpose, we tested the data for normality with Fig. 3 Peak wavelength (λ max ) (nm) and total irradiance (photons/s) were calculated with MATLAB R2016b (the MathWorks Inc., USA) from the Skewed Gaussian fitting curve (thick, broken curve). Average μW/nm values (thick, black curve; μW/nm is the default unit of the measurement software) were calculated and drawn from the six scans (thin, saw-toothed curves) taken from a single female glow-worm with 500-ms integration time. Small differences between the intensity curves are a result of the female moving slightly inside the integrating sphere during the measurement. Glow spectrum ranged typically between 500 and 700 nm, with circa 70-nm full width at half maximum Shapiro-Wilk tests, and applied logarithmic transformations when normality was not otherwise met. In addition, we visually checked model assumptions of normally distributed errors and homoscedasticity with quantile-quantile plots of residuals and scale-location plots. The first model included female glow intensity (logarithmic) as the response variable, and latitude, female body size (pronotum width), ambient light intensity (logarithmic), year, and temperature (categorical) as explanatory variables. The initial model also included all two-way interactions between the explanatory variables. We then simplified the model with backward elimination of nonsignificant interactions, having P = 0.10 as the threshold. The model objects produced by the final linear models were assessed with the ANOVA function in R. The second model included peak wavelength (λ max , its hue/colour) of female glow as the response variable, and latitude, female body size, year, and temperature category as explanatory variables. In other respects, the model was treated as above. To further assess any relationships between our key variables, we conducted a correlation test between the two response variables, latitude and female body size, and a partial correlation test between latitude and glow intensity while controlling for female body size.
We did not use dates of glow measurements and times of the day in the models because environmental conditions and developmental stages of glow-worms are likely to be very different between southern and northern study areas at a certain time. As much as possible, we attempted to collect data during the peak glowing period, independent of latitude. Measurement dates and times are provided in Table 1.

Light conditions along the latitudinal gradient
The ambient light conditions were highly different for females glowing at the different locations in Finland. The intensity of ambient light (over both study years) during glow measurements was circa nine times higher in the northernmost study area (N62.62°; 1.02 × 10 13 ± 1.94 × 10 12 photons/s, 2.40 ± 0.46 lx [mean ± SE], N = 13), than in the southernmost study area (N59.86°; 1.16 × 10 12 ± 1.59 × 10 11 photons/s, 0.27 ± 0.039 lx, N = 45). The intensity of the brightest ambient light, in which we found glowing females, was 2.29 × 10 13 photons/s (5.37 lx). Despite the challenges ambient light conditions pose to the visibility of female glow in higher latitudes, we found females glowing in the northernmost study area even on the 24th of June, the brightest time of summer (Table 1). Hence, the results suggest that females do not avoid glowing during the period of the brightest summer nights.

Female brightness
All the considered interaction terms were nonsignificant (P > 0.10). The final model with the main effects showed that female glow intensity increased significantly with both latitude and female size (Table 2; Figs. 4 and 5), indicating that larger females were brighter, and that females were brighter (relative to their size) in the north. Indeed, when controlling for female body size, there was a statistically significant positive partial correlation between latitude and female brightness (Pearson r partial = 0.276, t 78 = 2.536, P = 0.013). Female size also increased with latitude (Pearson correlation test: r = 0.367, N = 81, P < 0.001; Fig. 6). Over all study sites, the mean female glow intensity in 2017 was 3.50 × 10 12 ± 5.36 × 10 11 photons/s (range: from 2.57 × 10 11 to 1.50 × 10 13 photons/s, N = 34), and in 2018, it was 3.03 × 10 12 ± 2.79 × 10 11 photons/s (range: from 7.8 × 10 11 to 8.44 × 10 12 photons/ s, N = 51). Glow intensities did not significantly differ between the 2 years (Table 2). In addition, neither intensity of ambient light nor temperature significantly affected the intensity of female glow ( Table 2).

Discussion
We found that female glow-worms were glowing in our most northern study area (N62.62°) during midsummer nights, which suggests that they attract males even during the brightest time of the year, in the part of their range where summer nights are the least dark. Our results also show that glow intensity and body size of females increased with latitude (that is, towards north). Larger body size and greater glow intensity at higher latitudes presumably allow sufficient visibility to mate-searching males in the challenging signalling environment (relatively bright nights). Given that both female size and other factors that increase glow intensity are expected to increase mate attraction success (e.g. Hopkins et al. 2015;Lehtonen and Kaitala 2020), it is reasonable to ask what prevents females in the lower latitude sites from growing to a larger size and producing a brighter glow, similar to females at the northern sites. Most  Fig. 7 The average peak wavelength (λ max ) of female glow in each temperature category ((1) < 10°C, (2) 10-13.9°C, (3) 14-17.9°C, and (4) ≥ 18°C) with standard errors, including data from both study years importantly, larger size and other potentially costly factors that help to attain a brighter glow should provide a higher selection advantage in the brighter light conditions of higher latitudes. However, growing seasons are also shorter in the north, which makes larger body size potentially even costlier to reach there than in more southern regions. Hence, individuals in higher latitudes seem to experience a steeper trade-off between the costs of growing (due to the shorter growing seasons) and visibility (larger body size needed). Given that there is a clear relationship between development time and body size in many species, one way of reaching a larger adult body size is to prolong development time (Roff 1992(Roff , 2000, with development time and body size often correlating positively (Masaki 1967;Roff 1992). Therefore, insects are expected to prolong their development time, when feasible, to gain larger size and greater fecundity (Honêk 1993). At higher latitudes, prolonging the time available for growth may require spending at least one additional year in the larval stage, with an additional risky winter diapause. In other words, females presumably face a trade-off between the costs of postponing maturation (by at least one full year so that pupation occurs during the summer breeding season) and gaining a larger body size during that time, which ensures higher glow intensity and visibility to mate-searching males. Our results indicate that the way females respond to this trade-off depends on the latitude of their signalling site. Some animals can compensate for the shorter available development time at higher latitudes by growing at a faster rate (Abrams et al. 1996;Arendt 1997;Nylin and Gotthard 1998). Indeed, seasonal environments favour plasticity in developmental time (Nylin and Gotthard 1998). However, a shorter time available for the development does not typically result in greater body size (Abrams et al. 1996), and therefore does not explain our result of females being larger at higher latitudes. To gain further insights into the mechanisms that cause the positive association between adult female size and latitude, and to better understand the complex relationship between body size, development time, and growth rate, we encourage future studies that involve rearing glow-worms from populations located at a range of different latitudes under a common garden setup.
Similarly, the timing of maturation within a summer is likely to be subject to trade-offs, especially regarding the properties of the signalling environment, mate availability, and temperatures suitable for maturation and mate attraction. For instance, in many species, the first individuals to emerge as adults face decreased competition, but these individuals are often smaller in size (Morbey and Ydenberg 2001), which may be costly to attractiveness or fecundity. In glow-worms, a small body size has multiple costs to females. First, small females have a dimmer glow, and they therefore attract males less easily (Hopkins et al. 2015). Second, fecundity costs due to delays in mating seem to be higher for small females of these capital breeders (Hopkins 2018). Hence, small size poses a risk of not attracting a mate quickly, and this delay in mating may even result in a complete reproductive failure quicker than it does for larger females. Theoretically, these costs could, at least partly, be mitigated by postponing maturation to a time after the brightest midsummer nights. The fact that females did not do so, even at the higher latitude sites, suggests that it is not feasible. For instance, postponing maturation at higher latitudes might leave too small a timewindow to ensure successful mating, before the short season is over. Therefore, factors other than the challenges of glowing under bright ambient light conditions seem to be the main drivers of the timing of maturation.
Although female glow was brighter at higher latitudes, where nights were more illuminated, female glow intensity was not significantly linked to the intensity of ambient light at the measurement spot. In other words, within the study areas, individual females did not adjust their glow intensity to the immediate ambient brightness. Therefore, glow intensity does not seem to be a trait that is quick to control at an individual level. Previous findings are in accordance with this finding: females do not adjust their glow intensity between nights to respond to the time spent unmated, and competition does not affect the onset of glow (Baudry et al. unpublished data). Nevertheless, it remains unclear to which extent increasing glow intensity and body size towards higher latitudes are genetically based or represent phenotypic plasticity, as both types of responses can drive latitudinal body size differences in ectotherms (e.g. Mousseau and Roff 1989;James et al. 1997;Blanckenhorn and Demont 2004).
Interestingly, even when body size was controlled for, glow intensity increased with latitude, indicating that body size is not the only factor driving glow intensity, and that females have additional means to glow brighter at higher latitudes. Currently, the exact mechanism allowing females to do so is not known, providing fascinating avenues for further investigation, preferably using laboratory-reared individuals. We also note that, besides night brightness, other factors may differ with latitude, potentially contributing to the development and quality of the glow signal. We do not have data on the variation in soil quality or food availability, for instance. Therefore, while we show that the female phenotype (glow and size) matches with the requirements of the night-time conditions that are defined by latitude, further research is needed to investigate the full range of factors that might contribute to this pattern.
In contrast to the glow intensity, the peak wavelength (λ max , measure of hue/colour) of female glow did not covary with latitude. We offer the following two hypotheses that may explain the lack of a significant relationship. First, the effect of latitude on colour may become apparent only in a very large-scale study, that is, when measuring female glow spectra over a wider latitudinal range. In this regard, it is worth noting that although the latitude range used in the current study did not reveal a significant gradient in glow colour, it did so regarding female size and brightness. Second, if male mate searching is not sensitive to small differences in λ max values, the exact colour of female glow may not be under selective pressure, indicating that latitudedependent changes in ambient wavelength composition are less important in sexual signalling. Glow-worms may, indeed, adjust to ambient light more easily, or efficiently, in other ways than by changing glow colour, for example by altering some aspects of their signalling behaviour.
The temperature category did not significantly explain variation in female brightness. However, in accordance with glow intensity increasing with temperature to a certain optimum in in vitro experiments (Ueda et al. 1994), we did notice that females could not glow at their usual intensity during the coldest nights (approximately below 7°C). In contrast to brightness, the λ max of females' bioluminescent light was linked to temperature. This is also in accordance with earlier in vivo and in vitro experiments showing that an increase in temperature shifts the λ max towards longer wavelengths, that is, towards the yellow colour (Hastings 1996;Rabha et al. 2017). As with the colour of the glow-worm glow, firefly flash signalling patterns-which have a significant role in mate attraction and species recognition (Lewis and Cratsley 2008)-are also affected by temperature (Sharma et al. 2016;Rabha et al. 2017). Interestingly, despite the mean temperature being higher in 2018, the average λ max was significantly lower (just the opposite of what the temperature difference would have led us to expect). This could be a result of other factors besides temperature (potentially involving physiological restrictions on light production) affecting the λ max . We note that the small (~1 nm) between-year difference in λ max is unlikely to affect mate search or have other biological relevance. While some day-active insects may be able to discriminate such fine-tuned colour differences (Koshitaka et al. 2008;Telles et al. 2016), the ability for such precise colour discrimination is unlikely for night-active species. Moreover, the spectrum width of the female glow is so wide that small peak differences are even less likely to be detectable.
In conclusion, female sexual signalling was found to be related to its environmental light conditions. In addition, females in the most northern sites were found to attract males even during midsummer nights, despite the challengingly bright night-time conditions. Hence, our results emphasise how a species' signalling environment can promote regional differences in life-history traits and sexual signals, which, in turn, may contribute to population differences through sexual selection. These findings also highlight the potential of nocturnal animals that engage in visual signalling to partially adapt to the increasingly common light pollution over time.
Author contributions AMB and AK designed the experiments. AMB, assisted by TKL, performed the field experiments. PS purchased and prepared the measuring devices, developed the code for MATLAB, and assisted with the use of MATLAB. AMB and TKL carried out the statistical analyses. AMB wrote the first draft of the manuscript with all authors contributing to its preparation.
Funding Open access funding provided by University of Oulu including Oulu University Hospital. The work was funded by the Academy of Finland (grant number #294664 to AK, TKL and AMB) and Emil Aaltonen Foundation (AMB).
Data availability The data generated and analyzed for this study is available in the electronic supplementary materials.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.