Clock gene polymorphism and scheduling of migration: a geolocator study of the barn swallow Hirundo rustica.

Circannual rhythms often rely on endogenous seasonal photoperiodic timers involving 'clock' genes, and Clock gene polymorphism has been associated to variation in phenology in some bird species. In the long-distance migratory barn swallow Hirundo rustica, individuals bearing the rare Clock allele with the largest number of C-terminal polyglutamine repeats found in this species (Q8) show a delayed reproduction and moult later. We explored the association between Clock polymorphism and migration scheduling, as gauged by light-level geolocators, in two barn swallow populations (Switzerland; Po Plain, Italy). Genetic polymorphism was low: 91% of the 64 individuals tracked year-round were Q7/Q7 homozygotes. We compared the phenology of the rare genotypes with the phenotypic distribution of Q7/Q7 homozygotes within each population. In Switzerland, compared to Q7/Q7, two Q6/Q7 males departed earlier from the wintering grounds and arrived earlier to their colony in spring, while a single Q7/Q8 female was delayed for both phenophases. On the other hand, in the Po Plain, three Q6/Q7 individuals had a similar phenology compared to Q7/Q7. The Swiss data are suggestive for a role of genetic polymorphism at a candidate phenological gene in shaping migration traits, and support the idea that Clock polymorphism underlies phenological variation in birds.

The timing of seasonal activities has major fitness consequences and is subjected to intense natural selection12. In species experiencing seasonal peaks of resource availability, typically occurring at medium-high latitudes, proper matching of life-cycle events to such resource pulses is fundamental for achieving a high fitness. The timing of seasonal activities is often set by genetically controlled endogenous circadian and circannual rhythms3, which are synchronized with ecological conditions through external drivers4, among which photoperiod plays a pivotal role5. This seasonal photoperiodic timer involves ‘circadian clock’ genes, which are responsible for the onset and setting of circadian and circannual rhythmicity678. The circadian clock relies on an auto-regulated negative feedback loop, the ‘core circadian oscillator’ (CCO)910. In birds and mammals, the CLOCK transcription factor, encoded by the Circadian Locomotor Output Cycles Kaput (Clock) gene, plays a central role within the CCO, both by acting as a transcriptional activator of the CCO itself and by activating the transcription of several output genes1112. The C-terminal domain of the Clock gene contains a series of polyglutamine residues (Poly-Q) which may vary in length among species as well as among individuals within populations131415161718.

Latitudinal increase in the number of Poly-Q repeats across populations has been suggested to reflect adaptation to local regimes of circannual photoperiodic variation141619. Similarly, Poly-Q polymorphism may control variation in photoperiodic responses among individuals from the same population, causing variation in phenology at the individual level15172021. Though the evidence for an association between Clock polymorphism and timing of seasonal events at the population level is scant, two recent studies conducted on a long-distance migratory bird, the barn swallow Hirundo rustica, have shown that a very rare Clock genotype, the heterozygous Q7/Q8 (c.a. 1.5% of the population), which contains the largest Poly-Q stretch observed in this species (Q8), has a delayed timing compared to the most frequent (94%) Q7/Q7 genotype1720. These findings are consistent with the idea that larger Poly-Q stretches are associated with a delayed phenology15172021, and suggest that the Q7/Q8 genotype shows a delayed annual routine. Specifically, the first study, carried out on an Italian barn swallow population, showed that the onset of laying by Q7/Q8 females was delayed by 13 days relative to Q7/Q7 counterparts, corresponding to ca. 0.7 s.d. later than the mean of Q7/Q7 females17. In the second study, carried out on moulting barn swallows in the sub-Saharan wintering areas, wing feather moult of three Q7/Q8 individuals was found to be significantly delayed, by ca. 8/9 days, than Q7/Q7 birds20. In both studies, the phenology of the other rare genotype Q6/Q7 (4.8%)17, with a smaller Poly-Q stretch, did not differ from that of the common genotype Q7/Q71720.

