Space use and genetic structure do not maintain color polymorphism in a species with alternative behavioral strategies

Space use including territoriality and spatial arrangement within a population can reveal important information on the nature, dynamics, and evolutionary maintenance of alternative strategies in color polymorphic species. Despite the prevalence of color polymorphic species as model systems in evolutionary biology, the interaction between space use and genetic structuring of morphs within populations has rarely been examined. Here, we assess the spatial and genetic structure of male throat color morphs within a population of the tawny dragon lizard, Ctenophorus decresii. Male color morphs do not differ in morphology but differ in aggressive and antipredator behaviors as well as androgen levels. Despite these behavioral and endocrine differences, we find that color morphs do not differ in territory size, with their spatial arrangement being essentially random with respect to each other. There were no differences in genetic diversity or relatedness between morphs; however, there was significant, albeit weak, genetic differentiation between morphs, which was unrelated to geographic distance between individuals. Our results indicate potential weak barriers to gene flow between some morphs, potentially due to nonrandom pre‐ or postcopulatory mate choice or postzygotic genetic incompatibilities. However, space use, spatial structure, and nonrandom mating do not appear to be primary mechanisms maintaining color polymorphism in this system, highlighting the complexity and variation in alternative strategies associated with color polymorphism.     

Phylogeographic structure,demographic history and morph composition in a colour polymorphic lizard

In polymorphic species, population divergence in morph composition and frequency has the potential to promote speciation. We assessed the relationship between geographic variation in male throat colour polymorphism and phylogeographic structure in the tawny dragon lizard, Ctenophorus decresii. We identified four genetically distinct lineages, corresponding to two polymorphic lineages in the Northern Flinders Ranges and Southern Flinders Ranges/Olary Ranges regions respectively, and a monomorphic lineage in the Mt Lofty Ranges/Kangaroo Island region. The degree of divergence between these three lineages was consistent with isolation to multiple refugia during Pleistocene glacial cycles, whereas a fourth, deeply divergent (at the interspecific level) and monomorphic lineage was restricted to western New South Wales. The same four morphs occurred in both polymorphic lineages, although populations exhibited considerable variation in the frequency of morphs. By contrast, male throat coloration in the monomorphic lineages differed from each other and from the polymorphic lineages. Our results suggest that colour polymorphism has evolved once in the C. decresii species complex, with subsequent loss of polymorphism in the Mt Lofty Ranges/Kangaroo Island lineage. However, an equally parsimonious scenario, that polymorphism arose independently twice within C. decresii, could not be ruled out. We also detected evidence of a narrow contact zone with limited genotypic admixture between the polymorphic Olary Ranges and monomorphic Mt Lofty Ranges regions, yet no individuals of intermediate colour phenotype. Such genetic divergence and evidence for barriers to gene flow between lineages suggest incipient speciation between populations that differ in morph composition.     

Phenotypic divergence among threespine stickleback that differ in nuptial coloration

By studying systems in their earliest stages of differentiation, we can learn about the evolutionary forces acting within and among populations and how those forces could contribute to reproductive isolation. Such an understanding would help us to better discern and predict how selection leads to the maintenance of multiple morphs within a species, rather than speciation. The postglacial adaptive radiation of the threespine stickleback (Gasterosteus aculeatus) is one of the best‐studied cases of evolutionary diversification and rapid, repeated speciation. Following deglaciation, marine stickleback have continually invaded freshwater habitats across the northern hemisphere and established resident populations that diverged innumerable times from their oceanic ancestors. Independent freshwater colonization events have yielded broadly parallel patterns of morphological differences in freshwater and marine stickleback. However, there is also much phenotypic diversity within and among freshwater populations. We studied a lesser‐known freshwater “species pair” found in southwest Washington, where male stickleback in numerous locations have lost the ancestral red sexual signal and instead develop black nuptial coloration. We measured phenotypic variation in a suite of traits across sites where red and black stickleback do not overlap in distribution and at one site where they historically co‐occurred. We found substantial phenotypic divergence between red and black morphs in noncolor traits including shape and lateral plating, and additionally find evidence that supports the hypothesis of sensory drive as the mechanism responsible for the evolutionary switch in color from red to black. A newly described third “mixed” morph in Connor Creek, Washington, differs in head shape and size from the red and black morphs, and we suggest that their characteristics are most consistent with hybridization between anadromous and freshwater stickleback. These results lay the foundation for future investigation of the underlying genetic basis of this phenotypic divergence as well as the evolutionary processes that may drive, maintain, or limit divergence among morphs.     

