Growing in Cities: An Urban Penalty for Wild Birds? A Study of Phenotypic Differences between Urban and Rural Great Tit Chicks (Parus major)

Biard, Clotilde; Brischoux, François; Meillère, Alizée; Michaud, Bruno; Nivière, Manon; Ruault, Stéphanie; Vaugoyeau, Marie; Angelier, Frédéric

doi:10.3389/fevo.2017.00079

ORIGINAL RESEARCH article

Front. Ecol. Evol., 18 July 2017

Sec. Behavioral and Evolutionary Ecology

Volume 5 - 2017 | https://doi.org/10.3389/fevo.2017.00079

This article is part of the Research Topic Behavioural and Ecological Consequences of Urban Life in Birds View all 31 articles

Growing in Cities: An Urban Penalty for Wild Birds? A Study of Phenotypic Differences between Urban and Rural Great Tit Chicks (Parus major)

$\r\nClotilde Biard*$ Clotilde Biard¹^*

Bruno Michaud²

¹Institut d'Écologie et des Sciences de l'Environnement de Paris, Sorbonne universités, UPMC Université Paris 06, UPEC, Paris 7, Centre National de la Recherche Scientifique, Institut National de la Recherche Agronomique, IRD, Paris, France
²Centre d'Etudes Biologiques de Chizé, Centre National de la Recherche Scientifique-ULR, UMR 7372, Villiers en Bois, France

Urban sprawl is associated with deep and intense modifications of the natural habitats of wild vertebrates. Although, many species are unable to cope with such an environment, a few species can be found in cities and can help us assessing the impact of urbanization on wildlife. Urban-related environmental modifications are multiple and some of them seem beneficial while others seem rather detrimental to wild vertebrates. Moreover, the impact of these modifications on wild vertebrates is likely to vary depending on the phase of the annual life-cycle. Therefore, it is challenging to get a comprehensive picture of the impact of urbanization on wild vertebrates. Overall, urbanization is usually associated with reduced breeding performances in wild birds, but the impact of urbanization on the phenotype and quality of developing offspring has been less studied. In this study, we specifically investigated the impact of urbanization on several proxies of individual quality in great tits (Parus major). We concomitantly measured body size (tarsus length and body mass), plumage coloration, and telomere length in 14-days old chicks issued from 4 populations (two pairs of urban/rural populations located in two different geographical areas of France). First, rural chicks were significantly taller and heavier than urban birds although this impact of urbanization on body size/body mass appears only true for the most urbanized site. Interestingly, body size was also affected by the geographical area of capture, suggesting that regional environmental conditions may attenuate or exacerbate the influence of urbanization on nestling growth. Second, the carotenoid-based yellow plumage of rural nestlings was more colorful than that of urban birds, independently of the area of capture. This suggests that urban birds probably have a low-carotenoid diet relative to rural birds. Finally, telomere length did not differ between rural and urban chicks. These results suggest that urbanization probably imposes large developmental constraints in wild vertebrates and that this impact may primarily be related to constraining nutritional conditions.

Introduction

Urban sprawl is a worldwide phenomenon that deeply affects environmental conditions, and consequently, avian biodiversity (Grimm et al., 2008; Shanahan et al., 2014). Thus, cities are characterized by multiple abiotic and biotic modifications of the natural environment such as fragmentation (Crooks et al., 2004), human disturbance (Sauvajot et al., 1998; Fernández-Juricic, 2002), and increased pollutions: noise (Slabbekoorn and Peet, 2003; Meillère et al., 2015b), light (Dominoni et al., 2013; Spoelstra and Visser, 2014), electromagnetic (Engels et al., 2014), and chemical pollution (Roux and Marra, 2007; Isaksson, 2010; Bichet et al., 2013). Urbanization is also associated with altered predation risk (Fischer et al., 2012; Møller and Ibáñez-Álamo, 2012), inter and intra-specific competition (Duckworth, 2014), disease prevalence (Bichet et al., 2013; Giraudeau et al., 2014), food type and abundance (Robb et al., 2008; Harrison et al., 2010; Galbraith et al., 2015). All these modifications raise new important selective constraints for wild vertebrates, and as a consequence, biodiversity is overall impoverished in cities (Clergeau et al., 2006; McKinney, 2008). While some species seem unable to cope with this modified environment (urban avoiders) others can still be found in this specific environment (urban-tolerant species, Chace and Walsh, 2006). Comparative studies have allowed ecologists to determine the specific characteristics that may be compatible with an urban way of life in wild vertebrates (Møller, 2008, 2009; Evans et al., 2011; Sol et al., 2013; Jokimäki et al., 2016). In addition, an individual-based approach is also necessary to fully understand the ecological processes that may constrain the ability of wild populations to cope with urbanization (Gil and Brumm, 2014). Several studies have for instance demonstrated that living in an urban environment is associated with phenotypical changes in vertebrates (e.g., body size: Bókony et al., 2012; Meillère et al., 2017, secondary sexual signals: Fuller et al., 2007; Giraudeau et al., 2015, stress physiology: Bonier, 2012, immunity: Audet et al., 2015). Moreover, a large number of studies have also shown that urbanization affects reproductive performances, such as phenology or the number of young produced (Both et al., 2004; Chamberlain et al., 2009; Ausprey and Rodewald, 2011; Stracey and Robinson, 2012; Deviche and Davies, 2014; Møller et al., 2015).

