Structural plumage colour is one of the most enigmatic sexually selected traits. The information content of structural colour variation is debated, and the heterogeneity of the findings is hard to explain because the proximate background of within-species colour differences is very scarcely studied. We combined measurements of feather macrostructure and nanostructure to explain within-population variability in blue tit crown reflectance. We found that sexual dichromatism in aspects of crown reflectance was explained only by feather macrostructure, whereas nanostructural predictors accounted for some of the age-related differences in reflectance. Moreover, we found that both mean reflectance and spectral shape traits reflected a combination of quantity and regularity aspects in macrostructure and nanostructure. This rich proximate background provides ample scope for reflectance to convey various types of information on individual quality.

In the face of directional selection, sexual signal traits may be constrained by various costs and limitations. Production costs arise in connection with the development of the trait, such as the energy expenditure required for display (Hedenström, 1995) or the necessity of collecting and depositing pigments into some integumental colours at the expense of other physiological functions (Olson and Owens, 1998). In contrast, wearing costs are caused by the trait itself, for example, the additional predation risk imposed by visual conspicuousness (Slagsvold et al., 1995) or the negative effect of dragging a long ornament on flight efficiency (Thomas, 1993). Furthermore, intrinsic individual characteristics may limit the production of certain signal character states. For example, stress and impaired neural development early in life may limit song learning capability in birds (MacDonald et al., 2006), and body size may limit the amplitude and frequency of acoustic signals in amphibians (Hoskin et al., 2009).

