Genes and gene networks underlying spatial cognition in food-caching chickadees.

Substantial progress has been made in understanding the genetic architecture of phenotypes involved in a variety of evolutionary processes. Behavioral genetics remains, however, among the least understood. We explore the genetic architecture of spatial cognitive abilities in a wild passerine bird, the mountain chickadee (Poecile gambeli). Mountain chickadees cache thousands of seeds in the fall and require specialized spatial memory to recover these caches throughout the winter. We previously showed that variation in spatial cognition has a direct effect on fitness and has a genetic basis. It remains unknown which specific genes and developmental pathways are particularly important for shaping spatial cognition. To further dissect the genetic basis of spatial cognitive abilities, we combine experimental quantification of spatial cognition in wild chickadees with whole-genome sequencing of 162 individuals, a new chromosome-scale reference genome, and species-specific gene annotation. We have identified a set of genes and developmental pathways that play a key role in creating variation in spatial cognition and found that the mechanism shaping cognitive variation is consistent with selection against mildly deleterious non-coding mutations. Although some candidate genes were organized into connected gene networks, about half do not have shared regulation, highlighting that multiple independent developmental or physiological mechanisms contribute to variation in spatial cognitive abilities. A large proportion of the candidate genes we found are associated with synaptic plasticity, an intriguing result that leads to the hypothesis that certain genetic variants create antagonism between behavioral plasticity and long-term memory, each providing distinct benefits depending on ecological context.

Thermal acclimation in a non-migratory songbird occurs via changes to thermogenic capacity, but not conductance.

Thermoregulatory performance can be modified through changes in various subordinate traits, but the rate and magnitude of change in these traits is poorly understood. We investigated flexibility in traits that affect thermal balance between black-capped chickadees (Poecile atricapillus) acclimated for 6 weeks to cold (-5°C) or control (23°C) environments (n=7 per treatment). We made repeated measurements of basal and summit metabolic rates via flow-through respirometry and of body composition using quantitative magnetic resonance of live birds. At the end of the acclimation period, we measured thermal conductance of the combined feathers and skins. Cold-acclimated birds had a higher summit metabolic rate, reflecting a greater capacity for endogenous heat generation, and an increased lean mass. However, birds did not alter their thermal conductance. These results suggest that chickadees respond to cold stress by increasing their capacity for heat production rather than increasing heat retention, an energetically expensive strategy.

Quantitative Analyses of Coupling in Hybrid Zones.

In hybrid zones, whether barrier loci experience selection mostly independently or as a unit depends on the ratio of selection to recombination as captured by the coupling coefficient. Theory predicts a sharper transition between an uncoupled and coupled system when more loci affect hybrid fitness. However, the extent of coupling in hybrid zones has rarely been quantified. Here, we use simulations to characterize the relationship between the coupling coefficient and variance in clines across genetic loci. We then reanalyze 25 hybrid zone data sets and find that cline variances and estimated coupling coefficients form a smooth continuum from high variance and weak coupling to low variance and strong coupling. Our results are consistent with low rates of hybridization and a strong genome-wide barrier to gene flow when the coupling coefficient is much greater than 1, but also suggest that this boundary might be approached gradually and at a near constant rate over time.

The genetic basis of spatial cognitive variation in a food-caching bird.

Spatial cognition is used by most organisms to navigate their environment. Some species rely particularly heavily on specialized spatial cognition to survive, suggesting that a heritable component of cognition may be under natural selection. This idea remains largely untested outside of humans, perhaps because cognition in general is known to be strongly affected by learning and experience.1-4 We investigated the genetic basis of individual variation in spatial cognition used by non-migratory food-caching birds to recover food stores and survive harsh montane winters. Comparing the genomes of wild, free-living birds ranging from best to worst in their performance on a spatial cognitive task revealed significant associations with genes involved in neuron growth and development and hippocampal function. These results identify candidate genes associated with differences in spatial cognition and provide a critical link connecting individual variation in spatial cognition with natural selection. VIDEO ABSTRACT.

Variable Signatures of Selection Despite Conserved Recombination Landscapes Early in Speciation.

