Lamb T, Justice TC, Brewer MS, Moler PE, Hopkins H, Bond JE. A biogeographical profile of the sand cockroach *Arenivaga floridensis* and its bearing on origin hypotheses for Florida scrub biota. Ecol Evol. 2018;8:5254--5266. <https://doi.org/10.1002/ece3.3885>

1. INTRODUCTION {#ece33885-sec-0001}
===============

Florida scrub is a fragmented xeric ecosystem largely confined to peninsular Florida (Menges, [1999](#ece33885-bib-0032){ref-type="ref"}), where it is partitioned across a series of relict beach ridges that formed sequentially during the Miocene, Pliocene, and Pleistocene epochs (Scott, [1997](#ece33885-bib-0048){ref-type="ref"}; Figure [1](#ece33885-fig-0001){ref-type="fig"}). Older ridges occupy the central peninsula, whereas the youngest corresponds roughly with present shorelines (Opdyke, Spangler, Smith, Jones, & Lindquist, [1984](#ece33885-bib-0041){ref-type="ref"}). Following inception, most ridges experienced rounds of inundation and isolation associated with sea level fluctuation (Webb, [1990](#ece33885-bib-0058){ref-type="ref"}). Today, these ridges retain scrub and related sandhill ecosystems, which were considered to have been widespread on the peninsula in the late Pleistocene (Myers, [1990](#ece33885-bib-0037){ref-type="ref"}) but experienced significant contraction to the ridges proper under more mesic conditions of the Holocene (Watts & Hansen, [1994](#ece33885-bib-0057){ref-type="ref"}).

![Map showing the Florida peninsula\'s sand ridge system, with sampled ridges depicted in color, and images of male and female (wingless) *Arenivaga floridensis*, illustrating species sexual dimorphism](ECE3-8-5254-g001){#ece33885-fig-0001}

Florida\'s sand ridges are characterized by quartzipsamment soils, providing porosity and drainage necessary to support the xeromorphic plant community that defines scrub---a shrubland composed of small evergreen oaks (*Quercus chapmanii* Sargent, *Q. geminata* Small, *Q. inopina* Ashe, *Q. myrtifolia* Willdenow) interspersed with Florida rosemary (*Ceratiola ericoides* Michaux) and maintained by low‐frequency, high‐intensity fires (Menges, [1999](#ece33885-bib-0032){ref-type="ref"}). Scrub is inherently patchy, spatially and temporally, and its historical flux, with attendant opportunities for isolation, is considered to have promoted speciation within this ecosystem. Indeed, Florida scrub is distinguished by high levels of endemism, which includes some 40 species of plants (Christman & Judd, [1990](#ece33885-bib-0007){ref-type="ref"}), four vertebrates (Moler, [1992](#ece33885-bib-0035){ref-type="ref"}; Rodgers, Kale, & Smith, [1996](#ece33885-bib-0045){ref-type="ref"}), and over 50 arthropods (Deyrup, [1989](#ece33885-bib-0010){ref-type="ref"}). Scrub endemics vary widely in overall distribution, with some species being confined locally within a single ridge and others occupying multiple ridges.

The origin of Florida scrub biota has been variously ascribed to Pliocene or Pleistocene epochs (Hubbell, [1961](#ece33885-bib-0020){ref-type="ref"}; Kurz, [1942](#ece33885-bib-0024){ref-type="ref"}; Neill, [1957](#ece33885-bib-0038){ref-type="ref"}), involving a "combination of both historical and edaphic factors" (Huck et al., [1989](#ece33885-bib-0121){ref-type="ref"}). One long‐held hypothesis invoking eastward range expansions of arid‐adapted biota into peninsular Florida (Myers, [1990](#ece33885-bib-0037){ref-type="ref"}) is based on a rich Pliocene fossil record representing numerous extralimital species with western or tropical affinities (Meylan, [1982](#ece33885-bib-0033){ref-type="ref"}; Morgan & Emslie, [2010](#ece33885-bib-0036){ref-type="ref"}). Western taxa purportedly dispersed along a recurrent Gulf Coast corridor---a broad belt of thorn scrub and savanna established by the increased aridity and lowered sea levels associated with Pliocene glaciation (Morgan & Emslie, [2010](#ece33885-bib-0036){ref-type="ref"}). Conversely, interglacial periods reduced corridor width, generated more mesic climatic conditions, and initiated Gulf Coast vicariance with regard to arid‐adapted species. Sundered from their western lineages, Florida populations underwent allopatric speciation in xeric peninsular settings. An alternative scrub origin hypothesis, detailed in Germain‐Aubrey et al. ([2014](#ece33885-bib-0015){ref-type="ref"}), invokes southward dispersal events within eastern North America during Pleistocene glacial advances. Under this scenario, eastern taxa colonized the peninsula and, following local adaptation to xeric ridge environments, experienced ecological/spatial isolation sufficient for speciation (Swenson & Howard, [2005](#ece33885-bib-0054){ref-type="ref"}; Watts, [1975](#ece33885-bib-0056){ref-type="ref"}).

To test the two competing hypotheses, Germain‐Aubrey et al. ([2014](#ece33885-bib-0015){ref-type="ref"}) generated molecular phylogenies for four angiosperm genera, focusing on the topological placement and sister taxon relationships of respective scrub endemics. Their results were ambiguous for one plant, the scrub plum (*Prunus geniculata* Harper), but they identified eastern origins for the remaining three species: Lewton\'s milkwort (*Polygala lewtonii* Small), scrub holly (*Ilex opaca* var. *arenicola* (Ashe) Ashe), and silk bay (*Persea humilis* Nash). However, age estimates for all four angiosperms dated to the Pliocene (or Miocene), a timeframe inconsistent with a Pleistocene colonization predicated by the eastern hypothesis. Additional reports have offered distributional data supporting eastern biogeographical sources for certain scrub taxa (Huck et al., [1989](#ece33885-bib-0121){ref-type="ref"}) and western sources for others (Hubbell, [1961](#ece33885-bib-0020){ref-type="ref"}; Zona & Judd, [1986](#ece33885-bib-0061){ref-type="ref"}), which suggests that eastern versus western origin hypotheses need not be mutually exclusive. With some 90 endemic species yet to be examined, a prevailing colonization pattern for Florida scrub has yet to be vetted.

Here, we examine the biogeography of a scrub arthropod, the Florida sand cockroach, *Arenivaga floridensis* Caudell (family Corydiidae). This species is the eastern representative of a genus otherwise distributed in arid settings from central Texas westward into California and Mexico (Hopkins, [2014a](#ece33885-bib-0017){ref-type="ref"}). As is characteristic of the genus, *A. floridensis* is fossorial and sexually dimorphic (females are wingless; Figure [1](#ece33885-fig-0001){ref-type="fig"}). Inhabiting scrub as well as adjacent sandhills communities, the species requires patches of open sand (Deyrup, [1994](#ece33885-bib-0011){ref-type="ref"}) and demonstrates a strong preference for loose substrate beneath light leaf litter of sand live oaks, *Quercus geminata* (Lamb, Justice, & Justice, [2006](#ece33885-bib-0026){ref-type="ref"}). Males occasionally engage in low, erratic flight at dusk, but juveniles and females appear to be completely fossorial (Deyrup, [1994](#ece33885-bib-0011){ref-type="ref"}). With populations documented from 11 peninsular ridges, *A. floridensis* ranks as Florida scrub\'s most geographically widespread faunal endemic (Lamb, Justice, & Justice, [2006](#ece33885-bib-0026){ref-type="ref"}).

The overall geographic distribution of *Arenivaga* appears consistent with expectations of the western origin hypothesis, underscoring the potential of *A. floridensis* to provide additional biogeographical perspective on scrub colonization and endemism. Support for the western origin hypothesis (with concomitant rejection of the eastern origin hypothesis) would require *A. floridensis* to meet the following predictions: (1) Pliocene colonization of the Florida peninsula, (2) a derived topological placement within a phylogeny of the genus (its basal placement could potentially negate a western origin of *Arenivaga*), and (3) a nodal divergence estimate for *A*. *floridensis* and its sister species approximating a late Miocene--early Pliocene timeframe. To distinguish between hypotheses, we present a detailed intraspecific phylogeography for *A. floridensis* in conjunction with a phylogeny for the species and selected congeners. We compare spatiotemporal patterns of our results with the aforementioned predictions to assess support or refutation of the western origin hypothesis.

2. METHODS {#ece33885-sec-0002}
==========

2.1. Sampling and sequencing regimens {#ece33885-sec-0003}
-------------------------------------

To examine genetic variation in *Arenivaga floridensis,* we pursued a dense intra‐ and inter‐ridge sampling survey that yielded 237 roaches representing 65 localities throughout the species' range ([Appendix](#ece33885-app-0001){ref-type="app"}). Most specimens were captured by sifting sand samples through a two‐tier wire‐mesh (7.0 and 3.0 mm^2^) sieve, which retains all but the smallest nymphs. Roaches were preserved in 95% ethanol, and a rear leg of each was processed for genomic DNA using Qiagen\'s DNeasy kit. We selected the mitochondrial gene cytochrome oxidase II (*cox2*) as our initial genetic marker to assess population divergence and phylogeographical structure, using primers and amplification parameters listed in Table [1](#ece33885-tbl-0001){ref-type="table-wrap"}. Amplicons were cleaned using exoSAP‐IT (USB Corp.) prior to assay on an Applied Biosystems 3130 capillary sequencer. Resulting sequences, edited and assembled in SEQUENCER 4.9 (GeneCodes, Ann Arbor, MI), were aligned in CLUSTALX ver. 2.0 (Larkin et al., [2007](#ece33885-bib-0028){ref-type="ref"}) and translated to ensure correct reading frames.

###### 

Primers and cycling conditions used for DNA amplification

  Gene                                                     Primer                           Sequence (5′--3′)                Anneal. temp. (°C)   Cycles   Reference
  -------------------------------------------------------- -------------------------------- -------------------------------- -------------------- -------- ----------------------------------------------------------------
  *cox1*                                                   C1‐J‐2183                        CAA CAT TTA TTT TGA TTT TTT GG   50                   32       Simon et al. ([1994](#ece33885-bib-0050){ref-type="ref"})
  Roach‐t‐Leucine[a](#ece33885-note-0001){ref-type="fn"}   TCC ATT GCA CTA ATC TGC CA       This study                                                     
  *cox2*                                                   TL2‐J‐3037                       ATG GCA GAT TAG TGC AAT GG       50                   32       Simon et al. ([1994](#ece33885-bib-0050){ref-type="ref"})
  TK‐N‐3785                                                GTT TAA GAG ACC AGT ACT TG                                                                      
  *H3*                                                     HexAF                            ATG GCT ACC AAG CAG ACG GC       61.5                 45       Ogden and Whiting ([2003](#ece33885-bib-0039){ref-type="ref"})
  HexAR                                                    ATA TCC TTG GGC ATG ATG GTG AC                                                                  

A version of Simon et al.\'s ([1994](#ece33885-bib-0050){ref-type="ref"}) TL2‐N‐3014, modified slightly for roaches.

John Wiley & Sons, Ltd

To explore aspects of origin and dispersal, we also generated *cox2* sequences for five additional species of *Arenivaga*. One species, *Arenivaga erratica* Caudell, is the putative sister species to *A. floridensis*, whereas the remaining four, *Arenivaga tonkawa* Hebard, *Arenivaga gumperzae* Hopkins, *Arenivaga gaiophanes* Hopkins, and *Arenivaga bolianna* (Saussure), represent increasingly divergent clades within the genus (Hopkins & Giermakowski, [2014](#ece33885-bib-0019){ref-type="ref"}). We combined *cox2* data for the western species and 32 *A. floridensis* (representing all major and minor ridges in our sample) with sequence data for two additional loci, mitochondrial cytochrome oxidase I (*cox1*) and nuclear histone 3 (*H3*), amplified as detailed in Table [1](#ece33885-tbl-0001){ref-type="table-wrap"}. *Arenivaga erratica* served as the outgroup for the *cox2* dataset; *Eupolyphaga sinensis* (Walker) and *Ergaula capucina* (Brunner von Wattenwyl), representing additional corydiid genera, were outgroup taxa for multilocus analyses.

2.2. Phylogeographical and phylogenetic analysis {#ece33885-sec-0004}
------------------------------------------------

To examine phylogeographical structure in *A. floridensis*, we constructed *cox2* haplotype networks using the TCS algorithm (Clement, Posada, & Crandall, [2000](#ece33885-bib-0009){ref-type="ref"}) implemented in POPART (Leigh & Bryant, [2015](#ece33885-bib-0029){ref-type="ref"}). We also generated phylogenetic networks in SPLITSTREE v4 (Huson & Bryant, [2006](#ece33885-bib-0021){ref-type="ref"}) using the NeighborNet approach with default parameters. We used Bayesian inference (BI) analysis to estimate phylogenetic relationships within *A. floridensis* (*cox2* dataset) and, subsequently, among species (multilocus dataset). Best‐fit codon partitioning schemes and nucleotide substitution models were selected using the Bayesian information criterion (PartitionFinder v.1.1.0; Lanfear, Calcott, Ho, & Guindon, [2012](#ece33885-bib-0027){ref-type="ref"}). We used MRBAYES v.3.2.3 (Ronquist et al., [2012](#ece33885-bib-0046){ref-type="ref"}) to execute two concurrent runs involving four simultaneous chains of 20 million Markov chain Monte Carlo generations, sampling trees every 1,000 generations. The first 25% of the posterior distribution were discarded as burn‐in. Likelihood values for postanalysis trees and parameters were evaluated for convergence using the MRBAYES "sump" command and the program TRACER v. 1.6 (<http://evolve.zoo.ox.ac.uk/software.html?id=tracer>). We used the MRBAYES "sumt" command to generate a majority‐rule consensus tree and calculate posterior probabilities (PP) for consensus nodes.

2.3. Divergence time estimations {#ece33885-sec-0005}
--------------------------------

Divergence time estimates between the Florida lineages and between *A. floridensis* and western roach species were calculated using BEAST 1.8.2 (Drummond, Suchard, Xie, & Rambaut, [2012](#ece33885-bib-0013){ref-type="ref"}). Without fossil calibrations for *Arenivaga*, we used a substitution rate estimated across 15 beetle species for *cox2* (Pons, Ribera, Bertranpetit, & Balke, [2010](#ece33885-bib-0043){ref-type="ref"}) and an average of rates estimated for Hawaiian katydids (Shapiro, Strazanac, & Roderick, [2006](#ece33885-bib-0049){ref-type="ref"}) and beetles (Andújar, Serrano, & Gómez‐Zurita, [2012](#ece33885-bib-0002){ref-type="ref"}; Papadopoulou, Anastasiou, & Vogler, [2010](#ece33885-bib-0042){ref-type="ref"}; Pons et al., [2010](#ece33885-bib-0043){ref-type="ref"}) for *cox1*. Three partitions were analyzed (*cox1*,*cox2, and H3*; each separated by codon position) under PartitionFinder chosen substitution models and a birth--death speciation prior. Runs comprised 50 million generations and were sampled every 1,000 generations. To date the nodes, we used the following rates of substitution: *cox2 *=* *0.02610 substitutions per site per million yrs per lineage (subs/s/My/l) and *cox1 *=* *0.03605 subs/s/My/l. Under each model, the two nonfixed rates were estimated under a LogNormal relaxed clock. Two independent BEAST analyses were conducted and combined using LOGCOMBINER, with the first 10% of each run discarded as burn‐in. Convergence of parameters was accessed using TRACER, and trees were summed using TREEANNOTATOR.

In light of potential issues with concatenated gene data, we also produced a time‐calibrated species tree using StarBEAST2 (Ogilvie, Bouckaert, & Drummond, [2016](#ece33885-bib-0040){ref-type="ref"}) in BEAST2 v2.4.5 (Bouckaert et al., [2014](#ece33885-bib-0004){ref-type="ref"}). Site models, clock models, and topologies were unlinked between the mitochondrial and nuclear loci; ploidy assignments for *cox1* and *cox2* were designated haploid, and *H3* was designated diploid. The population model was set as "Analytical Population Size Integration." Topologies and divergence times were estimated as in BEAST analysis above, using the *cox2* calibration rate.

3. RESULTS {#ece33885-sec-0006}
==========

3.1. Haplotype distribution, phylogeography, and phylogenetics {#ece33885-sec-0007}
--------------------------------------------------------------

We obtained a 684‐bp fragment---effectively, the entire *cox2* gene---for all 237 specimens; sequence comparisons revealed 91 unique haplotypes ([Appendix](#ece33885-app-0001){ref-type="app"}). Haplotype distribution was extremely localized, with 88 of the 91 haplotypes (97%) being observed at single localities. The three haplotypes (*Fp2*,*SB1*,*Vi*) representing more than one locality were limited to proximate sites 2--12 km apart. No haplotypes were shared between ridges despite the fact that certain inter‐ridge populations were within 5 km of each other. Such localized spatial distribution, in part a function of scrub patchiness, may also reflect the species' limited vagility: *Arenivaga floridensis* is restricted to friable, sandy soils (Deyrup, [1994](#ece33885-bib-0011){ref-type="ref"}), and the wingless females may impose additional dispersal constraints.

The BI consensus tree for *cox 2* identified three strongly supported lineages (PP \> 0.99), with evident phylogeographical structure observed within and across lineages (Figures [2](#ece33885-fig-0002){ref-type="fig"} and [3](#ece33885-fig-0003){ref-type="fig"}). The largest lineage, genetically and geographically, is designated the Atlantic‐Central Lineage. Comprising 50 haplotypes, it encompasses four major (Atlantic Coastal, Mount Dora, Northern Brooksville, and Trail) and five minor (Bell, Crescent City, DeLand, Orlando, and Winter Haven) ridges as well as adjacent sandy uplands. This lineage also includes a five‐member haplogroup (*Cl1*,*Cl2*,*LL1*,*LL2*, and *PR*) from the northern end of Lake Wales Ridge. However, the 18 remaining Lake Wales localities, which span the southern 140 km of the ridge\'s 186 km north--south axis, constitute a second distinct lineage of 29 haplotypes, appropriately termed the Lake Wales Lineage. The third lineage is distributed across a single major (Southern Brooksville) and minor ridge (Cotton Plant) and two disjunct scrub parcels (Cedar Key, Tampa Bay). With just 13 haplotypes, this, the Southern Brooksville Lineage, is the smallest but most genetically distinct lineage, constituting the sister group to the Atlantic‐Central/Lake Wales clade. These three major lineages were also identified in both SPLITSTREE (Figure [3](#ece33885-fig-0003){ref-type="fig"}) and TCS haplotype networks (not illustrated), the latter providing finer details of phylogeographic structure less pertinent to the predictions being tested.

![BI consensus tree from *cox2* dataset illustrating the three major lineages of *Arenivaga floridensis*, coded as: orange = Atlantic‐Central Lineage, blue = Lake Wales Lineage, and purple = Southern Brooksville Lineage (Figure [S1](#ece33885-sup-0001){ref-type="supplementary-material"} depicts *cox2* tree with labeled terminals and complete PP listings)](ECE3-8-5254-g002){#ece33885-fig-0002}

![SPLITSTREE network and distribution map, illustrating the three major lineages of *Arenivaga floridensis*, coded as: orange = Atlantic‐Central, blue = Lake Wales, and purple = Southern Brooksville (see [Appendix](#ece33885-app-0001){ref-type="app"} for haplotype abbreviations and localities). Inset map depicts peninsular uplands exposed during the Pliocene\'s last glacial minimum (3.2--2.7 Ma)](ECE3-8-5254-g003){#ece33885-fig-0003}

The three lineages were again recovered (all PP = 1.0) in the concatenated BI analysis but with higher nodal support for the Atlantic‐Central/Lake Wales clade (PP = 0.98; Figure [4](#ece33885-fig-0004){ref-type="fig"}) than in the *cox2* tree (PP = 0.79). The concatenated tree also reveals strongly supported interspecific relationships (all PP = 1.0), which are topologically congruent with Hopkins and Giermakowski\'s ([2014](#ece33885-bib-0019){ref-type="ref"}) molecular--morphological phylogeny for the genus. Namely, *Arenivaga floridensis* is recovered as the sister taxon to an *A. erratica* + *A. tonkawa* clade, whereas the species *A. gumperzae*,*A. gaiophanes*, and *A. bolianna* represent increasingly divergent clades.

![BI consensus tree from concatenated dataset for a subset of *Arenivaga floridensis* and selected western congeners. The three major lineages of *Arenivaga floridensis* are coded as: orange = Atlantic‐Central, blue = Lake Wales, and purple = Southern Brooksville. Circles indicate PP = 1.0; PP values for within‐lineage nodes are not shown, and outgroup terminals are not labeled. (Figure [S2](#ece33885-sup-0001){ref-type="supplementary-material"} depicts concatenated tree with labeled terminals and complete PP listings)](ECE3-8-5254-g004){#ece33885-fig-0004}

3.2. Divergence dating {#ece33885-sec-0008}
----------------------

*Cox2* divergence time estimates within *A. floridensis* and among species are depicted in Figure [5](#ece33885-fig-0005){ref-type="fig"}. The nodal age between the Lake Wales and Atlantic--Central lineages is estimated to be 2.27 Ma (95% highest posterior density \[HPD\]: 2.92--1.71), whereas coalescence of all three roach lineages dates to 2.78 Ma (HPD: 3.61--2.08). The nodal age estimated for the *floridensis*--*erratica*+*tonkawa* clade is 6.36 Ma (HPD: 8.39--4.69). The species tree generated in StarBEAST2 yielded an identical topology (Figure [S3](#ece33885-sup-0001){ref-type="supplementary-material"}), though with weaker nodal support (a function of limited variation in *H3*) and a slightly younger age estimate for *floridensis*--*erratica*+*tonkawa* clade (*cox2 *=* *4.02 Ma; HPD: 5.40--2.26). In turn, *cox1* divergence time estimates are older (BEAST consensus tree; Figure [S4](#ece33885-sup-0001){ref-type="supplementary-material"}), with the Lake Wales and Atlantic--Central clade dating to 3.17 Ma (HPD: 3.98--2.47), lineage coalescence for *A*. *floridensis* to 3.88 Ma (HPD: 4.89--2.99) and the *floridensis*--*erratica*+*tonkawa* clade to 8.84 Ma (HPD: 11.33--6.61). Although the *Arenivaga* age estimates should be interpreted with caution in the absence of fossil calibration, we note that *cox2* and *cox1* HPD intervals exhibit substantial overlap for each of the three nodes. Furthermore, both genes offer general timeframes for peninsular colonization (late Miocene‐Pliocene) and ridge population isolation (Pliocene) that clearly predate the Pleistocene.

![BEAST chronogram for the major lineages of *Arenivaga floridensis* and western *Arenivaga* for the concatenated dataset using *cox2* calibration. Blue HPD divergence bars subtend selected nodes, and the gray panel delimits the Pliocene (5.333--2.58 Ma). With exception of two nodes within lineages, support throughout the tree is PP \>0.98; all HPD nodes have PP = 1. Lineages of *A. floridensis* are as follows: orange bar = Atlantic‐Central, blue = Lake Wales, and purple = Southern Brooksville](ECE3-8-5254-g005){#ece33885-fig-0005}

4. DISCUSSION {#ece33885-sec-0009}
=============

4.1. Phylogeographical structure and lineage ages {#ece33885-sec-0010}
-------------------------------------------------

Our genetic survey of *Arenivaga floridensis* reveals a scrub endemic characterized by pronounced phylogeographical structure, much of it partitioned by sand ridges. Genetic differentiation among the three major lineages is likely attributable to isolation imposed by episodic marine transgressions, which inundated substantial portions of Florida during the Pliocene and Pleistocene epochs. Estimated lineage ages correspond closely with peninsular fragmentation during the Pliocene warm period (3.2--2.7 Ma) when sea levels rose significantly, with a peak some 22 m higher than present day (Miller et al., [2012](#ece33885-bib-0034){ref-type="ref"}). The specifics of the Florida platform\'s inundation require incorporating concomitant geologic factors such as isostatic uplift (Adams, Opdyke, & Jaeger, [2010](#ece33885-bib-0001){ref-type="ref"}; Rowley et al., [2013](#ece33885-bib-0047){ref-type="ref"}). Nonetheless, peninsular uplands emergent at this time were likely limited to portions of the Lake Wales, Mount Dora, Trail, and Southern Brooksville ridges, which offer striking geographical congruence with the Lake Wales, Central‐Atlantic, and Southern Brooksville lineages (Figure [3](#ece33885-fig-0003){ref-type="fig"}). The older age observed for the Southern Brooksville Lineage suggests a longer isolation history, which is consistent with the presence of specific Southern Brooksville endemics (Deyrup, [2005](#ece33885-bib-0012){ref-type="ref"}; Neill, [1957](#ece33885-bib-0038){ref-type="ref"}; Squitier, Deyrup, & Capinera, [1998](#ece33885-bib-0053){ref-type="ref"}) and notable absence of more widespread scrub species, for example, the Florida sand skink, *Plestiodon reynoldsi* (Stejneger), and Florida scrub lizard, *Sceloporus woodi* Stejneger. Uniquely, the Lake Wales Ridge is inhabited by populations representing two major lineages: the Lake Wales Lineage, present from the ridge\'s central portion to its southern terminus, and the Atlantic‐Central Lineage, represented by a 5‐member haplogroup confined to the ridge\'s northernmost sector. The presence of Atlantic‐Central haplotypes on the Lake Wales Ridge could reflect dispersal from the adjacent Mount Dora Ridge. However, phylogeographical structuring in a second scrub endemic, the aforementioned Florida sand skink, provides a more compelling interpretation. Substantive mtDNA divergence (*cob*, 4.5%) between skink populations from northern Lake Wales and those further south (Richmond, Reid, Ashton, & Zamudio, [2009](#ece33885-bib-0044){ref-type="ref"}) exhibits precise geographical congruence and a similar level of divergence (*cox2*, 5.5%) with roach haplotypes. Their genealogical concordance suggests a shared history whereby codistributed ancestral populations experienced vicariant separation in the northern Lake Wales Ridge. Overall, the phylogeography of *A. floridensis* reveals the combined influences of relatively early colonization (2.78 ma), close edaphic ties to peninsular ridges, and an episodic flux of scrub habitat under the climatic oscillations of Plio‐Pleistocene glaciation. Divergence levels detected for the three major lineages clearly support prediction 1 of the western origin hypothesis: *Arenivaga floridensis* has been present on the Florida peninsula since the Pliocene. Moreover, the roach\'s fidelity to Florida scrub and persistence through rounds of ecosystem expansion/fragmentation are aptly reflected in its phylogeographical complexity.

4.2. A western origin for *Arenivaga floridensis* {#ece33885-sec-0011}
-------------------------------------------------

*Arenivaga* comprises a desert‐dwelling clade of roaches distributed largely within southwestern North America. The genus is now recognized as being far more species rich (48 vs. 9 spp) than traditionally perceived, with its highest levels of diversity centered in northern Mexico (Hopkins, [2014a](#ece33885-bib-0017){ref-type="ref"}). Although *A. floridensis* is considered to be the only species native to the southeastern United States, we note that two male specimens of *A. bolliana* have been collected in Florida (from the cities of Gainesville and Seminole). Otherwise, the eastern range boundary of *A. bolliana* terminates in Texas (some 1,400 km west); if viable populations do occur in Florida, they were probably recently introduced and remain extremely localized. Hopkins and Giermakowski\'s ([2014](#ece33885-bib-0019){ref-type="ref"}) phylogeny of *Arenivaga* (24 of the 48 species) identified *A. erratica* as the sister species to *A. floridensis*, providing consistency with an earlier morphological interpretation that placed *A. floridensis* in the "*erratica* group" (Hebard, [1920](#ece33885-bib-0016){ref-type="ref"}). Our phylogeny differs slightly, recovering a sister group relationship between *A. floridensis* and *A. erratica* + *A. tonkawa*. This topological discrepancy probably reflects the influence of morphological characters in Hopkins and Giermakowski\'s ([2014](#ece33885-bib-0019){ref-type="ref"}) integrated dataset. Hopkins' ([2014b](#ece33885-bib-0018){ref-type="ref"}) BI consensus trees for DNA sequence data, generated separately for *cox1*,*12S,* and *H3* genes, depict slightly different relationships among these three species, although all are invariably recovered in the same clade. Topological placement of *A. floridensis*, relative to other congeners, reveals a derived position within the genus, both in our BI tree and in Hopkins and Giermakowski\'s ([2014](#ece33885-bib-0019){ref-type="ref"}) phylogeny. As well, their phylogeny depicts a Baja California endemic*, Arenivaga diaphana* Hopkins, as the sister taxon to the remaining 23 species surveyed, which corroborates southwestern North America as the center of origin for *Arenivaga*. Our estimated dates for the most recent common ancestor of *A. floridensis* and *A. erratica* + *A. tonkawa* (6.36 Ma) provide a timeframe spanning the late Miocene to early Pliocene epochs. Arid conditions linking Florida and southwestern North America at this time facilitated biotic interchange and account for a strong western influence on the peninsula\'s Plio‐Pleistocene faunas (Krysko, Nuñez, Lippi, Smith, & Granatosky, [2016](#ece33885-bib-0023){ref-type="ref"}; Morgan & Emslie, [2010](#ece33885-bib-0036){ref-type="ref"}). Complementary patterns in diversity, distribution, and phylogenetic relationships for *Arenivaga*, together with nodal estimates for the *erratica*‐*floridensis*‐*tonkawa* clade*,* provide spatiotemporal consistency with an eastward migration from southwestern North America. These results satisfy predictions 2 (derived topological placement) and 3 (appropriate sister taxon divergence estimate) of the western origin hypothesis.

4.3. Reevaluating scrub origin hypotheses {#ece33885-sec-0012}
-----------------------------------------

Our phylogeographical and phylogenetic results for *Arenivaga* support a different scrub origin hypothesis than Germain‐Aubrey et al.\'s ([2014](#ece33885-bib-0015){ref-type="ref"}) angiosperm phylogenies but, importantly, indicate that the two hypotheses are not mutually exclusive. Although no other studies have actually tested either hypothesis, it is possible to infer eastern versus western origins for additional scrub endemics (or near endemics) from published phylogenies. For example, the evolutionary history of the Florida scrub‐jay exhibits remarkable parallels with that of *Arenivaga floridensis*. The Florida scrub‐jay---also the sole eastern representative of its genus---is sister taxon to four western species that form the 'scrub‐jay group' within *Aphelocoma* jays (McCormack et al., [2011](#ece33885-bib-0031){ref-type="ref"}). Sequence divergence (mitochondrial *cob* and *nad2*) between the Florida and western scrub‐jays, at 8.0%, with a nodal estimate of \~6.0 Ma (McCormack et al., [2011](#ece33885-bib-0031){ref-type="ref"}), corresponds with our observed *cox2* divergence and age estimate between *Arenivaga floridensis* and *A. erratica* + *A. tonkawa* (8.6%, 6.36 Ma). This case, together with phylogenies for additional taxa (Table [2](#ece33885-tbl-0002){ref-type="table-wrap"}), reveal unequivocal eastern and western contributions to Florida scrub biota. The eastern origin hypothesis advocates a Pleistocene inception of Florida scrub, invoking dispersal of eastern North American species onto the peninsula in response to glacial advances between 2.56 Ma and 10 ka. Compiling Germain‐Aubrey et al.\'s ([2014](#ece33885-bib-0015){ref-type="ref"}) phylogenetic results with those of others, we identified eight scrub endemics with sister species from eastern North America (Table [2](#ece33885-tbl-0002){ref-type="table-wrap"}). However, only two cases appear indicative of Pleistocene speciation. The first involves silk bay, *Persea humilis*, and closely related taxa, *Persea borbonia* (L.) Sprengel and *Persea palustris* (Rafinesque) Sargent, which form a clade dating to near the Pliocene--Pleistocene interface (Germain‐Aubrey et al., [2014](#ece33885-bib-0015){ref-type="ref"}). The second example involves the grasshopper *Schistocerca ceratiola* Hubbell and Walker, an obligate, sedentary specialist on Florida rosemary that exhibits minimal genetic distance across ridges (0.16%--0.3%; Lamb & Justice, [2005](#ece33885-bib-0025){ref-type="ref"}) and recent divergence (\~0.6 Ma) from its sister species (Song, Moulton, Hiatt, & Whiting, [2013](#ece33885-bib-0052){ref-type="ref"}). If the western origin cases are also considered, then age estimates for nine of ten scrub endemics date to peninsular residence times of Pliocene age or older. Although Pleistocene glaciation unquestionably influenced intraspecific phylogeography in eastern North America (Avise, [2000](#ece33885-bib-0003){ref-type="ref"}; Soltis et al., [2006](#ece33885-bib-0150){ref-type="ref"}), its role in scrub speciation appears to have been minimal. Germain‐Aubrey et al. ([2014](#ece33885-bib-0015){ref-type="ref"}) noted a Miocene cooling trend in eastern North American (Foster, Lunt, & Parrish, [2010](#ece33885-bib-0014){ref-type="ref"}) that may have elicited earlier southward migrations, which would provide greater temporal consistency with observed genetic differentiation between scrub endemics and their sister taxa. We have provided the first molecular phylogenetic evidence used to test and confirm the western origin hypothesis for a Florida scrub endemic and have compiled additional molecular phylogenies identifying both eastern and western progenitors to scrub species. These cases reveal more examples of eastern origin than western, and proximity alone would seem to favor eastern contribution potential over intermittent dispersal from the southwest. Still, we hesitate to endorse a prevailing eastern assembly when the evolutionary histories for most scrub species have yet to be determined. It is clear now, however, that arguments for a single regional source leading to the inception of this ecosystem can be rightly dismissed; rather, a composite biogeographical history involving disparate biotas better explains Florida scrub origins and endemism.

###### 

Geographic origins and epoch age assignments for selected Florida scrub endemics

  Taxon                                             Hypothesis support   Sister taxon range[a](#ece33885-note-0002){ref-type="fn"}   Epoch age assignment     Reference
  ------------------------------------------------- -------------------- ----------------------------------------------------------- ------------------------ -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  Plants                                                                                                                                                      
  Florida rosemary (*Ceratiola ericoides)*          Eastern              NE                                                          Miocene                  Li, Alexander, Ward, Del Tredici, and Nicholson ([2002](#ece33885-bib-0030){ref-type="ref"}) and Trapnell, Schmidt, Quintana‐Ascencio, and Hamrick ([2007](#ece33885-bib-0055){ref-type="ref"})
  Scrub holly (*Ilex opaca* var. *arenicola*)       Eastern              E                                                           Pliocene                 Germain‐Aubrey et al. ([2014](#ece33885-bib-0015){ref-type="ref"})
  Lewton\'s polygala (*Polygala lewtonii*)          Eastern              E                                                           Miocene                  Germain‐Aubrey et al. ([2014](#ece33885-bib-0015){ref-type="ref"})
  Silk bay (*Persea humilis*)                       Eastern              SE                                                          Pleistocene              Germain‐Aubrey et al. ([2014](#ece33885-bib-0015){ref-type="ref"})
  Scrub plum (*Prunus geniculata*)                  Ambiguous            NA                                                          Mio‐Pliocene interface   Germain‐Aubrey et al. ([2014](#ece33885-bib-0015){ref-type="ref"})
  Scrub hickory (*Carya floridana*)                 Eastern              SE                                                          Pliocene                 Zhang et al. ([2013](#ece33885-bib-0060){ref-type="ref"})
  Florida jujube (*Condalia celata*)                Western              SW                                                          \_                       Islam and Guralnick ([2015](#ece33885-bib-0022){ref-type="ref"})
  Animals                                                                                                                                                     
  Florida scrub‐jay (*Aphelocoma coerulescens*)     Western              SW                                                          Mio‐Pliocene interface   McCormack, Heled, Delaney, Peterson, and Knowles ([2011](#ece33885-bib-0031){ref-type="ref"})
  Florida sand skink (*Plestiodon reynoldsi*)       Eastern              SE                                                          Miocene                  Brandley et al. ([2011](#ece33885-bib-0006){ref-type="ref"})
  Florida scrub lizard (*Sceloporus woodi*)         Eastern              SE                                                          \_                       Wiens, Kuczynski, Arif, and Reeder ([2010](#ece33885-bib-0059){ref-type="ref"})
  Rosemary grasshopper (*Schistocerca ceratiola*)   Eastern              SE                                                          Pleistocene              Song ([2004](#ece33885-bib-0051){ref-type="ref"}) and Song et al. ([2013](#ece33885-bib-0052){ref-type="ref"})

North American sectors abbreviated as: E = eastern, NA = North America, NE = northeastern, SE = southeastern, SW = southwestern.

John Wiley & Sons, Ltd

CONFLICT OF INTEREST {#ece33885-sec-0014}
====================

None declared.

AUTHOR CONTRIBUTIONS {#ece33885-sec-0015}
====================

TL conceived the project; TJ, PM, TL, and JB collected specimens; TJ and TL generated/analyzed sequence data; MB and JB conducted network and phylogenetic analyses; HH provided material and sequence data for the western *Arenivaga*; TL wrote the manuscript; MB generated figures and contributed to the text.

DATA ACCESSIBILITY {#ece33885-sec-0016}
==================

DNA sequences: GenBank accession numbers---[DQ873966.1](DQ873966.1), [DQ874280.1](DQ874280.1), [FJ830540.1](FJ830540.1), [KP986405.1](KP986405.1), [MG761831](MG761831)--[MG761994](MG761994).

Supporting information
======================

###### 

 

###### 

Click here for additional data file.

We thank Mark Deyrup, Michael Evans, Michael Justice, Jackson Mosely, and Paul Skelley for their assistance in the field. Kenny Krysko graciously provided the inset map in Figure [2](#ece33885-fig-0002){ref-type="fig"}. Research was supported by contract NG01‐002 with the Florida Fish & Wildlife Conservation Commission to TL and TJ as well as by the American Museum of Natural History\'s Theodore Roosevelt Fund and an Archbold Biological Station internship to TJ.

Localities, coordinates, sample sizes, and *cox2* haplotype designations (derived from locality abbreviations) for populations of *Arenivaga floridensis*, partitioned by ridges, from north to south {#ece33885-sec-0017}
=====================================================================================================================================================================================================

RidgeCoordinates*n*Haplotype (n)Trail83Gold Head State Park29.82782, −81.949904*GH* (4)Ordway Preserve29.72840, −81.979902*Od* (2)Interlachen29.6243, −81.91982*It* (2)Atlantic Coastal3510Roseland27.82842, −80.478142*Ro1* (1), *Ro2* (1)Viking27.54285, −80.362953*Vi (3)*Ft. Pierce27.47035, −80.333815*Vi* (5)Savannas State Preserve27.29906, −80.258374*Sa* (4)Hobe Sound A27.10840, −80.168273*HS1* (3)Hobe Sound B27.05972, −80.140455*HS2*(4), *HS3* (1)Jonathan Dickinson State Park27.01693, −80.110153*JD* (3)Jupiter Ridge Natural Area26.91608, −80.073054*JR* (4)Juno Beach26.87847, −80.055406*JB* (6)Bell 42Bell29.79037, −82.85362 *Be1*(3), *Be2* (1)Northern Brooksville 156Archer29.50267, −82.572505*Ar1* (2), *Ar2* (3)Newberry29.6313, −82.70505*Ne* (5)Northern Brooksville Ridge (south)29.11872, −82.447271*NB*Williston29.3314, −82.53674*Wl1* (3), *Wl2* (1)Crescent City 115Pomona Park29.49041, −81.584345*PP1* (4), *PP2* (1)Como Lake A29.46995, −81.604615*Co1* (3), *Co2* (2)Como Lake B29.47275, −81.590631*Co3* (1)DeLand 41DeLand29.01872, −81.23567 *De* (3)Cotton Plant 63Morriston A29.1936, −82.40124*Mo1* (3), *Mo2* (1)Morriston B29.1891, −82.38832*Mo3* (2)Southern Brooksville 156Dunnellon29.03074, −82.531883*Du1* (2), *Du2* (1)Weeki Wachee28.52517, −82.569952*WW* (2)Weeki Wachee Springs28.51512, −82.573183*WS* (3)Withlacoochee State Forest28.79871, −82.381197*Wi1* (6), *Wi2* (1)Mount Dora 215Ocala Natl. Forest, Lake Kerr29.36503, −81.822044*Oc* (4)Ocala Natl. Forest, firetower29.17800, −81.777306*Oc* (6)Ocala Natl. Forest, Astor Park29.12346, −81.577892*AP* (2)Ocala Natl. Forest, Tomahawk Lake29.13381, −81.891075*TL* (5)Wekiwa Springs State Park28.71048, −81.466534*We1* (3), *We2* (1)Lake Apopka28.54625, −81.701002*LA* (2)Orlando 31Orlando28.54230, −81.31950 *Or* (3)Lake Wales 8335Clermont28.55583, −81.742333*Cl1* (2), *Cl2* (1)Palatlakaha River Park28.50334, −81.749501*PR* (1)Lake Louisa28.47346, −81.716005*LL1* (3), *LL2* (2)Lake Marion28.07444, −81.546594*LM* (4)Dundee28.02648, −81.634082*Dd* (2)Frostproof27.70455, −81.561203*Fp1* (2), *Fp2* (1)Snell Creek28.13136, −81.541613*SC* (3)Lake Wales27.92406, −81.605745*LW1* (3), *LW2* (1), *LW3* (1)Warner College27.70455, −81.561205*WC* (5)Hickory Lake Scrub Preserve27.69647, −81.539643*HL* (1), *Fp2* (1)Lake Streety27.69037, −81.563924*LS* (4)Saddle Blanket Scrub Preserve27.66950, −81.576334*SB1* (2), *SB2* (2)Silver Lake27.56431, −81.523155*SB1* (1), *SL1* (2), *SL2* (1), *SL3* (1)Carter Creek Scrub Preserve27.54008, −81.408604*CC* (4)Highlands Hammock State Park27.46728, −81.520694*HH* (4)Sebring27.44177, −81.418636*Se1* (4), *Se2* (2)Josephine Creek27.36860, −81.401734*JC* (4)Lake Placid27.24863, −81.300635*LP*1 (1), *LP2* (3), *LP3* (1)Placid Lakes27.24169, −81.401422*PL* (2)Archbold Biological Station27.18185, −81.352178*ABS1*(5), *ABS2* (1), *ABS3* (2)Venus27.10675, −81.332401*Ve* (1)Winter Haven 85Winter Haven28.04417, −81.727174*WH1* (2), *WH2* (2)Bartow27.96510, −81.773402*Ba* (2)McLeod Lake27.96309, −81.743142*Mc1* (1), *Mc2* (1)Unnamed ridges/uplandsDixie29.60855, −83.070155*Di1* (3), *Di2* (1), *Wl2* (1)Indian Lake State Forest29.2835, −82.05391*IL*Cedar Key29.18520, −83.019014*CK* (4)Silver Springs Shores A29.06142, −81.908235*SS1* (3) *SS2* (2)Silver Springs Shores B29.0415, −81.87102*SS3* (2)Silver Springs Shores C29.0558, −81.88461*SS4* (1)Tampa28.0698, −82.39486*Ta* (6)