Surprisingly, the influence of urbanization on the phenotype and quality of developing offspring has been much less studied (Heiss et al., 2009; Seress et al., 2012; Salmón et al., 2016). Early developmental conditions are however crucial to study because they are known to have important long-term consequences on individual phenotypes and adult performances (survival and reproduction, Lindström, 1999; Metcalfe and Monaghan, 2001; Monaghan, 2008). In terms of nutritional conditions, urbanization is very likely to be constraining for most bird species. Several studies have shown that the urban environment is associated with modifications of the arthropod community (Helden et al., 2012), and thus, with a lack of protein-rich diet (McIntyre et al., 2001; Shochat et al., 2004), vitamins, carotenoids (Isaksson and Andersson, 2007; Møller et al., 2010; Giraudeau and McGraw, 2014; Giraudeau et al., 2015), and calcium (Heiss et al., 2009). Similarly, urban-related disturbance is likely to affect the quality of parental care (e.g., incubation commitment or brood provisioning, Ibáñez-Álamo and Soler, 2010; Seress et al., 2011; Meillère et al., 2015a), and therefore the nutritional conditions of offspring. In addition, other urban-related environmental conditions could also constrain the development of offspring through direct or parental-mediated effects (e.g., chemical pollution, light pollution, intra- and inter-specific competition, predation risk, Shanahan et al., 2014).

In this study, we aimed to better understand the consequences of living in an urban environment for developing wild vertebrates. To do so, we propose to study concomitantly several complementary proxies of developmental conditions (growth, carotenoid-based coloration, and telomere length). First, growth is classically used to monitor nutritional developmental conditions in wild vertebrates and poor nutritional conditions are associated with delayed growth, an overall poor health status and a lower probability to fledge and to recruit into the population (Gebhardt-Henrich and Richner, 1998; Naef-Daenzer and Keller, 1999; Naef-Daenzer et al., 2001). Second, carotenoids are also relevant to study because plasma carotenoid levels are known to be positively associated with immunity (Saino et al., 2003b; Biard et al., 2005; Cucco et al., 2006), and growth (Tschirren et al., 2003; Biard et al., 2005, 2006; Berthouly et al., 2008). In addition to these physiological functions, carotenoids are also allocated into carotenoid-based signals (e.g., plumage) in developing animals (Fitze et al., 2003; Isaksson et al., 2006). Such signaling is crucial because it is used by parents to assess the viability/quality of the chicks and therefore to determine their parental investment (Saino et al., 2003a; Pike et al., 2007; Thorogood et al., 2011). Therefore, carotenoids mediate a trade-off between protective physiological functions (e.g., immunity, growth, and antioxidants) and signaling in developing animals. Because carotenoids cannot be synthetized by vertebrates and are only acquired through the diet (Isaksson, 2009), carotenoid-based plumage appear relevant to better assess the environmental constraints associated with an urban way of life in developing vertebrates (Giraudeau et al., 2015). Third, offspring survival is also classically used because it is a direct measure of reproductive success. However, it is often very challenging to follow the young after fledging in wild populations, and therefore, long-term survival is difficult to assess. In that respect, telomere length appears as a promising molecular marker to assess the survival probability of offspring after fledging (Monaghan, 2010, 2014). Telomeres are long repetitive sequences of non-coding DNA that are located at the end of chromosomes. They protect chromosomes from degradation and they also play a crucial role during the replication process (Blackburn, 2000). Telomeres usually shorten as individual age, and importantly, environmental stressors and the occurrence of oxidative stress have been shown to accelerate this shortening (von Zglinicki, 2002). Moreover, telomere length has also been linked with survival (Barrett et al., 2013; Fairlie et al., 2016) and it can be used as a reliable proxy of an individual's longevity (Heidinger et al., 2012).

Several studies have examined the influence of urbanization on these variables. For instance, urbanization often delays and impairs growth in bird species (Chamberlain et al., 2009; Heiss et al., 2009; Seress et al., 2012; Bailly et al., 2016). Similarly, urbanization is often associated with lower offspring productivity in wild vertebrates (Chamberlain et al., 2009; Bailly et al., 2016). However, these results appear somewhat inconsistent and sometimes differ between species or populations (reviewed in Chamberlain et al., 2009). Overall, carotenoid-based plumage color is paler in birds living in cities than in those living in rural areas (Eeva et al., 1998; Hõrak et al., 2000, 2001; Isaksson et al., 2005; Giraudeau and McGraw, 2014; Giraudeau et al., 2015). The influence of urbanization on telomere length has only been investigated a few times in wild animals (Salmón et al., 2016; Stauffer et al., 2017), and further studies are necessary to confirm that urbanization could lead to shortened telomeres in developing wild animals.

In this study, we investigated the influence of urbanization on growth, offspring survival, offspring carotenoid-based coloration and offspring telomere length in a wild bird species, the great tit (Parus major). To our knowledge, no study has concomitantly investigated the influence of urbanization on all these variables and their functional links remain unclear in an urbanization context. Moreover, most studies have only compared a rural and an urban population, preventing to assess the influence of the degree of urbanization on these variables. Here, we examined this question by comparing four populations of great tits. In addition to two rural sites, we studied two urban sites with contrasted degrees of urbanization: a very large city (Paris), and a medium-sized city (Niort) in France, Western Europe. Great tit is a relevant species for such a study because it is considered as an urban adapter species (Shochat et al., 2006; Kark et al., 2007), and it can be found in rural sites, but also all along the urbanization gradient (Croci et al., 2008). Moreover, as cavity-nesters, the development of the chicks is very easy to follow (Lambrechts et al., 2010). According to previous studies, we predict that urbanization will be detrimental to developing chicks. First, we predict that growth will be impaired in cities relative to rural areas. Second, we predict, that urban nestlings will be less colorful than rural ones because of a lack of carotenoid-rich diet in cities. Third, we predict that urban nestlings will have shorter telomeres due to increased developmental stress than rural ones. Importantly, we predict that these effects should be reinforced as the degree of urbanization increases. If telomere length and carotenoid-based coloration are two proxies of health status and individual quality, we finally predict that carotenoid-based coloration and telomere length will be positively correlated with growth.