To date, the analysis of individual variation in scheduling of migration events across the entire annual cycle in small-sized birds has been hampered by the technical difficulties of tracking large numbers of individuals from the same population. Here we exploit data retrieved from 64 individual barn swallows tracked along the entire annual migration cycle using miniaturized light-level geolocators22 to explore the association between Clock genotype polymorphism and scheduling of migration events in two barn swallow populations breeding in Southern Switzerland (Magadino) and Northern Italy (Po Plain). We expected that rare individuals bearing the heterozygous Q7/Q8 genotype should show a delayed scheduling of migration and wintering events compared to the phenotypic distribution of phenological traits of Q7/Q7 birds, complementing previous evidence concerning delayed timing of breeding and moult.

Results and Discussion

Most (90.6%) of the 64 individuals tracked for an entire annual cycle (Switzerland, n = 26; Po Plain, n = 38) were homozygous Q7/Q7, as expected. In the Swiss population, we detected two Q6/Q7 males and a single Q7/Q8 female, while in the Po Plain we detected three Q6/Q7 individuals (two males, one female) and no Q7/Q8 birds.

The migration and wintering phenology of Q7/Q7 birds showed a broad overlap between the two study populations, with the single exception of spring arrival date to the breeding colony, that was ca. 10 days later among Swiss birds compared to Po Plain ones (Fig. 1; Supplementary Table S1, Fig. S1; see also Liechti et al.22).

Figure 1

Frequency distribution of phenotypic values of the Q7/Q7 individuals and phenology of the Q6/Q7 and Q7/Q8 individuals in the Swiss population.

(a) Departure from the breeding colony; (b) arrival to the wintering area; (c) departure from the wintering area; (d) arrival to the breeding colony. The horizontal boxes above histograms show the mean and 95% non-parametric bootstrap confidence interval (CI) of the Q7/Q7 phenotypic distribution, while the horizontal lines show the standard deviation. Dots show the phenotypic values of the rare genotypes (open dot = Q6/Q7; filled dot = Q7/Q8). Histograms and dots refer to actual data, while asterisks above dots denote that the Q6/Q7 or Q7/Q8 values were outside the 95% CI of a given phenological trait for centred values (removing any between-year and between-sex variation in phenology; see Supplementary Table S2 and Methods for details).

In the Swiss population, the single Q7/Q8 female showed a consistently delayed phenology and was outside the 95% confidence interval (CI) compared to the phenotypic distribution of Q7/Q7 for all traits (Fig. 1), both when actual or within-group centered values (thus removing any between-year and between-sex variation in phenology) were considered, the single exception being departure date from the breeding colony for actual values (Fig. 1a; Supplementary Table S2). On the other hand, the two Q6/Q7 males were consistently delayed and outside the 95% CI compared to the phenotypic distribution of Q7/Q7 individuals for departure from the breeding colony, but were consistently early and outside the 95% CI compared to Q7/Q7 for traits related to pre-breeding migration (Fig. 1c,d; Supplementary Table S2). They were instead not consistently advanced for arrival date to wintering area (Fig 1b; Supplementary Table S2). Results were similar if the analyses were separated by sex (Supplementary Table S3) and if we accounted for the effect of individual variation in wintering latitude and longitude (stationary wintering positions, SWP; see Methods) on arrival date to the breeding colony (Supplementary Table S2).

In the Po Plain, the lack of Q7/Q8 individuals allowed only a partial evaluation of the association between Clock polymorphism and phenology. Although being somewhat delayed in some cases (Supplementary Fig. S1, Table S2), the scheduling of migration and wintering events of the three Q6/Q7 was not consistently different from the phenotypic distribution of Q7/Q7 birds, with the exception of colony departure, that was delayed in all three individuals compared to the 95% CI of Q7/Q7 birds, similarly to the Swiss population (Supplementary Fig. S1, Table S2).

It may be worth noting that the two Swiss Q6/Q7 males and the Q7/Q8 female were also outside the CIs (advanced and delayed, respectively) of the Q7/Q7 Po Plain birds for departure from the wintering area and arrival to the breeding colony (cf. Fig. 1 and Supplementary Fig. S1; Table S2).

