Parthenogenesis in an elasmobranch in the presence of conspecific males.

Parthenogenesis has been observed in several elasmobranch species, primarily in public aquaria. The majority of cases of parthenogenesis have occurred either when females were held without males or once a male was removed from a female's habitat. Here we report a second instance of parthenogenesis in a zebra shark female that was housed with conspecific mature males. This study calls into question the conditions under which elasmobranch females undergo parthenogenesis.

Phylogeography and population genetics of the cryptic bonnethead shark Sphyrna aff. tiburo in Brazil and the Caribbean inferred from mtDNA markers.

Resolving the identity, phylogeny and distribution of cryptic species within species complexes is an essential precursor to management. The bonnethead shark, Sphyrna tiburo, is a small coastal shark distributed in the Western Atlantic from North Carolina (U.S.A.) to southern Brazil. Genetic analyses based on mitochondrial markers revealed that bonnethead sharks comprise a species complex with at least two lineages in the Northwestern Atlantic and the Caribbean (S. tiburo and Sphyrna aff. tiburo, respectively). The phylogeographic and phylogenetic analysis of two mitochondrial markers [control region (mtCR) and cytochrome oxidase I (COI)] showed that bonnethead sharks from southeastern Brazil correspond to S. aff. tiburo, extending the distribution of this cryptic species >5000 km. Bonnethead shark populations are only managed in the U.S.A. and in the 2000s were considered to be regionally extinct or collapsed in southeast Brazil. The results indicate that there is significant genetic differentiation between S. aff. tiburo from Brazil and other populations from the Caribbean (ΦST  = 0.9053, P < 0.000), which means that collapsed populations in the former are unlikely to be replenished from Caribbean immigration. The species identity of bonnethead sharks in the Southwest Atlantic and their relationship to North Atlantic and Caribbean populations still remains unresolved. Taxonomic revision and further sampling are required to reevaluate the status of the bonnethead shark complex through its distribution range.

Microsatellite primer development in elasmobranchs using next generation sequencing of enriched libraries.

Microsatellites are useful in studies of population genetics, sibship, and parentage. Here, we screened for microsatellites from multiple elasmobranch genomic libraries using an enrichment protocol followed by sequencing on an Illumina platform. We concurrently screened five and then nine genomes and describe the number of potential loci from each respective round of sequencing. To validate the efficacy of the protocol, we developed and tested primers for the pelagic thresher shark, Alopias pelagicus. The method described here is a cost-effective protocol to increase the pool of potential useful loci and allows the concurrent screening of multiple libraries.

Annual breeding in a captive smalltooth sawfish, Pristis pectinata.

The critically endangered smalltooth sawfish Pristis pectinata reproduces biennially in central west Florida, U.S.A. Here we demonstrate that smalltooth sawfish are physiologically capable of reproducing annually in a captive environment. The smalltooth sawfish are held in an open system, with abiotic conditions that vary naturally with the surrounding environment in The Bahamas. This suggests wild smalltooth sawfish may also be capable of annual reproduction provided there are adequate prey resources, limited competition and mate availability.

No apparent negative tagging effects after 13 years at liberty for lemon shark, Negaprion brevirostris implanted with acoustic transmitter.

An intact and uncompromised internal acoustic transmitter was non-lethally recovered from a lemon shark Negaprion brevirostris, after 13 years at liberty. The shark, first tagged at an estimated age of 2 years old near South Bimini, Bahamas in 2004, was recaptured in 2017 with a total length of 264 cm. The tagged shark displayed typical growth rate, pregnancy, natal homing and pupping behaviour of other individuals in this population. This observation provides important evidence regarding the effects from long-term retention of implanted acoustic transmitters in a carcharhinid shark.

Species composition of the international shark fin trade assessed through a retail-market survey in Hong Kong.