Bird plumage colour is a conspicuous example of potential signal traits and has therefore attracted much research (for reviews of non-avian colour, see e.g. Kinoshita and Yoshioka, 2005; Caro et al., 2017). Although there is a classical categorization of plumage colour as pigment-based versus structural, it currently seems that, with very few exceptions, colour variation is generated by the interaction of pigments and structures (Shawkey and D'Alba, 2017). The keratin matrix always contributes to the expression of pigment-based colour (Shawkey et al., 2006a), and pigments play vital roles in the modulation of structural colour (Stavenga et al., 2011; Maia et al., 2012a; Tinbergen et al., 2013). The pigment-based component of colour usually involves carotenoid or melanin deposition, and the honesty of the resulting signal can be derived from the costs of obtaining, synthesizing or metabolically transforming the pigment (Griffith et al., 2006; Guindre-Parker and Love, 2014).

One type of structural colour is achromatic white reflectance produced by incoherent scattering of incident light (e.g. McGlothlin et al., 2007). In this colour type, the quantitative aspects of the feather keratin macrostructure may explain trait expression and honesty (Prum, 2006; Shawkey et al., 2006a). Another type is iridescent structural colour, where the dominant wavelength changes with the relative angles of incident light and observation (Osorio and Ham, 2002), and is typically produced by 2D photonic crystals or multi-layer thin-film interference in feather barbules (Yin et al., 2006; Maia et al., 2009; Eliason and Shawkey, 2012; Lee et al., 2015; Stavenga et al., 2017). Finally, many avian taxa exhibit non-iridescent chromatic structural colour. This is produced by constructive interference of short wavelengths by the spongy keratin layer of feather barbs, which shows short-distance structural regularity (Prum et al., 1999; Saranathan et al., 2012).

The honesty of chromatic structural colour may involve at least two limitations (Fitzpatrick, 1998). The first is the production cost of keratin and melanin required for the nanostructure (Meadows et al., 2012). The second is the accurate self-assembly of the ordered nanostructure necessary for constructive interference (Prum et al., 2009; Maia et al., 2012a). Therefore, predictions concerning the information content (e.g. condition dependence) of structural colouration depend on how feather structural properties translate to aspects of reflectance variation at the within-species level. However, despite extensive data at the species level (e.g. Saranathan et al., 2012), studies of within-species variation of colour and nanostructure are currently restricted to one species each for iridescent (satin bowerbird, Ptilonorhynchus violaceus; Doucet et al., 2006) and non-iridescent structural colour (eastern bluebird, Sialia sialis; Shawkey et al., 2003, 2005), which allows little generalization.

In addition to nanostructure, feather macrostructure, i.e. the size and density of coloured feather parts, may also contribute to colour expression. Macrostructure is certainly not the mechanistic basis of the overall colour, but it may still grossly modify the visual attributes of the given individual, principally through quantity and density aspects of the reflective surfaces. Macrostructure of contour feathers has been shown to be condition dependent (Vágási et al., 2012) and environment dependent (common garden experiment; Broggi et al., 2011). Moreover, in the only experimental study that measured feather nanostructure, manipulation of food access at moulting impaired feather macrostructure but not nanostructure (D'Alba et al., 2014). Still, only a few studies have correlated contour feather macrostructure and reflectance traits, and nearly all of these have focused on dark melanin-based colour produced by pigment-containing barbules in high densities (Galván, 2011; D'Alba et al., 2014; sex and morph differences in Enbody et al., 2017; but see one macrostructural variable in relation to structural colour in Shawkey et al., 2005). No general conclusion can therefore be drawn on the macrostructural mechanisms of colour determination but, according to the studies mentioned above, the surface (size, density, etc.) of the colour-producing feather parts may play an important role.

Here, we examine some potential nanostructural and macrostructural predictors of individual reflectance variation in a well-known structural plumage ornament, the crown of the blue tit (Cyanistes caeruleus). The UV–blue reflectance of this trait has been suggested to exhibit strong sexual dichromatism (Andersson et al., 1998) and age dependence (Delhey and Kempenaers, 2006), and be subject to intersexual (Hunt et al., 1999) and intrasexual (Rémy et al., 2010; Midamegbe et al., 2011) selection in both sexes. Furthermore, the expression of the trait in males appears to mediate vocal interactions among neighbours (Poesel et al., 2007), and perhaps also sex ratio adjustment by the female partner (Sheldon et al., 1999; Delhey et al., 2007). Finally, both sexes have been found to react to the partner's crown colour at parental care, although in varying directions (Limbourg et al., 2004, 2013; Mahr et al., 2012). Concerning the information content of blue tit crown colour, its age-related increase seems to reflect within-individual change between years (Delhey and Kempenaers, 2006), and crown colour expression has been found to be influenced by moult speed manipulation (Griggio et al., 2009) and brood size manipulation (Doutrelant et al., 2012), although not by mild nutritional manipulation during moulting (Peters et al., 2011). Therefore, given that different measures of blue tit crown colour have been shown to indicate phenotypic quality and to function as signals in both sexes, the species is ideal for examining the proximate background of variation in non-iridescent structural colour expression. In addition to individual differences, we also investigate whether the sex and age dependence of crown reflectance can be traced back to feather nanostructural properties.

Field methods

Among other hole-breeders, we have been studying a breeding population of blue tits [Cyanistes caeruleus (Linnaeus 1758)] in our nestbox plots in the Pilis–Visegrádi Mountains, Hungary. For the present study, we used data from the year with the largest sample size (2009). In this year, all parents were caught at nestling feeding, and approximately 10 crown feathers were collected in envelopes for spectral measurements in the laboratory. This study was conducted under a long-term research agreement with the Pilis Park Forestry (December 1988 and March 2007), and with a research permit from the regional nature conservation authority (KTVF 35464-1/2007). Blind to the reflectance values of their feathers, we chose 32 individuals for feather structure analyses. These included eight juvenile females, eight adult females, eight juvenile males and eight adult males.

Feather macrostructure

The crown feather was placed on a white sheet with a piece of millimetre scale paper attached for calibration, and it was gently pressed down with a microscope slide. Two feathers of each individual (one at a time) were photographed in RAW format under diffuse natural illumination using a Nikon D70 camera and AF Micro Nikkor 60 mm lens. The resulting photographs (Fig. 1A) and additional scanning electron microscope imaging (Fig. 1B) indicated that there are no or very few barbules (i.e. lower order branches) in the apical part of the feather. The UV–blue colour is probably produced principally by the bare distal sections of the apical barbs, which indeed look blue in the photographs. The rest of the feather is grey and densely covered by dark (likely melanin-containing) barbules. After size calibration, two objective measurements were made with a distance tool using SPIP (Image Metrology Inc., Hørsholm, Denmark). We measured the distance from the end of the rachis to the end of the most distal blue barb (hereafter, overhang length), and the distance (parallel to the final section of the rachis) between the most proximal bare point of the most proximal blue barb and the end of the feather (hereafter, blue section length). We also counted the number of blue barbs (with bare sections longer than 1 mm) in the feather (hereafter, number of barbs). All three measurements were repeatable for the two feathers of the same individual (Pearson’s r=0.763–0.815, P<0.001).

Fig. 1.

Structure of blue tit crown feathers. (A) Overall macrostructure, (B) SEM image of the base of the colour-generating bare barb sections (photograph courtesy of Károly Bóka) and (C) TEM image of nanostructure. a, air vacuole; c, barb cortex; m, melanosomes; s, spongy layer.

Fig. 1.

Structure of blue tit crown feathers. (A) Overall macrostructure, (B) SEM image of the base of the colour-generating bare barb sections (photograph courtesy of Károly Bóka) and (C) TEM image of nanostructure. a, air vacuole; c, barb cortex; m, melanosomes; s, spongy layer.

Feather nanostructure

One undamaged crown feather of each of the individuals was prepared for transmission electron microscopy (TEM). Feather barbs were cut from the upper 5 mm of contour feathers and incubated in 0.25 mol l−1 sodium hydroxide and 0.1% Tween-20 for 30 min on a bench-top shaker. These barbs were then changed into 2:3 (v/v) formic acid and ethanol for 2.5 h. Feather barbs were dehydrated by incubation in 100% ethanol for 2×15 min and 100% propylene oxide for 2×15 min, and infiltrated with Spurr Low-Viscosity resin (EM0300, Sigma-Aldrich) in successive concentrations of 25%, 50%, 75%, 90% and 100%. Each of these infiltration steps was performed for 24–48 h except the last one, which was only 6 h. Barbs were placed into moulds with the most distal tip of the barb at the top of the mould, and then the blocks were cured in an oven at 80°C for 48 h. Sections of 80 nm thickness were cut oriented perpendicularly to the barbs of the blue apical portion of the feather using a diamond knife on an Ultracut E ultramicrotome (Reichert-Jung, Leica Microsystems, Cambridge, UK). Sections were placed on a 150-mesh copper grid with Formvar support, post-stained in 2.5% water-based uranyl acetate for 10 min and Reynold's lead citrate for 3 min, and photographed on a Jeol JEM1011 transmission electron microscope (Jeol, Tokyo, Japan).

For one feather, the preparation for TEM was unsuccessful, so this individual was omitted from the dataset. Two undamaged barb cross-sections (4000× or 5000× magnification, see Fig. 1C) were chosen for analysis from each feather. The images were processed in SPIP, and the following nanostructural parameters were measured. Thickness of the homogeneous outer keratin matrix and the distance between the cortex and the melanosomes beneath the spongy layer (keratin–air layer beneath the cortex) were recorded five times per cross-section using a distance tool (yielding cortex thickness and spongy layer thickness, respectively). Spongy layer thickness was not recorded in a subset of the individuals (N=9) where the interior air vacuoles beneath the spongy layer were not filled by the resin, and therefore the melanosomes were displaced and dispersed by the ultramicrotome.

Two magnifications of a fixed scale (25,000×) and a fixed size (corresponding to 1.5 μm picture length and width) were prepared from the spongy layer of each cross-section. In these detail images, after noise reduction, we measured the density of the spongy layer (percentage of the picture area covered by keratin) using a detection tool (by manually adjusting the detection threshold value). We also manually measured the minimum diameter of 20 approximately circular keratin rods and air spaces using a distance tool. Twenty was the approximate number of circular rods in a detail image, and fewer than 20 rods were found in some images. We measured the minimum diameter as the spongy layer is a branching structure in which the exactly perpendicular intersection of rods or cavities is unlikely, and therefore the minimum diameter best approximates the true diameter. We only focused on circular rods and air spaces because branching ones probably correspond to the ‘joints’ of the air or keratin network, and are therefore not representative of gross nanostructure. From these data, we calculated mean and standard deviation (s.d.) values for keratin rod and air space diameter. Finally, ‘distance between scatterers’ (i.e. surfaces with a shift in refractive index) was estimated by summing the mean rod and air space diameters. This measure indicates periodicity in the spongy layer. Although it is a derived measure and therefore correlated with its component traits, we included it for comparison with the findings of Shawkey et al. (2005).

The repeatability of these measurements could be tested on different levels. Cortex and spongy layer thicknesses were first subjected to repeatability testing at the cross-section level (five data per cross-section) using the rptR procedure (Nakagawa and Schielzeth, 2011) in R (https://www.r-project.org/). Spongy area parameters were tested for repeatability between the two detail images of each cross-section. Repeatability was significant for all parameters, so the data were averaged and the repeatability testing was repeated at the level of cross-sections (two per individual). All of the measured structural variables were tested at this level and all repeatabilities were again significant, indicating that measurements were sufficiently accurate and that the two cross-sections were representative. All details on measurement repeatabilities are presented in Table S1. We then used overall averages of each parameter at the individual level to characterize feather nanostructure.

Reflectance

Reflectance measurements were made using at least eight collected feathers from each individual, using methods detailed and validated in Hegyi et al. (2015). Briefly, the feathers were put on top of one another in a parallel position on a piece of black velvet. Spectral data were collected using OOIBase software, a USB2000 spectrometer, a bifurcated fibre optic probe fitted with a black plastic tube to standardize measurement distance, and a DH2000 deuterium and tungsten halogen light source (Ocean Optics Europe, Ostfildern, Germany). Percent reflectance was given relative to a WS-1 spectralon white reflectance standard (Ocean Optics Europe). Three reflectance spectra were taken from each sample. We calculated three spectral variables from the obtained spectra: mean reflectance (average reflectance from 320 to 700 nm), UV chroma (UV reflectance from 320 to 400 nm, divided by total reflectance) and dominant wavelength (wavelength of maximum reflectance). All of these variables were repeatable (tested in rptR; rI=0.787–0.955, P<0.001), and were therefore averaged at the sample level. We illustrate variation in the three spectral variables in Fig. 2.

Fig. 2.

Averaged reflectance curves of the lower and upper halves of the feather sample after sorting the data by a given spectral variable. (A) Mean reflectance, (B) UV chroma and (C) dominant wavelength. Thin and thick lines refer to low and high values, respectively.

Fig. 2.

Averaged reflectance curves of the lower and upper halves of the feather sample after sorting the data by a given spectral variable. (A) Mean reflectance, (B) UV chroma and (C) dominant wavelength. Thin and thick lines refer to low and high values, respectively.

Five juvenile individuals (one male and four females) had dominant wavelength values greater than 400 nm. All of these samples had been collected late in the season. The reflectance spectra of these five samples showed a characteristic shape difference from the rest of the samples (Fig. S1), suggesting that an environmental impact on the crown feathers had caused a UV cutoff in these birds. We present the analyses without these individuals in the paper, and results including their data are given in Tables S2 and S3.

Statistical analyses

All statistical analyses were conducted in Statistica 5.5 (StatSoft Inc., Dell, Rock Round, TX, USA). To improve normality, we log transformed dominant wavelength, cortex thickness, spongy area thickness and mean keratin rod diameter before analysis. We first tested age and sex effects (and their interaction) on each reflectance and structural trait in general linear models with backward stepwise selection and reintroduction. In the second step, we correlated the nanostructural and macrostructural descriptors with one another and with the reflectance traits. Where we found significant age or sex effects on both structure and colour, we repeated the colour–structure correlations using the residuals of the affected variables.

Effects of age and sex

All three spectral traits differed between the sexes (see statistical details in Table 1). We found higher mean reflectance and UV chroma and lower dominant wavelength in males than in females. UV chroma and dominant wavelength (P=0.068) also varied with age, with higher UV chroma and lower dominant wavelength in adults than in juveniles.

Table 1.

Age- and sex-related variation in aspects of reflectance and feather structure in blue tits: general linear models with backward stepwise selection

Age- and sex-related variation in aspects of reflectance and feather structure in blue tits: general linear models with backward stepwise selection
Age- and sex-related variation in aspects of reflectance and feather structure in blue tits: general linear models with backward stepwise selection

Among the nanostructural traits, spongy area thickness and mean air space diameter (Fig. 3A) were greater in adults than in juveniles and there was a similar tendency (P=0.061) for distance between scatterers. No other nanostructural trait varied with age, and no nanostructural trait differed between sexes. By contrast, all three macrostructural traits differed between sexes, with males having a greater number of barbs (Fig. 3B) and greater overhang length (Fig. 3C) and blue tract length than females; however, no macrostructural trait differed between age classes.

Fig. 3.

Age and sex dependence of some structural parameters of blue tit crown feathers. (A) Age dependence of mean air space diameter, and sex dependence of (B) the number of blue barbs and (C) feather overhang length.

Fig. 3.

Age and sex dependence of some structural parameters of blue tit crown feathers. (A) Age dependence of mean air space diameter, and sex dependence of (B) the number of blue barbs and (C) feather overhang length.

Interrelations among structural variables

There were some strong correlations within groups of variables (details in Table 2). Distance between scatterers was positively correlated with mean rod and air space diameters, and all three macrostructural variables were strongly positively correlated. In contrast, relationships among the groups of structural variables were less prominent, with only a few exceptions. Spongy area thickness was positively correlated with mean air space diameter and distance between scatterers. Spongy area density was strongly negatively correlated with mean air space diameter and distance between scatterers. There was only one significant relationship between nanostructural and macrostructural traits: the number of barbs was negatively correlated with cortex thickness.

Table 2.

Correlations among descriptors of blue tit crown feather structure

Correlations among descriptors of blue tit crown feather structure
Correlations among descriptors of blue tit crown feather structure

Reflectance traits in relation to nanostructure and macrostructure

Among the nanostructural parameters, mean air space diameter (Fig. 4A) and distance between scatterers were negatively correlated with mean reflectance (see detailed results in Table 3). Cortex thickness was negatively correlated with UV chroma (Fig. 4B), and the pattern remained when using the residual reflectance traits. Spongy area thickness was positively correlated with UV chroma, but this pattern vanished when using standardized reflectance traits. Spongy area density and mean keratin rod diameter were not significantly correlated with reflectance parameters. Finally, the s.d. of rod diameter was negatively correlated with dominant wavelength (Fig. 4C). Among the macrostructural parameters, all three measures were positively correlated with mean reflectance and UV chroma, and negatively correlated with dominant wavelength (marginal between the number of barbs and dominant wavelength). When we corrected both macrostructure and reflectance traits for sex, overhang length continued to show a positive correlation with mean reflectance (raw relationship in Fig. 5A), and the number of barbs remained significantly correlated with UV chroma (raw relationship in Fig. 5B). No other correlation remained significant when corrected for sex.

Fig. 4.

Feather nanostructural predictors of blue tit crown reflectance. (A) Mean air space diameter and mean reflectance, (B) barb cortex thickness and UV chroma and (C) the s.d. of keratin rod diameter in the barb spongy layer and dominant wavelength.

Fig. 4.

Feather nanostructural predictors of blue tit crown reflectance. (A) Mean air space diameter and mean reflectance, (B) barb cortex thickness and UV chroma and (C) the s.d. of keratin rod diameter in the barb spongy layer and dominant wavelength.

Table 3.

Correlations (Pearson’s r) between aspects of feather structure and blue tit crown reflectance

Correlations (Pearson’s r) between aspects of feather structure and blue tit crown reflectance
Correlations (Pearson’s r) between aspects of feather structure and blue tit crown reflectance
Fig. 5.

Feather macrostructural predictors of blue tit crown reflectance. (A) Feather overhang length and mean reflectance and (B) the number of blue barbs and UV chroma.

Fig. 5.

Feather macrostructural predictors of blue tit crown reflectance. (A) Feather overhang length and mean reflectance and (B) the number of blue barbs and UV chroma.

Reflectance of the crown of blue tits is a sexually selected ornamental trait (Hunt et al., 1999) that has been found to show sexual dichromatism (Andersson et al., 1998) and an increased expression with age (Delhey and Kempenaers, 2006). Indeed, all reflectance traits we calculated here differed between the sexes, and spectral shape measures also differed between age classes, all in the expected direction (more exaggerated expression in males and in adults). Here, we aimed to find aspects of feather structure that may explain variation in the reflectance properties of the blue tit crown in three contexts: between sexes, between age groups and among individuals. In addition to explaining variation, the investigation of the pathways of reflectance modulation may also help to identify mechanisms that may limit the production of exaggerated character states (e.g. high mean reflectance, high UV chroma), and may thereby determine the information conveyed by crown reflectance as a sexual signal. We measured feather barb nanostructure that probably generates the UV–blue colour of the area. We also measured aspects of feather macrostructure (size and density traits) that may contribute to reflectance differences among individuals.

Strikingly, nanostructural traits showed some age differences, but no significant sex differences, while the reverse was true for macrostructural traits, all of which were strongly sex dependent but age independent. Based on the interrelations of age-corrected data, the s.d. of keratin rod diameter is age dependent and correlated with the dominant wavelength both among and within groups. Rod diameter may therefore explain some of the age difference in dominant wavelength. Following a similar logic, cortex thickness may explain some age difference in UV chroma. Analogously, relationships with macrostructure remaining after sex correction indicate that feather overhang length and number of barbs may explain some of the sexual dichromatism in mean reflectance and UV chroma, respectively. In summary, sexual dichromatism in the mean reflectance and chromaticity of blue tit crown colour seems to originate at least partly from feather macrostructure, and we found no sign that it originated from nanostructure. Sexual dichromatism represents a substantial proportion of variance in reflectance in the blue tit (Andersson et al., 1998; present study), so our findings highlight a striking discrepancy between the mechanisms generating species-specific overall colour (barb nanostructure) and those leading to variation in colour among individuals (feather macrostructure). As macrostructure should be prone to mechanical deterioration, it will also be very important to examine sex differences in crown macrostructure and reflectance immediately after the summer moult, and their changes until and during the breeding season. We need similar analyses as those done previously on seasonal patterns of reflectance (e.g. Delhey et al., 2010).

In contrast to sexual dichromatism, age differences in reflectance (see also Delhey and Kempenaers, 2006) were partly explained by barb nanostructure but were not significantly affected by feather macrostructure. Therefore, it is unlikely that age-related patterns in crown colour are influenced by external feather abrasion. Again, it will be very important to examine the within-individual changes of barb nanostructure and reflectance with subsequent moults (nanostructure is not very likely to substantially change with abrasion but this needs to be confirmed as well). Because of the well-known improvement of individual condition and intrinsic attributes with age (e.g. Forslund and Pärt, 1995), age-dependent nanostructural characters that predict reflectance (here particularly, barb cortex thickness) may be especially likely candidates for indicators of condition or stable aspects of individual quality.

Two previous studies have examined the structural basis of plumage sexual dichromatism. In eastern bluebirds, robust sex differences were found in nanostructural traits, while the only macrostructural measure (number of melanized barbules) was independent of sex (Shawkey et al., 2005). In two species of fairy-wren, sex- and morph-specific differences were suggested in both macrostructure and nanostructure (Enbody et al., 2017). We found no previous studies of age dependence in colour-producing feather nanostructures or macrostructures, so no comparison can be made at this point. It is nevertheless interesting that age and sex differences were produced at different levels of structural organization in our study species.

The most crucial characteristic of a sexually selected signal trait is the information it conveys to the receivers of the signal, which in turn depends on what creates signal variation. Unfortunately, the mechanism that generates within-species variation in ‘structural’ plumage colour signals is essentially unknown, with only one published study on iridescent colour (Doucet et al., 2006) and two on non-iridescent colour (Shawkey et al., 2003, 2005; but these two concerned the same species and population). Here, we examined a non-iridescent ornament and paid attention to both macrostructure and nanostructure of feathers.

Our macrostructural variables described the amount and density of colour-producing material in the feather. The positive relationship between overhang length and mean reflectance may reflect the increase of the total coloured surface in the crown plumage. The positive correlation of UV chroma with the number of blue barbs in the feather is also logical as there should be a positive link between the density of coloured tissue and the spectral purity of the reflected light. Similar relationships with macrostructure have previously been found for melanin-coloured plumage traits (Galván, 2011; D'Alba et al., 2014).

Concerning the nanostructural predictors of within-species variation in non-iridescent chromatic structural colour, most of our results were consistent with those of the only published research project, that of Shawkey et al. (2003, 2005) on the eastern bluebird. Mean reflectance was negatively correlated with two variables that were positively correlated with each other: mean air space diameter and the distance between scatterers. These patterns are straightforward as both variables probably correlate negatively with the reflective surface density of the spongy area (given the same thickness of the spongy layer), and a reduced surface density is expected to generate lower overall reflectance from the spongy layer (see e.g. simulations in Yu et al., 2017; see similar empirical results in Shawkey et al., 2005).

UV chroma was negatively correlated with cortex thickness (see also Shawkey et al., 2005). This is presumably due to the selective filtering of UV wavelengths by the keratin matrix of the cortex (Brink and van der Berg, 2004). Concerning dominant wavelength, its only nanostructural predictor was the s.d. of keratin rod diameter, and birds with more variable rod diameter showed more UV-shifted dominant wavelengths. Strikingly, mean values of keratin rod diameter, air space diameter and inter-scatterer distance did not predict dominant wavelength.

According to the accepted theory on the generation of non-iridescent structural plumage colour, the reflectance curve directly mirrors the distribution of light path lengths (i.e. inter-scatterer distances) in the spongy layer (Prum et al., 1998, 1999). Looking at the reflectance curves of the blue tit crown (see Fig. 2 for characteristic spectra), we see that these must originate from light path distributions of varying shapes. Assuming a direct correspondence with reflectance curve shape, light path distributions probably become more and more right-tailed as dominant wavelength decreases. This increasing tail may have caused the negative correlation of keratin rod s.d. with dominant wavelength in our dataset. There are multiple complementary explanations for the lack of correlation of dominant wavelength with mean particle diameters, pore diameters and estimated path lengths. For example, means may poorly reflect the true peaks of the non-normal distribution. The assumed short-range but not long-range order of the spongy layer (Noh et al., 2010) may also make it difficult to reliably derive the path length distribution of light from the observed distributions of particle and pore diameters in 2D planes (Shawkey et al., 2009). Even more sophisticated approaches than ours (2D Fourier analysis of TEM images) to directly estimate the path length distribution have previously failed to accurately estimate the dominant periodicity of the 3D spongy layer, with the result that calculated mean dominant wavelength was (roughly) accurate at the species level (e.g. Shawkey et al., 2006b), but individual differences could not be reliably discerned (Shawkey et al., 2003; see also Shawkey et al., 2009). To our knowledge, the only known method to reliably estimate the distributions of path lengths is the high-energy, high-resolution small-angle X-ray scattering (SAXS) screening of single, intact feather barbs (Saranathan et al., 2012). This method may give sufficient precision to capture individual differences in dominant wavelength, but it currently requires a synchrotron, limiting its application for evolutionary ecologists. Further technical developments are clearly needed in this area.

The information content of structural plumage reflectance, including its condition dependence, is unclear. If energetic constraints dominate, condition dependence may be most likely in reflectance traits that depend on the quantity of keratin invested (Meadows et al., 2012). It could be assumed that mean reflectance is a direct function of the quantity and density of scattering or absorbing elements at the macrostructural and nanostructural level (Prum, 2006), whereas the chromaticity of structural colour should instead more strongly depend on the regularity of nanostructure (Prum et al., 1999). Finally, the wavelength of maximum reflectance (dominant wavelength) is assumed to originate from the geometry of the spongy layer, i.e. the distances between scatterers that reinforce certain wavelengths by constructive interference (Prum et al., 1999). With these reflectance modulation pathways in mind, we would expect the greatest condition dependence for mean reflectance and the lowest for dominant wavelength, whereas regularity-dependent measures of spectral purity should be sensitive to stress during moulting.

However, the few existing studies of within-species variation show more complex interrelationships between nanostructure and reflectance. Mean reflectance indeed correlates negatively with quantity and/or density traits such as cortex thickness, keratin rod size and distance between scatterers, whereas UV chroma correlates with aspects of the regularity of the spongy layer (Shawkey et al., 2005; this study). However, UV chroma is also affected by quantity aspects of the nanostructure (cortex thickness: Shawkey et al., 2005; this study; number of circular keratin rods: Shawkey et al., 2003; possibly also spongy layer density: Yin et al., 2012). Finally, the main predictor of dominant wavelength in our study is apparently not one dominant structural parameter of the spongy layer (inter-scatterer distance; Shawkey et al., 2005) but rather its composition (s.d. of keratin rod diameter). Dominant wavelength has also been shown to correlate with quantity/density aspects of the nanostructure (spongy layer thickness: Shawkey et al., 2005; cortex thickness: Doucet et al., 2006). Concerning macrostructural traits, most previous quantitative studies of such traits examined melanin-based colour and are therefore not relevant (Galván, 2011; D'Alba et al., 2014). The only macrostructural trait measured by Shawkey et al. (2005), the number of melanin-containing barbules, reduced both mean reflectance and spectral purity of structural colour. Our present data also indicate that quantitative aspects of feather macrostructure may influence both mean reflectance (overhang length) and spectral shape (number of barbs).

It therefore currently seems that all descriptors of structural colour may depend on keratin quantity and density, and thus all of them may show condition dependence during moulting via this pathway. Experimental studies on condition dependence are few, and the results obtained were mixed. It has been suggested that structural plumage colour is dependent on developmental stress rather than body condition at moulting (Peters et al., 2011). Correlations with stress hormone levels and significant effects of various experimental manipulations (see examples in Table 4) indicate that structural colour may depend on both stress and true body condition. Mean reflectance was somewhat more frequently influenced by manipulations than dominant wavelength or chromaticity (Table 4). This suggests that condition or stress during moulting may indeed be more likely to affect the quantity than the regularity aspects of feather structure.

Table 4.

Some studies that provided correlative or experimental support for the phenotypic plasticity of structural plumage colour

Some studies that provided correlative or experimental support for the phenotypic plasticity of structural plumage colour
Some studies that provided correlative or experimental support for the phenotypic plasticity of structural plumage colour

In summary, our results surprisingly indicate that feather macrostructure may explain a similar amount of phenotypic variation in ‘structural’ plumage colour traits as nanostructure. Moreover, in our study population, macrostructure seemed to be an important determinant of sexual dichromatism in mean reflectance and spectral purity. Studies of feather macrostructure and structural colour are effectively lacking, although macrostructure is technically much more accessible than nanostructure. More research is urgently needed on this topic as the expected information content, vulnerability and seasonal expression profile of structural colour signals strongly depend on whether their main origin is macrostructural or nanostructural variation. In the blue tit, further observational and experimental studies are needed to elucidate the phenotypic plasticity and quality indicator value of the structural determinants of crown colour we identified here, and the consequences for signal information content.

We thank Mónika Duleba, Rita Főző, Rita Hargitai, Márton Herényi, Dorottya Kiss, Gábor Markó, Gergely Nagy and Anikó Szegedi for help with the fieldwork; Zsófia Kováts and Mónika Truszka for assistance in the laboratory; and Károly Bóka for preparing several SEM photographs, including the one shown in Fig. 1B.

Author contributions

Conceptualization: G.H., M.L., D.K., T.C., P.L., B.R., E.S., J.T.; Investigation: G.H., M.L., D.K., T.C., P.L., B.R., E.S., J.T.; Writing - original draft: G.H., M.L., D.K., T.C., P.L., B.R., E.S., J.T.; Writing - review & editing: G.H.; Supervision: G.H.; Funding acquisition: G.H., B.R.

Funding

This work was supported by the Nemzeti Kutatás-Fejlesztési és Innovációs Hivatal (NKFIH) (grant numbers F68295, K124443), the Pilis Park Forestry, and the Erdők a Közjóért Alapítvány. The project was also supported by the János Bolyai Research Scholarship of the Magyar Tudományos Akadémia (Hungarian Academy of Sciences) to B.R. and E.S. (grant numbers BO_050_17 and BO_663_17).

Andersson
,
S.
,
Örnborg
,
J.
and
Andersson
,
M.
(
1998
).
Ultraviolet sexual dimorphism and assortative mating in blue tits
.
Proc. R. Soc. B
265
,
445
-
450
.
Brink
,
D. J.
and
van der Berg
,
N. G.
(
2004
).
Structural colours from the feathers of the bird Bostrychia hagedash
.
J. Phys. D
37
,
813
-
818
.
Broggi
,
J.
,
Gamero
,
A.
,
Hohtola
,
E.
,
Orell
,
M.
and
Nilsson
,
J.-Å.
(
2011
).
Interpopulation variation in contour feather structure is environmentally determined in great tits
.
PLoS ONE
6
,
e24942
.
Caro
,
T.
,
Walker
,
H.
,
Santana
,
S. E.
and
Stankowich
,
T.
(
2017
).
The evolution of anterior coloration in carnivorans
.
Behav. Ecol. Sociobiol.
71
,
177
.
D'Alba
,
L.
,
Van Hemert
,
C.
,
Spencer
,
K. A.
,
Heidinger
,
B. J.
,
Gill
,
L.
,
Evans
,
N. P.
,
Monaghan
,
P.
,
Handel
,
C. M.
and
Shawkey
,
M. D.
(
2014
).
Melanin-based color of plumage: role of condition and of feathers’ microstructure
.
Integr. Comp. Biol.
54
,
633
-
644
.
Delhey
,
K.
and
Kempenaers
,
B.
(
2006
).
Age differences in blue tit Parus caeruleus plumage colour: within-individual changes or colour-biased survival?
J. Avian Biol.
37
,
339
-
348
.
Delhey
,
K.
,
Peters
,
A.
,
Johnsen
,
A.
and
Kempenaers
,
B.
(
2007
).
Brood sex ratio and male UV ornamentation in blue tits (Cyanistes caeruleus): correlational evidence and an experimental test
.
Behav. Ecol. Sociobiol.
61
,
853
-
862
.
Delhey
,
K.
,
Burger
,
C.
,
Fiedler
,
W.
and
Peters
,
A.
(
2010
).
Seasonal changes in colour: a comparison of structural, melanin- and carotenoid-based plumage colours
.
PLoS ONE
5
,
e11582
.
Doucet
,
S. M.
,
Shawkey
,
M. D.
,
Hill
,
G. E.
and
Montgomerie
,
R.
(
2006
).
Iridescent plumage in satin bowerbirds: structure, mechanisms and nanostructural predictors of individual variation in colour
.
J. Exp. Biol.
209
,
380
-
390
.
Doutrelant
,
C.
,
Grégoire
,
A.
,
Midamegbe
,
A.
,
Lambrechts
,
M.
and
Perret
,
P.
(
2012
).
Female plumage coloration is sensitive to the cost of reproduction. An experiment in blue tits
.
J. Anim. Ecol.
81
,
87
-
96
.
Doyle
,
A.
and
Siefferman
,
L.
(
2014
).
Supplemental food increases ornamentation of male nestling eastern bluebirds
.
J. Field Ornithol.
85
,
31
-
39
.
Eliason
,
C. M.
and
Shawkey
,
M. D.
(
2012
).
A photonic heterostructure produces diverse iridescent colours in duck wing patches
.
J. R. Soc. Interface
9
,
2279
-
2289
.
Enbody
,
E. D.
,
Lantz
,
S. M.
and
Karubian
,
J.
(
2017
).
Production of plumage ornaments among males and females of two closely related tropical passerine bird species
.
Ecol. Evol.
7
,
4024
-
4034
.
Fitzpatrick
,
S.
(
1998
).
Colour schemes for birds: structural coloration and signals of quality in feathers
.
Ann. Zool. Fenn.
35
,
67
-
77
.
Forslund
,
P.
and
Pärt
,
T.
(
1995
).
Age and reproduction in birds: hypotheses and tests
.
Trends Ecol. Evol.
10
,
374
-
378
.
Galván
,
I.
(
2011
).
Ultraviolet-blue plumage colouration can be perceived as an indicator of fluctuating asymmetry by blue tits (Cyanistes caeruleus)
.
J. Ornithol.
152
,
223
-
230
.
Griffith
,
S. C.
,
Parker
,
T. H.
and
Olson
,
V. A.
(
2006
).
Melanin-versus carotenoid-based sexual signals: is the difference really so black and red?
Anim. Behav.
71
,
749
-
763
.
Griggio
,
M.
,
Serra
,
L.
,
Licheri
,
D.
,
Campomori
,
C.
and
Pilastro
,
A.
(
2009
).
Moult speed affects structural feather ornaments in the blue tit
.
J. Evol. Biol.
22
,
782
-
792
.
Guindre-Parker
,
S.
and
Love
,
O. P.
(
2014
).
Revisiting the condition-dependence of melanin-based plumage
.
J. Avian Biol.
45
,
29
-
33
.
Hedenström
,
A.
(
1995
).
Song flight performance in the skylark Alauda arvensis
.
J. Avian Biol.
26
,
337
-
342
.
Hegyi
,
G.
,
Laczi
,
M.
,
Nagy
,
G.
,
Szász
,
E.
,
Kötél
,
D.
and
Török
,
J.
(
2015
).
Stable correlation structure among multiple plumage colour traits: can they work as a single signal?
Biol. J. Linn. Soc.
114
,
92
-
108
.
Henderson
,
L. J.
,
Heidinger
,
B. J.
,
Evans
,
N. P.
and
Arnold
,
K. E.
(
2013
).
Ultraviolet crown coloration in female blue tits predicts reproductive success and baseline corticosterone
.
Behav. Ecol.
24
,
1299
-
1305
.
Hoskin
,
C. J.
,
James
,
S.
and
Grigg
,
G. C.
(
2009
).
Ecology and taxonomy-driven deviations in the frog call-body size relationship across the diverse Australian frog fauna
.
J. Zool.
278
,
36
-
41
.
Hunt
,
S.
,
Cuthill
,
I. C.
,
Bennett
,
A. T. D.
and
Griffiths
,
R.
(
1999
).
Preferences for ultraviolet partners in the blue tit
.
Anim. Behav.
58
,
809
-
815
.
Jacot
,
A.
and
Kempenaers
,
B.
(
2007
).
Effects of nestling condition on UV plumage traits in blue tits: an experimental approach
.
Behav. Ecol.
18
,
34
-
40
.
Kinoshita
,
S.
and
Yoshioka
,
S.
(
2005
).
Structural colors in nature: The role of regularity and irregularity in the structure
.
Chemphyschem
6
,
1142
-
1459
.
Lee
,
C.-C.
,
Liao
,
S.-F.
and
Vukusic
,
P.
(
2015
).
Measuring and modelling the reflectance spectra of male Swinhoe's pheasant feather barbules
.
J. R. Soc. Interface
12
,
20141354
.
Limbourg
,
T.
,
Mateman
,
A. C.
,
Andersson
,
S.
and
Lessells
,
C. M.
(
2004
).
Female blue tits adjust parental effort to manipulated male UV attractiveness
.
Proc. R. Soc B.
271
,
1903
-
1908
.
Limbourg
,
T.
,
Mateman
,
A. C.
and
Lessells
,
C. M.
(
2013
).
Opposite differential allocation by males and females of the same species
.
Biol. Lett.
9
,
20120835
.
MacDonald
,
I. F.
,
Kempster
,
B.
,
Zanette
,
L.
and
MacDougall-Shackleton
,
S. A.
(
2006
).
Early nutritional stress impairs development of a song-control brain region in both male and female juvenile song sparrows (Melospiza melodia) at the onset of song learning
.
Proc. R. Soc. B.
273
,
2559
-
2564
.
Mahr
,
K.
,
Griggio
,
M.
,
Granatiero
,
M.
and
Hoi
,
H.
(
2012
).
Female attractiveness affects paternal investment: experimental evidence for male differential allocation in blue tits
.
Front. Zool.
9
,
14
.
Maia
,
R.
,
Caetano
,
J. V. O.
,
Bao
,
S. N.
and
Macedo
,
R. H.
(
2009
).
Iridescent structural colour production in male blue-black grassquit feather barbules: the role of keratin and melanin
.
J. R. Soc. Interface
6
,
S203
-
S211
.
Maia
,
R.
,
Macedo
,
R. H. F.
and
Shawkey
,
M. D.
(
2012a
).
Nanostructural self-assembly of iridescent feather barbules through depletion attraction of melanosomes during keratinisation
.
J. R. Soc. Interface
9
,
734
-
743
.
Maia
,
R.
,
Brasileiro
,
L.
,
Lacava
,
R. V.
and
Macedo
,
R. H.
(
2012b
).
Social environment affects acquisition and color of structural nuptial plumage in a sexually dimorphic tropical passerine
.
PLoS ONE
7
,
e47501
.
McGlothlin
,
J. W.
,
Duffy
,
D. L.
,
Henry-Freeman
,
J. L.
and
Ketterson
,
E. D.
(
2007
).
Diet quality affects an attractive white plumage pattern in dark-eyed juncos (Junco hyemalis)
.
Behav. Ecol. Sociobiol.
61
,
1391
-
1399
.
Meadows
,
M. G.
,
Roudybush
,
T. E.
and
McGraw
,
K. J.
(
2012
).
Dietary protein level affects iridescent coloration in Anna's hummingbirds, Calypte anna
.
J. Exp. Biol.
215
,
2742
-
2750
.
Midamegbe
,
A.
,
Gregoire
,
A.
,
Perret
,
P.
and
Doutrelant
,
C.
(
2011
).
Female-female aggressiveness is influenced by female coloration in blue tits
.
Anim. Behav.
82
,
245
-
253
.
Nakagawa
,
S.
and
Schielzeth
,
H.
(
2011
).
Repeatability for Gaussian and non-Gaussian data: a practical guide for biologists
.
Biol. Rev.
85
,
935
-
956
.
Noh
,
H.
,
Liew
,
S. F.
,
Saranathan
,
V.
,
Mochrie
,
S. G. J.
,
Prum
,
R. O.
,
Dufresne
,
E. R.
and
Cao
,
H.
(
2010
).
How noniridescent colors are generated by quasi-ordered structures of bird feathers
.
Adv. Mater.
22
,
2871
-
2880
.
Olson
,
V. A.
and
Owens
,
I. P. F.
(
1998
).
Costly sexual signals: are carotenoids rare, risky or required?
Trends Ecol. Evol.
13
,
510
-
514
.
Osorio
,
D.
and
Ham
,
A. D.
(
2002
).
Spectral reflectance and directional properties of structural coloration in bird plumage
.
J. Exp. Biol.
205
,
2017
-
2027
.
Peters
,
A.
,
Kurvers
,
R. H. J. M.
,
Roberts
,
M. L.
and
Delhey
,
K.
(
2011
).
No evidence for general condition-dependence of structural plumage colour in blue tits: an experiment
.
J. Evol. Biol.
24
,
976
-
987
.
Poesel
,
A.
,
Dabelsteen
,
T.
,
Darden
,
S.-K.
,
Delhey
,
K.
and
Peters
,
A.
(
2007
).
Territorial responses of male blue tits, Cyanistes caeruleus, to UV-manipulated neighbours
.
J. Ornithol.
148
,
179
-
187
.
Prum
,
R. O.
(
2006
).
Anatomy, physics and evolution of avian structural colors
. In
Bird Coloration
, Vol.
1
(ed.
G. E.
Hill
and
K. J.
McGraw
), pp.
295
-
353
.
Cambridge, MA
:
Harvard University Press
.
Prum
,
R. O.
,
Torres
,
R. H.
,
Williamson
,
S.
and
Dyck
,
J.
(
1998
).
Coherent light scattering by blue feather barbs
.
Nature
396
,
28
-
29
.
Prum
,
R. O.
,
Torres
,
R. H.
,
Williamson
,
S.
and
Dyck
,
J.
(
1999
).
Two-dimensional Fourier analysis of the spongy medullary keratin of structurally coloured feather barbs
.
Proc. R. Soc. B
266
,
13
-
22
.
Prum
,
R. O.
,
Dufresne
,
E. R.
,
Quinn
,
T.
and
Waters
,
K.
(
2009
).
Development of colour-producing β-keratin nanostructures in avian feather barbs
.
J. R. Soc. Interface
6
,
S253
-
S265
.
Rémy
,
A.
,
Grégoire
,
A.
,
Perret
,
P.
and
Doutrelant
,
C.
(
2010
).
Mediating male-male interactions: the role of the UV-blue crest coloration in blue tits
.
Behav. Ecol. Sociobiol.
64
,
1839
-
1847
.
Saranathan
,
V.
,
Forster
,
J. D.
,
Noh
,
H.
,
Liew
,
S.-F.
,
Mochrie
,
S. G. J.
,
Cao
,
H.
,
Dufresne
,
E. R.
and
Prum
,
R. O.
(
2012
).
Structure and optical function of amorphous photonic nanostructures from avian feather barbs: a comparative small-angle X-ray scattering (SAXS) analysis of 230 bird species
.
J. R. Soc. Interface
9
,
2563
-
2580
.
Shawkey
,
M. D.
and
D'Alba
,
L.
(
2017
).
Interactions between colour-producing mechanisms and their effects on the integumentary colour palette
.
Phil. Trans. R. Soc. B
372
,
20160536
.
Shawkey
,
M. D.
,
Estes
,
A. M.
,
Siefferman
,
L. M.
and
Hill
,
G. E.
(
2003
).
Nanostructure predicts intraspecific variation in ultraviolet-blue plumage colour
.
Proc. R. Soc. B
270
,
1455
-
1460
.
Shawkey
,
M. D.
,
Estes
,
A. M.
,
Siefferman
,
L.
and
Hill
,
G. E.
(
2005
).
The anatomical basis of sexual dichromatism in non-iridescent ultraviolet-blue structural coloration of feathers
.
Biol. J. Linn. Soc.
84
,
259
-
271
.
Shawkey
,
M. D.
,
Hill
,
G. E.
,
McGraw
,
K. J.
,
Hood
,
W. R.
and
Huggins
,
K.
(
2006a
).
An experimental test of the contributions and condition-dependence of microstructure and carotenoids in yellow plumage coloration
.
Proc. R. Soc. B
273
,
2985
-
2991
.
Shawkey
,
M. D.
,
Balenger
,
S. L.
,
Hill
,
G. E.
,
Johnson
,
L. S.
,
Keyser
,
A. J.
and
Siefferman
,
L.
(
2006b
).
Mechanisms of evolutionary change in structural plumage coloration among bluebirds (Sialia spp.)
.
J. R. Soc. Interface
3
,
527
-
532
.
Shawkey
,
M. D.
,
Saranathan
,
V.
,
Pálsdóttir
,
H.
,
Crum
,
J.
,
Ellisman
,
M. H.
,
Auer
,
M.
and
Prum
,
R. O.
(
2009
).
Electron tomography, three-dimensional Fourier analysis and colour prediction of a three-dimensional amorphous biophotonic nanostructure
.
J. R. Soc. Interface
6
,
S213
-
S220
.
Sheldon
,
B. C.
,
Andersson
,
S.
,
Griffith
,
S. C.
,
Örnborg
,
J.
and
Sendecka
,
J.
(
1999
).
Ultraviolet colour variation influences blue tit sex ratios
.
Nature
402
,
874
-
877
.
Siefferman
,
L.
and
Hill
,
G. E.
(
2005
).
Evidence for sexual selection on structural plumage coloration in female eastern bluebirds (Sialia sialis)
.
Evolution
59
,
1819
-
1828
.
Siefferman
,
L.
and
Hill
,
G. E.
(
2007
).
The effect of rearing environment on blue structural coloration of eastern bluebirds (Sialia sialis)
.
Behav. Ecol. Sociobiol.
61
,
1839
-
1846
.
Siefferman
,
L.
,
Liu
,
M.
,
Navara
,
K. J.
,
Mendonça
,
M. T.
and
Hill
,
G. E.
(
2013
).
Effect of prenatal and natal administration of testosterone on production of structurally based plumage coloration
.
Physiol. Biochem. Zool.
86
,
323
-
332
.
Slagsvold
,
T.
,
Dale
,
S.
and
Kruszewicz
,
A.
(
1995
).
Predation favors cryptic coloration of breeding male pied flycatchers
.
Anim. Behav.
50
,
1109
-
1121
.
Stavenga
,
D. G.
,
Tinbergen
,
J.
,
Leertouwer
,
H. L.
and
Wilts
,
B. D.
(
2011
).
Kingfisher feathers – colouration by pigments, spongy nanostructures and thin films
.
J. Exp. Biol.
214
,
3960
-
3967
.
Stavenga
,
D. G.
,
van der Kooi
,
C. J.
and
Wilts
,
B. D.
(
2017
).
Structural coloured feathers of mallards act by simple multilayer photonics
.
J. R. Soc. Interface
14
,
20170407
.
Thomas
,
A. L. R.
(
1993
).
The aerodynamics of birds’ tails
.
Phil. Trans. R. Soc. B
340
,
361
-
380
.
Tinbergen
,
J.
,
Wilts
,
B. D.
and
Stavenga
,
D. G.
(
2013
).
Spectral tuning of Amazon parrot feather coloration by psittacofulvin pigments and spongy structures
.
J. Exp. Biol.
216
,
4358
-
4364
.
Vágási
,
C. I.
,
Pap
,
P. L.
,
Vincze
,
O.
,
Benkő
,
Z.
,
Marton
,
A.
and
Barta
,
Z.
(
2012
).
Haste makes waste but condition matters: molt rate – feather quality trade-off in a sedentary songbird
.
PLoS ONE
7
,
e40651
.
Yin
,
H. W.
,
Shi
,
L.
,
Sha
,
J.
,
Li
,
Y. Z.
,
Qin
,
Y. H.
,
Dong
,
B. Q.
,
Meyer
,
S.
,
Liu
,
X. H.
,
Zhao
,
L.
and
Zi
,
J.
(
2006
).
Iridescence in the neck feathers of domestic pigeons
.
Phys. Rev. E
74
,
051916
.
Yin
,
H. W.
,
Dong
,
B. Q.
,
Liu
,
X. H.
,
Zhan
,
T. R.
,
Shi
,
L.
,
Zi
,
J.
and
Yablonovitch
,
E.
(
2012
).
Amorphous diamond-structured photonic crystal in the feather barbs of the scarlet macaw
.
Proc. Natl. Acad. Sci. USA
109
,
10798
-
10801
.
Yu
,
S. C.
,
Wang
,
C.
,
Zhang
,
Y. C.
,
Dong
,
B. Q.
,
Jiang
,
Z.
,
Chen
,
X. F.
,
Chen
,
W.
and
Sun
,
C.
(
2017
).
Design of non-deterministic quasi-random nanophotonic structures using Fourier space representations
.
Sci. Rep.
7
,
3752
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information