Speciation in Ficedula flycatchers

Regulatory and evolutionary impact of DNA methylation in two songbird species and their naturally occurring F1 hybrids

BACKGROUND

Regulation of transcription by DNA methylation in 5'-CpG-3' context is a widespread mechanism allowing differential expression of genetically identical cells to persist throughout development. Consequently, differences in DNA methylation can reinforce variation in gene expression among cells, tissues, populations, and species. Despite a surge in studies on DNA methylation, we know little about the importance of DNA methylation in population differentiation and speciation. Here we investigate the regulatory and evolutionary impact of DNA methylation in five tissues of two Ficedula flycatcher species and their naturally occurring F1 hybrids.

RESULTS

We show that the density of CpG in the promoters of genes determines the strength of the association between DNA methylation and gene expression. The impact of DNA methylation on gene expression varies among tissues with the brain showing unique patterns. Differentially expressed genes between parental species are predicted by genetic and methylation differentiation in CpG-rich promoters. However, both these factors fail to predict hybrid misexpression suggesting that promoter mismethylation is not a main determinant of hybrid misexpression in Ficedula flycatchers. Using allele-specific methylation estimates in hybrids, we also determine the genome-wide contribution of cis- and trans effects in DNA methylation differentiation. These distinct mechanisms are roughly balanced in all tissues except the brain, where trans differences predominate.

CONCLUSIONS

Overall, this study provides insight on the regulatory and evolutionary impact of DNA methylation in songbirds.

Sympatry in a nightingale contact zone has no effect on host-specific blood parasite prevalence and lineage diversity

Parasites are a key driving force behind many ecological and evolutionary processes. Prevalence and diversity of parasites, as well as their effects on hosts, are not uniform across host species. As such, the potential parasite spillover between species can significantly influence outcomes of interspecific interactions. We screened two species of Luscinia nightingales for haemosporidian blood parasites (Plasmodium, Leucocytozoon and Haemoproteus) along an approximately 3000 km transect in Europe, incorporating areas of host distant allopatry, close allopatry and sympatry. We found significant differences in infection rates between the two host species, with common nightingales having much lower parasite prevalence than thrush nightingales (36.7% versus 83.8%). This disparity was mostly driven by Haemoproteus prevalence, which was significantly higher in thrush nightingales while common nightingales had a small, but significantly higher, Plasmodium prevalence. Furthermore, we found no effect of proximity to the contact zone on infection rate in either host species. Despite having lower infection prevalence, common nightingales were infected with a significantly higher diversity of parasite lineages than thrush nightingales, and lineage assemblages differed considerably between the two species, even in sympatry. This pattern was mostly driven by the large diversity of comparatively rare lineages, while the most abundant lineages were shared between the two host species. This suggests that, despite the close evolutionary relationships between the two nightingales, there are significant differences in parasite prevalence and diversity, regardless of the distance from the contact zone. This suggests that spillover of haemosporidian blood parasites is unlikely to contribute towards interspecific interactions in this system.

Evidence that genetic drift not adaptation drives fast-Z and large-Z effects in Ficedula flycatchers

The sex chromosomes have been hypothesized to play a key role in driving adaptation and speciation across many taxa. The reason for this is thought to be the hemizygosity of the heteromorphic part of sex chromosomes in the heterogametic sex, which exposes recessive mutations to natural and sexual selection. The exposure of recessive beneficial mutations increases their rate of fixation on the sex chromosomes, which results in a faster rate of evolution. In addition, genetic incompatibilities between sex-linked loci are exposed faster in the genomic background of hybrids of divergent lineages, which makes sex chromosomes contribute disproportionately to reproductive isolation. However, in birds, which show a Z/W sex determination system, the role of adaptation versus genetic drift as the driving force of the faster differentiation of the Z chromosome (fast-Z effect) and the disproportionate role of the Z chromosome in reproductive isolation (large-Z effect) are still debated. Here, we address this debate in the bird genus Ficedula flycatchers based on population-level whole-genome sequencing data of six species. Our analysis provides evidence for both faster lineage sorting and reduced gene flow on the Z chromosome than the autosomes. However, these patterns appear to be driven primarily by the increased role of genetic drift on the Z chromosome, rather than an increased rate of adaptive evolution. Genomic scans of selective sweeps and fixed differences in fact suggest a reduced action of positive selection on the Z chromosome.

Species-specific song responses emerge as a by-product of tuning to the local dialect

Oscine birds preferentially respond to certain sounds over others from an early age, which focuses subsequent learning onto sexually relevant songs.^1,2,3 Songs vary both across species and, due to cultural evolution, among populations of the same species. As a result, early song responses are expected to be shaped by selection both to avoid the fitness costs of cross-species learning⁴ and to promote learning of population-typical songs.⁵ These sources of selection are not mutually exclusive but can result in distinct geographic patterns of song responses in juvenile birds: if the risks of interspecific mating are the main driver of early song discrimination, then discrimination should be strongest where closely related species co-occur.⁴ In contrast, if early discrimination primarily facilitates learning local songs, then it should be tuned to songs typical of the local dialect.^5,6,7 Here, we experimentally assess the drivers of song discrimination in nestling pied flycatchers (Ficedula hypoleuca). We first demonstrate that early discrimination against the songs of the closely related collared flycatcher (F. albicollis) is not strongly affected by co-occurrence. Second, across six European populations, we show that nestlings' early song responses are tuned to their local song dialect and that responses to the songs of collared flycatchers are similarly weak as to those of other conspecific dialects. Taken together, these findings provide clear experimental support for the hypothesis that cultural evolution, in conjunction with associated learning predispositions, drives the emergence of pre-mating reproductive barriers.