Materials and Methods

Study Sites and Data Collection

Data were collected in 2014 in two pairs of rural and urban populations of greattits breeding in France: (i) in the Paris area, and more specifically in city parks of Paris (48° 17′ N, 2° 38′ E) and in the forest surrounding the Foljuif field station, 80 km down south of Paris (CEREEP; 48° 16′ N, 2° 41′ E), and (ii) in the Niort area, and more specifically in the city of Niort (46°18′ N, 0°28′ W) and nearby forest of Chizé-CEBC, 30 km down south of Niort (46° 09′ N, 0° 24′ W). The rural study sites of Foljuif and Chizé respectively contained 120 and 100 nestboxes (Schwegler wood-concrete nest boxes 2 M, Valliance, Saint Pierre La Palud, France, floor area: 113 cm²) evenly distributed throughout homogenous deciduous woodlands composed mainly of oak (Quercus spp.), hornbeam (Carpinus betulus) and beech (Fagus sylvatica), and used by both blue and great tits (entrance hole diameter: 32 mm). The urban study sites of Paris and Niort contained respectively 60 and 64 nestboxes (Paris: 40 Schwegler nest boxes 2 M, and 20 handmade nestboxes in wooden panel, floor area: 105 cm², entrance hole diameter: 32 mm; Niort: 64 Schwegler nest boxes 2 M) installed in several urban parks or cemetery, separated from one another by 1–7 km as the crow flies. Paris is the biggest city in France (~2 million inhabitants), while Niort is a medium-sized city (~60 000 inhabitants), allowing us to investigate the potential impact of the degree of urbanization on nestlings' great tits. Paris is a very good example of historical European city with an extremely dense land occupation (housing development, industrial/commercial sites), and relatively small and isolated city Parks, while Niort is a classical “province” city less densely built with relatively larger city vegetation cover. Accordingly, these differences translated in highly divergent urbanization scores between our study sites as measured in another study (Meillère et al., 2017).

Nests were regularly inspected from March to June to determine laying date, clutch size, incubation date, hatching date of the first egg(s), number of hatchlings and unhatched eggs, brood size (at day 14 post-hatch), and number of successfully fledged nestlings. Only first clutches were used in this study. Data were collected for a total of 424 nestlings, originating from 28 nests in the rural sites (14 in Foljuif and 14 in Chizé) and from 26 nests in the urban sites (10 in Paris and 16 in Niort). Nestlings were captured at the nest when aged 14-days old, except for 2 broods aged 13 (one nest from Paris and one nest from Foljuif) and 4 broods aged 15 (two nests from Paris and two nests from Foljuif). Variation in nestling age being confounded with a pair of urban/rural sites and very low (mean nestling age ± s.e.: 14.02 ± 0.34), nestling age was not entered as a covariate in statistical models. Nestling capture spanned from May, 22th to June 10th (mean capture date ± s.e.: 123.63 ± 7.51; with 121 = May, 1st), there was no significant difference in phenology among sites [F_{(3, 50)} = 1.31 p = 0.28]. Tarsus length and wing length were measured to the nearest 0.1 mm with a caliper and a ruler, respectively, and body mass to the nearest 0.25 g with a Pesola spring balance. A sample of 5–8 yellow feathers was plucked from the center of the yellow breast for each bird, and stored in individual plastic bags in the dark until later color analysis. A blood sample (20–100 μl) was taken from the brachial vein in heparinized micro-haematocrit tubes, and centrifuged when back in the lab. Packed red blood cells were separated from plasma and stored at −20° C until later DNA extraction.

Telomere Measurements

Genomic DNA was extracted from frozen red blood cells using DNeasy Blood and Tissue Kit (Qiagen) according to the manufacturer's protocol. DNA yield and purity were checked using a spectrophotometer (Nanodrop ND-1000; Thermo Scientific, USA) according to previous recommendations (Nussey et al., 2014). All sample yields were >20 ng μl⁻¹. Moreover, ranges for absorption of all samples were within an acceptable range (between 1.8 and 2.0 for 260 nm/280 nm ratio and between 1.9 and 2.2 for 260 nm/230 nm ratio). Telomeres were measured using real-time quantitative PCR, following a protocol previously validated for birds and successfully used in great tits (Salmón et al., 2016). Briefly, a single-copy gene (glyceraldehyde-3-phosphate dehydrogenase, GAPDH) was amplified using the primers specifically designed for great tits (Atema et al., 2013). Telomeres primers were Tel1b (5′-CGGTTTGTTTGGGTTTGGGTTTGGGTTTGGGTTTGGGTT-3′) and Tel2b (5′-GGCTTGCCTTACCCTTACCCTTACCCTTACCCTTACCCT-3′). qPCR for both GAPDH and telomeres was performed using 5 ng of DNA per reaction. After preliminary optimization procedures, the telomere and GAPDH primers were respectively used at a final concentration of 1,000 and 200 nM. To control for the amplifying efficiency of the qPCR, each plate included a standard curve, which consisted of various concentrations of a great tit DNA pool (50–0.08 ng mL⁻¹). The efficiencies of the qPCR were within the acceptable range for both GAPDH and Telomere (Mean ± SD; GAPDH: 101.9 ± 1.09; Telomere: 96.7 ± 3.13). Telomere length is expressed relative to the single-copy gene (GAPDH) measured on the same sample of DNA (i.e., T/S ratio, Cawthon, 2002; Salmón et al., 2016). All standard DNA samples were run in duplicates. Moreover, a common sample was run on every plate and the inter-plate CV in TS ratio was 2.5%. Telomeres measurements were obtained for a subsample of 111 randomly chosen nestlings [number of nestlings (number of nests): Foljuif: n = 36(14), Paris: n = 24(10), Chizé: n = 30(12), Niort: n = 24(9)].

