Ecology and genetics of speciation inFicedulaflycatchers

Trait Variation and Spatiotemporal Dynamics across Avian Secondary Contact Zones

A secondary contact zone (SCZ) is an area where incipient species or divergent populations may meet, mate, and hybridize. Due to the diverse patterns of interspecific hybridization, SCZs function as field labs for illuminating the on-going evolutionary processes of speciation and the establishment of reproductive isolation. Interspecific hybridization is widely present in avian populations, making them an ideal system for SCZ studies. This review exhaustively summarizes the variations in unique traits within avian SCZs (vocalization, plumage, beak, and migratory traits) and the various movement patterns of SCZs observed in previous publications. It also highlights several potential future research directions in the genomic era, such as the relationship between phenotypic and genomic differentiation in SCZs, the genomic basis of trait differentiation, SCZs shared by multiple species, and accurate predictive models for forecasting future movements under climate change and human disturbances. This review aims to provide a more comprehensive understanding of speciation processes and offers a theoretical foundation for species conservation.

A taste for the familiar: explaining the inbreeding paradox

The negative consequences of inbreeding have led animal biologists to assume that mate choice is generally biased against relatives. However, inbreeding avoidance is highly variable and by no means the rule across animal taxa. Even when inbreeding is costly, there are numerous examples of animals failing to avoid inbreeding or even preferring to mate with close kin. We argue that selective and mechanistic constraints interact to limit the evolution of inbreeding avoidance, notably when there is a risk of mating with heterospecifics and losing fitness through hybridization. Further, balancing inbreeding avoidance with conspecific mate preference may drive the evolution of multivariate sexual communication. Studying different social and sexual decisions within the same species can illuminate trade-offs among mate-choice mechanisms.

Assortative Mating in an Ecological Context: Effects of Mate Choice Errors and Relative Species Abundance on the Frequency and Asymmetry of Hybridization

AbstractThe frequency and asymmetry of mixed-species mating set the initial stage for the ecological and evolutionary implications of hybridization. How such patterns of mixed-species mating, in turn, are influenced by the combination of mate choice errors and relative species abundance remains largely unknown. We develop a mathematical model that generates predictions for how relative species abundances and mate choice errors affect hybridization patterns. When mate choice errors are small (<5%), the highest frequency of hybridization occurs when one of the hybridizing species is at low abundance, but when mate choice errors are high (>5%), the highest hybridization frequency occurs when species occur in equal proportions. Furthermore, females of the less abundant species are overrepresented in mixed-species matings. We compare our theoretical predictions with empirical data on naturally hybridizing Ficedula flycatchers and find that hybridization is highest when the two species occur in equal abundance, implying rather high mate choice errors. We discuss ecological and evolutionary implications of our findings and encourage future work on hybrid zone dynamics that take demographic aspects, such as relative species abundance, into account.

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.

Avian Population Genomics Taking Off: Latest Findings and Future Prospects

Birds are one of the most recognizable and diverse groups of organisms on earth. This group has played an important role in many fields, including the development of methods in behavioral ecology and evolutionary theory. The use of population genomics took off following the advent of high-throughput sequencing in various taxa. Several features of avian genomes make them particularly amenable for work in this field, including their nucleated red blood cells permitting easy DNA extraction and small, compact genomes. We review the latest findings in the population genomics of birds here, emphasizing questions related to behavior, ecology, evolution, and conservation. Additionally, we include insights in trait mapping and the ability to obtain accurate estimates of important summary statistics for conservation (e.g., genetic diversity and inbreeding). We highlight roadblocks that will need to be overcome in order to advance work on the population genomics of birds and prospects for future work. Roadblocks include the assembly of more contiguous reference genomes using long-reads and optical mapping. Prospects include the integration of population genomics with additional fields (e.g., landscape genetics, phylogeography, and genomic mapping) along with studies beyond genetic variants (e.g., epigenetics).

The Genomic Landscape of Divergence Across the Speciation Continuum in Island-Colonising Silvereyes (Zosterops lateralis)