Recently diverged taxa often exhibit heterogeneous landscapes of genomic differentiation, characterized by regions of elevated differentiation on an otherwise homogeneous background. While divergence peaks are generally interpreted as regions responsible for reproductive isolation, they can also arise due to background selection, selective sweeps unrelated to speciation, and variation in recombination and mutation rates. To investigate the association between patterns of recombination and landscapes of genomic differentiation during the early stages of speciation, we generated fine-scale recombination maps for six southern capuchino seedeaters (Sporophila) and two subspecies of White Wagtail (Motacilla alba), two recent avian radiations in which divergent selection on pigmentation genes has likely generated peaks of differentiation. We compared these recombination maps to those of Collared (Ficedula albicollis) and Pied Flycatchers (Ficedula hypoleuca), non-sister taxa characterized by moderate genomic divergence and a heterogenous landscape of genomic differentiation shaped in part by background selection. Although recombination landscapes were conserved within all three systems, we documented a weaker negative correlation between recombination rate and genomic differentiation in the recent radiations. All divergence peaks between capuchinos, wagtails, and flycatchers were located in regions with lower-than-average recombination rates, and most divergence peaks in capuchinos and flycatchers fell in regions of exceptionally reduced recombination. Thus, co-adapted allelic combinations in these regions may have been protected early in divergence, facilitating rapid diversification. Despite largely conserved recombination landscapes, divergence peaks are specific to each focal comparison in capuchinos, suggesting that regions of elevated differentiation have not been generated by variation in recombination rate alone.

Asymmetric introgression reveals the genetic architecture of a plumage trait.

Genome-wide variation in introgression rates across hybrid zones offers a powerful opportunity for studying population differentiation. One poorly understood pattern of introgression is the geographic displacement of a trait implicated in lineage divergence from genome-wide population boundaries. While difficult to interpret, this pattern can facilitate the dissection of trait genetic architecture because traits become uncoupled from their ancestral genomic background. We studied an example of trait displacement generated by the introgression of head plumage coloration from personata to alba subspecies of the white wagtail. A previous study of their hybrid zone in Siberia revealed that the geographic transition in this sexual signal that mediates assortative mating was offset from other traits and genetic markers. Here we show that head plumage is associated with two small genetic regions. Despite having a simple genetic architecture, head plumage inheritance is consistent with partial dominance and epistasis, which could contribute to its asymmetric introgression.

Migratory divides coincide with reproductive barriers across replicated avian hybrid zones above the Tibetan Plateau.

Migratory divides are proposed to be catalysts for speciation across a diversity of taxa. However, it is difficult to test the relative contributions of migratory behaviour vs. other divergent traits to reproductive isolation. Comparing hybrid zones with and without migratory divides offers a rare opportunity to directly examine the contribution of divergent migratory behaviour to reproductive barriers. We show that across replicate sampling transects of two pairs of barn swallow (Hirundo rustica) subspecies, strong reproductive isolation coincided with a migratory divide spanning 20 degrees of latitude. A third subspecies pair exhibited no evidence for a migratory divide and hybridised extensively. Within migratory divides, overwintering habitats were associated with assortative mating, implicating a central contribution of divergent migratory behaviour to reproductive barriers. The remarkable geographic coincidence between migratory divides and genetic breaks supports a long-standing hypothesis that the Tibetan Plateau is a substantial barrier contributing to the diversity of Siberian avifauna.

Unifying Theoretical and Empirical Perspectives on Genomic Differentiation.

Differentiation is often heterogeneous across the genomes of diverging populations. Despite substantial recent progress, much work remains to improve our abilities to connect genomic patterns to underlying evolutionary processes. Crosstalk between theoretical and empirical research has shaped the field of evolutionary genetics since its foundation and needs to be greatly enhanced for modern datasets. We leverage recent insights from theoretical and empirical studies to identify existing gaps and suggest pathways across them. We stress the importance of reporting empirical data in standardized ways to enable meta-analyses and to facilitate parameterization of analyses and models. Additionally, a more comprehensive view of potential mechanisms - especially considering variable recombination rates and ubiquitous background selection - and their interactions should replace common, oversimplified assumptions.

Using Networks to Connect Individual-Level Reproductive Behavior to Population Patterns.

We propose an application of network analysis to determine which traits and behaviors predict fertilizations within and between populations. This approach quantifies how reproductive behavior between individuals shapes patterns of selection and gene flow, filling an important gap in our understanding of the connection between evolutionary processes and emergent patterns.

Extensive phenotypic diversification coexists with little genetic divergence and a lack of population structure in the White Wagtail subspecies complex (Motacilla alba).

