After 60 years, an answer to the question: what is the Karner blue butterfly?

Matthew L. Forister, Zachariah Gompert, James A. Fordyce, Chris C. Nice

Abstract

The Karner blue butterfly (KBB), Lycaeides melissa samuelis, is a federally protected taxon whose relationship to the Melissa blue, Lycaeides melissa, has been a point of contention during the 66 years since the KBB was first described. Using a large population-genomic dataset and a model of population divergence with migration, we investigated the relationship between the KBB and L. melissa, as well as the relationship between L. melissa and a third taxon, Lycaeides idas. We report that gene flow between the KBB and L. melissa is low, and comparable to gene flow between L. melissa and L. idas. Considering this population-genetic evidence, we conclude that the KBB is a unique evolutionary lineage that should be recognized as Lycaeides samuelis.

1. Introduction

Estimates of gene flow among natural populations have aided in the discovery, demarcation and prioritization of units for conservation, from evolutionarily significant units (ESUs) to species [1]. Typically, efforts in conservation genetics have relied on a small number of genetic markers, and in many cases a single marker, mitochondrial DNA, which has limited the ability of conservationists to address key issues [2]. Technological advances now make it possible to sample an unprecedented number of genetic regions distributed across the genome of non-model species, though these tools have yet to be widely applied to species of conservation concern [3,4]. Here, we focus on the Karner blue butterfly (KBB), Lycaeides melissa samuelis Nabokov, and use a population-genomic dataset to address questions related to the genetic distinctness of this taxon.

Populations of the KBB occur from Minnesota to New Hampshire in association with a single larval host plant, Lupinus perennis, in habitats that include oak savannah and pine barrens [5]. Suitable KBB habitat decreased sharply in the twentieth century as a result of fire suppression and land conversion, resulting in range-wide declines in population size, ultimately leading to placement of the KBB on the US list of endangered species in 1992 [5]. Based on variation in wing pattern and genitalic morphology, the KBB was described as a subspecies of Lycaeides melissa 66 years ago by Vladimir Nabokov in 1944, who subsequently and informally revised his opinion to consider it a distinct species [6,7]. Population-genetic evidence bearing on the distinctness of the KBB was reported by Packer et al. [8], who concluded that the KBB was indistinguishable from L. melissa based on 34 allozyme loci. Subsequent molecular studies raised the possibility, contra Packer et al. [8], that the KBB might be a genetically distinct lineage, albeit with introgression of mitochondrial DNA from L. melissa [9,10].

We used a population-genomic dataset described by Gompert et al. [3] to investigate the relationship between the KBB and the closely related L. melissa. To provide a context for that comparison, we also investigated the relationship between L. melissa and a third closely related taxon, Lycaeides idas. Specifically, we asked the following questions: (i) when did the KBB and L. melissa diverge relative to L. melissa and L. idas? (ii) how do population sizes vary among the three taxa? and (iii) has there been a history of gene flow between the KBB and L. melissa? If so, how do rates of gene flow compare with those between L. melissa and L. idas? Answers to these questions will provide a more complete picture of the evolutionary history of the KBB.

2. Material and methods

The population-genomic dataset used in this study was generated using 454 pyrosequencing of individuals from 12 Lycaeides populations [3]. Here, we focus on four populations: one KBB population, two L. melissa populations (‘west’ and ‘east’) and one L. idas population (figure 1a). The ‘populations’ consist of pooled DNA of 15 individuals from three sampling localities per population. Pyrosequencing of these four populations produced 43 618 sequences, which assembled into a total of 9955 contiguous alignments (contigs), with an average of 329 base pairs per alignment. The KBB population consisted of individuals from Fort McCoy (Wisconsin), Necedah (Wisconsin) and Saratoga (New York). Further locality information, and details of data generation are given in the electronic supplementary material.

Figure 1.

(a) Map of sampled locations and taxa in North America and posterior probability distributions for parameters as follows: (b) splitting times, (c,e and g) populations sizes and (d,f and h) migration rates. Height of curves corresponds to the estimated probability that a given parameter value is true, given the data (95% confidence intervals are reported in table 1).