Inferring the evolutionary dynamics at play during the process of speciation by analyzing the genomic landscape of divergence is a major pursuit in population genomics. However, empirical assessments of genomic landscapes under varying evolutionary scenarios that are known a priori are few, thereby limiting our ability to achieve this goal. Here we combine RAD-sequencing and individual-based simulations to evaluate the genomic landscape of divergence in the silvereye (Zosterops lateralis). Using pairwise comparisons that differ in divergence timeframe and the presence or absence of gene flow, we document how genomic patterns accumulate along the speciation continuum. In contrast to previous predictions, our results provide limited support for the idea that divergence accumulates around loci under divergent selection or that genomic islands widen with time. While a small number of genomic islands were found in populations diverging with and without gene flow, in few cases were SNPs putatively under selection tightly associated with genomic islands. The transition from localized to genome-wide levels of divergence was captured using individual-based simulations that considered only neutral processes. Our results challenge the ubiquity of existing verbal models that explain the accumulation of genomic differences across the speciation continuum and instead support the idea that divergence both within and outside of genomic islands is important during the speciation process.

Differential divergence in autosomes and sex chromosomes is associated with intra-island diversification at a very small spatial scale in a songbird lineage

Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome-wide markers to investigate the diversification of the Reunion grey white-eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping-by-sequencing and pooled RAD-seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z-linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process.

The Past and Future of Experimental Speciation

Speciation is the result of evolutionary processes that generate barriers to gene flow between populations, facilitating reproductive isolation. Speciation is typically studied via theoretical models and snapshot tests in natural populations. Experimental speciation enables real-time direct tests of speciation theory and has been long touted as a critical complement to other approaches. We argue that, despite its promise to elucidate the evolution of reproductive isolation, experimental speciation has been underutilised and lags behind other contributions to speciation research. We review recent experiments and outline a framework for how experimental speciation can be implemented to address current outstanding questions that are otherwise challenging to answer. Greater uptake of this approach is necessary to rapidly advance understanding of speciation.

Ecological divergence plays an important role in strong but complex reproductive isolation in campions (Silene)

New species arise through the evolution of reproductive barriers between formerly interbreeding lineages. Yet, comprehensive assessments of potential reproductive barriers, which are needed to make inferences on processes driving speciation, are only available for a limited number of systems. In this study, we estimated individual and cumulative strengths of seven prezygotic and six postzygotic reproductive barriers between the recently diverged taxa Silene dioica (L.) Clairv. and S. latifolia Poiret using both published and new data. A combination of multiple partial reproductive barriers resulted in near-complete reproductive isolation between S. dioica and S. latifolia, consistent with earlier estimates of gene flow between the taxa. Extrinsic barriers associated with adaptive ecological divergence were most important, while intrinsic postzygotic barriers had moderate individual strength but contributed only little to total reproductive isolation. These findings are in line with ecological divergence as driver of speciation. We further found extensive variation in extrinsic reproductive isolation, ranging from sites with very strong selection against migrants and hybrids to intermediate sites where substantial hybridization is possible. This situation may allow for, or even promote, heterogeneous genetic divergence.

Behavioral Isolation and Incipient Speciation in Birds

Behavioral changes, such as those involved in mating, foraging, and migration, can generate reproductive barriers between populations. Birds, in particular, are known for their great diversity in these behaviors, and so behavioral isolation is often proposed to be the major driver of speciation. Here, we review empirical evidence to evaluate the importance of behavioral isolation in the early stages of avian speciation. Experimentally measured mating preferences indicate that changes in mating behavior can result in premating barriers, with their strength depending on the extent of divergence in mating signals. Differences in migratory and foraging behavior also can play important roles in generating reproductive barriers in the early stages of speciation. However, because premating behavioral isolation is imperfect, extrinsic postzygotic barriers, in the form of selection against hybrids having intermediate phenotypes, also play an important role in avian diversification, especially in completing the speciation process.

Fine-scale geographic patterns of gene flow and reproductive character displacement in Drosophila subquinaria and Drosophila recens

