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.

Anthropogenic Change and the Process of Speciation.

Anthropogenic impacts on the environment alter speciation processes by affecting both geographical contexts and selection patterns on a worldwide scale. Here we review evidence of these effects. We find that human activities often generate spatial isolation between populations and thereby promote genetic divergence but also frequently cause sudden secondary contact and hybridization between diverging lineages. Human-caused environmental changes produce new ecological niches, altering selection in diverse ways that can drive diversification; but changes also often remove niches and cause extirpations. Human impacts that alter selection regimes are widespread and strong in magnitude, ranging from local changes in biotic and abiotic conditions to direct harvesting to global climate change. Altered selection, and evolutionary responses to it, impacts early-stage divergence of lineages, but does not necessarily lead toward speciation and persistence of separate species. Altogether, humans both promote and hinder speciation, although new species would form very slowly relative to anthropogenic hybridization, which can be nearly instantaneous. Speculating about the future of speciation, we highlight two key conclusions: (1) Humans will have a large influence on extinction and "despeciation" dynamics in the short term and on early-stage lineage divergence, and thus potentially speciation in the longer term, and (2) long-term monitoring combined with easily dated anthropogenic changes will improve our understanding of the processes of speciation. We can use this knowledge to preserve and restore ecosystems in ways that promote (re-)diversification, increasing future opportunities of speciation and enhancing biodiversity.

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.

Survey of haemosporidian parasite infections in an endangered high alpine bird.

The Brown-capped Rosy-Finch is an endangered high alpine specialist that breeds on cliffs in the Rocky Mountains of North America. We know little about the parasites that infect these birds but filling this knowledge gap is important for understanding their population decline. The aim of this study was to survey haemosporidian parasite infections in Brown-capped Rosy-Finches. We sampled 104 Brown-capped Rosy-Finches during their breeding season at six sampling sites spanning the Colorado Rocky Mountains where they are nearly endemic. We used nested PCR methods to screen birds for parasite infections, and Sanger sequencing data were used to identify parasite lineages. Four of the sampling sites had birds with infections. Females were more often infected than males (prevalence = 15.6% and 9.7%, respectively). We observed twice as many infected birds sampled in July compared to birds sampled in June (prevalence = 12.5% and 6.25%, respectively). The percent of infected birds by sampling site ranged from 0 to over 20%. In total we identified eight different genetic lineages of haemosporidian parasites infecting Brown-capped Rosy-Finches-seven were in the genus Leucocytozoon, and one was in the genus Haemoproteus. Network analysis clearly separates Haemoproteus from Leucocytozoon, with Leucocytozoon lineages comprising two major clusters. Based on reports made to the MalAvi database, all lineages in Cluster 1 have the same reported insect vector, Simulium silvestre. We report that Brown-capped Rosy-Finches experience infections with haemosporidian parasites and that the percentage of infected birds differed markedly between sampling sites. We hypothesize that vector ecology and associated variation in climate at sampling sites explain our observations.

The genetic basis of plumage coloration and elevation adaptation in a clade of recently diverged alpine and arctic songbirds.

Trait genetic architecture plays an important role in the probability that variation in that trait leads to divergence and speciation. In some cases, speciation may be driven by the generation of novel phenotypes through the recombination of genes associated with traits that are important for local adaptation or sexual selection. Here, we investigate the genetic basis of three plumage color traits, and one ecological trait, breeding elevation, in a recent avian radiation, the North American rosy-finches (Leucosticte spp.). We identify unique genomic regions associated with each trait and highlight 11 candidate genes. Among these are well-characterized melanogenesis genes, including Mitf and Tyrp1, and previously reported hypoxia-related genes including Egln1. Additionally, we use mitochondrial data to date the divergence of rosy-finch clades which appear to have diverged within the past 250 ky. Given the low levels of genome-wide differentiation among rosy-finch taxa, and evidence for extensive introgression in North America, plumage coloration and adaptation to high elevations have likely played large roles in generating the observed patterns of lineage divergence. The relative independence of these candidate regions across the genome suggests that recombination might have led to multiple phenotypes, and subsequent rosy-finch speciation, over short periods of time.

Hybridization between closely related songbirds is related to human habitat disturbance.