For the Swiss Q7/Q8 female, delayed arrival could lead to a delayed reproduction17. This female started egg laying on May 28th, 2012 (day 148), but unfortunately we have no other breeding data from geolocator females in the same year (nor from the same female in the previous year) for comparison. However, data from Q7/Q7 geolocator females from the previous year (spring 2011) revealed a mean first egg laying date that was ca. 20 days earlier (day 126, CI 120 to 136, n = 5) than the Q7/Q8 female. We are confident that the delayed laying of the Q7/Q8 female compared to Q7/Q7 females in the previous year was robust to between-year differences in egg laying dates: in fact, the Q7/Q8 female was a delayed phenodeviant for arrival date (Table 1; Fig. 1), and over the entire sample of geolocator females (both Swiss and Po Plain individuals) there was a strict association between arrival date and breeding date (r = 0.90, n = 16; Zr = 1.47, CI 0.93 to 2.01; Fig. 2). On the other hand, the two Swiss Q6/Q7 males did not have a clearly advanced breeding date compared to Q7/Q7 birds (details not shown for brevity). Indeed, the correlation between breeding date (start of egg laying by the social mate) and arrival date in geolocator males was weaker than in females (r = 0.31, n = 37; Zr = 0.32, CI −0.02 to 0.66), suggesting that timing of reproduction is mainly under female control17.

Figure 2

Arrival date predicts onset of reproduction among female barn swallows equipped with geolocators (day 1 = January 1st).

The fitted regression equation is: first egg laying date = 0.88 (0.12 s.e.m.) × arrival date + 30.48 (12.51 s.e.m.) (R2 = 0.80). The filled dot denotes the single Q7/Q8 female from the Swiss population.

These findings, based on the first individual year-round tracking data for the barn swallow, complement previous evidence that the rare Q7/Q8 barn swallow genotype has a delayed timing1720. Importantly, this evidence originates from the analysis of a further independent dataset, as no individuals from the Swiss population had been included in previous studies1720. The delay of the Swiss Q7/Q8 individual was most pronounced for phenological variables related to timing of spring migration (departure from the wintering area and arrival date to the breeding colony) (Fig. 1).

Intriguingly, the two Swiss Q6/Q7 males showed an apparently advanced timing of spring migration compared to Q7/Q7 birds in the same population. In previous studies of the barn swallow, no statistically significant differences in timing emerged between Q6/Q7 and Q7/Q7 birds1720. The Po Plain data show instead an inconsistent tendency towards a delayed pre-breeding migration phenology of the three Q6/Q7 birds compared to Q7/Q7 ones. Albeit not conclusively due to the small sample size of the ‘rare’ genotypes, the picture emerging from this and our earlier work may suggest that ‘clock’ genes differently affect phenology in different barn swallow populations, as previously documented by an analysis of the songbird genus Junco18. The most likely explanation for such among-population variation in ‘clock’ gene-phenology associations has been proposed to be variation in genetic background among different populations of the same species, with possibly other genes related to phenology variably interacting with Clock in different geographical populations18.

The proximate mechanisms leading to an association between length of Poly-Q stretches and phenology remain unknown, but may involve a different response to photoperiod or to other astronomical cues while in Africa. Despite photoperiodic changes around the Equator are minor, birds are able to detect them and tune the timing of their seasonal activities (including timing of moult) accordingly23. Moreover, it has recently been suggested that tropical birds can use variation in time of sunrise and sunset to adjust the timing of life-history events along the annual cycle24. Our findings suggest that the previously documented reproductive delay in Q7/Q8 females was likely to be mainly a consequence of carry-over effects of events occurring during wintering (such as moult) and at the onset of spring migration. Elucidating whether the Clock-phenology association causally stems from variable photoperiodic responses while in the African wintering grounds, triggering variable onset of migration activities with carry-over effects on timing of reproduction, would be a challenge for future research.

The fact that the association between Clock genotype and timing of autumn migration was somewhat less clear may be due to the difficulties in identifying the onset of autumn migration. Geolocator data allow only an estimate of the timing of departure from the breeding colony, which may not be indicative of variation in onset of post-breeding migration. In fact, the timing of departure from the breeding colony may be heavily affected by involvement in reproductive activities (e.g. nest failure, laying of multiple clutches and extent of post-fledging parental care) besides timing of reproduction. In addition, departures are very concentrated and apparently highly synchronized with stochastic events such as rainfall bouts, which may vary locally between study colonies22.