When two species are incompletely isolated, strengthening premating isolation barriers in response to the production of low fitness hybrids may complete the speciation process. Here, we use the sister species Drosophila subquinaria and Drosophila recens to study the conditions under which this reinforcement of species boundaries occurs in natural populations. We first extend the region of known sympatry between these species, and then we conduct a fine-scale geographic survey of mate discrimination coupled with estimates of gene flow within and admixture between species. Within D. subquinaria, reinforcement is extremely effective: we find variation in mate discrimination both against D. recens males and against conspecific allopatric males on the scale of a few kilometres and in the face of gene flow both from conspecific populations and introgression from D. recens. In D. recens, we do not find evidence for increased mate discrimination in sympatry, even where D. recens is rare, consistent with substantial gene flow throughout the species' range. Finally, we find that introgression between species is asymmetric, with more from D. recens into D. subquinaria than vice versa. Within each species, admixture is highest in the geographic region where it is rare relative to the other species, suggesting that when hybrids are produced they are of low fitness. In sum, reinforcement within D. subquinaria is effective at maintaining species boundaries, but even when reinforcing selection is strong it may not always result in a pattern of strong reproductive character displacement due to variation in the frequency of hybridization and gene flow from neighbouring populations.

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.

Correlated patterns of genetic diversity and differentiation across an avian family

Comparative studies of closely related taxa can provide insights into the evolutionary forces that shape genome evolution and the prevalence of convergent molecular evolution. We investigated patterns of genetic diversity and differentiation in stonechats (genus Saxicola), a widely distributed avian species complex with phenotypic variation in plumage, morphology and migratory behaviour, to ask whether similar genomic regions have become differentiated in independent, but closely related, taxa. We used whole-genome pooled sequencing of 262 individuals from five taxa and found that levels of genetic diversity and divergence are strongly correlated among different stonechat taxa. We then asked whether these patterns remain correlated at deeper evolutionary scales and found that homologous genomic regions have become differentiated in stonechats and the closely related Ficedula flycatchers. Such correlation across a range of evolutionary divergence and among phylogenetically independent comparisons suggests that similar processes may be driving the differentiation of these independently evolving lineages, which in turn may be the result of intrinsic properties of particular genomic regions (e.g. areas of low recombination). Consequently, studies employing genome scans to search for areas important for reproductive isolation or adaptation should account for corresponding regions of differentiation, as these regions may not necessarily represent speciation islands or evidence of local adaptation.

Changing climate mediates sapsucker (Aves: Sphyrapicus) hybrid zone movement

Hybrid zones, where two divergent taxa meet and interbreed, offer unique opportunities to investigate how climate contributes to reproductive isolation between closely related taxa and how these taxa may respond to climatic changes. Red-naped (Sphyrapicus nuchalis) and Red-breasted (Sphyrapicus ruber) sapsuckers (Aves: Picidae) hybridize along a narrow contact zone that stretches from northern California to British Columbia. The hybrid zone between these species has been studied extensively for more than 100 years and represents an excellent system for investigations of the evolution of reproductive isolation. Shifts in the proportions of phenotypes at hybrid localities since 1910 that were inferred using specimens from museum collections were confirmed using species distribution models. We predicted the historical, current, and future distributions of parental and hybrid sapsuckers using Random Forests models to quantify how climate change is affecting hybrid zone movement in the Pacific Northwest. We found observed distribution shifts of parental sapsuckers were likely the result of climate change over the past 100 years, with these shifts predicted to continue for both sapsuckers over the next 80 years. We found Red-breasted Sapsuckers are predicted to continue to expand, while Red-naped Sapsuckers are predicted to contract substantially under future climate scenarios. As a result of the predicted changes, the amount of overlap in the distribution of these sapsuckers may decrease. Using hybrid phenotypes, we found the climate niche occupied by the hybrid zone is predicted to disappear under future conditions. The disappearance of this climate niche where the two parental species come into contact and hybridize may lead to a substantial reduction in genetic introgression. Understanding the impacts of global climate change on hybrid zones may help us to better understand how speciation has been shaped by climate in the past, as well as how evolution may respond to climate change in the future.

Haemosporidian blood parasite diversity and prevalence in the semi-collared flycatcher (Ficedula semitorquata) from the eastern Balkans

We present the first survey on prevalence and diversity of haemosporidian parasites of the genera Plasmodium and Haemoproteus in a poorly studied migratory passerine, the semi-collared flycatcher (Ficedula semitorquata). In total, 110 individuals were sampled during two breeding seasons in Eastern Bulgaria. We collected both blood samples for PCR identification and blood films for microscopic identification of haemosporidians. We found six distinctive parasite cyt b lineages present in the blood of the semi-collared flycatcher (three Haemoproteus and three Plasmodium). Two of the lineages, i.e. pWW3 of Plasmodium sp. and hCCF2 of Haemoproteus sp., are recorded for the first time in the family Muscicapidae. The overall prevalence ranged between 12.2 and 15.9% and we did not find co-infections. We hypothesize that the low prevalence of haemosporidians in this species might be linked to its small population size and relatively restricted geographical range.

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.