Feather Color Analysis

Yellow breast feather color was analyzed in the laboratory, using a spectrophotometer (Ocean Optics, IDIL Fibers Optiques, Lannion, France) following previously published methods (Andersson and Prager, 2006), samples being randomly chosen with respect to origin (study site and nest), and among nestlings from the same nest. Feathers were illuminated at an angle of 90° with a deuterium-tungsten lamp (DH2000), and reflected light was measured at the same angle, with a 200 μm fiber-optic reflection probe (QR200-7-SR-BX) connected to the spectrometer (USB2000+) operated using the Spectra Suite software. To exclude ambient light and standardize measuring distance (3 mm), an opaque black tube was fitted at the end of the measuring probe. Feather color was measured against a black short velvet background regularly cleaned and changed to keep its reflectance constant and equal to 0. A set of four feathers was randomly chosen among each sample, and piled up to reflect the bird's plumage (Quesada and Senar, 2006), with their rachis superimposed and parallel to each other. Preliminary tests showed that four feathers were sufficient to get reliable and repeatable spectra while keeping the number of feathers collected to a minimum for ethical reasons. Two sets of four feathers were analyzed for each individual with four independent (probe lift up and placed again on the sample) measures per set. Percentage of reflectance at each 1 nm interval was calculated between 300 and 700 nm, with respect to white (Spectralon diffuse reflectance Standard, WS-1-SL) and dark references, as R_(λ) = 100 × [(sample-white)/(white-dark)]. From these values, we calculated two independent and objective parameters to characterize carotenoid-based plumage color. Brightness is the total intensity of light reflected by feathers (spectral intensity), and was calculated as the integral of reflectance over [300, 700] nm (e.g., Endler, 1990; Andersson et al., 1998). Brightness reflects the structural quality of feathers and captures the UV component of carotenoid-based feather reflectance (Shawkey and Hill, 2005; Andersson and Prager, 2006). Plumage brightness has been shown to reflect body condition in this species (Galván, 2010). Carotenoid chroma reflects the carotenoid content of feathers and was calculated as the relative difference between the maximal (plateau at wavelength above 500 nm) and minimal (between 445 and 455 nm, corresponding to wavelengths of maximum absorbance of carotenoids) reflectance in the visible part of the spectrum: [R₇₀₀ − R₄₅₀]/R₇₀₀ (Örnborg et al., 2002; Andersson and Prager, 2006; Peters et al., 2007). Feather reflectance spectra were obtained by a single observer (MN) for 403 nestlings. Repeatability of measurements calculated as the intra-class correlation coefficient (Lessells and Boag, 1987) was always highly significant (p < 0.0001) with the following values: (a) repeatability within sets of feathers: brightness: 0.92, carotenoid chroma: 0.75, (b) and repeatability within individuals: brightness: 0.85, carotenoid chroma: 0.63. Average values for the 8 measures per nestling were used in subsequent statistical analyses.

Statistical Analyses

Statistical analyses were run using SAS v9.3 (SAS Institute Inc., Cary, NC, USA). Values are reported as mean ± s. e.

Generalized linear models were used to investigate variation in clutch size, number of hatchlings and brood size (number of nestlings surviving to the day of capture), as a function of site, including date as a covariate (GENMOD procedure with a Poisson distribution and a log link function).

Mixed linear models (MIXED procedure) were used to analyze variation in nestling phenotype (tarsus length, body mass, wing length, plumage color, and telomere length) with nest as a random effect to account for non-independence of nestlings raised in the same nest and sharing genetic and common rearing environment effects, using an unstructured covariance and REML (Restricted/Residual Maximum Likelihood) estimation method. The Wald test of the covariance parameter assessed whether the variance was significantly structured by the random nest effect. The fixed part of the models initially included the effects of site, date and the interaction between site and date. Models were compared with Akaike's Information Criterion (AIC), and the most parsimonious was retained (lowest AIC, Burnham and Anderson, 1998). Final models including significant differences between sites were followed by comparisons of means or least square means with adjusted values of p, using the Tukey-Kramer method. Tests of the residuals for normality and homoscedasticity were used to check the validity of the models; feather brightness and carotenoid chroma were log-transformed to meet these constraints. Finally, mixed linear models were also used to investigate whether plumage color and telomere length were related to nestling body mass and growth conditions (date and brood size) and whether telomere length was related to nestling plumage color. As there were significant differences between sites in brood size, nestling body condition and feather carotenoid chroma (see Section Results), these models were run separately for each site to avoid multi-collinearity issues (Graham, 2003).

Results

Clutch size differed significantly among sites [F_{(3, 50)} = 8.76, p = 0.03] and was significantly greater in Chizé as compared to Paris (Table 1), but there was no difference among sites in the number of hatchlings [F_{(1, 50)} = 5.35, p = 0.15]. Date and its interaction with site were not retained in the final models for clutch size and the number of hatchlings. Brood size, i.e., the number of nestlings surviving to 14-days old, decreased with advancing season [F_{(1, 49)} = 4.38, p = 0.04, slope estimate ± s.e.: −0.008 ± 0.004] and differed among sites [F_{(1, 49)} = 3.19, p = 0.03], with great tits in Foljuif rearing more nestlings to fledging than in Niort and Paris, and in Chizé as compared to Paris (Table 1).