Range limits,large‐scale biogeographic variation,and localized evolutionary dynamics in a polymorphic damselfly

Studies of heritable colour polymorphisms allow investigators to track the genetic dynamics of natural populations. By comparing polymorphic populations over large geographic areas and across generations, issues about both morph stability and evolutionary dynamics can be addressed, increasing our understanding of the potential mechanisms maintaining genetic polymorphisms. In the present study, we investigated population morph frequencies in a sex‐limited heritable colour polymorphic damselfly (Ischnura elegans, Vander Linden), with three discrete female morphs. We compared the frequencies of these three female morphs in 120 different populations from ten European countries at differing latitudes and longitudes. There were pronounced differences in morph frequencies both across the entire European biogeographic range, as well as at a smaller scale within regions. We also found considerable between‐population variation at the local scale within regions, particularly at the edges of the range of this species. We discuss these findings in the context of recent models of adaptive population divergence along the range of a species. This polymorphism is thus highly dynamic, with stable morph frequencies at the core of the species range but fluctuating morph dynamics at the range limits. We finish with a discussion of how local interactions and climatic factors can be expected to have a strong influence on the biogeographic patterns in this species and other sexually selected polymorphisms. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 775–785.     

Convergence and divergence in lizard colour polymorphisms

Colour polymorphic species are model systems for examining the evolutionary processes that generate and maintain discrete phenotypic variation in natural populations. Lizards have repeatedly evolved strikingly similar polymorphic sexual signals in distantly related lineages, providing an opportunity to examine convergence and divergence in colour polymorphism, correlated traits and associated evolutionary processes. Herein, we synthesise the extensive literature on lizard colour polymorphisms in both sexes, including recent advances in understanding of the underlying biochemical, cellular and genetic mechanisms, and correlated behavioural, physiological and life-history traits. Male throat, head or ventral colour morphs generally consist of red/orange, yellow and white/blue morphs, and sometimes mixed morphs with combinations of two colours. Despite these convergent phenotypes, there is marked divergence in correlated behavioural, physiological and life-history traits. We discuss the need for coherence in morph classification, particularly in relation to ‘mixed’ morphs. We highlight future research directions such as the genetic basis of convergent phenotypes and the role of environmental variation in the maintenance of polymorphism. Research in this very active field promises to continue to provide novel insights with broad significance to evolutionary biologists.     

Predator-prey relationships and the evolution of colour polymorphism: a comparative analysis in diurnal raptors

Genetically based variation in coloration occurs in populations of many organisms belonging to various taxa, including birds, mammals, frogs, molluscs, insects and plants. Colour polymorphism has evolved in raptors more often than in any other group of birds, suggesting that predator–prey relationships was a driving evolutionary force. Individuals displaying a new invading colour morph may enjoy an initial foraging advantage because prey have difficulties in learning the colour of a rare morph (apostatic selection), or because morphs provide alternative foraging benefits allowing differently coloured individuals to exploit distinct food niches (disruptive selection). Plumage polymorphism should therefore have evolved in species that prey upon animals having the physiological ability to distinguish between differently coloured predators but also to flee once a predator has been detected. From this assumption, we can predict that closely related polymorphic and monomorphic species prey upon different animals. They may also differ in morphology, because foraging upon different prey may require different foraging modes, and in turn different morphological structures. We tested these two predictions in a comparative study of raptors. As expected, polymorphic and monomorphic species had a different diet, and there was a difference in wing length between polymorphic and monomorphic species within two genera ( Buteo and Accipiter ). Across all raptors for which phylogenetic relationships are known, polymorphic species preyed more often upon mammals than did monomorphic ones. These two types of raptor did not differ in the frequency of birds, insects and reptiles in their diets. We discuss these results in the light of the hypothesis that predator–prey relationships played a role in the evolution of colour polymorphism. © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 81 , 565–578.     

Patterns of genetic diversity in the polymorphic ground snake (Sonora semiannulata)