Recombination Rate Evolution and the Origin of Species

A recipe for dissolving incipient species into a continuum of phenotypes is to recombine their genetic material. Therefore, students of speciation have become increasingly interested in the mechanisms by which recombination between locally adapted lineages is reduced. Evidence abounds that chromosomal rearrangements, via their suppression of recombination during meiosis in hybrids, play a major role in adaptation and speciation. By contrast, genic modifiers of recombination rates have been largely ignored in studies of speciation. We show how both types of reduction in recombination rates facilitate divergence in the face of gene flow, including the early stages of adaptive divergence, the persistence of species after secondary contact, and reinforcement.

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.

Linked selection and recombination rate variation drive the evolution of the genomic landscape of differentiation across the speciation continuum of Ficedula flycatchers

Speciation is a continuous process during which genetic changes gradually accumulate in the genomes of diverging species. Recent studies have documented highly heterogeneous differentiation landscapes, with distinct regions of elevated differentiation (“differentiation islands”) widespread across genomes. However, it remains unclear which processes drive the evolution of differentiation islands; how the differentiation landscape evolves as speciation advances; and ultimately, how differentiation islands are related to speciation. Here, we addressed these questions based on population genetic analyses of 200 resequenced genomes from 10 populations of four Ficedula flycatcher sister species. We show that a heterogeneous differentiation landscape starts emerging among populations within species, and differentiation islands evolve recurrently in the very same genomic regions among independent lineages. Contrary to expectations from models that interpret differentiation islands as genomic regions involved in reproductive isolation that are shielded from gene flow, patterns of sequence divergence (d_xy and relative node depth) do not support a major role of gene flow in the evolution of the differentiation landscape in these species. Instead, as predicted by models of linked selection, genome-wide variation in diversity and differentiation can be explained by variation in recombination rate and the density of targets for selection. We thus conclude that the heterogeneous landscape of differentiation in Ficedula flycatchers evolves mainly as the result of background selection and selective sweeps in genomic regions of low recombination. Our results emphasize the necessity of incorporating linked selection as a null model to identify genome regions involved in adaptation and speciation.

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.

Evolutionary analysis of the female-specific avian W chromosome

The typically repetitive nature of the sex-limited chromosome means that it is often excluded from or poorly covered in genome assemblies, hindering studies of evolutionary and population genomic processes in non-recombining chromosomes. Here, we present a draft assembly of the non-recombining region of the collared flycatcher W chromosome, containing 46 genes without evidence of female-specific functional differentiation. Survival of genes during W chromosome degeneration has been highly non-random and expression data suggest that this can be attributed to selection for maintaining gene dose and ancestral expression levels of essential genes. Re-sequencing of large population samples revealed dramatically reduced levels of within-species diversity and elevated rates of between-species differentiation (lineage sorting), consistent with low effective population size. Concordance between W chromosome and mitochondrial DNA phylogenetic trees demonstrates evolutionary stable matrilineal inheritance of this nuclear–cytonuclear pair of chromosomes. Our results show both commonalities and differences between W chromosome and Y chromosome evolution.

The evolution of non-recombining chromosomes is poorly understood. Here, the authors sequence the collared flycatcher female-specific W chromosome and show nonrandom survival of genes during W chromosome degeneration which is due to selection for maintaining gene dose and expression levels of essential genes.

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.

Speciation in Western Scrub-Jays, Haldane's rule, and genetic clines in secondary contact

Background

Haldane’s Rule, the tendency for the heterogametic sex to show reduced fertility in hybrid crosses, can obscure the signal of gene flow in mtDNA between species where females are heterogametic. Therefore, it is important when studying speciation and species limits in female-heterogametic species like birds to assess the signature of gene flow in the nuclear genome as well. We studied introgression of microsatellites and mtDNA across a secondary contact zone between coastal and interior lineages of Western Scrub-Jays (Aphelocoma californica) to test for a signature of Haldane’s Rule: a narrower cline of introgression in mtDNA compared to nuclear markers.