The shark fin trade is a major driver of shark exploitation in fisheries all over the world, most of which are not managed on a species-specific basis. Species-specific trade information highlights taxa of particular concern and can be used to assess the efficacy of management measures and anticipate emerging threats. The species composition of the Hong Kong Special Administrative Region of China, one of the world's largest fin trading hubs, was partially assessed in 1999-2001. We randomly selected and genetically identified fin trimmings (n = 4800), produced during fin processing, from the retail market of Hong Kong in 2014-2015 to assess contemporary species composition of the fin trade. We used nonparametric species estimators to determine that at least 76 species of sharks, batoids, and chimaeras supplied the fin trade and a Bayesian model to determine their relative proportion in the market. The diversity of traded species suggests species substitution could mask depletion of vulnerable species; one-third of identified species are threatened with extinction. The Bayesian model suggested that 8 species each comprised >1% of the fin trimmings (34.1-64.2% for blue [Prionace glauca], 0.2-1.2% for bull [Carcharhinus leucas] and shortfin mako [Isurus oxyrinchus]); thus, trade was skewed to a few globally distributed species. Several other coastal sharks, batoids, and chimaeras are in the trade but poorly managed. Fewer than 10 of the species we modeled have sustainably managed fisheries anywhere in their range, and the most common species in trade, the blue shark, was not among them. Our study and approach serve as a baseline to track changes in composition of species in the fin trade over time to better understand patterns of exploitation and assess the effects of emerging management actions for these animals.

Risky business for a juvenile marine predator? Testing the influence of foraging strategies on size and growth rate under natural conditions.

Mechanisms driving selection of body size and growth rate in wild marine vertebrates are poorly understood, thus limiting knowledge of their fitness costs at ecological, physiological and genetic scales. Here, we indirectly tested whether selection for size-related traits of juvenile sharks that inhabit a nursery hosting two dichotomous habitats, protected mangroves (low predation risk) and exposed seagrass beds (high predation risk), is influenced by their foraging behaviour. Juvenile sharks displayed a continuum of foraging strategies between mangrove and seagrass areas, with some individuals preferentially feeding in one habitat over another. Foraging habitat was correlated with growth rate, whereby slower growing, smaller individuals fed predominantly in sheltered mangroves, whereas larger, faster growing animals fed over exposed seagrass. Concomitantly, tracked juveniles undertook variable movement behaviours across both the low and high predation risk habitat. These data provide supporting evidence for the hypothesis that directional selection favouring smaller size and slower growth rate, both heritable traits in this shark population, may be driven by variability in foraging behaviour and predation risk. Such evolutionary pathways may be critical to adaptation within predator-driven marine ecosystems.

Global population genetic dynamics of a highly migratory, apex predator shark.

Knowledge of genetic connectivity dynamics in the world's large-bodied, highly migratory, apex predator sharks across their global ranges is limited. One such species, the tiger shark (Galeocerdo cuvier), occurs worldwide in warm temperate and tropical waters, uses remarkably diverse habitats (nearshore to pelagic) and possesses a generalist diet that can structure marine ecosystems through top-down processes. We investigated the phylogeography and the global population structure of this exploited, phylogenetically enigmatic shark by using 10 nuclear microsatellites (n = 380) and sequences from the mitochondrial control region (CR, n = 340) and cytochrome oxidase I gene (n = 100). All three marker classes showed the genetic differentiation between tiger sharks from the western Atlantic and Indo-Pacific ocean basins (microsatellite FST  > 0.129; CR ΦST  > 0.497), the presence of North vs. southwestern Atlantic differentiation and the isolation of tiger sharks sampled from Hawaii from other surveyed locations. Furthermore, mitochondrial DNA revealed high levels of intraocean basin matrilineal population structure, suggesting female philopatry and sex-biased gene flow. Coalescent- and genetic distance-based estimates of divergence from CR sequences were largely congruent (dcorr  = 0.0015-0.0050), indicating a separation of Indo-Pacific and western Atlantic tiger sharks <1 million years ago. Mitochondrial haplotype relationships suggested that the western South Atlantic Ocean was likely a historical connection for interocean basin linkages via the dispersal around South Africa. Together, the results reveal unexpectedly high levels of population structure in a highly migratory, behaviourally generalist, cosmopolitan ocean predator, calling for management and conservation on smaller-than-anticipated spatial scales.

Surgeons and suture zones: Hybridization among four surgeonfish species in the Indo-Pacific with variable evolutionary outcomes.