TABLE 1

Table 1. Mean ± s.e. great tit clutch size (total number of eggs laid), number of hatchlings and brood size (number of nestlings surviving to 14 days old) in the four study sites.

Nestling growth, i.e., tarsus length, body mass, and wing length, significantly differed among sites (Table 2A, Figure 1): Nestlings were significantly taller (longer tarsi) in the rural site of Foljuif as compared to the other rural site of Chizé, and both urban sites of Niort and Paris. Nestlings were also smaller in Chizé and Niort as compared to Paris (Figure 1A). Nestlings were heavier in the rural site of Foljuif as compared to Chizé, Niort, and Paris (Figure 1B). Wing length was greater both in Foljuif and Paris as compared to Chizé and Niort (Figure 1C). Regarding plumage color, while yellow feather brightness did not differ among sites, there were significant differences in feather carotenoid chroma (Table 2B, Figure 2): yellow feathers were significantly more chromatic in both woodland sites of Foljuif and Chizé as compared to both urban sites of Paris and Niort. Finally, telomere length did not differ among sites (Table 2B, Figure 3). Nestling growth, plumage color, and telomere length did not vary with season in any site (date and its interaction with site were not retained in the final models).

TABLE 2

Table 2. Mixed linear models investigating variation in great tit nestlings' phenotype as a function of site and season.

FIGURE 1

Figure 1. Variation of (A) tarsus length, (B) body mass, and (C) wing length, among breeding sites in great tit nestlings. Nestling tarsus length was significantly greater in the rural site of Foljuif as compared to the other rural site of Chizé, and to both urban sites of Niort and Paris [differences of least squares means: t_(df:370) = 11.67, adj. p < 0.0001; t_(df:370) = 11.58, adj. p < 0.0001; and t_(df:370) = 4.38, adj. p < 0.0001, respectively]. Nestlings were also smaller in Chizé and Niort as compared to Paris [t_(df:370) = −6.20, adj. p < 0.0001; and t_(df:370) = −5.95, adj. p < 0.0001, respectively]. Nestlings were significantly heavier in the rural site of Foljuif as compared to Chizé, Niort, and Paris [differences of least squares means: t_(df:370) = 3.56, adj. p = 0.002; t_(df:370) = 4.03, adj. p = 0.0004; and t_(df:370) = 2.54, adj. p = 0.05, respectively]. Wing length was significantly greater both in Foljuif and Paris as compared to Chizé [t_(df:369) = 4.19, adj. p = 0.0002; and t_(df:369) = 3.18, adj. p = 0.01] and Niort [t_(df:369) = 3.84, adj. p = 0.001; t_(df:369) = 2.82, adj. p = 0.03], respectively. Different letters indicate significant differences.

FIGURE 2

Figure 2. Variation of yellow feather color among breeding sites in great tit nestlings: (A) brightness and (B) carotenoid chroma. Nestling feathers were significantly more chromatic in both woodland sites of Foljuif and Chizé as compared to both urban sites of Paris [t_(df:349) = 3.68, adj. p = 0.001; t_(df:349) = 4.76, adj. p < 0.0001] and Niort [t_(df:349) = 3.32, adj. p = 0.005; t_(df:349) = 4.57, adj. p < 0.0001], respectively. Different letters indicate significant differences.

FIGURE 3

Figure 3. Variation of telomere length among breeding sites in great tit nestlings reared in Chizé (woodland) and Niort (urban), Foljuif (woodland), and Paris (urban).

Yellow feather brightness was not related to nestling growth conditions or body mass in any site (all F < 1.87 and p > 0.17). Similarly, carotenoid chroma was not related to nestling growth conditions or body mass in Niort and Foljuif (all F < 2.16 and p > 0.14, Figures 4B,C). However, carotenoid chroma was positively related to nestling body mass [F_{(1, 40)} = 8.02, p = 0.007, estimate ± s.e. = 0.59 ± 0.02, Figure 4D], while a marginal negative relationship was found in Chizé [F_{(1, 98)} = 3.51, p = 0.06, estimate ± s.e. = −0.02 ± 0.01, Figure 4A]. In Foljuif and Niort, telomere length was not related to nestling body mass or growth conditions (all F < 0.59 and p > 0.45, Figures 5B,C). However, telomere length was positively related to nestling body mass in Chizé [F_{(1, 17)} = 7.71, p = 0.01, estimate ± s.e. = 0.02 ± 0.01], while the opposite pattern was found in Paris [F_{(1, 13)} = 5.45, p = 0.04, estimate ± s.e. = −0.02 ± 0.01; Figures 5A,D]. Finally, telomere length was not related to yellow feather color in any site (all F < 4.19 and p > 0.07).

FIGURE 4

Figure 4. Relationships between feather carotenoid chroma and body mass in great tit nestlings reared in (A) Chizé and (B) Niort, (C) Foljuif, and (D) Paris. For information, the marginally negative relationship in Chizé (see text for details), turned significant when the seemingly outlier data point is excluded [F_{(1, 97)} = 8.38, p = 0.005, estimate ± s.e. = −0.03 ± 0.01). There was however no biological or methodological reason to exclude a posteriori this nestling from the data set.

FIGURE 5

Figure 5. Relationships between telomere length and body mass in great tit nestlings reared in (A) Chizé and (B) Niort, (C) Foljuif, and (D) Paris (see text for details of the statistics).

Discussion

In this study, we aimed to assess the impact of urbanization on the phenotype and quality of developing offspring by comparing four populations of great tits. To do so, we focused on complementary morphological and molecular components (growth, plumage coloration, and telomere length). According to our predictions, we found that the body size and the carotenoid-based plumage of great tits chicks were affected by urbanization, although the impact of urbanization on body size was only found in the most urbanized site. However, we did not find any strong evidence for an impact of urbanization on telomere length.