Geographically clustered phenotypes often demonstrate consistent patterns in molecular markers, particularly mitochondrial DNA (mtDNA) traditionally used in phylogeographic studies. However, distinct evolutionary trajectories among traits and markers can lead to their discordance. First, geographic structure in phenotypic traits and nuclear molecular markers can be co-aligned but inconsistent with mtDNA (mito-nuclear discordance). Alternatively, phenotypic variation can have little to do with patterns in neither mtDNA nor nuclear markers. Disentangling between these distinct patterns can provide insight into the role of selection, demography and gene flow in population divergence. Here, we examined a previously reported case of strong inconsistency between geographic structure in mtDNA and plumage traits in a widespread polytypic bird species, the White Wagtail (Motacilla alba). We tested whether this pattern is due to mito-nuclear discordance or discrepancy between morphological evolution and both nuclear and mtDNA markers. We analysed population differentiation and structure across six out of nine commonly recognized subspecies using 17 microsatellite loci and a combination of microsatellites and plumage indices in a comprehensively sampled region of a contact between two subspecies. We did not find support for the mito-nuclear discordance hypothesis: nuclear markers indicated a subtle signal of genetic clustering only partially consistent with plumage groups, similar to previous findings that relied on mtDNA. We discuss evolutionary factors that could have shaped the intricate patterns of phenotypic diversification in the White wagtail and the role that repeated selection on plumage 'hotspots' and hybridization may have played.

Effects of assortative mate choice on the genomic and morphological structure of a hybrid zone between two bird subspecies.

Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between the alba and personata subspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437 SNP loci on a 3,000 km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome-wide divergence. Variation in only one trait-head plumage patterning-was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100 km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome-wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.

Hybridization Associated with Cycles of Ecological Succession in a Passerine Bird.

Identifying the diversity of contexts that can lead to hybridization is important for understanding its prevalence and dynamics in natural populations. Despite the potential of ecological succession to dramatically alter species co-occurrence and abundances, it is unknown whether it directly promotes hybridization and, if so, has long-lasting consequences. Here, we summarize 30 years of survey data across 10 populations to show that in western and mountain bluebirds, heterospecific pairing occurs during repeatable and transient colonization events at the early stages of species turnover. Despite mixed pairing occurring only during early succession, genetic data showed presence of hybrids at both early and late successional stages. Moreover, hybrids showed novel patterns of variation in morphology and behavior, emphasizing that even ephemeral contexts for hybridization can have important evolutionary consequences. Our results suggest that because ecological succession often brings together closely related competitors in disparate numbers but lasts for only a brief period of time, it may be a widespread but underappreciated context for hybridization.

Genomic variation across two barn swallow hybrid zones reveals traits associated with divergence in sympatry and allopatry.

Hybrid zones are geographic regions where isolating barriers between divergent populations are challenged by admixture. Identifying factors that facilitate or inhibit hybridization in sympatry can illuminate the processes that maintain those reproductive barriers. We analysed patterns of hybridization and phenotypic variation across two newly discovered hybrid zones between three subspecies of barn swallow (Hirundo rustica). These subspecies differ in ventral coloration and wing length, traits that are targets of sexual and natural selection, respectively, and are associated with genome-wide differentiation in allopatry. We tested the hypothesis that the degree of divergence in these traits is associated with the extent of hybridization in secondary contact. We applied measures of population structure based on >23,000 SNPs to confirm that named subspecies correspond to distinct genomic clusters, and assessed coincidence between geographic clines for ancestry and phenotype. Although gene flow was ongoing across both hybrid zones and pairwise FST between subspecies was extremely low, we found striking differences in the extent of hybridization. In the more phenotypically differentiated subspecies pair, clines for ancestry, wing length and ventral coloration were steep and coincident, suggestive of strong isolation and, potentially, selection associated with phenotype. In the less phenotypically differentiated pair, gene flow and phenotypic variation occurred over a wide geographic span, indicative of weaker isolation. Traits associated with genome-wide differentiation in allopatry may thus also contribute to isolation in sympatry. We discuss potentially important additional roles for evolutionary history and ecology in shaping variation in the extent hybridization between closely related pairs of subspecies.

Effects of asymmetric nuclear introgression, introgressive mitochondrial sweep, and purifying selection on phylogenetic reconstruction and divergence estimates in the Pacific clade of Locustella warblers.