We evaluated the genetic diversity of a snake species with color polymorphism to understand the evolutionary processes that drive genetic structure across a large geographic region. Specifically, we analyzed genetic structure of the highly polymorphic ground snake, Sonora semiannulata, (1) among populations, (2) among color morphs (3) at regional and local spatial scales, using an amplified fragment length polymorphism dataset and multiple population genetic analyses, including F_ST-based and clustering analytical techniques. Based upon these methods, we found that there was moderate to low genetic structure among populations. However, this diversity was not associated with geographic locality at either spatial scale. Similarly, we found no evidence for genetic divergence among color morphs at either spatial scale. These results suggest that despite dramatic color polymorphism, this phenotypic diversity is not a major driver of genetic diversity within or among populations of ground snakes. We suggest that there are two mechanisms that could explain existing genetic diversity in ground snakes: recent range expansion from a genetically diverse founder population and current or recent gene flow among populations. Our findings have further implications for the types of color polymorphism that may generate genetic diversity in snakes.     

Disassortative mating prevails in style‐dimorphic Narcissus papyraceus despite low reciprocity and compatibility of morphs

Evolution to reduce inbreeding can favor disassortative (intermorph) over assortative (intramorph) mating in hermaphroditic sexually polymorphic plant species. Heterostyly enhances disassortative pollination through reciprocal placement of stigmas and anthers of morphs and appropriate pollinators. Stylar dimorphism in which there is not reciprocal anther placement may compromise disassortative mating, particularly when there is not intramorph incompatibility. Variable rates of disassortative mating along with differential female fecundity or siring success among floral morphs could lead to variation in morph ratio. We investigated mating patterns, female fecundity, and siring success of style‐length morphs in Narcissus papyraceus, a self‐incompatible but morph‐compatible species with dimorphic (long‐ and short‐styled) and monomorphic (long‐styled) populations in central and north regions of its range, respectively. We established experimental populations in both regions and exposed them to ambient pollinators. Using paternity analysis, we found similar siring success of morphs and high disassortative mating in most populations. Female fecundity of morphs was similar in all populations. Although these results could not completely explain the loss of dimorphism in the species’ northern range, they provided evidence for the evolutionary stability of stylar dimorphism in N. papyraceus in at least some populations. Our findings support the hypothesis that prevailing intermorph mating is key for the maintenance of stylar dimorphism.     

Population genetics of speciation in nonmodel organisms: I. Ancestral polymorphism in mangroves

The level of DNA polymorphism in the ancestral species at the time of speciation can be estimated using DNA sequences from many loci sampled from 2 or more extant species. The comparison between ancestral and extant polymorphism can be informative about the population genetics of speciation. In this study, we collected and analyzed DNA sequences of approximately 60 genes from 4 species of Sonneratia, a common genus of mangroves on the Indo-Pacific coasts. We found that the 3 ancestral species were comparable to each other in terms of level of polymorphism. However, the ancestral species at the time of speciation were substantially more polymorphic than the extant geographical populations. This ancestral polymorphism is in fact larger than, or at least equal to, the level of polymorphism of the entire species across extant geographical populations. The observations are not fully compatible with speciation by strict allopatry. We suggest that, at the time of speciation, the ancestral species consisted of interconnected but strongly divided geographical populations. This population structure would give rise to high level of polymorphism across species range. This approach of studying the speciation history by genomic means should be applicable to nonmodel organisms.     

Multiple cases of striking genetic similarity between alternate electric fish signal morphs in sympatry

Striking trait polymorphisms are worthy of study in natural populations because they can often shed light on processes of phenotypic divergence and specialization, adaptive evolution, and (in some cases) the early stages of speciation. We examined patterns of genetic variation within and between populations of mormyrid fishes that are morphologically cryptic in sympatry but produce alternate types of electric organ discharge (EOD). Other species in a large group containing a clade of these morphologically cryptic EOD types produce stereotyped, species-typical EOD waveforms thought to function in mate recognition. First, for six populations from Gabon's Brienomyrus species flock, we confirm that forms of electric fish that exhibit distinctive morphologies and unique EOD waveforms (i.e., good reference species) are reproductively isolated from coexisting congeners. These sympatric species deviate from genetic panmixia across five microsatellite loci. Given this result, we examined three focal pairs of syntopic and morphologically cryptic EOD waveform types that are notable exceptions to the pattern of robust genetic partitioning among unique waveform classes within assemblages. These exceptional pairs constitute a monophyletic group within the Brienomyrus flock known as the magnostipes complex. One member of each pair (type I) produces a head-negative EOD, while the other member (either type II or type III, depending on location) produces a longer duration EOD differing in waveform from type I. We show that signal development in these pairs begins with juveniles of all magnostipes-complex morphs emitting head-positive EODs resembling those of type II adults. Divergence of EOD waveforms occurs with growth such that there are two discrete and fixed signal types in morphologically indistinguishable adults at each of several localities. Strong microsatellite partitioning between allopatric samples of any of these morphologically cryptic signal types suggests that geographically isolated populations are genetically decoupled from one another. By contrast, sympatric morphs appear genetically identical across microsatellite loci in Mouvanga Creek and the Okano River and only very weakly diverged, if at all, in the Ivindo River. Our results for the magnostipes complex fail to detect species boundaries between the focal morphs and are, instead, fully consistent with the existence of relatively stable signal dimorphisms at each of several different localities. No mechanism for the maintenance of this electrical polymorphism is suggested by the known natural history of the magnostipes complex. Despite a lack of evidence for genetic differentiation, the possibility of incipient sympatric speciation between morphs (especially type I and type II within the Ivindo River) merits further testing due to behavioral and neurobiological lines of evidence implying a general role for stereotyped EOD waveforms in species recognition. We discuss alternative hypotheses concerning the origins, stability, and evolutionary significance of these intriguing electrical morphs in light of geographical patterns of population structure and signal variation.     