Individual differences in song plasticity in response to social stimuli and singing position

Individual animals can react to the changes in their environment by exhibiting behaviors in an individual‐specific way leading to individual differences in phenotypic plasticity. However, the effect of multiple environmental factors on multiple traits is rarely tested. Such a complex approach is necessary to assess the generality of plasticity and to understand how among‐individual differences in the ability to adapt to changing environments evolve. This study examined whether individuals adjust different song traits to varying environmental conditions in the collared flycatcher (Ficedula albicollis), a passerine with complex song. We also aimed to reveal among‐individual differences in behavioral responses by testing whether individual differences in plasticity were repeatable. The presence of general plasticity across traits and/or contexts was also tested. To assess plasticity, we documented (1) short‐scale temporal changes in song traits in different social contexts (after exposition to male stimulus, female stimulus or without stimuli), and (2) changes concerning the height from where the bird sang (singing position), used as a proxy of predation risk and acoustic transmission conditions. We found population‐level relationships between singing position and both song length (SL) and complexity, as well as social context‐dependent temporal changes in SL and maximum frequency (MF). We found among‐individual differences in plasticity of SL and MF along both the temporal and positional gradients. These among‐individual differences in plasticity were repeatable. Some of the plastic responses correlated across different song traits and environmental gradients. Overall, our results show that the plasticity of bird song (1) depends on the social context, (2) exists along different environmental gradients, and (3) there is evidence for trade‐offs between the responses of different traits to different environmental variables. Our results highlight the need to consider individual differences and to investigate multiple traits along multiple environmental axes when studying behavioral plasticity.

Divergent selection on behavioural and chemical traits between reproductively isolated populations of Drosophila melanogaster

Speciation is driven by traits that can act to prevent mating between nascent lineages, including male courtship and female preference for male traits. Mating barriers involving these traits evolve quickly because there is strong selection on males and females to maximize reproductive success, and the tight co-evolution of mating interactions can lead to rapid diversification of sexual behaviour. Populations of Drosophila melanogaster show strong asymmetrical reproductive isolation that is correlated with geographic origin. Using strains that capture natural variation in mating traits, we ask two key questions: which specific male traits are females selecting, and are these traits under divergent sexual selection? These questions have proven extremely challenging to answer, because even in closely related lineages males often differ in multiple traits related to mating behaviour. We address these questions by estimating selection gradients for male courtship and cuticular hydrocarbons for two different female genotypes. We identify specific behaviours and particular cuticular hydrocarbons that are under divergent sexual selection and could potentially contribute to premating reproductive isolation. Additionally, we report that a subset of these traits are plastic; males adjust these traits based on the identity of the female genotype they interact with. These results suggest that even when male courtship is not fixed between lineages, ongoing selection can act on traits that are important for reproductive isolation.

Asymmetric allelic introgression across a hybrid zone of the coal tit (Periparus ater) in the central Himalayas

In the Himalayas, a number of secondary contact zones have been described for vicariant vertebrate taxa. However, analyses of genetic divergence and admixture are missing for most of these examples. In this study, we provide a population genetic analysis for the coal tit (Periparus ater) hybrid zone in Nepal. Intermediate phenotypes between the distinctive western "spot-winged tit" (P. a. melanolophus) and Eastern Himalayan coal tits (P. a. aemodius) occur across a narrow range of <100 km in western Nepal. As a peculiarity, another distinctive cinnamon-bellied form is known from a single population so far. Genetic admixture of western and eastern mitochondrial lineages was restricted to the narrow zone of phenotypically intermediate populations. The cline width was estimated 46 km only with a center close to the population of the cinnamon-bellied phenotype. In contrast, allelic introgression of microsatellite loci was asymmetrical from eastern P. a. aemodius into far western populations of phenotypic P. a. melanolophus but not vice versa. Accordingly, the microsatellite cline was about 3.7 times wider than the mitochondrial one.

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.

Should females prefer old males?

Whether females should prefer to mate with old males is controversial. Old males may sire offspring of low quality because of an aging germline, but their proven ability to reach an old age can also be an excellent indicator of superior genetic quality, especially in natural populations. These genetic effects are, however, hard to study in nature, because they are often confounded with direct benefits offered by old males to the female, such as experience and high territory quality. We, therefore, used naturally occurring extra‐pair young to disentangle different aspects of male age on female fitness in a natural population of collared flycatchers because any difference between within‐ and extra‐pair young within a nest should be caused by paternal genetic effects only. Based on 18 years of long‐term data, we found that females paired with older males as social partners experienced an overall reproductive advantage. However, offspring sired by old males were of lower quality as compared to their extra‐pair half‐siblings, whereas the opposite was found in nests attended by young males. These results imply a negative genetic effect of old paternal age, given that extra‐pair males are competitive middle‐age males. Thus, offspring may benefit from being sired by young males but raised by old males, to maximize both genetic and direct effects. Our results show that direct and genetic benefits from pairing with old males may act in opposing directions and that the quality of the germline may deteriorate before other signs of senescence become obvious.

Low Heritability but Significant Early Environmental Effects on Resting Metabolic Rate in a Wild Passerine

