<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Ecol. Evol.</journal-id>
<journal-title>Frontiers in Ecology and Evolution</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Ecol. Evol.</abbrev-journal-title>
<issn pub-type="epub">2296-701X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fevo.2017.00064</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Ecology and Evolution</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Elevated Immune Gene Expression Is Associated with Poor Reproductive Success of Urban Blue Tits</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Capilla-Lasheras</surname> <given-names>Pablo</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/379929/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Dominoni</surname> <given-names>Davide M.</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Babayan</surname> <given-names>Simon A.</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/59379/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>O&#x00027;Shaughnessy</surname> <given-names>Peter J.</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Mladenova</surname> <given-names>Magdalena</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/438299/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Woodford</surname> <given-names>Luke</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/445608/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Pollock</surname> <given-names>Christopher J.</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Barr</surname> <given-names>Tom</given-names></name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Baldini</surname> <given-names>Francesco</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Helm</surname> <given-names>Barbara</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/426456/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Centre for Ecology and Conservation, College of Life and Environmental Sciences, University of Exeter</institution> <country>Penryn, United Kingdom</country></aff>
<aff id="aff2"><sup>2</sup><institution>Institute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow</institution> <country>Glasgow, United Kingdom</country></aff>
<aff id="aff3"><sup>3</sup><institution>Department of Animal Ecology, Netherlands Institute of Ecology</institution> <country>Wageningen, Netherlands</country></aff>
<aff id="aff4"><sup>4</sup><institution>Institute of Infection, Immunity and Inflammation, University of Glasgow</institution> <country>Glasgow, United Kingdom</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Caroline Isaksson, Lund University, Sweden</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Adam Michael Fudickar, Indiana University Bloomington, United States; Marta Szulkin, University of Warsaw, Poland</p></fn>
<fn fn-type="corresp" id="fn001"><p>&#x0002A;Correspondence: Pablo Capilla-Lasheras <email>p.capilla&#x00040;exeter.ac.uk</email></p></fn>
<fn fn-type="other" id="fn002"><p>This article was submitted to Behavioral and Evolutionary Ecology, a section of the journal Frontiers in Ecology and Evolution</p></fn></author-notes>
<pub-date pub-type="epub">
<day>16</day>
<month>06</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="collection">
<year>2017</year>
</pub-date>
<volume>5</volume>
<elocation-id>64</elocation-id>
<history>
<date date-type="received">
<day>25</day>
<month>02</month>
<year>2017</year>
</date>
<date date-type="accepted">
<day>31</day>
<month>05</month>
<year>2017</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2017 Capilla-Lasheras, Dominoni, Babayan, O&#x00027;Shaughnessy, Mladenova, Woodford, Pollock, Barr, Baldini and Helm.</copyright-statement>
<copyright-year>2017</copyright-year>
<copyright-holder>Capilla-Lasheras, Dominoni, Babayan, O&#x00027;Shaughnessy, Mladenova, Woodford, Pollock, Barr, Baldini and Helm</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract><p>Urban and forest habitats differ in many aspects that can lead to modifications of the immune system of wild animals. Altered parasite communities, pollution, and artificial light at night in cities have been associated with exacerbated inflammatory responses, with possibly negative fitness consequences, but few data are available from free-living animals. Here, we investigate how urbanization affects major immune pathways and experimentally test potentially contributing factors in blue tits (<italic>Cyanistes caeruleus</italic>) from an urban and forest site. We first compared breeding adults by quantifying the mRNA transcript levels of proteins associated with anti-bacterial, anti-malarial (TLR4, LY86) and anti-helminthic (Type 2 transcription factor GATA3) immune responses. Adult urban and forest blue tits differed in gene expression, with significantly increased <italic>TLR4</italic> and <italic>GATA3</italic>, but not <italic>LY86</italic>, in the city. We then experimentally tested whether these differences were environmentally induced by cross-fostering eggs between the sites and measuring mRNA transcripts in nestlings. The populations differed in reduced reproductive success, with a lower fledging success and lower fledgling weight recorded at the urban site. This mirrors the findings of our twin study reporting that the urban site was severely resource limited when compared to the forest. Because of low urban survival, robust gene expression data were only obtained from nestlings reared in the forest. Transcript levels in these nestlings showed no (<italic>TLR4, LY86</italic>), or weak (<italic>GATA3</italic>), differences according to their origin from forest or city nests, suggesting little genetic or maternal contribution to nestling immune transcript levels. Lastly, to investigate differences in parasite pressure between urban and forest sites, we measured the prevalence of malaria in adult and nestling blood. Prevalence was invariably high across environments and not associated with the transcript levels of the studied immune genes. Our results support the hypothesis that inflammatory pathways are activated in an urban environment and suggest that these differences are most likely induced by environmental factors.</p></abstract>
<kwd-group>
<kwd>urban ecology</kwd>
<kwd>inflammation</kwd>
<kwd>immunity</kwd>
<kwd>gene expression</kwd>
<kwd>blue tits</kwd>
<kwd>TLR4</kwd>
<kwd>LY86</kwd>
<kwd>GATA3</kwd>
</kwd-group>
<counts>
<fig-count count="3"/>
<table-count count="5"/>
<equation-count count="0"/>
<ref-count count="69"/>
<page-count count="13"/>
<word-count count="10373"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>Urban areas are the fastest growing land cover globally and are projected to increase a further 30% by 2030 (Seto et al., <xref ref-type="bibr" rid="B60">2012</xref>). The modification of the natural landscape caused by urbanization poses novel challenges to wildlife. For example, urbanization is associated with altered food webs, parasite communities, and chemical, light and noise pollution (Grimm et al., <xref ref-type="bibr" rid="B28">2008</xref>; Alberti, <xref ref-type="bibr" rid="B1">2015</xref>; Isaksson, <xref ref-type="bibr" rid="B33">2015</xref>). To thrive in these novel environments requires changes in behavior and physiology at the phenological and, possibly, genetic level. Such changes have, indeed, been found in urban populations of several species, indicating that they accommodate some of the challenges of city life (Isaksson, <xref ref-type="bibr" rid="B32">2010</xref>; Munshi-South and Kharchenko, <xref ref-type="bibr" rid="B44">2010</xref>; Dominoni et al., <xref ref-type="bibr" rid="B21">2013a</xref>, <xref ref-type="bibr" rid="B20">2016</xref>; Nemeth et al., <xref ref-type="bibr" rid="B47">2013</xref>; Atwell et al., <xref ref-type="bibr" rid="B4">2014</xref>; Gil et al., <xref ref-type="bibr" rid="B26">2014</xref>; Giraudeau et al., <xref ref-type="bibr" rid="B27">2014</xref>; Davies et al., <xref ref-type="bibr" rid="B18">2016</xref>; Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>). Nonetheless, urban environments are linked to reductions in fitness and health in a substantial number of organisms (Chamberlain et al., <xref ref-type="bibr" rid="B17">2009</xref>; Murray et al., <xref ref-type="bibr" rid="B45">2015</xref>). Hence, urban habitats are sometimes considered ecological traps which attract organisms, for example by anthropogenic food availability, but are insufficient for supporting successful rearing of offspring (Plummer et al., <xref ref-type="bibr" rid="B55">2013</xref>; Sumasgutner et al., <xref ref-type="bibr" rid="B63">2014</xref>; Lamb et al., <xref ref-type="bibr" rid="B34">2017</xref>). There is a need, therefore, to understand the mechanisms by which urbanization reduces reproductive success.</p>
<p>Despite their relevance for understanding urban adaptation, the proximate physiological mechanisms that underpin differences in behavior, life histories and fitness between urban and rural populations of wild animals still remain largely elusive (but see Dominoni et al., <xref ref-type="bibr" rid="B22">2013b</xref>; Atwell et al., <xref ref-type="bibr" rid="B4">2014</xref>; Davies et al., <xref ref-type="bibr" rid="B18">2016</xref>; Salm&#x000F3;n et al., <xref ref-type="bibr" rid="B58">2016</xref>). Whereas some physiological responses in urban animals have been investigated thoroughly, such as reproductive (Schoech et al., <xref ref-type="bibr" rid="B59">2004</xref>; Partecke et al., <xref ref-type="bibr" rid="B54">2006</xref>; Dominoni et al., <xref ref-type="bibr" rid="B22">2013b</xref>) and stress physiology (Bonier, <xref ref-type="bibr" rid="B13">2012</xref>), others have received little attention. Recently, a study using transcriptome analysis (RNA-seq) has broadly examined the ways physiological systems differ between an urban and rural population of a songbird, the great tit (<italic>Parus major</italic>) (Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>). Among the systems highlighted in this study, immunity showed particularly clear differences, corroborating reports of substantially reduced health of urban populations in several species (Martin et al., <xref ref-type="bibr" rid="B39">2010</xref>; Isaksson, <xref ref-type="bibr" rid="B33">2015</xref>; Murray et al., <xref ref-type="bibr" rid="B45">2015</xref>). Given the central role of the immune system for determining fitness and given its sensitivity to environmental factors, differences in immunity are thus strong candidates for linking urbanization to reduced reproductive success (Martin et al., <xref ref-type="bibr" rid="B39">2010</xref>; Isaksson, <xref ref-type="bibr" rid="B33">2015</xref>; Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>).</p>
<p>Urban environments can affect the immune system by modulating the structure of the parasite community to which animals are exposed and potentially respond to. Cities have been associated with altered parasite communities and altered infection dynamics (Fokidis et al., <xref ref-type="bibr" rid="B24">2008</xref>; Evans et al., <xref ref-type="bibr" rid="B23">2009</xref>; Giraudeau et al., <xref ref-type="bibr" rid="B27">2014</xref>; Neiderud, <xref ref-type="bibr" rid="B46">2015</xref>). For example, foraging behavior in cities is thought to enhance the prevalence of intestinal coccidians (Giraudeau et al., <xref ref-type="bibr" rid="B27">2014</xref>), potentially driven by increasing contact rates at feeders (Becker et al., <xref ref-type="bibr" rid="B10">2015</xref>). Several other factors associated with urban environments can affect the activation of particular immune pathways of wild animals. Prevalent urban conditions, such as chemical pollution and artificial light at night, can act as environmental stressors influencing the immune system and promoting pro-inflammatory immune profiles (Halliwell and Gutteridge, <xref ref-type="bibr" rid="B29">2002</xref>; Bedrosian et al., <xref ref-type="bibr" rid="B11">2011</xref>; Fonken et al., <xref ref-type="bibr" rid="B25">2013</xref>; Isaksson, <xref ref-type="bibr" rid="B33">2015</xref>; Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>). Low-quality and changed diet in urban areas can also enhance pro-inflammatory immune processes, which in turn can impose behavioral and physiological costs on animals (McGraw et al., <xref ref-type="bibr" rid="B40">2006</xref>; Andersson et al., <xref ref-type="bibr" rid="B3">2015</xref>; Becker et al., <xref ref-type="bibr" rid="B10">2015</xref>; Isaksson, <xref ref-type="bibr" rid="B33">2015</xref>). Furthermore, nutritional stress and resource restriction can exacerbate trade-offs between different biological processes, for example through resource allocation to growth vs. immunity (Sheldon and Verhulst, <xref ref-type="bibr" rid="B61">1996</xref>; Norris and Evans, <xref ref-type="bibr" rid="B49">2000</xref>). The findings of the transcriptomics study by Watson et al. (<xref ref-type="bibr" rid="B69">2017</xref>) showed that in an urban songbird population, where some of the above factors are effective, inflammatory immune genes were indeed upregulated. Experimental studies that measure candidate genes are now needed to clarify the origin of such differences, and links to fitness need to be established.</p>
<p>Here we examined how urban life affects avian immunity and fitness traits. As a study species, we used another parid songbird that is common in city and forest environments, the blue tit (<italic>Cyanistes caeruleus</italic>). Blue tits reportedly show reduced reproductive investment (e.g., clutch size and egg size) and reduced reproductive success in urban compared to rural populations (Chamberlain et al., <xref ref-type="bibr" rid="B17">2009</xref>; Bailly et al., <xref ref-type="bibr" rid="B7">2016b</xref>). We experimentally studied blue tits at two sites, one in an urban park area in Glasgow (UK), and one in a National Park 40 km away. The same populations were simultaneously investigated in a twin study on food resource availability, food provisioning and stable isotope signatures in eggs and blood samples (Pollock et al., <xref ref-type="bibr" rid="B56">2017</xref>). This parallel study design allowed us to link, for the same sites, comparisons of immunity with resource restriction. We first tested the hypothesis that adult urban blue tits show increased expression of genes associated with inflammation compared to blue tits from a forest habitat. Specifically, using RT-qPCR, we measured the transcript levels of TLR4 and LY86, which are involved in anti-bacterial and anti-malarial responses (Medzhitov, <xref ref-type="bibr" rid="B41">2001</xref>). We also tested the type 2 transcription factor GATA3, which is central to innate and adaptive immunity against parasitic helminths and immune homeostasis (Wang et al., <xref ref-type="bibr" rid="B67">2011</xref>; Tindemans et al., <xref ref-type="bibr" rid="B64">2014</xref>) and showed elevated gene expression in urban great tits (Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>).</p>
<p>We then investigated whether any urban-forest differences in immune gene transcript levels are already present during the nestling stage. Urban-forest immune differences in nestlings could arise either through environmental effects during early life, genetic or maternal effects. Thus, we cross-fostered eggs between and within the urban and forest sites to distinguish whether any difference in immune activation of nestlings was induced by the environment or controlled by genetic or maternal mechanisms (Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>). Finally, from the cross-fostering experiment, we quantified the reproductive success at the two sites and, thus, investigated possible associations between environment, immunity and reproductive success. We also examined prevalence of a common avian disease, avian malaria, as a potential driver of differences in expression of our immune markers (Martin et al., <xref ref-type="bibr" rid="B38">2014</xref>; Videvall et al., <xref ref-type="bibr" rid="B66">2015</xref>).</p>
<p>Increased levels of <italic>TLR4</italic> and <italic>LY86</italic> in adult urban birds are predicted on the basis of the hypothesis that the urban environment enhances inflammation. We also predicted increased levels of <italic>GATA3</italic> in adult urban birds as a consequence of high activation of immune pathways against intestinal parasites, whose transmission is thought to be exacerbated by anthropogenic food provisioning at bird feeders (Becker et al., <xref ref-type="bibr" rid="B10">2015</xref>). Expression of these three genes was also reported to be increased in at least one tissue in the transcriptome study on great tits (Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>). For nestlings blue tits, we predicted that urban-raised offspring might also show heightened expression of immune genes, indicating a direct effect of the environment on their immune profiles within their short postnatal life-span. We further predicted that cross-fostered nestlings should also differ in their immune gene expression profiles by origin, either due to micro-evolutionary change or based on differences in maternal investment. Assuming that higher expression levels in immune genes were selected for in urban birds, we expected higher expression levels in nestlings originating from eggs laid the urban environment regardless of their rearing environment. We also expected lower reproductive investment, fledging success and fledgling body mass in blue tits from the urban compared to the forest environment. The differences in reproductive success of urban and forest nests were indeed striking: 88% fledgling success in the forest compared to only 30% in the city. Consequently, in the cross-fostering experiment we obtained robust sample sizes on gene expression only for nestlings reared in the forest. We acknowledge limitations of our study arising from low sample sizes and from a design that involves only two sites (one urban and one rural) in a single year. Our experimental approach, using cross-fostering, addresses some of these concerns by effectively providing within-site replication. However, urban ecology will require studies in multiple sites and years, and on a broad range of organisms, before major, generalizable advances can be achieved (Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>).</p>
</sec>
<sec sec-type="materials and methods" id="s2">
<title>Materials and methods</title>
<sec>
<title>Experimental design and field protocol</title>
<p>Field work was carried out in one urban and one forest location in Scotland in May and June of 2014. In both sites, existing nest-box study systems were used (woodcrete boxes: 260 &#x000D7; 170 &#x000D7; 180 mm, &#x000D8; &#x0003D; 32 mm, Schwegler, Germany). The urban site was located in Kelvingrove Park, Glasgow (55&#x000B0; 52.18N 4&#x000B0; 17.22W), with a total of 60 bird nest-boxes. The forest location was situated in oak woodland near the Scottish Centre for Ecology and the Natural Environment (SCENE; 56&#x000B0; 7.73N 4&#x000B0; 36.79W), with a total of 143 bird nest-boxes. Twenty clutches of blue tits in the city and 22 in the forest were manipulated prior to clutch completion: 10 nests in the city and 12 in the forest were swapped within locations, representing control nests; and 10 nests were swapped across locations, representing experimentally cross-fostered nests (in total, 42 manipulated nests). Clutches within and across locations were matched based on sixth-egg laying date. Before swapping clutches, every egg involved in the experiment was individually marked, weighed (&#x000B1;0.01grams) and kept at 4&#x000B0;C overnight. Clutch size was reduced to six viable eggs at both sites in order to control for possible inter-habitat differences in clutch size. After females laid the sixth egg, clutches were swapped as explained above. When more than six eggs were present in a nest on the day of swapping, the six experimental eggs were randomly chosen. After swapping clutches, nests were checked every other day and newly laid eggs were replaced by dummy eggs. Using dummy eggs, we always kept the original number of eggs laid by females (each nest contained six real eggs plus a variable number of dummy eggs depending on the number of eggs that females actually laid). Dummy eggs were removed from nests after real eggs hatched. The total number of eggs laid by each female was recorded and termed &#x0201C;natural clutch size.&#x0201D;</p>
<p>After the 10th day of incubation, nests were monitored daily and hatch date was precisely assigned for every nest. On day 13 after first-egg hatching, blue tit nestlings were weighed (&#x000B1;0.05 grams), ringed with a unique metal ring and between 20 and 75 &#x003BC;L of blood were collected from their brachial veins (stored in 250 &#x003BC;L of RNAlater&#x000AE; for gene expression analysis). Additionally, between 20 and 75 &#x003BC;L of whole blood were stored in &#x0003E;99% ethanol for molecular screening of malaria parasites. Between 10 and 12 days post-hatching, we aimed to capture one or two of the respective parents while provisioning and to sample adults as described for chicks. However, when we realized the low overall breeding success in our urban population we greatly reduced efforts of catching parents to avoid additional risks to the broods. In total, we obtained samples of 24 adult birds (Table <xref ref-type="table" rid="T1">1</xref>). Nest-boxes were checked 16&#x02013;20 days after expected fledging dates in search of dead nestlings. Since clutch size was reduced to six eggs, hatching success was defined as number of hatchlings divided by six. Number of fledglings over number of hatchlings represented fledging success. Table <xref ref-type="table" rid="T1">1</xref> summarizes sample sizes throughout different breeding stages between our urban and forest site, which became progressively disparate because of high chick mortality in the city.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Summary of sample sizes that entered the analysis: number of clutches, nestlings and adults of each experimental group at every breeding stage for our two study sites.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left" colspan="5"><bold>Initial number of nests</bold></th>
</tr>
</thead>
<tbody>
<tr style="border-bottom: thin solid #000000;">
<td/>
<td/>
<td valign="top" align="center"><bold>Urban</bold></td>
<td valign="top" align="center"><bold>Forest</bold></td>
<td valign="top" align="center"><bold>Total</bold></td>
</tr>
<tr>
<td valign="top" align="left" colspan="2">Cross-fostered group</td>
<td valign="top" align="center">10</td>
<td valign="top" align="center">10</td>
<td valign="top" align="center">20</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left" colspan="2">Control group</td>
<td valign="top" align="center">10</td>
<td valign="top" align="center">12</td>
<td valign="top" align="center">22</td>
</tr>
 <tr>
<td valign="top" align="left" colspan="2">Total</td>
<td valign="top" align="center">20</td>
<td valign="top" align="center">22</td>
<td valign="top" align="center">42</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left" colspan="5"><bold>Number of clutches successfully hatched</bold><xref ref-type="table-fn" rid="TN1"><sup>a</sup></xref></td>
</tr>
 <tr>
<td valign="top" align="left" colspan="2">Cross-fostered group</td>
<td valign="top" align="center">7</td>
<td valign="top" align="center">6</td>
<td valign="top" align="center">13</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left" colspan="2">Control group</td>
<td valign="top" align="center">8</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">17</td>
</tr>
 <tr>
<td valign="top" align="left" colspan="2">Total</td>
<td valign="top" align="center">15</td>
<td valign="top" align="center">15</td>
<td valign="top" align="center">30</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left" colspan="5"><bold>Number of hatchlings</bold><xref ref-type="table-fn" rid="TN2"><sup>b</sup></xref></td>
</tr>
 <tr>
<td valign="top" align="left" colspan="2">Cross-fostered group</td>
<td valign="top" align="center">32</td>
<td valign="top" align="center">25</td>
<td valign="top" align="center">57</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left" colspan="2">Control group</td>
<td valign="top" align="center">30</td>
<td valign="top" align="center">39</td>
<td valign="top" align="center">69</td>
</tr>
 <tr>
<td valign="top" align="left" colspan="2">Total</td>
<td valign="top" align="center">62</td>
<td valign="top" align="center">64</td>
<td valign="top" align="center">126</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left" colspan="5"><bold>Number of fledglings</bold></td>
</tr>
 <tr>
<td valign="top" align="left" colspan="2">Cross-fostered group</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">20</td>
<td valign="top" align="center">29</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left" colspan="2">Control group</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">36</td>
<td valign="top" align="center">45</td>
</tr>
 <tr>
<td valign="top" align="left" colspan="2">Total</td>
<td valign="top" align="center">18</td>
<td valign="top" align="center">56</td>
<td valign="top" align="center">74</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left" colspan="5"><bold>Number of individuals sampled for RT-qPCR analysis</bold><xref ref-type="table-fn" rid="TN3"><sup>c</sup></xref></td>
</tr>
 <tr>
<td valign="top" align="left">Nestlings</td>
<td valign="top" align="left">Cross-fostered group</td>
<td valign="top" align="center">4</td>
<td valign="top" align="center">16</td>
<td valign="top" align="center">20</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td/>
<td valign="top" align="left">Control group</td>
<td valign="top" align="center">4</td>
<td valign="top" align="center">24</td>
<td valign="top" align="center">29</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td/>
<td valign="top" align="left">Total</td>
<td valign="top" align="center">8</td>
<td valign="top" align="center">40</td>
<td valign="top" align="center">48</td>
</tr>
 <tr>
<td valign="top" align="left">Adults</td>
<td/>
<td valign="top" align="center">13</td>
<td valign="top" align="center">11</td>
<td valign="top" align="center">24</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>Cross-fostered clutches were swapped between sites, control clutches were swapped within sites</italic>.</p>
<fn id="TN1">
<label>a</label>
<p><italic>Clutches with at least one egg hatched</italic>.</p></fn>
<fn id="TN2">
<label>b</label>
<p><italic>Maximum number of hatchlings per clutches was six since every clutch was experimentally reduced to six eggs</italic>.</p></fn>
<fn id="TN3">
<label>c</label>
<p><italic>With successful gene expression data for both reference genes used in data normalization</italic>.</p></fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec>
<title>Ethics statement</title>
<p>All bird sampling was conducted following the directions and legislations of UK Home Office (project license: 70/7899 to BH, and personal licenses to DD and BH), Scottish Natural Heritage (52463 to BH) and British Trust for Ornithology (Scientific C and T licenses to BH and DD, respectively).</p>
</sec>
<sec>
<title>Expression of immune and reference genes</title>
<sec>
<title>Primer design</title>
<p>Blue tit sequences of <italic>TLR4</italic> (toll-like receptor 4), <italic>LY86</italic> (lymphocyte antigen 86, also known as MD1), <italic>GATA3</italic> (GATA Binding Protein 3) and of the candidate reference genes, <italic>HPRT</italic> (Hypoxanthine-guanine phosphoribosyl transferase), <italic>PMM1</italic> (Phosphomannomutase 1), <italic>SDHA</italic> (succinate dehydrogenase complex, subunit A) and <italic>TFRC</italic> (transferrin receptor protein 1) (Olias et al., <xref ref-type="bibr" rid="B53">2014</xref>) were obtained from the blue tit genome (Mueller et al., <xref ref-type="bibr" rid="B42">2016</xref>). Gene sequences were compared in BLAST against the zebra finch (<italic>Taeniopygia guttata</italic>) genome (Warren et al., <xref ref-type="bibr" rid="B68">2010</xref>) in order to design primers on the correct gene regions. Primers were designed using Primer Express&#x02122; 2.0.0 (Applied Biosystems) as described previously (O&#x00027;Shaughnessy et al., <xref ref-type="bibr" rid="B51">2008</xref>). In order to avoid genomic DNA (gDNA) amplification, every primer pair was designed to flank an intron of more than 1,000 base pairs whenever possible. All primer sequences (Table <xref ref-type="supplementary-material" rid="SM1">S1</xref>), as well as the validation of reference genes (Figure <xref ref-type="supplementary-material" rid="SM1">S1</xref>), are detailed as Supplementary Material.</p>
</sec>
<sec>
<title>RNA extraction, reverse transcription and RT-qPCR protocol</title>
<p>Blood samples stored in RNAlater&#x000AE; were centrifuged for 5 min at 5,000 RPM, and RNA in the cell fraction isolated using TRIzol&#x000AE; reagent (Life Technologies, Thermo Fisher Scientific 2015). Extracted RNA was reverse transcribed using random hexamers and Moloney murine leukemia virus reverse transcriptase (Superscript III, Invitrogen Ltd.) (for detailed protocols see O&#x00027;Shaughnessy and Murphy, <xref ref-type="bibr" rid="B52">1993</xref>; O&#x00027;Shaughnessy et al., <xref ref-type="bibr" rid="B50">1994</xref>). Levels of mRNA were measured by real-time quantitative PCR using the SYBR Green method (O&#x00027;Shaughnessy and Murphy, <xref ref-type="bibr" rid="B52">1993</xref>; O&#x00027;Shaughnessy et al., <xref ref-type="bibr" rid="B50">1994</xref>). The efficiency of the seven employed primer pairs was between 91 and 120% in all the cases. None of the primer pairs amplified gDNA. One unique gene was run in each RT-qPCR 96-well plate along with a non-template control (NTC). Every sample was always run in duplicate and only samples with similar duplicate results were considered for final analysis. The vast majority of samples included in analyses had C<sub>t</sub> differences between replicates below one. Five samples showed C<sub>t</sub> differences between one and two, and two additional samples in the final analysis had replicate C<sub>t</sub> differences larger than two (2.07 and 2.18). The exclusion of the samples with C<sub>t</sub> differences above 2 or the samples with C<sub>t</sub> differences above 1 did not change the results. For details of RT-qPCR data normalization, see Supplementary Material.</p>
</sec>
</sec>
<sec>
<title>Screening of malaria parasites</title>
<p>Avian malaria collectively refers to blood parasites of the three taxa <italic>Leucocytozoon, Haemoproteus</italic> and/or <italic>Plasmodium</italic> [PMID:15357072]. We tested for the presence or absence of any of these parasites using a nested PCR approach which identified <italic>Leucocytozoon</italic> in one reaction, and either of <italic>Haemoproteus</italic> and <italic>Plasmodium</italic> in another reaction (Hellgren et al., <xref ref-type="bibr" rid="B31">2004</xref>). In short, DNA from blood samples stored in &#x0003E;99%-ethanol was extracted using a commercial kit (DNeasy whole-blood extraction kit, Qiagen). Eluted DNA was amplified by 20 cycles at 94&#x000B0;C for 30 s, 50&#x000B0;C for 30 s and 72&#x000B0;C for 45 s. A 2 &#x003BC;l-aliquot from the result of this reaction was further amplified for 35 cycles at 94&#x000B0;C for 30 s, 54&#x000B0;C for 30 s, and 72&#x000B0;C for 45 s. Both PCRs were performed in 20-&#x003BC;l reaction mixture using 10 &#x003BC;l of GoTaq&#x000AE; universal PCR master mix (Promega) and 0.6 &#x003BC;M of each primer. Products of the second PCR reaction were examined on a 1% agarose gel. The presence of <italic>Leucocytozoon</italic> and/or <italic>Haemoproteus</italic>/<italic>Plasmodium</italic> parasites was further checked by sequencing 24 PCR positive bands (12 for <italic>Leucocytozoon</italic> and 12 for <italic>Haemoproteus</italic>/<italic>Plasmodium</italic>, corresponding to samples positive for both taxa) (Eurofins Genomics). In all the positive samples that we tested, sequencing confirmed the presence of at least one avian malaria parasite, <italic>Leucocytozoon</italic> (GenBank Accession Numbers <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="KY981756-KY981765">KY981756-KY981765</ext-link>) and/or <italic>Haemoproteus</italic> (GenBank Accession Numbers <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="KY981753-KY981755">KY981753-KY981755</ext-link>). No <italic>Plasmodium</italic> was detected. In nine out of 12 samples that were positive for both <italic>Leucocytozoon</italic> and <italic>Haemoproteus</italic> by PCR, sequencing revealed that only one of the two parasites was actually present, suggesting that the observed double infections might not always be true but possibly an artifact of the PCR procedure due to cross-reaction of the primers.</p>
</sec>
<sec>
<title>Statistical analysis</title>
<sec>
<title>Model selection approach</title>
<p>For every modeled response variable, data analysis started with a global model including every biologically important predictor (see below). Then, model selection based on Akaike&#x00027;s Information Criterion (AIC, Burnham and Anderson, <xref ref-type="bibr" rid="B15">2002</xref>; Burnham, <xref ref-type="bibr" rid="B14">2004</xref>; Burnham et al., <xref ref-type="bibr" rid="B16">2011</xref>) was applied to investigate alternative hypotheses on the effects of urban&#x02013;forests conditions and of our experimental manipulation on immune gene expression, breeding investment and breeding success. Models were ranked based on their AICc and considered similarly supported if their &#x00394;AICc value was &#x0003C;2. Poisson and binomial models were checked for over-dispersion by comparing residual deviance and residual degrees of freedom. Linear models and linear mixed model residuals were inspected to check that they met the assumption of normality, and collinearity between pairs of explanatory variables was checked before accepting the results of any statistical model. Data analysis was carried out in R 3.3.2 (R Core Team, <xref ref-type="bibr" rid="B57">2016</xref>) using the <italic>lme4</italic> (Bates et al., <xref ref-type="bibr" rid="B9">2014</xref>) and <italic>MuMin</italic> (Barton, <xref ref-type="bibr" rid="B8">2016</xref>) packages.</p>
</sec>
<sec>
<title>Immune transcript levels and malaria prevalence in adult blue tits</title>
<p>After RT-qPCR data normalization, adult levels of <italic>TLR4, LY86</italic>, and <italic>GATA3</italic> were analyzed using linear models (LMs). Habitat (a factor with two levels: urban or forest) and adult body mass 10&#x02013;12 days after hatching of their clutch were originally included as predictors in every model. Weight information was missed for two adult birds with successful gene expression data. Weight appeared to have little importance in explaining variation in gene expression and we, therefore, present statistical results for models containing weight in Table <xref ref-type="supplementary-material" rid="SM1">S2</xref>. In the main text, we report results for models without weight but include the two additional birds to increase our sample size and, hence, the confidence in our estimates of the effect of habitat. Malaria prevalence was compared across habitats using Fisher&#x00027;s exact test. Because little variation in infection status was observed in adult birds, we were unable to investigate its association with immune transcript levels.</p>
</sec>
<sec>
<title>Immune transcript levels and malaria prevalence in cross-fostered nestlings</title>
<p>Originally, the main purpose of this analysis was to investigate the effect of the original and rearing environment on immune transcript levels. However, because of the high mortality in our urban environment (see Section Results), sample sizes from urban raised nestlings were very low. Hence, we focus our main analysis only on nestlings raised in the forest and present preliminary results for the full cross-fostering design as Supplementary Material.</p>
<p>Nestling TLR4, LY86, and GATA3 transcript levels were analyzed using linear-mixed models (LMMs). Original environment (a factor with two levels: urban or forest) and weight of chicks on day 13 were included as fixed effects. Malaria infection status (factor with two levels, YES/NO) was also included as a predictor. Nest ID was included as a random factor in every model as several nestlings were measured per nest.</p>
</sec>
<sec>
<title>Reproductive investment and breeding success of urban and forest cross-fostered nests</title>
<p>We compared several measures of reproductive investment and breeding success of adult females as follows. Female investment was quantified by number (i.e., natural clutch size) and size of eggs laid. Natural clutch size was modeled as a Poisson variable in a generalized linear model (GLM) including habitat, and clutch completion date (as a linear and a quadratic term) to account for temporal trends. We also included experimental group (cross-fostered or control) and its interaction with site to test whether our experimental design differed across habitats. Individual egg weight was modeled using LMMs including as explanatory variables original environment, experimental group, their interaction, natural clutch size and 1st egg laying date (as a linear and a quadratic term). Nest ID was always kept as random factor.</p>
<p>We then examined reproductive success of the urban and forest nests. Hatching and fledging success were analyzed with generalized linear mixed models (GLMMs) using binomial distributions with logit link functions. Due to over-dispersed model residuals, observation-level random factors were employed and yielded good model fit (Harrison, <xref ref-type="bibr" rid="B30">2015</xref>). Original and rearing environment, their interaction, and clutch completion date&#x02014;for hatching success&#x02014;or hatch date (linear and quadratic terms)&#x02014;for fledgling success&#x02014;were included as explanatory variables. Weight of fledglings on day 13 after hatching was analyzed by a LMM, keeping Nest ID as a random factor and using original and rearing environment, their interaction, hatch date (linear and quadratic terms), brood size on day 13, and the interaction between rearing environment and brood size as explanatory variables.</p>
</sec>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<sec>
<title>Immune transcript levels and malaria prevalence of urban and forest adults</title>
<p>Adult urban blue tits showed higher levels of <italic>TLR4</italic> and <italic>GATA3</italic> than forest blue tits (Figure <xref ref-type="fig" rid="F1">1</xref>). Including habitat as a predictor generated models for these two immune genes that were superior in AICc than intercept-only models (Table <xref ref-type="table" rid="T2">2</xref>). Removing a high-value outlier in GATA3 transcript levels of urban birds (Figure <xref ref-type="fig" rid="F1">1</xref>) did not qualitatively change the statistical outcome. After removal of this outlier, the model containing habitat as a predictor was still the most supported by the data, with a decrease in AICc of 2.4 compared to the intercept-only model. For <italic>LY86</italic> expression, we did not find statistical evidence for an alteration across habitats (Table <xref ref-type="table" rid="T2">2</xref>, Figure <xref ref-type="fig" rid="F1">1B</xref>). Body weight of adults 10&#x02013;12 days after hatching of their clutch was not an important predictor for transcript levels of any immune gene (Table <xref ref-type="supplementary-material" rid="SM1">S2</xref>).</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Relative expression of <bold>(A)</bold> <italic>TLR4</italic>, <bold>(B)</bold> <italic>LY86</italic>, and <bold>(C)</bold> <italic>GATA3</italic> for adult, breeding blue tits. Y-axis values refer to levels of gene expression relative to the expression of reference genes. Main bars illustrate raw data mean values and dark error bars illustrate &#x000B1; 2 standard errors. Dots represent raw data points. Sample sizes are illustrated at the base of each bar.</p></caption>
<graphic xlink:href="fevo-05-00064-g0001.tif"/>
</fig>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>Summarizing table of statistical models employed to explain variation in TLR4, LY86, and GATA3 transcript levels for adult blue tits.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th/>
<th valign="top" align="center" colspan="6" style="border-bottom: thin solid #000000;"><bold>Estimates of model coefficients (standard error)</bold></th>
</tr>
<tr>
<th valign="top" align="left"><bold>Response term</bold></th>
<th valign="top" align="center"><bold>Intercept</bold></th>
<th valign="top" align="center"><bold>Habitat &#x02013; urban<xref ref-type="table-fn" rid="TN4"><sup>a</sup></xref></bold></th>
<th valign="top" align="center"><bold><italic>k</italic></bold></th>
<th valign="top" align="left"><bold>AICc</bold></th>
<th valign="top" align="left"><bold>&#x00394;AICc</bold></th>
<th valign="top" align="left"><bold><italic>w</italic></bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">TLR4 gene expression (<italic>n</italic> &#x0003D; 16; 14 nest-boxes)</td>
<td valign="top" align="center">0.055 (0.033)</td>
<td valign="top" align="center">0.090 (0.041) <xref ref-type="table-fn" rid="TN5"><sup>&#x0002A;</sup></xref></td>
<td valign="top" align="center">3</td>
<td valign="top" align="left">&#x02212;29.5</td>
<td valign="top" align="left">0.0</td>
<td valign="top" align="left">0.69</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td/>
<td valign="top" align="center">0.111 (0.022)</td>
<td/>
<td valign="top" align="center">2</td>
<td valign="top" align="left">&#x02212;27.9</td>
<td valign="top" align="left">1.6</td>
<td valign="top" align="left">0.31</td>
</tr>
<tr>
<td valign="top" align="left">LY86 gene expression (<italic>n</italic> &#x0003D; 20; 17 nest-boxes)</td>
<td valign="top" align="center">0.464 (0.055)</td>
<td/>
<td valign="top" align="center">2</td>
<td valign="top" align="left">4.1</td>
<td valign="top" align="left">0.0</td>
<td valign="top" align="left">0.79</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td/>
<td valign="top" align="center">0.486 (0.083)</td>
<td valign="top" align="center">&#x02212;0.041 (0.112)</td>
<td valign="top" align="center">3</td>
<td valign="top" align="left">6.7</td>
<td valign="top" align="left">2.7</td>
<td valign="top" align="left">0.21</td>
</tr>
<tr>
<td valign="top" align="left">GATA3 gene expression (<italic>n</italic> &#x0003D; 22, 17 nest-boxes)</td>
<td valign="top" align="center">0.206 (0.077)</td>
<td valign="top" align="center">0.229 (0.106) <xref ref-type="table-fn" rid="TN5"><sup>&#x0002A;</sup></xref></td>
<td valign="top" align="center">3</td>
<td valign="top" align="left">6.5</td>
<td valign="top" align="left">0.0</td>
<td valign="top" align="left">0.72</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.321 (0.058)</td>
<td/>
<td valign="top" align="center">2</td>
<td valign="top" align="left">8.4</td>
<td valign="top" align="left">1.9</td>
<td valign="top" align="left">0.28</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>Models are presented along with their coefficient estimates and standard errors [estimate (SE)]. k, number of model parameters; w &#x0003D; Akaike&#x00027;s weight, defined as the probability of a model given the data and the candidate set of alternative models</italic>.</p>
<fn id="TN4">
<label>a</label>
<p><italic>&#x0201C;Habitat &#x02013; forest&#x0201D; set as reference level and, therefore, equals zero. When top models contained &#x0201C;habitat,&#x0201D; this predictor was assessed by a likelihood-ratio test against the intercept-only model and</italic></p></fn>
<fn id="TN5">
<label>&#x0002A;</label>
<p><italic>illustrates a p-value &#x0003C; 0.05. See Table <xref ref-type="supplementary-material" rid="SM1">S2</xref> for model estimates including adult weight</italic>.</p></fn>
</table-wrap-foot>
</table-wrap>
<p>Every screened adult urban bird scored positive for malaria infection (<italic>n</italic> &#x0003D; 6). Out of 9 forest birds tested, 8 were positive for malaria parasites. These results indicated no association between malaria prevalence and habitat in adult blue tits (Fisher&#x00027;s exact test, <italic>p</italic> &#x0003E; 0.90), suggesting that differential pressure of malaria parasites across environments does not underlay our results of immune gene expression.</p>
</sec>
<sec>
<title>Immune transcript levels and malaria prevalence of forest-reared nestlings</title>
<p>Our data suggested no differences in transcript levels of TLR4 and LY86 in forest-reared nestlings due to habitat of origin, malaria parasite infection or weight 13 days after hatching (Figure <xref ref-type="fig" rid="F2">2</xref>). For the expression of these two genes, intercept-only models were most supported by the data (Table <xref ref-type="table" rid="T3">3</xref>). The presence of malaria parasites did not predict <italic>TLR4</italic> or <italic>LY86</italic> levels. Malaria infection status was only kept in models with &#x00394;AICc values larger than 2. For <italic>GATA3</italic>, we found some evidence for an effect of the original environment on nestling transcript levels. Urban-originated nestlings reared in the forest had higher levels of <italic>GATA3</italic> than forest-originated forest-reared birds (Figure <xref ref-type="fig" rid="F2">2C</xref>). Four competing models scored very similar AICc values, two of them containing original environment and weight on day 13 as predictors (Table <xref ref-type="table" rid="T3">3</xref>). However, after removal of an urban-originated outlier with high <italic>GATA3</italic> (Figure <xref ref-type="fig" rid="F2">2C</xref>), the effect of weight and of the original environment lost much importance and the intercept-only model became the most supported one (Table <xref ref-type="supplementary-material" rid="SM1">S3</xref>). Results for the full cross-fostering experiment, including urban-reared forest birds, are shown in Figure <xref ref-type="supplementary-material" rid="SM1">S2</xref>.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Relative expression of <bold>(A)</bold> <italic>TLR4</italic>, <bold>(B)</bold> <italic>LY86</italic>, and <bold>(C)</bold> <italic>GATA3</italic> for forest-reared nestlings hatched from cross-fostered eggs of the urban or forest site. Y-axis values refer to transcript levels relative to the expression of reference genes. Main bars illustrate raw data mean values and dark error bars illustrate &#x000B1; 2 standard errors. Dots represent raw data points. Sample sizes are illustrated at the base of each bar. Removal of the urban outlier in <italic>GATA3</italic> analysis <bold>(C)</bold> led to qualitative and quantitative changes in the results regarding this gene (Table <xref ref-type="supplementary-material" rid="SM1">S3</xref>). Differences in sample size between Figure <xref ref-type="fig" rid="F2">2</xref> and Table <xref ref-type="table" rid="T1">1</xref> are due to two nestlings with missing information for body weight. The inclusion of these additional data points does not change the results regarding the effect of the rearing habitat.</p></caption>
<graphic xlink:href="fevo-05-00064-g0002.tif"/>
</fig>
<table-wrap position="float" id="T3">
<label>Table 3</label>
<caption><p>Summarizing table of statistical models employed to explain variation in TLR4, LY86, and GATA3 transcript levels for forest-reared nestlings.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th/>
<th valign="top" align="center" colspan="5" style="border-bottom: thin solid #000000;"><bold>Estimates of model coefficients (standard error)</bold></th>
<th/>
<th/>
<th/>
<th/>
</tr>
<tr>
<th valign="top" align="left"><bold>Response term</bold></th>
<th valign="top" align="center"><bold>Intercept</bold></th>
<th valign="top" align="center"><bold>Original Habitat &#x02013; urban<xref ref-type="table-fn" rid="TN6"><sup>a</sup></xref></bold></th>
<th valign="top" align="center"><bold>Malaria parasites<xref ref-type="table-fn" rid="TN7"><sup>b</sup></xref></bold></th>
<th valign="top" align="center"><bold>Weight</bold></th>
<th valign="top" align="center"><bold>Original Habitat &#x02013; urban<xref ref-type="table-fn" rid="TN6"><sup>a</sup></xref> X weight</bold></th>
<th valign="top" align="center"><bold><italic>k</italic></bold></th>
<th valign="top" align="center"><bold>AICc</bold></th>
<th valign="top" align="center"><bold>&#x00394;AICc</bold></th>
<th valign="top" align="center"><bold><italic>w</italic></bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">TLR4 gene expression (<italic>n</italic> &#x0003D; 28; 13 nest-boxes)</td>
<td valign="top" align="center">0.036 (0.004)</td>
<td/>
<td/>
<td/>
<td/>
<td valign="top" align="center">3</td>
<td valign="top" align="center">&#x02212;128.9</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.39</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.104 (0.060)</td>
<td/>
<td/>
<td valign="top" align="center">&#x02212;0.006 (0.005)</td>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">&#x02212;127.4</td>
<td valign="top" align="center">1.46</td>
<td valign="top" align="center">0.19</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.039 (0.005)</td>
<td valign="top" align="center">&#x02212;0.007 (0.008)</td>
<td/>
<td/>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">&#x02212;126.9</td>
<td valign="top" align="center">1.98</td>
<td valign="top" align="center">0.14</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.036 (0.010)</td>
<td/>
<td valign="top" align="center">0.001 (0.011)</td>
<td/>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">&#x02212;126.2</td>
<td valign="top" align="center">2.73</td>
<td valign="top" align="center">0.10</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.115 (0.059)</td>
<td valign="top" align="center">&#x02212;0.008 (0.008)</td>
<td/>
<td valign="top" align="center">&#x02212;0.007 (0.005)</td>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">&#x02212;125.5</td>
<td valign="top" align="center">3.36</td>
<td valign="top" align="center">0.07</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.106 (0.060)</td>
<td/>
<td valign="top" align="center">0.004 (0.011)</td>
<td valign="top" align="center">&#x02212;0.007 (0.005)</td>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">&#x02212;124.6</td>
<td valign="top" align="center">4.33</td>
<td valign="top" align="center">0.05</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td/>
<td valign="top" align="center">0.040 (0.011)</td>
<td valign="top" align="center">&#x02212;0.001 (0.011)</td>
<td valign="top" align="center">&#x02212;0.001 (0.008)</td>
<td/>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">&#x02212;123.9</td>
<td valign="top" align="center">4.96</td>
<td valign="top" align="center">0.03</td>
</tr>
<tr>
<td valign="top" align="left">LY86 gene expression (<italic>n</italic> &#x0003D; 32; 13 nest-boxes)</td>
<td valign="top" align="center">0.418 (0.039)</td>
<td/>
<td/>
<td/>
<td/>
<td valign="top" align="center">3</td>
<td valign="top" align="center">1.7</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.41</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.453 (0.050)</td>
<td valign="top" align="center">&#x02212;0.086 (0.079)</td>
<td/>
<td/>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">3.1</td>
<td valign="top" align="center">1.45</td>
<td valign="top" align="center">0.20</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.392 (0.091)</td>
<td/>
<td valign="top" align="center">0.032 (0.101)</td>
<td/>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">4.2</td>
<td valign="top" align="center">2.52</td>
<td valign="top" align="center">0.12</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.357 (0.597)</td>
<td/>
<td/>
<td valign="top" align="center">0.005 (0.053)</td>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">4.3</td>
<td valign="top" align="center">2.61</td>
<td valign="top" align="center">0.11</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.434 (0.098)</td>
<td valign="top" align="center">&#x02212;0.084 (0.079)</td>
<td valign="top" align="center">0.023 (0.010)</td>
<td/>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">5.9</td>
<td valign="top" align="center">4.23</td>
<td valign="top" align="center">0.05</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.504 (0.601)</td>
<td valign="top" align="center">&#x02212;0.087 (0.080)</td>
<td/>
<td valign="top" align="center">&#x02212;0.005 (0.053)</td>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">5.9</td>
<td valign="top" align="center">4.27</td>
<td valign="top" align="center">0.05</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.389 (0.605)</td>
<td/>
<td valign="top" align="center">0.032 (0.106)</td>
<td valign="top" align="center">0.0002 (0.056)</td>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">7.0</td>
<td valign="top" align="center">5.35</td>
<td valign="top" align="center">0.03</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td/>
<td valign="top" align="center">1.183 (0.761)</td>
<td valign="top" align="center">&#x02212;1.710 (1.168)</td>
<td/>
<td valign="top" align="center">&#x02212;0.065 (0.067)</td>
<td valign="top" align="center">0.146 (0.105)</td>
<td valign="top" align="center">6</td>
<td valign="top" align="center">7.1</td>
<td valign="top" align="center">5.44</td>
<td valign="top" align="center">0.03</td>
</tr>
<tr>
<td valign="top" align="left">GATA3 gene expression (<italic>n</italic> &#x0003D; 32, 14 nest-boxes)</td>
<td valign="top" align="center">0.548 (0.228)</td>
<td/>
<td/>
<td valign="top" align="center">&#x02212;0.037 (0.020)</td>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">&#x02212;57.6</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.21</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.109 (0.020)</td>
<td valign="top" align="center">0.059 (0.032)</td>
<td/>
<td/>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">&#x02212;57.6</td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.21</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.131 (0.017)</td>
<td/>
<td/>
<td/>
<td/>
<td valign="top" align="center">3</td>
<td valign="top" align="center">&#x02212;57.2</td>
<td valign="top" align="center">0.37</td>
<td valign="top" align="center">0.18</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.458 (0.229)</td>
<td valign="top" align="center">0.047 (0.031)</td>
<td/>
<td valign="top" align="center">&#x02212;0.030 (0.020)</td>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">&#x02212;57.0</td>
<td valign="top" align="center">0.65</td>
<td valign="top" align="center">0.15</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.574 (0.230)</td>
<td/>
<td valign="top" align="center">0.021 (0.044)</td>
<td valign="top" align="center">&#x02212;0.041 (0.021)</td>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">&#x02212;55.0</td>
<td valign="top" align="center">2.59</td>
<td valign="top" align="center">0.06</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.113 (0.040)</td>
<td valign="top" align="center">0.059 (0.032)</td>
<td valign="top" align="center">&#x02212;0.005 (0.041)</td>
<td/>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">&#x02212;54.8</td>
<td valign="top" align="center">2.85</td>
<td valign="top" align="center">0.05</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.137 (0.040)</td>
<td/>
<td valign="top" align="center">&#x02212;0.008 (0.044)</td>
<td/>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">&#x02212;54.6</td>
<td valign="top" align="center">2.96</td>
<td valign="top" align="center">0.05</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.343 (0.292)</td>
<td valign="top" align="center">0.333 (0.454)</td>
<td/>
<td valign="top" align="center">&#x02212;0.020 (0.025)</td>
<td valign="top" align="center">&#x02212;0.025 (0.040)</td>
<td valign="top" align="center">6</td>
<td valign="top" align="center">&#x02212;54.3</td>
<td valign="top" align="center">3.31</td>
<td valign="top" align="center">0.04</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">0.478 (0.232)</td>
<td valign="top" align="center">0.046 (0.031)</td>
<td valign="top" align="center">0.019 (0.042)</td>
<td valign="top" align="center">&#x02212;0.034 (0.021)</td>
<td/>
<td valign="top" align="center">6</td>
<td valign="top" align="center">&#x02212;54.1</td>
<td valign="top" align="center">3.51</td>
<td valign="top" align="center">0.04</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>Models within a &#x00394;ICc value of six are presented along with their coefficient estimates and standard errors [estimate (SE)]. k &#x0003D; number of model parameters; w &#x0003D; Akaike&#x00027;s weight, defined as the probability of a model given the data and the candidate set of alternative models</italic>.</p>
<fn id="TN6">
<label>a</label>
<p><italic>&#x0201C;Original habitat &#x02013; forest&#x0201D; set as reference level and, thus, equals zero</italic>.</p></fn>
<fn id="TN7">
<label>b</label>
<p><italic>Malaria parasite &#x0003D; &#x0201C;NO&#x0201D; fixed as reference level; therefore, the given coefficients represent the change in gene transcript levels associated with the presence of malaria parasites. The statistical importance of weight and habitat in the top two models for GATA3 was further assessed by a likelihood-ratio test comparing such models against the intercept-only one. In both cases, these tests yielded a p-value of 0.08</italic>.</p></fn>
</table-wrap-foot>
</table-wrap>
<p>Prevalence of avian malaria in forest-reared nestlings was 79.16% (<italic>n</italic> &#x0003D; 48, <italic>SE</italic> &#x0003D; 5.92%), whereas every screened urban-reared nestling was found positive (<italic>n</italic> &#x0003D; 4). Within forest-reared birds, malaria prevalence did not vary based on nestling origin (Fisher&#x00027;s exact test, <italic>p</italic> &#x0003D; 0.468, <italic>n</italic> &#x0003D; 38. Prevalence &#x000B1; standard error: forest-originated birds &#x0003D; 0.83 &#x000B1; 0.38, urban-originated birds &#x0003D; 0.72 &#x000B1; 0.46).</p>
</sec>
<sec>
<title>Breeding success in urban and forest cross-fostered nests</title>
<p>Most metrics of breeding investment and success were lower in the urban compared to the forest site (Figure <xref ref-type="fig" rid="F3">3</xref>). Despite observed differences in clutch size between habitats (mean clutch size &#x000B1; SE: urban clutches &#x0003D; 7.70 &#x000B1; 0.34, <italic>n</italic> &#x0003D; 20; forest clutches &#x0003D; 8.63 &#x000B1; 0.48, <italic>n</italic> &#x0003D; 22), they did not received strong statistical support (Table <xref ref-type="table" rid="T4">4</xref>). Urban-originated eggs were slightly lighter than forest-originated eggs and habitat of origin appeared in three out of four models with &#x00394;AICc &#x0003C;2 for egg weight (the most supported one amongst them); however, the effect size was small, with urban eggs only 0.045 g lighter than forest ones. Experimental group was also kept in the set of top models showing that cross-fostered eggs were 0.05 grams heavier than control eggs (Figure <xref ref-type="fig" rid="F3">3A</xref>). No predictor was retained in the most supported model for hatching success (Table <xref ref-type="table" rid="T5">5</xref>, Figure <xref ref-type="fig" rid="F3">3B</xref>) and habitat of origin only appeared in a model featuring a &#x00394;AICc value of 1.89 (ranked third in support). These results suggested no differences in hatching success due to habitat of origin or rearing habitat (Table <xref ref-type="table" rid="T5">5</xref>). The urban rearing environment, however, had a very strong negative effect on fledging success as well as on nestling weight (Table <xref ref-type="table" rid="T5">5</xref>). Regardless of their origin, forest-reared nestlings were on average 1.33 grams heavier and more than twice as likely to fledge as urban-reared birds (Figures <xref ref-type="fig" rid="F3">3C,D</xref>). Rearing habitat appeared in every model within the &#x00394;2 set for fledging success and nestling weight (Table <xref ref-type="table" rid="T5">5</xref>).</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Breeding investment and breeding success for urban and forest blue tits. <bold>(A)</bold> Egg size across original habitats (urban and forest) and experimental groups. Cross-fostered eggs were swapped between sites and were reared in the fostering habitat; control eggs were swapped between nests within site. In the urban rearing habitat, forest- and urban-originated eggs differed in weight, whereas in the forest rearing habitat eggs in both experimental groups had similar weights. <bold>(B)</bold> Hatching success, <bold>(C)</bold> fledging success, and <bold>(D)</bold> weight of 13-day old nestlings across original and rearing habitats. Black dots illustrate raw data mean values and black bars illustrate &#x000B1; 2 standard errors. Raw data points are represented as partially transparent dots (see legend). Sample sizes are given beside mean values.</p></caption>
<graphic xlink:href="fevo-05-00064-g0003.tif"/>
</fig>
<table-wrap position="float" id="T4">
<label>Table 4</label>
<caption><p>Summarizing table of statistical models employed to explain variation in natural clutch size and egg weight (in g).</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th/>
<th valign="top" align="center" colspan="6" style="border-bottom: thin solid #000000;"><bold>Estimates of model coefficients (standard error)</bold></th>
<th/>
<th/>
<th/>
<th/>
</tr>
<tr>
<th valign="top" align="left"><bold>Response term</bold></th>
<th valign="top" align="center"><bold>Intercept</bold></th>
<th valign="top" align="center"><bold>Habitat &#x02013; urban<xref ref-type="table-fn" rid="TN8"><sup>a</sup></xref></bold></th>
<th valign="top" align="center"><bold>Experimental Group<xref ref-type="table-fn" rid="TN8"><sup>a</sup></xref> &#x02013; Cross-fostered</bold></th>
<th valign="top" align="center"><bold>Completion Date &#x02013; linear</bold></th>
<th valign="top" align="center"><bold>Habitat &#x000D7; experimental group<xref ref-type="table-fn" rid="TN8"><sup>a</sup></xref></bold></th>
<th valign="top" align="center"><bold>Clutch size</bold></th>
<th valign="top" align="center"><bold><italic>k</italic></bold></th>
<th valign="top" align="center"><bold>AICc</bold></th>
<th valign="top" align="center"><bold>&#x00394;AICc</bold></th>
<th valign="top" align="center"><bold><italic>w</italic></bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Clutch size (<italic>n</italic> &#x0003D; 42)</td>
<td valign="top" align="center">2.103 (0.054)</td>
<td/>
<td/>
<td/>
<td/>
<td valign="top" align="center">Not included</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">186.0</td>
<td valign="top" align="center">0</td>
<td valign="top" align="center">0.29</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">2.156 (0.073)</td>
<td valign="top" align="center">&#x02212;0.115 (0.108)</td>
<td/>
<td/>
<td/>
<td valign="top" align="center">Not included</td>
<td valign="top" align="center">2</td>
<td valign="top" align="center">187.1</td>
<td valign="top" align="center">1.08</td>
<td valign="top" align="center">0.17</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td/>
<td valign="top" align="center">1.785 (0.406)</td>
<td/>
<td/>
<td valign="top" align="center">0.008 (0.011)</td>
<td/>
<td valign="top" align="center">Not included</td>
<td valign="top" align="center">2</td>
<td valign="top" align="center">187.6</td>
<td valign="top" align="center">1.58</td>
<td valign="top" align="center">0.13</td>
</tr>
<tr>
<td valign="top" align="left">Egg weight (<italic>n</italic> &#x0003D; 248)</td>
<td valign="top" align="center">1.185 (0.020)</td>
<td valign="top" align="center">&#x02212;0.045 (0.024)</td>
<td valign="top" align="center">0.050 (0.024)<xref ref-type="table-fn" rid="TN9"><sup>&#x0002A;</sup></xref></td>
<td/>
<td/>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">&#x02212;632.3</td>
<td valign="top" align="center">0</td>
<td valign="top" align="center">0.17</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">1.164 (0.017)</td>
<td/>
<td valign="top" align="center">0.048 (0.025)</td>
<td/>
<td/>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">&#x02212;631.0</td>
<td valign="top" align="center">1.33</td>
<td valign="top" align="center">0.09</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">1.151 (0.054)</td>
<td valign="top" align="center">&#x02212;0.043 (0.024)</td>
<td valign="top" align="center">0.048 (0.024)</td>
<td/>
<td/>
<td valign="top" align="center">0.004 (0.006)</td>
<td valign="top" align="center">6</td>
<td valign="top" align="center">&#x02212;630.7</td>
<td valign="top" align="center">1.66</td>
<td valign="top" align="center">0.07</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">1.178 (0.022)</td>
<td valign="top" align="center">&#x02212;0.031 (0.033)</td>
<td valign="top" align="center">0.064 (0.033)</td>
<td/>
<td valign="top" align="center">&#x02212;0.029 (0.047)</td>
<td/>
<td valign="top" align="center">6</td>
<td valign="top" align="center">&#x02212;630.6</td>
<td valign="top" align="center">1.73</td>
<td valign="top" align="center">0.07</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>Only models within &#x00394;AICc &#x0003C;2 are presented along with their estimates and standard errors [estimate (SE)]. k &#x0003D; number of model parameters; w &#x0003D; Akaike&#x00027;s weight, defined as the probability of a model given the data and the candidate set of alternative models</italic>.</p>
<fn id="TN8">
<label>a</label>
<p><italic>&#x0201C;Forest&#x0201D; and &#x0201C;Control&#x0201D; categories are set as reference levels for &#x0201C;Habitat&#x0201D; and &#x0201C;Experimental Group,&#x0201D; respectively, and, therefore fixed to zero. Correspondingly, estimates for &#x0201C;Habitat &#x000D7; Group&#x0201D; refer to the urban habitat and cross-fostered experimental group. The additional term &#x0201C;Completion Date&#x02014;quadratic&#x0201D; was included in candidate models but estimates for this variable are not shown as they do not appear in any model within a &#x00394;AICc value of 2. Raw data mean values &#x000B1; SE for clutch size: urban clutches &#x0003D; 7.70 &#x000B1; 0.34; forest clutches &#x0003D; 8.63 &#x000B1; 0.48. Raw data mean values &#x000B1; SE for egg weight: urban-originated eggs &#x0003D; 1.17 &#x000B1; 0.01; forest-originated eggs &#x0003D; 1.21 &#x000B1; 0.01. Model coefficients for clutch size are shown in the scale of the link function (log). The statistical importance of variables in the top model for egg weight was further assessed by a likelihood-ratio test, dropping one predictor at a time, and</italic></p></fn>
<fn id="TN9">
<label>&#x0002A;</label>
<p><italic>illustrates p-values &#x0003C; 0.05</italic>.</p></fn>
</table-wrap-foot>
</table-wrap>
<table-wrap position="float" id="T5">
<label>Table 5</label>
<caption><p>Summarizing table of statistical models employed to explain variation in hatching success, fledging success and nestling weight (in g).</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th/>
<th valign="top" align="center" colspan="5" style="border-bottom: thin solid #000000;"><bold>Estimates of model coefficients (standard error)</bold></th>
<th/>
<th/>
<th/>
<th/>
</tr>
<tr>
<th valign="top" align="left"><bold>Response term</bold></th>
<th valign="top" align="center"><bold>Intercept</bold></th>
<th valign="top" align="center"><bold>Original Habitat &#x02013; urban<xref ref-type="table-fn" rid="TN10"><sup>a</sup></xref></bold></th>
<th valign="top" align="center"><bold>Rearing Habitat &#x02013; urban<xref ref-type="table-fn" rid="TN10"><sup>a</sup></xref></bold></th>
<th valign="top" align="center"><bold>Completion &#x02013; linear</bold></th>
<th valign="top" align="center"><bold>Completion &#x02013; quadratic</bold></th>
<th valign="top" align="center"><bold><italic>k</italic></bold></th>
<th valign="top" align="center"><bold>AICc</bold></th>
<th valign="top" align="center"><bold>&#x00394;AICc</bold></th>
<th valign="top" align="center"><bold><italic>w</italic></bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Hatching Success (<italic>n</italic> &#x0003D; 42)</td>
<td valign="top" align="center">&#x02212;0.126 (0.422)</td>
<td/>
<td/>
<td/>
<td/>
<td valign="top" align="center">2</td>
<td valign="top" align="center">166.4</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.33</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">2.187 (3.142)</td>
<td/>
<td/>
<td valign="top" align="center">&#x02212;0.060 (0.082)</td>
<td/>
<td valign="top" align="center">3</td>
<td valign="top" align="center">168.2</td>
<td valign="top" align="center">1.80</td>
<td valign="top" align="center">0.13</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td/>
<td valign="top" align="center">0.146 (0.578)</td>
<td valign="top" align="center">&#x02212;0.565 (0.838)</td>
<td/>
<td/>
<td/>
<td valign="top" align="center">3</td>
<td valign="top" align="center">168.3</td>
<td valign="top" align="center">1.89</td>
<td valign="top" align="center">0.13</td>
</tr>
<tr>
<td valign="top" align="left">Fledgling Success (<italic>n</italic> &#x0003D; 30)</td>
<td valign="top" align="center">4.652 (1.876)</td>
<td/>
<td valign="top" align="center">&#x02212;7.046 (3.082)<xref ref-type="table-fn" rid="TN11"><sup>&#x0002A;</sup></xref></td>
<td/>
<td/>
<td valign="top" align="center">3</td>
<td valign="top" align="center">74.3</td>
<td valign="top" align="center">0</td>
<td valign="top" align="center">0.35</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">&#x02212;9.252 (11.493)</td>
<td/>
<td valign="top" align="center">&#x02212;6.503 (2.935)</td>
<td valign="top" align="center">0.246 (0.213)</td>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">75.3</td>
<td valign="top" align="center">1.07</td>
<td valign="top" align="center">0.20</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td/>
<td valign="top" align="center">5.112 (2.003)</td>
<td valign="top" align="center">&#x02212;1.335 (1.615)</td>
<td valign="top" align="center">&#x02212;6.780 (2.828)</td>
<td/>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">76.2</td>
<td valign="top" align="center">1.95</td>
<td valign="top" align="center">0.13</td>
</tr>
<tr>
<td valign="top" align="left">Nestling Weight (<italic>n</italic> &#x0003D; 68)</td>
<td valign="top" align="center">11.190 (0.180)</td>
<td/>
<td valign="top" align="center">&#x02212;1.328 (0.341) <xref ref-type="table-fn" rid="TN11"><sup>&#x0002A;</sup></xref></td>
<td/>
<td/>
<td valign="top" align="center">4</td>
<td valign="top" align="center">178.9</td>
<td valign="top" align="center">0</td>
<td valign="top" align="center">0.27</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">8.822 (2.145)</td>
<td/>
<td valign="top" align="center">&#x02212;1.318 (0.330)</td>
<td valign="top" align="center">0.042 (0.038)</td>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">180.0</td>
<td valign="top" align="center">1.15</td>
<td valign="top" align="center">0.15</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">&#x02212;32.250 (26.282)</td>
<td/>
<td valign="top" align="center">&#x02212;1.075 (0.350)</td>
<td valign="top" align="center">1.526 (0.947)</td>
<td valign="top" align="center">&#x02212;0.013 (0.009)</td>
<td valign="top" align="center">6</td>
<td valign="top" align="center">180.1</td>
<td valign="top" align="center">1.22</td>
<td valign="top" align="center">0.15</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">11.300 (0.210)</td>
<td valign="top" align="center">&#x02212;0.278 (0.310)</td>
<td valign="top" align="center">&#x02212;1.322 (0.335)</td>
<td/>
<td/>
<td valign="top" align="center">5</td>
<td valign="top" align="center">180.4</td>
<td valign="top" align="center">1.54</td>
<td valign="top" align="center">0.12</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>Only models within &#x00394;AICc &#x0003C;2 are presented along with their estimates and standard errors (estimate (SE)). k &#x0003D; number of model parameters; w &#x0003D; Akaike&#x00027;s weight, defined as the probability of a model given the data and the candidate set of alternative models. The additional terms &#x0201C;Original habitat &#x000D7; Rearing habitat,&#x0201D; &#x0201C;Brood size&#x0201D; and &#x0201C;Rearing habitat &#x000D7; Brood size&#x0201D; were included in global models but estimates for these variables are not shown as they do not appear in any model within a &#x00394;AICc value of 2</italic>.</p>
<fn id="TN10">
<label>a</label>
<p><italic>&#x0201C;Original Habitat &#x02013; forest&#x0201D; and &#x0201C;Rearing habitat &#x02013; forest&#x0201D; were set as reference levels for coefficient estimation. Model coefficients for hatching and fledgling success are illustrated in link function scale (logit). The statistical importance of rearing habitat in the top model for fledgling success and nestling weight was further assessed by a likelihood-ratio test against the intercept-only model, and</italic></p></fn>
<fn id="TN11">
<label>&#x0002A;</label>
<p><italic>illustrates p-values &#x0003C; 0.05</italic>.</p></fn>
</table-wrap-foot>
</table-wrap>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>Reproductive success of blue tits was dramatically reduced in the city compared to the forest. Our study adds to existing evidence that urban environments commonly impair fitness and health in wild populations (Chamberlain et al., <xref ref-type="bibr" rid="B17">2009</xref>; Murray et al., <xref ref-type="bibr" rid="B45">2015</xref>; Bailly et al., <xref ref-type="bibr" rid="B7">2016b</xref>). In most avian studies, fledging success and nestling weight were lower in urban than in more rural habitats (Chamberlain et al., <xref ref-type="bibr" rid="B17">2009</xref>; Bailly et al., <xref ref-type="bibr" rid="B7">2016b</xref>; Sprau et al., <xref ref-type="bibr" rid="B62">2017</xref>). These findings are matched by observations in this study and in the parallel study of effects of food availability on other nests at our sites (Pollock et al., <xref ref-type="bibr" rid="B56">2017</xref>). In addition to immediate effects of the environment on developing offspring, such differences could also reflect prenatal maternal investment (e.g., differences in egg composition, Toledo et al., <xref ref-type="bibr" rid="B65">2016</xref>) or genetic differences between urban and rural populations (Mueller et al., <xref ref-type="bibr" rid="B43">2013</xref>). Our cross-fostering experiment does not support these latter hypotheses as we found no strong effects of natal origin on reproductive and fitness traits (Figure <xref ref-type="fig" rid="F3">3</xref>). Slightly larger cross-fostered eggs from both sites were a spurious outcome of our alternating, experimental swaps. Overall, we provide experimental support for the existence of a significant negative postnatal effect of the urban environmental on breeding success in blue tits, although we acknowledge that further replication will be needed to consolidate these results.</p>
<p>These findings converge with a recent study on closely related great tits which also used a cross-fostering experiment and showed that negative implications of urban nesting (in this case, shortened telomeres) arose from the raising environment and not from genetic or maternal factors (Salm&#x000F3;n et al., <xref ref-type="bibr" rid="B58">2016</xref>); however, this study did not investigate the possible mechanistic causes of such differences. Our parallel study on food availability and provisioning behavior provided clear evidence of resource limitation and nutritional stress in our urban site compared to our forest location (Pollock et al., <xref ref-type="bibr" rid="B56">2017</xref>). Whereas nestlings in the forest received predominantly caterpillars, the proportion of this preferred, high-quality diet was substantially lower for nestlings in the city. Pollock&#x00027;s data further suggest that parents partly resorted to other food, possibly including anthropogenic sources, to feed their nestlings and presumably themselves (Pollock et al., <xref ref-type="bibr" rid="B56">2017</xref>). Shifts in diet, which were also reported for urban populations in other species (e.g., Murray et al., <xref ref-type="bibr" rid="B45">2015</xref>), can have important effects on the nutritional state of urban animals, and consequently, also on their health (Andersson et al., <xref ref-type="bibr" rid="B3">2015</xref>). Thus, the findings that we present here on immunity have to be interpreted against the backdrop of resource restriction in the urban site.</p>
<p>Our results from adult blue tits on transcript levels of immune genes support some of our initial hypotheses. Urban blue tits showed increased transcript levels of TLR4, a marker of inflammatory processes, although these findings were not paralleled for <italic>LY86</italic>. Several factors have been proposed to explain how the immune system of urban organisms is expected to change in response to the environment (Isaksson, <xref ref-type="bibr" rid="B33">2015</xref>). Malaria parasite infection has been shown to affect the expression of <italic>TLR4</italic> (Martin et al., <xref ref-type="bibr" rid="B38">2014</xref>) and also <italic>LY86</italic> in passerines birds (Videvall et al., <xref ref-type="bibr" rid="B66">2015</xref>). In our study, malaria infection status probably had little importance for <italic>TLR4</italic> and <italic>LY86</italic> levels across habitats. In contrast to other studies (Evans et al., <xref ref-type="bibr" rid="B23">2009</xref>), malaria prevalence was consistently high at both of our sites; however, our statistical power was small and the lack of differences across habitat needs to be considered cautiously. Other intra-cellular pathogens found to vary in prevalence between urban and rural sites (Giraudeau et al., <xref ref-type="bibr" rid="B27">2014</xref>) or in association with anthropogenic food provisioning (Becker et al., <xref ref-type="bibr" rid="B10">2015</xref>), may also be important determinants of <italic>TLR4</italic> and <italic>LY86</italic> levels. Because we lack information on the wider pathogen assembly in our study sites, we cannot discard differences in other inflammatory pathogens as a cause of the observed gene expression patterns. Additional environmental factors could also explain differences in immunity between urban and forest populations. Low-quality and restricted diets are known to promote pro-inflammatory immune processes (Blount et al., <xref ref-type="bibr" rid="B12">2003</xref>; McGraw et al., <xref ref-type="bibr" rid="B40">2006</xref>; Isaksson, <xref ref-type="bibr" rid="B33">2015</xref>; Nettle et al., <xref ref-type="bibr" rid="B48">2017</xref>). Given our parallel findings of significant diet differences between our study populations (Pollock et al., <xref ref-type="bibr" rid="B56">2017</xref>), it is likely that resource limitation has contributed to elevated inflammation in the city (Larsson et al., <xref ref-type="bibr" rid="B35">2004</xref>; Isaksson, <xref ref-type="bibr" rid="B33">2015</xref>). Costs of inflammation under resource limitation could have contributed to the birds&#x00027; low reproductive success via exacerbated physiological trade-offs.</p>
<p>In contrast to <italic>TLR4</italic>, transcript levels of <italic>LY86</italic> in adult blue tits did not match our predictions. Because TLR4 and LY86 interact (Lee et al., <xref ref-type="bibr" rid="B36">2012</xref>), we expected to find a correlated pattern of expression between their gene expression. However, the interaction between these molecules is complex and can vary between cell types (Divanovic et al., <xref ref-type="bibr" rid="B19">2005</xref>). The lack of differences in <italic>LY86</italic> across habitats, in contrast to our findings regarding <italic>TLR4</italic>, indicates that these two genes might respond differentially to the urban environment. <italic>TLR4</italic> expression may be particularly sensitive to urban-related environmental stressors. As well as acting as a receptor for pathogens (gram negative bacteria molecular patterns), TLR4 is implicated in the recognition of damage-associated molecular patterns that follow tissue damage or cellular apoptosis (Liu et al., <xref ref-type="bibr" rid="B37">2014</xref>). Conceivably, <italic>TLR4</italic> could have been additionally enhanced by urban-associated environmental factors that promote oxidative stress and tissue damage, for example, air pollution and artificial light at night (Isaksson, <xref ref-type="bibr" rid="B32">2010</xref>; Fonken et al., <xref ref-type="bibr" rid="B25">2013</xref>). Our findings on <italic>TLR4</italic> and <italic>LY86</italic> transcript levels can be compared to recent findings of the transcriptomic comparison between urban and rural populations of another parid, the great tit (Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>). In fully grown great tits, expression of <italic>TLR4</italic> (gene ID ENSTGUG00000003342) tended to be elevated in blood and liver also in the city compared to the forest (Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>, Supplementary Datasets 1 and 2). Expression of <italic>LY86</italic> (gene ID ENSTGUG00000002305) also tended to be higher in liver of urban great tits, but not in blood. Largely, therefore, patterns were similar in the two parid species.</p>
<p>The findings of higher <italic>GATA3</italic> transcript levels in adult urban compared to forest blue tits also followed our initial predictions. The type 2 transcription factor GATA3 constitutes a major regulatory component of the immunity against helminth parasites (Tindemans et al., <xref ref-type="bibr" rid="B64">2014</xref>). Although we were unable to monitor intestinal parasite load in our study, our results could be explained by potentially higher helminth pressure in our urban location. Aggregation in feeding stations causes high rates of horizontal transmission of directly-transmitted parasites and, hence, affects parasite distribution (Becker et al., <xref ref-type="bibr" rid="B10">2015</xref>). If increased activation of GATA3-induced immunity in our urban birds was, indeed, caused by increased helminth infection, it might thus reflect a parasitological cost imposed by supplementary feeding and contribute to some of the negative effects of supplementary feeding on fitness (Plummer et al., <xref ref-type="bibr" rid="B55">2013</xref>). Results for <italic>GATA3</italic> expression (gene ID ENSTGUG00000002134) from an urban great tit study in Sweden were comparable for liver, with higher levels in the city, whereas for blood expression did not differ between sites (Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>).</p>
<p>For the nestling stage, our findings from the cross-fostering experiment provide no evidence that urban as compared to forest origin was associated with higher immune transcript levels in nestlings (Figure <xref ref-type="fig" rid="F2">2</xref>). Unfortunately, our sparse data from the urban rearing environment do not allow comprehensive conclusions from the full cross-fostering experiment (Figure <xref ref-type="supplementary-material" rid="SM1">S2</xref>). However, overall, transcript levels of nestling immune genes appeared to be similar across habitats (Figure <xref ref-type="supplementary-material" rid="SM1">S2</xref>). Thus, pending confirmation by a larger dataset, our findings suggest a differential impact of urbanization on the immune system, becoming more pronounced in later life-history stages. The contrasting patterns among age groups might be explained by chronic effects of long-term exposure to urban-specific environmental factors in adults (e.g., Bedrosian et al., <xref ref-type="bibr" rid="B11">2011</xref>). Genetic or prenatal maternal causes for increased <italic>TLR4</italic> and <italic>GATA3</italic> transcript levels in adult blue tits are made unlikely by our findings that among nestlings reared in the forest, origin (city or forest) did not strongly affect transcript levels of any gene (Figure <xref ref-type="fig" rid="F2">2</xref>). Hence, the cross-fostering experiment indicates that urban-forest differences in transcript levels of adult blue tits were most likely driven by the environment (Salm&#x000F3;n et al., <xref ref-type="bibr" rid="B58">2016</xref>). However, we cannot fully exclude the possibility that genetic or maternal factors are only identifiable in adulthood and not in nestlings. In contrast to our findings, a recent study on great tits did report differences in immune status between urban and rural nestlings (Bailly et al., <xref ref-type="bibr" rid="B6">2016a</xref>). Interestingly, the differences were opposite to those we found for adults blue tits: using physiological assays, rather than gene expression studies, Bailly et al. (<xref ref-type="bibr" rid="B6">2016a</xref>) show that urban nestlings produce less haptoglobin (a marker of inflammation) than forest birds, and offer as an explanation that haptoglobin production is compromised by food resource availability in the urban site. Although caution is needed when comparing different species and inflammatory markers, the opposite findings for the age groups might reflect different trade-offs between immunity and development: under resource restriction, growing nestlings might not be able to mount costly responses in the same way as fully grown adults (Sheldon and Verhulst, <xref ref-type="bibr" rid="B61">1996</xref>; Norris and Evans, <xref ref-type="bibr" rid="B49">2000</xref>; Alonso-Alvarez and Tella, <xref ref-type="bibr" rid="B2">2001</xref>).</p>
<p>Our results contribute to a body of evidence suggesting that urban living has reproductive costs and can impact the health of wild animals by altering their immune system (Audet et al., <xref ref-type="bibr" rid="B5">2016</xref>; Bailly et al., <xref ref-type="bibr" rid="B6">2016a</xref>; Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>). We acknowledge that caution is needed when interpreting our result given our sample size and the existence of only one urban-forest study pair and one study year. However, our experimental approach, with within-forest and within-city controls, allows us to draw conclusions on the causal links between immunity, fitness and the urban environment in our study system. Furthermore, the results generally confirm our original predictions and are in line with widespread evidence of reductions in fitness in urban environments (Chamberlain et al., <xref ref-type="bibr" rid="B17">2009</xref>; Bailly et al., <xref ref-type="bibr" rid="B7">2016b</xref>) and with recent discoveries on the effect of the urban environment on gene expression profiles in wild birds (Watson et al., <xref ref-type="bibr" rid="B69">2017</xref>). We find evidence for elevated expression of immune genes in adults of our urban compared to our forest population, and data from our cross-fostering experiment indicate that such changes are best explained by environmental factors. In combination with our twin study on reduced food availability and modified provisioning in the city (Pollock et al., <xref ref-type="bibr" rid="B56">2017</xref>), we tentatively identify links between resource limitation and altered immunity. To better understand the impact of urbanization on the immune system of wild animals, we need fine-scale characterization of the urban environment including parasite assemblies, coupled with a broader immune assessment of wild populations and information on fitness across life stages. Such information would reveal the role of the immune system for adaptation to urban life, as well as the long-term demographic consequences of altered immunity for urban-dwelling species.</p>
</sec>
<sec id="s5">
<title>Author contributions</title>
<p>PC-L, BH, DD, SB, and TB conceived the original study and designed experimental procedures. PC-L, DD, CP, and BH carried out field work, and PC-L, PO, MM, LW, and FB performed laboratory analysis. Statistical analyses were done by PC-L with advice from DD and BH. PC-L, BH, and DD wrote the manuscript with input from all other co-authors.</p>
<sec>
<title>Conflict of interest statement</title>
<p>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.</p>
</sec>
</sec>
</body>
<back>
<ack>
<p>We would like to thank Dan Haydon, Barbara Mable, and the IBAHCM for their support of the project. We would also like to thank Jessica Clark, Stephen Larcombe, Steve Duncan, Bernard Lundie, Paul Baker, Stewart White, Kim Mortega, Ruedi Nager, Iain Malzer, Paul Jerem, Ana Monteiro, and Yoana Ivanova for their help in the field and in the laboratory. We acknowledge the Trustees of the RSFS Forest Trust for access to Cashel Forest. Two reviewers provided helpful advice that greatly improved an earlier version of the manuscript.</p>
</ack><sec sec-type="supplementary-material" id="s6">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="http://journal.frontiersin.org/article/10.3389/fevo.2017.00064/full#supplementary-material">http://journal.frontiersin.org/article/10.3389/fevo.2017.00064/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Presentation1.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Alberti</surname> <given-names>M.</given-names></name></person-group> (<year>2015</year>). <article-title>Eco-evolutionary dynamics in an urbanizing planet</article-title>. <source>Trends Ecol. Evol.</source> <volume>30</volume>, <fpage>114</fpage>&#x02013;<lpage>126</lpage>. <pub-id pub-id-type="doi">10.1016/j.tree.2014.11.007</pub-id><pub-id pub-id-type="pmid">25498964</pub-id></citation></ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Alonso-Alvarez</surname> <given-names>C.</given-names></name> <name><surname>Tella</surname> <given-names>J. L.</given-names></name></person-group> (<year>2001</year>). <article-title>Effects of experimental food restriction and body-mass changes on the avian T-cell-mediated immune response</article-title>. <source>Can. J. Zool.</source> <volume>79</volume>, <fpage>101</fpage>&#x02013;<lpage>105</lpage>. <pub-id pub-id-type="doi">10.1139/z00-190</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Andersson</surname> <given-names>M. N.</given-names></name> <name><surname>Wang</surname> <given-names>H.-L.</given-names></name> <name><surname>Nord</surname> <given-names>A.</given-names></name> <name><surname>Salm&#x000F3;n</surname> <given-names>P.</given-names></name> <name><surname>Isaksson</surname> <given-names>C.</given-names></name></person-group> (<year>2015</year>). <article-title>Composition of physiologically important fatty acids in great tits differs between urban and rural populations on a seasonal basis</article-title>. <source>Front. Ecol. Evol.</source> <volume>3</volume>:<fpage>93</fpage>. <pub-id pub-id-type="doi">10.3389/fevo.2015.00093</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Atwell</surname> <given-names>J. W.</given-names></name> <name><surname>Cardoso</surname> <given-names>C.</given-names></name> <name><surname>Whittaker</surname> <given-names>D. J.</given-names></name> <name><surname>Price</surname> <given-names>T. D.</given-names></name> <name><surname>Ketterson</surname> <given-names>E. D.</given-names></name></person-group> (<year>2014</year>). <article-title>Correlated shifts in relation to population establishment in a novel environment</article-title>. <source>Am. Nat.</source> <volume>184</volume>, <fpage>E147</fpage>&#x02013;<lpage>E160</lpage>. <pub-id pub-id-type="doi">10.1086/678398</pub-id><pub-id pub-id-type="pmid">25438184</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Audet</surname> <given-names>J. N.</given-names></name> <name><surname>Ducatez</surname> <given-names>S.</given-names></name> <name><surname>Lefebvre</surname> <given-names>L.</given-names></name></person-group> (<year>2016</year>). <article-title>The town bird and the country bird: problem solving and immunocompetence vary with urbanization</article-title>. <source>Behav. Ecol.</source> <volume>27</volume>, <fpage>637</fpage>&#x02013;<lpage>644</lpage>. <pub-id pub-id-type="doi">10.1093/beheco/arv201</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bailly</surname> <given-names>J.</given-names></name> <name><surname>Scheifler</surname> <given-names>R.</given-names></name> <name><surname>Belvalette</surname> <given-names>M.</given-names></name> <name><surname>Garnier</surname> <given-names>S.</given-names></name> <name><surname>Boissier</surname> <given-names>E.</given-names></name> <name><surname>Anne</surname> <given-names>V.</given-names></name> <etal/></person-group>. (<year>2016a</year>). <article-title>Negative impact of urban habitat on immunity in the great tit Parus major</article-title>. <source>Oecologia</source> <volume>182</volume>, <fpage>1053</fpage>&#x02013;<lpage>1062</lpage>. <pub-id pub-id-type="doi">10.1007/s00442-016-3730-2</pub-id><pub-id pub-id-type="pmid">27646717</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bailly</surname> <given-names>J.</given-names></name> <name><surname>Scheifler</surname> <given-names>R.</given-names></name> <name><surname>Berthe</surname> <given-names>S.</given-names></name> <name><surname>Cl&#x000E9;ment-Demange</surname> <given-names>V. A.</given-names></name> <name><surname>Leblond</surname> <given-names>M.</given-names></name> <name><surname>Pasteur</surname> <given-names>B.</given-names></name> <etal/></person-group>. (<year>2016b</year>). <article-title>From eggs to fledging: Negative impact of urban habitat on reproduction in two tit species</article-title>. <source>J. Ornithol.</source> <volume>157</volume>, <fpage>377</fpage>&#x02013;<lpage>392</lpage>. <pub-id pub-id-type="doi">10.1007/s10336-015-1293-3</pub-id></citation></ref>
<ref id="B8">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Barton</surname> <given-names>K.</given-names></name></person-group> (<year>2016</year>). <source>MuMIn: Multi-Model Inference. R package version 1.15.6</source>. Available online at: <ext-link ext-link-type="uri" xlink:href="http://cran.r-project.org/package=MuMIn">http://cran.r-project.org/package=MuMIn</ext-link></citation></ref>
<ref id="B9">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Bates</surname> <given-names>D.</given-names></name> <name><surname>Maechler</surname> <given-names>M.</given-names></name> <name><surname>Bolker</surname> <given-names>B.</given-names></name> <name><surname>Walker</surname> <given-names>S.</given-names></name></person-group> (<year>2014</year>). <source>lme4: Linear Mixed-Effects Models using Eigen and S4. R Package Version 1.1-7</source>. Available online at: <ext-link ext-link-type="uri" xlink:href="http://cran.r-project.org/package=lme4">http://cran.r-project.org/package=lme4</ext-link></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Becker</surname> <given-names>D. J.</given-names></name> <name><surname>Streicker</surname> <given-names>D. G.</given-names></name> <name><surname>Altizer</surname> <given-names>S.</given-names></name></person-group> (<year>2015</year>). <article-title>Linking anthropogenic resources to wildlife-pathogen dynamics: a review and meta-analysis</article-title>. <source>Ecol. Lett.</source> <volume>18</volume>, <fpage>483</fpage>&#x02013;<lpage>495</lpage>. <pub-id pub-id-type="doi">10.1111/ele.12428</pub-id><pub-id pub-id-type="pmid">25808224</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bedrosian</surname> <given-names>T. A.</given-names></name> <name><surname>Fonken</surname> <given-names>L. K.</given-names></name> <name><surname>Walton</surname> <given-names>J. C.</given-names></name> <name><surname>Nelson</surname> <given-names>R. J.</given-names></name></person-group> (<year>2011</year>). <article-title>Chronic exposure to dim light at night suppresses immune responses in <italic>Siberian hamsters</italic></article-title>. <source>Biol. Lett.</source> <volume>7</volume>, <fpage>468</fpage>&#x02013;<lpage>471</lpage>. <pub-id pub-id-type="doi">10.1098/rsbl.2010.1108</pub-id><pub-id pub-id-type="pmid">21270021</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Blount</surname> <given-names>J. D.</given-names></name> <name><surname>Metcalfe</surname> <given-names>N. B.</given-names></name> <name><surname>Birkhead</surname> <given-names>T. R.</given-names></name> <name><surname>Surai</surname> <given-names>P. F.</given-names></name></person-group> (<year>2003</year>). <article-title>Carotenoid modulation of immune function and sexual attractiveness in zebra finches</article-title>. <source>Science</source> <volume>300</volume>, <fpage>125</fpage>&#x02013;<lpage>127</lpage>. <pub-id pub-id-type="doi">10.1126/science.1082142</pub-id><pub-id pub-id-type="pmid">12677066</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bonier</surname> <given-names>F.</given-names></name></person-group> (<year>2012</year>). <article-title>Hormones in the city: Endocrine ecology of urban birds</article-title>. <source>Horm. Behav.</source> <volume>61</volume>, <fpage>763</fpage>&#x02013;<lpage>772</lpage>. <pub-id pub-id-type="doi">10.1016/j.yhbeh.2012.03.016</pub-id><pub-id pub-id-type="pmid">22507448</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Burnham</surname> <given-names>K. P.</given-names></name></person-group> (<year>2004</year>). <article-title>Multimodel inference: understanding AIC and BIC in model selection</article-title>. <source>Sociol. Methods Res.</source> <volume>33</volume>, <fpage>261</fpage>&#x02013;<lpage>304</lpage>. <pub-id pub-id-type="doi">10.1177/0049124104268644</pub-id></citation></ref>
<ref id="B15">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Burnham</surname> <given-names>K. P.</given-names></name> <name><surname>Anderson</surname> <given-names>D. R.</given-names></name></person-group> (<year>2002</year>). <source>Model Selection and Multimodel Inference, 2nd Edn</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Springer-Verlag</publisher-name> <publisher-loc>New York</publisher-loc>.</citation></ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Burnham</surname> <given-names>K. P.</given-names></name> <name><surname>Anderson</surname> <given-names>D. R.</given-names></name> <name><surname>Huyvaert</surname> <given-names>K. P.</given-names></name></person-group> (<year>2011</year>). <article-title>AIC model selection and multimodel inference in behavioral ecology : some background, observations, and comparisons</article-title>. <source>Behav. Ecol. Sociobiol.</source> <volume>65</volume>, <fpage>23</fpage>&#x02013;<lpage>35</lpage>. <pub-id pub-id-type="doi">10.1007/s00265-010-1029-6</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chamberlain</surname> <given-names>D. E.</given-names></name> <name><surname>Cannon</surname> <given-names>A. R.</given-names></name> <name><surname>Toms</surname> <given-names>M. P.</given-names></name> <name><surname>Leech</surname> <given-names>D. I.</given-names></name></person-group> (<year>2009</year>). <article-title>Avian productivity in urban landscapes : a review and meta-analysis</article-title>. <source>Ibis</source> <volume>151</volume>, <fpage>1</fpage>&#x02013;<lpage>18</lpage>. <pub-id pub-id-type="doi">10.1111/j.1474-919X.2008.00899.x</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Davies</surname> <given-names>S.</given-names></name> <name><surname>Lane</surname> <given-names>S.</given-names></name> <name><surname>Meddle</surname> <given-names>S. L.</given-names></name> <name><surname>Tsutsui</surname> <given-names>K.</given-names></name> <name><surname>Deviche</surname> <given-names>P.</given-names></name></person-group> (<year>2016</year>). <article-title>The ecological and physiological bases of variation in the phenology of gonad growth in an urban and desert songbird</article-title>. <source>Gen. Comp. Endocrinol.</source> <volume>230</volume>, <fpage>17</fpage>&#x02013;<lpage>25</lpage>. <pub-id pub-id-type="doi">10.1016/j.ygcen.2016.03.013</pub-id><pub-id pub-id-type="pmid">26972152</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Divanovic</surname> <given-names>S.</given-names></name> <name><surname>Trompette</surname> <given-names>A.</given-names></name> <name><surname>Atabani</surname> <given-names>S. F.</given-names></name> <name><surname>Madan</surname> <given-names>R.</given-names></name> <name><surname>Golenbock</surname> <given-names>D. T.</given-names></name> <name><surname>Visintin</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2005</year>). <article-title>Negative regulation of Toll-like receptor 4 signaling by the Toll-like receptor homolog RP105</article-title>. <source>Nat. Immunol.</source> <volume>6</volume>, <fpage>571</fpage>&#x02013;<lpage>578</lpage>. <pub-id pub-id-type="doi">10.1038/ni1198</pub-id><pub-id pub-id-type="pmid">15852007</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dominoni</surname> <given-names>D. M.</given-names></name> <name><surname>Greif</surname> <given-names>S.</given-names></name> <name><surname>Nemeth</surname> <given-names>E.</given-names></name> <name><surname>Brumm</surname> <given-names>H.</given-names></name></person-group> (<year>2016</year>). <article-title>Airport noise predicts song timing of European birds</article-title>. <source>Ecol. Evol.</source> <volume>6</volume>, <fpage>6151</fpage>&#x02013;<lpage>6159</lpage>. <pub-id pub-id-type="doi">10.1002/ece3.2357</pub-id><pub-id pub-id-type="pmid">27648232</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dominoni</surname> <given-names>D. M.</given-names></name> <name><surname>Helm</surname> <given-names>B.</given-names></name> <name><surname>Lehmann</surname> <given-names>M.</given-names></name> <name><surname>Dowse</surname> <given-names>H. B.</given-names></name> <name><surname>Partecke</surname> <given-names>J.</given-names></name></person-group> (<year>2013a</year>). <article-title>Clocks for the city: circadian differences between forest and city songbirds</article-title>. <source>Proc. R. Soc. B</source> <volume>280</volume>:<fpage>20130593</fpage>. <pub-id pub-id-type="doi">10.1098/rspb.2013.0593</pub-id><pub-id pub-id-type="pmid">23740778</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dominoni</surname> <given-names>D.</given-names></name> <name><surname>Quetting</surname> <given-names>M.</given-names></name> <name><surname>Partecke</surname> <given-names>J.</given-names></name></person-group> (<year>2013b</year>). <article-title>Artificial light at night advances avian reproductive physiology</article-title>. <source>Proc. R. Soc. B</source> <volume>280</volume>:<fpage>20123017</fpage>. <pub-id pub-id-type="doi">10.1098/rspb.2012.3017</pub-id><pub-id pub-id-type="pmid">23407836</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Evans</surname> <given-names>K. L.</given-names></name> <name><surname>Gaston</surname> <given-names>K. J.</given-names></name> <name><surname>Sharp</surname> <given-names>S. P.</given-names></name> <name><surname>McGowan</surname> <given-names>A.</given-names></name> <name><surname>Simeoni</surname> <given-names>M.</given-names></name> <name><surname>Hatchwell</surname> <given-names>B. J.</given-names></name></person-group> (<year>2009</year>). <article-title>Effects of urbanisation on disease prevalence and age structure in blackbird Turdus merula populations</article-title>. <source>Oikos</source> <volume>118</volume>, <fpage>774</fpage>&#x02013;<lpage>782</lpage>. <pub-id pub-id-type="doi">10.1111/j.1600-0706.2008.17226.x</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fokidis</surname> <given-names>B. H.</given-names></name> <name><surname>Greiner</surname> <given-names>E. C.</given-names></name> <name><surname>Deviche</surname> <given-names>P.</given-names></name></person-group> (<year>2008</year>). <article-title>Interspecific variation in avian blood parasites and haematology associated with urbanization in a desert habitat</article-title>. <source>J. Avian Biol.</source> <volume>39</volume>, <fpage>300</fpage>&#x02013;<lpage>310</lpage>. <pub-id pub-id-type="doi">10.1111/j.0908-8857.2008.04248.x</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fonken</surname> <given-names>L. K.</given-names></name> <name><surname>Weil</surname> <given-names>Z. M.</given-names></name> <name><surname>Nelson</surname> <given-names>R. J.</given-names></name></person-group> (<year>2013</year>). <article-title>Mice exposed to dim light at night exaggerate inflammatory responses to lipopolysaccharide</article-title>. <source>Brain. Behav. Immun.</source> <volume>34</volume>, <fpage>159</fpage>&#x02013;<lpage>163</lpage>. <pub-id pub-id-type="doi">10.1016/j.bbi.2013.08.011</pub-id><pub-id pub-id-type="pmid">24012645</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gil</surname> <given-names>D.</given-names></name> <name><surname>Honarmand</surname> <given-names>M.</given-names></name> <name><surname>Pascual</surname> <given-names>J.</given-names></name> <name><surname>Perez-Mena</surname> <given-names>E.</given-names></name> <name><surname>Macias Garcia</surname> <given-names>C.</given-names></name></person-group> (<year>2014</year>). <article-title>Birds living near airports advance their dawn chorus and reduce overlap with aircraft noise</article-title>. <source>Behav. Ecol.</source> <volume>26</volume>, <fpage>435</fpage>&#x02013;<lpage>443</lpage>. <pub-id pub-id-type="doi">10.1093/beheco/aru207</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Giraudeau</surname> <given-names>M.</given-names></name> <name><surname>Mousel</surname> <given-names>M.</given-names></name> <name><surname>Earl</surname> <given-names>S.</given-names></name> <name><surname>McGraw</surname> <given-names>K.</given-names></name></person-group> (<year>2014</year>). <article-title>Parasites in the city: degree of urbanization predicts poxvirus and coccidian infections in house finches (<italic>Haemorhous mexicanus</italic>)</article-title>. <source>PLoS ONE</source> <volume>9</volume>:<fpage>e86747</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0086747</pub-id><pub-id pub-id-type="pmid">24503816</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Grimm</surname> <given-names>N. B.</given-names></name> <name><surname>Faeth</surname> <given-names>S. H.</given-names></name> <name><surname>Golubiewski</surname> <given-names>N. E.</given-names></name> <name><surname>Redman</surname> <given-names>C. L.</given-names></name> <name><surname>Wu</surname> <given-names>J.</given-names></name> <name><surname>Bai</surname> <given-names>X.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title>Global change and the ecology of cities</article-title>. <source>Science</source> <volume>319</volume>, <fpage>756</fpage>&#x02013;<lpage>760</lpage>. <pub-id pub-id-type="doi">10.1126/science.1150195</pub-id><pub-id pub-id-type="pmid">18258902</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Halliwell</surname> <given-names>B.</given-names></name> <name><surname>Gutteridge</surname> <given-names>J. M.</given-names></name></person-group> (<year>2002</year>). <source>Free Radicals in Biology and Medicine</source>. <publisher-loc>Oxford</publisher-loc>: <publisher-name>Oxford University Press</publisher-name>.</citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Harrison</surname> <given-names>X.</given-names></name></person-group> (<year>2015</year>). <article-title>A comparison of observation-level random effect and Beta-Binomial models for modelling overdispersion in Binomial data in ecology and evolution</article-title>. <source>PeerJ</source> <volume>3</volume>:<fpage>e1114</fpage>. <pub-id pub-id-type="doi">10.7717/peerj.1114</pub-id><pub-id pub-id-type="pmid">26244118</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hellgren</surname> <given-names>O.</given-names></name> <name><surname>Waldenstr&#x000F6;m</surname> <given-names>J.</given-names></name> <name><surname>Bensch</surname> <given-names>S.</given-names></name></person-group> (<year>2004</year>). <article-title>A new PCR assay for simultaneous studies of Leucocytozoon, Plasmodium, and Haemoproteus from avian blood</article-title>. <source>J. Parasitol.</source> <volume>90</volume>, <fpage>797</fpage>&#x02013;<lpage>802</lpage>. <pub-id pub-id-type="doi">10.1645/GE-184R1</pub-id><pub-id pub-id-type="pmid">15357072</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Isaksson</surname> <given-names>C.</given-names></name></person-group> (<year>2010</year>). <article-title>Pollution and its impact on wild animals: a meta-analysis on oxidative stress</article-title>. <source>Ecohealth</source> <volume>7</volume>, <fpage>342</fpage>&#x02013;<lpage>350</lpage>. <pub-id pub-id-type="doi">10.1007/s10393-010-0345-7</pub-id><pub-id pub-id-type="pmid">20865439</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Isaksson</surname> <given-names>C.</given-names></name></person-group> (<year>2015</year>). <article-title>Urbanization, oxidative stress and inflammation: a question of evolving, acclimatizing or coping with urban environmental stress</article-title>. <source>Funct. Ecol.</source> <volume>29</volume>, <fpage>913</fpage>&#x02013;<lpage>923</lpage>. <pub-id pub-id-type="doi">10.1111/1365-2435.12477</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lamb</surname> <given-names>C. T.</given-names></name> <name><surname>Mowat</surname> <given-names>G.</given-names></name> <name><surname>Mclellan</surname> <given-names>B. N.</given-names></name> <name><surname>Nielsen</surname> <given-names>S. E.</given-names></name> <name><surname>Boutin</surname> <given-names>S.</given-names></name></person-group> (<year>2017</year>). <article-title>Forbidden fruit : human settlement and abundant fruit create an ecological trap for an apex omnivore</article-title>. <source>J. Anim. Ecol.</source> <volume>86</volume>, <fpage>55</fpage>&#x02013;<lpage>65</lpage>. <pub-id pub-id-type="doi">10.1111/1365-2656.12589</pub-id><pub-id pub-id-type="pmid">27677529</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Larsson</surname> <given-names>S. C.</given-names></name> <name><surname>Kumlin</surname> <given-names>M.</given-names></name> <name><surname>Ingelman-Sundberg</surname> <given-names>M.</given-names></name> <name><surname>Wolk</surname> <given-names>A.</given-names></name></person-group> (<year>2004</year>). <article-title>Dietary long-chain n-3 fatty acids for the prevention of cancer: a review of potential mechanisms</article-title>. <source>Am. J. Clin. Nutr</source>. <volume>79</volume>, <fpage>935</fpage>&#x02013;<lpage>945</lpage>. <pub-id pub-id-type="pmid">15159222</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>C. C.</given-names></name> <name><surname>Avalos</surname> <given-names>A. M.</given-names></name> <name><surname>Ploegh</surname> <given-names>H. L.</given-names></name></person-group> (<year>2012</year>). <article-title>Accessory molecules for Toll-like receptors and their function</article-title>. <source>Nat. Rev. Immunol.</source> <volume>12</volume>, <fpage>168</fpage>&#x02013;<lpage>179</lpage>. <pub-id pub-id-type="doi">10.1038/nri3151</pub-id><pub-id pub-id-type="pmid">22301850</pub-id></citation></ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>Y.</given-names></name> <name><surname>Yin</surname> <given-names>H.</given-names></name> <name><surname>Zhao</surname> <given-names>M.</given-names></name> <name><surname>Lu</surname> <given-names>Q.</given-names></name></person-group> (<year>2014</year>). <article-title>TLR2 and TLR4 in autoimmune diseases: a comprehensive review</article-title>. <source>Clin. Rev. Allerg. Immunol.</source> <volume>47</volume>, <fpage>136</fpage>&#x02013;<lpage>147</lpage>. <pub-id pub-id-type="doi">10.1007/s12016-013-8402-y</pub-id><pub-id pub-id-type="pmid">24352680</pub-id></citation></ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Martin</surname> <given-names>L. B.</given-names></name> <name><surname>Coon</surname> <given-names>C. A. C.</given-names></name> <name><surname>Liebl</surname> <given-names>A. L.</given-names></name> <name><surname>Schrey</surname> <given-names>A. W.</given-names></name></person-group> (<year>2014</year>). <article-title>Surveillance for microbes and range expansion in house sparrows</article-title>. <source>Proc. R. Soc. B</source> <volume>281</volume>:<fpage>20132690</fpage>. <pub-id pub-id-type="doi">10.1098/rspb.2013.2690</pub-id><pub-id pub-id-type="pmid">24258722</pub-id></citation></ref>
<ref id="B39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Martin</surname> <given-names>L. B.</given-names></name> <name><surname>Hopkins</surname> <given-names>W. A.</given-names></name> <name><surname>Mydlarz</surname> <given-names>L. D.</given-names></name> <name><surname>Rohr</surname> <given-names>J. R.</given-names></name></person-group> (<year>2010</year>). <article-title>The effects of anthropogenic global changes on immune functions and disease resistance</article-title>. <source>Ann. N. Y. Acad. Sci.</source> <volume>1195</volume>, <fpage>129</fpage>&#x02013;<lpage>148</lpage>. <pub-id pub-id-type="doi">10.1111/j.1749-6632.2010.05454.x</pub-id><pub-id pub-id-type="pmid">20536821</pub-id></citation></ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>McGraw</surname> <given-names>K. J.</given-names></name> <name><surname>Crino</surname> <given-names>O. L.</given-names></name> <name><surname>Medina-Jerez</surname> <given-names>W.</given-names></name> <name><surname>Nolan</surname> <given-names>P. M.</given-names></name></person-group> (<year>2006</year>). <article-title>Effect of dietary carotenoid supplementation on food intake and immune function in a songbird with no carotenoid coloration</article-title>. <source>Ethology</source> <volume>112</volume>, <fpage>1209</fpage>&#x02013;<lpage>1216</lpage>. <pub-id pub-id-type="doi">10.1111/j.1439-0310.2006.01280.x</pub-id></citation></ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Medzhitov</surname> <given-names>R.</given-names></name></person-group> (<year>2001</year>). <article-title>Toll-like receptors and innate immunity</article-title>. <source>Nat. Rev. Immunol.</source> <volume>1</volume>, <fpage>135</fpage>&#x02013;<lpage>145</lpage>. <pub-id pub-id-type="doi">10.1038/35100529</pub-id><pub-id pub-id-type="pmid">11905821</pub-id></citation></ref>
<ref id="B42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mueller</surname> <given-names>J. C.</given-names></name> <name><surname>Kuhl</surname> <given-names>H.</given-names></name> <name><surname>Timmermann</surname> <given-names>B.</given-names></name> <name><surname>Kempenaers</surname> <given-names>B.</given-names></name></person-group> (<year>2016</year>). <article-title>Characterization of the genome and transcriptome of the blue tit <italic>Cyanistes caeruleus</italic> : polymorphisms, sex-biased expression and selection signals</article-title>. <source>Mol. Ecol. Resour.</source> <volume>16</volume>, <fpage>549</fpage>&#x02013;<lpage>561</lpage>. <pub-id pub-id-type="doi">10.1111/1755-0998.12450</pub-id><pub-id pub-id-type="pmid">26220359</pub-id></citation></ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mueller</surname> <given-names>J. C.</given-names></name> <name><surname>Partecke</surname> <given-names>J.</given-names></name> <name><surname>Hatchwell</surname> <given-names>B. J.</given-names></name> <name><surname>Gaston</surname> <given-names>K. J.</given-names></name> <name><surname>Evans</surname> <given-names>K. L.</given-names></name></person-group> (<year>2013</year>). <article-title>Candidate gene polymorphisms for behavioural adaptations during urbanization in blackbirds</article-title>. <source>Mol. Ecol.</source> <volume>22</volume>, <fpage>3629</fpage>&#x02013;<lpage>3637</lpage>. <pub-id pub-id-type="doi">10.1111/mec.12288</pub-id><pub-id pub-id-type="pmid">23495914</pub-id></citation></ref>
<ref id="B44">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Munshi-South</surname> <given-names>J.</given-names></name> <name><surname>Kharchenko</surname> <given-names>K.</given-names></name></person-group> (<year>2010</year>). <article-title>Rapid, pervasive genetic differentiation of urban white-footed mouse (<italic>Peromyscus leucopus</italic>) populations in New York City</article-title>. <source>Mol. Ecol.</source> <volume>19</volume>, <fpage>4242</fpage>&#x02013;<lpage>4254</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-294X.2010.04816.x</pub-id><pub-id pub-id-type="pmid">20819163</pub-id></citation></ref>
<ref id="B45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Murray</surname> <given-names>M.</given-names></name> <name><surname>Edwards</surname> <given-names>M. A.</given-names></name> <name><surname>Abercrombie</surname> <given-names>B.</given-names></name> <name><surname>Cassady</surname> <given-names>C.</given-names></name> <name><surname>Clair</surname> <given-names>S.</given-names></name></person-group> (<year>2015</year>). <article-title>Poor health is associated with use of anthropogenic resources in an urban carnivore</article-title>. <source>Proc. R. Soc. B</source> <volume>282</volume>:<fpage>20150009</fpage>. <pub-id pub-id-type="doi">10.1098/rspb.2015.0009</pub-id><pub-id pub-id-type="pmid">25876843</pub-id></citation></ref>
<ref id="B46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Neiderud</surname> <given-names>C.-J.</given-names></name></person-group> (<year>2015</year>). <article-title>How urbanization affects the epidemiology of emerging infectious diseases</article-title>. <source>Infect. Ecol. Epidemiol.</source> <volume>5</volume>:<fpage>27060</fpage>. <pub-id pub-id-type="doi">10.3402/iee.v5.27060</pub-id><pub-id pub-id-type="pmid">26112265</pub-id></citation></ref>
<ref id="B47">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nemeth</surname> <given-names>E.</given-names></name> <name><surname>Pieretti</surname> <given-names>N.</given-names></name> <name><surname>Zollinger</surname> <given-names>S. A.</given-names></name> <name><surname>Geberzahn</surname> <given-names>N.</given-names></name> <name><surname>Partecke</surname> <given-names>J.</given-names></name> <name><surname>Miranda</surname> <given-names>A. C.</given-names></name> <etal/></person-group>. (<year>2013</year>). <article-title>Bird song and anthropogenic noise: vocal constraints may explain why birds sing higher-frequency songs in cities</article-title>. <source>Proc. R. Soc. B</source> <volume>280</volume>:<fpage>20122798</fpage>. <pub-id pub-id-type="doi">10.1098/rspb.2012.2798</pub-id><pub-id pub-id-type="pmid">23303546</pub-id></citation></ref>
<ref id="B48">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nettle</surname> <given-names>D.</given-names></name> <name><surname>Andrews</surname> <given-names>C.</given-names></name> <name><surname>Reichert</surname> <given-names>S.</given-names></name> <name><surname>Bedford</surname> <given-names>T.</given-names></name> <name><surname>Kolenda</surname> <given-names>C.</given-names></name> <name><surname>Parker</surname> <given-names>C.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Early-life adversity accelerates cellular ageing and affects adult inflammation: experimental evidence from the European starling</article-title>. <source>Sci. Rep.</source> <volume>7</volume>:<fpage>40794</fpage>. <pub-id pub-id-type="doi">10.1038/srep40794</pub-id><pub-id pub-id-type="pmid">28094324</pub-id></citation></ref>
<ref id="B49">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Norris</surname> <given-names>K.</given-names></name> <name><surname>Evans</surname> <given-names>M. R.</given-names></name></person-group> (<year>2000</year>). <article-title>Ecological immunology: life history trade-offs and immune defense in birds</article-title>. <source>Behav. Ecol.</source> <volume>11</volume>, <fpage>19</fpage>&#x02013;<lpage>26</lpage>. <pub-id pub-id-type="doi">10.1093/beheco/11.1.19</pub-id></citation></ref>
<ref id="B50">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>O&#x00027;Shaughnessy</surname> <given-names>P. J.</given-names></name> <name><surname>Marsh</surname> <given-names>P.</given-names></name> <name><surname>Dudley</surname> <given-names>K.</given-names></name></person-group> (<year>1994</year>). <article-title>Follicle-stimulating hormone receptor mRNA in the mouse ovary during post-natal development in the normal mouse and in the adult hypogonadal (hpg) mouse: structure of alternate transcripts</article-title>. <source>Mol. Cell. Endocrinol.</source> <volume>101</volume>, <fpage>197</fpage>&#x02013;<lpage>201</lpage>. <pub-id pub-id-type="doi">10.1016/0303-7207(94)90235-6</pub-id><pub-id pub-id-type="pmid">9397953</pub-id></citation></ref>
<ref id="B51">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>O&#x00027;Shaughnessy</surname> <given-names>P. J.</given-names></name> <name><surname>Morris</surname> <given-names>I. D.</given-names></name> <name><surname>Baker</surname> <given-names>P. J.</given-names></name></person-group> (<year>2008</year>). <article-title>Leydig cell re-generation and expression of cell signaling molecules in the germ cell-free testis</article-title>. <source>Reproduction</source> <volume>135</volume>, <fpage>851</fpage>&#x02013;<lpage>858</lpage>. <pub-id pub-id-type="doi">10.1530/REP-07-0529</pub-id><pub-id pub-id-type="pmid">18502897</pub-id></citation></ref>
<ref id="B52">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>O&#x00027;Shaughnessy</surname> <given-names>P. J.</given-names></name> <name><surname>Murphy</surname> <given-names>L.</given-names></name></person-group> (<year>1993</year>). <article-title>Cytochrome P-450 17 -hydroxylase protein and mRNA in the testis of the testicular feminized (Tfm) mouse</article-title>. <source>J. Mol. Endocrinol.</source> <volume>11</volume>, <fpage>77</fpage>&#x02013;<lpage>82</lpage>. <pub-id pub-id-type="doi">10.1677/jme.0.0110077</pub-id><pub-id pub-id-type="pmid">8240674</pub-id></citation></ref>
<ref id="B53">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Olias</surname> <given-names>P.</given-names></name> <name><surname>Adam</surname> <given-names>I.</given-names></name> <name><surname>Meyer</surname> <given-names>A.</given-names></name> <name><surname>Scharff</surname> <given-names>C.</given-names></name> <name><surname>Gruber</surname> <given-names>A. D.</given-names></name></person-group> (<year>2014</year>). <article-title>Reference genes for quantitative gene expression studies in multiple avian species</article-title>. <source>PLoS ONE</source> <volume>9</volume>:<fpage>e99678</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0099678</pub-id><pub-id pub-id-type="pmid">24926893</pub-id></citation></ref>
<ref id="B54">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Partecke</surname> <given-names>J.</given-names></name> <name><surname>Schwabl</surname> <given-names>I.</given-names></name> <name><surname>Gwinner</surname> <given-names>E.</given-names></name></person-group> (<year>2006</year>). <article-title>Stress and the city: urbanization and its effects on the stress physiology in European blackbirds</article-title>. <source>Ecology</source> <volume>87</volume>, <fpage>1945</fpage>&#x02013;<lpage>1952</lpage>. <pub-id pub-id-type="doi">10.1890/0012-9658(2006)87[1945:SATCUA]2.0.CO;2</pub-id><pub-id pub-id-type="pmid">16937632</pub-id></citation></ref>
<ref id="B55">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Plummer</surname> <given-names>K. E.</given-names></name> <name><surname>Bearhop</surname> <given-names>S.</given-names></name> <name><surname>Leech</surname> <given-names>D. I.</given-names></name> <name><surname>Chamberlain</surname> <given-names>D. E.</given-names></name> <name><surname>Blount</surname> <given-names>J. D.</given-names></name></person-group> (<year>2013</year>). <article-title>Winter food provisioning reduces future breeding performance in a wild bird</article-title>. <source>Sci. Rep.</source> <volume>3</volume>:<fpage>2002</fpage>. <pub-id pub-id-type="doi">10.1038/srep02002</pub-id><pub-id pub-id-type="pmid">23788126</pub-id></citation></ref>
<ref id="B56">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pollock</surname> <given-names>C. J.</given-names></name> <name><surname>Capilla-Lasheras</surname> <given-names>P.</given-names></name> <name><surname>McGill</surname> <given-names>R. A. R.</given-names></name> <name><surname>Helm</surname> <given-names>B.</given-names></name> <name><surname>Dominoni</surname> <given-names>D.</given-names></name></person-group> (<year>2017</year>). <article-title>Integrated behavioural and stable isotope data reveal altered diet linked to low breeding success in urban-dwelling blue tits (<italic>Cyanistes caeruleus</italic>)</article-title>. <source>Sci. Rep</source>.</citation></ref>
<ref id="B57">
<citation citation-type="web"><person-group person-group-type="author"><collab>R Core Team</collab></person-group> (<year>2016</year>). <source>R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing</source>. <publisher-loc>Vienna</publisher-loc>. Available online at: <ext-link ext-link-type="uri" xlink:href="http://www.r-project.org/">http://www.r-project.org/</ext-link></citation></ref>
<ref id="B58">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Salm&#x000F3;n</surname> <given-names>P.</given-names></name> <name><surname>Nilsson</surname> <given-names>J. F.</given-names></name> <name><surname>Nord</surname> <given-names>A.</given-names></name> <name><surname>Bensch</surname> <given-names>S.</given-names></name> <name><surname>Isaksson</surname> <given-names>C.</given-names></name></person-group> (<year>2016</year>). <article-title>Urban environment shortens telomere length in nestling great tits, Parus major</article-title>. <source>Biol. Lett.</source> <volume>12</volume>:<fpage>20160155</fpage>. <pub-id pub-id-type="doi">10.1098/rsbl.2016.0155</pub-id><pub-id pub-id-type="pmid">27303051</pub-id></citation></ref>
<ref id="B59">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schoech</surname> <given-names>S. J.</given-names></name> <name><surname>Bowman</surname> <given-names>R.</given-names></name> <name><surname>Reynolds</surname> <given-names>S. J.</given-names></name></person-group> (<year>2004</year>). <article-title>Food supplementation and possible mechanisms underlying early breeding in the Florida Scrub-Jay (<italic>Aphelocoma coerulescens</italic>)</article-title>. <source>Horm. Behav.</source> <volume>46</volume>, <fpage>565</fpage>&#x02013;<lpage>573</lpage>. <pub-id pub-id-type="doi">10.1016/j.yhbeh.2004.06.005</pub-id><pub-id pub-id-type="pmid">15555498</pub-id></citation></ref>
<ref id="B60">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Seto</surname> <given-names>K. C.</given-names></name> <name><surname>G&#x000FC;neralp</surname> <given-names>B.</given-names></name> <name><surname>Hutyra</surname> <given-names>L. R.</given-names></name></person-group> (<year>2012</year>). <article-title>Global forecasts of urban expansion to 2030 and direct impacts on biodiversity and carbon pools</article-title>. <source>Proc. Natl. Acad. Sci. U.S.A.</source> <volume>109</volume>, <fpage>16083</fpage>&#x02013;<lpage>16088</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1211658109</pub-id><pub-id pub-id-type="pmid">22988086</pub-id></citation></ref>
<ref id="B61">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sheldon</surname> <given-names>B. C.</given-names></name> <name><surname>Verhulst</surname> <given-names>S.</given-names></name></person-group> (<year>1996</year>). <article-title>Ecological immunology: costly parasite defences and trade-offs in evolutionary ecology</article-title>. <source>Trends Ecol. Evol.</source> <volume>11</volume>, <fpage>317</fpage>&#x02013;<lpage>321</lpage>. <pub-id pub-id-type="doi">10.1016/0169-5347(96)10039-2</pub-id><pub-id pub-id-type="pmid">21237861</pub-id></citation></ref>
<ref id="B62">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sprau</surname> <given-names>P.</given-names></name> <name><surname>Mouchet</surname> <given-names>A.</given-names></name> <name><surname>Dingemanse</surname> <given-names>N. J.</given-names></name></person-group> (<year>2017</year>). <article-title>Multidimensional environmental predictors of variation in avian forest and city life histories</article-title>. <source>Behav. Ecol</source>. <volume>28</volume>, <fpage>59</fpage>&#x02013;<lpage>68</lpage>. <pub-id pub-id-type="doi">10.1093/beheco/arw130</pub-id></citation></ref>
<ref id="B63">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sumasgutner</surname> <given-names>P.</given-names></name> <name><surname>Nemeth</surname> <given-names>E.</given-names></name> <name><surname>Tebb</surname> <given-names>G.</given-names></name> <name><surname>Krenn</surname> <given-names>H. W.</given-names></name> <name><surname>Gamauf</surname> <given-names>A.</given-names></name></person-group> (<year>2014</year>). <article-title>Hard times in the city - attractive nest sites but insufficient food supply lead to low reproduction rates in a bird of prey</article-title>. <source>Front. Zool.</source> <volume>11</volume>:<fpage>48</fpage>. <pub-id pub-id-type="doi">10.1186/1742-9994-11-48</pub-id><pub-id pub-id-type="pmid">24872836</pub-id></citation></ref>
<ref id="B64">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tindemans</surname> <given-names>I.</given-names></name> <name><surname>Serafini</surname> <given-names>N.</given-names></name> <name><surname>DiSanto</surname> <given-names>J. P.</given-names></name> <name><surname>Hendriks</surname> <given-names>R. W.</given-names></name></person-group> (<year>2014</year>). <article-title>GATA-3 function in innate and adaptive immunity</article-title>. <source>Immunity</source> <volume>41</volume>, <fpage>191</fpage>&#x02013;<lpage>206</lpage>. <pub-id pub-id-type="doi">10.1016/j.immuni.2014.06.006</pub-id><pub-id pub-id-type="pmid">25148023</pub-id></citation></ref>
<ref id="B65">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Toledo</surname> <given-names>A.</given-names></name> <name><surname>Andersson</surname> <given-names>M. N.</given-names></name> <name><surname>Wang</surname> <given-names>H.-L.</given-names></name> <name><surname>Salm&#x000F3;n</surname> <given-names>P.</given-names></name> <name><surname>Watson</surname> <given-names>H.</given-names></name> <name><surname>Burdge</surname> <given-names>G. C.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Fatty acid profiles of great tit (<italic>Parus major</italic>) eggs differ between urban and rural habitats, but not between coniferous and deciduous forests</article-title>. <source>Sci. Nat</source>. <volume>103</volume>:<fpage>55</fpage>. <pub-id pub-id-type="doi">10.1007/s00114-016-1381-0</pub-id></citation></ref>
<ref id="B66">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Videvall</surname> <given-names>E.</given-names></name> <name><surname>Cornwallis</surname> <given-names>C. K.</given-names></name> <name><surname>Palinauskas</surname> <given-names>V.</given-names></name> <name><surname>Valki&#x0016B;nas</surname> <given-names>G.</given-names></name> <name><surname>Hellgren</surname> <given-names>O.</given-names></name></person-group> (<year>2015</year>). <article-title>The avian transcriptome response to malaria infection</article-title>. <source>Mol. Biol. Evol.</source> <volume>32</volume>, <fpage>1255</fpage>&#x02013;<lpage>1267</lpage>. <pub-id pub-id-type="doi">10.1093/molbev/msv016</pub-id><pub-id pub-id-type="pmid">25636457</pub-id></citation></ref>
<ref id="B67">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>Y.</given-names></name> <name><surname>Su</surname> <given-names>M. A.</given-names></name> <name><surname>Wan</surname> <given-names>Y. Y.</given-names></name></person-group> (<year>2011</year>). <article-title>an essential role of the transcription factor GATA-3 for the function of regulatory T cells</article-title>. <source>Immunity</source> <volume>35</volume>, <fpage>337</fpage>&#x02013;<lpage>348</lpage>. <pub-id pub-id-type="doi">10.1016/j.immuni.2011.08.012</pub-id><pub-id pub-id-type="pmid">21924928</pub-id></citation></ref>
<ref id="B68">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Warren</surname> <given-names>W. C.</given-names></name> <name><surname>Clayton</surname> <given-names>D. F.</given-names></name> <name><surname>Ellegren</surname> <given-names>H.</given-names></name> <name><surname>Arnold</surname> <given-names>A. P.</given-names></name> <name><surname>Hillier</surname> <given-names>L. W.</given-names></name> <name><surname>K&#x000FC;nstner</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2010</year>). <article-title>The genome of a songbird</article-title>. <source>Nature</source> <volume>464</volume>, <fpage>757</fpage>&#x02013;<lpage>762</lpage>. <pub-id pub-id-type="doi">10.1038/nature08819</pub-id><pub-id pub-id-type="pmid">20360741</pub-id></citation></ref>
<ref id="B69">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Watson</surname> <given-names>H.</given-names></name> <name><surname>Videvall</surname> <given-names>E.</given-names></name> <name><surname>Andersson</surname> <given-names>M. N.</given-names></name> <name><surname>Isaksson</surname> <given-names>C.</given-names></name></person-group> (<year>2017</year>). <article-title>Transcriptome analysis of a wild bird reveals physiological responses to the urban environment</article-title>. <source>Sci. Rep.</source> <volume>7</volume>:<fpage>44180</fpage>. <pub-id pub-id-type="doi">10.1038/srep44180</pub-id><pub-id pub-id-type="pmid">28290496</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn fn-type="financial-disclosure"><p><bold>Funding.</bold> PC-L was funded by a postgraduate scholarship from Iberdrola Foundation. Funding of DD and of gene analysis was provided by a Marie-Curie Career Integration Grant to BH [EC CIG (618578) Wildclocks].</p>
</fn>
</fn-group>
</back>
</article>