Urbanization and Body Size in Nestling Great Tits

Overall and according to our prediction, we found some support for an effect of urbanization on nestling body size and body mass. In the area of the biggest city (Paris), we specifically found that urban nestlings were smaller and lighter than rural ones (Foljuif). However, we did not find any difference in body size and body mass between urban (Niort) and rural (Chizé) nestlings in the area of the smallest city. Supporting these contradictory patterns, some studies previously found that urban great tit nestlings are smaller than rural ones (Riddington and Gosler, 1995; Eeva et al., 2009; Geens et al., 2009) while others did not report such a difference (Isaksson and Andersson, 2007). All together, these results suggest that the influence of urbanization on growth and body size may only be apparent when the degree of urbanization reaches an upper threshold, or/and when urbanization is associated with specific anthropogenic perturbations that vary from one locality to another (e.g., pollution, Eeva et al., 2009; Geens et al., 2009). Thus, the effect of urbanization on body size may only be minor and non-significant for great tits reared in small to medium-sized city, such as Niort (~70,000 inhabitants) while it may become important and significant when they grow in large cities, such as Paris (over 2 million inhabitants). Similarly, some studies found an effect of urbanization on body size or growth in other bird species (Heiss et al., 2009; Bókony et al., 2012; Seress et al., 2012), supporting the idea that urbanization can have detrimental effects on growth and development. Interestingly, the influence of urbanization on body size also seems to depend on the degree of urbanization in some of these studies (Bókony et al., 2012; Meillère et al., 2017). For example, Meillère et al. (2017) found that urban house sparrows were not significantly smaller than rural ones when living in a medium-sized city such as Niort, while their body size was significantly reduced when living in bigger cities. To our knowledge, most studies have so far focused on a single pair of sites when comparing urban and rural populations of wild birds. Our study highlights the importance of studying multiple populations all over the urbanization gradient in order to better assess the influence of the degree of urbanization on the morphology of wild vertebrates (see Evans et al., 2009; Bókony et al., 2012; Meillère et al., 2017 for some examples).

At the proximate levels, our results support the idea that developmental conditions may be impoverished in big cities relative to forest areas. Overall, cities are characterized by multiple severe environmental modifications that certainly raise new important developmental constraints for wild vertebrates (Clergeau et al., 2006; McKinney, 2008; Gil and Brumm, 2014). Although, we did not quantify abiotic and biotic constraints in our populations, our study suggests that food abundance could explain the patterns we report because the growth of great tit nestlings has been directly correlated with food supply (Naef-Daenzer and Keller, 1999; Mägi et al., 2009, but see Tremblay et al., 2005). Supporting this interpretation, we also found that clutch size and brood size (number of nestlings reaching 14-days old) were generally reduced in cities relative to woodlands. Interestingly, we found that brood sizes differ between Paris and Foljuif, but not between Niort and Chizé, supporting further the idea that the degree of urbanization may have an impact on the ability of parent great tits to rear numerous young. During the chick-rearing period, great tits rely mainly on arthropods to feed their young and arthropod biomass and diversity is overall reduced in cities (McIntyre, 2000; Helden et al., 2012). Because nestlings rely on protein-rich diet during their development, this arthropod impoverishment is probably a major nutritional constraint for urban developing great tits (McIntyre et al., 2001; Shochat et al., 2004). Supporting this idea, food supplementation experiments have been shown to result in increased body size and nestling survival in urban birds (Zanette et al., 2003; Schoech et al., 2007; Heiss et al., 2009; Meyrier et al., 2017) and cross-fostering experiments demonstrated that reduced nestling body size in urban areas is indeed related to the rearing urban environment (e.g., Seress et al., 2012).

Urbanization and Plumage Characteristics

According to our prediction, we found a strong difference between the plumages of urban and rural nestlings: nestlings from Paris and Niort had less chromatic plumages, indicating their feathers contained less carotenoids (e.g., Peters et al., 2007) than those from Foljuif and Chizé. However, we did not report any effect of urbanization on plumage brightness, i.e., feather structural quality. These results are supported by previous studies that have reported similar results in great tits (Hõrak et al., 2000; Isaksson et al., 2005; Isaksson and Andersson, 2007) and other wild bird species (e.g., Jones et al., 2010; Giraudeau et al., 2015). In wild passerines, the diet seems to be the main factor affecting plumage coloration (Tschirren et al., 2003; Isaksson et al., 2007; Eeva et al., 2009; Giraudeau et al., 2015), and therefore, nestling plumage coloration certainly depends on the ability of parents to deposit carotenoids into their eggs (Biard et al., 2005, 2007) and to provide their chicks with carotenoid-rich diet (Biard et al., 2006; Isaksson et al., 2006). Our results suggest therefore that carotenoid availability and/or carotenoid uptake is reduced in the urban environment (Isaksson and Andersson, 2007; Isaksson et al., 2007). This interpretation is supported by eco-toxicological studies that were conducted in great tit nestlings: Eeva et al. (1998) and Dauwe and Eens (2008) found that nestlings had a paler plumage when they were raised in more polluted areas, and this effect was mediated by a reduced availability of carotenoid-rich diet in polluted areas (Sillanpää et al., 2008; Eeva et al., 2009). Furthermore, we found that the nestlings from Paris had a similar plumage color as those from Niort despite a large difference in the degree of urbanization between these two cities. Interestingly, this suggests that nestling plumage coloration is not strongly dependent on the degree of urbanization of a city, and more importantly that carotenoid availability or/and uptake by nestlings may not differ between medium-sized and large cities. The diet of nestlings surely differs to a large extent between urban and rural habitats (Isaksson and Andersson, 2007) as shown by contrasted body composition between urban and rural wild birds (e.g., yolk and plasma fatty acids, Andersson et al., 2015; Toledo et al., 2016). However, our results also suggest that the carotenoid availability and/or uptake may only slightly differ within a given habitat type (e.g., between different urban habitats, between different forests). To our knowledge, our study is the first to compare the plumage coloration of two urban populations, but another study did not report any difference in plumage coloration between woodlands differing in size (Ferns and Hinsley, 2008). Similarly, egg yolk fatty acid composition did not significantly differ between different types of woodlands, confirming therefore that the diet may not dramatically differ between different types of woodlands (Toledo et al., 2016).