AbstractPredicting the impact of climate change on biodiversity requires understanding the adaptation potential of wild organisms. Evolutionary responses depend on the additive genetic variation associated with the phenotypic traits targeted by selection. We combine 5 years of cross-fostering experiments, measurements of resting metabolic rate (RMR) on nearly 200 wild collared flycatcher (Ficedula albicollis) nestlings, and animal models using a 17-year pedigree to evaluate the potential for an evolutionary response to changing environmental conditions. Contrary to other avian studies, we find no significant heritability of whole-organism, mass-independent, or mass-specific RMR, but we report a strong effect of nest environment instead. We therefore conclude that variation in nestling RMR is explained by variation in the early-life environment provided by the parents. We discuss possible underlying specific parental effects and the importance of taking different mechanisms into account to understand how animals phenotypically adapt (or fail to adapt) to climate change.

The effect of social environment on bird song: listener-specific expression of a sexual signal

Animal signals should consistently differ among individuals to convey distinguishable information about the signalers. However, behavioral display signals, such as bird song are also loaded with considerable within-individual variance with mostly unknown function. We hypothesized that the immediate social environment may play a role in mediating such variance component, and investigated in the collared flycatcher (Ficedula albicollis) if the identity and quality of listeners could affect song production in signalers. After presenting territorial males with either a female or male social stimulus, we found in the subsequent song recordings that the among-stimulus effects corresponded to non-zero variance components in several acoustic traits indicating that singing males are able to plastically adjust their songs according to stimulus identity. Male and female stimuli elicited different responses as the identity of the female stimuli affected song complexity only, while the identity of male stimuli altered also song length, maximum frequency, and song rate. The stimulus-specific effect on song in some cases decreased with time, being particularly detectable right after the removal of the stimulus and ceasing later, but this pattern varied across the sex of the stimulus and the song traits. We were able to identify factors that can explain the among-stimulus effects (e.g., size and quality of the stimuli) with roles that also varied among song traits. Our results confirm that the variable social environment can raise considerable variation in song performance, highlighting that within-individual plasticity of bird song can play important roles in sexual signaling.

Positive selection plays a major role in shaping signatures of differentiation across the genomic landscape of two independent Ficedula flycatcher species pairs

A current debate within population genomics surrounds the relevance of patterns of genomic differentiation between closely related species for our understanding of adaptation and speciation. Mounting evidence across many taxa suggests that the same genomic regions repeatedly develop elevated differentiation in independent species pairs. These regions often coincide with high gene density and/or low recombination, leading to the hypothesis that the genomic differentiation landscape mostly reflects a history of background selection, and reveals little about adaptation or speciation. A comparative genomics approach with multiple independent species pairs at a timescale where gene flow and ILS are negligible permits investigating whether different evolutionary processes are responsible for generating lineage-specific versus shared patterns of species differentiation. We use whole-genome resequencing data of 195 individuals from four Ficedula flycatcher species comprising two independent species pairs: collared and pied flycatchers, and red-breasted and taiga flycatchers. We found that both shared and lineage-specific F_ST peaks could partially be explained by selective sweeps, with recurrent selection likely to underlie shared signatures of selection, whereas indirect evidence supports a role of recombination landscape evolution in driving lineage-specific signatures of selection. This work therefore provides evidence for an interplay of positive selection and recombination to genomic landscape evolution.

Adaptive coloration in pied flycatchers (Ficedula hypoleuca)-The devil is in the detail

Understanding the origin and persistence of phenotypic variation within and among populations is a major goal in evolutionary biology. However, the eagerness to find unadulterated explanatory models in combination with difficulties in publishing replicated studies may lead to severe underestimations of the complexity of selection patterns acting in nature. One striking example is variation in plumage coloration in birds, where the default adaptive explanation often is that brightly colored individuals signal superior quality across environmental conditions and therefore always should be favored by directional mate choice. Here, we review studies on the proximate determination and adaptive function of coloration traits in male pied flycatchers (Ficedula hypoleuca). From numerous studies, we can conclude that the dark male color phenotype is adapted to a typical northern climate and functions as a dominance signal in male-male competition over nesting sites, and that the browner phenotypes are favored by relaxed intraspecific competition with more dominant male collared flycatchers (Ficedula albicollis) in areas where the two species co-occur. However, the role of avoidance of hybridization in driving character displacement in plumage between these two species may not be as important as initially thought. The direction of female choice on male coloration in pied flycatchers is not simply as opposite in direction in sympatry and allopatry as traditionally expected, but varies also in relation to additional contexts such as climate variation. While some of the heterogeneity in the observed relationships between coloration and fitness probably indicate type 1 errors, we strongly argue that environmental heterogeneity and context-dependent selection play important roles in explaining plumage color variation in this species, which probably also is the case in many other species studied in less detail.

Tissue-specific patterns of regulatory changes underlying gene expression differences among Ficedula flycatchers and their naturally occurring F1 hybrids

Changes in interacting cis- and trans-regulatory elements are important candidates for Dobzhansky-Muller hybrid incompatibilities and may contribute to hybrid dysfunction by giving rise to misexpression in hybrids. To gain insight into the molecular mechanisms and determinants of gene expression evolution in natural populations, we analyzed the transcriptome from multiple tissues of two recently diverged Ficedula flycatcher species and their naturally occurring F₁ hybrids. Differential gene expression analysis revealed that the extent of differentiation between species and the set of differentially expressed genes varied across tissues. Common to all tissues, a higher proportion of Z-linked genes than autosomal genes showed differential expression, providing evidence for a fast-Z effect. We further found clear signatures of hybrid misexpression in brain, heart, kidney, and liver. However, while testis showed the highest divergence of gene expression among tissues, it showed no clear signature of misexpression in F₁ hybrids, even though these hybrids were found to be sterile. It is therefore unlikely that incompatibilities between cis-trans regulatory changes explain the observed sterility. Instead, we found evidence that cis-regulatory changes play a significant role in the evolution of gene expression in testis, which illustrates the tissue-specific nature of cis-regulatory evolution bypassing constraints associated with pleiotropic effects of genes.