Closely related species can provide valuable insights into evolutionary processes through comparison of their ecology, geographic distribution and the history recorded in their genomes. In the Indo-Pacific, many reef fishes are divided into sister species that come into secondary contact at biogeographic borders, most prominently where Indian Ocean and Pacific Ocean faunas meet. It is unclear whether hybridization in this contact zone represents incomplete speciation, secondary contact, an evolutionary dead-end (for hybrids) or some combination of the above. To address these issues, we conducted comprehensive surveys of two widely-distributed surgeonfish species, Acanthurus leucosternon (N=141) and A. nigricans (N=412), with mtDNA cytochrome b sequences and ten microsatellite loci. These surgeonfishes are found primarily in the Indian and Pacific Oceans, respectively, but overlap at the Christmas and Cocos-Keeling Islands hybrid zone in the eastern Indian Ocean. We also sampled the two other Pacific members of this species complex, A. achilles (N=54) and A. japonicus (N=49), which are known to hybridize with A. nigricans where their ranges overlap. Our results indicate separation between the four species that range from the recent Pleistocene to late Pliocene (235,000-2.25million years ago). The Pacific A. achilles is the most divergent (and possibly ancestral) species with mtDNA dcorr≈0.04, whereas the other two Pacific species (A. japonicus and A. nigricans) are distinguishable only at a population or subspecies level (ΦST=0.6533, P<0.001). Little population structure was observed within species, with evidence of recent population expansion across all four geographic ranges. We detected sharing of mtDNA haplotypes between species and extensive hybridization based on microsatellites, consistent with later generation hybrids but also the effects of allele homoplasy. Despite extensive introgression, 98% of specimens had concordance between mtDNA lineage and species identification based on external morphology, indicating that species integrity may not be eroding. The A. nigricans complex demonstrates a range of outcomes from incomplete speciation to secondary contact to decreasing hybridization with increasing evolutionary depth.

Population structure and dispersal of the coral-excavating sponge Cliona delitrix.

Some excavating sponges of the genus Cliona compete with live reef corals, often killing and bioeroding entire colonies. Important aspects affecting distribution of these species, such as dispersal capability and population structure, remain largely unknown. Thus, the aim of this study was to determine levels of genetic connectivity and dispersal of Cliona delitrix across the Greater Caribbean (Caribbean Sea, Bahamas and Florida), to understand current patterns and possible future trends in their distribution and effects on coral reefs. Using ten species-specific microsatellite markers, we found high levels of genetic differentiation between six genetically distinct populations: one in the Atlantic (Florida-Bahamas), one specific to Florida and four in the South Caribbean Sea. In Florida, two independent breeding populations are likely separated by depth. Gene flow and ecological dispersal occur among other populations in the Florida reef tract, and between some Florida locations and the Bahamas. Similarly, gene flow occurs between populations in the South Caribbean Sea, but appears restricted between the Caribbean Sea and the Atlantic (Florida-Bahamas). Dispersal of C. delitrix was farther than expected for a marine sponge and favoured in areas where currents are strong enough to transport sponge eggs or larvae over longer distances. Our results support the influence of ocean current patterns on genetic connectivity, and constitute a baseline to monitor future C. delitrix trends under climate change.

Genetic Diversity of White Sharks, Carcharodon carcharias, in the Northwest Atlantic and Southern Africa.

The white shark, Carcharodon carcharias, is both one of the largest apex predators in the world and among the most heavily protected marine fish. Population genetic diversity is in part shaped by recent demographic history and can thus provide information complementary to more traditional population assessments, which are difficult to obtain for white sharks and have at times been controversial. Here, we use the mitochondrial control region and 14 nuclear-encoded microsatellite loci to assess white shark genetic diversity in 2 regions: the Northwest Atlantic (NWA, N = 35) and southern Africa (SA, N = 131). We find that these 2 regions harbor genetically distinct white shark populations (Φ ST = 0.10, P < 0.00001; microsatellite F ST = 0.1057, P < 0.021). M-ratios were low and indicative of a genetic bottleneck in the NWA (M-ratio = 0.71, P < 0.004) but not SA (M-ratio = 0.85, P = 0.39). This is consistent with other evidence showing a steep population decline occurring in the mid to late 20th century in the NWA, whereas the SA population appears to have been relatively stable. Estimates of effective population size ranged from 22.6 to 66.3 (NWA) and 188 to 1998.3 (SA) and evidence of inbreeding was found (primarily in NWA). Overall, our findings indicate that white population dynamics within NWA and SA are determined more by intrinsic reproduction than immigration and there is genetic evidence of a population decline in the NWA, further justifying the strong domestic protective measures that have been taken for this species in this region. Our study also highlights how assessment of genetic diversity can complement other sources of information to better understand the status of threatened marine fish populations.