We did not find any correlation between body mass and plumage coloration in Niort, Foljuif and Chizé (although the relationship was marginally significant in Chizé for feather carotenoid chroma). This suggests that growth and plumage color are relatively independent and may be governed by different pathways, as previously suggested by Hõrak et al. (2000). Actually, carotenoid availability appears as the main driver of plumage coloration whereas structural growth is probably rather related to nutritional conditions that could be independent of carotenoid contents. In contrast to Niort, Chizé, and Foljuif, body mass was positively correlated with plumage carotenoid chroma in Paris (i.e., the most urbanized population in our study). In Paris, the largest chicks were also the most colorful. This suggests that growth and plumage coloration might be simultaneously reduced when nutritional conditions are especially constraining (see also Senar et al., 2003). Supporting the idea that plumage coloration can be related to the nutritional status of individuals under some circumstances, previous studies found that the plumage of infected individuals was paler than that of healthy ones (Møller et al., 2000; Hõrak et al., 2001). Similarly, Tschirren et al. (2003) experimentally demonstrated a positive relationship between growth conditions and plumage coloration in great tit nestlings. Therefore, our study highlights a complex and context-dependent relationship between growth and plumage coloration in developing nestlings (Hõrak et al., 2000; Tschirren et al., 2003; Isaksson and Andersson, 2007). Under constraining environmental conditions, parents may decide to allocate most of the food and most of the carotenoid-rich items to the healthier nestlings, resulting therefore in a positive relationship between growth and plumage coloration.

In the highly-urbanized site (Paris), we clearly found that the development of great tit nestlings was detrimentally impacted (smaller structural size and paler plumage). Although structural growth was not slowed down in the moderately-urbanized site (Niort), we found that plumage color was paler relative to rural nestlings. Because both structural size and plumage color have been related to subsequent performances (e.g., Hõrak et al., 2001; Naef-Daenzer et al., 2001; Blount et al., 2003), our study highlights that urbanization may have multiple effects on developing wild birds. Further studies are now needed to better assess how these urban-related phenotypic modifications translate into fitness consequences.

Urbanization and Telomere Length

Contrary to our prediction, we did not find any effect of urbanization on telomere length in 14-days old great tit nestlings. Thus, we did not detect any difference in telomere length between the four populations (Paris, Niort, Foljuif, Chizé), suggesting no impact of urbanization on telomere length in nestling great tits. Conversely, Salmón et al. (2016) found that telomeres of urban great tit nestlings were shorter than those of rural nestlings and these findings were recently supported by another ecotoxicological study (Stauffer et al., 2017). However, Stauffer et al. (2017) also found that telomere length was not affected by urban-like pollution in nestlings from another passerine species, the pied flycatcher (Ficedula hypoleuca). To our knowledge, the impact of urbanization on telomere dynamics has so far rarely been studied in wild vertebrates, and taken all together, these studies suggest that the influence of urbanization on telomere length may depend on the species but also on the environmental characteristics of a given city. It also indicates that urbanization does not necessarily translate into significantly shorter telomeres in wild passerines.

Other hypotheses could also explain why urbanization affects telomere length in some studies but not others. First, these differences could result from regional differences in environmental conditions. In our study, we captured nestlings from four sites (two cities and two forests) that were located in Western Europe. In contrast, Salmón et al. (2016) and Stauffer et al. (2017) studied Northern populations that were located in Scandinavia. Regional differences in life-history strategies or environmental conditions could therefore potentially account for these contrasted results. Second, Niort is a medium-sized city, and thus urbanization may not have been intense enough to induce a strong effect on telomere length in that site. This hypothesis is however quite unlikely because we also did not detect any difference in telomere length in Paris, which is one of the most densely inhabited city in Europe, as compared to the other three sites studied. Moreover, Salmón et al., 2016) sampled great tits in Malmö, which is much smaller than Paris. Thirdly, selection processes could have masked a potential effect of urbanization on telomere length (Haussmann and Mauck, 2008). Urbanization is known to be associated with intense brood competition, and thus with brood reduction (reviewed in Chamberlain et al., 2009). Accordingly, we found that brood sizes were overall smaller in cities (Paris and Niort) relative to woodlands (Foljuif and Chizé), although the difference was marginally significant between Niort and Chizé. This brood reduction means that most of the costs are probably borne by a few nestlings while the others are only slightly disadvantaged. Supporting this hypothesis, Stier et al. (2015) found in great tits that the last-hatched nestling has shorter telomeres than its siblings at fledging. Therefore, reduced telomere length could have only occurred in the weakest nestlings that did not reach 14-days old, and were therefore not sampled in our study. This would however mean that urbanization translates into the production of fewer offspring and not necessarily into the production of fledglings with shorter telomeres. Lastly, we may also raise the hypothesis that the telomere length of nestlings may have been preserved at the cost of growth and plumage coloration in the cities we studied (Paris and Niort). Telomere attrition is known to result from oxidative stress (Monaghan, 2014) that is associated with rapid growth (Geiger et al., 2012) and reduced antioxidant protection in wild birds (Monaghan et al., 2009). Because the nestlings from Paris were smaller relative to those from Foljuif, the growth of the nestlings from Paris could have been reduced to limit the oxidative damages. Similarly, most carotenoids could have been allocated to antioxidant functions, and thus telomere protection, at the cost of plumage coloration in urbanized sites (Niort and Paris) although the antioxidant role of carotenoids is debated (Costantini, 2008; Isaksson and Andersson, 2008).