Heterospecific mating interactions as an interface between ecology and evolution

Reproductive interference (costly interspecific sexual interactions) is well-understood to promote divergence in mating-relevant traits (i.e. reproductive character displacement: RCD), but it can also reduce population growth, eventually leading to local extinction of one of the species. The ecological and evolutionary processes driven by reproductive interference can interact with each other. These interactions are likely to influence whether the outcome is coexistence or extinction, but remain little studied. In this paper, we first develop an eco-evolutionary perspective on reproductive interference by integrating ecological and evolutionary processes in a common framework. We also present a simple model to demonstrate the eco-evolutionary dynamics of reproductive interference. We then identify a number of factors that are likely to influence the relative likelihoods of extinction or RCD. We discuss particularly relevant factors by classifying them into four categories: the nature of the traits responding to selection, the mechanisms determining the expression of these traits, mechanisms of reproductive interference and the ecological background. We highlight previously underappreciated ways in which these factors may influence the relative likelihoods of RCD and local extinction. By doing so, we also identify questions and future directions that will increase our holistic understanding of the outcomes of reproductive interference.

Symmetrical discrimination despite weak song differentiation in 2 suboscine bird sister species

Song mediates territorial competition and mate choice in birds and population divergence in this signal can have important evolutionary consequences. For example, divergent songs can act in specific recognition and limit gene flow and, hence, have a fundamental role on the origin and/or integrity of evolutionary lineages. Especially interesting systems to test the role of song in specific recognition are species pairs that present small structural differences in this signal. Here, we perform song play-back experiments on males of a long-diverged sister pair of Neotropical Suboscine species, the squamate antbird (Myrmoderus squamosus) and the white-bibbed antbird (Myrmoderus loricatus), which occur in parapatry in the Atlantic Forest and that overlap extensively in song variation. Previous evidence indicates that genetic introgression between these species is either absent or negligible, suggesting that vocal discrimination or other mechanisms function as effective barriers to gene flow. Our results show that responses to heterospecific songs were symmetrical and intermediary compared with responses to conspecific songs in both species. A stronger response to conspecific territorial songs suggests that conspecific individuals pose greater competitive threat than heterospecifics. An important implication of our study is that even small song differences can play an important role in specific recognition.

Differential introgression of a female competitive trait in a hybrid zone between sex-role reversed species

Mating behavior between recently diverged species in secondary contact can impede or promote reproductive isolation. Traditionally, researchers focus on the importance of female mate choice and male-male competition in maintaining or eroding species barriers. Although female-female competition is widespread, little is known about its role in the speciation process. Here, we investigate a case of interspecific female competition and its influence on patterns of phenotypic and genetic introgression between species. We examine a hybrid zone between sex-role reversed, Neotropical shorebird species, the northern jacana (Jacana spinosa) and wattled jacana (J. jacana), in which female-female competition is a major determinant of reproductive success. Previous work found that females of the more aggressive and larger species, J. spinosa, disproportionately mother hybrid offspring, potentially by monopolizing breeding territories in sympatry with J. jacana. We find a cline shift of female body mass relative to the genetic center of the hybrid zone, consistent with asymmetric introgression of this competitive trait. We suggest that divergence in sexual characteristics between sex-role reversed females can influence patterns of gene flow upon secondary contact, similar to males in systems with more typical sex roles.

Interspecific transfer of parasites following a range-shift in Ficedula flycatchers

Human-induced climate change is expected to cause major biotic changes in species distributions and thereby including escalation of novel host-parasite associations. Closely related host species that come into secondary contact are especially likely to exchange parasites and pathogens. Both the Enemy Release Hypothesis (where invading hosts escape their original parasites) and the Novel Weapon Hypothesis (where invading hosts bring new parasites that have detrimental effects on native hosts) predict that the local host will be most likely to experience a disadvantage. However, few studies evaluate the occurrence of interspecific parasite transfer by performing wide-scale geographic sampling of pathogen lineages, both within and far from host contact zones. In this study, we investigate how haemosporidian (avian malaria) prevalence and lineage diversity vary in two, closely related species of passerine birds; the pied flycatcher Ficedula hypoleuca and the collared flycatcher F. albicollis in both allopatry and sympatry. We find that host species is generally a better predictor of parasite diversity than location, but both prevalence and diversity of parasites vary widely among populations of the same bird species. We also find a limited and unidirectional transfer of parasites from pied flycatchers to collared flycatchers in a recent contact zone. This study therefore rejects both the Enemy Release Hypothesis and the Novel Weapon Hypothesis and highlights the complexity and importance of studying host-parasite relationships in an era of global climate change and species range shifts.

Maladaptive learning and memory in hybrids as a reproductive isolating barrier

Selection against hybrid offspring, or postzygotic reproductive isolation, maintains species boundaries in the face of gene flow from hybridization. In this review, we propose that maladaptive learning and memory in hybrids is an important, but overlooked form of postzygotic reproductive isolation. Although a role for learning in premating isolation has been supported, whether learning deficiencies can contribute to postzygotic isolation has rarely been tested. We argue that the novel genetic combinations created by hybridization have the potential to impact learning and memory abilities through multiple possible mechanisms, and that any displacement from optima in these traits is likely to have fitness consequences. We review evidence supporting the potential for hybridization to affect learning and memory, and evidence of links between learning abilities and fitness. Finally, we suggest several avenues for future research. Given the importance of learning for fitness, especially in novel and unpredictable environments, maladaptive learning and memory in hybrids may be an increasingly important source of postzygotic reproductive isolation.