Two decades of genetic profiling yields first evidence of natal philopatry and long-term fidelity to parturition sites in sharks.

Sharks are a globally threatened group of marine fishes that often breed in their natal region of origin. There has even been speculation that female sharks return to their exact birthplace to breed ('natal philopatry'), which would have important conservation implications. Genetic profiling of lemon sharks (Negaprion brevirostris) from 20 consecutive cohorts (1993-2012) at Bimini, Bahamas, showed that certain females faithfully gave birth at this site for nearly two decades. At least six females born in the 1993-1997 cohorts returned to give birth 14-17 years later, providing the first direct evidence of natal philopatry in the chondrichthyans. Long-term fidelity to specific nursery sites coupled with natal philopatry highlights the merits of emerging spatial and local conservation efforts for these threatened predators.

Comparative population genetics and evolutionary history of two commonly misidentified billfishes of management and conservation concern.

BACKGROUND: Misidentifications between exploited species may lead to inaccuracies in population assessments, with potentially irreversible conservation ramifications if overexploitation of either species is occurring. A notable showcase is provided by the realization that the roundscale spearfish (Tetrapturus georgii), a recently validated species, has been historically misidentified as the morphologically very similar and severely overfished white marlin (Kajikia albida) (IUCN listing: Vulnerable). In effect, no information exists on the population status and evolutionary history of the enigmatic roundscale spearfish, a large, highly vagile and broadly distributed pelagic species. We provide the first population genetic evaluation of the roundscale spearfish, utilizing nuclear microsatellite and mitochondrial DNA sequence markers. Furthermore, we re-evaluated existing white marlin mitochondrial genetic data and present our findings in a comparative context to the roundscale spearfish. RESULTS: Microsatellite and mitochondrial (control region) DNA markers provided mixed evidence for roundscale spearfish population differentiation between the western north and south Atlantic regions, depending on marker-statistical analysis combination used. Mitochondrial DNA analyses provided strong signals of historical population growth for both white marlin and roundscale spearfish, but higher genetic diversity and effective female population size (1.5-1.9X) for white marlin. CONCLUSIONS: The equivocal indications of roundscale spearfish population structure, combined with a smaller effective female population size compared to the white marlin, already a species of concern, suggests that a species-specific and precautionary management strategy recognizing two management units is prudent for this newly validated billfish.

Insights into the nuclear and mitochondrial genome of the Lemon shark Negaprion brevirostris using low-coverage sequencing: Genome size, repetitive elements, mitochondrial genome, and phylogenetic placement.

The Lemon shark Negaprion brevirostris is an important species experiencing conservation issues that is in need of genomic resources. Herein, we conducted a genome survey sequencing in N. brevirostris and determined genome size, explored repetitive elements, assembled and annotated the 45S rRNA DNA operon, and assembled and described in detail the mitochondrial genome. Lastly, the phylogenetic position of N. brevirostris in the family Carcharhinidae was examined using translated protein coding genes. The estimated haploid genome size ranged between 2.29 and 2.58 Gbp using a k-mer analysis, which is slightly below the genome size estimated for other sharks belonging to the family Carcharhinidae. Using a k-mer analysis, approx. 64-71 % of the genome of N. brevirostris was composed of repetitive elements. A relatively large proportion of the 'repeatome' could not be annotated. Taking into account only annotated repetitive elements, Class I - Long Interspersed Nuclear Element (LINE) were the most abundant repetitive elements followed by Class I - Penelope and Satellite DNA. The nuclear ribosomal operon was fully assembled. The AT-rich complete mitochondrial genome was 16,703 bp long and encoded 13 protein coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Negaprion brevirostris is closely related to the genera Carcharhinus, Glyphis and Lamiopsis in the family Carcharinidae. This new genomic resources will aid with the development of conservation plans for this large coastal shark.

Accounting for unobserved population dynamics and aging error in close-kin mark-recapture assessments.