Multivariate phenotypes and the potential for alternative phenotypic optima in wall lizard (Podarcis muralis) ventral colour morphs

A major goal in evolutionary biology is to determine how phenotypic variation arises and is maintained in natural populations. Recent studies examining the morphological, physiological and behavioural differences among discrete colour morphotypes (morphs) have revealed several mechanisms that maintain discrete variation within populations, including frequency‐dependence, density‐dependence and correlational selection. For example, trade‐offs over resource allocation to morphological, physiological and behavioural traits can drive correlational selection for morph‐specific phenotypic optima. Here, we describe a ventral colour polymorphism in the wall lizard (Podarcis muralis) and test the hypothesis that morphs differ along multivariate axes defined by trade‐offs in morphological, physiological, and immunological traits. We show that ventral colour is a discrete trait and that morphs differ in body size, prevalence of infection by parasites and infection intensity. We also find that morphs differ along multivariate phenotypic axes and experience different multivariate selection pressures. Our results suggest that multivariate selection pressures may favour alternative optimal morph‐specific phenotypes in P. muralis.     

Correlated morphological and colour differences among females of the damselfly Ischnura elegans

Abstract 1. The female‐limited colour polymorphic damselfly Ischnura elegans has proven to be an interesting study organism both as an example of female sexual polymorphism, and in the context of the evolution of colour polymorphism, as a model of speciation processes. 2. Previous research suggests the existence of correlations between colour morph and other phenotypic traits, and the different female morphs in I. elegans may be pursuing alternative phenotypically integrated strategies. However, previous research on morphological differences in southern Swedish individuals of this species was only carried out on laboratory‐raised offspring from a single population, leaving open the question of how widespread such differences are. 3. The present study therefore analysed multi‐generational data from 12 populations, investigating morphological differences between the female morphs in the field, differences in the pattern of phenotypic integration between morphs, and quantified selection on morphological traits. 4. It was found that consistent morphological differences indeed existed between the morphs across populations, confirming that the previously observed differences were not simply a laboratory artefact. It was also found, somewhat surprisingly, that despite the existence of sexual dimorphism in body size and shape, patterns of phenotypic integration differed most between the morphs and not between the sexes. Finally, linear selection gradients showed that female morphology affected fecundity differently between the morphs. 5. We discuss the relevance of these results to the male mimicry hypothesis and to the existence of potential ecological differences between the morphs.     

Recent speciation between sympatric Tanganyikan cichlid colour morphs

Lake Tanganyika, Africa's oldest lake, harbours an impressive diversity of cichlid fishes. Although diversification in its radiating groups is thought to have been initially rapid, cichlids from Lake Tanganyika show little evidence for ongoing speciation. In contrast, examples of recent divergence among sympatric colour morphs are well known in haplochromine cichlids from Lakes Malawi and Victoria. Here, we report genetic evidence for recent divergence between two sympatric Tanganyikan cichlid colour morphs. These Petrochromis morphs share mitochondrial haplotypes, yet microsatellite loci reveal that their sympatric populations form distinct genetic groups. Nuclear divergence between the two morphs is equivalent to that which arises geographically within one of the morphs over short distances and is substantially smaller than that among other sympatric species in this genus. These patterns suggest that these morphs diverged only recently, yet that barriers to gene flow exist which prevent extensive admixture despite their sympatric distribution. The morphs studied here provide an unusual example of active diversification in Lake Tanganyika's generally ancient cichlid fauna and enable comparisons of speciation processes between Lake Tanganyika and other African lakes.     