Human habitat disturbances can promote hybridization between closely related, but typically reproductively isolated, species. We explored whether human habitat disturbances are related to hybridization between two closely related songbirds, black-capped and mountain chickadees, using both genomic and citizen science data sets. First, we genotyped 409 individuals from across both species' ranges using reduced-representation genome sequencing and compared measures of genetic admixture to a composite measure of human landscape disturbance. Then, using eBird observations, we compared human landscape disturbance values for sites where phenotypically diagnosed hybrids were observed to locations where either parental species was observed to determine whether hybrid chickadees are reported in more disturbed areas. We found that hybridization between black-capped and mountain chickadees positively correlates with human habitat disturbances. From genomic data, we found that (1) hybrid index (HI) significantly increased with habitat disturbance, (2) more hybrids were sampled in disturbed habitats, (3) mean HIs were higher in disturbed habitats versus wild habitats, and (4) hybrids were detected in habitats with significantly higher disturbance values than parentals. Using eBird data, we found that both hybrid and black-capped chickadees were significantly more disturbance-associated than mountain chickadees. Surprisingly, we found that nearly every black-capped chickadee we sampled contained some proportion of hybrid ancestry, while we detected very few mountain chickadee backcrosses. Our results highlight that hybridization between black-capped and mountain chickadees is widespread, but initial hybridization is rare (few F1s were detected). We conclude that human habitat disturbances can erode pre-zygotic reproductive barriers between chickadees and that post-zygotic isolation is incomplete. Understanding what becomes of recently hybridizing species following large-scale habitat disturbances is a new, but pressing, consideration for successfully preserving genetic biodiversity in a rapidly changing world.

Sympatry leads to reduced body condition in chickadees that occasionally hybridize.

Both abiotic and biotic drivers influence species distributions. Abiotic drivers such as climate have received considerable attention, even though biotic drivers such as hybridization often interact with abiotic drivers. We sought to explore the (1) costs of co-occurrence for ecologically similar species that hybridize and (2) associations between ecological factors and condition to understand how abiotic and biotic factors influence species distributions. For two closely related and ecologically similar songbirds, black-capped and mountain chickadees, we characterized body condition, as a proxy for fitness, using a 1358-individual range-wide dataset. We compared body condition in sympatry and allopatry with several abiotic and biotic factors using species-specific generalized linear mixed models. We generated genomic data for a subset of 217 individuals to determine the extent of hybridization-driven admixture in our dataset. Within this data subset, we found that ~11% of the chickadees had hybrid ancestry, and all hybrid individuals had typical black-capped chickadee plumage. In the full dataset, we found that birds of both species, independent of demographic and abiotic factors, had significantly lower body condition when occurring in sympatry than birds in allopatry. This could be driven by either the inclusion of cryptic, likely poor condition, hybrids in our full dataset, competitive interactions in sympatry, or range edge effects. We are currently unable to discriminate between these mechanisms. Our findings have implications for mountain chickadees in particular, which will encounter more black-capped chickadees as black-capped chickadee ranges shift upslope and could lead to local declines in mountain chickadee populations.

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.

Estimating hybridization in the wild using citizen science data: A path forward.

Genomic evidence of introgression in natural populations has reinvigorated the study of hybridization in recent years. Still, it is largely unknown how frequently individual organisms mate across species lines. Recently, Justyn et al. suggested that eBird, one of the world's largest citizen science databases, may supply adequate data for estimating hybridization rates. Here, we compare Justyn et al.'s estimates-and their conclusions that hybridization is rare-with estimates from museum and molecular data. We also estimate hybridization using eBird observations from areas and times when hybridization is possible, namely, in contact zones during the breeding season. These estimates are all considerably higher than those reported in Justyn et al., emphasizing that inferences from multiple datasets can differ radically. Finally, we demonstrate an approach for predicting the location of hybrid zones using eBird data, which can be done with high confidence and with unprecedented resolution. We show that citizen science data, far from settling the question of how frequently bird species hybridize, instead offer a promising step toward more focused study of hybrid zones.

A supergene underlies linked variation in color and morphology in a Holarctic songbird.