Why such rare genotypes persist in barn swallow populations is open to speculation. It is known that occasional spring cold spells, occurring after arrival of barn swallows to their breeding areas, can cause large mortality2526. Clearly, the rare, delayed Q7/Q8 individuals may accrue a survival advantage under such intense natural selection episodes, and this may contribute to maintain Clock polymorphism in barn swallow populations. On the other hand, early arriving Q6/Q7 males may accrue a high fitness advantage by mating earlier if environmental conditions are favourable, even if the association between timing of spring arrival and reproduction was not as strong as in females in the study populations, since timing of reproduction predicts seasonal fitness in barn swallows2527. Finally, we wish to point out that both here and in our previous studies we focused on the association between Clock gene polymorphism and phenology, but future studies will need to address whether the huge phenological variance observed among the most common Q7/Q7 genotype is mostly due to environmental effects or whether it can be explained by genetic variation at other phenological candidate genes.

To conclude, capitalizing on a unique dataset of individual year-round tracking data of a small long-distance migratory bird obtained by miniaturized light-level geolocators, the Swiss data are suggestive for a role of genetic polymorphism at a candidate phenological gene in shaping variation in migration traits. Despite the very low genetic variation and the consequently small sample of ‘rare’ genotypes, our findings foster the idea that Clock gene polymorphism may underlie phenological variation in avian species.

Methods

Study species, study sites and field procedures

The barn swallow is a small (ca. 20 g) Holarctic semi-colonial insectivorous passerine bird. European populations are largely migratory28. Birds move southwards in September-October, and winter in sub-Saharan Africa, where they undergo their annual complete plumage moult. During wintering, birds are mostly stationary22. They then leave for spring migration in February-March and reach their breeding sites mid-March to late May22. Females lay one to three clutches (3–7 eggs) per season28. In Europe, barn swallows mostly breed in small colonies settled in rural buildings (farms, cowsheds, stables)29.

The study was carried out at three study sites, one in southern Switzerland, in an Alpine valley floor (Magadino, 46°09’ N, 8°55’ E, 211 m a.s.l.) and two in the Po Plain (northern Italy; Lombardy, 45°19’ N, 9°40’ E, 60 m a.s.l.; Piedmont, 45°33’ N, 8°44’ E, 160 m a.s.l.) (details in Scandolara et al.2930 and Liechti et al.22). We captured barn swallows at several colonies during two breeding seasons (2010 and 2011). A very large sample (640 individuals; 2010, n = 310; 2011, n = 330) of adult birds (i.e. at least 1 year old) was equipped with miniaturized light-level geolocators weighting on average 0.68–0.77 g (3.74–4.14% of body mass)30. During the subsequent breeding season we retrieved 124 geolocators (2010, n = 78; 2011, n = 46)22. We collected a small blood sample from birds that successfully returned to their breeding colony with the geolocator. Blood was collected from the brachial vein into heparinized capillary tubes, stored in a cool bag and subsequently kept at −20/−80 °C. Birds were marked with individual colour rings upon capture and were assigned to their nests by direct observation1730. Laying date (Julian date of first egg laying) was determined for most individual by means of frequent inspections of nest content (details in Scandolara et al.30).

All procedures were performed in accordance with the Swiss and Italian regulations concerning scientific investigations on bird species in the wild, and approved by the Office fédéral de l’environnement (OFEV, Division Espèces, écosystèmes, paysages; Switzerland) (permit n. F044-0799), by Regione Lombardia (auth. n. 329 issued on January 21, 2009 and n. 2141 issued on March 9, 2011) and by the Provincia di Novara (auth. n. 905 issued on March 21, 2011).

Geolocator data analysis