Difference in plasticity of resting metabolic rate - the proximate explanation to different niche breadth in sympatric Ficedula flycatchers

Variation in relative fitness of competing recently formed species across heterogeneous environments promotes coexistence. However, the physiological traits mediating such variation in relative fitness have rarely been identified. Resting metabolic rate (RMR) is tightly associated with life history strategies, thermoregulation, diet use, and inhabited latitude and could therefore moderate differences in fitness responses to fluctuations in local environments, particularly when species have adapted to different climates in allopatry. We work in a long‐term study of collared (Ficedula albicollis) and pied flycatchers (Ficedula hypoleuca) in a recent hybrid zone located on the Swedish island of Öland in the Baltic Sea. Here, we explore whether differences in RMR match changes in relative performance of growing flycatcher nestlings across environmental conditions using an experimental approach. The fitness of pied flycatchers has previously been shown to be less sensitive to the mismatch between the peak in food abundance and nestling growth among late breeders. Here, we find that pied flycatcher nestlings have lower RMR in response to higher ambient temperatures (associated with low food availability). We also find that experimentally relaxed nestling competition is associated with an increased RMR in this species. In contrast, collared flycatcher nestlings did not vary their RMR in response to these environmental factors. Our results suggest that a more flexible nestling RMR in pied flycatchers is responsible for the better adaptation of pied flycatchers to the typical seasonal changes in food availability experienced in this hybrid zone. Generally, subtle physiological differences that have evolved when species were in allopatry may play an important role to patterns of competition, coexistence, or displacements between closely related species in secondary contact.

Heterogeneous Patterns of Genetic Diversity and Differentiation in European and Siberian Chiffchaff (Phylloscopus collybita abietinus/P. tristis)

Identification of candidate genes for trait variation in diverging lineages and characterization of mechanistic underpinnings of genome differentiation are key steps toward understanding the processes underlying the formation of new species. Hybrid zones provide a valuable resource for such investigations, since they allow us to study how genomes evolve as species exchange genetic material and to associate particular genetic regions with phenotypic traits of interest. Here, we use whole-genome resequencing of both allopatric and hybridizing populations of the European (Phylloscopus collybita abietinus) and the Siberian chiffchaff (P. tristis)—two recently diverged species which differ in morphology, plumage, song, habitat, and migration—to quantify the regional variation in genome-wide genetic diversity and differentiation, and to identify candidate regions for trait variation. We find that the levels of diversity, differentiation, and divergence are highly heterogeneous, with significantly reduced global differentiation, and more pronounced differentiation peaks in sympatry than in allopatry. This pattern is consistent with regional differences in effective population size and recurrent background selection or selective sweeps reducing the genetic diversity in specific regions prior to lineage divergence, but the data also suggest that postdivergence selection has resulted in increased differentiation and fixed differences in specific regions. We find that hybridization and backcrossing is common in sympatry, and that phenotype is a poor predictor of the genomic composition of sympatric birds. The combination of a differentiation scan approach with identification of fixed differences pinpoint a handful of candidate regions that might be important for trait variation between the two species.

Song discrimination by nestling collared flycatchers during early development

Pre-zygotic isolation is often maintained by species-specific signals and preferences. However, in species where signals are learnt, as in songbirds, learning errors can lead to costly hybridization. Song discrimination expressed during early developmental stages may ensure selective learning later in life but can be difficult to demonstrate before behavioural responses are obvious. Here, we use a novel method, measuring changes in metabolic rate, to detect song perception and discrimination in collared flycatcher embryos and nestlings. We found that nestlings as early as 7 days old respond to song with increased metabolic rate, and, by 9 days old, have increased metabolic rate when listening to conspecific when compared with heterospecific song. This early discrimination between songs probably leads to fewer heterospecific matings, and thus higher fitness of collared flycatchers living in sympatry with closely related species.

Explosive radiation and spatial expansion across the cold environments of the Old World in an avian family

Our objective was to elucidate the biogeography and speciation patterns in an entire avian family, which shows a complex pattern of overlapping and nonoverlapping geographical distributions, and much variation in plumage, but less in size and structure. We estimated the phylogeny and divergence times for all of the world's species of Prunella based on multiple genetic loci, and analyzed morphometric divergence and biogeographical history. The common ancestor of Prunella was present in the Sino‐Himalayan Mountains or these mountains and Central Asia–Mongolia more than 9 million years ago (mya), but a burst of speciations took place during the mid‐Pliocene to early Pleistocene. The relationships among the six primary lineages resulting from that differentiation are unresolved, probably because of the rapid radiation. A general increase in sympatry with increasing time since divergence is evident. With one exception, species in clades younger than c. 3.7 my are allopatric. Species that are widely sympatric, including the most recently diverged (2.4 mya) sympatric sisters, are generally more divergent in size/structure than allo‐/parapatric close relatives. The distributional pattern and inferred ages suggest divergence in allopatry and substantial waiting time until secondary contact, likely due to competitive exclusion. All sympatrically breeding species are ecologically segregated, as suggested by differences in size/structure and habitat. Colonizations of new areas were facilitated during glacial periods, followed by fragmentation during interglacials—contrary to the usual view that glacial periods resulted mainly in fragmentations.

