Search

Sperm morphology varies considerably among species. Sperm traits may contribute to speciation if they diverge fast in allopatry and cause conspecific sperm precedence upon secondary contact. However, their role in driving prezygotic isolation has been poorly investigated. Here we test the hypothesis that, early in the speciation process, female promiscuity promotes a reduction in overlap in sperm length distributions among songbird populations. We assembled a data set of 20 pairs of populations with known sperm length distributions, a published estimate of divergence time, and an index of female promiscuity derived from extrapair paternity rates or relative testis size. We found that sperm length distributions diverged more rapidly in more promiscuous species. Faster divergence between sperm length distributions was caused by the lower variance in the trait in more promiscuous species, and not by faster divergence of the mean sperm lengths. The reduced variance is presumably due to stronger stabilizing selection on sperm length mediated by sperm competition. If divergent sperm length optima in allopatry causes conspecific sperm precedence in sympatry, which remains to be shown empirically, female promiscuity may promote prezygotic isolation and rapid speciation in songbirds., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE).)

The intertidal gastropod Littorina saxatilis is a model system to study speciation and local adaptation. The repeated occurrence of distinct ecotypes showing different levels of genetic divergence makes L. saxatilis particularly suited to study different stages of the speciation continuum in the same lineage. A major finding is the presence of several large chromosomal inversions associated with the divergence of ecotypes and, specifically, the species offers a system to study the role of inversions in this divergence. The genome of L. saxatilis is 1.35 Gb and composed of 17 chromosomes. The first reference genome of the species was assembled using Illumina data, was highly fragmented (N50 of 44 kb), and was quite incomplete, with a BUSCO completeness of 80.1% on the Metazoan dataset. A linkage map of one full-sibling family enabled the placement of 587 Mbp of the genome into 17 linkage groups corresponding to the haploid number of chromosomes, but the fragmented nature of this reference genome limited the understanding of the interplay between divergent selection and gene flow during ecotype formation. Here, we present a newly generated reference genome that is highly contiguous, with a N50 of 67 Mb and 90.4% of the total assembly length placed in 17 super-scaffolds. It is also highly complete with a BUSCO completeness of 94.1% of the Metazoa dataset. This new reference will allow for investigations into the genomic regions implicated in ecotype formation as well as better characterization of the inversions and their role in speciation., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Chromosomal inversions have been identified in many natural populations and can be responsible for novel traits and rapid adaptation. In zebra finch, a large region on the Z chromosome has been subject to multiple inversions, which have pleiotropic effects on multiple traits but especially on sperm phenotypes, such as midpiece and flagellum length. To understand the effect, the Z inversion has on these traits, we examined testis and liver transcriptomes of young males at different maturation times. We compared gene expression differences among three inversion karyotypes: AA, B*B* and AB*, where B* denotes the inverted regions on Z with respect to A. In testis, 794 differentially expressed genes were found and most of them were located on chromosome Z. They were functionally enriched for sperm-related traits. We also identified clusters of co-expressed genes that matched with the inversion-related sperm phenotypes. In liver, there were some enriched functions and some overrepresentation on chromosome Z with similar location as in testis. In both tissues, the overrepresented genes were located near the distal end of Z but also in the middle of the chromosome. For the heterokaryotype, we observed several genes with one allele being dominantly expressed, similar to expression patterns in one or the other homokaryotype. This was confirmed with SNPs for three genes, and interestingly one gene, DMGDH, had allele-specific expression originating mainly from one inversion haplotype in the testis, yet both inversion haplotypes were expressed equally in the liver. This karyotype-specific difference in tissue-specific expression suggests a pleiotropic effect of the inversion and thus suggests a mechanism for divergent phenotypic effects resulting from an inversion., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

When populations become geographically isolated, they begin to diverge in various traits and at variable rates. The dynamics of such trait divergences are relevant for understanding evolutionary processes such as local adaptation and speciation. Here we examine divergences in sperm and body structures in a polygynandrous songbird, the alpine accentor (Prunella collaris) between two allopatric high-altitude populations, in Morocco and Spain. The populations diverged around 82,000 years ago, as estimated with a coalescence-based phylogenetic analysis of genome-wide single-nucleotide polymorphisms. We found that birds in the two areas had nonoverlapping sperm lengths, which suggests adaptation to divergent female reproductive tract environments. Sperm length also showed an exceptionally low coefficient of among-male variation, a signal of strong stabilizing selection imposed by sperm competition. The evolutionary rate of sperm length was almost twice the rates for the most divergent morphological traits and more than three times higher than expected from literature data over a similar generational timescale. This rapid evolution of a key reproductive trait has implications for reproductive isolation and ultimately for speciation. Strong selection for different sperm length optima in allopatry predicts conspecific sperm precedence and disruptive selection in sympatry, hence a possible postcopulatory prezygotic barrier to gene flow., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE).)

