# Serotonin transporter (SERT) polymorphisms, personality and problem-solving in urban great tits

Andrea S. Grunst, Melissa L. Grunst, Nicky Staes, Bert Thys, Rianne Pinxten, Marcel Eens

► To cite this version:

Andrea S. Grunst, Melissa L. Grunst, Nicky Staes, Bert Thys, Rianne Pinxten, et al. Serotonin transporter (SERT) polymorphisms, personality and problem-solving in urban great tits. Scientific Reports, 2021, 11, 10.1038/s41598-021-03466-7. insu-03708694

HAL Id: insu-03708694

https://insu.hal.science/insu-03708694

Submitted on 29 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Serotonin transporter (SERT) polymorphisms, personality and problem-solving in urban great titls

Andrea S. Grunst<sup>1,2,5✉</sup>, Melissa L. Grunst<sup>1,2,5</sup>, Nicky Staes<sup>1,3</sup>, Bert Thys<sup>1</sup>, Rianne Pinxten<sup>1,4</sup> & Marcel Eens<sup>1</sup>

Understanding underlying genetic variation can elucidate how diversity in behavioral phenotypes evolves and is maintained. Genes in the serotonergic signaling pathway, including the serotonin transporter gene (SERT), are candidates for affecting animal personality, cognition and fitness. In a model species, the great tit (Parus major), we reevaluated previous findings suggesting relationships between SERT polymorphisms, neophobia, exploratory behavior and fitness parameters, and performed a first test of the relationship between single nucleotide polymorphisms (SNPs) in SERT and problem-solving in birds. We found some evidence for associations between SERT SNPs and neophobia, exploratory behavior and laying date. Furthermore, several SNPs were associated with behavioral patterns and success rates during obstacle removal problem-solving tests performed at nest boxes. In females, minor allele homozygotes (AA) for nonsynonymous SNP226 in exon 1 made fewer incorrect attempts and were more likely to problem-solve. In both sexes, there was some evidence that minor allele homozygotes (CC) for SNP84 in exon 9 were more likely to problem-solve. Only one SNP-behavior relationship was statistically significant after correcting for multiple comparisons, but several were associated with substantial effect sizes. Our study provides a foundation for future research on the genetic basis of behavioral and cognitive variation in wild animal populations.

Animal personality is defined by consistent, among individual variation in behavioral traits such as neophobia, exploratory behavior and aggressiveness. Among individual variation in behavior was once discounted as noise around an adaptive mean. However, behavioral ecologists now recognize that personality variation describes an integral part of the variation in behavioral phenotypes, and is shaped by selective pressures<sup>1-6</sup>. Personality variation reflects alternative behavioral strategies that involve differential balance of life-history tradeoffs, and affects how individuals respond to ecological and social challenges, with ultimate fitness implications<sup>1-6</sup>. For example, individuals that are less neophobic, bolder and more exploratory are often more successful in invading new habitats, but can also be susceptible to risks arising from novel predators and poisoning<sup>7,8</sup>. In addition, personality variation may be related to variation in cognitive strategies and innovative problem-solving performance, wherein individuals solve novel problems, independent of overall cognitive ability<sup>9</sup>. For instance, less neophobic and more exploratory individuals may learn more rapidly and perform better on innovative problem-solving tasks, because they are more likely to attempt novel tasks and interact with novel environmental features<sup>9-11</sup>. Thus, cognitive performance traits and personality may be interrelated, and might be affected by shared physiological and genetic mechanisms.

Personality and cognitive traits are often heritable and show genetic variance. Quantitative genetic studies present estimates as high as 50% for the amount of personality variation attributed to genetics<sup>12-14</sup>. Quantitative genetics and artificial selection studies also suggest that cognitive traits, including general cognitive ability and cognitive performance within specific domains, can show substantial genetic variation<sup>15</sup>. However, other studies

<footnote n="1"><sup>1</sup>Department of Biology, Behavioural Ecology and Ecophysiology Group, University of Antwerp, Antwerp, Belgium. <sup>2</sup>Littoral Environnement Et Sociétés, La Rochelle Université, La Rochelle, France. <sup>3</sup>Centre for Research and Conservation, Royal Zoological Society of Antwerp, Antwerp, Belgium. <sup>4</sup>Faculty of Social Sciences, Antwerp School of Education, University of Antwerp, Antwerp, Belgium. <sup>5</sup>These authors contributed equally: Andrea S. Grunst and Melissa L. Grunst. <sup>6</sup>email: andrea.grunst@uantwerpen.be</footnote>
have yielded low estimates for the heritability of cognitive performance traits, including innovative problem-solving¹⁶. Developmental, environmental and genetic factors combine to determine behavioral and cognitive traits, and multiple, interacting genes affect these complex phenotypes, with most having small effects¹⁴,¹⁵. The extent to which genetics underlie personality and cognitive variation is therefore challenging to elucidate, and specific genes underlying this variation are challenging to identify. For instance, personality traits, including boldness and aggressiveness¹⁷⁻²⁰, and cognitive traits including problem-solving performance²¹, often differ between urban and rural bird populations. However, whether these differences are genetic in origin, and which genes might be involved, remains debated²²,²³. Nonetheless, there are a number of exciting cases in which specific genes underlying personality and cognitive variation have been identified, via approaches including candidate gene and genome wide association studies. For example, candidate gene studies in species ranging from humans²⁴⁻²⁶ to wild great tits²⁷⁻²⁹ have linked polymorphisms in the dopamine receptor D4 (DRD4) gene to personality variation in novelty seeking and exploratory behavior. Moreover, genes linked to variation in cognitive traits have also been identified in some cases. To give a few examples, catechol-O-methyltransferase (COMT) gene variants influence insight problem-solving in humans³⁰, and polymorphism in the arginine vasopressin V1a receptor gene (AVPR1A) affects social cognition in chimpanzees³¹.

Though more limited in scope than genome wide association studies, candidate gene studies allow researchers to access specific hypotheses regarding relationships between behavioral traits and genes of interest¹³,³². Genes related to the dopaminergic and serotonergic neurotransmitter systems are good candidates for affecting variation in personality and cognitive performance patterns²³⁻²⁹,³³⁻⁴¹. Serotonin is a monoamine hormone that has multifaceted physiological, behavioral and cognitive effects in both vertebrates and invertebrates. Serotonin acts as a critical neurotransmitter, but is also prevalent outside of the central nervous system and plays central roles in a diversity of physiological functions, including cardiovascular function and gut motility. Thus, effects of serotonin function on behavioral patterns may conceivably arise through modulation of physiological patterns, as well as via important effects on the central nervous system⁴². Effects of serotonin function on behavior and cognition have been most extensively studied in humans, wherein serotonin levels are inversely related to anxiety, impulsiveness and aggressiveness⁴³,⁴⁴. Serotonin dysfunction in nonhuman primates and model laboratory rodents has also been linked to similar behavioral traits⁴⁵,⁴⁶, and shown to interfere with cognitive processes⁴⁶⁻⁴⁸. The serotonin transporter gene (SERT), which was the focus of our study, is a member of the neurotransmitter sodium symporter family. The serotonin transporter (SERT) protein selectively cotransports serotonin (5-HT) out of the synaptic cleft and into nerve cells along with one Na⁺ and one Cl⁻ ion, while one K⁺ is transported in the opposite direction. Uptake of extracellular serotonin (5-HT) by SERT terminates the action of serotonin at receptor sites after its release, and thus serves as a key regulator of serotonergic signaling⁴⁹. SERT polymorphism in humans⁵⁰⁻⁵³, nonhuman primates⁵⁴⁻⁵⁶ and model rodents⁵⁷ has been implicated in an array of behavioral and cognition-related traits, especially sensitivity to stress, anxiety and social cognition. Variation in SERT function may thus affect both personality and cognitive traits in animals via changes in serotonergic signaling.

Avian studies that suggest a role for *SERT* polymorphisms in mediating genetic differences in behavior are growing in scope. Some avian studies have found differences in *SERT* polymorphisms between urban and rural bird populations³⁷. Thus, *SERT* polymorphisms that differ between urban and rural populations could underlie some of the aforementioned behavioral differences between the two¹⁷⁻²², including differences in neophobia¹⁸, exploratory behavior¹⁹, and innovativeness²¹. A few other avian studies have linked *SERT* polymorphisms to individual differences in personality traits. In great tits (*Parus major*), *SERT* SNPs were related to neophobia during parental care⁴⁰ and hissing behavior during nest defense⁴¹. Research that replicates and expands upon such studies is needed to elucidate whether *SERT* polymorphism is implicated in personality variation across multiple great tit populations, and to what extent the same polymorphisms are involved. Our research team recently failed to replicate the previous finding that *SERT* SNP187 in exon 1 is related to hissing behavior in great tits⁴¹, but found some association between two other *SERT* SNPs and hissing behavior (SNP66 in exon 13 and SNP144 in exon 12)⁵⁸. In the same study, we found no associations between female-female aggression and *SERT* SNPs, even when taking age-related plasticity in aggression into account⁵⁸. Thus, different *SERT* polymorphisms might be responsible for behavioral variation across different populations and behavioral contexts.

Avian studies regarding the genetics of problem-solving are few in number, and have not yet investigated effects of *SERT* polymorphisms on cognitive traits. This is despite a growing body of evidence linking serotonin and cognitive processes in other animal taxa⁴⁶,⁵³,⁵⁶, although note that documented relationships between serotonin and cognitive processes are often indirect effects mediated via individual differences in stress-sensitivity⁴⁸,⁵³. In addition, although differences in problem-solving performance have been documented between urban and rural bird populations²¹, potential genetic underpinnings of these differences remain largely unidentified. However, a few previous avian studies in free-living populations have investigated the genetic underpinnings of problem-solving performance. For example, expression patterns of the N-methyl-D-aspartate (NMDA) glutamate receptor were linked to differences in innovative problem-solving between two species of closely related finches⁵⁹. On the other hand, a study in a pedigreed population of great tits suggested little heritability of problem-solving performance¹⁶.

Using great tits as a study species, we aimed to expand upon past avian studies involving relationships between *SERT* polymorphisms and personality traits (boldness, neophobia and exploratory behavior), and perform a first test of relationships between *SERT* polymorphisms and problem-solving in birds. We previously found that birds occupying more disturbed territories, closer to roads and paths, are less neophobic and more likely to problem solve⁶⁰,⁶¹. We thus specifically aimed to determine whether *SERT* polymorphisms could provide a genetic basis for differences in neophobia and problem-solving between birds on more and less disturbed territories. We assessed associations between SNPs in the promoter region and 13 exons of the great tit *SERT* gene and behaviors measured during three standardized tests: (1) boldness and novel object (neophobia) tests (2) novel environment exploration tests, and (3) obstacle removal problem-solving tests. Neophobia and exploratory
behavior show moderate repeatability in our study population ($r \ge 0.30$), and thus constitute putative personality traits$^{60-64}$. We also evaluated whether *SERT* SNPs were related to reproductive fitness (nestlings fledged, laying date), as suggested by a previous study in great tits$^{41}$. Understanding genetic underpinnings of fitness traits may grant insight into how behavioral variation is maintained in populations.

## Materials and methods