Results

Our initial phylogeographic analysis revealed that there is only one major area of contact between coastal and interior lineages and identified five genetic clusters with strong spatial structuring: Pacific Slope, Interior US, Edwards Plateau (Texas), Northern Mexico, and Southern Mexico. Consistent with predictions from Haldane’s Rule, mtDNA showed a narrower cline than nuclear markers across a transect through the hybrid zone. This result is not being driven by female-biased dispersal because neutral diffusion analysis, which included estimates of sex-specific dispersal rates, also showed less diffusion of mtDNA. Lineage-specific plumage traits were associated with nuclear genetic profiles for individuals in the hybrid zone, indicating that these differences are under genetic control.

Conclusions

This study adds to a growing list of studies that support predictions of Haldane’s Rule using cline analysis of multiple loci of differing inheritance modes, although alternate hypotheses like selection on different mtDNA types cannot be ruled out. That Haldane’s Rule appears to be operating in this system suggests a measure of reproductive isolation between the Pacific Slope and interior lineages. Based on a variety of evidence from the phenotype, ecology, and genetics, we recommend elevating three lineages to species level: A. californica (Pacific Slope); A. woodhouseii (Interior US plus Edwards Plateau plus Northern Mexico); A. sumichrasti (Southern Mexico). The distinctive Edwards Plateau population in Texas, which was monophyletic in mtDNA except for one individual, should be studied in greater detail given habitat threat.

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.

Genetic, morphological, and acoustic evidence reveals lack of diversification in the colonization process in an island bird

Songbirds with recently (i.e., early Holocene) founded populations are suitable models for studying incipient differentiation in oceanic islands. On such systems each colonization event represents a different evolutionary episode that can be studied by addressing sets of diverging phenotypic and genetic traits. We investigate the process of early differentiation in the spectacled warbler (Sylvia conspicillata) in 14 populations separated by sea barriers from three Atlantic archipelagos and from continental regions spanning from tropical to temperate latitudes. Our approach involved the study of sexual acoustic signals, morphology, and genetic data. Mitochondrial DNA did not provide clear population structure. However, microsatellites analyses consistently identified two genetic groups, albeit without correspondence to subspecies classification and little correspondence to geography. Coalescent analyses showed significant evidence for gene flow between the two genetic groups. Discriminant analyses could not correctly assign morphological or acoustic traits to source populations. Therefore, although theory predicting that in isolated populations genetic, morphological, or acoustic traits can lead to radiation, we have strikingly failed to document differentiation on these attributes in a resident passerine throughout three oceanic archipelagos.

Genomics and the origin of species

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

Reproductive character displacement of epicuticular compounds and their contribution to mate choice in Drosophila subquinaria and Drosophila recens

Interactions between species can alter selection on sexual displays used in mate choice within species. Here we study the epicuticular pheromones of two Drosophila species that overlap partially in geographic range and are incompletely reproductively isolated. Drosophila subquinaria shows a pattern of reproductive character displacement against Drosophila recens, and partial behavioral isolation between conspecific sympatric versus allopatric populations, whereas D. recens shows no such variation in mate choice. First, using manipulative perfuming experiments, we show that females use pheromones as signals for mate discrimination both between species and among populations of D. subquinaria. Second, we show that patterns of variation in epicuticular compounds, both across populations and between species, are consistent with those previously shown for mating probabilities: pheromone compositions differ between populations of D. subquinaria that are allopatric versus sympatric with D. recens, but are similar across populations of D. recens regardless of overlap with D. subquinaria. We also identify differences in pheromone composition among allopatric regions of D. subquinaria. In sum, our results suggest that epicuticular compounds are key signals used by females during mate recognition, and that these traits have diverged among D. subquinaria populations in response to reinforcing selection generated by the presence of D. recens.

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.

Geographic variation in animal colour polymorphisms and its role in speciation