Obtaining robust estimates of population abundance is a central challenge hindering the conservation and management of many threatened and exploited species. Close-kin mark-recapture (CKMR) is a genetics-based approach that has strong potential to improve the monitoring of data-limited species by enabling estimates of abundance, survival, and other parameters for populations that are challenging to assess. However, CKMR models have received limited sensitivity testing under realistic population dynamics and sampling scenarios, impeding the application of the method in population monitoring programs and stock assessments. Here, we use individual-based simulation to examine how unmodeled population dynamics and aging uncertainty affect the accuracy and precision of CKMR parameter estimates under different sampling strategies. We then present adapted models that correct the biases that arise from model misspecification. Our results demonstrate that a simple base-case CKMR model produces robust estimates of population abundance with stable populations that breed annually; however, if a population trend or non-annual breeding dynamics are present, or if year-specific estimates of abundance are desired, a more complex CKMR model must be constructed. In addition, we show that CKMR can generate reliable abundance estimates for adults from a variety of sampling strategies, including juvenile-focused sampling where adults are never directly observed (and aging error is minimal). Finally, we apply a CKMR model that has been adapted for population growth and intermittent breeding to two decades of genetic data from juvenile lemon sharks (Negaprion brevirostris) in Bimini, Bahamas, to demonstrate how application of CKMR to samples drawn solely from juveniles can contribute to monitoring efforts for highly mobile populations. Overall, this study expands our understanding of the biological factors and sampling decisions that cause bias in CKMR models, identifies key areas for future inquiry, and provides recommendations that can aid biologists in planning and implementing an effective CKMR study, particularly for long-lived data-limited species.

Gene expression, evolution, and the genetics of electrosensing in the smalltooth sawfish, Pristis pectinata.

Sawfishes (Pristidae) are large, highly threatened rays named for their tooth-studded rostrum, which is used for prey sensing and capture. Of all five species, the smalltooth sawfish, Pristis pectinata, has experienced the greatest decline in range, currently found in only ~20% of its historic range. To better understand the genetic underpinnings of these taxonomically and morphologically unique animals, we collected transcriptomic data from several tissue types, mapped them to the recently completed reference genome, and contrasted the patterns observed with comparable data from other elasmobranchs. Evidence of positive selection was detected in 79 genes in P. pectinata, several of which are involved in growth factor/receptor tyrosine kinase signaling and body symmetry and may be related to the unique morphology of sawfishes. Changes in these genes may impact cellular responses to environmental conditions such as temperature, dissolved oxygen, and salinity. Data acquired also allow for examination of the molecular components of P. pectinata electrosensory systems, which are highly developed in sawfishes and have likely been influential in their evolutionary success.

Vectored dispersal of Symbiodinium by larvae of a Caribbean gorgonian octocoral.

The ability of coral reefs to recover from natural and anthropogenic disturbance is difficult to predict, in part due to uncertainty regarding the dispersal capabilities and connectivity of their reef inhabitants. We developed microsatellite markers for the broadcast spawning gorgonian octocoral Eunicea (Plexaura) flexuosa (four markers) and its dinoflagellate symbiont, Symbiodinium B1 (five markers), and used them to assess genetic connectivity, specificity and directionality of gene flow among sites in Florida, Panama, Saba and the Dominican Republic. Bayesian analyses found that most E. flexuosa from the Florida reef tract, Saba and the Dominican Republic were strongly differentiated from many E. flexuosa in Panama, with the exception of five colonies from Key West that clustered with colonies from Panama. In contrast, Symbiodinium B1 was more highly structured. At least seven populations were detected that showed patterns of isolation by distance. The symbionts in the five unusual Key West colonies also clustered with symbionts from Panama, suggesting these colonies are the result of long-distance dispersal. Migration rate tests indicated a weak signal of northward immigration from the Panama population into the lower Florida Keys. As E. flexuosa clonemates only rarely associated with the same Symbiodinium B1 genotype (and vice versa), these data suggest a dynamic host-symbiont relationship in which E. flexuosa is relatively well dispersed but likely acquires Symbiodinium B1 from highly structured natal areas prior to dispersal. Once vectored by host larvae, these symbionts may then spread through the local population, and/or host colonies may acquire different local symbiont genotypes over time.

The behavioural and genetic mating system of the sand tiger shark, Carcharias taurus, an intrauterine cannibal.