Insights into the genetic architecture of morphological traits in two passerine bird species

Knowledge about the underlying genetic architecture of phenotypic traits is needed to understand and predict evolutionary dynamics. The number of causal loci, magnitude of the effects and location in the genome are, however, still largely unknown. Here, we use genome-wide single-nucleotide polymorphism (SNP) data from two large-scale data sets on house sparrows and collared flycatchers to examine the genetic architecture of different morphological traits (tarsus length, wing length, body mass, bill depth, bill length, total and visible badge size and white wing patches). Genomic heritabilities were estimated using relatedness calculated from SNPs. The proportion of variance captured by the SNPs (SNP-based heritability) was lower in house sparrows compared with collared flycatchers, as expected given marker density (6348 SNPs in house sparrows versus 38 689 SNPs in collared flycatchers). Indeed, after downsampling to similar SNP density and sample size, this estimate was no longer markedly different between species. Chromosome-partitioning analyses demonstrated that the proportion of variance explained by each chromosome was significantly positively related to the chromosome size for some traits and, generally, that larger chromosomes tended to explain proportionally more variation than smaller chromosomes. Finally, we found two genome-wide significant associations with very small-effect sizes. One SNP on chromosome 20 was associated with bill length in house sparrows and explained 1.2% of phenotypic variation (V_P), and one SNP on chromosome 4 was associated with tarsus length in collared flycatchers (3% of V_P). Although we cannot exclude the possibility of undetected large-effect loci, our results indicate a polygenic basis for morphological traits.

Population genetic isolation and limited connectivity in the purple finch (Haemorhous purpureus)

Using a combination of mitochondrial and z-linked sequences, microsatellite data, and spatio-geographic modeling, we examined historical and contemporary factors influencing the population genetic structure of the purple finch (Haemorhous purpureus). Mitochondrial DNA data show the presence of two distinct groups corresponding to the two subspecies, H. p. purpureus and H. p. californicus. The two subspecies likely survived in separate refugia during the last glacial maximum, one on the Pacific Coast and one east of the Rocky Mountains, and now remain distinct lineages with little evidence of gene flow between them. Southwestern British Columbia is a notable exception, as subspecies mixing between central British Columbia and Vancouver Island populations suggests a possible contact zone in this region. Z-linked data support two mitochondrial groups; however, Coastal Oregon and central British Columbia sites show evidence of mixing. Contemporary population structure based on microsatellite data identified at least six genetic clusters: three H. p. purpureus clusters, two H. p. californicus clusters, and one mixed cluster, which likely resulted from high site fidelity and isolation by distance, combined with sexual selection on morphological characters reinforcing subspecies differences.

Hybrid Dysfunction Expressed as Elevated Metabolic Rate in Male Ficedula Flycatchers

Studies of ecological speciation are often biased towards extrinsic sources of selection against hybrids, resulting from intermediate hybrid morphology, but the knowledge of how genetic incompatibilities accumulate over time under natural conditions is limited. Here we focus on a physiological trait, metabolic rate, which is central to life history strategies and thermoregulation but is also likely to be sensitive to mismatched mitonuclear interactions. We measured the resting metabolic rate of male collared, and pied flycatchers as well as of naturally occurring F1 hybrid males, in a recent hybrid zone. We found that hybrid males had a higher rather than intermediate metabolic rate, which is indicative of hybrid physiological dysfunction. Fitness costs associated with elevated metabolic rate are typically environmentally dependent and exaggerated under harsh conditions. By focusing on male hybrid dysfunction in an eco-physiological trait, our results contribute to the general understanding of how combined extrinsic and intrinsic sources of hybrid dysfunction build up under natural conditions.

Competition-driven build-up of habitat isolation and selection favoring modified dispersal patterns in a young avian hybrid zone

Competition-driven evolution of habitat isolation is an important mechanism of ecological speciation but empirical support for this process is often indirect. We examined how an on-going displacement of pied flycatchers from their preferred breeding habitat by collared flycatchers in a young secondary contact zone is associated with (a) access to an important food resource (caterpillar larvae), (b) immigration of pied flycatchers in relation to habitat quality, and (c) the risk of hybridization in relation to habitat quality. Over the past 12 years, the estimated access to caterpillar larvae biomass in the habitat surrounding the nests of pied flycatchers has decreased by a fifth due to shifted establishment possibilities, especially for immigrants. However, breeding in the high quality habitat has become associated with such a high risk of hybridization for pied flycatchers that overall selection currently favors pied flycatchers that were forced to immigrate into the poorer habitats (despite lower access to preferred food items). Our results show that competition-driven habitat segregation can lead to fast habitat isolation, which per se caused an opportunity for selection to act in favor of future "voluntarily" altered immigration patterns and possibly strengthened habitat isolation through reinforcement.

Genome-wide association mapping in a wild avian population identifies a link between genetic and phenotypic variation in a life-history trait

Understanding the genetic basis of traits involved in adaptation is a major challenge in evolutionary biology but remains poorly understood. Here, we use genome-wide association mapping using a custom 50 k single nucleotide polymorphism (SNP) array in a natural population of collared flycatchers to examine the genetic basis of clutch size, an important life-history trait in many animal species. We found evidence for an association on chromosome 18 where one SNP significant at the genome-wide level explained 3.9% of the phenotypic variance. We also detected two suggestive quantitative trait loci (QTLs) on chromosomes 9 and 26. Fitness differences among genotypes were generally weak and not significant, although there was some indication of a sex-by-genotype interaction for lifetime reproductive success at the suggestive QTL on chromosome 26. This implies that sexual antagonism may play a role in maintaining genetic variation at this QTL. Our findings provide candidate regions for a classic avian life-history trait that will be useful for future studies examining the molecular and cellular function of, as well as evolutionary mechanisms operating at, these loci.