When isolated but reproductively compatible populations expand geographically and meet, simulations predict asymmetric introgression of neutral loci from a local to invading taxon. Genetic introgression may affect phylogenetic reconstruction by obscuring topology and divergence estimates. We combined phylogenetic analysis of sequences from one mtDNA and 12 nuDNA loci with analysis of gene flow among 5 species of Pacific Locustella warblers to test for presence of genetic introgression and its effects on tree topology and divergence estimates. Our data showed that nuDNA introgression was substantial and asymmetrical among all members of superspecies groups whereas mtDNA showed no introgression except a single species pair where the invader's mtDNA was swept by mtDNA of the local species. This introgressive sweep of mtDNA had the opposite direction of the nuDNA introgression and resulted in the paraphyly of the local species' mtDNA haplotypes with respect to those of the invader. Тhe multilocus nuDNA species tree resolved all inter- and intraspecific relationships despite substantial introgression. However, the node ages on the species tree may be underestimated as suggested by the differences in node age estimates based on non-introgressing mtDNA and introgressing nuDNA. In turn, the introgressive sweep and strong purifying selection appear to elongate internal branches in the mtDNA gene tree.

Limited phylogeographic signal in sex-linked and autosomal loci despite geographically, ecologically, and phenotypically concordant structure of mtDNA variation in the Holarctic avian genus Eremophila.

Phylogeographic studies of Holarctic birds are challenging because they involve vast geographic scale, complex glacial history, extensive phenotypic variation, and heterogeneous taxonomic treatment across countries, all of which require large sample sizes. Knowledge about the quality of phylogeographic information provided by different loci is crucial for study design. We use sequences of one mtDNA gene, one sex-linked intron, and one autosomal intron to elucidate large scale phylogeographic patterns in the Holarctic lark genus Eremophila. The mtDNA ND2 gene identified six geographically, ecologically, and phenotypically concordant clades in the Palearctic that diverged in the Early-Middle Pleistocene and suggested paraphyly of the horned lark (E. alpestris) with respect to the Temminck's lark (E. bilopha). In the Nearctic, ND2 identified five subclades which diverged in the Late Pleistocene. They overlapped geographically and were not concordant phenotypically or ecologically. Nuclear alleles provided little information on geographic structuring of genetic variation in horned larks beyond supporting the monophyly of Eremophila and paraphyly of the horned lark. Multilocus species trees based on two nuclear or all three loci provided poor support for haplogroups identified by mtDNA. The node ages calculated using mtDNA were consistent with the available paleontological data, whereas individual nuclear loci and multilocus species trees appeared to underestimate node ages. We argue that mtDNA is capable of discovering independent evolutionary units within avian taxa and can provide a reasonable phylogeographic hypothesis when geographic scale, geologic history, and phenotypic variation in the study system are too complex for proposing reasonable a priori hypotheses required for multilocus methods. Finally, we suggest splitting the currently recognized horned lark into five Palearctic and one Nearctic species.

Geographic mode of speciation in a mountain specialist Avian family endemic to the Palearctic.

Mountains host greater avian diversity than lowlands at the same latitude due to their greater diversity of habitats stratified along an elevation gradient. Here we test whether this greater ecological heterogeneity promotes sympatric speciation. We selected accentors (Prunellidae), an avian family associated with mountains of the Palearctic, as a model system. Accentors differ in their habitat/elevation preferences and south-central Siberia and Himalayan regions each host 6 of the 13 species in the family. We used sequences of the mtDNA ND2 gene and the intron 9 of the Z chromosome specific ACO1 gene to reconstruct a complete species-level phylogeny of Prunellidae. The tree based on joint analysis of both loci was used to reconstruct the family's biogeographic history and to date the diversification events. We also analyzed the relationship between the node age and sympatry, to determine the geographic mode of speciation in Prunellidae. Our data suggest a Miocene origin of Prunellidae in the Himalayan region. The major division between alpine species (subgenus Laiscopus) and species associated with shrubs (subgenus Prunella) and initial diversification events within the latter happened within the Himalayan region in the Miocene and Pliocene. Accentors colonized other parts of the Palearctic during the Pliocene-Pleistocene transition. This spread across the Palearctic resulted in rapid diversification of accentors. With only a single exception dating to 0.91 Ma, lineages younger than 1.5 Ma are allopatric. In contrast, sympatry values for older nodes are >0. There was no relationship between node age and range symmetry. Allopatric speciation (not to include peripatric) is the predominant geographic mode of speciation in Prunellidae despite the favorable conditions for ecological diversification in the mountains and range overlaps among species.