We found that telomere length was not related to growth in Niort and Foljuif, suggesting that developmental nutritional deficit does not necessarily translate into shorter telomeres. However, we found a positive relationship between growth and telomere length in Chizé, suggesting that higher quality nestlings also had longer telomeres. This finding is supported by several studies that have reported a positive relationship between body size and telomere length in chicks or young individuals (Caprioli et al., 2013; Angelier et al., 2015; Parolini et al., 2015; Mizutani et al., 2016) and it suggests that long telomeres may be a proxy of good nutritional conditions, at least under some circumstances. Surprisingly, we found the opposite pattern in Paris where larger individuals had shorter telomeres. This suggests that structural growth may be made at the expense of telomere protection in highly urbanized populations. Accordingly, previous studies have found that accelerated growth can be associated with oxidative stress and telomere attrition in wild birds (Geiger et al., 2012; Stier et al., 2015). Overall, it remains however unclear why the relationship between telomere length and body size/growth differs to such extent between sites. It emphasizes that the link between urbanization, environmental conditions and telomere dynamics is complex and deserves further studies.

Regional Differences in Nestlings' Morphology and Reproductive Output

Interestingly, we reported a large regional difference in body size between the Paris (Paris and Foljuif) and the Niort geographical areas (Niort and Chizé). Thus, the nestlings from the Paris area (Paris and Foljuif) were overall larger and bigger than those from the Niort area (Niort and Chizé). Although, this difference could result from a different way of measuring the nestlings between sites, this is quite unlikely because all the ringers in France are trained to use similar measurement techniques. Moreover, we found exactly the same pattern for tarsus length and wing length, but also for body mass that is not subject to a handler effect (body mass is measured with an electronic or spring scale). This supports therefore the idea that these differences do not result from different measurements techniques. However, we cannot totally exclude the possibility that some differences in morphological data have to do with variation in measurement techniques Despite this potential bias, our result is supported by the Bergman's rule (i.e., selection of larger body size as temperature decreases, and thus as latitude increases, Ashton, 2002) and by previous studies that have also reported an increase in body size as latitude increases in this species (Snow, 1954; Encabo et al., 2002). Alternatively, other explanations could also be involved in these structural differences between the two areas. For instance, predation risk, food resources and more generally environmental conditions are known to be associated with regional variations in body size (reviewed in Blanckenhorn, 2000). In our study, such hypothesis could explain the differences in body size between populations. For example, environmental conditions (food resources) could have been particularly favorable to nestlings' development (e.g., growth) in Foljuif relative to our other study populations.

More generally, this variation in body size might have important implications when focusing on the impact of urbanization on growth and nestlings' phenotypes because growing a larger structural size requires necessarily more energy and this may be especially challenging in an impoverished habitat, such as the urban one. Consequently, it suggests that the influence of urbanization on reproduction and nestlings' phenotypes might differ between populations of great tits. So far, the influence of urbanization on nestlings' phenotypes of great tits has mainly been studied in Northern environments (e.g., Riddington and Gosler, 1995; Isaksson and Andersson, 2007; Eeva et al., 2009; Geens et al., 2009). Further studies should therefore focus on multiple environmental conditions (populations spread over a large geographic range) in order to better assess the impact of urbanization not only on nestlings and reproductive performances, but also on adult wild birds.

Author Contributions

CB: Conceived the research, conducted field work, analyzed data, wrote the manuscript. FB: Conceived the research, conducted field work. AM: Conducted field work, provided comments on the manuscript. BM: Conducted field work. MN: Analyzed feather color. SR: Conducted molecular lab work. MV: Conducted field work, provided comments on the manuscript. FA: Conceived the research, conducted field work, wrote the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are grateful to the ANR (URBASTRESS, ANR-16-CE02-0004-01), the Centre National de la Recherche Scientifique (CNRS), the University of La Rochelle (ULR), the Région Poitou-Charentes, the Conseil Départemental des Deux-Sèvres (DC 79), the Contrat Plan Etat region (CPER ECONAT), and Sorbonne Universités—University Pierre et Marie Curie Paris 06 (Émergence-UPMC EME 1110) for funding this research. This work conforms to French legal requirements, including those relating to conservation and welfare, and was conducted under ringing licenses and authorizations for blood and feather sampling from the CRBPO (National Museum of Natural History) to CB and FA (research programs # 537 and 664, respectively). The protocol was approved by the COMETHEA ethical committee (CE2013-03). We would also like to acknowledge the support of the National Forest Office (ONF) and the cities of Paris and Niort, in allowing us to survey tit populations breeding in state-owned forest and city parks, respectively. We also thank all the students for their help in conducting fieldwork and Janice Duru for her help in the lab with feather color analyses.

References

Andersson, M. N., Wang, H.-L., Nord, A., Salmon, P., and Isaksson, C. (2015). Composition of physiologically important fatty acids in great tits differs between urban and rural populations on a seasonal basis. Front. Ecol. Evol. 3:93. doi: 10.3389/fevo.2015.00093