**Study site.** We researched great tits from 5 study sites: Umicore (UM; 51°09'58.28"N, 4°21'23.76"E), Fort 7 (F7; 51°09'52.55"N, 4°22'40.46"E), Fort 6 (F6/Campus Drie Eiken), Fort 5 (F5; 51°10'02.90"N, 4°26'00.26"E) and Fort 4 (F4; 51°10'23.80"N, 4°27'38.46"E) located in the greater metropolis area of Antwerp, Belgium. Only novel environment exploration testing was performed at F6, to avoid conflicts with on-going studies. Study populations are not genetically isolated, with regular instances of dispersal occurring between sites. Most adults at all sites are color banded and females are equipped with passive integrated transponder (PIT) tags.

Study sites are located at varying distances from the Umicore smelter and metal refinery facility, a metal pollution point source, in Hoboken (south of Antwerp city center). In previous studies, we found no effect of proximity to Umicore on neophobia or problem-solving$^{60,61}$, but found a relationship between metal exposure and exploratory behavior$^{63,64}$. We do not consider gene-by-environment interactions mediated by metal exposure$^{65,66}$ due to lack of individual metal exposure measurements and because testing for genotype-by-site interactions would require a larger sample size. Rather, we include site as a random effect in statistical models to account for potential site-level differences in behavior. Within study sites, territories differ in anthropogenic disturbance levels, with territories closer to site edges being closer to roads, paths and buildings$^{60,61,67}$. As a proxy for disturbance level, we measured the distance of nest boxes from the nearest road or path using Google Earth measurement tools$^{60,61,67}$.

**Behavioral tests.** We tested 93 genotyped, incubating females on novel object tests$^{60}$, 156 genotyped birds (92 females, 62 males) on novel environment exploration tests$^{62-64}$, and 93 genotyped birds (52 females, 41 males) for problem-solving$^{61}$. Novel object tests consisted of two trial types, in random order, zero to three days apart (2018 and 2019, March to April). During novel object trials, we flushed females from the nest, placed a novel object on the nest box, and video-recorded for 1-h to determine the female's latency to return and enter$^{60}$. Baseline trials followed the same protocol, but with no novel object. We determined the female's boldness following disturbance at the nest from baseline trials, and additionally determined neophobia from novel object trials. However, note that return latencies during novel object trials are also affected by boldness as measured during baseline trials. Return latencies are individually repeatable across novel object and baseline trials, and across repeated novel object trials, suggesting repeatability of both boldness and neophobia$^{60}$.

Novel environment exploration tests occurred during winter (2017 to 2019, December to early March). Birds sleeping in nest boxes were captured at night, transported to the University of Antwerp's Campus Drie Eiken, and housed overnight in individual cages with food and water. The next morning (0800-1100), we performed tests using standard protocol$^{68}$. Birds were individually released into a novel environment room containing 5 artificial trees via sliding doors that connect cages to the room. We counted hops and flights made within 2-min, and calculated exploration scores as the total number of movements. After all individuals were tested, birds were weighed, measured, aged, sexed and banded (if not already banded), and released at capture sites before noon. Novel environment exploration behavior is repeatable in our populations$^{62-64,68}$.

Obstacle removal tests occurred in 2019 at UM, F7 and F4 during the nestling stage (nestlings 8-12 days old, April 26-May 22, 0930-1630)$^{61}$. We first placed a barrier (a modified nest box trap) on top of the nest box secured it in an open position, and video-recorded nestling provisioning for 1-h. This period allowed parent birds to habituate to the barrier, to minimize effects of neophobia on problem-solving. After one hour, we returned to the nest box and lowered the barrier. To create a mechanism whereby birds could lift the barrier, a string was looped around the metal hinge at the top of the trap and a wooden popsicle stick attached to the string. Birds could raise the barrier by landing and pulling on the stick, while lifting the barrier with their heads. We video-recorded nest boxes for 1.25 h. with the barrier lowered, and determined latency to contact the apparatus, time on the nest box, incorrect attempts directed at the barrier (pecks at the barrier) and number of entries. Note that birds that entered the nest box multiple times during the problem-solving test mastered the problem-solving task faster, or more thoroughly, than those that entered only once. Thus, number of entries grants additional information about individual problem-solving performance. Latency to contact the barrier might be influenced by both the motivation of birds to enter the nest box, and neophobia. However, effects of neophobia are expected to be minimized, because we habituated birds to the test apparatus before lowering the barrier$^{61}$. Obstacle removal tests were mostly performed with the barrier elevated 1.5 cm from the nest box face, to facilitate entry. A subset of 31 pairs (for 28 pairs at least one pair member was successful on the first test) was tested a second time with the barrier flat against the nest box$^{61}$. Trials for which an individual was never observed by the nest box were dropped from analyses, and nests with similar brood sizes were selected whenever possible. Testing protocols are described in detail in previous publications$^{60-64}$.

**Fitness-related breeding parameters.** During 2018 and 2019, we monitored reproduction via nest box checks every other day (late March through May). We recorded the date that the first egg was laid, clutch size, brood size, and fledging success. We analyzed relationships between SNPs and female laying date, and number of nestlings fledged for both sexes, because these parameters were the most variable and least redundant.

**Genotyping.** To identify SNPs, we first amplified the great tit *SERT* promoter region and thirteen exons using PCR. Most primers were as in Timm et al. (2018)$^{41}$, but those for exon 1 and 12 were redesigned using
<table><caption>**Table 1.** Twenty-six SNPs in the *SERT* gene and associated amino acid (AA) changes for exonic SNPs found in great tits in Antwerp, Belgium. Loci are named according to nucleotide positions within the promoter region or exons. The coordinate column gives the position within the great tit genome on chromosome 19. M = major allele, m = minor allele, %m = minor allele frequency. Chi-square tests were used to assess the null hypothesis of Hardy-Weinberg (HW) equilibrium. P values in red italics indicate significant deviations from HW (Holm's corrected P<sub>critical</sub> = 0.002). Sample sizes differ between loci due to sequencing failures. SNPs considered in statistical analyses are shown in bold text.</caption><thead><tr><th>Locus</th><th>Coordinate</th><th>Location</th><th>M/m</th><th>mm</th><th>Mm</th><th>MM</th><th>%m</th><th>X<sub>1</sub><sup>2</sup></th><th>P</th><th>Protein coding</th><th>AA change</th></tr></thead><tbody><tr><td>SNP030</td><td>chr19:5979530</td><td>Promoter</td><td>c/t</td><td>8</td><td>30</td><td>145</td><td>12.57</td><td>11.81</td><td>0.001</td><td></td><td></td></tr><tr><td>SNP031</td><td>chr19:5979529</td><td>Promoter</td><td>a/t</td><td>9</td><td>61</td><td>114</td><td>21.47</td><td>0.051</td><td>0.820</td><td></td><td></td></tr><tr><td>SNP100</td><td>chr19:5979460</td><td>Promoter</td><td>a/g</td><td>4</td><td>32</td><td>147</td><td>10.93</td><td>1.898</td><td>0.168</td><td></td><td></td></tr><tr><td>SNP115</td><td>chr19:5979445</td><td>Promoter</td><td>c/a</td><td>34</td><td>54</td><td>98</td><td>32.80</td><td>21.67</td><td>&lt;0.001</td><td></td><td></td></tr><tr><td>SNP250</td><td>chr19:5979310</td><td>Promoter</td><td>g/a</td><td>12</td><td>54</td><td>125</td><td>20.42</td><td>3.230</td><td>0.072</td><td></td><td></td></tr><tr><td>SNP288</td><td>chr19:5979272</td><td>Promoter</td><td>g/t</td><td>15</td><td>28</td><td>148</td><td>15.18</td><td>35.45</td><td>&lt;0.001</td><td></td><td></td></tr><tr><td>SNP415</td><td>chr19:5979145</td><td>Promoter</td><td>c/t</td><td>3</td><td>37</td><td>150</td><td>11.32</td><td>0.168</td><td>0.682</td><td></td><td></td></tr><tr><td>SNP438</td><td>chr19:5979122</td><td>Promoter</td><td>c/t</td><td>12</td><td>73</td><td>106</td><td>25.39</td><td>0.014</td><td>0.904</td><td></td><td></td></tr><tr><td>SNP592</td><td>chr19:5978968</td><td>Promoter</td><td>t/c</td><td>16</td><td>46</td><td>129</td><td>20.42</td><td>12.80</td><td>&lt;0.001</td><td></td><td></td></tr><tr><td>SNP663</td><td>chr19:5978897</td><td>Promoter</td><td>g/a</td><td>1</td><td>44</td><td>144</td><td>12.17</td><td>1.498</td><td>0.221</td><td></td><td></td></tr><tr><td>SNP163</td><td>chr19:5978840</td><td>Exon 1</td><td>g/a</td><td>4</td><td>47</td><td>145</td><td>14.03</td><td>0.007</td><td>0.933</td><td>Synonymous</td><td></td></tr><tr><td>SNP187</td><td>chr19:5978816</td><td>Exon 1</td><td>c/t</td><td>4</td><td>58</td><td>134</td><td>16.84</td><td>0.630</td><td>0.427</td><td>Synonymous</td><td></td></tr><tr><td>SNP226</td><td>chr19:5978777</td><td>Exon 1</td><td>t/a</td><td>20</td><td>72</td><td>104</td><td>28.57</td><td>1.96</td><td>0.162</td><td>Non-synonymous</td><td>E26D</td></tr><tr><td>SNP32</td><td>chr19:5977679</td><td>Exon 2</td><td>c/a</td><td>0</td><td>19</td><td>176</td><td>4.87</td><td>0.511</td><td>0.475</td><td>Synonymous</td><td></td></tr><tr><td>SNP100</td><td>chr19:5976873</td><td>Exon 3</td><td>c/t</td><td>2</td><td>28</td><td>168</td><td>8.08</td><td>0.457</td><td>0.499</td><td>Non-synonymous</td><td>A231V</td></tr><tr><td>SNP101</td><td>chr19:5976872</td><td>Exon 3</td><td>g/a</td><td>0</td><td>17</td><td>181</td><td>4.29</td><td>0.398</td><td>0.528</td><td>Synonymous</td><td></td></tr><tr><td>SNP125</td><td>chr19:5976812</td><td>Exon 3</td><td>g/a</td><td>3</td><td>32</td><td>163</td><td>9.60</td><td>0.929</td><td>0.335</td><td>Synonymous</td><td></td></tr><tr><td>SNP170</td><td>chr19:5976794</td><td>Exon 3</td><td>c/g</td><td>6</td><td>33</td><td>159</td><td>11.36</td><td>5.902</td><td>0.015</td><td>Synonymous</td><td></td></tr><tr><td>SNP187</td><td>chr19:5976777</td><td>Exon 3</td><td>t/a</td><td>2</td><td>21</td><td>175</td><td>6.31</td><td>2.116</td><td>0.146</td><td>Non-synonymous</td><td>L260Q</td></tr><tr><td>SNP51</td><td>chr19:5975781</td><td>Exon 5</td><td>c/t</td><td>0</td><td>16</td><td>178</td><td>4.12</td><td>0.358</td><td>0.549</td><td>Synonymous</td><td></td></tr><tr><td>SNP36</td><td>chr19:5973968</td><td>Exon 6</td><td>t/c</td><td>50</td><td>91</td><td>57</td><td>48.23</td><td>1.256</td><td>0.262</td><td>Synonymous</td><td></td></tr><tr><td>SNP48</td><td>chr19:5971868</td><td>Exon 9</td><td>g/a</td><td>0</td><td>23</td><td>175</td><td>5.81</td><td>0.752</td><td>0.386</td><td>Synonymous</td><td></td></tr><tr><td>SNP51</td><td>chr19:5971865</td><td>Exon 9</td><td>c/t</td><td>10</td><td>49</td><td>139</td><td>17.42</td><td>3.881</td><td>0.049</td><td>Synonymous</td><td></td></tr><tr><td>SNP84</td><td>chr19:5971832</td><td>Exon 9</td><td>t/c</td><td>41</td><td>100</td><td>57</td><td>45.96</td><td>0.055</td><td>0.814</td><td>Synonymous</td><td></td></tr><tr><td>SNP144</td><td>chr19:5968682</td><td>Exon 12</td><td>c/t</td><td>4</td><td>21</td><td>172</td><td>7.36</td><td>9.390</td><td>0.002</td><td>Synonymous</td><td></td></tr><tr><td>SNP66</td><td>chr19:5967914</td><td>Exon 13</td><td>c/t</td><td>14</td><td>64</td><td>101</td><td>25.70</td><td>0.727</td><td>0.394</td><td>Synonymous</td><td></td></tr></tbody></table>

Geneious Prime 2020.0.4 (Table S1). Thermocycling conditions and reaction mixes have been previously described⁵⁸. Direct sequencing of PCR amplicons was performed at the University of Antwerp's genomics core facility using a sanger sequencing platform. We used Geneious Prime 2020.0.4 to align sequences to the great tit reference genome (parus_major1.1:ENSPMIGJ00000019668, Ensembl genome browser) and identified SNPs via manual inspection. We initially genotyped between 46 and 48 individuals representing a wide range in behavioral phenotypes across all loci to determine the presence and frequencies of SNP. We then genotyped SNPs exceeding 5% minor allele frequency in the rest of the population (N = 181). SNPs with minor allele frequencies of <10% in the whole population were excluded from analyses. SNPs for which there were <4 individuals with the minor allele (minor allele homozygotes or heterozygotes) in behavioral datasets were also excluded from analyses due to model convergence problems. SNPs that deviated from Hardy-Weinberg equilibrium (assessed via Chi-square tests with a Holm's correction, Table 1) were additionally excluded.

Linkage disequilibrium was assessed with the web-based application SNPStats, which uses matrices of genotype data to calculate haplotypes, linkage disequilibrium statistics (*D*, *D'*, Pearson's *r*) and *P* values based on Chi-square tests⁶⁹. For exonic SNPs, we also determined whether SNPs were synonymous (not changing the amino acid sequence) or nonsynonymous (changing the sequence). Sample sizes varied slightly across loci as not all SNPs were successfully sequenced for all individuals.

To infer putative functional consequences of coding variants, we used SNAP2 to predict effects of variants on protein function⁷⁰. SNAP2 is a trained classifier based on a machine-learning device called neural network, which distinguishes between nonsynonymous SNPs that cause functional effects and neutral variants. The effect of a variant is believed to be important to native protein function and structure if the SNAP2 score exceeds 50; neutral if the score is below -50; and not possible to reliably determine when the score is between 50 and -50⁷⁰. SNAP2 results were verified using a different function prediction tool, PROVEAN⁷¹. In addition, we also used the web-based tool Protter⁷² to visualize the location of nonsynonymous SNPs within SERT.

**Statistical analyses.** We performed statistical analyses in R 3.6.1⁷³. We first used linear (LMMs; for modeling log-transformed return latencies) and generalized linear mixed models (GLMMs; binomial for problem-
solving success, Poisson for incorrect attempts)⁷⁴,⁷⁵ to investigate associations among behavioral and cognitive traits. Specifically, we tested the relationship between neophobia and boldness (log-transformed return latencies, dependent variable), and the first measured exploration score of females (predictor). Relationships between traits measured during problem-solving tests (dependent) and the putative personality traits were tested using results from first tests of personality traits as predictor variables. Relationships among traits measured during problem-solving tests were also tested. For all models, individual and study site were entered as random effects, and trial number as a covariate. When testing the above relationships, neophobia was adjusted for boldness by extracting residuals of an LMM predicting return latencies during novel object trials from return latencies during baseline trials. We also report the relationship between return latencies during novel object and baseline trials.

We assessed the repeatability of putative personality traits (neophobia, boldness, exploratory behavior) in the current dataset using R package rptR. Package rptR estimates repeatability based on variance components extracted from mixed models⁷⁶. Models used to estimate repeatability were linear mixed models (LMMs)⁷⁵ for boldness and neophobia, and a Poisson GLMM⁷⁴ for exploratory behavior. Repeatability estimates were adjusted for effects of year, test date and test number for exploratory behavior, and year and test number for boldness and neophobia. For neophobia, we assessed repeatability both when adjusting for baseline return latencies (see above), and when not making this adjustment.

To model effects of SNPs on behaviors, problem-solving and breeding parameters, we constructed general generic models, wherein genotypes were entered as three level factors. If results from general generic models suggested dominance (e.g. CT = TT ≠ CC) or overdominance (heterozygote ≠ homozygotes), we reran models combining genotypes. We analyzed the relationship between each SNP (in separate models to avoid collinearity and over complexity) and six behaviors: (1) female latency to return to the nest during novel object tests, wherein latency during baseline trials measured boldness and latency during novel object tests measured neophobia, (2) exploratory behavior, (3) problem-solving success, (4) number of entries during obstacle removal tests (results only reported where significant relationships were found with problem-solving success), (5) incorrect attempts on the barrier, and (6) latency to contact the obstacle removal test apparatus. For behavioral patterns during obstacle removal tests (problem-solving success, number of entries, incorrect attempts and latency to contact), we performed separate models for males and females, because random effects were overparameterized with both sexes in models. For males, a low number of problem-solving individuals in our dataset made testing associations with problem-solving success difficult for some SNPs (18 genotyped females, but only 9 genotyped males were successful). We also constructed models predicting two fitness-related parameters: female laying date (log-transformed) and number of fledglings in both sexes.

We had repeated measures for many individuals, and thus adopted a mixed modeling approach⁷⁴. Return latencies (log-transformed) during novel object tests, latency to contact during obstacle removal tests, and laying date (log-transformed) were analyzed using linear mixed effects models (LMMs) with Satterthwaite approximations for degrees of freedom (R package lmerTest)⁷⁵. Success on obstacle removal tests was analyzed using a generalized linear mixed effects model (GLMM) with a binomial distribution. Exploratory behavior, nest box entries, incorrect attempts on the barrier, and fledgling number were analyzed using Poisson GLMMs. For most models, we used individual identity and study site as random effects, and added year as a random when predicting behaviors measured across multiple years (neophobia, exploratory behavior). However, for problem-solving success in males, we used only study site as a random effect, because attempting to include an individual-level random effect led to convergence problems. Thus, only results from first problem-solving tests on individual males were included in analyses. Binomial and Poisson GLMMs included an observation level random effect to control for overdispersion, where necessary⁷⁷.

Some models contained additional interactions and covariates. Models predicting return latencies during novel object tests included two-way interactions between trial type (baseline, novel object) and genotype, and trial type and distance to the path/road (we previously found this interaction to be significant)⁶⁰. If neophobia were affected by SNPs, we expected statistically significant relationships only within novel object trials, whereas if boldness were affected by SNPs, we expected that relationships would be present within baseline trials or across both trial types. Trial number (first versus second test on a bird) was included as a covariate, whereas clutch size, trial date and female age were not considered, as these variables were nonsignificant in a previous study⁶⁰. Models for exploratory behavior included an interaction between sex and genotype, test date and trial number. Models for problem-solving success, nest box entries, incorrect attempts on the barrier, and latency to contact the barrier contained trial number and distance to paths/roads (square-root transformed). Brood size, nestling age, time and date were not included in these models, because these variables were found to be unrelated to problem-solving performance in a previous study⁶¹. However, to better understand potential behavioral mechanisms underlying relationships between SNPs and problem-solving success, we repeated models that included significant SNP-problem-solving associations with incorrect attempts on the barrier and latency to contact the test apparatus as additional predictor variables. Models involving incorrect attempts on the barrier as the dependent variable additionally used time on the nest box (log-transformed) as an offset variable. Lay date was initially included in models involving fledgling numbers, but was statistically nonsignificant, so was removed.

We performed posthoc tests for comparisons among genotypes using R package emmeans (Tukey method)⁷⁸, and applied Bonferroni corrections to account for the number of SNPs (N = 12; significance threshold α = 0.004). Although Bonferroni corrections are stringent, applying a different method, such as false discovery rate⁷⁹, would not qualitatively change conclusions. To obtain a measure of effect size for SNPs that showed substantial associations with behavior, we used marginal R² (R²m), calculated using the method of Nakagawa et al. (2013)⁸⁰ (function r2 in R package performance)⁸¹. R²m reflects the proportion of variance explained by fixed effects in a mixed effects model. We also report conditional R² (R²c), which reflects the proportion of variance explained by random plus fixed effects. Models were reduced using a stepwise process (α = 0.05) involving removing nonsignificant
interaction terms first, and then nonsignificant main effects when significant predictors were present. Model diagnostics were performed using R package DHARMA⁸². Sample sizes varied among the behavioral tests and among SNPs due to sequencing failures. We thus indicate sample sizes for all models in the results section.

**Ethical statement.** Animal experiments were approved by the University of Antwerp's ethical committee (ID 2016-71), and conducted in accordance with Belgian and Flemish laws. Methodology adhered to the ASAB/ABS guidelines for use of animals in behavioral research, and we made all possible efforts to minimize the stress experienced by birds during handling and behavioral testing. The Belgian Royal Institute for Natural Sciences (Koninklijk Belgisch Instituut voor Natuurwetenschappen; KBIN) provided banding licenses for authors and technical personnel.

## Results

**Relationships among traits and repeatability of behaviors.** Novel object neophobia (adjusted for baseline return latency) was related to exploration behavior ($\beta$ = -0.016 ± 0.007, $t_{78}$ = -2.092, *P* = 0.039), but did not predict problem-solving success ($\beta$ = 0.296 ± 0.303, *Z* = 0.978, *P* = 0.328), incorrect attempts ($\beta$ = -0.288 ± 0.975, *Z* = -295, *P* = 0.768), or latency to contact the test apparatus ($\beta$ = 0.075 ± 0.141, $t_{66}$ = 0.532, *P* = 0.596). Problem solving success ($\beta$ = -0.255 ± 0.321, *Z* = -0.794, *P* = 0.427), incorrect attempts ($\beta$ = 0.119 ± 0.482, *Z* = 0.247, *P* = 0.805), and latency to contact the test apparatus ($\beta$ = -0.867 ± 2.43, $t_{87}$ = -0.356, *P* = 0.722) were also not predicted by exploratory behavior. Return latency during novel object trials (unadjusted for baseline behavior) was related to return latency during baseline trials (boldness) ($\beta$ = 0.018 ± 0.004, $t_{146}$ = 4.281, *P* < 0.001), but boldness was not significantly related to exploratory behavior ($\beta$ = -0.001 ± 0.009, $t_{75}$ = -0.113, *P* = 0.91) or traits measured during problem-solving tests (*P* > 0.30 in all cases). Problem-solving success was negatively associated with incorrect attempts ($\beta$ = -0.905 ± 0.279, *Z* = -2.121, *P* = 0.033) and latency to contact the test apparatus ($\beta$ = -0.960 ± 0.427, *Z* = -3.430, *P* < 0.001).

Exploratory behavior ($r$ ± SE [95% CI] = 0.581 ± 0.041 [0.319, 0.481], *P* < 0.001) and novel object neophobia either adjusted ($r$ ± SE [95% CI] = 0.432 ± 0.103 [0.221, 0.609], *P* < 0.001) or unadjusted ($r$ ± SE [95% CI] = 0.473 ± 0.091 [0.290, 0.635], *P* < 0.001) for boldness showed statistically significant repeatability. Boldness was also significantly repeatable ($r$ ± SE [95% CI] = 0.436 ± 0.176 [0.052, 0.758], *P* = 0.047).

**Genetic polymorphisms.** Of 26 SNPs identified at a minor allele frequency of > 5% in our initial sample, 5 deviated from Hardy-Weinberg equilibrium, and 9 displayed a minor allele frequency of < 10% in the entire sample or too few individuals (< 4) with the minor allele in behavioral datasets, resulting in 12 SNPs retained for analysis (Table 1). Four retained SNPs were in the promoter region. Of the 8 exonic SNPs retained for analysis, only SNP226 in exon 1 was nonsynonymous, causing the amino acid change E26D (glutamic acid replaced by aspartic acid). Two other nonsynonymous SNPs located in exon 3 (SNP100, A231V and SNP187, L260Q) were present at < 10% minor allele frequency. None of the SNPs had been previously reported in great tits, although Timm et al. (2018, 2019)⁴⁰,⁴¹ identified a different nucleotide change at one of the positions (SNP187 in exon 1; A/T in previous study, C/T in current study). There was no complete linkage disequilibrium, but substantial linkage disequilibrium existed between some loci within and between the promoter region and exon 1 (Supplementary Table S2).

SNAP2 analysis could not reliably predict the functional significance of the amino acid change caused by SNP226 in exon 1 (E26D: strength: -42, reliability: 72%), the only nonsynonymous SNP retained in our analyses. For the other two amino acid changes identified (Table 1), SNAP2 analysis revealed a likely neutral effect for A231V caused by SNP100 in exon 3 (strength: -70, reliability: 82%), and failed to reliably predict functionality for the change L260Q associated with SNP187 in exon 3 (strength: -19, reliability: 57%). Visualization of the location of nonsynonymous SNPs using Protter demonstrated that E26D (SNP226 in exon 1) is located within the N terminus of SERT, whereas A231V (SNP100 exon 3) and L260Q (SNP187 exon 3) are located within extracellular loop 2 (Supplementary Fig. S1).

**Relationships between SNPs and behaviors.** *Novel object neophobia, boldness and exploratory behavior.* Genotype at SNP36 in exon 6 and trial type interacted to predict return latency during novel object tests ($F_{2,173}$ = 3.570, *P* = 0.030). Females heterozygous for SNP36 (CT) returned to the nest box faster during novel object trials than homozygous birds, suggesting an effect on neophobia (genotype CC or TT; $\beta$ = 0.402 ± 0.166, $t_{104}$ = 2.418, *P* = 0.017; $R^2m$ = 0.048, $R^2c$ = 0.49), whereas heterozygotes and homozygotes did not differ in return latencies during baseline trials, suggesting no effect on boldness following disturbance at the nest ($\beta$ = -0.049 ± 0.177, $t_{83}$ = -0.276, *P* = 0.783) (Fig. 1a). However, the interaction was nonsignificant after accounting for multiple comparisons, as was also the case for all other SNP-behavior relationships examined, with one exception, noted below.

Promoter region SNP438 was associated with exploratory behavior ($\chi^2$ = 6.830, *P* = 0.033), with heterozygotes (CT) tending to have higher exploration scores than homozygotes ($\beta$ = -0.592 ± 0.246, *Z* = -2.405, *P* = 0.016; $R^2m$ = 0.035, $R^2c$ = 0.463) (Fig. 1b). We found no other relationships between promoter region or exonic SNPs and either neophobia/boldness (Table S3, S4) or exploratory behavior (Table S5, S6). Including distance to paths/roads in models did not alter conclusions.

**Problem-solving performance.** Incorrect attempts and latency to contact. In females, but not males, there was a relationship between SNP36 in exon 6 and incorrect attempts during obstacle removal tests (Table 2a; Fig. 2a). As was also the case for novel object neophobia (see above), heterozygous females (CT) differed from homozygotes, with heterozygotes pecking at the barrier more (Table 2a; Fig. 2a). The overall relationship between SNP36
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Relationships between SERT SNPs, neophobia/boldness and exploratory behavior. Mean values and standard errors (error bars) for groups were extracted using R package emmeans and back transformed from the log scale; (a) exon 6 SNP36 and return latency (minutes) during novel object tests. Triangles indicate baseline trials (boldness) and circles novel object trials (neophobia). N (tests, females) = baseline: CC: 24, 20; CT: 48, 41; TT: 37, 29; novel object: CC: 31, 21; CT: 74, 43; TT: 50, 28. (b) Promoter region SNP438 and exploratory behavior (movements/2 min). N (tests, birds) = CC: 133, 87; CT: 81; TT: 12, 9.</figcaption></figure>

and incorrect attempts on the barrier was significant at the Bonferroni threshold level, although this was not the case for pairwise comparisons (Table 2a).

Nonsynonymous SNP226 (Table 2b; Fig. 2b) and synonymous SNP187 in exon 1 (Table 2c; Fig. 2c) were also associated with incorrect attempts by females. Females with rare genotype AA at SNP226 pecked the barrier less than females with AT or TT, although the *P* value for the AA-AT contrast was > 0.05 after posthoc correction (Table 2b; Fig. 2b). Females with genotype CC at SNP187 pecked the barrier more than birds with CT. Behavior did not differ between females with genotype TT at SNP187 and those with the other genotypes (Table 2c; Fig. 2c), but few (three) females exhibited genotype TT, making these comparisons difficult.

No other promoter region or exonic SNP was related to incorrect attempts on the barrier in either females (Table S7, S8) or males (Table S9, S10). We also found no statistically significant relationships between any of the SNPs and latency to contact the test apparatus in either females (S11, S12) or males (S13, S14).

**Problem-solving success.** Nonsynonymous SNP226 in exon 1 was related to female problem-solving success. Pairwise contrasts from the general generic model were nonsignificant after posthoc corrections for multiple comparisons (Table 3a), but in a dominance effect model, females with genotype AA were more likely to problem-solve than females with genotype AT or TT at $α=0.05$ (Table 3b; Fig. 3a), and entered the nest box more times (dominance effect model, $β=1.714±0.855, Z=2.128, P=0.045$). With distance to paths/roads in the dominance model, there was an effect of both genotype ($β=1.694±0.787, Z=2.150, P=0.031$) and distance from paths/roads ($β=-1.239±0.588, Z=-2.106, P=0.035$) on female problem-solving success. We found no statistical support for an effect of SNP226 on male problem-solving success (Table S18). However, it was difficult to test the effect of SNP226 on male problem-solving, because only two males had genotype AA.

In addition, we found some evidence for an effect of SNP84 in exon 9 on problem-solving success, with this association being stronger in males than females (Table 4; Fig. 3b,c). In males, pairwise contrasts from the general generic model were nonsignificant after posthoc corrections for multiple comparisons (Table 4a), but in a dominance model males with genotype CC were more likely to problem-solve than males with CT or TT ($β=3.504±0.1526, Z=2.296, P=0.021$; Table 4b), and tended to enter the nest box more times ($β=3.160±1.807, Z=1.749, P=0.080$). The effect size associated with the SNP84-problem-solving success association in males was large ($R^2m=0.292$, general generic model). In females, the relationship between SNP84 in exon 9 and problem-solving success was nonsignificant, but was in the same direction as the effect observed in males, and was associated with a substantial effect size, comparable to that for SNP226 in exon 1 (Table 4c, d). With distance to paths/roads in the dominance models, there was some evidence for an effect of both SNP84 (male: $β=3.727±1.435, Z=2.597, P=0.009$; females: $β=1.296±0.7314, Z=-1.757, P=0.087$) and distance from paths/roads (males: $β=-0.348±0.205, Z=-1.694, P=0.090$; females: $β=-0.190±0.093, Z=-2.040, P=0.041$) on problem-solving success in both sexes. No other promoter region or exonic SNP was associated with problem-solving success in either females (Table S15, S16) or males (Table S17, S18).

We added incorrect attempts and latency to contact the test apparatus to dominance models predicting female problem-solving success from SNP226 and SNP84, to investigate how considering these behavioral traits might modify the relationship between SNPs and problem-solving. In genotyped males, incorrect attempts and latency to contact unexpectedly did not have significant effects on problem-solving success (*P* > 0.20), so we did not add these variables to the model involving SNP84 and problem-solving in males. The effect of SNP226 on problem-solving success in females was no longer statistically significant ($β=-1.066±0.808, Z=-1.320, P=0.186$) when including these variables in the model, while both incorrect attempts ($β=-0.048±0.025, Z=-1.925, P=0.054$) and latency to contact ($β=-0.01±0.005, Z=-1.928, P=0.054$) had near significant effects. A similar pattern was observed when adding latency to contact and incorrect attempts to the model involving SNP84 and female problem-solving success.
(a) SNP36 exon 6

<table><thead><tr><th colspan="2">General generic model</th><th colspan="2">R<sup>2</sup>m = 0.300, 0.181<sup>a</sup></th><th colspan="2">R<sup>2</sup>c = 0.381</th></tr><tr><th>Fixed effects</th><th>β ± SE</th><th>Z</th><th>P</th><th>χ<sub>2</sub><sup>2</sup></th><th>P</th></tr></thead><tbody><tr><td>Intercept</td><td>0.697 ± 1.235</td><td>0.565</td><td>0.572</td><td></td><td></td></tr><tr><td>SNP36 CC</td><td>-2.811 ± 1.128</td><td>2.491</td><td>0.034<sup>b</sup></td><td></td><td></td></tr><tr><td>SNP36 TT</td><td>-1.570 ± 0.608</td><td>-2.581</td><td>0.026<sup>b</sup></td><td>11.34</td><td>0.003</td></tr><tr><td>Trial number</td><td>-1.895 ± 0.641</td><td>-2.954</td><td>0.003</td><td></td><td></td></tr><tr><th>Random effects</th><th>Variance</th><th>SD</th><th>N</th><th></th><th></th></tr><tr><td>Observation</td><td>2.754</td><td>1.659</td><td>73</td><td></td><td></td></tr><tr><td>Individual</td><td>0.510</td><td>0.714</td><td>52</td><td></td><td></td></tr><tr><td>Study site</td><td>&lt; 0.001</td><td>0.002</td><td>3</td><td></td><td></td></tr></tbody></table>

(b) SNP226 exon 1

<table><thead><tr><th colspan="2">General generic model</th><th colspan="2">R<sup>2</sup>m = 0.214, 0.145<sup>a</sup></th><th colspan="2">R<sup>2</sup>c = 0.416</th></tr><tr><th>Fixed Effects</th><th>β ± SE</th><th>Z</th><th>P</th><th>χ<sub>2</sub><sup>2</sup></th><th>P</th></tr></thead><tbody><tr><td>Intercept</td><td>0.494 ± 1.175</td><td>0.421</td><td>0.6739</td><td></td><td></td></tr><tr><td>SNP226 AT</td><td>1.849 ± 1.002</td><td>1.845</td><td>0.155<sup>c</sup></td><td></td><td></td></tr><tr><td>SNP226 TT</td><td>2.268 ± 0.952</td><td>2.382</td><td>0.045<sup>c</sup></td><td>6.335</td><td>0.042</td></tr><tr><td>Trial number</td><td>-1.549 ± 0.612</td><td>-2.53</td><td>0.011</td><td></td><td></td></tr><tr><th>Random effects</th><th>Variance</th><th>SD</th><th>N</th><th></th><th></th></tr><tr><td>Observation</td><td>2.385</td><td>1.544</td><td>73</td><td></td><td></td></tr><tr><td>Individual</td><td>1.218</td><td>1.104</td><td>52</td><td></td><td></td></tr><tr><td>Study site</td><td>&lt; 0.001</td><td>&lt; 0.001</td><td>3</td><td></td><td></td></tr></tbody></table>

<table><thead><tr><th colspan="2">Dominance effect model</th><th colspan="2">R<sup>2</sup>m = 0.212, 0.141<sup>a</sup></th><th colspan="2">R<sup>2</sup>c = 0.397</th></tr><tr><th>Fixed effects</th><th>β ± SE</th><th>Z</th><th>P</th><th>χ<sub>2</sub><sup>2</sup></th><th>P</th></tr></thead><tbody><tr><td>Intercept</td><td>0.510 ± 1.182</td><td>0.432</td><td>0.665</td><td></td><td></td></tr><tr><td>SNP226 AT/TT</td><td>2.114 ± 0.913</td><td>2.315</td><td>0.020</td><td></td><td></td></tr><tr><td>Trial number</td><td>-1.555 ± 0.620</td><td>-2.508</td><td>0.012</td><td></td><td></td></tr><tr><th>Random effects</th><th>Variance</th><th>SD</th><th>N</th><th></th><th></th></tr><tr><td>Observation</td><td>2.510</td><td>1.580</td><td>73</td><td></td><td></td></tr><tr><td>Individual</td><td>1.110</td><td>1.060</td><td>52</td><td></td><td></td></tr><tr><td>Study site</td><td>&lt; 0.001</td><td>&lt; 0.001</td><td>3</td><td></td><td></td></tr></tbody></table>

(c) SNP187 exon 1

<table><thead><tr><th colspan="2">General generic model</th><th colspan="2">R<sup>2</sup>m = 0.210, 0.126<sup>a</sup></th><th colspan="2">R<sup>2</sup>c = 0.378</th></tr><tr><th>Fixed effects</th><th>β ± SE</th><th>Z</th><th>P</th><th>χ<sub>2</sub><sup>2</sup></th><th>P</th></tr></thead><tbody><tr><td>Intercept</td><td>2.959 ± 0.799</td><td>3.703</td><td>&lt; 0.001</td><td></td><td></td></tr><tr><td>SNP187 TC</td><td>-1.814 ± 0.700</td><td>-2.59</td><td>0.026<sup>d</sup></td><td></td><td></td></tr><tr><td>SNP187 TT</td><td>-0.665 ± 1.235</td><td>-0.538</td><td>0.852<sup>d</sup></td><td>7.130</td><td>0.028</td></tr><tr><td>Trial number</td><td>-1.660 ± 0.625</td><td>-2.656</td><td>0.007</td><td></td><td></td></tr><tr><th>Random effects</th><th>Variance</th><th>SD</th><th>N</th><th></th><th></th></tr><tr><td>Observation</td><td>2.543</td><td>1.595</td><td>73</td><td></td><td></td></tr><tr><td>Individual</td><td>0.993</td><td>0.996</td><td>52</td><td></td><td></td></tr><tr><td>Study site</td><td>&lt; 0.001</td><td>&lt; 0.001</td><td>3</td><td></td><td></td></tr></tbody></table>

**Table 2.** General generic and dominance effect GLMMs for relationships between exonic *SERT* SNPs and incorrect attempts on the barrier by female great titls during obstacle removal problem-solving tests. Time on the nest box was entered as an offset variable. (a) SNP36 exon 6, (b) SNP226 exon 1, (c) SNP187 exon 1. <sup>a</sup>Marginal R<sup>2</sup> for all fixed effects, and the SNP alone. <sup>b</sup>P values for post-hoc Tukey tests with CT as reference, P = 0.534 for the CC–TT contrast. <sup>c</sup>P values for post-hoc Tukey tests, P = 0.843 for the AT–TT contrast. <sup>d</sup>P values from posthoc Tukey tests, P = 0.657 for the TT–TC contrast. N = 73 observations, 52 females; exon 6 SNP36: 9 CC, 28 CT, 15 TT (a); exon 1 SNP226: 9 AA, 16 AT, 27 (b) TT; exon 1 SNP187: 33 CC, 16 CT, 3 TT (c).

**SERT** polymorphisms, distribution patterns and breeding parameters. We found no statistically significant relationships between any SNP and distance from paths/roads (Table S19, S20). SNP170 in exon 3
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Relationships between SERT SNPs and incorrect attempts directed at the barrier during obstacle removal tests. Mean values and standard errors (error bars) for groups were extracted using R package emmeans and back transformed from the log scale. Values are incorrect attempts per minute on the nest box; (a) Exon 6 SNP36, N (tests, females) = CC: 12, 9; CT: 39, 28; TT: 22, 15. (b) Exon 1 SNP226, N = AA: 13, 9; AT: 20, 16; TT: 40, 27, (c) exon 1 SNP187, N = CC: 47, 34; CT: 21, 15; TT: 5, 3.</figcaption></figure>

<table><caption>**Table 3.** Binomial GLMM for the relationship between nonsynonymous SNP226 in exon 1 and success on obstacle removal tests in female great tits. (a) General generic model, (b) dominance effect model with AT = TT. <sup>a</sup>P values from posthoc Tukey test. P = 0.974 for the AT–TT contrast. N = 73 observations on 52 females; 9 AA, 16 TA, 27 TT.</caption><thead><tr><th colspan="2">(a) General generic model</th><th colspan="2">R<sup>2</sup>m = 0.083</th><th colspan="2">R<sup>2</sup>c = 0.280</th></tr><tr><th>Fixed effects</th><th>β ± SE</th><th>Z</th><th>P</th><th>χ<sup>2</sup></th><th>P</th></tr></thead><tbody><tr><td>Intercept</td><td>0.532 ± 0.851</td><td>0.625</td><td>0.532</td><td></td><td></td></tr><tr><td>SNP226 AT vs AA</td><td>-1.693 ± 0.872</td><td>-1.942</td><td>0.127<sup>a</sup></td><td></td><td></td></tr><tr><td>SNP226 TT vs AA</td><td>-1.551 ± 0.810</td><td>-1.914</td><td>0.134<sup>a</sup></td><td>5.238</td><td>0.072</td></tr><tr><th>Random effects</th><th>Variance</th><th>SD</th><th>N</th><th></th><th></th></tr><tr><td>Individual</td><td>0.064</td><td>0.253</td><td>52</td><td></td><td></td></tr><tr><td>Study site</td><td>0.832</td><td>0.912</td><td>3</td><td></td><td></td></tr></tbody></table><table><thead><tr><th colspan="2">(b) Dominance effect model</th><th colspan="2">R<sup>2</sup>m = 0.084</th><th colspan="2">R<sup>2</sup>c = 0.280</th></tr><tr><th>Fixed effects</th><th>β ± SE</th><th>Z</th><th>P</th><th></th><th></th></tr></thead><tbody><tr><td>Intercept</td><td>0.535 ± 0.849</td><td>0.630</td><td>0.528</td><td></td><td></td></tr><tr><td>SNP226 (AT/TT vs AA)</td><td>-1.609 ± 0.775</td><td>-2.077</td><td>0.037</td><td></td><td></td></tr><tr><th>Random effects</th><th>Variance</th><th>SD</th><th>N</th><th></th><th></th></tr><tr><td>Individual</td><td>0.085</td><td>0.292</td><td>52</td><td></td><td></td></tr><tr><td>Study site</td><td>0.811</td><td>0.901</td><td>3</td><td></td><td></td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Relationships between SERT SNPs and problem-solving success. Mean values and standard errors (error bars) for groups were extracted using R package emmeans and back transformed from the logit scale; (a) exon 1 SNP226 in females, N (test, females) = AA: 13, 9; AT: 20, 16; TT: 40, 27, (b) exon 9 SNP84 in males, N (males) = CC: 5; CT: 20; TT: 16, (c) exon 9 SNP84 in females. N (tests, females) = CC: 19, 11; CT: 37, 28; TT: 17, 13.</figcaption></figure>

(R<sup>2</sup>m = 0.022, general generic model) and SNP84 in exon 9 (R<sup>2</sup>m = 0.016, general generic model) showed associations with laying date, but effect sizes were relatively small (Table 5a,b; Fig. 4; Table S21, S22 show models for all SNPs). In general generic models, posthoc pairwise comparisons were nonsignificant at α = 0.05 (Table 5). However, in dominance effects models, females with genotypes CG and CC at SNP170 in exon 3 lay earlier than females with GG (β = 0.041 ± 0.019, t<sub>114</sub> = 2.134, P = 0.034; Table 5a), and females with genotypes CT and TT at
<table><caption>**Table 4.** Binomial GLMM showing the relationship between synonymous SNP84 in exon 9 and success on the obstacle removal test in male and female great tits; (a) general generic model, males (b) dominance effect model, males with CT = TT, (c) general generic model, males (d) dominance effect model, females with CT = TT. <sup>a</sup>P values from posthoc Tukey test. <sup>P</sup> = 0.764, 0.717 for the CT–TT contrast in males and females, respectively. <sup>N</sup> = 41 males, 5 CC, 20 CT, 16 TT (a, b); 73 observations on 52 females, 11 CC, 28 CT, 13 TT (c, d).</caption><thead><tr><th colspan="2">(a) General generic model, males</th><th colspan="2">R<sup>2</sup>m = 0.292</th><th colspan="2">R<sup>2</sup>c = 0.324</th></tr></thead><tbody><tr><th colspan="2">Fixed effects</th><th>Z</th><th>P</th><th>χ<sub>2</sub><sup>2</sup></th><th>P</th></tr><tr><td colspan="2">Intercept</td><td>1.531 ± 1.307</td><td>1.172</td><td>0.241</td><td></td></tr><tr><td colspan="2">SNP84 CT vs CC</td><td>-3.771 ± 1.585</td><td>-2.380</td><td>0.045<sup>a</sup></td><td></td></tr><tr><td colspan="2">SNP84 TT vs CC</td><td>-1.551 ± 0.810</td><td>-1.914</td><td>0.134<sup>a</sup></td><td>9.77</td></tr><tr><th colspan="2">Random effects</th><th>Variance</th><th>SD</th><th>N</th><th></th></tr><tr><td colspan="2">Study site</td><td>0.155</td><td>0.394</td><td>3</td><td></td></tr></tbody><thead><tr><th colspan="2">(b) Dominance effect, males</th><th colspan="2">R<sup>2</sup>m = 0.277</th><th colspan="2">R<sup>2</sup>c = 0.324</th></tr></thead><tbody><tr><th colspan="2">Fixed effects</th><th>Z</th><th>P</th><th>χ<sub>2</sub><sup>2</sup></th><th>P</th></tr><tr><td colspan="2">Intercept</td><td>1.595 ± 1.327</td><td>1.202</td><td>0.229</td><td></td></tr><tr><td colspan="2">SNP84 (CT/TT vs CC)</td><td>-3.504 ± 1.526</td><td>-2.296</td><td>0.021</td><td></td></tr><tr><th colspan="2">Random effects</th><th>Variance</th><th>SD</th><th>N</th><th></th></tr><tr><td colspan="2">Study site</td><td>0.228</td><td>0.478</td><td>3</td><td></td></tr></tbody><thead><tr><th colspan="2">(c) General generic model, females</th><th colspan="2">R<sup>2</sup>m = 0.092</th><th colspan="2">R<sup>2</sup>c = 0.310</th></tr></thead><tbody><tr><th colspan="2">Fixed Effects</th><th>Z</th><th>P</th><th>χ<sub>2</sub><sup>2</sup></th><th>P</th></tr><tr><td colspan="2">Intercept</td><td>0.036 ± 0.752</td><td>0.049</td><td>0.960</td><td></td></tr><tr><td colspan="2">SNP84 CT vs CC</td><td>-1.536 ± 0.801</td><td>-1.917</td><td>0.133<sup>a</sup></td><td></td></tr><tr><td colspan="2">SNP84 TT vs CC</td><td>-0.624 ± 0.803</td><td>-0.777</td><td>0.443<sup>a</sup></td><td>4.969</td></tr><tr><th colspan="2">Random effects</th><th>Variance</th><th>SD</th><th>N</th><th></th></tr><tr><td colspan="2">Individual</td><td>0.288</td><td>0.536</td><td>52</td><td></td></tr><tr><td colspan="2">Study site</td><td>0.753</td><td>0.868</td><td>3</td><td></td></tr></tbody><thead><tr><th colspan="2">(d) Dominance model, females</th><th colspan="2">R<sup>2</sup>m = 0.059</th><th colspan="2">R<sup>2</sup>c = 0.296</th></tr></thead><tbody><tr><th colspan="2">Fixed Effects</th><th>Z</th><th>P</th><th>χ<sub>2</sub><sup>2</sup></th><th>P</th></tr><tr><td colspan="2">Intercept</td><td>0.036 ± 0.752</td><td>0.049</td><td>0.960</td><td></td></tr><tr><td colspan="2">SNP84 (TT/CT vs CC)</td><td>-1.193 ± 0.708</td><td>-1.684</td><td>0.092</td><td></td></tr><tr><th colspan="2">Random effects</th><th>Variance</th><th>SD</th><th>N</th><th></th></tr><tr><td colspan="2">Individual</td><td>0.306</td><td>0.553</td><td>52</td><td></td></tr><tr><td colspan="2">Study site</td><td>0.800</td><td>0.894</td><td>3</td><td></td></tr></tbody></table>

SNP84 in exon 9 lay later than females with CC (β = 0.022 ± 0.009, t<sub>74</sub> = 2.371, P = 0.020; Table 5b). We found no statistical support for relationships between SNPs and fledgling numbers (Tables S23–S26).

Discussion

We found relationships between *SERT* SNPs and multiple personality traits, problem-solving performance and breeding parameters. Only one of these relationships (between SNP36 in exon 6 and incorrect attempts in females) remained significant after adjusting for multiple comparisons, but several were associated with substantial effect sizes (R<sup>2</sup> = 0.03–0.30). These relationships warrant further investigation, given the large sample sizes needed to robustly test for associations between behaviors and SNPs, each of which is expected to have a relatively small independent effect on behavioral or cognitive variation<sup>13–15</sup>.

Past studies have suggested that *SERT* polymorphisms might be involved in variation in neophobia, boldness and exploratory behavior, which are often interrelated traits<sup>37,40</sup>. We found some additional support for this hypothesis in the current study. With respect to neophobia behavior, the strongest gene-behavior association involved synonymous SNP36 in exon 6, which explained 4.8% of the variance in return latencies during novel object, but not baseline, trials. Heterozygous females (CT) were quicker to return to the nest box during novel object trials than homozygotes, suggesting that heterozygous females are less neophobic. However, SNP36 was not related to return latencies during baseline trials (boldness following disturbance at the nest), suggesting that SNP36 may be specifically related to neophobia, rather than affecting boldness-related responses across contexts (i.e. both neophobia and boldness following disturbance). The largest SNP-behavior association for exploratory behavior was with promoter region SNP438, which explained 3.5% of the behavioral variance, with heterozygotes (CT) having higher exploration scores than homozygotes. Our sample sizes for neophobia and exploratory behavior were similar to those in previous studies on great tits<sup>37,40,41</sup>. We thus had similar statistical power to detect a given effect size, and also found effect sizes similar to those reported by past studies<sup>37,40,41</sup>. However, past studies that found statistically significant relationships involved fewer SNPs, and thus lesser corrections for multiple comparisons. In sum, research on relationships between *SERT* polymorphisms, neophobia
<table><caption>**Table 5.** General generic and dominance effect LMMs showing the relationship between, (a) SNP170 in exon 3, and (b) SNP84 in exon 9 and log-transformed laying date in female great tits. <sup>a</sup>P values from posthoc Tukey test. P=0.613 for the CC-CG contrast. <sup>b</sup>P values from posthoc Tukey test. P=0.864 for the CT-TT contrast. N=117 observations on 92 females; SNP170: 77 CC, 11 CG, 3 GG; SNP84: 20 CC, 47 CT, 24 TT.</caption><thead><tr><th colspan="9">(a) SNP170 exon 3</th></tr><tr><th colspan="2">General generic model</th><th colspan="2">R<sup>2</sup>m=0.017</th><th colspan="2">R<sup>2</sup>c=0.653</th><th colspan="2"></th></tr></thead><tbody><tr><th colspan="2">Fixed effects</th><th>t</th><th>d.f</th><th>P</th><th>F</th><th>d.f</th><th>P</th></tr><tr><td colspan="2">Intercept</td><td>4.66±0.039</td><td>118.5</td><td>1.069</td><td>0.003</td><td></td><td></td></tr><tr><td colspan="2">SNP170 CG</td><td>-0.039±0.019</td><td>-1.966</td><td>92.7</td><td>0.131<sup>a</sup></td><td></td><td></td></tr><tr><td colspan="2">SNP170 CC</td><td>-0.051±0.022</td><td>2.285</td><td>58.5</td><td>0.065<sup>a</sup></td><td>2.712</td><td>2.113</td><td>0.070</td></tr><tr><th colspan="2">Random effects</th><th>Variance</th><th>SD</th><th>N</th><th colspan="2"></th><th colspan="2"></th></tr><tr><td colspan="2">Individual</td><td>&lt;0.001</td><td>&lt;0.001</td><td>91</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Study site</td><td>&lt;0.001</td><td>0.013</td><td>4</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Year</td><td>0.003</td><td>0.054</td><td>2</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Residual</td><td>0.002</td><td>0.041</td><td>117</td><td colspan="2"></td><td colspan="2"></td></tr><tr><th colspan="2">Dominance effect</th><th colspan="2">R<sup>2</sup>m=0.014</th><th colspan="2">R<sup>2</sup>c=0.654</th><th colspan="2"></th></tr><tr><th colspan="2">Fixed Effects</th><th>t</th><th>d.f</th><th>P</th><th>F</th><th>d.f</th><th>P</th></tr><tr><td colspan="2">Intercept</td><td>4.661±0.039</td><td>118.6</td><td>1.075</td><td>0.003</td><td></td><td></td></tr><tr><td colspan="2">SNP170 (CG/CC vs GG)</td><td>0.041±0.019</td><td>2.134</td><td>113.6</td><td>0.034</td><td></td><td></td></tr><tr><th colspan="2">Random effects</th><th>Variance</th><th>SD</th><th>N</th><th colspan="2"></th><th colspan="2"></th></tr><tr><td colspan="2">Individual</td><td>&lt;0.001</td><td>&lt;0.001</td><td>91</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Study site</td><td>&lt;0.001</td><td>0.014</td><td>4</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Year</td><td>0.002</td><td>0.054</td><td>2</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Residual</td><td>0.002</td><td>0.041</td><td>117</td><td colspan="2"></td><td colspan="2"></td></tr><tr><th colspan="9">(b) SNP84 exon 9</th></tr><tr><th colspan="2">General generic model</th><th colspan="2">R<sup>2</sup>m=0.016</th><th colspan="2">R<sup>2</sup>c=0.508</th><th colspan="2"></th></tr><tr><th colspan="2">Fixed effects</th><th>t</th><th>d.f</th><th>P</th><th>F</th><th>d.f</th><th>P</th></tr><tr><td colspan="2">Intercept</td><td>4.646±0.038</td><td>120.0</td><td>1.12</td><td>0.002</td><td></td><td></td></tr><tr><td colspan="2">SNP84 CT</td><td>-0.024±0.010</td><td>-2.330</td><td>77.8</td><td>0.057<sup>b</sup></td><td></td><td></td></tr><tr><td colspan="2">SNP84 TT</td><td>-0.019±0.011</td><td>-1.664</td><td>80.3</td><td>0.225<sup>b</sup></td><td>2.916</td><td>2.73</td><td>0.060</td></tr><tr><th colspan="2">Random effects</th><th>Variance</th><th>SD</th><th>N</th><th colspan="2"></th><th colspan="2"></th></tr><tr><td colspan="2">Individual</td><td>&lt;0.001</td><td>0.008</td><td>91</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Study site</td><td>&lt;0.001</td><td>0.010</td><td>4</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Year</td><td>0.002</td><td>0.052</td><td>2</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Residual</td><td>0.002</td><td>0.040</td><td>117</td><td colspan="2"></td><td colspan="2"></td></tr><tr><th colspan="2">Dominance model</th><th colspan="2">R<sup>2</sup>m=0.016</th><th colspan="2">R<sup>2</sup>c=0.509</th><th colspan="2"></th></tr><tr><th colspan="2">Fixed effects</th><th>t</th><th>d.f</th><th>P</th><th>F</th><th>d.f</th><th>P</th></tr><tr><td colspan="2">Intercept</td><td>4.646±0.038</td><td>120.1</td><td>1.12</td><td>0.002</td><td></td><td></td></tr><tr><td colspan="2">SNP84 (CT/TT vs CC)</td><td>0.022±0.009</td><td>2.371</td><td>73.9</td><td>0.020</td><td></td><td></td></tr><tr><th colspan="2">Random effects</th><th>Variance</th><th>SD</th><th>N</th><th colspan="2"></th><th colspan="2"></th></tr><tr><td colspan="2">Individual</td><td>&lt;0.001</td><td>0.007</td><td>91</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Study site</td><td>&lt;0.001</td><td>0.010</td><td>4</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Year</td><td>0.002</td><td>0.052</td><td>2</td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Residual</td><td>0.002</td><td>0.040</td><td>117</td><td colspan="2"></td><td colspan="2"></td></tr></tbody></table>

and exploratory behavior in great tits suggests no consistent sources of genetic control involving *SERT*, but rather that different *SERT* SNPs may be associated with behavioral variation in different populations. Also, despite interrelationships between neophobia, boldness and exploratory behavior, these traits were not related to the same SNPs, likely reflecting the fact that complex behavioral traits can be influenced by many different environmental and genetic effects<sup>13–16</sup>.

Our study represents a first attempt to characterize associations between problem solving and *SERT* polymorphisms in birds, and we found some evidence for such associations in the case of nonsynonymous SNP226 in exon 1 and synonymous SNP84 in exon 9. Nonsynonymous SNP226 in exon 1 was associated with 8% percent of the variance in female problem-solving success and 14% percent of the variance in female incorrect attempts during problem-solving. Females with rare genotype AA at SNP226 directed fewer incorrect attempts at the barrier and succeeded more on the test compared to those with genotype AT or TT. A higher rate of incorrect attempts directed at the barrier during the obstacle removal test (persistent pecking at the barrier) is associated with lower rates of success<sup>61</sup>. Therefore, SNP226 could affect problem-solving success by modulating the rate of incorrect attempts. Indeed, when including both incorrect attempts and latency to contact the test apparatus in
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Relationships between SERT SNPs and lay date (days since January 1). Mean values and standard errors (error bars) for groups were extracted using R package emmeans and back transformed from the log scale; (a) exon 3 SNP170, *N* (observations, females) = CC: 99, 77; CG: 13, 11; GG: 5, 3; (b) exon 9 SNP84. *N* (observations, females) = CC: 25, 20; CT: 61, 47; TT: 31, 24.</figcaption></figure>

statistical models, the effect of SNP226 on female problem-solving success was no longer statistically significant. This suggests that variance in both of these behavioral patterns could contribute to explaining the relationship between SNP226 and problem-solving success, although there was not a statistically significant relationship between SNP226 and latency to contact the test apparatus. The specific aspect of cognitive function which is reflected by incorrect attempts remains unclear, whereas latency to contact could be affected by both motivation and neophobia.

SNP226 is a nonsynonymous variant, and could thus affect serotonergic signaling, cognition and behavior during problem solving through direct effects on SERT protein structure and function. The E26D amino acid change caused by SNP226 is located within the intracellular N terminus of the SERT protein (Fig. S1). Extensive research has focused on the role of transmembrane helices and their connecting intracellular and extracellular loops in mediating substrate recognition and transport dynamics of transmembrane transport proteins, with relatively little work involving the function of transport protein tails (i.e. the N and C termini)⁸². However, emerging evidence suggests that the N and C termini can also play important roles in transporter function⁸², with research specifically indicating that the N terminus is involved in regulating transformational changes in SERT during substrate exchange⁸³,⁸⁴. However, SNAP2 analysis could not reliably predict functionality for the change in amino acid sequence from glutamic acid to aspartic acid caused by SNP226. Why relationships between SNP226, incorrect attempts and problem-solving success were present in females, but not males, also remains to be elucidated, but note that the low number of males with genotype AA at SNP226 made the effect of this SNP difficult to assess in males.

In addition to SNP226, synonymous SNP84 in exon 9 was also associated with problem-solving success, and there was some evidence for an effect of this SNP on problem-solving in both females and males. Specifically, 29% of the variance in male problem-solving and 9% of the variance in female problem-solving was associated with this SNP, although this effect was statistically nonsignificant in females at $α = 0.05$. SNP84 was not significantly related to the other behavioral patterns that we measured during the obstacle removal test, and might thus affect problem-solving through effects on traits not quantified through our protocol.

Synonymous SNP36 in exon 6 and synonymous SNP187 in exon 1 were also associated with female incorrect attempts directed at the barrier, although neither of these SNPs was significantly related to problem-solving success. The relationship between incorrect attempts on the barrier and SNP36 was associated with the largest effect size ($R^2 = 0.30$) found in our study, and remained significant after accounting for multiple comparisons. As was also the case for neophobia behavior, females heterozygous at SNP36 differed in behavior from homozygotes, pecking the barrier more. Thus, one might expect that SNP36 could mediate an association between incorrect attempts and neophobia behavior. However, we did not find a statistically significant association between incorrect attempts and neophobia behavior. The relationship between SNP187 in exon 1 and incorrect attempts was associated with a smaller, but non-negligible, effect size ($R^2 = 0.108$). A different nucleotide change at position 187 in exon 1 (A/T instead of C/T) was previously related to hissing behavior during nest defense in female great tits⁴¹. However, recent work in our study population showed no relationship between SNP187 and either hissing behavior or female-female aggression³⁸. Thus, whether SNPs at this locus have significant behavioral effects across multiple traits remains equivocal. In our dataset, linkage between SNP226 and SNP187 in exon 1 made the effects of the two difficult to disentangle.

We hypothesized that genetic differences between birds could explain reduced neophobia and enhanced problem-solving among individuals occupying highly disturbed territories. However, we found no association between SERT SNPs and distance from roads/paths. Also, neophobia and problem-solving were not significantly associated. Rather, distance from roads/paths and SERT SNPs had independent effects on neophobia (for SNP36 exon 6) and problem-solving (for SNP226 exon 1 and SNP84 exon 9). Thus, multiple environmental and genetic effects may combine to determine behavior and problem-solving. Genotype-by-environment interactions also remain possible, but we had an insufficient sample size to test this possibility.

SNP170 in exon 3 and SNP84 in exon 9 showed associations with egg laying date, suggesting that SERT polymorphisms could affect this reproductive trait. However, these relationships had relatively small effect sizes (explaining ~2% of the variance), and there was no evidence for a relationship between SNPs and fledgling
numbers. Egg laying date in great tits can be important to matching timing of the nestling period with peak caterpillar abundance, and can affect reproductive success<sup>41,85</sup>. One past study in great tits reported that two different *SERT* SNPs in exon 8 (not present in our population) correlated with laying date and nestlings hatched<sup>41</sup>. Relationships between *SERT* polymorphisms and laying date could be mediated through multiple pathways, including effects of serotonin on feeding behavior and endocrine systems, both of which affect the time at which individuals attain breeding condition<sup>41</sup>.

Besides SNP226 in exon 1, the SNPs for which we found associations with behavioral and fitness traits were synonymous SNPs. Synonymous SNPs do not affect the amino acid sequence of the SERT protein, and thus appear unlikely to affect SERT's structure and enzymatic function. However, recent research suggests that synonymous SNPs can affect the expression level, structure and function of proteins, through mechanisms including effects on transcription, splicing and folding of proteins<sup>86</sup>. Alternatively, synonymous SNPs may be linked to functional variants. Synonymous SNPs are as frequently linked to disease in humans as are nonsynonymous SNPs<sup>86</sup>, and in birds nonsynonymous and synonymous SNPs have been linked to behavioral variation at similar rates<sup>35,36</sup>. In great tits, all of the past studies on *SERT*<sup>37,40,41</sup> (and *DRD4*<sup>23,27-29,36</sup>) polymorphisms that have linked SNPs to behavioral variation have involved synonymous SNPs. These past results and our own findings suggest that synonymous SNPs may be linked to a nontrivial amount of genetic variation in behavioral traits.

## Conclusions

Our study revealed associations between *SERT* polymorphisms, personality traits, problem-solving and laying date. Associations were mostly non-significant after corrections for multiple comparisons, but some were associated with substantial effect sizes. Relationships between *SERT*, behavior and cognition in birds warrant further investigation, given the complex genetic control of behavior and large sample sizes required to robustly test for effects of any one genetic variant. Our study serves as a basis for future research, including meta-analyses aimed at synergizing evidence for effects of *SERT* on behavior, cognition and fitness.

## Data availability

Data will be available through the Dryad Digital Repository<sup>87</sup>.

Received: 25 May 2021; Accepted: 2 November 2021

Published online: 20 December 2021

## References

1. Dingemanse, N. J. & Wolf, M. Recent models for adaptive personality differences: A review. *Phil. Trans. R. Soc. B* **365**, 3947-3958. https://doi.org/10.1098/rstb.2010.0221 (2010).

2. Wolf, M., van Doorn, G., Leimar, O. & Weissing, F. J. Life-history trade-offs favour the evolution of animal personalities. *Nature* **447**, 581–584. https://doi.org/10.1038/nature05835 (2007).

3. Dingemanse, N. J., Both, C., Drent, P. J. & Tinbergen, J. M. Fitness consequences of avian personalities in a fluctuating environment. *Proc. R. Soc. B.* **271**, 847–852. https://doi.org/10.1098/rspb.2004.2680 (2004).

4. Sih, A. & Bell, A. M. Insights for behavioral ecology from behavioral syndromes. *Adv. Study Behav.* **38**, 227–281. https://doi.org/10.1016/S0065-3454(08)00005-3 (2008).

5. Sih, A., Bell, A. M. & Johnson, J. C. Behavioral syndromes: An ecological and evolutionary overview. *Trends Ecol. Evol.* **19**, 372–378. https://doi.org/10.1016/j.tree.2004.04.009 (2004).

6. Drent, P. J., van Oers, K. & van Noordwijk, A. J. Realized heritability of personalities in the great tit (*Parus major*). *Proc. R. Soc. B.* **270**, 45–51. https://doi.org/10.1098/rspb.2002.2168 (2003).

7. Sol, D., Griffin, A. S., Bartomeus, I. & Boyce, H. Exploring or avoiding novel food resources? The novelty conflict in an invasive bird. *PLoS ONE* **6**, e19535. https://doi.org/10.1371/journal.pone.0019535 (2011).

8. Dammhahn, M., Mazza, V., Schirmer, A., Göttliche, C. & Eccard, J. C. Of city and village mice: Behavioural adjustments of striped field mice to urban environments. *Sci. Rep.* **10**, 13056. https://doi.org/10.1038/s41598-020-69998-6 (2020).

9. Sih, A. & Del Giudice, M. Linking behavioural syndromes and cognition: A behavioural ecological perspective. *Phil. Trans. R. Soc. B* **367**, 2762–2772. https://doi.org/10.1098/rstb.2012.0216 (2012).

10. Stoewe, M. & Kotrschal, K. Behavioural phenotypes may determine whether social context facilitates or delays novel object exploration in ravens (*Corvus corax*). *J. Ornithol.* **148**, S179-S184. https://doi.org/10.1007/s10336-007-0145-1 (2007).

11. Guillette, L. M., Reddon, A. R., Hoeschele, M. & Sturdy, C. B. Sometimes slower is better: Slow-exploring birds are more sensitive to changes in a vocal discrimination task. *Proc. R. Soc. B* **278**, 767-773. https://doi.org/10.1098/rspb.2010.1669 (2011).

12. Dochtermann, N. A., Schwab, T. & Sih, A. The contribution of additive genetic variation to personality variation: Heritability of personality. *Proc. R. Soc. B* **282**, 20142201. https://doi.org/10.1098/rspb.2014.2201 (2015).

13. Van Oers, K., De Jong, G., Van Noordwijk, A. J., Kempenaers, B. & Drent, P. J. Contribution of genetics to the study of animal personalities: A review of case studies. *Behaviour* **142**, 1185-1206. https://doi.org/10.1163/156853905774539364 (2005).

14. Van Oers, K. & Mueller, J. C. Evolutionary genomics of animal personality. *Phil. Trans. R. Soc. B* **365**, 3991–4000. https://doi.org/10.1098/rstb.2010.0178 (2010).

15. Croston, R., Branch, C. L., Kozlovsky, D. Y., Dukas, R. & Pravosudov, V. V. Heritability and the evolution of cognitive traits. *Behav. Ecol.* **26**, 1447-1459. https://doi.org/10.1093/beheco/arv088 (2015).

16. Quinn, J. L., Cole, E. F., Reed, T. E. & Morand-Ferron, J. Environmental and genetic determinants of innovativeness in a natural population of birds. *Phil. Trans. R. Soc. B* **371**, 20150184. https://doi.org/10.1098/rstb.2015.0184 (2016).

17. Evans, J., Boudreau, K. & Hyman, J. Behavioural syndromes in urban and rural populations of song sparrows. *Ethology* **116**, 588-595. https://doi.org/10.1111/j.1439-0310.2010.01771.x (2010).

18. Bókony, V., Kulcsár, A., Tóth, Z. & Liker, A. Personality traits and behavioral syndromes in differently urbanized populations of house sparrows (*Passer domesticus*). *PLoS ONE* **7**, 36639. https://doi.org/10.1371/journal.pone.0036639 (2007).

19. Charmantier, A., Deyeyrier, V., Lambrechts, M., Perret, S. & Grégoire, A. Urbanization is associated with divergence in pace-of-life in great tits. *Front. Ecol. Evol.* https://doi.org/10.3389/fevo.2017.00053 (2017).

20. Isaksson, C., Rodewald, A. D. & Gil, D. Editorial: Behavioural and ecological consequences of urban life in birds. *Front. Ecol. Evol.* https://doi.org/10.3389/fevo.2018.00050 (2018).

21. Audet, J.-N., Ducatez, S. & Lefebvre, L. The town bird and the country bird: Problem solving and immunocompetence vary with urbanization. *Behav. Ecol.* **27**, 637-644. https://doi.org/10.1093/beheco/arv201 (2016).
22. Miranda, A. C., Schielzeth, H., Sonntag, T. & Partecke, J. Urbanization and its effects on personality traits: A result of microevolution or phenotypic plasticity. *Glob. Change Biol.* **19**, 2634–2644. https://doi.org/10.1111/gcb.12258 (2013).

23. Riyahi, S., Björklund, M., Mateos-Gonzalez, F. & Senar, J. C. Personality and urbanization: Behavioural traits and DRD4 SNP830 polymorphisms in great tits in Barcelona city. *J. Ethol.* **35**, 101–108. https://doi.org/10.1007/s10164-016-0496-2 (2017).

24. Schinka, J. A., Letsch, E. A. & Crawford, F. C. DRD4 and novelty seeking: Results of meta-analyses. *Am. J. Med. Genet.* **114**, 643–648. https://doi.org/10.1002/ajmg.10649 (2002).

25. Chen, C. S., Burton, M., Greenberger, E. & Dmitrieva, J. Population migration and the variation of Dopamine D4 Receptor (DRD4) allele frequencies around the globe. *Evol. Hum. Behav.* **20**, 309–324. https://doi.org/10.1016/S1090-5138(99)00015-X (1999).

26. Shimada, M. K. et al. Polymorphism in the second intron of dopamine receptor D4 gene in humans and apes. *Biochem. Biophys. Res. Commun.* **316**, 1186–1190. https://doi.org/10.1016/j.bbcrc.2004.03.006 (2004).

27. Fidler, A. E. et al. Drd4 gene polymorphisms are associated with personality variation in a passerine bird. *Proc. R. Soc. B.* **274**, 1685–1691. https://doi.org/10.1098/rspb.2007.0337 (2007).

28. Mueller, J. C. et al. Haplotype structure, adaptive history and associations with exploratory behaviour of the DRD4 gene region in four great tit (*Parus major*) populations. *Mol. Ecol.* **22**, 2797–2809. https://doi.org/10.1111/mec.12282 (2013).

29. Korsten, P. et al. Association between DRD4 gene polymorphism and personality variation in great tits: A test across four wild populations. *Mol. Ecol.* **19**, 832–843. https://doi.org/10.1111/j.1365-294X.2009.04518.x (2010).

30. Jiang, W., Shang, S. & Su, Y. Genetic influences on insight problem solving: The role of catechol-O-methyltransferase polymorphisms. *Front. Psychol.* **6**, 1569. https://doi.org/10.3389/fpsyg.2015.01569 (2015).

31. Hopkins, W. et al. Genetic influences on receptive joint attention in chimpanzees (*Pan troglodytes*). *Sci. Rep.* **4**, 3774. https://doi.org/10.1038/srep03774 (2014).

32. Fitzpatrick, M. J. et al. Candidate genes for behavioural ecology. *Trends Ecol. Evol.* **20**, 96–104. https://doi.org/10.1016/j.tree.2004.11.017 (2005).

33. Munafo, M. R., Brown, S. M. & Harkless, K. C. Serotonin transporter (5-HTTLPR) genotype and amygdala activation: A meta-analysis. *Biol. Psychiatry* **63**, 852–857. https://doi.org/10.1016/j.biopsych.2007.08.016 (2008).

34. Staes, N. et al. Serotonin receptor 1A variation is associated with anxiety and agonistic behavior in chimpanzees. *Mol. Biol. Evol.* **36**, 1418–1429. https://doi.org/10.1093/molbev/msz061 (2019).

35. Mueller, J. C. et al. Behaviour-related DRD4 polymorphisms in invasive bird populations. *Mol. Ecol.* **23**, 2876–2885. https://doi.org/10.1111/mec.12763 (2014).

36. Timm, K., Tilgar, V. & Saag, P. DRD4 gene polymorphism in great tits: Gender-specific association with behavioural variation in the wild. *Behav. Ecol. Sociobiol.* **69**, 729–735. https://doi.org/10.1007/s00265-015-1887-z (2015).

37. Riyahi, S., Sánchez-Delgado, M., Calafell, F., Monk, D. & Senar, J. C. Combined epigenetic and intraspecific variation of the DRD4 and SERT genes influence novelty seeking behaviour in great tit *Parus major*. *Epigenetics* **10**, 516–525. https://doi.org/10.1080/15592294.2015.1046027 (2015).

38. Holtmann, B. et al. Population differentiation and behavioural association of the two 'personality' genes DRD4 and SERT in dunnocks (*Prunella modularis*). *Mol. Ecol.* **25**, 706–722. https://doi.org/10.1111/mec.13514 (2016).

39. Krause, E. T., Kjaer, J. B., Lüders, C. & van Phi, L. A polymorphism in the 5'-flanking region of the serotonin transporter (5-HTT) gene affects fear-related behaviors of adult domestic chickens. *Behav. Brain Res.* **14**, 92–96. https://doi.org/10.1016/j.bbr.2017.04.051 (2017).

40. Timm, K., van Oers, K. & Tilgar, V. SERT gene polymorphisms are associated with risk-taking behaviour and breeding parameters in wild great tits. *J. Exp. Biol.* **221**, jeb171595. https://doi.org/10.1242/jeb.171595 (2018).

41. Timm, K., Koosa, K. & Tilgar, V. The serotonin transporter gene could play a role in anti-predator behaviour in a forest passerine. *J. Ethol.* **37**, 221–227. https://doi.org/10.1007/s10164-019-00593-7 (2019).

42. Berger, M., Gray, J. A. & Roth, B. L. The expanded biology of serotonin. *Annu. Rev. Med.* **60**, 355–366. https://doi.org/10.1146/annurev.med.60.042307.110802 (2009).

43. Lesch, K. P. & Merschdorf, U. Impulsivity, aggression, and serotonin: A molecular psychobiological perspective. *Behav. Sci. Law* **18**, 581–604 (2000).

44. Duke, A. A., Bègue, L., Bell, R. & Eisenlohr-Moul, T. Revisiting the serotonin-aggression relation in humans: A meta-analysis. *Psychol. Bull.* **139**, 1148–1172. https://doi.org/10.1037/a0031544 (2013).

45. Ferrari, P. F., Palanza, P., Parmigiani, S., de Almeida, R. M. & Miczek, K. A. Serotonin and aggressive behavior in rodents and nonhuman primates: Predispositions and plasticity. *Eur. J. Pharmacol.* **526**, 259–273. https://doi.org/10.1016/j.ejphar.2005.10.002 (2005).

46. Bacqué-Cazenave, J. et al. Serotonin in animal cognition and behavior. *Int. J. Mol. Sci.* **21**, 1649. https://doi.org/10.3390/ijms2105169 (2020).

47. Walker, S. C. et al. Selective prefrontal serotonin depletion impairs acquisition of a detour-reaching task. *Eur. J. Neurosci.* **23**, 3119–3123. https://doi.org/10.1111/j.1460-9568.2006.04826.x (2006).

48. Cools, R., Roberts, A. C. & Robbins, T. W. Serotoninergic regulation of emotional and behavioural control processes. *Trends Cogn. Sci.* **12**, 31–40. https://doi.org/10.1016/j.tics.2007.10.011 (2008).

49. Rudnick, G. & Sandtner, W. Serotonin transport in the 21st century. *J. Gen. Physiol.* **151**, 1248–1264. https://doi.org/10.1085/jgp.201812066 (2018).

50. Lesch, K. P. et al. Association of anxiety-related traits with a polymorphism in the serotonin transporter gene regulatory region. *Science* **274**, 1527–1531. https://doi.org/10.1126/science.274.5292.1527 (1996).

51. Sen, S., Burmeister, M. & Ghosh, D. Meta-analysis of the association between a serotonin transporter promoter polymorphism (5-HTTLPR) and anxiety-related personality traits. *Am. J. Med. Genet.* **127**, 85–89. https://doi.org/10.1002/ajmg.b.20158 (2004).

52. Karg, K., Burmeister, M., Shedden, K. & Sen, S. The serotonin transporter promoter variant (5-HTTLPR), stress, and depression meta-analysis revisited: Evidence of genetic moderation. *Arch. Gen. Psychiatry* **68**, 444–454. https://doi.org/10.1001/archgenpsychotherapy.2010.189 (2011).

53. Beversdorf, D. Q. et al. Influence of serotonin transporter *SLC6A4* genotype on the effect of psychosocial stress on cognitive performance: An exploratory pilot study. *Cogn. Behav. Neurol.* **31**, 79–85. https://doi.org/10.1097/WNN.0000000000000153 (2018).

54. Canli, T. & Lesch, P.-K. Long story short: The serotonin transporter in emotion regulation and social cognition. *Nat. Neurosci.* **10**, 1103–1109. https://doi.org/10.1038/nn1964 (2007).

55. Jarrell, H. et al. Polymorphisms in the serotonin reuptake transporter gene modify the consequences of social status on metabolic health in female rhesus monkeys. *Physiol. Behav.* **93**, 807–819. https://doi.org/10.1016/j.physbeh.2007.11.042 (2008).

56. Bennett, A. et al. Early experience and serotonin transporter gene variation interact to influence primate CNS function. *Mol. Psychiatry* **7**, 118–122. https://doi.org/10.1038/smp.4000949 (2002).

57. Golebiowska, J. et al. Serotonin transporter deficiency alters socioemotional ultrasonic communication in rats. *Sci. Rep.* **9**, 20283. https://doi.org/10.1038/s41598-019-56629-y (2019).

58. Thys, B. et al. The serotonin transporter gene and female personality variation in a free-living passerine. *Sci. Rep.* **11**, 8577. https://doi.org/10.1038/s41598-021-88225-4 (2021).

59. Audet, J.-N. et al. Divergence in problem-solving skills is associated with differential expression of glutamate receptors in wild finches. *Sci. Adv.* **4**, eaa06369. https://doi.org/10.1126/sciadv.aao6369 (2018).
60. Grunst, A. S., Grunst, M. L., Pinxten, R. & Eens, M. Personality and plasticity in neophobia levels vary with anthropogenic disturbance but not toxic metal exposure in urban great tits. *Sci. Total Environ.* **656**, 997–1009. https://doi.org/10.1016/j.scitotenv.2018.11.383 (2019).

61. Grunst, A. S., Grunst, M. L., Pinxten, R. & Eens, M. Sources of individual variation in problem-solving performance in urban great tits (*Parus major*): Exploring effects of metal pollution, urban disturbance and personality. *Sci. Tot. Environ.* **749**, 141436. https://doi.org/10.1016/j.scitotenv.2020.141436 (2020).

62. Thys, B. et al. The female perspective of personality in a wild songbird: Repeatable aggressiveness relates to exploration behavior. *Sci. Rep.* **7**, 7656. https://doi.org/10.1038/s41598-017-08001-1 (2017).

63. Grunst, A. S. et al. An important personality trait varies with blood and plumage metal concentrations in a free-living songbird. *Environ. Sci. Technol.* **53**, 10487–10496. https://doi.org/10.1021/acs.est.9b03548 (2019).

64. Grunst, A. S. et al. Variation in personality traits across a metal pollution gradient in a free-living songbird. *Sci. Total Environ.* **630**, 668–678. https://doi.org/10.1016/j.scitotenv.2018.02.19 (2018).

65. Laucht, M. et al. Interaction between the 5-HTTLPR serotonin transporter polymorphism and environmental adversity for mood and anxiety psychopathology: Evidence from a high-risk community sample of young adults. *Int. J. Neuropharmacol.* **12**, 737–747. https://doi.org/10.1017/S1461145708009875 (2009).

66. Wang, Z. et al. Genome-wide gene by lead exposure interaction analysis identifies UNC5D as a candidate gene for neurodevelopment. *Environ. Health* **16**, 81. https://doi.org/10.1186/s12940-017-0288-3 (2017).

67. Grunst, A. S., Grunst, M. L., Pinxten, R. & Eens, M. Proximity to roads, but not exposure to metal pollution, is associated with accelerated developmental telomere shortening in nestling great tits. *Environ. Pollut.* **256**, 113373. https://doi.org/10.1016/j.envpol.2019.113373 (2019).

68. Dingemanse, N. J. et al. Repeatability and heritability of exploratory behaviour in great tits from the wild. *Anim. Behav.* **64**, 929–937. https://doi.org/10.1006/anbe.2002.2006 (2002).

69. Solé, X. et al. SNPStats: A web tool for the analysis of association studies. *Bioinformatics* **22**, 1928–1929. https://doi.org/10.1093/bioinformatics/bti283 (2005).

70. Hecht, M., Bromberg, Y. & Rost, B. Better prediction of functional effects for sequence variants from VarI-SIG 2014: Identification and annotation of genetic variants in the context of structure, function and disease. *BMC Genom.* **16**, S1. https://doi.org/10.1186/1471-2164-16-S8-S1 (2015).

71. Choi, Y. & Chan, A. PROVEAN web server: A tool to predict the functional effect of amino acid substitutions and indels. *Bioinformatics* **31**, 2745–2747. https://doi.org/10.1093/bioinformatics/btv195 (2015).

72. Omasits, Ü., Ahrens, C. H., Müller, S. & Wollsccheid, B. Protter: Interactive protein feature visualization and integration with experimental proteomic data. *Bioinformatics* **30**(6), 884–886. https://doi.org/10.1093/bioinformatics/btt607 (2014).

73. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria (2019). URL https://www.R-project.org/.

74. Bates, D., Maechler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. *J. Stat. Softw.* **67**, 1–48. https://doi.org/10.18637/jss.v067.i01 (2014).

75. Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest Package: Tests in linear mixed effects models. *J. Stat. Softw.* **82**, 1–26. https://doi.org/10.18637/jss.v082.i13 (2017).

76. Stoffel, M. A., Nakagawa, S. & Schielzeth, H. rptR: repeatability estimation and variance decomposition by generalized linear mixed-effects models. *Methods Ecol. Evol.* **8**, 1639–1644. https://doi.org/10.1111/2041-210X.12797 (2017).

77. Harrison, X. A. Using observation-level random effects to model overdispersion in count data in ecology and evolution. *PeerJ* **2**, e616. https://doi.org/10.7717/peerj.616 (2014).

78. Lenth, R. emmeans: Estimated Marginal Means, aka Least-Squares Means. R package version 1.4.3.01 (2019). https://CRAN.R-project.org/package=emmeans.

79. Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. *J. R. Stat. Soc. B* **57**, 289–300. https://doi.org/10.2307/2346101 (1995).

80. Nakagawa, S. & Schielzeth, H. A general and simple method for obtaining R² from generalized linear mixed-effects models. *Methods Ecol. Evol.* **4**, 133–142. https://doi.org/10.1111/j.2041-210x.2012.00261.x (2013).

81. Lüdecke, D., Makowski, D., Waggoner, P. & Patil, I. performance: Assessment of Regression Models Performance. R package version 0.4.6 (2020). https://CRAN.R-project.org/package=performance.

82. Hartig, F. DHARMA: Residual Diagnostics for Hierarchical (Multi-Level/Mixed) Regression Models. R package version 0.2.6 (2019). https://CRAN.R-project.org/package=DHARMA.

83. Mikros, E. & Diallinas, G. Tales of tails in transporters. *Open Biol.* **9**, 190083. https://doi.org/10.1098/rsob.190083 (2019).

84. Kern, C. et al. The N terminus specifies the switch between transporter modes of the human serotonin transporter. *J. Biol. Chem.* **292**, 3603–3613. https://doi.org/10.1074/jbc.M116.771360 (2017).

85. Visser, M. E., Van Noordwijk, A. J., Tinbergen, J. M. & Lessells, C. M. Warmer springs lead to mistimed reproduction in great tits (*Parus major*). *Proc. R. Soc. B* **265**, 1867–1870. https://doi.org/10.1098/rspb.1998.0514 (1998).

86. Hunt, R., Sauna, Z. E., Ambudkar, S. V., Gottesman, M. M. & Kimchi-Sarfaty, C. Silent (Synonymous) SNPs: Should we care about them? In *Single Nucleotide Polymorphisms Methods in Molecular Biology (Methods and Protocols)* Vol. 578 (ed. Komar, A.) (Humana Press, 2009). https://doi.org/10.1007/978-1-60327-411-1_2.

87. Grunst, A. S., Grunst, M. L. & Staes, N., Bert, T., Pinxten, R., Eens, M. Data for: Serotonin Transporter (SERT) Polymorphisms, Personality and Problem-Solving in Urban Great Tits. (Dryad Digital Repository, 2021).

## Acknowledgements

We thank Geert Eens and Peter Scheys for providing support during fieldwork, Jasmijn Daans for help with genetic work in the laboratory, and members of the Behavioural Ecology and Ecophysiology group at the University of Antwerp for intellectual support. This study received financial support from the FWO Flanders (to ASG, Grant ID: 1.2135.17N) and European Commission (Marie Skłodowska-Curie fellowship 799667 to MLG). The authors declare that they have no conflicts of interest.

## Author contributions

All authors contributed to study design. A.S.G., M.L.G. and B.T. performed behavioral testing and other fieldwork. A.S.G., M.L.G. and N.S. performed laboratory work. N.S. designed primers, aligned sequences and performed SNAP2 analysis. A.S.G. and M.L.G. analyzed data. A.S.G. wrote the first draft of the manuscript, M.L.G. made a substantial revision, and all authors contributed to reviewing and finalizing the paper.

## Competing interests

The authors declare no competing interests.
**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-03466-7.

**Correspondence** and requests for materials should be addressed to A.S.G.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

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

© The Author(s) 2021