Climate drives community‐wide divergence within species over a limited spatial scale: evidence from an oceanic island

Geographic isolation substantially contributes to species endemism on oceanic islands when speciation involves the colonisation of a new island. However, less is understood about the drivers of speciation within islands. What is lacking is a general understanding of the geographic scale of gene flow limitation within islands, and thus the spatial scale and drivers of geographical speciation within insular contexts. Using a community of beetle species, we show that when dispersal ability and climate tolerance are restricted, microclimatic variation over distances of only a few kilometres can maintain strong geographic isolation extending back several millions of years. Further to this, we demonstrate congruent diversification with gene flow across species, mediated by Quaternary climate oscillations that have facilitated a dynamic of isolation and secondary contact. The unprecedented scale of parallel species responses to a common environmental driver for evolutionary change has profound consequences for understanding past and future species responses to climate variation.     

Colour polymorphism is associated with lower extinction risk in birds

Colour polymorphisms have played a major role in enhancing current understanding of how selection and demography can impact phenotypes. Because different morphs often display alternative strategies and exploit alternative ecological niches, colour polymorphism can be expected to promote adaptability to environmental changes. However, whether and how it could influence populations' and species' response to global changes remains debated. To address this question, we built an up‐to‐date and complete database on avian colour polymorphism based on the examination of available data from all 10,394 extant bird species. We distinguished between true polymorphism (where different genetically determined morphs co‐occur in sympatry within the same population) and geographic variation (parapatric or allopatric colour variation), because these two patterns of variation are expected to have different consequences on populations' persistence. Using the IUCN red list, we then showed that polymorphic bird species are at lesser risk of extinction than nonpolymorphic ones, after controlling for a range of factors such as geographic range size, habitat breadth, life history, and phylogeny. This appears consistent with the idea that high genetic diversity and/or the existence of alternative strategies in polymorphic species promotes the ability to adaptively respond to changing environmental conditions. In contrast, polymorphic species were not less vulnerable than nonpolymorphic ones to specific drivers of extinction such as habitat alteration, direct exploitation, climate change, and invasive species. Thus, our results suggest that colour polymorphism acts as a buffer against environmental changes, although further studies are now needed to understand the underlying mechanisms. Developing accurate quantitative indices of sensitivity to specific threats is likely a key step towards a better understanding of species response to environmental changes.     

Population genetic structure of Arctic charr, Salvelinus alpinus from northwest Europe on large and small spatial scales

To examine the population genetic structure of lake-resident Arctic charr, Salvelinus alpinus from northwest Europe on multiple spatial scales, 2367 individuals from 43 lakes located in three geographical regions (Iceland, the British Isles and Scandinavia) were genotyped at six microsatellite loci. On a large scale, data provided little evidence to support clustering of populations according to geographical region. Hierarchical analysis of molecular variance indicated that, although statistically significant, only 2.17% of the variance in allelic frequencies was partitioned at the among-region level. Within regions, high levels of genetic differentiation were typically found between lakes regardless of the geographical distance separating them. These results are consistent with the hypothesis of rapid postglacial recolonization of all of northwest Europe from a single charr lineage, with subsequent restriction of gene flow. On a smaller scale, there was evidence for close genetic relationships among lakes from within common drainage basins in Scotland. Thus, interlake genetic structure reflects localized patterns of recent (or contemporary) gene flow superimposed onto a larger scale structure that is largely a result of historical processes. There was also evidence for widespread genetic structuring at the within-lake level, with sympatric populations detected in 10 lakes, and multilocus heterozygote deficits found in 23 lakes. This evidence of the Wahlund effect was found in all lakes known to contain discrete phenotypic morphs, as well as many others, suggesting that morphs may often represent separate breeding populations, and also that the phenomenon of polymorphism in this species may be more widespread than is currently realized.     

The dynamics of sexually antagonistic coevolution and the complex influences of mating system and genetic correlation