Sand tiger sharks (Carcharias taurus) have an unusual mode of reproduction, whereby the first embryos in each of the paired uteri to reach a certain size ('hatchlings') consume all of their smaller siblings during gestation ('embryonic cannibalism' or EC). If females commonly mate with multiple males ('behavioural polyandry') then litters could initially have multiple sires. It is possible, however, that EC could exclude of all but one of these sires from producing offspring thus influencing the species genetic mating system ('genetic monogamy'). Here, we use microsatellite DNA profiling of mothers and their litters (n = 15, from two to nine embryos per litter) to quantify the frequency of behavioural and genetic polyandry in this system. We conservatively estimate that nine of the females we examined (60%) were behaviourally polyandrous. The genetic mating system was characterized by assessing sibling relationships between hatchlings and revealed only 40 per cent genetic polyandry (i.e. hatchlings were full siblings in 60% of litters). The discrepancy stemmed from three females that were initially fertilized by multiple males but only produced hatchlings with one of them. This reveals that males can be excluded even after fertilizing ova and that some instances of genetic monogamy in this population arise from the reduction in litter size by EC. More research is needed on how cryptic post-copulatory and post-zygotic processes contribute to determining paternity and bridging the behavioural and genetic mating systems of viviparous species.

Phylogeography of two closely related Indo-Pacific butterflyfishes reveals divergent evolutionary histories and discordant results from mtDNA and microsatellites.

Marine biogeographic barriers can have unpredictable consequences, even among closely related species. To resolve phylogeographic patterns for Indo-Pacific reef fauna, we conducted range-wide surveys of sister species, the scrawled butterflyfish (Chaetodon meyeri; N = 134) and the ornate butterflyfish (Chaetodon ornatissimus; N = 296), using mitochondrial DNA cytochrome b sequences and 10 microsatellite loci. The former is distributed primarily in the Indian Ocean but also extends to the Line Islands in the Central Pacific, whereas the latter is distributed primarily in the Central-West Pacific (including Hawaii and French Polynesia) but extends to the eastern margin of the Indian Ocean. Analyses of molecular variance and Bayesian STRUCTURE results revealed 1 range-wide group for C. meyeri and 3 groups for C. ornatissimus: 1) eastern Indian Ocean and western Pacific, 2) Central Pacific, and 3) Hawaii. Estimates of the last population expansion were much more recent for C. meyeri (61 500 to 95 000 years) versus C. ornatissimus (184 700 to 286 300 years). Despite similarities in ecology, morphology, life history, and a broadly overlapping distribution, these sister species have divergent patterns of dispersal and corresponding evolutionary history. The mtDNA and microsatellite markers did not provide concordant results within 1 of our study species (C. meyeri), or in 7 out of 12 other cases of marine fishes in the published literature. This discordance renews caution in relying on one or a few markers for reconstructing historical demography.

Deep ecomorphological and genetic divergence in Steller's Jays (Cyanocitta stelleri, Aves: Corvidae).

The relationship between ecology and morphology is a cornerstone of evolutionary biology, and quantifying variation across environments can shed light on processes that give rise to biodiversity. Three morphotypes of the Steller's Jay (Cyanocitta stelleri) occupy different ecoregions in western North America, which vary in climate and landcover. These morphotypes (Coastal, Interior, Rocky Mountain) differ in size, plumage coloration, and head pattern. We sampled 1080 Steller's Jays from 68 populations (plus 11 outgroups) to address three main questions using data on morphology, plumage, genetics (mtDNA, microsatellites), and ecological niches: (1) How do phenotypic and genetic traits vary within and among populations, morphotypes, and ecoregions? (2) How do population-level differences in Steller's Jays compare with other sister species pairs of North American birds? (3) What can we infer about the population history of Steller's Jays in relation to past climates, paleoecology, and niche evolution? We found substantial morphological, genetic, and ecological differentiation among morphotypes. The greatest genetic divergence separated Coastal and Interior morphotypes from the Rocky Mountain morphotype, which was associated with warmer, drier, and more open habitats. Microsatellites revealed additional structure between Coastal and Interior groups. The deep mtDNA split between Coastal/Interior and Rocky Mountain lineages of Steller's Jay (ND2 ~ 7.8%) is older than most North American avian sister species and dates to approximately 4.3 mya. Interior and Rocky Mountain morphotypes contact across a narrow zone with steep clines in traits and reduced gene flow. The distribution of the three morphotypes coincides with divergent varieties of ponderosa pine and Douglas fir. Species distribution models support multiple glacial refugia for Steller's Jays. Our integrative dataset combined with extensive geographic sampling provides compelling evidence for recognizing at least two species of Steller's Jay.