Polymorphic species, in which multiple variants coexist within a population, are often used as model systems in evolutionary biology. Recent research has been dominated by the hypothesis that polymorphism can be a precursor to speciation. To date, the majority of research regarding polymorphism and speciation has focused on whether polymorphism is maintained within a population or whether morphs within populations may diverge to form separate species (sympatric speciation); however, the geographical context of speciation in polymorphic systems is likely to be both diverse and complex. In this review, we draw attention to the geographic variation in morph composition and frequencies that characterises many, if not most polymorphic species. Recent theoretical and empirical developments suggest that such variation in the number, type and frequency of morphs present among populations can increase the probability of speciation. Thus, the geographical context of a polymorphism requires a greater research focus. Here, we review the prevalence, causes and evolutionary consequences of geographic variation in polymorphism in colour-polymorphic animal species. The prevalence and nature of geographic variation in polymorphism suggests that polymorphism may be a precursor to and facilitate speciation more commonly than appreciated previously. We argue that a better understanding of the processes generating geographic variation in polymorphism is vital to understanding how polymorphism can promote speciation.

Heterospecific female mimicry in Ficedula flycatchers

Mimicry is a widespread phenomenon. Vertebrate visual mimicry often operates in an intraspecific sexual context, with some males resembling conspecific females. Pied flycatcher (Ficedula hypoleuca) dorsal plumage varies from the ancestral black to female-like brown. Experimental studies have shown that conspecific and heterospecific (collared flycatcher, F. albicollis) individuals of both sexes respond, at least initially, to brown individuals as if they were female. We quantified the perceptual and biochemical differences between brown feathers and found that brown pied flycatcher males are indistinguishable from heterospecific, but not from conspecific, females in both aspects. To our knowledge, this is the first evidence of a visual mimetic signalling system in a sexual context where the model is heterospecific to the mimic. By only mimicking heterospecific females, brown pied flycatcher males can establish territories next to the more dominant collared flycatcher in sympatry, suffer less aggression by darker conspecifics in allopatry and preserve within-species sexual recognition throughout the breeding range. A closer look at the evolutionary history and ecology of these two species illustrates how such a mimetic system can evolve. Although likely rare, this phenomenon might not be unique to Ficedula flycatchers.

Demographic divergence history of pied flycatcher and collared flycatcher inferred from whole-genome re-sequencing data

Profound knowledge of demographic history is a prerequisite for the understanding and inference of processes involved in the evolution of population differentiation and speciation. Together with new coalescent-based methods, the recent availability of genome-wide data enables investigation of differentiation and divergence processes at unprecedented depth. We combined two powerful approaches, full Approximate Bayesian Computation analysis (ABC) and pairwise sequentially Markovian coalescent modeling (PSMC), to reconstruct the demographic history of the split between two avian speciation model species, the pied flycatcher and collared flycatcher. Using whole-genome re-sequencing data from 20 individuals, we investigated 15 demographic models including different levels and patterns of gene flow, and changes in effective population size over time. ABC provided high support for recent (mode 0.3 my, range <0.7 my) species divergence, declines in effective population size of both species since their initial divergence, and unidirectional recent gene flow from pied flycatcher into collared flycatcher. The estimated divergence time and population size changes, supported by PSMC results, suggest that the ancestral species persisted through one of the glacial periods of middle Pleistocene and then split into two large populations that first increased in size before going through severe bottlenecks and expanding into their current ranges. Secondary contact appears to have been established after the last glacial maximum. The severity of the bottlenecks at the last glacial maximum is indicated by the discrepancy between current effective population sizes (20,000-80,000) and census sizes (5-50 million birds) of the two species. The recent divergence time challenges the supposition that avian speciation is a relatively slow process with extended times for intrinsic postzygotic reproductive barriers to evolve. Our study emphasizes the importance of using genome-wide data to unravel tangled demographic histories. Moreover, it constitutes one of the first examples of the inference of divergence history from genome-wide data in non-model species.

Evidence for asymmetrical divergence-gene flow of nuclear loci, but not mitochondrial loci, between seabird sister species: blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies

Understanding the process of speciation requires understanding how gene flow influences divergence. Recent analyses indicate that divergence can take place despite gene flow and that the sex chromosomes can exhibit different levels of gene flow than autosomes and mitochondrial DNA. Using an eight marker dataset including autosomal, z-linked, and mitochondrial loci we tested the hypothesis that blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies diverged from their common ancestor with gene flow, paying specific attention to the differences in gene flow estimates from nuclear and mitochondrial markers. We found no gene flow at mitochondrial markers, but found evidence from the combined autosomal and z-linked dataset that blue-footed and Peruvian boobies experienced asymmetrical gene flow during or after their initial divergence, predominantly from Peruvian boobies into blue-footed boobies. This gene exchange may have occurred either sporadically between periods of allopatry, or regularly throughout the divergence process. Our results add to growing evidence that diverging species can remain distinct but exchange genes.