Deriving accurate geographical positions from light-level data obtained from geolocators constitutes a major methodological challenge31. The full details and rationale of the analyses of barn swallow geolocator data are reported in Liechti et al.22. Hereafter we briefly summarize the interpretation of the phenological traits of interest for defining individual variation in the timing of the annual routine. First, based on standardized analyses of daily changes in sunrise and sunset time derived from light-level data (latitudinal estimates were not calculated during migration because the barn swallow migration periods coincide with the period around the Equinoxes, when data derived from geolocators do not allow to accurately calculate latitudes2231), we assigned each daily position to a stationary or movement (migration) period. We then calculated the following phenological traits at the individual level: departure from the breeding colony (day of departure from the breeding colony, determined by means of visual inspection of variation in daily light-level profiles); arrival to the wintering area (day of the first stationary period south of the Sahara, i.e. south of 23.5°N); departure from the wintering area (day of the last stationary period south of the Sahara); arrival to the breeding colony (date of arrival to the breeding colony, determined by means of visual inspection of variation in daily light-level profiles). Among the 124 individuals that returned with the geolocator (see above), we could obtain complete phenological information (i.e. we could obtain phenological data until spring arrival date to the breeding colony) for 68 individuals only (2010, n = 37; 2011, n = 31), due to total or partial failure of the geolocators (e.g. battery failure), which did not allow recovering the entire migration path (see Table 1 in Liechti et al.22). For each of these complete tracks we determined the stationary wintering position (SWP) as the centre of density (mode) of longitude and latitude of all the daily positions, taking into account the stationary periods south of the Sahara22.

Genetic analyses

Genomic DNA extraction was performed by alkaline lysis using 6 μl of blood in 100 μl of a 50 mM NaOH at 100 °C for 20 minutes. DNA was quantified using a spectrophotometer and diluted to a final concentration of 50–100 ng/μl. PCR amplification was performed using primers designed on the barn swallow genomic sequence: forward primer (5′-labelled with 6-FAM dye) [6FAM]GGGACAGGTGGTGACAGCTTATC and reverse primer CTGCTGATGGTCCTGCTGACT (Sigma-Aldrich)17; PCR fragments were screened for length polymorphism at ClkpolyQcds by resolution and detection on a conventional DNA sequencing machine using a commercial fragment analysis service (Macrogen Inc., Seoul, Republic of Korea)17. The reliability of molecular data for the barn swallow ClkpolyQcds locus has previously been confirmed by independently repeating the fragment analysis of several individuals1720. Alleles were named according to the number of glutamine repeats predicted in the mature protein as Q6 to Q8, after sequencing the most common allele (Q7, 112 bp long, coding for a CLOCK protein containing a stretch of 7 glutamine repeats)32. We genotyped the vast majority (64 out of 68 individuals; 2010, n = 34; 2011, n = 30) of the birds for which we obtained complete migration tracks (missing genotypes were due to failed genotyping or missing samples).

Statistical analyses

Since the genotype-phenology associations may differ between populations18, tests of associations between Clock polymorphism and phenology were conducted within population. First, we calculated the 95% non-parametric bootstrap CIs (BCa method, based on 5000 replicates; R-package boot, ver. 1.3–933) of each phenological trait of Q7/Q7 individuals. We then checked whether the phenological values of the Q6/Q7 and Q7/Q8 individuals were included or not within the CIs. If a phenological value was outside the Q7/Q7 CI of a given trait, that individual was considered phenodeviant for that trait (i.e. it had a phenotype that significantly deviated from the phenotypic distribution of the reference genotype). Male and female data were pooled for analyses as there were no statistically significant differences between the sexes for most phenological traits, with the exception of departure from the breeding colony, that was 2 days later on average in males compared to females22. In addition, to completely remove between-year and between-sex differences22, we adopted a within-group centring procedure, subtracting from a given value the mean of its corresponding group (i.e. 2010 males, 2010 females, 2011 males, 2011 females). CIs were calculated for centred values as well. Data from the Lombardy and Piedmont study areas were pooled because both were in the same geographical area (Po Plain) and because of the small sample of rare genotypes at each study site (see Results and Discussion).

Additional Information

How to cite this article: Bazzi, G. et al. Clock gene polymorphism and scheduling of migration: a geolocator study of the barn swallow Hirundo rustica. Sci. Rep. 5, 12443; doi: 10.1038/srep12443 (2015).