PSMC analysis of effective population sizes in molecular ecology and its application to black-and-white Ficedula flycatchers

Climatic fluctuations during the Quaternary period governed the demography of species and contributed to population differentiation and ultimately speciation. Studies of these past processes have previously been hindered by a lack of means and genetic data to model changes in effective population size (Ne ) through time. However, based on diploid genome sequences of high quality, the recently developed pairwise sequentially Markovian coalescent (PSMC) can estimate trajectories of changes in Ne over considerable time periods. We applied this approach to resequencing data from nearly 200 genomes of four species and several populations of the Ficedula species complex of black-and-white flycatchers. Ne curves of Atlas, collared, pied and semicollared flycatcher converged 1-2 million years ago (Ma) at an Ne of ≈ 200 000, likely reflecting the time when all four species last shared a common ancestor. Subsequent separate Ne trajectories are consistent with lineage splitting and speciation. All species showed evidence of population growth up until 100-200 thousand years ago (kya), followed by decline and then start of a new phase of population expansion. However, timing and amplitude of changes in Ne differed among species, and for pied flycatcher, the temporal dynamics of Ne differed between Spanish birds and central/northern European populations. This cautions against extrapolation of demographic inference between lineages and calls for adequate sampling to provide representative pictures of the coalescence process in different species or populations. We also empirically evaluate criteria for proper inference of demographic histories using PSMC and arrive at recommendations of using sequencing data with a mean genome coverage of ≥18X, a per-site filter of ≥10 reads and no more than 25% of missing data.

Increasing the power of genome wide association studies in natural populations using repeated measures - evaluation and implementation

Genomewide association studies (GWAS) enable detailed dissections of the genetic basis for organisms' ability to adapt to a changing environment. In long‐term studies of natural populations, individuals are often marked at one point in their life and then repeatedly recaptured. It is therefore essential that a method for GWAS includes the process of repeated sampling. In a GWAS, the effects of thousands of single‐nucleotide polymorphisms (SNPs) need to be fitted and any model development is constrained by the computational requirements. A method is therefore required that can fit a highly hierarchical model and at the same time is computationally fast enough to be useful.

Our method fits fixed SNP effects in a linear mixed model that can include both random polygenic effects and permanent environmental effects. In this way, the model can correct for population structure and model repeated measures. The covariance structure of the linear mixed model is first estimated and subsequently used in a generalized least squares setting to fit the SNP effects. The method was evaluated in a simulation study based on observed genotypes from a long‐term study of collared flycatchers in Sweden.

The method we present here was successful in estimating permanent environmental effects from simulated repeated measures data. Additionally, we found that especially for variable phenotypes having large variation between years, the repeated measurements model has a substantial increase in power compared to a model using average phenotypes as a response.

The method is available in the r package RepeatABEL. It increases the power in GWAS having repeated measures, especially for long‐term studies of natural populations, and the R implementation is expected to facilitate modelling of longitudinal data for studies of both animal and human populations.

Resolving Evolutionary Relationships in Closely Related Species with Whole-Genome Sequencing Data

Using genetic data to resolve the evolutionary relationships of species is of major interest in evolutionary and systematic biology. However, reconstructing the sequence of speciation events, the so-called species tree, in closely related and potentially hybridizing species is very challenging. Processes such as incomplete lineage sorting and interspecific gene flow result in local gene genealogies that differ in their topology from the species tree, and analyses of few loci with a single sequence per species are likely to produce conflicting or even misleading results. To study these phenomena on a full phylogenomic scale, we use whole-genome sequence data from 200 individuals of four black-and-white flycatcher species with so far unresolved phylogenetic relationships to infer gene tree topologies and visualize genome-wide patterns of gene tree incongruence. Using phylogenetic analysis in nonoverlapping 10-kb windows, we show that gene tree topologies are extremely diverse and change on a very small physical scale. Moreover, we find strong evidence for gene flow among flycatcher species, with distinct patterns of reduced introgression on the Z chromosome. To resolve species relationships on the background of widespread gene tree incongruence, we used four complementary coalescent-based methods for species tree reconstruction, including complex modeling approaches that incorporate post-divergence gene flow among species. This allowed us to infer the most likely species tree with high confidence. Based on this finding, we show that regions of reduced effective population size, which have been suggested as particularly useful for species tree inference, can produce positively misleading species tree topologies. Our findings disclose the pitfalls of using loci potentially under selection as phylogenetic markers and highlight the potential of modeling approaches to disentangle species relationships in systems with large effective population sizes and post-divergence gene flow.

Climate adaptation and speciation: particular focus on reproductive barriers in Ficedula flycatchers