Relative performance of hybrid nestlings in Ficedula flycatchers: a translocation experiment

Ecological speciation predicts that hybrids should experience relatively low fitness in the local environments of their parental species. In this study, we performed a translocation experiment of nestling hybrids between collared and pied flycatchers into the nests of conspecific pairs of their parental species. Our aim was to compare the performance of hybrids with purebred nestlings. Nestling collared flycatchers are known to beg and grow faster than nestling pied flycatchers under favorable conditions, but to experience higher mortality than nestling pied flycatchers under food limitation. The experiment was performed relatively late in the breeding season when food is limited. If hybrid nestlings have an intermediate growth potential and begging intensity, we expected them to beg and grow faster, but also to experience lower survival than pied flycatchers. In comparison with nestling collared flycatchers, we expected them to beg and grow slower, but to survive better. We found that nestling collared flycatchers indeed begged significantly faster and experienced higher mortality than nestling hybrids. Moreover, nestling hybrids had higher weight and tended to beg faster than nestling pied flycatchers, but we did not detect a difference in survival between the latter two groups of nestlings. We conclude that hybrid Ficedula nestlings appear to have a better intrinsic adaptation to food limitation late in the breeding season compared with nestling collared flycatchers. We discuss possible implications for gene flow between the two species.

Fluctuating selection and immigration as determinants of the phenotypic composition of a population

It is important to identify the factors that affect the evolutionary potential of populations to respond to environmental changes. Such processes are for example the ones affecting the amount of heritable phenotypic variation in a population. We examined factors explaining the wide phenotypic variation in the genetically determined black-brown dorsal colouration of male pied flycatchers (Ficedula hypoleuca) during a period of >50 years in a northern European breeding population. We demonstrate that the temperature-dependent relative breeding success of brown males predicts the inter-annual change in the proportion of the brown male phenotype. The proportion of brown males also appears to reflect immigration from Central Europe, where the brown type prevails due to local selection pressure. Warm springs in northern Central Europe had a positive effect on the proportion of the brown phenotype in the north in the early part of the study period, which suggests prolonged migration in favourable conditions. However, the association between warm springs and a high proportion of brown males has weakened from the 1950s to the present, which may explain why the proportion of the brown males in our study area decreased by a third during the period 1954 to 2008. This is likely a result of decreasing population size in Central Europe. These results demonstrate that temporal variation in environmental conditions is maintaining variation in the pied flycatcher male phenotype. They also indicate that climate warming has the potential to change the population composition both through temperature-dependent selection and environmental factors affecting long-distance immigration.

Ecological speciation along an elevational gradient in a tropical passerine bird?

Local adaptation of populations along elevational gradients is well known, but conclusive evidence that such divergence has resulted in the origin of distinct species in parapatry remains lacking. We integrated morphological, vocal, genetic and behavioural data to test predictions pertaining to the hypothesis of parapatric ecological speciation associated with elevation in populations of a tropical montane songbird, the Grey-breasted Wood-wren (Henicorhina leucophrys: Troglodytidae), from the Sierra Nevada de Santa Marta, Colombia. We confirmed that two distinct populations exist along the elevational gradient. Phylogenetic analyses tentatively indicate that the two populations are not sister taxa, suggesting they did not differentiate from a single ancestor along the gradient, but rather resulted from separate colonization events. The populations showed marked divergence in morphometrics, vocalizations and genetic variation in mitochondrial and nuclear loci, and little to no evidence of hybridization. Individuals of both populations responded more strongly to their own local songs than to songs from another elevation. Although the two forms do not appear to have differentiated locally in parapatry, morphological and vocal divergence along the elevational gradient is consistent with adaptation, suggesting a possible link between adaptive evolution in morphology and songs and the origin of reproductive isolation via a behavioural barrier to gene flow. The adaptive value of phenotypic differences between populations requires additional study.