Abstract
Phylogeographic patterns in phytophagous organisms are often contextualized in light of geographic isolation and ecological (host, habitat) specialization. However, assessing the relative impact of these phenomena is not straightforward, even in areas where phylogeography is well-studied, such as the California Floristic Province. Here, we use genome-wide markers to elucidate population genomic and phylgeographic patterns for a group of monophytophagous butterflies in southern California. This group is of high conservation interest because it includes the El Segundo blue, Euphilotes battoides allyni, one of the first insects listed under the U.S. Endangered Species Act, and a newly discovered population putatively assigned to E. b. allyni. Despite using the same unique host and coastal habitat, our results indicate that the newly discovered populations are not E. b. allyni and are more closely related to geographically proximate populations of the E. battoides group using a different habitat/host combination. Aside from E. b. allyni and the newly discovered populations, the rest of the group shows only fine-scale structure and apparently maintains genetic connectivity throughout southern California, across a vast range of habitats and climates, and on multiple hosts. Thus, habitat and host specialization did not elicit genetic isolation in neighboring populations suggesting that: (1) other phenomena are needed to explain the remarkable and idiosyncratic divergence of these highly restricted, proximate, taxa, and (2) fine-scale genomic markers suggest broader implications for understanding the mechanisms of speciation and reinvestigation of phylogeographic patterns in regions like the California Floristic Province.
Similar content being viewed by others
Data availability
Raw sequence data for ddRAD dataset is available on NCBI, BioProject PRJNA555218 (SRA accessions: SRR9707610-SRR9707926), and COI data is available on GenBank (accessions: MN224680-MN224984). Data files are available as a Dryad Data repository (repository: https://doi.org/10.5061/dryad.cvdncjt15). Data citation: Dupuis et al. (2019).
References
Allendorf FW (2017) Genetics and the conservation of natural populations: allozymes to genomes. Mol Ecol 26:420–430
Allendorf FW, Hohenlohe PA, Luikart G (2010) Genomics and the future of conservation genetics. Nat Rev Genet 11:697–709. https://doi.org/10.1038/nrg2844
Andrews KR, Good JM, Miller MR, Luikart G, Hohenlohe PA (2016) Harnessing the power of RADseq for ecological and evolutionary genomics. Nat Rev Genet 17:81–92. https://doi.org/10.1038/nrg.2015.28
Bowden R et al (2012) Genomic tools for evolution and conservation in the chimpanzee: Pan troglodytes ellioti is a genetically distinct population. PLoS Genet 8:e1002504. https://doi.org/10.1371/journal.pgen.1002504
Butts CT (2016) sna: tools for social network analysis. R package version 2.4. https://CRAN.R-project.org/package=sna
Calsbeek R, Thompson JN, Richardson JE (2003) Patterns of molecular evolution and diversification in a biodiversity hotspot: the California Floristic Province. Mol Ecol 12:1021–1029. https://doi.org/10.1046/j.1365-294X.2003.01794.x
Catchen JM, Amores A, Hohenlohe P, Cresko W, Postlethwait JH (2011) Stacks: building and genotyping loci de novo from short-read sequences. G3 (Bethesda) 1:171–182. https://doi.org/10.1534/g3.111.000240
Catchen JM, Hohenlohe PA, Bassham S, Amores A, Cresko WA (2013) Stacks: an analysis tool set for population genomics. Mol Ecol 22:3124–3140. https://doi.org/10.1111/mec.12354
Chatzimanolis S, Caterino MS (2007) Toward a better understanding of the "Transverse Range break": lineage diversification in southern California. Evolution 61:2127–2141. https://doi.org/10.1111/j.1558-5646.2007.00186.x
Clement M, Posada D, Crandall KA (2000) TCS: a computer program to estimate gene genealogies. Mol Ecol 9:1657–1659
Crandall KA, Bininda-Emonds ORP, Mace GM, Wayne RK (2000) Considering evolutionary processes in conservation biology. Trends Ecol Evol 15:290–295
Danecek P et al (2011) The variant call format and VCFtools. Bioinformatics 27:2156–2158. https://doi.org/10.1093/bioinformatics/btr330
DeAngelis MM, Wang DG, Hawkins TL (1995) Solid-phase reversible immobilization for the isolation of PCR products. Nucleic Acids Res 23:4742–4743
Decker J et al (2009) Resolving the evolution of extant and extinct ruminants with high-throughput phylogenomics. Proc Natl Acad Sci USA 106:18644–18649
Dobson AP, Rodriguez JP, Roberts WM, Wilcove DS (1997) Geographic distribution of endangered species in the United States. Science 275:550–553. https://doi.org/10.1126/science.275.5299.550
DOD (2008) Diverse communications project at Vandenberg Air Force Base Supplemental Environmental Assessment. Department of Defense Missile Defence Agency, Washington DC. https://www.mda.mil/global/documents/pdf/env_gmd_vafb_sea.pdf
Dupuis JR, Geib SM, Osborne KH, Rubinoff D (2019) Data from: genomics confirms surprising ecological divergence and isolation in an endangered butterfly
Dupuis JR, Oliver JC, Brunet BMT, Longcore T, Johnson JJ, Sperling FAH (2018a) Genomic data indicate ubiquitous evolutionary distinctiveness among populations of California metalmark butterflies. Conserv Genet 19:1097–1108. https://doi.org/10.1007/s10592-018-1081-8
Dupuis JR, Peigler RS, Geib SM, Rubinoff D (2018b) Phylogenomics supports incongruence between ecological specialization and taxonomy in a charismatic clade of buck moths. Mol Ecol. https://doi.org/10.1111/mec.14883
Durham JW, Allison EC (1960) The geologic history of Baja California and its marine faunas. Syst Zool 9:47–91
Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14:2611–2620. https://doi.org/10.1111/j.1365-294X.2005.02553.x
Excoffier L, Smouse PE, Quattro JM (1992) Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131:479–491
Falush D, Stephens M, Pritchard JK (2003) Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164:1567–1587
Feldman CR, Spicer GS (2006) Comparative phylogeography of woodland reptiles in California: repeated patterns of cladogenesis and population expansion. Mol Ecol 15:2201–2222. https://doi.org/10.1111/j.1365-294X.2006.02930.x
Feldman TD, Jonas AEG (2000) Sage scrub revolution? property rights, political fragmentation, and conservation planning in southern California under the federal Endangered Species Act. Ann Assoc Am Geogr 90:256–292
Fitzpatrick BM, Ryan ME, Johnson JR, Corush J, Carter ET (2015) Hybridization and the species problem in conservation. Curr Zool 61:206–216
Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R (1994) DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol 3:294–299
Fraser DJ, Bernatchez L (2001) Adaptive evolutionary conservation: towards a unified concept for defining conservation units. Mol Ecol 10:2741–2752
Gould B, Moeller DA, Eckhart VM, Tiffin P, Fabio E, Geber MA, Whitney K (2014) Local adaptation and range boundary formation in response to complex environmental gradients across the geographical range ofClarkia xantianassp. xantiana. J Ecol 102:95–107. https://doi.org/10.1111/1365-2745.12188
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol 59:307–321. https://doi.org/10.1093/sysbio/syq010
Haddad NM, Tilman D, Haarstad J, Ritchie M, Knops JMH (2001) Contrasting effects of plant richness and composition on insect communities: a field experiment. Am Nat 158:17–35
Hamilton JA, Miller JM (2016) Adaptive introgression as a resource for management and genetic conservation in a changing climate. Conserv Biol 30:33–41. https://doi.org/10.1111/cobi.12574
Hare MP (2001) Prospects for nuclear gene phylogeography. Trends Ecol Evol 16:700–706
Harrison SP (2013) Plant and animal endemism in California. University of California Press, Berkeley
Hawkins BA, Porter EE (2003) Does herbivore diversity depend on plant diversity? The case of California butterflies. Am Nat 161:40–49
Hijmans RJ (2017) geosphere: spherical trigonometry. R package version 1.5-7. https://CRAN.R-project.org/package=geosphere
Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS (2018) UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol 35:518–522. https://doi.org/10.1093/molbev/msx281
Hodel RGJ, Chen S, Payton AC, McDaniel SF, Soltis P, Soltis DE (2017) Adding loci improves phylogeographic resolution in red mangroves despite increased missing data: comparing microsatellites and RAD-Seq and investigating loci filtering. Sci Rep 7:17598. https://doi.org/10.1038/s41598-017-16810-7
Huson DH, Bryant D (2006) Application of phylogenetic networks in evolutionary studies. Mol Biol Evol 23:254–267
Jeffries DL, Copp GH, Lawson Handley L, Olsen KH, Sayer CD, Hanfling B (2016) Comparing RADseq and microsatellites to infer complex phylogeographic patterns, an empirical perspective in the Crucian carp, Carassius carassius, L. Mol Ecol 25:2997–3018. https://doi.org/10.1111/mec.13613
Jombart T (2008) adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24:1403–1405. https://doi.org/10.1093/bioinformatics/btn129
Jombart T, Devillard S, Balloux F (2010) Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet 11:94. https://doi.org/10.1186/1471-2156-11-94
Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS (2017) ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods 14:587–589. https://doi.org/10.1038/nmeth.4285
Kamvar ZN, Brooks JC, Grunwald NJ (2015) Novel R tools for analysis of genome-wide population genetic data with emphasis on clonality. Front Genet 6:208. https://doi.org/10.3389/fgene.2015.00208
Kamvar ZN, Tabima JF, Grunwald NJ (2014) Poppr: an R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ 2:e281. https://doi.org/10.7717/peerj.281
Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I (2015) Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour 15:1179–1191. https://doi.org/10.1111/1755-0998.12387
Kovach RP, Luikart G, Lowe WH, Boyer MC, Muhlfeld CC (2016) Risk and efficacy of human-enabled interspecific hybridization for climate-change adaptation: response to Hamilton and Miller (2016). Conserv Biol 30:428–430. https://doi.org/10.1111/cobi.12678
Lanfear R, Calcott B, Kainer D, Mayer C, Stamatakis A (2014) Selecting optimal partitioning schemes for phylogenomic datasets. BMC Evol Biol 14:82
Lapointe F-J, Rissler LJ (2005) Congruence, consensus, and the comparative phylogeography of codistributed species in California. Am Nat 166:290–299
Law JH, Crespi BJ (2002) The evolution of geographic parthenogenesis in Timema walking-sticks. Mol Ecol 11:1471–1489
Legendre P, Fortin M-J, Borcard D, Peres-Neto P (2015) Should the Mantel test be used in spatial analysis? Methods Ecol Evol 6:1239–1247. https://doi.org/10.1111/2041-210x.12425
Leigh JW, Bryant D, Nakagawa S (2015) popart: full-feature software for haplotype network construction. Methods Ecol Evol 6:1110–1116. https://doi.org/10.1111/2041-210x.12410
Lewis PO (2001) A likelihood approach to estimating phylogeny from discrete morphological character data. Syst Biol 50:913–925
Li YL, Liu JX (2017) StructureSelector: A web-based software to select and visualize the optimal number of clusters using multiple methods. Mol Ecol Resour. https://doi.org/10.1111/1755-0998.12719
Longcore T, Osborne KH (2015) Butterflies are not grizzly bears: lepidoptera conservation in practice. In: Daniels JC (ed) Butterfly conservation in North America. Springer, Netherlands, pp 161–192
Luikart G, England PR, Tallmon D, Jordan S, Taberlet P (2003) The power and promise of population genomics: from genotyping to genome typing. Nat Rev Genet 4:981–994. https://doi.org/10.1038/nrg1226
Maddison DR, Maddison WP (2018a) Chromaseq: a Mesquite package for analyzing sequence chromatograms. Version 1.31. https://chromaseq.mesquiteproject.org
Maddison WP, Maddison DR (2018b) Mesquite: a modular system for evolutionary analysis. Version 3.51. https://www.mesquiteproject.org
Marschalek DA, Klein MW (2010) Distribution, ecology, and conservation of Hermes copper (Lycaenidae: Lycaena [Hermelycaena] hermes). J Insect Conserv 14:721–730. https://doi.org/10.1007/s10841-010-9302-6
Mattoni R, Longcore T, Zonneveld C, Novotny V (2001) Analysis of transect counts to monitor population size in endangered insects: The case of the El Segundo blue butterfly, Euphilotes bernardino allyni. J Insect Conserv 5:197–206
Meirmans PG, Van Tienderen PH (2004) genotype and genodive: two programs for the analysis of genetic diversity of asexual organisms. Mol Ecol Notes 4:792–794. https://doi.org/10.1111/j.1471-8286.2004.00770.x
Michalakis Y, Excoffier L (1996) A generic estimation of population subdivision using distances between alleles with special reference for microsatellite loci. Genetics 142:1061–1064
Miller JM, Hamilton JA (2016) Interspecies hybridization in the conservation toolbox: response to. Conserv Biol 30:431–433. https://doi.org/10.1111/cobi.12677
Moritz C (1994a) Applications of mitochondrial DNA analysis in conservation: a critical review. Mol Ecol 3:401–411
Moritz C (1994b) Defining "evolutionarily significant units" for conservation. Trends Ecol Evol 9:373–375
Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GA, Kent J (2000) Biodiversity hotspots for conservation priorities. Nature 403:853–858. https://doi.org/10.1038/35002501
Narum SR, Buerkle CA, Davey JW, Miller MR, Hohenlohe PA (2013) Genotyping-by-sequencing in ecological and conservation genomics. Mol Ecol 22:2841–2847. https://doi.org/10.1111/mec.12350
Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ (2015) IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 32:268–274. https://doi.org/10.1093/molbev/msu300
Nosil P et al (2018) Natural selection and the predictability of evolution in Timema stick insects. Science 359:765–770
Ogden R et al (2013) Sturgeon conservation genomics: SNP discovery and validation using RAD sequencing. Mol Ecol 22:3112–3123. https://doi.org/10.1111/mec.12234
Opler PA, Wright AB (1999) Peterson field guides: a field guide to western butterflies. Houghton Mifflin, New York, New York
Pelham JP (2008) A catalogue of the butterflies of the United States and Canada, with a complete bibliography of the descriptive and systematic literature. J Res Lepidoptera 40:1–672
Peterson BK, Weber JN, Kay EH, Fisher HS, Hoekstra HE (2012) Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS ONE 7:e37135. https://doi.org/10.1371/journal.pone.0037135
Pfeifer B, Wittelsburger U, Ramos-Onsins SE, Lercher MJ (2014) PopGenome: an efficient Swiss army knife for population genomic analyses in R. Mol Biol Evol 31:1929–1936. https://doi.org/10.1093/molbev/msu136
Pratt GF (1988) The evolution and biology of Euphilotes biotypes. University of California at Riverside, Riverside
Pratt GF, Emmel JF (1998) Revision of the Euphilotes enoptes and E. battoides complexes (Lepidoptera: Lycaenidae). In: Emmel TC (ed) Systematics of Western North American Butterflies. Mariposa Press, Gainesville, Florida, pp 207–270
Pratt GF, Stouthamer R (2008) The genetic relationships between the El Segundo Blues from Los Angeles County and Santa Barbara County. Report. University of California at Riverside, Riverside
Preston KL, Rotenberry JT, Redak RA, Allen MF (2008) Habitat shifts of endangered species under altered climate conditions: importance of biotic interactions. Glob Change Biol. https://doi.org/10.1111/j.1365-2486.2008.01671.x
Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155:945–959
Pritchard JK, Wen X (2004) Documentation for structure software: version 2. University of Chicago Press, Chicago
Puechmaille SJ (2016) The program structure does not reliably recover the correct population structure when sampling is uneven: subsampling and new estimators alleviate the problem. Mol Ecol Resour 16:608–627. https://doi.org/10.1111/1755-0998.12512
Pyron RA, Costa GC, Patten MA, Burbrink FT (2015) Phylogenetic niche conservatism and the evolutionary basis of ecological speciation. Biol Rev Camb Philos Soc 90:1248–1262. https://doi.org/10.1111/brv.12154
R Core Team (2018) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-projectorg/
Rambaut A, Drummond AJ (2010) FigTree v1.4.4. Institute of Evolutionary Biology, University of Edinburgh. https://www.treebioedacuk/software/figtree
Reynolds J, Weir BS, Cockerham CC (1983) Estimation of the coancestry coefficient: basis for a short-term genetic distance. Genetics 105:767–779
Rissler LJ, Hijmans RJ, Graham CH, Moritz C, Wake DB (2006) Phylogeographic lineages and species comparisons in conservation analyses: a case studi of California herpetofauna. Am Nat 167:655–666
Rogers R, Mattoni R (1993) Observations on the natural history and conservation biology of the giant flower loving flies, Rhaphiomidas (Diptera: Apioceridae). Dipterol Res 4:21–34
Rubinoff D (2001) Evaluating the California gnatcatcher as an umbrella species for conservation of southern California coastal sage scrub. Conserv Biol 15:1374–1383
Rubinoff D, Cameron S, Will K (2006) A genomic perspective on the shortcomings of mitochondrial DNA for "barcoding" identification. J Hered 97:581–594. https://doi.org/10.1093/jhered/esl036
Rubinoff D, Sperling FAH (2004) Mitochondrial DNA sequence, morphology and ecology yield contrasting conservation implications for two threatened buckmoths (Hemileuca: Saturniidae). Biol Conserv 118:341–351. https://doi.org/10.1016/j.biocon.2003.09.013
Sandoval C, Carmean DA, Crespi BJ (1998) Molecular phylogenetics of sexual and parthenogenic Timema walking-sticks. Proc R Soc B 265:589–595
Savolainen O, Lascoux M, Merila J (2013) Ecological genomics of local adaptation. Nat Rev Genet 14:807–820. https://doi.org/10.1038/nrg3522
Schierenbeck KA (2014) Phylogeography of California: an introduction. University of California Press, Oakland, California
Scott JA (1986) The Butterflies of North America: a natural history and field guide. Stanford University Press, Stanford, CA
Shields O, Reveal JL (1988) Sequential evolution of Euphilotes (Lycaenidae: Scolitantidini) on their plant host Eriogonum (Polygonaceae: Eriogonoideae). Biol J Linn Soc 33:51–93
Siemann E, Tilman D, Haarstad J, Ritchie M (1998) Experimental tests of the dependence of Arthropod diversity on plant diversity. Am Nat 152:738–750
Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P (1994) Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am 87:651–701
Sork VL, Gugger PF, Chen JM, Werth S (2016) Evolutionary lessons from California plant phylogeography. Proc Natl Acad Sci USA 113:8064–8071. https://doi.org/10.1073/pnas.1602675113
Spinks PQ, Thomson RC, Shaffer HB (2010) Nuclear gene phylogeography reveals the historical legacy of an ancient inland sea on lineages of the western pond turtle, Emys marmorata in California. Mol Ecol 19:542–556. https://doi.org/10.1111/j.1365-294X.2009.04451.x
Templeton AR, Crandall KA, Sing CF (1992) A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics 132:619–633
The Inkscape Team (2017) Inkscape v0.91. https://www.inkscape.org
Toews DP, Brelsford A (2012) The biogeography of mitochondrial and nuclear discordance in animals. Mol Ecol 21:3907–3930. https://doi.org/10.1111/j.1365-294X.2012.05664.x
Ugelvig LV, Vila R, Pierce NE, Nash DR (2011) A phylogenetic revision of the Glaucopsyche section (Lepidoptera: Lycaenidae), with special focus on the Phengaris-Maculinea clade. Mol Phylogenet Evol 61:237–243. https://doi.org/10.1016/j.ympev.2011.05.016
USFWS (2008) El Segundo Blue Butterfly (Euphilotes battoides allyni) 5-year review: summary and evaluation. United States Fish and Wildlife Service, Carlsbad Fish and Wildlife Office, Carlsbad, CA. https://www.fws.gov/carlsbad/SpeciesStatusList/5YR/20080331_5YR_ESB.pdf
Vähä JP, Erkinaro J, Niemela E, Primmer CR (2007) Life-history and habitat features influence the within-river genetic structure of Atlantic salmon. Mol Ecol 16:2638–2654. https://doi.org/10.1111/j.1365-294X.2007.03329.x
Wagner CE et al (2013) Genome-wide RAD sequence data provide unprecedented resolution of species boundaries and relationships in the Lake Victoria cichlid adaptive radiation. Mol Ecol 22:787–798. https://doi.org/10.1111/mec.12023
Wan QH, Wu H, Fujihara T, Fang SG (2004) Which genetic marker for which conservation genetics issue? Electrophoresis 25:2165–2176. https://doi.org/10.1002/elps.200305922
Zedler JB (1996) Coastal mitigation in southern California: the need for a regional restoration strategy. Ecol Appl 6:84–93
Acknowledgements
This project would not have been possible without the enthusiasm, assistance, and advice of Nancy Ferguson (United States Fish and Wildlife Service), and support by the USFWS and VAFB. Eric Porter and William Miller (USFWS), and Travis Longcore (Urban Wildlands Group) provided advice and support necessary for our success. We thank the following individuals for assistance in collecting specimens: Alice Abela (ManTech), Sam Fisher, Robert Fisher, Marie Fisher, Carolyn Lin (Los Angeles World Airports), Adrienne Mohan (Palos Verdes Peninsula Land Conservancy), Katrina Olthof (Wildlands Conservation Science), Nancy Price (Los Angeles World Airports), and Patrick Tyrrell (Friends of Balona Wetlands). We also thank M. Wells, Torrey Pines State Reserve; D. Lydy U.S. Navy, Pt. Loma; T. Walker of the Tri-canyon Parks (San Diego), Environmental Protection Services, U.S. Marine Corps Camp Pendelton; and especially S. Weber, Chief of Natural Resources at Cabrillo National Monument. D.R. thanks J. and P. Brown and family for support and invaluable advice. J. Vanderweir and E. Hein, U.S. Fish and Wildlife Service provided essential technical advice and assistance for collection of San Diego County samples. ddRAD library preparation and wet lab procedures were performed at the United States Department of Agriculture—Agricultural Research Service (USDA-ARS) Daniel K. Inouye US PBARC Genomics facility by Angela Kauwe and Kimberley Morris, and HiSeq sequencing was conducted at the Vincent J. Coates Genomics Sequencing Laboratory at University of California at Berkeley, supported by National Institutes of Health S10 Instrumentation Grants S10RR029668 and S10RR027303. Phylogenetic data file with flanking sequence produced using https://github.com/jrdupuis/stacks_full_fasta_out, and figures produced using Inkscape v0.91 (The Inkscape Team 2017) and R (R Core Team 2018). USDA is an equal opportunity employer. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the USDA.
Author information
Authors and Affiliations
Contributions
KHO and DR conceptualized study with input from JRD and SMG. KHO and DR collected specimens and JRD conducted data processing and analyses with input from SMG. JRD and DR wrote the manuscript with input from SMG and KHO.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Ethical approval
All ethical standards have been followed.
Additional information
Communicated by Nigel E. Stork.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
10531_2020_1950_MOESM1_ESM.pdf
Figure S1. Full results of STRUCTURE analyses for all individuals (a), and the allyni (b), VAFB (c), and battoides (d) clusters, including plots of LnPr(K|X), ΔK, and the Peuchmaille statistics, and barplots for all values of K. Individuals labeled at an optimal value of K, if appropriate. Alternative sort order (by general area within the VAFB (northern, southern, off-base) and by longitude) provided for VAFB cluster below (c). (PDF 2316 kb)
10531_2020_1950_MOESM2_ESM.pdf
Figure S2. Results of DAPC for all individuals (a), the allyni and VAFB clusters (b), and the battoides cluster (c). Inset graphs show plots of number of clusters versus BIC and the relative contribution of each discriminant axis to the analysis. In (b), find.clusters predicted K = 1, so DAPC plots not shown; however in (c), K = 2 is forced in the final iterative hierarchical analysis of the battoides cluster to show partial separation of San Diego populations from San Bernardino/Riverside populations. (PDF 324 kb)
10531_2020_1950_MOESM3_ESM.pdf
Figure S3. ML consensus tree from phylogenetic dataset (260 individuals, 54,305 SNPs), with SH-aLRT/ufBS branch support displayed (strong support = >0.8/>0.95). (PDF 16 kb)
10531_2020_1950_MOESM4_ESM.pdf
Figure S4. ML consensus tree from supplementary phylogenetic dataset 1 (260 individuals, 1,882 SNPs matching the population genetic dataset), with SH-aLRT/ufBS branch support displayed (strong support = >0.8/>0.95). (PDF 16 kb)
10531_2020_1950_MOESM5_ESM.pdf
Figure S5. ML consensus tree from supplementary phylogenetic dataset 2 (260 individuals, 189,837 characters, matching the population genetic dataset but with flanking sequence around each SNP), with SH-aLRT/ufBS branch support displayed (strong support = >0.8/>0.95). (PDF 16 kb)
10531_2020_1950_MOESM6_ESM.pdf
Figure S6. ML consensus tree from COI dataset (305 individuals, 650 bp), with SH-aLRT/ufBS branch support displayed (strong support = >0.8/>0.95). (PDF 14 kb)
Rights and permissions
About this article
Cite this article
Dupuis, J.R., Geib, S.M., Osborne, K.H. et al. Genomics confirms surprising ecological divergence and isolation in an endangered butterfly. Biodivers Conserv 29, 1897–1921 (2020). https://doi.org/10.1007/s10531-020-01950-6
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10531-020-01950-6