Genetic evaluation of the evolutionary distinctness of a federally endangered butterfly, Lange’s Metalmark
© Proshek et al.; licensee BioMed Central. 2015
Received: 8 December 2014
Accepted: 20 April 2015
Published: 25 April 2015
The Mormon Metalmark (Apodemia mormo) species complex occurs as isolated and phenotypically variable colonies in dryland areas across western North America. Lange’s Metalmark, A. m. langei, one of the 17 subspecies taxonomically recognized in the complex, is federally listed under the U.S. Endangered Species Act of 1973. Metalmark taxa have traditionally been described based on phenotypic and ecological characteristics, and it is unknown how well this nomenclature reflects their genetic and evolutionary distinctiveness. Genetic variation in six microsatellite loci and mitochondrial cytochrome oxidase subunit I sequence was used to assess the population structure of the A. mormo species complex across 69 localities, and to evaluate A. m. langei’s qualifications as an Evolutionarily Significant Unit.
We discovered substantial genetic divergence within the species complex, especially across the Continental Divide, with population genetic structure corresponding more closely with geographic proximity and local isolation than with taxonomic divisions originally based on wing color and pattern characters. Lange’s Metalmark was as genetically divergent as several other locally isolated populations in California, and even the unique phenotype that warranted subspecific and conservation status is reminiscent of the morphological variation found in some other populations.
This study is the first genetic treatment of the A. mormo complex across western North America and potentially provides a foundation for reassessing the taxonomy of the group. Furthermore, these results illustrate the utility of molecular markers to aid in demarcation of biological units below the species level. From a conservation point of view, Apodemia mormo langei’s diagnostic taxonomic characteristics may, by themselves, not support its evolutionary significance, which has implications for its formal listing as an Endangered Species.
Global loss of biodiversity, fueled by unprecedented anthropogenic influences, has elevated the importance of conservation biology in mainstream public consciousness. Despite this increased attention, one of the most fundamental challenges of conservation remains unsolved and contentious: how do we accurately identify and delimit manageable units of biodiversity [1-9]? Since species are one of the fundamental units of biology, this challenge is rooted in taxonomy [10-13]. Successful conservation, however, also relies on our understanding of diversity below the species level , as evolutionary and ecological potential is often recognized at the population level [6,15].
Evolutionarily (or Evolutionary) Significant Units (ESUs) were originally developed to facilitate objective prioritization of conservation units (CUs) below the species level . At a time when conservation managers (particularly in mammalian systems) were grappling with sometimes-trivial subspecific classifications , ESUs aimed to refocus conservation resources on populations exhibiting the most distinct evolutionary characteristics. However, as with species concepts in taxonomy (e.g. ), alternate definitions and operational criteria of ESUs have proliferated [1,3-8,18-20], as have alternate qualifiers for CUs (e.g. management units (MUs)/demographic units (DUs) [2,14], discrete population segments (DPSs) , and service-providing units (SPUs) ). Starting with relatively conceptual, integrative origins [1,16] that were criticized for their subjectivity , various redefinitions of the ESU have focused on more objective methods for discriminating populations that have evolutionary potential, such as: 1) the use of consistently congruent gene phylogenies , 2) reciprocal monophyly of mitochondrial DNA (mtDNA) or allele frequency divergence at nuclear loci , and 3) diagnostic characters (including ecological, behavioral, etc.) that exclusively cluster individuals or populations using the phylogenetic species concept . While alleviating some subjectivity, these redefinitions met criticism focused on inconsistent phylogenetic reconstruction , the stringency of reciprocal monophyly and diagnostic characters [6,7], and the utility of the phylogenetic species concept .
In the midst of this operational debate, Crandall et al.  have argued that a dichotomous designation (“ESU or not”) betrays the goal of ESUs and undermines biological complexity. They presented eight categories to discern population distinctiveness based on genetic and ecological exchangeability sensu , and proposed that dichotomous use of the term ESU be abandoned in favor of a more holistic approach . Along similar lines, Fraser & Bernatchez  contended that all ESU concepts share the same “fundamental essence” and goal, but differ in specific optimality criteria. Following the unified species concept , these authors described a framework for adaptive evolutionary conservation that recognizes the situational strengths and weaknesses of each ESU definition and integrates them to define biologically meaningful ESUs anywhere along the evolutionary continuum . In the past 15 years, conceptual redefinitions of the ESU have been replaced by new alternate classifications and subdivisions of CUs. de Guia & Saitoh  proposed the use of full and partial ESUs to distinguish populations described based on knowledge of both neutral and adaptive genetic variation from those described using only one aspect of their variation, respectively. Demographic or management units (DUs or MUs, referred to as DUs from this point on) that describe demographically independent populations (generally within ESUs), have also become commonplace alongside the more commonly discussed ESU [14,15,26].
Balancing the complexity of the evolutionary continuum with the practical inventory needs of conservation management is an Augean task, and it can be instructive to reevaluate old decisions. Here, we reassess the conservation status of an endangered butterfly in California in light of its range-wide population structure and phylogeny. Lange’s Metalmark butterfly, Apodemia mormo langei , is an endangered subspecies characterized by a constrained range and a unique phenotype. Described from the banks of the Sacramento-San Joaquin River system downstream of Sacramento, California, this single site remains the only recognized location for this taxon. Due to habitat loss from sand mining and other activities, and these butterflies’ low vagility and constrained distribution, the population is extremely geographically restricted, and was placed on the US Endangered Species List in 1976 under the Endangered Species Act (ESA) (Federal Register 41:22044, 1976). The Antioch Dunes National Wildlife Refuge (NWR) was established in 1980 to protect A. m. langei as well as two species of wildflowers (the Antioch Dunes evening-primrose, Oenothera deltoides subsp. howellii [Munz] W. Klein and the Contra Costa wallflower, Erysimum capitatum var. angustatum [Greene] G. Rossb.). Lange’s Metalmark has undergone several boom-and-bust cycles, but is being maintained through extensive conservation efforts by several organizations, including a captive rearing program [27-31].
Apodemia mormo langei is a member of a variable species group, the A. mormo species complex. Three species are currently recognized in this complex: A. mormo , A. virgulti , and A. mejicanus [33,34]. Apodemia mormo occurs across western North America, from Mexico to Canada, and is by far the widest ranging metalmark species (Riodinidae) in North America. Apodemia virgulti and A. mejicanus are found in the American Southwest and western Mexico , and A. mejicanus also has an isolated population in Colorado . The species complex shows considerable variation in wing markings, voltinism, flight periods, host use, and oviposition behavior . Although most feed exclusively on plants in the genus Eriogonum (wild buckwheat, Polygonaceae), some also feed on Krameria (Rhatany, Krameriaceae) . Due to this variability, most of which is found in the Southwest region of the USA, there is significant taxonomic interest in the group. Currently 17 described subspecies are recognized in the complex , but the number and status of these taxa is far from settled (e.g. [39-41]).
In this study, we investigate the identification and demarcation of biological units in the A. mormo species complex, from both an evolutionary biology and a conservation biology point of view. We use mitochondrial gene sequence data and microsatellite markers to describe range-wide population structure and phylogenetic relationships in the A. mormo species complex, including the endangered subspecies A. m. langei. To objectively relate our genetic assessment of the group to preexisting taxonomy, we compare a selection of morphological characters based on taxonomic designations to the genetic diversity across the species complex. In light of our new genetic data, we then reassess the conservation status of A. m. langei under various definitions of CUs.
Within the eastern lineage, two major clades are observed: One consisting of the Mexican haplotypes, the one sample from Nevada, and the outlier from San Bernardino, CA (clade E2: Figure 2); and the other containing the remaining eastern haplotypes (clade E1). With the exception of one small clade of haplotypes from Wyoming and Montana with 94% bootstrap support and two distinct monophyletic clades comprised of Montana (localities 51 and 52) and Saskatchewan & Montana haplotypes (localities 43, 44, and 45), all other clades showed shallow relationships with little geographic pattern (see Additional file 1: Figure S1 for complete trees).
The Western lineage was composed of four major clades, labeled W1-W4 (Figure 2, Additional file 1: Figure S1). Clade W1 was composed of specimens from the five northernmost collection locations in California (localities 8–12), including the Antioch Dunes population of A. m. langei. Although branch lengths within this clade were short, divergence of 2.51% was observed relative to the rest of the western lineage. Relationships between the three remaining western clades were unresolved, but each was moderately divergent from the others. Excluding one outlier from Santa Barbara Co., CA (haplotype 1189, locality 31), clade W2 included specimens from four geographically intermediate populations in central California (localities 13–16) as well as one in south central California (locality 31). Clade W3 was comprised of the majority of the individuals from south-central and southern California. Finally, clade W4 included all individuals from the Pacific Northwest, but interestingly, also included several individuals from south-central Californian populations (localities 20, 21, 23, 25). The only western lineage to exhibit significant substructure was Clade W3, which displayed four main internal clades (W3a-W3d). Divergence of these clades ranged from 0.94% to 1.10%, but exhibited little geographic pattern.
Descriptive statistics and tests of isolation by distance: Descriptive statistics, population differentiation values, and isolation by distance (IBD) results for all microsatellite data (Global) and for localities separated regionally (PNW: Pacific Northwest; East: east of the Continental Divide; CA: California)
avg. H obs (95% CI)
IBD r 2
Fu’s F s
Six microsatellite loci were genotyped for 447 specimens. Across all localities, observed heterozygosity ranged from 0.1905 (locality 7) to 0.7578 (locality 10), and strong population differentiation was observed overall (FST = 0.2496: Table 1). In pair-wise tests of differentiation, most population pairs were significantly differentiated (only 13 out of 595 pair-wise comparisons were not significantly differentiated). When considering sets of populations regionally (Pacific Northwest, east of the divide, and California), the Pacific Northwest showed statistically lower heterozygosity (pair-wise tests: PNW vs. CA: W = 176, p = <0.00001; PNW vs. East: W = 0, p = <0.00001; CA vs. East: W = 86, p = 0.14163), while populations east of the Continental divide showed the lowest genetic differentiation. Considering all localities, there was a weak, but significant signature of IBD (Table 1). Regionally, significant signatures were observed in the Pacific Northwest and east of the Continental Divide, but not in California. Observed heterozygosity and pair-wise FST values for all populations are provided in Additional file 1: Table S5.
DAPC and STRUCTURE were also run on regional datasets in an attempt to elucidate population structure at that scale. Regional STRUCTURE analyses did not provide any additional resolution, and all regions mirrored the global analysis (K = 1, 3, and 2 for the Pacific Northwest, California, and east of the Continental Divide; results not shown). DAPC, however, provided resolution for the Pacific Northwest and Californian populations, but not for those east of the Continental Divide (Figure 3b-d). In the Pacific Northwest, the two most geographically distant populations (localities 06 and 07) were separated from a main cluster including the remaining populations (Figure 3c). In California, DAPC clustering roughly matched ancestry estimates from STRUCTURE, although as in the global DAPC analysis, more genetic clusters were estimated (K = 5), which graphically resembled over-splitting of the main groups. Similar to the clustering in the Pacific Northwest, in California, a main cluster was observed that showed some geographic substructure (i.e. geographically proximate populations tend to cluster more closely with each other). Several populations from northern California (localities 08, 09, 10, and 11) were separated from this main cluster, however, including the population of A. m. langei at Antioch Dunes (Figure 3b).
Morphology-based taxonomic assignment and concordance to genetic data
Comparison of genetic clades/clusters and taxonomic determinations: Main mitochondrial DNA clades and microsatellite clusters from STRUCTURE exhibited by specimens assigned to each taxonomic entity that was sampled from the A. mormo species group
Mitochondrial DNA clades
A. mormo mormo
A. mormo cythera
A. mormo langei
A. mormo tuolumensis
A. nr. mormo
A. m. nr. langei
A. virgulti virgulti
A. virgulti nigrescens
A. mejicanus mejicanus
A. mejicanus pueblo
The taxonomy and systematics of the A. mormo complex is unsettled (e.g. [34,39-41]), and we do not assert that these morphological designations are an answer to the delimitation of units within this taxonomic complexity. Rather, we use these designations, based on multiple interpretations and sources, to provide an objective link between our genetic data and the preexisting taxonomic nomenclature. For consistency’s sake, we follow the names used by , despite his warning that his arrangement of taxa within the A. mormo complex “should be considered tentative.” As discussed below, we hope this study can provide insight for a taxonomic revision in the future, which should include a more thorough treatment of the morphology of the group (e.g. including ventral and/or quantitative morphometric characters).
Our study represents the first DNA-based treatment of the Apodemia mormo species complex. Both maternally inherited mtDNA and biparentally inherited microsatellites show deep divergences across the Continental divide, as well as along the west coast. Within these regions, fine-scale population structure varied from being almost absent east of the Continental Divide and in the Pacific Northwest, to being highly structured in California. Mitochondrial DNA and microsatellite markers roughly agreed on population relationships, although several discontinuities were observed when comparing the two marker sets. This data provides a starting point for a reassessment of the taxonomy of the group. Accordingly, we reevaluate support for A. m. langei as an ESU under multiple definitions, and as a federally listed endangered species.
Range-wide population genetic structure/phylogeography
Across its range, the A. mormo complex includes several distinct lineages. The most pronounced of these is the mitochondrial divergence between haplotypes 0366 and h1472, assigned to A. mejicanus pueblo, and the rest of the A. mormo complex, including other specimens assigned to A. mejicanus. This is the largest intraspecific divergence observed in the complex, and these populations merit further attention and focused collecting to elucidate their phylogenetic and taxonomic placement. Apart from these divergent haplotypes, both mitochondrial and microsatellite loci support deep divergences across the Continental Divide and between the Pacific Northwest and California. Although some genetic variation exists east of the Continental Divide (particularly with microsatellites), the nexus of diversity for both marker types is found in California. Here however, the concordance between mtDNA and microsatellites breaks down. Microsatellites show some populations of pure ancestry (inferred via STRUCTURE), that tend to cluster geographically, but also many admixed populations and individuals. Mitochondrial DNA also shows several “pure” clusters of geographically proximate populations, in which all individuals belong to the same clade. Some of these roughly match clusters displayed by the microsatellite dataset, although the geographic positions of the breaks between clusters overlap. California also contains many populations of mixed mitochondrial ancestry, including individuals exhibiting haplotypes from the eastern lineage (localities 37, 61, 64, and 65) and the Pacific Northwest (localities 20, 21, 23, and 25).
Populations in California exhibit higher regional population differentiation (FST) than both other regions, and populations east of the Continental Divide show the lowest, despite having comparable levels of genetic diversity (observed heterozygosity) to California. Despite low vagility in these butterflies, only weak signatures of IBD are observed in the Pacific Northwest and east of the Continental Divide. Virtually no signatures of IBD were observed in California, which is unsurprising given the lack of geographic structure in microsatellite groupings. Populations in the Pacific Northwest were not only distinct, but exhibited statistically lower genetic diversity (heterozygosity) than the other regions, as well as fewer average alleles per locus. This confirms previous findings  but mtDNA haplotypes matching those found in the Pacific Northwest were also found in central California. Given the lower genetic diversity compared to other regions, and the presence of these mtDNA haplotypes, it is likely that the Pacific Northwest region was populated by post-glacial migrants originating in California, who experienced a population bottleneck during this process. Similar patterns have been observed in other species occupying recently glaciated areas (e.g. ), and this pattern is known as the leading edge hypothesis . Tajima’s D and Fu’s F s statistics potentially support recent population bottlenecks in the Pacific Northwest, but these results were not significant. East of the Continental Divide, a model of non-neutral sequence evolution was supported (i.e. population growth or purifying selection ), and interestingly, we also observed eastern mtDNA haplotypes in southern California, Nevada, and Mexico. The presence of these eastern haplotypes west of the Continental Divide may indicate a southwestern origin for the eastern lineage, departing from the previously hypothesized Great Plains origin presented by . If the signature of non-neutral evolution of this mtDNA is due to recent population expansion in this region (which assumes neutrality for mtDNA, but see [46,47]), this would likely make this origin/divergence older than that for the Pacific Northwest. More detailed sampling between California and these two regions will help to resolve these phylogeographic hypotheses, and may reveal finer scale population structure that our sampling omits.
California’s high haplotypic and genetic diversity could be attributed to multiple phenomena. First, introgressive hybridization between neighboring populations can create patterns of admixture and facilitate movement of mtDNA haplotypes between populations/species (e.g. [48,49]), which is possible given the close relationships of members within this complex and their geographic proximity. However, under this scenario we would expect adjacent localities to share more haplotypes than distant ones (i.e. IBD), particularly in organisms that have low vagility, such as A. mormo . Isolation by distance was virtually nonexistent and not statistically significant in California; furthermore, sites that shared identical mtDNA haplotypes were as frequently tens to hundreds of miles apart as they were adjacent. Therefore, a more likely explanation for this diversity may be retained ancestral polymorphism (i.e. incomplete lineage sorting e.g. [44,51,52]). Apodemia mormo’s low vagility and high habitat specificity [35,50] restricts its dispersal ability, which can lead to high levels of phylogeographic subdivision via limited gene flow . This scenario, combined with high genetic diversity and differentiation observed in California, supports the hypothesis of retained ancestral polymorphism.
Comparing genetic data and taxonomy
We found little to no concordance between our taxonomic designations and genetic placement with either mtDNA or microsatellites. The nominate subspecies of the group, A. m. mormo, exemplifies this point at the largest geographic scale considered in this study. From California to British Columbia and from New Mexico to Saskatchewan these butterflies can be more-or-less phenotypically identical, but exhibit deep genetic divergences in both their mitochondrial and nuclear genomes. Regionally, geographic proximity is a better predictor of genetic similarity than phenotypic appearance. However, within the three broad regions considered here, some areas exhibited stronger taxonomic/genetic/geographic correspondence than others. Northern Californian populations, for example, display higher correspondence between genetic variation and geography than in southern California (Figure 2), and even some correspondence between taxonomic assignment and mtDNA subclades (Additional file 1: Figure S1). Conversely, in southern California there is no taxonomic/genetic correspondence, and some mtDNA haplotypes (e.g. h1161, h1163, h1176) are shared between sites and species (Additional file 1: Figure S1).
It is important to reiterate that the taxonomy of the A. mormo complex is far from resolved, and we do not present our morphological taxonomic treatment as a resolution to this situation. Rather, this treatment aimed to provide a relatively objective comparison between our genetic data and the preexisting taxonomy. Our focus on characters associated with A. m. langei, and use of dorsal characters only, narrows the scope of this treatment. Given the morphological complexity within the group, we hope this genetic data can provide a novel lens with which to reevaluate the entire A. mormo complex, although such a reevaluation is beyond the scope of this study.
Apodemia mormo langei
Due to its federally endangered status, the single population of A. m. langei at Antioch Dunes provides a pertinent focus for this genetic assessment of the A. mormo complex. Within northern California (mtDNA clade W1), individuals from Antioch Dunes form a monophyletic mtDNA subclade (Additional file 1: Figure S1), and are most closely related (<0.5% sequence divergence) to individuals from geographically proximate populations (Mt. Diablo and Vallejo, localities 10 and 12, respectively). Genetic distance estimates generated with microsatellites, on the other hand, place Antioch Dunes in closer relationship to Tumey Hills, Arroyo Bayo, and Del Puerto (localities 16, 14, and 13, respectively: Figure 4) to the south. While discordance between maternally- and biparentally-inherited markers is not uncommon [47,54], it does suggest complexity in the evolutionary history of a lineage, which can affect conservation implications. Despite their morphological similarity, A. m. langei is quite genetically divergent from individuals assigned to A. m. nr. langei, collected in the Mendota area (locality 15). Several individuals from this population shared some of A. m. langei’s distinct phenotypic characters, but exhibited divergent mtDNA (clade W2) and high levels of admixture with microsatellite loci. The habitats of these populations are not particularly similar, and the use of different Eriogonum hosts (Additional file 1: Table S1) makes it unlikely that convergent evolution is responsible for the phenotypic similarity.
Placed in the context of the entire species complex, A. m. langei is no more genetically distinctive than most populations in California, and other populations exist in this region that exhibit higher mitochondrial and nuclear divergence. Additionally, some of the morphological characteristics that earned A. m. langei its subspecific status are found in other populations in California. In this light, does this current data support ESU status for A. m. langei and, if so, under what ESU criteria? Along with many populations in northern California, the population at Antioch Dunes is reciprocally monophyletic with respect to mtDNA (Additional file 1: Figure S1), satisfying half of Moritz’s  criteria. This mitochondrial pattern may also support criteria such as lack of “genetic exchangeability” (no recent gene flow ), “long term isolation” , and “highly-restricted gene flow from other such lineages within the higher organization level of the species” , although the lack of quantitative determination of these criteria highlights their subjective nature. Microsatellite loci support A. m. langei’s genetic distinctiveness, but only with regional, ordination-based methods (DAPC), not with individual-based Bayesian clustering methods (STRUCTURE). Ecologically, A. m. langei is potentially distinctive, as its host Eriogonum nudum var. psychicola Reveal 2007 is also an endemic of the Antioch Dunes . Group-wide information on A. mormo hosts, however, is relatively sparse (Additional file 1: Table S1), so we are reluctant to draw broad conclusions about the variability of host plants as criteria for “ecological exchangeability”  or general ecological distinctiveness. Below the level of an ESU, A. m. langei’s low vagility combined with its moderate genetic distinctiveness likely demonstrates demographic independence, a criterion of DUs [2,14]. In fact, our sole use of potentially neutral genetic loci (rather than both neutral and adaptive loci) may be more appropriate for defining DUs , and some argue that without both types of loci, our conclusion is inherently limited to defining a “partial ESU” rather than a “full ESU” . Regardless, defining DUs quantitatively (e.g. ) would require increased sampling and finer-scale genetic data.
Collectively, these data provide conflicting support for A. m. langei’s evolutionary significance. While more data (particularly that of host plant) will provide insight into this issue, we expect that the multifarious nature of the A. mormo complex will undermine future determinations that A. m. langei is more evolutionarily divergent than other isolated populations in northern California. However, this study’s focus on evolutionary significance is not an effort to describe or compare all members of the A. mormo complex evolutionarily, but to address whether genetic data support A. m. langei’s federal listing. As we have illustrated, this question can be answered differently depending on the ESU criteria considered. While one side of this answer is that genetic data do not support A. m. langei’s federal listing, we believe that such an interpretation oversimplifies the role of endangered species with regard to ecosystems, particularly in the case of the Antioch Dunes. Balancing the evolutionary significance of a single species against conservation efforts for an entire ecosystem is an exceptionally difficult task, particularly when valuation of those species to the “success” of conservation is not tied to obvious economic value [56,57]. Apodemia mormo langei is undeniably a “flagship species” [58,59] for the Antioch Dunes ecosystem (e.g. [60,61]), but how can a value be placed on public perception or engagement? Additionally, should the momentum of several decades of conservation effort towards this ecosystem be included in a valuation? These conceptual questions are much farther-reaching than the scope of this study, however for these reasons we are not arguing for the delisting of A. m. langei. We hope our genetic assessment of the complex provides a foundation to reevaluate the taxonomy of the group as well as current and future conservation efforts, and we stress that these reassessments should occur in that order. Without a thorough understanding of the group’s taxonomy and systematics, conservation prioritization may not be maximally efficient. This is particularly true for invertebrates, as federally recognized units of conservation below the subspecies level are generally only applied to vertebrates (“distinct population segments” ).
Here we used several genetic markers to assess range-wide population genetic structure of the A. mormo species complex. We then used this information to reassess the evolutionary significance of a federally endangered subspecies within the group. We found highly divergent lineages across the range of the A. mormo complex, indicating a complex evolutionary history. The nexus of genetic diversity was observed in California, where delimitation of highly structured populations agreed poorly with taxonomic designations. Apodemia mormo langei is no more genetically unique than various other populations of A. mormo in California, and even some of the morphological characteristics that earned it subspecific status are not unique. These results indicate that both the taxonomy and conservation prioritization of the A. mormo complex should be reassessed, both with a more fine-grained genetic survey and a greater focus on locally adapted phenotypic traits.
A total of 548 specimens of Apodemia and outgroups were obtained from six principal sources (Additional file 1: Table S8). Specimens of A. mormo langei from Antioch Dunes NWR were collected in 1997 under US Fish & Wildlife Service collection permit PRT-832200. Collection, vouchering, and preservation of specimens differed among sources (Additional file 1: Table S1). We sequenced the full mitochondrial cytochrome oxidase subunit I (COI) gene (1498 base pairs) and up to five microsatellite loci were scored from as many specimens as possible. The principal exceptions were 82 specimens from Opler, Davenport et al. (Source #3, Additional file 1: Table S8). The only genetic data available to us from those samples were the 648 base pairs of the “barcode region” of the COI gene, which were obtained by PO from the Barcode of Life Database (BOLD). In total, sampling of the A. mormo species complex spanned 69 geographic locations in 12 states and two provinces in three countries (Additional file 1: Table S1), comprising most of the known range.
All vouchered specimens were photographed, with the exception of 38 specimens from Grasslands National Park, SK, which were preserved in ethanol. Dorsal-view photographs of most specimens were taken by BP with an 8.0 megapixel Nikon Coolpix 8400 mounted on an Olympus SZX16 dissecting microscope illuminated with a fiber-optic light source. Images are available at the DNA voucher site of the Strickland Museum of Entomology (http://www.biology.ualberta.ca/uasm/Vouchers/index.html). Dorsal images of the specimens from Opler, Davenport et al. (Source #3, Additional file 1: Table S8) were obtained by PO from BOLD (www.boldsystems.org).
Several methods of DNA extraction were used. DNA was extracted from the samples collected by Proshek et al., Powell, and Davenport (sources #1, 4, 5, Additional file 1: Table S8) from two legs (or leg fragments and antennae, if the specimen was in poor condition) using the DNeasy Tissue Extraction Kit (Qiagen, Valencia, CA). DNA from the specimens collected by Sperling, Powell et al. (Additional file 1: Table S8) was extracted from the thorax using a phenol-chloroform method as outlined in . Sequences for the samples collected by Opler, Davenport et al. (Additional file 1: Table S8) were obtained from the Canadian Centre for DNA Barcoding (Guelph, ON) (www.dnabarcoding.ca) , but we did not have access to these DNA extractions. DNA was extracted from the samples collected by Crawford and Desjardins (Additional file 1: Table S8) from wing clips as in .
The mitochondrial gene COI was sequenced in its entirety for as many specimens as possible. In total, 469 sequences of 1498 base pairs in length were obtained. For the samples collected by Sperling, Powell et al. (Additional file 1: Table S8), 398 base pairs of the gene were initially sequenced using the primers Jerry (C1-J-2183)  and K741 (C1-N-2578a) , following . The rest of the gene was later amplified in two fragments using the primer pairs LCO1490 to HCO2198  and BrianXXVII to Pat . For all other samples, the COI gene was sequenced in two fragments: LCO to HCO and Jerry to Pat, unless chromatogram signal was poor, in which case the internal primers Jerry and Mila (MilaX, GATAGTCCTGTAAATAATGG, for samples from west of the Rocky Mountains and MilaXI, GATAATCCTGTAAATAATGG, for samples from east of the Rocky Mountains) and BrianXXVII and Pat were used. The polymerase chain reaction and cycle sequencing protocols are given in detail in . Chromatograms were checked for signal quality in Lasergene (DNASTAR, Madison, WI). Priming sites were manually removed and sequences were manually aligned in Mesquite 2.72 .
Microsatellite development, amplification and genotyping
We isolated and characterized six novel microsatellite loci from two libraries. Details of library development and locus amplification are given in . Genotyping was carried out in GeneMapper (Applied Biosystems, Foster City, CA). We obtained genotype scores for a total of 447 samples from all sampling sources except Source #3 (Additional file 1: Table S8). Amplification success was not consistent across sampling areas. Locus E7 did not amplify in individuals west of the Rockies, and locus M2 did not amplify for individuals east of the Rockies (one exception: the two samples from Sonora, Mexico did not amplify at E7 but one did amplify at M2). All samples, therefore, were genotyped at a maximum of five loci. Sixty samples were genotyped at four loci and 52 at three loci. Samples that amplified less than three loci were not used.
There was a total of 205 unique COI haplotypes: 157 haplotypes were 1498 base pairs in length, and 48 “barcode” haplotypes 648 base pairs in length. Haplotypes were only considered unique if there was at least one base substitution relative to all other haplotypes. Missing base pairs were scored as “N” (missing). All 1498-base pair COI haplotypes from the specimens from Sperling, Powell et al. (Additional file 1: Table S8) had 18 missing base pairs from the middle of the haplotypes where the internal primers overlapped; all other 1498-base pair COI haplotypes had 11 missing base pairs at the same location.
A maximum-likelihood phylogenetic analysis of the COI sequences was performed in Garli 1.0  under the TPM2uf + I + G model, which was selected by jModelTest 0.1.1 [71,72] as the most likely model for our data under the AIC, AICc, and BIC model selection criteria. For the best-tree analysis, rates were constrained so that r[AC] = r[AT], r[AG] = r[CT], and r[CG] = r[GT]. The rate parameters, base frequencies, proportion of invariable sites, and gamma shape parameter were estimated during analysis. Twenty-five search replicates were performed to find the best tree. Two hundred fifty bootstrap replicates were also performed under the same model, except with parameters fixed at the following values: r[AC] = r[AT] = 4.6640; r[AG] = r [CT] = 36.8988; r[CG] = r [GT] = 1.000; eqA = 0.3248, eqC = 0.1293, eqG = 0.1129, eqT = 0.4330; proportion invariable sites = 0.4780; and gamma shape parameter = 0.3120. Calephelis wrighti and Emesis emesia were selected as outgroups, both of which are members of the subfamily Riodinidae; the former in the Riodinini and the latter, like Apodemia, is incertae sedis . We tested for neutral sequence evolution by calculating Tajima’s D  and Fu’s F s statistics  in Arlequin v3.5 . Positive values for these statistics indicate potential population bottlenecks or balancing/overdominant selection, while negative values indicate potential population size increase or purifying selection . These statistics were calculated globally, regionally (Pacific Northwest, east of the divide, and California), and for each population in California (Tajima’s D and Fu’s F s were calculated for the other two regions in ).
The program STRUCTURE  was used to determine the smallest number of genetic clusters that maximized Hardy-Weinberg equilibrium. We tested K values (number of genetic clusters) between 2 and 20 with seven replications each, using the admixture model and correlated allele frequencies. 40,000 burn-in generations and 240,000 post burn-in generations were run. ∆K was calculated following the Evanno et al.  method using STRUCTURE HARVESTER , and CLUMPP was used to average multiple runs of each K value . STRUCTURE was run with and without population information as a prior, but results did not differ significantly with the addition of population information. To assess microsatellite clustering without assuming Hardy-Weinberg and gametic equilibrium, we conducted discriminant analysis of principal components (DAPC ). This method maximizes between- and minimizes within-group variability by conducting a principal components analysis on genetic data, before submitting those principal components to a discriminant analysis. DAPC was implemented in R v3.0.1  using adegenet v1.3.1 . adegenet’s find.clusters function was used to estimate the ideal value of K (default parameters, retaining all principal components), and optim.a.score was used to estimate the ideal number of principal components to retain in the final discriminant analysis (using default settings and 25 full simulations).
Pair-wise population differentiation (FST) and summary statistics (heterozygosity, number of alleles per locus) were calculated using GENEPOP v4.2 [85,86] and the Excel Microsatellite Toolkit v3.1.1 , respectively. To test for signatures of isolation by distance (IBD ), matrices of population-level, standardized genetic differentiation [FST/(1-FST)]  and geographic distance between localities were constructed using GENEPOP v4.2 and the Geographic Distance Matrix Generator v1.2.3 , respectively. The degree of correlation between these matrices was evaluated with a Mantel test  implemented in GENEPOP v4.2 using 9,999 randomizations. All measures of population differentiation, descriptive statistics, and tests for IBD were evaluated for the entire microsatellite dataset, as well as for each main geographic region (Pacific Northwest, east of the divide, and California) separately. To compare regional values of heterozygosity, non-parametric Wilcox-rank sum tests were conducted in R . A Bonferroni correction was used for all multiple pair-wise tests.
DNA sequence & microsatellite comparison
In order to provide a more direct comparison between microsatellites and mtDNA, a phylogenetic tree was generated using 154 unique mtDNA haplotypes from samples for which we also had microsatellite genotypes. We generated a maximum-likelihood tree in Garli 1.0  under the GTR + I + G model selected by jModelTest 0.1.1 [71,72] as the most likely for our data under both the AIC and hLRT criteria. All parameters were estimated during analysis. Twenty-five search replicates were performed to find the best tree. Two hundred fifty bootstrap replicates were also performed under the same conditions. All trees were unrooted. The terminal tips were manually condensed into simplified groupings that approximated the populations from which the haplotypes were sampled, in order to compare the topology of the tree to a tree generated by analysis of microsatellite genetic distances in a priori populations. The program Poptree2  was used to generate a neighbor-joining tree based on the Da genetic distances  of population microsatellite allele frequencies within a priori populations. 1000 bootstrap replicates were performed.
Eleven wing characters were selected to objectively relate the genetic data produced in this study to the existing taxonomic nomenclature of the 17 subspecies in the A. mormo species complex recognized by  (Additional file 1: Table S2). All subspecies in the A. mormo complex (as documented by ) were included in this morphological treatment, despite our genetic focus on only part of the complex. Characters were selected based on the original taxonomic descriptions (Additional file 1: Table S3), further descriptions of character variation in , examination of photographs in the Butterflies of America website , and geographic considerations. Images of type specimens were examined when available; otherwise images of several representative specimens were used. Characters were limited to the dorsal side only and chosen to be independent of specimen size or interpretation of shades of color. This allowed specimens to be scored with these characters based only on a dorsal-view photograph, irrespective of light source or camera settings, and without need of a scale bar.
Availability of supporting data
The sequence dataset supporting the results of this article are available in the GenBank repository (Additional file 1: Table S9), and other data files are available upon request.
We thank Jim Brock, Evi Buckner, Chuck and Cindy Harp, Dan Rubinoff, Eric Runquist, and Ray Stanford for additional collections, and the ‘All Leps’ campaign of the ‘Barcode of Life’ initiative of the University of Guelph, Ontario, for sequencing the ‘barcode’ fragment of mtDNA in these specimens. We thank Jennifer Frei, Betsy Radke, John Steiner and David Wright for permits to collect A. mormo langei; and Anthony Cognato, Becky Miller, and Robert Reed for field assistance at Antioch. We also thank all the people who contributed to the northern collections that are documented in detail in Proshek et al. , and Corey Davis for his advice on microsatellite genotyping of the additional material in the current study. Lab and field activities were funded by an NSERC Discovery Grant and a Hatch grant to Felix Sperling. Finally, we thank Boyd Mori for analytical assistance, and Bryan Brunet, Jason Dombroskie, Lisa Lumley, and Thomas Simonsen for advice and assistance in preparing this manuscript.
- Waples RS. Pacific salmon, Oncorhyncus spp., and the definition of ‘species’ under the Endangered Species Act. Mar Fish Rev. 1991;53:11–21.Google Scholar
- Moritz C. Defining ‘Evolutionarily Significant Units’ for conservation. Trends Ecol Evol. 1994;9:373–5.View ArticlePubMedGoogle Scholar
- Moritz C. Applications of mitochondrial DNA analysis in conservation: a critical review. Mol Ecol. 1994;3:401–11.View ArticleGoogle Scholar
- Vogler AP, DeSalle R. Diagnosing units of conservation management. Conserv Biol. 1994;8:354–63.View ArticleGoogle Scholar
- Bowen BW. What is wrong with ESUs? the gap between evolutionary theory and conservation principles. J Shellfish Res. 1998;17:1355–8.Google Scholar
- Crandall KA, Bininda-Emonds ORP, Mace G, Wayne RK. Considering evolutionary processes in conservation biology. Trends Ecol Evol. 2000;15:290–5.View ArticlePubMedGoogle Scholar
- Fraser DJ, Bernatchez L. Adaptive evolutionary conservation: towards a unified concept for defining conservation units. Mol Ecol. 2001;10:2741–52.View ArticlePubMedGoogle Scholar
- de Guia APO, Saitoh T. The gap between the concept and definitions in the evolutionarily significant unit: the need to integrate neutral genetic variation and adaptive variation. Ecol Res. 2007;22:604–12.View ArticleGoogle Scholar
- Forister ML, Nice CC, Fordyce JA, Gompert Z, Shapiro AM. Considering evolutionary processes in the use of single-locus genetic data for conservation, with examples from Lepidoptera. J Insect Conserv. 2008;12:37–51.View ArticleGoogle Scholar
- Hey J, Waples RS, Arnold ML, Butlin RK, Harrison RG. Understanding and confronting species uncertainty in biology and conservation. Trends Ecol Evol. 2003;18:597–603.View ArticleGoogle Scholar
- Mace GM. The role of taxonomy in species conservation. Philos Trans R Soc Lond B Biol Sci. 2004;359:711–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Frankham R, Ballou JD, Dudash MR, Eldridge MDB, Fenster CB, Lacy RC, et al. Implications of different species concepts for conserving biodiversity. Biol Conserv. 2012;153:25–31.View ArticleGoogle Scholar
- Sackett LC, Seglund A, Guralnick RP, Mazella MN, Wagner DM, Busch JD, et al. Evidence for two subspecies of Gunnison’s prairie dogs (Cynomys gunnisoni), and the general importance of the subspecies concept. Biol Conserv. 2014;174:1–11.View ArticleGoogle Scholar
- Luck GW, Daily GC, Ehrlich PR. Population diversity and ecosystem services. Trends Ecol Evol. 2003;18:331–6.View ArticleGoogle Scholar
- Funk WC, McKay JK, Hohenlohe PA, Allendorf FW. Harnessing genomics for delineating conservation units. Trends Ecol Evol. 2012;27:489–96.View ArticlePubMed CentralPubMedGoogle Scholar
- Ryder OA. Species conservation and systematics: the dilemma of subspecies. Trends Ecol Evol. 1986;1:9–10.View ArticleGoogle Scholar
- de Queiroz K. Species concepts and species delimitation. Syst Biol. 2007;56:879–86.View ArticlePubMedGoogle Scholar
- Dizon AE, Lockyer C, Perrin WF, Demaster DP, Sisson J. Rethinking the stock concept—a phylogeographic approach. Conserv Biol. 1992;6:24–36.View ArticleGoogle Scholar
- Avise JC. Molecular markers, natural history and evolution. New York: Chapman and Hall; 1994.View ArticleGoogle Scholar
- Casacci LP, Barbero F, Balletto E. The “evolutionarily significant unit” concept and its applicability in biological conservation. Ital J Zool. 2013;81:182–93.View ArticleGoogle Scholar
- (USFWS and NMFS) US Fish and Wildlife Service and National Marine Fisheries Service. Policy regarding the recognition of distinct vertebrate population segments under the Endangered Species Act. Fed Regist. 1996;61:4721–5.Google Scholar
- Moritz C, Lavery S, Slade R. Using allele frequency and phylogeny to define units for conservation and management. In: Nielsen JL, Powers GA, editors. Evolution and the Aquatic Ecosystem: Defining Unique Units in Population Conservation. Maryland: American Fisheries Society; 1995. p. 249–62.Google Scholar
- Waples RS. Evolutionary significant units and the conservation of biological diversity under the Endangered Species Act. In: Nielsen JL, Powers GA, editors. Evolution and the Aquatic Ecosystem: Defining Unique Units in Population Conservation. Maryland: American Fisheries Society; 1995. p. 8–27.Google Scholar
- Templeton AR. The meaning of species and speciation: a genetic perspective. In: Otte D, Endler JA, editors. Speciation and its Consequences. Sunderland: Sinauer; 1989. p. 3–27.Google Scholar
- de Queiroz. The general lineage concept of species, species criteria, and the process of speciation: A conceptual unification and terminological recommendations. In: Howard DJ, Berlocher SH, editors. Endless Forms: Species and Speciation. New York: Oxford University Press; 1998. p. 57–75.Google Scholar
- Palsbøll PJ, Bérubé M, Allendorf FW. Identification of management units using population genetic data. Trends Ecol Evol. 2007;22:11–6.View ArticlePubMedGoogle Scholar
- (USFWS) US Fish and Wildlife Service. Revised Recovery Plan for Three Endangered Species Endemic to the Antioch Dunes, California (Lange’s Metalmark Butterfly, Contra Costa Wallflower, and Antioch Dunes Evening Primrose). Portland, OR: U.S. Fish and Wildlife Service; 1984.Google Scholar
- (USFWS) US Fish and Wildlife Service. Antioch Dunes National Wildlife Refuge Comprehensive Conservation Plan. Sacramento, CA: US Fish and Wildlife Service, California/Nevada Planning Office; 2002.Google Scholar
- (USFWS) US Fish and Wildlife Service: Species account: Lange’s Metalmark butterfly, Apodemia mormo langei. [http://ecos.fws.gov/docs/five_year_review/doc1927.pdf]; 2008.
- Johnson JJ, Longcore T, Clause A, Pratt G, Dunn J, Osborne K. Propagation handbook for Lange’s metalmark butterfly, Apodemia mormo langei. 1st ed. The Urban Wildlands Group: Los Angeles; 2007.Google Scholar
- Powell JA, Parker MW. Lange’s Metalmark, Apodemia mormo langei Comstock. In: Conservation Biology of Lycaenidae. New TR: 119; 2010. p. 116.Google Scholar
- Felder C, Felder R. Lepidopterologische Fragmente. Wien Entomol Monatschr. 1859;3:271.Google Scholar
- Behr H. On California Lepidoptera. Proc Cal Acad Sci. 1865;3:178–9.Google Scholar
- Pelham JP. A catalogue of the butterflies of the United States and Canada with a complete bibliography of the descriptive and systematic literature. J Res Lepid. 2008;40:1–652.Google Scholar
- Opler PA, Powell JA. Taxonomic and distributional studies on the Western components of the Apodemia mormo complex (Riodinidae). J Lep Soc. 1961;15:145–71.Google Scholar
- Scott JA. The butterflies of North America: A natural history and field guide. Stanford: Stanford University Press; 1986.Google Scholar
- Wick AA, Janelle J, Pruss S, Erbilgin N. First observations of Mormon Metalmark (Apodemia mormo) oviposition behaviour in Canada. Can Field Nat. 2012;126:34–7.Google Scholar
- Opler PA, Wright AB. A Field Guide to Western Butterflies. Boston: Houghton Mifflin Company; 1999.Google Scholar
- Pratt G, Ballmer G. Three biotypes of Apodemia mormo (Riodinidae) in the Mojave Desert. J Lep Soc. 1991;45:46–57.Google Scholar
- Davenport K. A concise update of the information provided in The Butterflies of Southern California (1973) by Thomas C Emmel and John F Emmel. The Taxonomic Report of the International Lepidoptera Survey. 2004;4:1–23.Google Scholar
- Pratt G, Emmel JF, Bernard G. The buckwheat metalmarks, American Butterflies. 2011. p. 4–31.Google Scholar
- Proshek B, Crawford LA, Davis CS, Desjardins S, Henderson AE, Sperling F. Apodemia mormo in Canada: population genetic data support prior conservation ranking. J Insect Conserv. 2013;17:155–70.View ArticleGoogle Scholar
- Comstock JA. A new Apodemia from California (Lepidopt). Bull South Calif Acad Sci. 1939;37:129–31.Google Scholar
- Patricelli D, Sielezniew M, Poniwicka-Tyszko D, Ratkiewicz M, Bonelli S, Barbero F, et al. Contrasting genetic structure of rear edge and continuous range populations of a parasitic butterfly infected with Wolbachia. BMC Evol Biol. 2013;13:14.View ArticlePubMed CentralPubMedGoogle Scholar
- Hewitt G. The genetic legacy of the Quaternary ice ages. Nature. 2000;405:907–13.View ArticlePubMedGoogle Scholar
- Ballard JWO, Whitlock MC. The incomplete natural history of the mitochondria. Mol Ecol. 2004;13:729–44.View ArticlePubMedGoogle Scholar
- Toews DPL, Brelsford A. The biogeography of mitochondrial and nuclear discordance in animals. Mol Ecol. 2012;21:3907–30.View ArticlePubMedGoogle Scholar
- McGuire JA, Linkem CW, Koo MS, Hutchison DW, Lappin AK, Orange DI, et al. Mitochondrial introgression and incomplete lineage sorting through space and time: phylogenetics of crotaphytid lizards. Evolution. 2007;61:2879–97.View ArticlePubMedGoogle Scholar
- Schmidt BC, Sperling FAH. Widespread decoupling of mtDNA variation and species integrity in Grammia tiger moths (Lepidoptera: Noctuidae). Syst Ent. 2008;33:613–34.View ArticleGoogle Scholar
- Arnold RA, Powell JA. Apodemia mormo langei: Ecological Studies of Six Endangered Butterflies (Lepidoptera, Lycaenidae): Island Biogeography, Patch Dynamics, and the Design of Habitat Preserves. University of California Publications in Entomology. 1983;99:99–128.Google Scholar
- Bulgin NL, Gibbs L, Vickery P, Baker AJ. Ancestral polymorphisms in genetic markers obscure detection of evolutionary distinct population in the endangered Florida grasshopper sparrow (Ammodramus savannarum floridanus). Mol Ecol. 2003;12:831–44.View ArticlePubMedGoogle Scholar
- Steiner CC, Mitelberg A, Tursi R, Ryder OA. Molecular phylogeny of extant equids and effects of ancestral polymorphism in resolving species-level phylogenies. Mol Phylogenet Evol. 2012;65:573–81.View ArticlePubMedGoogle Scholar
- Papadopoulou A, Anastasiou I, Keskin B, Vogler AP. Comparative phylogeography of tenebrionid beetles in the Aegean archipeligo: the effect of dispersal ability and habitat preference. Mol Ecol. 2009;18:2503–17.View ArticlePubMedGoogle Scholar
- Hammouti N, Schmitt T, Seitz A, Kosuch J, Veith M. Combining mitochondrial and nuclear evidences: a refined evolutionary history of Erebia medusa (Lepidoptera: Nymphalidae: Satyrinae) in Central Europe based on the COI gene. J Zool Sys Evol Res. 2009;48:115–25.View ArticleGoogle Scholar
- Reveal JL. A new variety of Eriogonum nudum (Polygonaceae) from California. Phytologia. 2007;89:287–9.Google Scholar
- Heal G. Valuing ecosystem services. Ecosystems. 2000;3:24–30.View ArticleGoogle Scholar
- Mace GM. Whose conservation? Science. 2014;345:1558–60.View ArticlePubMedGoogle Scholar
- Caro TM, O’Doherty G. On the use of surrogate species in conservation biology. Conserv Biol. 1999;13:805–14.View ArticleGoogle Scholar
- Walpole MJ, Leader-Williams N. Tourism and flagship species in conservation. Biodivers Conserv. 2002;11:543–7.View ArticleGoogle Scholar
- Bandy C: Antioch Dunes National Wildlife Refuge. Tidelines Newsletter 2002, 22 [http://www.fws.gov/refuge/san_pablo_bay/TidelineArchives.html].
- Nickles J: Fluttering to extinction? preventing the Lange’s metalmark butterfly from floating away. Tidelines Newsletter 2007, 27 [http://www.fws.gov/refuge/san_pablo_bay/TidelineArchives.html].
- Waples RS, Nammack M, Cochrane JF, Hutchings JA. A tale of two acts: endangered species listing practices in Canada and the United States. Bioscience. 2013;63:723–34.View ArticleGoogle Scholar
- Sperling FAH, Harrison RG. Mitochondrial DNA variation within and between species of the Papilio machaon group of swallowtail butterflies. Evolution. 1994;48:408–22.View ArticleGoogle Scholar
- Hajibabaei M, de Waard JR, Ivanova NV, Ratnasingham S, Dooh RT, Kirk SL, et al. Critical factors for assembling a high volume of DNA barcodes. Philos Trans R Soc Biol B Biol Sci. 2005;360:1959–67.View ArticleGoogle Scholar
- Keyghobadi N, Crawford L, Maxwell S. Successful analysis of AFLPs from non-lethally sampled wing tissues in butterflies. Conserv Genet. 2009;10:2021–4.View ArticleGoogle Scholar
- Simon C, Frati F, Beckenbach A, Crespi BJ, Liu H, Flook P. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am. 1994;87:651–701.View ArticleGoogle Scholar
- Caterino MS, Sperling FAH. Papilio phylogeny based on mitochondrial cytochrome oxidase I and II genes. Mol Phylogenet Evol. 1999;11:122–37.View ArticlePubMedGoogle Scholar
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol and Biotech. 1994;3:294–9.Google Scholar
- Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. [http://mesquiteproject.org]; 2009.
- Zwickl DJ: Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. Ph.D. thesis, The University of Texas at Austin; 2006.Google Scholar
- Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704.View ArticlePubMedGoogle Scholar
- Brower, AVZ: Riodinidae Grote 1895. Metalmarks. The Tree of Life Web Project. [http://tolweb.org/Riodinidae/12174/2008.01.01]; 2008.
- Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.PubMed CentralPubMedGoogle Scholar
- Fu YX. New statistical teszts of neutrality for DNA samples from a population. Genetics. 1996;143:557–70.PubMed CentralPubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005;1:47–50.PubMed CentralGoogle Scholar
- Ramos-Onsins SE, Rozas J. Statistical properties of new neutrality tests against population growth. Mol Biol Evol. 2002;19:2092–100.View ArticlePubMedGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.PubMed CentralPubMedGoogle Scholar
- Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.View ArticlePubMedGoogle Scholar
- Earl DA, Von Holdt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4:359–61.View ArticleGoogle Scholar
- Jakobson M, Rosenberg NA. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007;23:1801–6.View ArticleGoogle Scholar
- Jombart T, Devillard S, Balloux F. Discriminate analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010;11:94.View ArticlePubMed CentralPubMedGoogle Scholar
- R Core Team: R: a language and environment for statistical computing. [http://www.r-project.org/]; 2013.
- Jombart T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24:1403–5.View ArticlePubMedGoogle Scholar
- Raymond M, Rousset F. GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995;86:248–9.Google Scholar
- Rousset F. Genepop’007: a complete reimplementation of the Genepop software for Windows and Linux. Mol Ecol Resourc. 2008;8:103–6.View ArticleGoogle Scholar
- Park SDE: Trypanotolerance in West African Cattle and the Population Genetic Effects of Selection. Ph.D. thesis. University of Dublin; 2001.Google Scholar
- Wright S. Isolation by distance. Genetics. 1943;28:114–38.PubMed CentralPubMedGoogle Scholar
- Rousset F. Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics. 1997;145:1219–28.PubMed CentralPubMedGoogle Scholar
- Erst PJ: Geographic distance matrix generator (version 1.2.3). [http://biodiversityinformatics.amnh.org/open_source/gdmg/]; 2014.
- Mantel N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967;27:207–20.Google Scholar
- Takezaki N, Nei M, Tamura K. POPTREE2: Software for constructing population trees from allele frequency data and computing other population statistics with Windows interface. Mol Biol Evol. 2010;27:747–52.View ArticlePubMed CentralPubMedGoogle Scholar
- Nei M, Tajima F, Tateno Y. Accuracy of estimated phylogenetic trees from molecular data. J Mol Evol. 1983;19:153–70.View ArticlePubMedGoogle Scholar
- Warren AD, Davis KJ, Grishin NV, Pelham JP, Stangeland EM: Interactive Listing of American Butterflies. [http://www.butterfliesofamerica.com]; 2010.
- Brown JW, Real HG, Faulkner DK. Butterflies of Baja California: faunal survey, natural history, conservation biology. Beverly Hills: The Lepidoptera Research Foundation; 1992.Google Scholar
- Brock JP, Kaufman K. Butterflies of North America. New York: Houghton Mifflin; 2003.Google Scholar
- Fisher MS. The Butterflies of Colorado: Riodinidae and Lycaenidae – Part 4: The Metalmarks, Coppers, Hairstreaks and Blues. Contributions of the Gillette Museum of Arthropod Diversity: Lepidoptera of North America title 74. Hansen Brothers Printing: Littleton, Colorado; 2009.Google Scholar
- Opler PA, Lotts K, Naberhaus T: Butterflies and Moths of North America. [http://www.butterfliesandmoths.org]; 2010.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.