Climate adaptation is surprisingly rarely reported as a cause for the build‐up of reproductive isolation between diverging populations. In this review, we summarize evidence for effects of climate adaptation on pre‐ and postzygotic isolation between emerging species with a particular focus on pied (Ficedula hypoleuca) and collared (Ficedula albicollis) flycatchers as a model for research on speciation. Effects of climate adaptation on prezygotic isolation or extrinsic selection against hybrids have been documented in several taxa, but the combined action of climate adaptation and sexual selection is particularly well explored in Ficedula flycatchers. There is a general lack of evidence for divergent climate adaptation causing intrinsic postzygotic isolation. However, we argue that the profound effects of divergence in climate adaptation on the whole biochemical machinery of organisms and hence many underlying genes should increase the likelihood of genetic incompatibilities arising as side effects. Fast temperature‐dependent co‐evolution between mitochondrial and nuclear genomes may be particularly likely to lead to hybrid sterility. Thus, how climate adaptation relates to reproductive isolation is best explored in relation to fast‐evolving barriers to gene flow, while more research on later stages of divergence is needed to achieve a complete understanding of climate‐driven speciation.

Sympatric divergence and clinal variation in multiple coloration traits of Ficedula flycatchers

Geographic variation in phenotypes plays a key role in fundamental evolutionary processes such as local adaptation, population differentiation and speciation, but the selective forces behind it are rarely known. We found support for the hypothesis that geographic variation in plumage traits of the pied flycatcher Ficedula hypoleuca is explained by character displacement with the collared flycatcher Ficedula albicollis in the contact zone. The plumage traits of the pied flycatcher differed strongly from the more conspicuous collared flycatcher in a sympatric area but increased in conspicuousness with increasing distance to there. Phenotypic differentiation (PST ) was higher than that in neutral genetic markers (FST ), and the effect of geographic distance remained when statistically controlling for neutral genetic differentiation. This suggests that a cline created by character displacement and gene flow explains phenotypic variation across the distribution of this species. The different plumage traits of the pied flycatcher are strongly to moderately correlated, indicating that they evolve non-independently from each other. The flycatchers provide an example of plumage patterns diverging in two species that differ in several aspects of appearance. The divergence in sympatry and convergence in allopatry in these birds provide a possibility to study the evolutionary mechanisms behind the highly divergent avian plumage patterns.

A high-density linkage map enables a second-generation collared flycatcher genome assembly and reveals the patterns of avian recombination rate variation and chromosomal evolution

Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high-density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best-order map contained 4215 markers, with a total distance of 3132 cm and a mean genetic distance between markers of 0.12 cm. Facilitated by the array being designed to include markers from most scaffolds, we obtained a second-generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super-scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event (6.6–7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cm/Mb and correlated closely with chromosome size, from 2 cm/Mb for chromosomes >100 Mb to >10 cm/Mb for chromosomes <10 Mb. This size dependence seemed entirely due to an obligate recombination event per chromosome; if 50 cm was subtracted from the genetic lengths of chromosomes, the rate per physical unit DNA was constant across chromosomes. Flycatcher recombination rate showed similar variation along chromosomes as chicken but lacked the large interior recombination deserts characteristic of zebra finch chromosomes.

Estimation of linkage disequilibrium and interspecific gene flow in Ficedula flycatchers by a newly developed 50k single-nucleotide polymorphism array

With the access to draft genome sequence assemblies and whole-genome resequencing data from population samples, molecular ecology studies will be able to take truly genome-wide approaches. This now applies to an avian model system in ecological and evolutionary research: Old World flycatchers of the genus Ficedula, for which we recently obtained a 1.1 Gb collared flycatcher genome assembly and identified 13 million single-nucleotide polymorphism (SNP)s in population resequencing of this species and its sister species, pied flycatcher. Here, we developed a custom 50K Illumina iSelect flycatcher SNP array with markers covering 30 autosomes and the Z chromosome. Using a number of selection criteria for inclusion in the array, both genotyping success rate and polymorphism information content (mean marker heterozygosity = 0.41) were high. We used the array to assess linkage disequilibrium (LD) and hybridization in flycatchers. Linkage disequilibrium declined quickly to the background level at an average distance of 17 kb, but the extent of LD varied markedly within the genome and was more than 10-fold higher in ‘genomic islands’ of differentiation than in the rest of the genome. Genetic ancestry analysis identified 33 F₁ hybrids but no later-generation hybrids from sympatric populations of collared flycatchers and pied flycatchers, contradicting earlier reports of backcrosses identified from much fewer number of markers. With an estimated divergence time as recently as <1 Ma, this suggests strong selection against F₁ hybrids and unusually rapid evolution of reproductive incompatibility in an avian system.

Malaria infections reinforce competitive asymmetry between two Ficedula flycatchers in a recent contact zone

Parasites may influence the outcome of interspecific competition between closely related host species through lower parasite virulence in the host with which they share the longer evolutionary history. We tested this idea by comparing the prevalence of avian malaria (Haemosporidia) lineages and their association with survival in pied and collared flycatchers (Ficedula hypoleuca and F. albicollis) breeding in a recent contact zone on the Swedish island of Öland. A nested PCR protocol amplifying haemosporidian fragments of mtDNA was used to screen the presence of malaria lineages in 1048 blood samples collected during 6 years. Competitively inferior pied flycatchers had a higher prevalence of blood parasites, including the lineages that were shared between the two flycatcher species. Multistate mark-recapture models revealed a lower survival of infected versus uninfected female pied flycatchers, while no such effects were detected in male pied flycatchers or in collared flycatchers of either sex. Our results show that a comparatively new host, the collared flycatcher, appears to be less susceptible to a local northern European malarial lineage where the collared flycatchers have recently expanded their distribution. Pied flycatchers experience strong reproductive interference from collared flycatchers, and the additional impact of species-specific blood parasite effects adds to this competitive exclusion. These results support the idea that parasites can strongly influence the outcome of interspecific competition between closely related host species, but that the invading species need not necessarily be more susceptible to local parasites.