The genetic architecture of a phenotype can have considerable effects on the evolution of a trait or species. Characterizing genetic architecture provides insight into the complexity of a given phenotype and, potentially, the role of the phenotype in evolutionary processes like speciation. We use genome sequences to investigate the genetic basis of phenotypic variation in redpoll finches (Acanthis spp.). We demonstrate that variation in redpoll phenotype is broadly controlled by a ~55-Mb chromosomal inversion. Within this inversion, we find multiple candidate genes related to melanogenesis, carotenoid coloration, and bill shape, suggesting the inversion acts as a supergene controlling multiple linked traits. A latitudinal gradient in ecotype distribution suggests supergene driven variation in color and bill morphology are likely under environmental selection, maintaining supergene haplotypes as a balanced polymorphism. Our results provide a mechanism for the maintenance of ecotype variation in redpolls despite a genome largely homogenized by gene flow.

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.

A highly invasive malaria parasite has expanded its range to non-migratory birds in North America.

Parasite range expansions are a direct consequence of globalization and are an increasing threat to biodiversity. Here, we report a recent range expansion of the SGS1 strain of a highly invasive parasite, Plasmodium relictum, to two non-migratory passerines in North America. Plasmodium relictum is considered one of the world's most invasive parasites and causes the disease avian malaria: this is the first reported case of SGS1 in wild non-migratory birds on the continent. Using a long-term database where researchers report avian malaria parasite infections, we summarized our current understanding of the geographical range of SGS1 and its known hosts. We also identified the most likely geographical region of this introduction event using the MSP1 allele. We hypothesize that this introduction resulted from movements of captive birds and subsequent spillover to native bird populations, via the presence of competent vectors and ecological fitting. Further work should be conducted to determine the extent to which SGS1 has spread following its introduction in North America.

Extensive historical and contemporary hybridization suggests premating isolation in Vermivora warblers is not strong: A reply to Confer et al.

We present comments on an article published by Confer et al. (Ecology and Evolution, 10, 2020). Confer et al. (2020) aggregate data from multiple studies of social pairing between Vermivora chrysoptera and V. cyanoptera, two wood warblers in the family Parulidae that hybridize extensively where they co-occur. From analysis of these data, they conclude there is near-complete reproductive isolation between these two species. In our reply, we show that this finding is not supported by other lines of evidence, and significant drawbacks of their study design preclude such strong conclusions. In our critique, we show that (a) coarse-scale plumage classifications cannot be used to accurately estimate hybrid ancestry in Vermivora; (b) extra-pair paternity is very high in Vermivora and is likely facilitating hybridization, yet was not considered by Confer et al. (2020), and we suggest this will have a substantial influence on the interpretation of reproductive isolation in the system; and (c) the central finding of strong total reproductive isolation is not compatible with the results of other long-term studies, which demonstrate low isolation and high gene flow. We conclude with a more comprehensive interpretation of hybridization and reproductive isolation in Vermivora warblers.

An isolated white-tailed deer (Odocoileus virginianus) population on St. John, US Virgin Islands shows low inbreeding and comparable heterozygosity to other larger populations.

This is the first study to document the genetic diversity of the white-tailed deer population on St. John, US Virgin Islands. The island population was founded by a small number of animals, has very limited hunting or predation, and recently experienced a reduction in size following an extended drought in 2015. DNA samples were collected from hair from 23 anesthetized adult deer (13 males, 10 females) ranging in age from 1 to 8 years (3.36 ± 1.9 years) and also from fecal DNA samples, for a total of 42 individuals analyzed for genetic diversity. The St. John deer data set averaged 4.19 alleles per marker and demonstrates the second lowest number of alleles (A) when compared to other populations of Odocoileus virginianus (4.19). Heterozygosity was similar to the other studies (0.54) with little evidence of inbreeding. To explain the level of heterozygosity and level of inbreeding within the St. John population, three hypotheses are proposed, including the effect of intrinsic biological traits within the population, a recent infusion of highly heterogeneous loci from North American populations, and a consistent level of immigration from a nearby island. Additional work is needed to further understand the genetic history of the St. John and regional deer populations.

Rapid speciation via the evolution of pre-mating isolation in the Iberá Seedeater.