To investigate the population-genetic history of the KBB and related taxa, we used the isolation with migration model (IM) of Nielsen & Wakeley [11] and Hey & Nielsen [12], as implemented in the software IMa [13]. IMa uses Markov chain Monte Carlo (MCMC) sampling to obtain maximum-likelihood estimates of six parameters: ancestral population size (ΘA), current population sizes (Θ1 and Θ2), rates of migration between two populations subsequent to an initial split (m1 and m2) and the timing of divergence (t). We used separate runs of IMa to investigate the histories of the following pairs of populations: KBB and L. melissa east; L. melissa east and L. melissa west; and L. melissa west and L. idas (figure 1a). Note that the comparisons involving L. melissa are made by geographical proximity, with L. melissa east being closest to the KBB, and L. melissa west closest to L. idas. The number of genetic regions (contigs) involved in the comparisons are as follows: 317 regions for KBB and L. melissa, 256 for the comparison within L. melissa and 188 for L. melissa and L. idas (a total of 440 unique regions across all comparisons). These analyses included contigs with a minimum of three sequences per population; the maximum number per population was 198, and the mean was 9.16 sequences per contig per population. For each pair of populations, 15 independent MCMC simulations were run for 500 000 generations, with a burn in of 100 000 generations; results were combined prior to maximum-likelihood parameter estimation. Additional details of the IM model are discussed in the electronic supplementary material.

The parameters estimated by IMa are in units that incorporate mutation rate (see electronic supplemental material). Given mutation rates (per gene), these parameters can be converted into years or individuals. We do not have an estimate of genome-wide mutation rate for Lycaeides. Instead, we use the parameters directly to draw comparative conclusions about the histories of our taxa. We note also that our taxa have differing numbers of generations per year (KBB has two, L. melissa has two to three and L. idas has one), which would complicate a conversion of the splitting parameter (t) into years. Instead, t can be thought of as the evolutionarily effective amount of time that has passed since divergence.

3. Results

MCMC simulations produced posterior probability distributions for most of our focal parameters that were sharp, suggesting high confidence in our estimates (figure 1; point estimates and 95% confidence intervals are given in table 1). Estimated splitting times for the three comparisons were clustered (figure 1b). Estimates of population sizes found small ancestral populations in all cases, larger contemporary populations for L. melissa and L. idas, and a small contemporary effective population size for the KBB (figure 1c,e and g). Estimated rates of migration were asymmetrical in all cases (figure 1d,f and h). For example, gene flow from L. melissa east to KBB was higher than the reverse direction (figure 1d). The highest rates of gene flow were observed in the intraspecific comparison of the two L. melissa populations: note the point estimates as well as the higher upper bounds of the confidence intervals relative to the comparisons involving the KBB and L. idas (table 1).

View this table:
Table 1.

Estimates of divergence times (t), migration rates (m), Θ (four times the effective population size times μ), and 95% confidence intervals (see supplementary text for more details). Lowercase letters after population names refer to panels in figure 1.

4. Discussion

The ecological and morphological distinctness of the KBB have been well established. Differences in wing pattern and genitalic variation distinguish the KBB from other Lycaeides [6,14]. Ecologically, the KBB uses a single larval host, whereas the closely related L. melissa uses a large number of legumes [15]. Previous genetic work has suggested varying levels of hybridization with L. melissa, or very recent derivation from a common ancestor [3,9]. We find evidence for limited gene flow between L. melissa and the KBB. In particular, we find rates of gene flow between those taxa that are lower than rates of gene flow between populations of L. melissa, and comparable to rates of gene flow between L. melissa and L. idas. It is difficult to make direct comparisons between these results and previous work in this system, given the data involved here which is different in both quantity and quality. However, times of divergence and population size estimates agree with previous work in this system: the split between L. idas and L. melissa, for example, encompasses a greater amount of evolutionary time than the split between L. melissa and the KBB, and the estimated population size for the KBB is smaller than for L. melissa.

Although questions of taxonomy are of secondary interest to the investigation of evolutionary processes, we conclude by asking if subspecific status for the KBB is appropriate given what has been learned in the last 66 years. There are no thresholds for multi-locus genetic distinctness or levels of gene flow that are widely agreed upon as being consistent with taxonomic species status [2,16]. Given the ubiquity and diversity of introgression that has been discovered in recent years for many taxa [17], it is hard to imagine how a universal, taxonomic-genetic threshold would be possible or even desirable. Instead, we suggest that sufficiently powerful genetic data can be used in a comparative sense to address localized taxonomic problems. We report levels of gene flow between the KBB and L. melissa that are comparable to levels of gene flow between L. melissa and L. idas, which are recognized as distinct taxa. Rates of gene flow between the KBB and L. melissa are also lower than rates detected between conspecific populations of L. melissa. In light of the comparative picture of gene flow presented here based on hundreds of loci, coupled with what we know about the distinctiveness of the KBB in terms of host plant association, morphology and life-history variation, we conclude that the current subspecific designation for the KBB is inappropriate, and that the KBB should be elevated to Lycaeides samuelis Nabokov.

  • Received November 15, 2010.
  • Accepted November 26, 2010.

References

View Abstract