The Centre for Marine Evolutionary Biology (CeMEB) at the University of Gothenburg, Sweden, was established in 2008 through a 10-year research grant of 8.7 m€ to a team of senior researchers. Today, CeMEB members have contributed >500 scientific publications, 30 PhD theses and have organised 75 meetings and courses, including 18 three-day meetings and four conferences. What are the footprints of CeMEB, and how will the centre continue to play a national and international role as an important node of marine evolutionary research? In this perspective article, we first look back over the 10 years of CeMEB activities and briefly survey some of the many achievements of CeMEB. We furthermore compare the initial goals, as formulated in the grant application, with what has been achieved, and discuss challenges and milestones along the way. Finally, we bring forward some general lessons that can be learnt from a research funding of this type, and we also look ahead, discussing how CeMEB's achievements and lessons can be used as a springboard to the future of marine evolutionary biology., Competing Interests: We declare we have no competing interests., (© 2022 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Species invasions are a global problem of increasing concern, especially in highly connected aquatic environments. Despite this, salinity conditions can pose physiological barriers to their spread, and understanding them is important for management. In Scandinavia's largest cargo port, the invasive round goby ( Neogobius melanostomus ) is established across a steep salinity gradient. We used 12,937 SNPs to identify the genetic origin and diversity of three sites along the salinity gradient and round goby from western, central and northern Baltic Sea, as well as north European rivers. Fish from two sites from the extreme ends of the gradient were also acclimated to freshwater and seawater, and tested for respiratory and osmoregulatory physiology. Fish from the high-salinity environment in the outer port showed higher genetic diversity, and closer relatedness to the other regions, compared to fish from lower salinity upstream the river. Fish from the high-salinity site also had higher maximum metabolic rate, fewer blood cells and lower blood Ca 2+ . Despite these genotypic and phenotypic differences, salinity acclimation affected fish from both sites in the same way: seawater increased the blood osmolality and Na + levels, and freshwater increased the levels of the stress hormone cortisol. Our results show genotypic and phenotypic differences over short spatial scales across this steep salinity gradient. These patterns of the physiologically robust round goby are likely driven by multiple introductions into the high-salinity site, and a process of sorting, likely based on behaviour or selection, along the gradient. This euryhaline fish risks spreading from this area, and seascape genomics and phenotypic characterization can inform management strategies even within an area as small as a coastal harbour inlet., Competing Interests: The authors declare no conflict of interest., (© 2022 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

In species with alternative reproductive tactics, there is much empirical support that parasitically spawning males have larger testes and greater sperm numbers as an evolved response to a higher degree of sperm competition, but support for higher sperm performance (motility, longevity and speed) by such males is inconsistent. We used the sand goby ( Pomatoschistus minutus ) to test whether sperm performance differed between breeding-coloured males (small testes, large mucus-filled sperm-duct glands; build nests lined with sperm-containing mucus, provide care) and parasitic sneaker-morph males (no breeding colouration, large testes, rudimentary sperm-duct glands; no nest, no care). We compared motility (per cent motile sperm), velocity, longevity of sperm, gene expression of testes and sperm morphometrics between the two morphs. We also tested if sperm-duct gland contents affected sperm performance. We found a clear difference in gene expression of testes between the male morphs with 109 transcripts differentially expressed between the morphs. Notably, several mucin genes were upregulated in breeding-coloured males and two ATP-related genes were upregulated in sneaker-morph males. There was a partial evidence of higher sperm velocity in sneaker-morph males, but no difference in sperm motility. Presence of sperm-duct gland contents significantly increased sperm velocity, and nonsignificantly tended to increase sperm motility, but equally so for the two morphs. The sand goby has remarkably long-lived sperm, with only small or no decline in motility and velocity over time (5 min vs. 22 h), but again, this was equally true for both morphs. Sperm length (head, flagella, total and flagella-to-head ratio) did not differ between morphs and did not correlate with sperm velocity for either morph. Thus, other than a clear difference in testes gene expression, we found only modest differences between the two male morphs, confirming previous findings that increased sperm performance as an adaptation to sperm competition is not a primary target of evolution., Competing Interests: The authors declare no conflict of interest., (© 2022 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

For externally fertilising organisms in the aquatic environment, the abiotic fertilisation medium can be a strong selecting force. Among bony fishes, sperm are adapted to function in a narrow salinity range. A notable exception is the family Gobiidae, where several species reproduce across a wide salinity range. The family also contains several wide-spread invasive species. To better understand how these fishes tolerate such varying conditions, we measured sperm performance in relation to salinity from a freshwater and a brackish population within their ancestral Ponto-Caspian region of the round goby, Neogobius melanostomus . These two ancestral populations were then compared to nine additional invaded sites across northern Europe, both in terms of their sperm traits and by using genomic SNP markers. Our results show clear patterns of ancestral adaptations to freshwater and brackish salinities in their sperm performance. Population genomic analyses show that the ancestral ecotypes have generally established themselves in environments that fit their sperm adaptations. Sites close to ports with intense shipping show that both outbreeding and admixture can affect the sperm performance of a population in a given salinity. Rapid adaptation to local conditions is also supported at some sites. Historical and contemporary evolution in the traits of the round goby sperm cells is tightly linked to the population and seascape genomics as well as biogeographic processes in these invasive fishes. Since the risk of a population establishing in an area is related to the genotype by environment match, port connectivity and the ancestry of the round goby population can likely be useful for predicting the species spread.

Temperature is a key environmental parameter affecting both the phenotypes and distributions of organisms, particularly ectotherms. Rapid organismal responses to thermal environmental changes have been described for several ectotherms; however, the underlying molecular mechanisms often remain unclear. Here, we studied whole genome cytosine methylation patterns of European grayling ( Thymallus thymallus ) embryos from five populations with contemporary adaptations of early life history traits at either 'colder' or 'warmer' spawning grounds. We reared fish embryos in a common garden experiment using two temperatures that resembled the 'colder' and 'warmer' conditions of the natal natural environments. Genome-wide methylation patterns were similar in populations originating from colder thermal origin subpopulations, whereas single nucleotide polymorphisms uncovered from the same data identified strong population structure among isolated populations, but limited structure among interconnected populations. This was surprising because the previously studied gene expression response among populations was mostly plastic, and mainly influenced by the developmental temperature. These findings support the hypothesis of the magnified role of epigenetic mechanisms in modulating plasticity. The abundance of consistently changing methylation loci between two warmer-to-colder thermal origin population pairs suggests that local adaptation has shaped the observed methylation patterns. The dynamic nature of the methylomes was further highlighted by genome-wide and site-specific plastic responses. Our findings support both the presence of a plastic response in a subset of CpG loci, and the evolutionary role of methylation divergence between populations adapting to contrasting thermal environments.

Studies of colonization of new habitats that appear from rapidly changing environments are interesting and highly relevant to our understanding of divergence and speciation. Here, we analyse phenotypic and genetic variation involved in the successful establishment of a marine fish (sand goby, Pomatoschistus minutus) over a steep salinity drop from 35 PSU in the North Sea (NE Atlantic) to two PSU in the inner parts of the post-glacial Baltic Sea. We first show that populations are adapted to local salinity in a key reproductive trait, the proportion of motile sperm. Thereafter, we show that genome variation at 22,190 single nucleotide polymorphisms (SNPs) shows strong differentiation among populations along the gradient. Sequences containing outlier SNPs and transcriptome sequences, mapped to a draft genome, reveal associations with genes with relevant functions for adaptation in this environment but without overall evidence of functional enrichment. The many contigs involved suggest polygenic differentiation. We trace the origin of this differentiation using demographic modelling and find the most likely scenario is that at least part of the genetic differentiation is older than the Baltic Sea and is a result of isolation of two lineages prior to the current contact over the North Sea-Baltic Sea transition zone., (© 2020 The Authors. Journal of Evolutionary Biology published by John Wiley & Sons Ltd on behalf of European Society for Evolutionary Biology.)

In recent years, the field of sexual selection has exploded, with advances in theoretical and empirical research complementing each other in exciting ways. This perspective piece is the product of a "stock-taking" workshop on sexual selection and sexual conflict. Our aim is to identify and deliberate on outstanding questions and to stimulate discussion rather than provide a comprehensive overview of the entire field. These questions are organized into four thematic sections we deem essential to the field. First we focus on the evolution of mate choice and mating systems. Variation in mate quality can generate both competition and choice in the opposite sex, with implications for the evolution of mating systems. Limitations on mate choice may dictate the importance of direct vs. indirect benefits in mating decisions and consequently, mating systems, especially with regard to polyandry. Second, we focus on how sender and receiver mechanisms shape signal design. Mediation of honest signal content likely depends on integration of temporally variable social and physiological costs that are challenging to measure. We view the neuroethology of sensory and cognitive receiver biases as the main key to signal form and the 'aesthetic sense' proposed by Darwin. Since a receiver bias is sufficient to both initiate and drive ornament or armament exaggeration, without a genetically correlated or even coevolving receiver, this may be the appropriate 'null model' of sexual selection. Thirdly, we focus on the genetic architecture of sexually selected traits. Despite advances in modern molecular techniques, the number and identity of genes underlying performance, display and secondary sexual traits remains largely unknown. In-depth investigations into the genetic basis of sexual dimorphism in the context of long-term field studies will reveal constraints and trajectories of sexually selected trait evolution. Finally, we focus on sexual selection and conflict as drivers of speciation. Population divergence and speciation are often influenced by an interplay between sexual and natural selection. The extent to which sexual selection promotes or counteracts population divergence may vary depending on the genetic architecture of traits as well as the covariance between mating competition and local adaptation. Additionally, post-copulatory processes, such as selection against heterospecific sperm, may influence the importance of sexual selection in speciation. We propose that efforts to resolve these four themes can catalyze conceptual progress in the field of sexual selection, and we offer potential avenues of research to advance this progress., Competing Interests: Scott Edwards is an Academic Editor for PeerJ., (©2019 Lindsay et al.)

The climate-change-driven increase in temperature is occurring rapidly and decreasing the predictability of seasonal rhythms at high latitudes. It is therefore urgent to understand how a change in the relationship between photoperiod and temperature can affect ectotherms in these environments. We tested whether temperature affects daily rhythms of transcription in a cold-adapted salmonid using high-throughput RNA sequencing. Arctic char ( Salvelinus alpinus ) from a subarctic population were reared at a high and a low temperature (15 and 8°C) for 1 month under natural, decreasing day length during late summer. Liver transcriptomes were compared between samples collected in the middle and towards the end of the light period and in the middle of the dark period. Daily variation in transcription was lower in fish from the low temperature compared with strong daily variation in warm-acclimated fish, suggesting that cold temperatures dampen the cycling of transcriptional rhythms under a simultaneously decreasing day length. Different circadian clock genes had divergent expression patterns, responding either by decreased expression or by increased rhythmicity at 15°C compared with 8°C. The results point out mechanisms that can affect the ability of fish to adapt to increasing temperatures caused by climate change., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2018. Published by The Company of Biologists Ltd.)

The distributions of species are not only determined by where they can survive - they must also be able to reproduce. Although immigrant inviability is a well-established concept, the fact that immigrants also need to be able to effectively reproduce in foreign environments has not been fully appreciated in the study of adaptive divergence and speciation. Fertilization and reproduction are sensitive life-history stages that could be detrimentally affected for immigrants in non-native habitats. We propose that "immigrant reproductive dysfunction" is a hitherto overlooked aspect of reproductive isolation caused by natural selection on immigrants. This idea is supported by results from experiments on an externally fertilizing fish (sand goby, Pomatoschistus minutus). Growth and condition of adults were not affected by non-native salinity whereas males spawning as immigrants had lower sperm motility and hatching success than residents. We interpret these results as evidence for local adaptation or acclimation of sperm, and possibly also components of paternal care. The resulting loss in fitness, which we call "immigrant reproductive dysfunction," has the potential to reduce gene flow between populations with locally adapted reproduction, and it may play a role in species distributions and speciation., (© 2017 The Author(s). Evolution © 2017 The Society for the Study of Evolution.)

Much adaptive evolutionary change is underlain by mutational variation in regions of the genome that regulate gene expression rather than in the coding regions of the genes themselves. An understanding of the role of gene expression variation in facilitating local adaptation will be aided by an understanding of underlying regulatory networks. Here, we characterize the genetic architecture of gene expression variation in the threespine stickleback (Gasterosteus aculeatus), an important model in the study of adaptive evolution. We collected transcriptomic and genomic data from 60 half-sib families using an expression microarray and genotyping-by-sequencing, and located expression quantitative trait loci (eQTL) underlying the variation in gene expression in liver tissue using an interval mapping approach. We identified eQTL for several thousand expression traits. Expression was influenced by polymorphism in both cis- and trans-regulatory regions. Trans-eQTL clustered into hotspots. We did not identify master transcriptional regulators in hotspot locations: rather, the presence of hotspots may be driven by complex interactions between multiple transcription factors. One observed hotspot colocated with a QTL recently found to underlie salinity tolerance in the threespine stickleback. However, most other observed hotspots did not colocate with regions of the genome known to be involved in adaptive divergence between marine and freshwater habitats., (Copyright © 2017 Pritchard et al.)

Phthalate esters are plasticizers frequently found in wastewater effluents. Previous studies on phthalates have reported anti-androgenic activity in mammals, causing concerns of their potential effects on the reproduction of aquatic organisms. Another group of environmental endocrine disrupters, steroidal estrogens, are known to inhibit steroid biosynthesis in the gonads, but the effects related to spermatogenesis are not well understood in fish. In this study, three-spined sticklebacks were exposed to di-n-butyl phthalate (DBP) and 17α ethinyl-oestradiol (EE2) at nominal concentrations 35μg/L and 40ng/L, respectively, for four days. The aim of the study was to obtain insight into the acute transcriptional responses putatively associated with endocrine disruption. RNA samples from eight individual male fish per treatment (including controls) were used in microarray analysis, covering the expression of approximately 21,000 genes. In the EE2 treatment the results show transcriptional downregulation of genes associated with steroid biosynthesis pathway and up-regulation of genes involved in pathways related to epidermal growth factor signaling and xenobiotic metabolism. The transcriptional response to DBP was in general weaker than to EE2, but based on enrichment analysis, we suggest adverse effects on retinoid metabolism, creatine kinase activity and cell adhesion. Among the genes showing highest fold changes after DBP treatment compared to control was the teleost fish -specific cytochrome P450 17A2. Overall, this study promotes our understanding on molecular responses to anti-androgens and estrogens in fish testes., (Copyright © 2015 Elsevier Inc. All rights reserved.)

In this study, we explored the hypothesis that killifish acclimate to thermal extremes through regulation of genes involved in stress and metabolism. We examined the liver and gonadal transcription of heat shock proteins (hsp70, hsp90a, hsp90b), glucokinase (gck), and high mobility group b1 (hmgb1) protein in wild killifish species from hot springs and rivers using quantitative real-time PCR. Moreover, we exposed a river killifish species to a long-term thermal regime of hot spring (37-40°C) and examined the liver transcription of the heat shock genes. Our results showed that hot spring killifish showed a significant, strong upregulation of liver hsp90a. Moreover, the testicular transcript levels of hsp90a, hsp90b, and hsp70 were higher in hot spring killifish than the river ones. The results of the common garden experiments showed that the transcripts of hsp70, hsp90b, and hmgb1 were mildly induced (> twofold) at the time when temperature reached to 37-40°C, while the transcripts of hsp90a were strongly induced (17-fold increase). The level of hsp90a was dramatically more upregulated when fish were maintained in thermal extreme (42-fold change higher than in ambient temperature). Moreover, a significant downregulation of gck transcripts was observed at the time when temperature was raised to 37-40°C (80-fold decrease) and during exposure to long-term thermal extreme (56-fold decrease). It can be concluded that the regulation of heat shock genes particularly hsp90a might be a key factor of the acclimation of fish to high temperature environments like hot springs., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Understanding how populations adapt to changing environmental conditions is a long-standing theme in evolutionary biology. Gene expression changes have been recognized as an important driver of local adaptation, but relatively little is known regarding the direction of change and in particular, about the interplay between plastic and evolutionary gene expression. We have previously shown that the gene expression profiles of European grayling (Thymallus thymallus) populations inhabiting different thermal environments include both plastic and evolutionary components. However, whether the plastic and evolutionary responses were in the same direction was not investigated in detail, nor was the identity of the specific genes involved. In this study, we show that the plastic changes in protein expression in response to different temperatures are highly correlated with the evolutionary response in grayling subpopulations adapted to different thermal environments. This finding provides preliminary evidence that the plastic response most likely facilitates adaptation during the early phases of colonization of thermal environments. The proteins that showed significant changes in expression level between warm and cold temperature treatments were mostly related to muscle development, which is consistent with earlier findings demonstrating muscle mass differentiation between cold and warm grayling populations., (© The American Genetic Association. 2015. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Evidence implicating differential gene expression as a significant driver of evolutionary novelty continues to accumulate, but our understanding of the underlying sources of variation in expression, both environmental and genetic, is wanting. Heritability in particular may be underestimated when inferred from genetic mapping studies, the predominant "genetical genomics" approach to the study of expression variation. Such uncertainty represents a fundamental limitation to testing for adaptive evolution at the transcriptomic level. By studying the inheritance of expression levels in 10,495 genes (10,527 splice variants) in a threespine stickleback pedigree consisting of 563 individuals, half of which were subjected to a thermal treatment, we show that 74-98% of transcripts exhibit significant additive genetic variance. Dominance variance is also prevalent (41-99% of transcripts), and genetic sources of variation seem to play a more significant role in expression variance in the liver than a key environmental variable, temperature. Among-population comparisons suggest that the majority of differential expression in the liver is likely due to neutral divergence; however, we also show that signatures of directional selection may be more prevalent than those of stabilizing selection. This predominantly aligns with the neutral model of evolution for gene expression but also suggests that natural selection may still act on transcriptional variation in the wild. As genetic variation both within- and among-populations ultimately defines adaptive potential, these results indicate that broad adaptive potential may be found within the transcriptome., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The loss of Y-linked genes during sex chromosome evolution creates a potentially deleterious low gene dosage in males. Recent studies have reported different strategies of dosage compensation. Unfortunately, most of these studies investigated taxa with comparatively old sex chromosome systems, which may limit insights into the evolution of dosage compensation and thus into the causes of different compensation strategies. Using deep RNA sequencing, we investigate differential expression patterns along the young XY chromosomes of threespine sticklebacks. Our strata-specific analyses provide new insights into the spatial patterns during the early stages of the evolution of dosage compensation. In particular, our results indicate systematic upregulation of male gene expression in stratum II, which in turn causes female hypertranscription in the same stratum. These findings are consistent with theoretical predictions that selection during early stages of sex chromosome evolution is stronger for a compensating upregulation in males than for the countercompensation of female hyperexpression. In contrast, no elevated gene expression is detectable in stratum I. We argue that strata-specific differences in compensating male gene expression may evolve in response to differences in the prevailing mechanism of Y chromosome degeneration., (© The Author(s) 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

A common challenge in phylogenetic reconstruction is to find enough suitable genomic markers to reliably trace splitting events with short internodes. Here, we present phylogenetic analyses based on genomewide single-nucleotide polymorphisms (SNPs) of an enigmatic avian radiation, the subspecies complex of Afrocanarian blue tits (Cyanistes teneriffae). The two sister species, the Eurasian blue tit (Cyanistes caeruleus) and the azure tit (Cyanistes cyanus), constituted the out-group. We generated a large data set of SNPs for analysis of population structure and phylogeny. We also adapted our protocol to utilize degraded DNA from old museum skins from Libya. We found strong population structuring that largely confirmed subspecies monophyly and constructed a coalescent-based phylogeny with full support at all major nodes. The results are consistent with a recent hypothesis that La Palma and Libya are relic populations of an ancient Afrocanarian blue tit, although a small data set for Libya could not resolve its position relative to La Palma. The birds on the eastern islands of Fuerteventura and Lanzarote are similar to those in Morocco. Together they constitute the sister group to the clade containing the other Canary Islands (except La Palma), in which El Hierro is sister to the three central islands. Hence, extant Canary Islands populations seem to originate from multiple independent colonization events. We also found population divergences in a key reproductive trait, viz. sperm length, which may constitute reproductive barriers between certain populations. We recommend a taxonomic revision of this polytypic species, where several subspecies should qualify for species rank., (© 2014 John Wiley & Sons Ltd.)