Behavioral isolation can catalyze speciation and permit the slow accumulation of additional reproductive barriers between co-occurring organisms. We illustrate how this process occurs by examining the genomic and behavioral bases of pre-mating isolation between two bird species (Sporophila hypoxantha and the recently discovered S. iberaensis) that belong to the southern capuchino seedeaters, a recent, rapid radiation characterized by variation in male plumage coloration and song. Although these two species co-occur without obvious ecological barriers to reproduction, we document behaviors indicating species recognition by song and plumage traits and strong assortative mating associated with genomic regions underlying male plumage patterning. Plumage differentiation likely originated through the reassembly of standing genetic variation, indicating how novel sexual signals may quickly arise and maintain species boundaries.

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.

Phylogenomic Data Reveal Widespread Introgression Across the Range of an Alpine and Arctic Specialist.

Understanding how gene flow affects population divergence and speciation remains challenging. Differentiating one evolutionary process from another can be difficult because multiple processes can produce similar patterns, and more than one process can occur simultaneously. Although simple population models produce predictable results, how these processes balance in taxa with patchy distributions and complicated natural histories is less certain. These types of populations might be highly connected through migration (gene flow), but can experience stronger effects of genetic drift and inbreeding, or localized selection. Although different signals can be difficult to separate, the application of high-throughput sequence data can provide the resolution necessary to distinguish many of these processes. We present whole-genome sequence data for an avian species group with an alpine and arctic tundra distribution to examine the role that different population genetic processes have played in their evolutionary history. Rosy-finches inhabit high elevation mountaintop sky islands and high-latitude island and continental tundra. They exhibit extensive plumage variation coupled with low levels of genetic variation. Additionally, the number of species within the complex is debated, making them excellent for studying the forces involved in the process of diversification, as well as an important species group in which to investigate species boundaries. Total genomic variation suggests a broadly continuous pattern of allele frequency changes across the mainland taxa of this group in North America. However, phylogenomic analyses recover multiple distinct, well supported, groups that coincide with previously described morphological variation and current species-level taxonomy. Tests of introgression using D-statistics and approximate Bayesian computation reveal significant levels of introgression between multiple North American taxa. These results provide insight into the balance between divergent and homogenizing population genetic processes and highlight remaining challenges in interpreting conflict between different types of analytical approaches with whole-genome sequence data. [ABBA-BABA; approximate Bayesian computation; gene flow; phylogenomics; speciation; whole-genome sequencing.].

Leveraging eDNA to expand the study of hybrid zones.

Hybrid zones are important windows into ecological and evolutionary processes. Our understanding of the significance and prevalence of hybridization in nature has expanded with the generation and analysis of genome-spanning data sets. That said, most hybridization research still has restricted temporal and spatial resolution, which limits our ability to draw broad conclusions about evolutionary and conservation related outcomes. Here, we argue that rapidly advancing environmental DNA (eDNA) methodology could be adopted for studies of hybrid zones to increase temporal sampling (contemporary and historical), refine and geographically expand sampling density, and collect data for taxa that are difficult to directly sample. Genomic data in the environment offer the potential for near real-time biological tracking of hybrid zones, and eDNA provides broad, but as yet untapped, potential to address eco-evolutionary questions.

Genomic regions underlying metabolic and neuronal signaling pathways are temporally consistent in a moving avian hybrid zone.

The study of hybrid zones can provide insight into the genetic basis of species differences that are relevant for the maintenance of reproductive isolation. Hybrid zones can also provide insight into climate change, species distributions, and evolution. The hybrid zone between black-capped chickadees (Poecile atricapillus) and Carolina chickadees (Poecile carolinensis) is shifting northward in response to increasing winter temperatures but is not increasing in width. This pattern indicates strong selection against chickadees with admixed genomes. Using high-resolution genomic data, we identified regions of the genomes that are outliers in both time points and do not introgress between the species; these regions may be involved in the maintenance of reproductive isolation. Genes involved in metabolic regulation processes were overrepresented in this dataset. Several gene ontology categories were also temporally consistent-including glutamate signaling, synaptic transmission, and catabolic processes-but the nucleotide variants leading to this pattern were not. Our results support recent findings that hybrids between black-capped and Carolina chickadees have higher basal metabolic rates than either parental species and suffer spatial memory and problem-solving deficits. Metabolic breakdown, as well as spatial memory and problem-solving, in hybrid chickadees may act as strong postzygotic isolation mechanisms in this moving hybrid zone.