Sexual conflict has been proposed to be a mediator of speciation but recent theoretical work, as well as empirical studies, suggests that sexual conflict may also be able to prevent speciation and to preserve genetic polymorphism within a species. Here, we develop a population genetic model and study the effects of sexual conflict in a polymorphic population. The morphs mate assortatively based on different sexually antagonistic traits and females are assumed to suffer a cost when the proportion of matching males is high. We consider the model in two different mating systems; promiscuity and polygyny. Our results show that genetic polymorphism may be maintained through negative frequency dependent selection established by assortative mating and female conflict costs. However, the outcome significantly differs between mating systems. Furthermore, we show that indirect selection may have profound effects on the evolutionary dynamics of a sexual conflict.     

Correlated evolution of floral morphology and mating-type frequencies in a sexually polymorphic plant

In sexually polymorphic species, reproductive morphology governs mating patterns and the character of negative frequency-dependent selection. If local environmental conditions cause sexual morphs to differ between populations, then frequency-dependent selection should create corresponding geographic variation in morph frequencies. We investigate this relation with a model of morph-ratio evolution and analysis of geographic variation in the heterostylous plant Narcissus triandrus. Unlike other tristylous species, N. triandrus possesses both imperfect reciprocity among morphs in sex-organ position and a self-incompatibility system that permits outcrossing within and between morphs. We sampled 137 populations throughout the Iberian Peninsula for floral-morph ratios, and measured floral morphology in 31 populations. Morph ratios exhibited three atypical features: (1) predominance of the long-styled (L) morph; (2) absence of the mid-styled (M) morph from 17.5% of populations; and (3) a negative relation between the frequencies of the L and M morphs among populations. Morph ratios varied geographically, with decreasing frequency of the M morph from the southeast to the northwest of the species' range. Much of this variation accompanied allometric change in the positions of sex organs, especially the mid-level organs, with the M morph declining in frequency and ultimately being lost in large-flowered populations. Using multivariate multiple regression, we demonstrate that variation in floral morphology among populations predicts this geographic variation in morph frequencies. Our theoretical analysis illustrates that patterns of pollen transfer governed by imperfect sex-organ reciprocity can select for unequal equilibrium morph ratios like those observed for N. triandrus. We interpret the L-biased morph ratios and the unusual morphology of N. triandrus as a consequence of its atypical intramorph compatibility system.     

Loss of floral polymorphism in heterostylous Luculia pinceana (Rubiaceae): a molecular phylogeographic perspective

Both deterministic and stochastic forces determine the representation and frequency of floral morphs in heterostylous plant populations. Phylogeographic analysis of molecular variation can provide information on the role of historical factors, including founder events, in affecting population morph structure. Here, we investigate geographical patterns of floral morph variation in a distylous shrub Luculia pinceana (Rubiaceae) by examining the relations between floral polymorphism and molecular (cpDNA and microsatellite) variation in 25 populations sampled throughout the distribution of the species in southwest China and adjacent countries. In 19 of the 25 populations, the frequency of floral morphs was not significantly different from the expected 1:1 ratio. The remaining populations were either L‐morph biased (2) or monomorphic (4) for this form and were morphologically differentiated from the remaining populations in several floral traits, that is, corolla tube length, sex organ position and stigma‐anther separation. Phylogeographic analysis supports the hypothesis that L. pinceana was initially split into west‐central and eastern lineages in the Early Pleistocene (～1.982 Mya). A centrally located lineage composed of morph‐biased and monomorphic populations appears to have been subsequently derived from the west‐central lineage, perhaps by a founder event after the last glacial maximum. Hypotheses to explain why these populations have not returned to equilibrium morph frequencies are considered.     

Delimiting shades of gray: phylogeography of the Northern Fulmar,Fulmarus glacialis

The Northern Fulmar (Fulmarus glacialis) is a common tube‐nosed seabird with a disjunct Holarctic range. Taxonomic divisions within the Northern Fulmar have historically been muddled by geographical variation notably including highly polymorphic plumage. Recent molecular analyses (i.e., DNA barcoding) have suggested that genetic divergence between Atlantic and Pacific populations could be on par with those typically observed between species. We employ a multigene phylogenetic analysis to better explore the level of genetic divergence between these populations and to test an old hypothesis on the origin of the modern distribution of color morphs across their range. Additionally, we test whether mutations in the melanocortin‐1 receptor gene (MC1R) are associated with dark plumage in the Northern Fulmar. We confirmed that mitochondrial lineages in the Atlantic and Pacific populations are highly divergent, but nuclear markers revealed incomplete lineage sorting. Genetic divergence between these populations is consistent with that observed between many species of Procellariiformes and we recommend elevating these two forms to separate species. We also find that MC1R variation is not associated with color morph but rather is better explained by geographical divergence.