Gene expression in the brain of a migratory songbird during breeding and migration
© Boss et al. 2016
Received: 8 September 2015
Accepted: 28 January 2016
Published: 15 February 2016
We still have limited knowledge about the underlying genetic mechanisms that enable migrating species of birds to navigate the globe. Here we make an attempt to get insight into the genetic architecture controlling this complex innate behaviour. We contrast the gene expression profiles of two closely related songbird subspecies with divergent migratory phenotypes. In addition to comparing differences in migratory strategy we include a temporal component and contrast patterns between breeding adults and autumn migrating juvenile birds of both subspecies. The two willow warbler subspecies, Phylloscopus trochilus trochilus and P. t. acredula, are remarkably similar both in phenotype and genotype and have a narrow contact zone in central Scandinavia. Here we used a microarray gene chip representing 23,136 expressed sequence tags (ESTs) from the zebra finch Taeniopygia guttata to identify mRNA level differences in willow warbler brain tissue in relation to subspecies and season.
Out of the 22,109 EST probe sets that remained after filtering poorly binding probes, we found 11,898 (51.8 %) probe sets that could be reliably and uniquely matched to a total of 6,758 orthologous zebra finch genes. The two subspecies showed very similar levels of gene expression with less than 0.1 % of the probe sets being significantly differentially expressed. In contrast, 3,045 (13.8 %) probe sets were found to be differently regulated between samples collected from breeding adults and autumn migrating juvenile birds. The genes found to be differentially expressed between seasons appeared to be enriched for functional roles in neuronal firing and neuronal synapse formation.
Our results show that only few genes are differentially expressed between the subspecies. This suggests that the different migration strategies of the subspecies might be governed by few genes, or that the expression patterns of those genes are time-structured or tissue-specific in ways, which our approach fails to uncover. Our findings will be useful in the planning of new experiments designed to unravel the genes involved in the migratory program of birds.
KeywordsMigration Microarray Birds Gene expression Behavior Calcium transport Synapse Phylloscopus trochilus
Bird migration is a complex behavior that involves a suite of traits evolved to optimize the movements between areas for reproduction, staging and winter survival . Therefore, the combined behavioral, morphological and physiological traits associated with flight and navigation are characterized as migration phenotypes . Comparative analyses of migratory and sedentary populations within passerine bird species, or between closely related species that differ in their migratory behaviors, have demonstrated strikingly parallel patterns of adaptive evolution on the phenotype level [3, 4]. The coordination of a set of complex behavioral adaptations has been referred to as ‘behavioral syndrome’ , a concept that has been used to discuss the evolution of migration [6, 7].
Quantitative genetic approaches have demonstrated heritable variation and strong selection for several key components of the migratory phenotype, such as migratory propensity, timing, distance and migratory orientation in birds [8, 9], migratory propensity in fish  and wing morphology in insects  and birds [12, 13]. In many species of long distance migratory songbirds, the core components of migration (timing, distance and direction) likely have a genetic basis since juveniles typically migrate on their own without influence of experienced adults . Hence, we may expect differentially migrating bird populations to differ in genes encoding for several adaptations directly or indirectly linked to modulation of the migratory phenotype.
Whereas quantitative genetic studies demonstrate the presence or absence of within-species genetic variation associated with migratory traits, they do not provide data required for testing whether the same or different genes are involved in the migratory phenotypes of different species . The work of identifying the genes that influence the behavioral adaptations associated with migration is still in its infancy and among the so far proposed candidates there is none that has consistently shown to be influencing migratory traits across species. The most extensive study comparing differential expression patterns between phenotypes used microarray analyses of monarch butterflies Danaus plexippus to compare summer and migratory phenotypes. This study identified 40 differentially expressed genes between seasons using brain tissue . Two genes of particular interest were the clock gene vrille and the locomotion-relevant TβH (tyramine beta hydroxylase) gene as their presumed functions directly link to migratory components (timing and movement respectively). A recent genome re-sequencing study of monarch butterflies identified a gene (collagen IV) influencing flight muscle functioning with the migratory alleles resulting in improved flight efficiency . One study on European blackcaps Sylvia atricapilla  analyzed six genes with functions suggesting potential association to the migratory program. A microsatellite polymorphism in the 3′-UTR of the gene encoding Adenylate Cyclase Activating Polypeptide 1 (ADCYAP1) was found to explain a small but significant proportion of migratory restlessness within populations, as well as variation in migratory distance across populations .
In Northern Europe, two geographically separated subspecies of the willow warbler (Phylloscopus trochilus) have been used extensively to study differential migration behavior with regard to migration direction and wintering locations in Africa [18–20]. The subspecies P. t. trochilus, breeding in southern Scandinavia, Britain and western parts of Central Europe, migrates towards the south-west in autumn to wintering grounds in West Africa. In contrast, the subspecies P. t. acredula, breeds in northern Scandinavia, Finland and Russia, migrates south-south-east to eastern and southern Africa [18–21]. Whereas the contact zone between the two species east of the Baltic Sea is complex, with different positions of cline centres and cline widths for size, coloration and migration strategy, they form a well-defined and 200 km wide contact zone in central Scandinavia [21, 22]. The subspecies are inseparable with regard to neutral genetic markers (e.g., mitochondrial DNA, microsatellites and AFLPs; ) and breeding habitat  suggesting that the few divergent genomic regions are maintained by selection associated with their different migration phenotypes .
To identify molecular components of the brain that are potentially involved in migration direction behaviour, we used a microarray chip designed from expressed sequence tags (ESTs) in brain tissue of zebra finches  to compare genome-wide expression of mRNA transcripts in brain tissue collected from two willow warbler subspecies with different migratory strategies during two crucial life history stages, i.e. breeding in adults as well as the early phase of migration in juvenile birds. Our rational behind choosing the brain as the target for the analyses is because it is the tissue where migration-related decision-making gets integrated. We decided to use the whole brains rather than separate regions because we are lacking precise information and characterisation of focal brain regions that process migration-related decision making genes. Based on the low degree of genetic differentiation between the subspecies [21, 24] and the different challenges the birds are facing during breeding and migration, we expect to find more difference in gene expression when comparing breeding and autumn samples within the subspecies than when comparing the subspecies within seasons. If genes involved in migration-direction decision-making are included in our microarray chip, we expected these to be differentially expressed when comparing the autumn samples of the two subspecies. However, given the lack of knowledge as to whether expression of such genes is restricted to the migration period, we included the overall contrast between the subspecies (within seasons) in our analysis.
Fieldwork and sample collection
RNA preparation and analysis
Nine months post collection of the samples (May and June 2009) brains were removed from the RNAlater buffer and homogenized in 1 ml QIA-zol Lysis Reagent per 100 mg tissue (Qiagen) using a TissueLyser (Qiagen). The samples were processed in batches of four individuals (one from each group) using the RNeasy Lipid Tissue Mini Kit (Qiagen) according to the manufacturer’s instructions. Samples were quality-controlled by electrophoresis on formaldehyde agarose gels and confirmed to have an acceptable 260 nm/280 nm absorption ratio (all values between 1.95 and 2.05) using an Ultraspec 3000 spectrophotometer (Pharmacia).
For this project, we had 29 Lund-zf2 microarray chips available and therefore RNA samples with the highest RNA Integrity Numbers  were selected (8 samples from each group except for P. t. acredula in autumn for which only 5 samples passed the RNA quality criteria). The Lund-zf2 array is a custom Affymetrix array with probe sets representing 23,136 zebra finch ESTs expressed in brain tissue. Each EST probe set is represented by eleven 25 bases long oligonucleotides (probes) and together the probe sets represent about 15,800 non-redundant gene models, for details see . The microarrays were hybridized and analyzed at the Swegene Center for Integrative Biology at Lund University (SCIBLU genomics, http://www.lth.se/sciblu) using standard procedures (for details, see ). Across-species array hybridization is challenging as it may result in measured expression differences that are methodological artefacts that needs to be controlled for. Sequence mismatches between the probe and the target may for example underestimate the expression level. We have taken several measures to minimize such biases: First, all probes that had failed to hybridize against genomic DNA from willow warblers were removed , see details below. Second, this microarray chip has previously been thoroughly evaluated for across-species hybridization in the common whitethroat Sylvia communis [28–30], a species positioned at the same phylogenetic distance to the zebra finch as the willow warbler. In these studies, the patterns of sex-biased gene expression were very similar for the zebra finch (within-species array hybridization) and the whitethroat (across-species array hybridization), both regarding the specific genes found being differentially regulated between the sexes and the magnitude of sex-biased expression.
All the birds were sexed by standard molecular methods  and the collected wings of the birds captured during the autumn were carefully checked to confirm that they had been correctly aged . Of the 29 analyzed samples (P. t. trochilus, 8 summer and 8 autumn; P. t. acredula, 5 summer and 8 autumn) we excluded two samples from the autumn birds (one from each subspecies) that were subsequently shown not to be juvenile males (one juvenile female and one adult male). All statistical analyses of the microarray dataset were carried out using R version 2.15.1 (R Core Development Team 2012). Quality of the array results was evaluated using the affyPLM package. A background correction for each probe was performed using the robust multi-array average (RMA) method as implemented in the affy package. Subsequently, the probe signals were normalized using RMA normalization. We evaluated the normalizations by inspecting the MA, NUSE and RLE plots before and after RMA normalization. The expression data was further evaluated by hierarchical clustering (average clustering method using Euclidean distances), and principal component analysis to identify possible outlier arrays (data not shown). One array (from a breeding acredula male) was marked as an outlier and was therefore removed from the subsequent analyses. Data from the remaining 26 arrays was re-normalized followed by a new round of quality control. Consequently, our reported analyses were carried out on data from eight breeding P.t. trochilus males, four breeding P. t. acredula males, seven migrating P. t. trochilus males, and seven migrating P. t. acredula males.
Since low probe hybridization signals can result not only from a low gene expression level but also from a poor sequence match of willow warbler transcripts to the probes designed from zebra finch ESTs, we made use of a previous genomic DNA hybridization data using the Lund-zf array  to determine cut-off parameters. This experiment contained data for two arrays each of which had been hybridized with pooled DNA from eight P. t. trochilus and eight P. t. acredula birds, respectively. Probes that had a signal at least two standard deviations above the mean background expression level on at least one of these arrays were flagged as present (i.e., the willow warbler sequence matches the zebra finch probe), while probes yielding signals below two standard deviations were flagged as absent (i.e. poor cross-species sequence match). Out of the 23,139 probe sets on the zf-2 array, 1,030 (4.4 %) had less than eight functional probes (Affymetrix standard cut-off; ) after subtraction of poorly matching probes and were excluded from the downstream analyses (Additional file 1: Figure S1). For the remaining 22,109 probe sets, confirmed to show significant hybridization signals to willow warbler DNA, we calculated the mean intensities for probes in each probe set as estimates of gene expression levels.
Probe sets were linked to genes using blastn 2.2.24+ in the toolset of NCBIs standalone blast (http://blast.ncbi.nlm.nih.gov/Blast.cgi?CMD=Web&PAGE_TYPE=BlastNews#1). The database (zebra finch gene set version 3.2.4) was downloaded from Ensembl (www.ensembl.org) using the biomart tool and included 18,581 annotated genes together with potentially uncharacterized untranslated region 1,000 bp up- and downstream of the first and last exons, respectively. Flanking sequences were included to allow identification of EST sequences corresponding to untranslated regions of gene transcripts. These blast results were then evaluated using an in-house Python script (supplementary file) in order to remove probes with more than one possible hit (second hit within an e-value range of 103 of the first) and alignments shorter than 50 bp. The four sample groups were compared in order to identify gene expression differences associated with the effects of season and subspecies. Two additional comparisons were made by combining sample groups according to subspecies or season. To avoid artefacts associated with identification of false positives due to multiple testing (i.e. erroneously rejecting the null hypothesis of no difference between groups) we used the false discovery rate correction (FDR) with a statistical cut-off value of 5 % .
Significant enrichment of differentially regulated genes associated with gene ontology terms (GO terms) was determined using the CORNA R package . We selected the annotated and significantly differentially regulated genes between breading and migratory season (1,640) for gene ontology analysis. We did not carry out GO term analyses with respect to subspecies/migratory strategy as the numbers for significantly differentially regulated genes between subspecies were too low for meaningful subsequent analyses. Gene descriptions were gathered from Ensembl and statistical significance of the enrichment was evaluated using a Fisher’s exact test with correction for false discovery rate (FDR). The inbound relationships between GO terms were examined using the web tool Revigo (http://revigo.irb.hr/) and the method SimRel . SimRel calculates a differentiation score, pairing all GO terms against each other. Revigo was used to create scatterplots of the physiological differences between GO terms of biological processes below p < 0.05.
Estimated number of differentially expressed genes
Gene expression differences between breeding adults and juvenile migrants
Top list of the 20 annotated genes that differed most significantly between samples from breeding adult and juvenile migrating male willow warblers
Ensembl gene ID
ATPase, Na+/K+ transporting, alpha 1 polypeptide
human immunodeficiency virus type I enhancer binding protein 2
fascin actin-bundling protein 1
cAMP-regulated phosphoprotein, 21 kDa
family with sequence similarity 49, member A
CREB3 regulatory factor
collagen, type I, alpha 2
bone morphogenetic protein/retinoic acid inducible neural-specific 3
SRSF protein kinase 2
complement component 8, gamma polypeptide
peripheral myelin protein 2
NHL repeat containing E3 ubiquitin protein ligase 1
DCN1, defective in cullin neddylation 1, domain containing 5
transducin-like enhancer of split 4
pleckstrin homology domain containing, family A member 5
Taeniopygia guttata parvalbumin-like
Gene expression differences between subspecies
The overall clustering of the expression profiles (Fig. 2) showed that the difference in the number of differentially expressed genes between subspecies was small. In line with this, within-season comparisons of gene expression differences between the subspecies did not provide support for any differences between subspecies either in summer or in autumn (Fig. 4). In order to increase the power of the comparison, we compared probe data for all birds from each subspecies regardless of season. Expression data from 11 P. t. acredula birds were thus compared to data from 15 P. t. trochilus individuals. This comparison identified 14 differentially expressed probe sets (FDR adjusted p-value <0.05, 0.6 % of the total) (Fig. 4). Including the interaction term subspecies x season only marginally changed the results (12 differentially expressed probe sets of which 10 overlapped with the 14 detected without the interaction). Seven of the 14 significant genes did not have Ensembl annotations or were annotated as uncharacterized. Nonetheless, a significant and annotated gene of potential interest is the PPP3CA gene, which encodes a calcium-regulated (via Calmodulin) protein phosphatase. PPP3CA is the only gene that was shown to be significantly changed in comparisons of both the seasons and subspecies, perhaps due to a role of calcium-related signalling in both contexts. PPP3CA mRNA levels were significantly higher in migrating birds compared to breeding birds, and higher in the southern subspecies, P. t. trochilus, compared to the northern subspecies, P. t. acredula.
Differentially regulated genes are enriched in specific functional gene classes
Lists of the ten most significant results from the gene ontology enrichment analyses between breeding adult and juvenile migrating willow warblers
positive regulation of gene expression
calcium ion transport
calcium ion transmembrane transport
structural constituent of ribosome
RNA polymerase II regulatory region sequence-specific DNA binding
protein serine/threonine kinase activity
protein domain specific binding
regulation of mRNA stability
Expression differences of genes included in the ontology term referred to as “calcium ion transport”
ATP2B2 ATPase, Ca++ transporting, plasma membrane 2
ATPase, Ca++ transporting, plasma membrane 1
Transient receptor potential cation channel, subfamily M, member 7
Sodium channel, voltage-dependent, type III, alpha subunit
Calcium channel, voltage-dependent, N type, alpha-1B subunit
Inositol 1,4,5-triphosphate receptor, type 1
Protein phosphatase 3, catalytic subunit, alpha isozyme
Calcium channel, voltage-dependent, gamma-2 subunit
Glutamate receptor, ionotropic, N-methyl D-aspartate 1
ATPase, Ca++ transporting, type 2C member 1
Calcium channel, voltage-dependent, gamma-3 subunit
Transient receptor potential cation channel, subfamily C, member 5
The two closely related willow warbler subspecies with clearly different migratory strategy provide a good system to study changes in gene expression associated with migration. As expected, the overall gene expression differed substantially between breeding adults and migrating juveniles, whereas the differences between the subspecies were much smaller. This is in line with the general assumption that differences in the migratory phenotype, such as migratory direction, are likely governed by a limited number of genes. This might also be due to the fact that the expressions of the genes involved in these processes are time-structured or tissue-specific in ways which our approach failed to uncover. The low level of differential gene expression between the subspecies is in line with previous studies examining sequence-based differences [21, 38]. In these analyses, the subspecies were inseparable based on mitochondrial DNA, microsatellites and at >99 % of 500 investigated AFLP markers. Analyzes of sequence differences in the transcriptome  has confirmed the main absence of genomic divergence between the two subspecies of willow warblers, with the exception of three regions (each of a few Mbp) on the putative chromosomes 1, 3 and 5. However as this data was based on pooled RNA samples in combination with a normalization procedure to equalize transcripts across the transcriptome , we could not use it for evaluating expression differences between the subspecies.
Large gene expression differences between breeding adults and migrating juveniles
The highly skewed raw p-value distribution (Fig. 3a) suggested that the actual number of genes involved in seasonal or age dependent differences of expression is higher than those 3,045 determined using the Benjamini and Hochberg false discovery rate . Our data suggest that as many as 22 % of EST probe sets (4,953 out of 22,109) may be differentially expressed between season or age classes. To correctly set the tradeoff between Type I and Type II errors is a well-known problem in all gene expression studies  and it has been suggested that well under 50 % of truly differently expressed genes fail detection even in well-powered studies . Furthermore, this number could also be an underestimation because of the experimental setup where the whole brain was investigated. Differences in gene expression restricted to small regions of the brain might be swamped by mRNA levels from larger brain regions, or possibly due to upregulation in one region and down-regulation in others .
Many of the gene expression differences seen in this study are probably due to the physiological requirements of different seasons and associated life history stages, but still there are alternative explanations that cannot be ruled out. For example, our sampling protocol does not allow discriminating between season and age differences of the population. The analyzed migrating male willow warblers were all in their first year (2–3 months old), while the breeding individuals consisted of males that were one year or older. Thus, any putative age-dependent difference shaping the pattern will be picked up as season differences in our experimental design.
Inspection of the 3,045 probe sets with the lowest p-values (Table 1) reveals a high number of neuron and calcium flux related genes. Regulation of calcium concentrations have been linked to multiple functions within brain cells such as neuronal transmission, excitatory synapses, cell growth, proliferation and synaptic plasticity . One of the significantly differentiating expressed gene, the adenylate cyclase activation polypeptide 1 (ADCYAP1, Additional file 1: Table S1), has previously been associated with migration behavior. It has been shown that the length of the 3′ untranslated region (UTR) on ADCYAP1 alleles partly predicts migratory restlessness within some bird species [17, 43]. Because we find ADCYAP1 mRNA levels to be lower during migration, our study suggests that ADCYAP1 expression might inhibit migratory restlessness. But since the breeding birds were adults and the migrating birds were juveniles on their first migration we cannot disentangle the effects of season from any putatively confounding age effect. The gene ADCYAP1 is highly conserved among birds and mammals, and encodes the protein PACAP that together with its selective receptor PAC1 is integral to anxiety behavior and stress responses in the limbic system of the brain .
The gene ontology enrichment analysis identified the terms “positive regulation of gene expression” and “calcium ion transport” as up-regulated in juvenile migrants compared to adult breeding willow warblers. The up-regulation of genes that increase gene expression agrees with our observed pattern the most differentially expressed genes being up-regulated in migrating juveniles compared to breeding adults.
Our background set included all genes that we could identify as matching the array probes. Since the sample list contains predominantly brain-expressed genes while our background set contains both brain expressed and non-brain expressed genes, one concern might be that typical nervous-tissue related genes become significance just because our samples consisted of brain tissue. However, we only found one of the possible brain-specific GO categories (calcium ion transport) as being significantly enriched; this would not be expected if we were just generally finding brain GO classes. Our data do not allow us to say why calcium ion transport is up-regulated in the brain during migration, but one potential explanation could involve the proliferation of neurons and creation of new synaptic connections needed during migration. Migration flights require the ability to use biological compasses and maps for navigation [2, 45–47], which are expected to be associated with high neuronal activity and a capacity to integrate information from different sensory systems [48, 49]. Another possibility is that the biological processes differing between the two seasons are mainly linked to the breeding period. For example, earlier studies suggest that genes involved in ion transport are repressed during singing behavior .
Small subspecies-dependent gene expression differences
The raw p-value distribution suggests that up to around 1,300 probe sets (6 %) might potentially involve genes differentially expressed in both subspecies, even those that cannot be individually confirmed with statistical confidence. Only 14 genes (0.06 %) were statistically supported (adjusted p < 0.05) to have differential expression between the two subspecies after adjusting for multiple testing. A study investigating subspecies differences between the hooded crow , Corvus corone cornix and the carrion crow, Corvus corone corone, that mainly differ in plumage color, found a much higher number of differentially expressed genes (106). Although no gene overlapped with this study, the cellular functions of neural regulation were differentially expressed in both studies.
The strongest supported probe set was similar to adenylate kinase isoenzyme 5 (AK5) which had a higher expression in P. t. trochilus than in P. t. acredula (FC 1.52). This protein is predicted to sense cellular energy levels . Another significant gene is PPP3CA (higher expression in P. t. acredula than in P. t. trochilus; FC 1.90), which encodes the catalytic subunit of serine/threonine-protein phosphatase (alpha isozyme), a catalytic subunit of the phosphatase calcineurin. PPP3CA was the only gene with significant expression change in both the seasonal (autumn > summer) and subspecies (P. t. trochilus > P. t. acredula) comparisons.
Differentially regulated genes involved in neuronal signaling
The experimental approach used for this first study on gene expression differences between two contrasting migratory strategies in willow warblers is limited by both methodological and conceptual factors, and is not sufficient for identifying all the specific genes involved. Our tentative list of genes differentially expressed in the brain between breeding adults and migrating juveniles as well as between subspecies indicate that specific functional areas and classes might be important for adaptation to migration and associated traits. However, more experiments are needed to confirm this as our seasonal classification criterion is confounded by age. With these limitations in mind, our study suggests an importance of genes involved in establishment of neuron networks (e.g., synapse formation) and neuron signaling (calcium ion transport). Both these functions might potentially play important roles involved in the cognitive/behavioral differences that characterize birds in different seasons, age classes and/or different subspecies. Calcium signaling is involved both in neuron firing and more long-term plasticity changes in neurons that fire repeatedly, leading to calcium-mediated signaling to the nucleus in cells in which calcium levels are elevated during longer times [51, 52]. Based on the present results, it will be possible to focus future studies on which calcium transporter subunits are differentially expressed in different brain regions. In such studies, markers for different classes of neurons should also be studied. Measurements of mRNA levels in different brain regions should also be complemented by in situ hybridization studies to determine whether the differentially expressed genes are predominantly expressed in neurons.
Our study suggests that relatively few genes are involved in shaping the phenotypic differences in migratory strategies between the willow warbler subspecies. However, our snapshot approach (one tissue, summer and autumn) cannot rule out that major differences in gene expression exist in other tissues or at other time points. For example, gene expression differences restricted to small regions of the brain would likely not be detected in analyses of total brain mRNA extracts. Critical expression differences might also be timed to when the birds execute the directed flights, i.e. after sunset, whereas the samples we analyzed were collected in mornings. A tighter coupling of differential gene expression to migration direction choice will require development of experimental migration-related endpoints that can be measured in controlled experiments using birds in captivity where relevant cues (e.g. magnetic field, polarized light, direction of sun set) can be manipulated [53, 54]. Such studies should of course make use of the advancements in next generation sequencing technologies such as RNA-seq which will be more powerful to detect differences in gene expression than our microarray approach.
The sample collection was performed with an ethical permission from Malmö/Lund Committee for Animal Experiment Ethics (no. M22-05). The study was funded by grants from the Swedish research council (SB: 621-2010-5277; SÅ: 621-2007-5930, 621-2010-5584; AW: K2013-66X-14027-10-3; DH: 621-2010-5399) and by The Foundation for Baltic and East European Studies (MG, AW). Support was received from the Centre for Animal Movement Research (CAnMove) financed by a Linnaeus grant (349-2007-8690) from the Swedish Research Council and Lund University.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Alerstam T, Hedenström A, Åkesson S. Long-distance migration: evolution and determinants. Oikos. 2003;103:247–60.View ArticleGoogle Scholar
- Åkesson S, Hedenström A. How migrants get there: Migratory performance and orientation. Bioscience. 2007;57:123–33.View ArticleGoogle Scholar
- Forschler MI, Bairlein F. Morphological shifts of the external flight apparatus across the range of a passerine (Northern Wheatear) with diverging migratory behaviour. Plos One. 2011;6:9.View ArticleGoogle Scholar
- Pulido F. The genetics and evolution of avian migration. Bioscience. 2007;57:165–74.View ArticleGoogle Scholar
- Sih A, Bell AM, Johnson JC, Ziemba RE. Behavioral syndromes: an integrative overview. Quart Rev Biol. 2004;79:241–77.View ArticleGoogle Scholar
- Piersma T, Pérez-Tris J, Mouritsen H, Bauchinger U, Bairlein F. Is there a “migratory syndrome” common to all migrant birds? Bird Hormones and Bird Migrations: Analyzing Hormones in Droppings and Egg Yolks and Assessing Adaptations in Long-Distance Migration. 2005;1046:282–93.Google Scholar
- Dingle H. Animal migration: is there a common migratory syndrome? J Ornithol. 2006;147:212–20.View ArticleGoogle Scholar
- Helbig AJ. Inheritance of migratory direction in a bird species - a cross-breeding experiment with SE-migrating and SW-migrating blackcaps (Sylvia-Atricapilla). Behav Ecol Sociobiol. 1991;28:9–12.View ArticleGoogle Scholar
- Pulido F, Berthold P, van Noordwijk AJ. Frequency of migrants and migratory activity are genetically correlated in a bird population: Evolutionary implications. Proc Natl Acad Sci USA. 1996;93:14642–7.View ArticleGoogle Scholar
- Paez DJ, Brisson-Bonenfant C, Rossignol O, Guderley HE, Bernatchez L, Dodson JJ. Alternative developmental pathways and the propensity to migrate: a case study in the Atlantic salmon. J Evol Biol. 2011;24:245–55.View ArticleGoogle Scholar
- Roff DA, Tucker J, Stirling G, Fairbairn DJ. The evolution of threshold traits: effects of selection on fecundity and correlated response in wing dimorphism in the sand cricket. J EvolBiol. 1999;12:535–46.Google Scholar
- Tarka M, Åkesson M, Beraldi D, Hernandez-Sanchez J, Hasselquist D, Bensch S, et al. A strong quantitative trait locus for wing length on chromosome 2 in a wild population of great reed warblers. Proc R Soc Lond B. 2010;277:2361–9.View ArticleGoogle Scholar
- Tarka M, Åkesson M, Hasselquist D, Hansson B. Intralocus sexual conflict over wing length in a wild migratory bird. Am Nat. 2014;183:62–73.View ArticleGoogle Scholar
- Liedvogel M, Åkesson S, Bensch S. The genetics of migration on the move. Trends Ecol Evol. 2011;26:561–9.View ArticleGoogle Scholar
- Zhu HS, Gegear RJ, Casselman A, Kanginakudru S, Reppert SM. Defining behavioral and molecular differences between summer and migratory monarch butterflies. BMC Biol. 2009;7. doi:10.1186/1741-7007-7-14.
- Zhan S, Zhang W, Niitepold K, Hsu J, Haeger JF, Zalucki MP, et al. The genetics of monarch butterfly migration and warning colouration. Nature. 2014;514:317–21.View ArticleGoogle Scholar
- Mueller JC, Pulido F, Kempenaers B. Identification of a gene associated with avian migratory behaviour. Proc R Soc Lond B. 2011;278:2848–56.View ArticleGoogle Scholar
- Hedenström A, Pettersson J. Migration routes and wintering areas of willow warblers Phylloscopus trochilus (L) ringed in Fennoscandia. Ornis Fenn. 1987;64:137–43.Google Scholar
- Chamberlain CP, Bensch S, Feng X, Åkesson S, Andersson T. Stable isotopes examined across a migratory divide in Scandinavian willow warblers (Phylloscopus trochilus trochilus and Phylloscopus trochilus acredula) reflect their African winter quarters. Proc R Soc Lond B. 2000;267:43–8.View ArticleGoogle Scholar
- Bensch S, Bengtsson G, Åkesson S. Patterns of stable isotope signatures in willow warbler Phylloscopus trochilus feathers collected in Africa. J Avian Biol. 2006;37:323–30.View ArticleGoogle Scholar
- Bensch S, Grahn M, Müller N, Gay L, Åkesson S. Genetic, morphological, and feather isotope variation of migratory willow warblers show gradual divergence in a ring. Mol Ecol. 2009;18:3087–96.View ArticleGoogle Scholar
- Bensch S, Andersson T, Åkesson S. Morphological and molecular variation across a migratory divide in willow warblers, Phylloscopus trochilus. Evolution. 1999;53:1925–35.View ArticleGoogle Scholar
- Larson KW, Liedvogel M, Addison B, Kleven O, Laskemoen T, Lifjeld JT, et al. Allelic variation in a willow warbler genomic region is associated with climate clines. PloS ONE. 2014;9(5):e95252.View ArticleGoogle Scholar
- Lundberg M, Boss J, Canbäck B, Liedvogel M, Larson KW, Grahn M, et al. Characterisation of a transcriptome to find sequence differences between two differentially migrating subspecies of the willow warbler Phylloscopus trochilus. BMC Genomics. 2013;14:11.View ArticleGoogle Scholar
- Naurin S, Bensch S, Hansson B, Johansson T, Clayton DF, Albrekt AS, et al. A microarray for large-scale genomic and transcriptional analyses of the zebra finch (Taeniopygia guttata) and other passerines. Mol Ecol Res. 2008;8:275–81.View ArticleGoogle Scholar
- Bensch S, Lindström A. The age of young willow warblers Phylloscopus trochilus estimated from different stages of post-juvenile moult. Ornis Svecica. 1992;2:23–8.Google Scholar
- Pettersson J, Hasselquist D. Fat deposition and migrationcapacity of robins Erithacus rubecula and goldcrests at Ottenby, Sweden. Ring Migr. 1985;6:66–76.View ArticleGoogle Scholar
- Schroeder A, Mueller O, Stocker S, Salowsky R, Leiber M, Gassmann M, et al. The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. 2006;7:3. doi:10.1186/1471-2199-7-3.
- Naurin S, Hansson B, Hasselquist D, Kim YH, Bensch S. The sex-biased brain: sexual dimorphism in gene expression in two species of songbirds. BMC Genomics. 2011;12:37. http://www.biomedcentral.com/1471-2164/12/37.
- Naurin S, Hasselquist D, Bensch S, Hansson B. Sex-biased gene expression on the avian Z chromosome: highly expressed genes show higher male-biased expression. Plos One. 2012;7:12.View ArticleGoogle Scholar
- Fridolfsson AK, Ellegren H. A simple and universal method for molecular sexing of non-ratite birds. J Avian Biol. 1999;30:116–21.View ArticleGoogle Scholar
- Larson KW, Liedvogel M, Bensch S, Akesson S, Wassenaar LI, Hobson KA. Inferring the ecology of willow warblers during their winter moult by sequential stable isotope analyses of remiges. J Avian Biol. 2013;44:561–6.View ArticleGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate - a practical and powerful approach to multiple testing. J Royal Statist Soc Series B-Method. 1995;57:289–300.Google Scholar
- Wu X, Watson M. CORNA: testing gene lists for regulation by microRNAs. Bioinformatics. 2009;25:832–3.View ArticleGoogle Scholar
- Schlicker A, Domingues FS, Rahnenfuhrer J, Lengauer T. A new measure for functional similarity of gene products based on Gene Ontology. BMC Bioinformatics. 2006;7:302. doi:10.1186/1471-2105-7-302.
- Dean C, Dresbach T. Neuroligins and neurexins: linking cell adhesion, synapse formation and cognitive function. Trends Neurosciences. 2006;29:21–9.Google Scholar
- Hu H, Gan J, Jonas P. Interneurons. Fast-spiking, parvalbumin+ GABAergic interneurons: from cellular design to microcircuit function. Science. 2014;345:doi: 10.1126/science.1255263.
- Bensch S, Åkesson S, Irwin DE. The use of AFLP to find an informative SNP: genetic differences across a migratory divide in willow warblers. Mol Ecology. 2002;11:2359–66.View ArticleGoogle Scholar
- De Smet F, Moreau Y, Engelen K, Timmerman D, Vergote I, De Moor B. Balancing false positives and false negatives for the detection of differential expression in malignancies. Br J Cancer. 2004;91:1160–5.Google Scholar
- Robles JA, Qureshi SE, Stephen SJ, Wilson SR, Burden CJ, Taylor JM. Efficient experimental design and analysis strategies for the detection of differential expression using RNA-Sequencing. BMC Genomics. 2012;13:484.View ArticleGoogle Scholar
- Lein ES, Hawrylycz MJ, Ao N, Ayres M, Bensinger A, Bernard A, et al. Genome-wide atlas of gene expression in the adult mouse brain. Nature. 2007;445:168–76.View ArticleGoogle Scholar
- Cooper EC, Jan LY. Ion channel genes and human neurological disease: Recent progress, prospects, and challenges. Proc Natl Acad Sci USA. 1999;96:4759–66.View ArticleGoogle Scholar
- Peterson MP, Abolins-Abols M, Atwell JW, Rice RJ, Mila B, Ketterson ED. Variation in candidate genes CLOCK and ADCYAP1 does not consistently predict differences in migratory behavior in the songbird genus Junco. F1000Research. 2013;2:115.Google Scholar
- Ressler KJ, Mercer KB, Bradley B, et al. Post-traumatic stress disorder is associated with PACAP and the PAC1 receptor. Nature. 2011;470:492–7.View ArticleGoogle Scholar
- Åkesson S, Boström J, Liedvogel M, Muheim R. Animal Navigation. In: Akesson S, Hansson L-A, editors. Animal movements across scales. Oxford: Oxford University Press; 2014. p. 151-78.Google Scholar
- Muheim R, Phillips JB, Åkesson S. Polarized light cues underlie compass calibration in migratory songbirds. Science. 2006;313:837–9.View ArticleGoogle Scholar
- Muheim R, Åkesson S, Boström J, Liedvogel M. Sensory mechanisms of animal orientation and navigation. In: Hansson L-A, Akesson S, editors. Animal Movement Across Scales. Oxford University Press; 2014.Google Scholar
- Warren WC, Clayton DF, Ellegren H, Arnold AP, Hillier LW, Kunstner A, et al. The genome of a songbird. Nature. 2010;464:757–62.View ArticleGoogle Scholar
- Wolf JBW, Bayer T, Haubold B, Schilhabel M, Rosenstiel P, Tautz D. Nucleotide divergence vs. gene expression differentiation: comparative transcriptome sequencing in natural isolates from the carrion crow and its hybrid zone with the hooded crow. Mol Ecol. 2010;19:162–75.View ArticleGoogle Scholar
- Dzeja P, Terzic A. Adenylate Kinase and AMP Signaling Networks: Metabolic Monitoring, Signal Communication and Body Energy Sensing. Int J Mol Sciences. 2009;10:1729–72.View ArticleGoogle Scholar
- Simms BA, Zamponi GW. Neuronal Voltage-Gated Calcium Channels: Structure, Function, and Dysfunction. Neuron. 2014;82:24–45.View ArticleGoogle Scholar
- Bading H. Nuclear calcium signalling in the regulation of brain function. Nat Rev Neurosci. 2013;14:593–608.View ArticleGoogle Scholar
- Wiltschko W, Wiltschko R. Magnetic compass of European robins. Science. 1972;176:62–4.View ArticleGoogle Scholar
- Åkesson S, Bäckman J. Orientation in pied flycatchers: the relative importance of magnetic and visual information at dusk. Anim Behav. 1999;57:819–28.View ArticleGoogle Scholar