CALIBRATING DIVERGENCE TIMES ON SPECIES TREES VERSUS GENE TREES: IMPLICATIONS FOR SPECIATION HISTORY OF APHELOCOMA JAYS

Mitogenomics clarifies the position of the Nearctic magpies (Pica hudsonia and Pica nuttalli) within the Holarctic magpie radiation

Partial separation of a peripheral population may lead to its divergence and, potentially, speciation due to genetic drift followed by selection and geographic isolation. This process may cause taxonomic uncertainty because reproductive isolation in allopatry cannot be verified directly. The two Nearctic allopatric species of magpies (Aves, Corvidae: Pica) serve as a good example of these problems. The Black-billed magpie Pica hudsonia is widely distributed in North America, whereas the Yellow-billed Magpie Pica nuttalli is endemic to a restricted range in California. Their relationships with Palearctic species have been little studied. We obtained complete mitochondrial genomes of both Nearctic magpie species, along with the Eurasian Magpie (Pica pica) and the Oriental Magpie (Pica serica), 20 mitogenomes in total. Phylogenetic analysis reveals a basal position of P. serica, and P. pica as a sister clade to the two Nearctic species. P. hudsonia and P. nuttalli form reciprocal monophyletic subclades, showing recent divergence between and within them. Our data show that the Nearctic magpie lineage diverged from the common ancestor with P. pica, with a single migration wave via the Beringia. Within the Nearctic, we hypothesize a peripatric mode of speciation among Pica taxa due to the divergence and separation of the small marginal population in California below the Sierra-Nevada mountains. Diversifying amino acid substitutions in ND4-ND5-ND6 genes along the branch leading to the New World clade may indicate selection for heat-tolerance. Considering the clear phenotypic differences between P. hudsonia and P. nuttalli, our data, showing their reciprocal monophylies and genetic distinctness, is consistent with the two-species taxonomy.

Museomics help resolving the phylogeny of snowfinches (Aves, Passeridae, Montifringilla and allies)

Historical specimens from museum collections provide a valuable source of material also from remote areas or regions of conflict that are not easily accessible to scientists today. With this study, we are providing a taxon-complete phylogeny of snowfinches using historical DNA from whole skins of an endemic species from Afghanistan, the Afghan snowfinch, Pyrgilauda theresae. To resolve the strong conflict between previous phylogenetic hypotheses, we generated novel mitogenome sequences for selected taxa and genome-wide SNP data using ddRAD sequencing for all extant snowfinch species endemic to the Qinghai-Tibet Plateau (QTP) and for an extended intraspecific sampling of the sole Central and Western Palearctic snowfinch species (Montifringilla nivalis). Our phylogenetic reconstructions unanimously refuted the previously suggested paraphyly of genus Pyrgilauda. Misplacement of one species-level taxon (Onychostruthus tazcanowskii) in previous snowfinch phylogenies was undoubtedly inferred from chimeric mitogenomes that included heterospecific sequence information. Furthermore, comparison of novel and previously generated sequence data showed that the presumed sister-group relationship between M. nivalis and the QTP endemic M. henrici was suggested based on flawed taxonomy. Our phylogenetic reconstructions based on genome-wide SNP data and on mitogenomes were largely congruent and supported reciprocal monophyly of genera Montifringilla and Pyrgilauda with monotypic Onychostruthus being sister to the latter. The Afghan endemic P. theresae likely originated from a rather ancient Pliocene out-of-Tibet dispersal probably from a common ancestor with P. ruficollis. Our extended trans-Palearctic sampling for the white-winged snowfinch, M. nivalis, confirmed strong lineage divergence between an Asian and a European clade dated to 1.5 - 2.7 million years ago (mya). Genome-wide SNP data suggested subtle divergence among European samples from the Alps and from the Cantabrian mountains.

Molecular phylogeny and morphology of Sporodiniella sinensis sp. nov. (Syzygitaceae, Mucorales), an invertebrate-associated species from Yunnan, China

During investigations of invertebrate-associated fungi in Yunnan Province of China, a new species, Sporodiniella sinensis sp. nov., was collected. Morphologically, S. sinensis is similar to Sporodiniella umbellata; however, it is distinguished from S. umbellata by its greater number of sporangiophore branches, longer sporangiophores, larger sporangiospores, and columellae. The novel species exhibits similarities of 91.62 % for internal transcribed spacer (ITS), 98.66-99.10 % for ribosomal small subunit (nrSSU), and 96.36-98.22 % for ribosomal large subunit (nrLSU) sequences, respectively, compared to S. umbellata. Furthermore, phylogenetic analyses based on combined sequences of ITS, nrLSU and nrSSU show that it forms a separate clade in Sporodiniella, and clusters closely with S. umbellata with high statistical support. The phylogenetic and morphological evidence support S. sinensis as a distinct species. Here, it is formally described and illustrated, and compared with other relatives.

Unfamiliarity generates costly aggression in interspecific avian dominance hierarchies

Dominance hierarchies often form between species, especially at common feeding locations. Yet, relative to work focused on the factors that maintain stable dominance hierarchies within species, large-scale analyses of interspecific dominance hierarchies have been comparatively rare. Given that interspecific behavioral interference mediates access to resources, these dominance hierarchies likely play an important and understudied role in community assembly and behavioral evolution. To test alternative hypotheses about the formation and maintenance of interspecific dominance hierarchies, we employ an large, participatory science generated dataset of displacements observed at feeders in North America in the non-breeding season. Consistent with the hypothesis that agonistic interference can be an adaptive response to exploitative competition, we find that species with similar niches are more likely to engage in costly aggression over resources. Among interacting species, we find broad support for the hypothesis that familiarity (measured as fine-scale habitat overlap) predicts adherence to the structure of the dominance hierarchy and reduces aggression between species. Our findings suggest that the previously documented agonistic hierarchy in North American birds emerges from species-level adaptations and learned behaviors that result in the avoidance of costly aggression.

De novo assembly of a chromosome-level reference genome for the California Scrub-Jay, Aphelocoma californica

We announce the assembly of the first de novo reference genome for the California Scrub-Jay (Aphelocoma californica). The genus Aphelocoma comprises four currently recognized species including many locally adapted populations across Mesoamerica and North America. Intensive study of Aphelocoma has revealed novel insights into the evolutionary mechanisms driving diversification in natural systems. Additional insights into the evolutionary history of this group will require continued development of high-quality, publicly available genomic resources. We extracted high molecular weight genomic DNA from a female California Scrub-Jay from northern California and generated PacBio HiFi long-read data and Omni-C chromatin conformation capture data. We used these data to generate a de novo partially phased diploid genome assembly, consisting of two pseudo-haplotypes, and scaffolded them using inferred physical proximity information from the Omni-C data. The more complete pseudo-haplotype assembly (arbitrarily designated "Haplotype 1") is 1.35 Gb in total length, highly contiguous (contig N50 = 11.53 Mb), and highly complete (BUSCO completeness score = 97%), with comparable scaffold sizes to chromosome-level avian reference genomes (scaffold N50 = 66.14 Mb). Our California Scrub-Jay assembly is highly syntenic with the New Caledonian Crow reference genome despite ~10 million years of divergence, highlighting the temporal stability of the avian genome. This high-quality reference genome represents a leap forward in publicly available genomic resources for Aphelocoma, and the family Corvidae more broadly. Future work using Aphelocoma as a model for understanding the evolutionary forces generating and maintaining biodiversity across phylogenetic scales can now benefit from a highly contiguous, in-group reference genome.

Crossing the Rift valley: using complete mitogenomes to infer the diversification and biogeographic history of ethiopian highlands Ptychadena (anura: Ptychadenidae)

The Ethiopian Highlands are considered a biodiversity hotspot, harboring a high number of endemic species. Some of the endemic species probably diversified in situ; this is, for example, the case of a monophyletic clade containing 12 known species of grass frogs of the genus Ptychadena. The different species occur at elevations ranging from 1,500 to above 3,400 m and constitute excellent models to study the process of diversification in the highlands as well as adaptations to high elevations. In this study, we sampled 294 specimens across the distribution of this clade and used complete mitogenomes and genome-wide SNP data to better understand how landscape features influenced the population structure and dispersal of these grass frogs across time and space. Using phylogenetic inference, population structure analyses, and biogeographic reconstructions, we found that the species complex probably first diversified on the south-east side of the Great Rift Valley. Later on, species dispersed to the north-west side, where more recent diversification occurred. We further demonstrate that Ptychadena species have dispersed across the Great Rift Valley at different times. Our analyses allowed for a more complete understanding of the contribution of geological events, biogeographic barriers and climatic changes as drivers of species diversification and adaptation in this important biogeographic region.

Revisiting phylogeny, systematics, and biogeography of a Pleistocene radiation

PREMISE

Pilosocereus (Cactaceae) is an important dry forest element in all subregions and transitional zones of the neotropics, with the highest diversity in eastern Brazil. The genus is subdivided into informal taxonomic groups; however, most of these are not supported by recent molecular phylogenetic inferences. This lack of confidence is probably due to the use of an insufficient number of loci and the complexity of cactus diversification. Here, we explored the species relationships in Pilosocereus in more detail, integrating multilocus phylogenetic approaches with the assessment of the ancestral range and the effect of geography on diversification shifts.

METHODS

We used 28 nuclear, plastid, and mitochondrial loci from 54 plant samples of 31 Pilosocereus species for phylogenetic analyses. We used concatenated and coalescent phylogenetic trees and Bayesian models to estimate the most likely ancestral range and diversification shifts.

RESULTS

All Pilosocereus species were clustered in the same branch, except P. bohlei. The phylogenetic relationships were more associated with the geographic distribution than taxonomic affinities among taxa. The genus began diversifying during the Plio-Pleistocene transition in the Caatinga domain and experienced an increased diversification rate during the Calabrian age.

CONCLUSIONS

We recovered a well-supported multispecies coalescent phylogeny. Our results refine the pattern of rapid diversification of Pilosocereus species across neotropical drylands during the Pleistocene and highlight the need for taxonomic rearrangements in the genus. We recovered a pulse of diversification during the Pleistocene that was likely driven by multiple dispersal and vicariance events within and among the Caatinga, Cerrado, and Atlantic Forest domains.

The Implications of Incongruence between Gene Tree and Species Tree Topologies for Divergence Time Estimation

Phylogenetic analyses are increasingly being performed with data sets that incorporate hundreds of loci. Due to incomplete lineage sorting, hybridization, and horizontal gene transfer, the gene trees for these loci may often have topologies that differ from each other and from the species tree. The effect of these topological incongruences on divergence time estimation has not been fully investigated. Using a series of simulation experiments and empirical analyses, we demonstrate that when topological incongruence between gene trees and the species tree is not accounted for, the temporal duration of branches in regions of the species tree that are affected by incongruence is underestimated, whilst the duration of other branches is considerably overestimated. This effect becomes more pronounced with higher levels of topological incongruence. We show that this pattern results from the erroneous estimation of the number of substitutions along branches in the species tree, although the effect is modulated by the assumptions inherent to divergence time estimation, such as those relating to the fossil record or among-branch-substitution-rate variation. By only analyzing loci with gene trees that are topologically congruent with the species tree, or only taking into account the branches from each gene tree that are topologically congruent with the species tree, we demonstrate that the effects of topological incongruence can be ameliorated. Nonetheless, even when topologically congruent gene trees or topologically congruent branches are selected, error in divergence time estimates remains. This stems from temporal incongruences between divergence times in species trees and divergence times in gene trees, and more importantly, the difficulty of incorporating necessary assumptions for divergence time estimation. [Divergence time estimation; gene trees; species tree; topological incongruence.].

Environmental correlates of phenotypic evolution in ecologically diverse Liolaemus lizards

Evolutionary correlations between phenotypic and environmental traits characterize adaptive radiations. However, the lizard genus Liolaemus, one of the most ecologically diverse terrestrial vertebrate radiations on earth, has so far shown limited or mixed evidence of adaptive diversification in phenotype. Restricted use of comprehensive environmental data, incomplete taxonomic representation and not considering phylogenetic uncertainty may have led to contradictory evidence. We compiled a 26-taxon dataset for the Liolaemus gracilis species group, representing much of the ecological diversity represented within Liolaemus and used environmental data to characterize how environments occupied by species' relate to phenotypic evolution. Our analyses, explicitly accounting for phylogenetic uncertainty, suggest diversification in phenotypic traits toward the present, with body shape evolution rapidly evolving in this group. Body shape evolution correlates with the occupation of different structural habitats indicated by vegetation axes suggesting species have adapted for maximal locomotory performance in these habitats. Our results also imply that the effects of phylogenetic uncertainty and model misspecification may be more extensive on univariate, relative to multivariate analyses of evolutionary correlations, which is an important consideration in analyzing data from rapidly radiating adaptive radiations.

Combining Species Delimitation, Species Trees, and Tests for Gene Flow Clarifies Complex Speciation in Scrub-Jays

Complex speciation, involving rapid divergence and multiple bouts of post-divergence gene flow, can obfuscate phylogenetic relationships and species limits. In North America, cases of complex speciation are common, due at least in part to the cyclical Pleistocene glacial history of the continent. Scrub-jays in the genus Aphelocoma provide a useful case study in complex speciation because their range throughout North America is structured by phylogeographic barriers with multiple cases of secondary contact between divergent lineages. Here, we show that a comprehensive approach to genomic reconstruction of evolutionary history, i.e., synthesizing results from species delimitation, species tree reconstruction, demographic model testing, and tests for gene flow, is capable of clarifying evolutionary history despite complex speciation. We find concordant evidence across all statistical approaches for the distinctiveness of an endemic southern Mexico lineage (A. w. sumichrasti), culminating in support for the species status of this lineage under any commonly applied species concept. We also find novel genomic evidence for the species status of a Texas endemic lineage A. w. texana, for which equivocal species delimitation results were clarified by demographic modeling and spatially explicit models of gene flow. Finally, we find that complex signatures of both ancient and modern gene flow between the non-sister California Scrub-Jay (A. californica) and Woodhouse's Scrub-Jay (A. woodhouseii), result in discordant gene trees throughout the species' genomes despite clear support for their overall isolation and species status. In sum, we find that a multi-faceted approach to genomic analysis can increase our understanding of complex speciation histories, even in well-studied groups. Given the emerging recognition that complex speciation is relatively commonplace, the comprehensive framework that we demonstrate for interrogation of species limits and evolutionary history using genomic data can provide a necessary roadmap for disentangling the impacts of gene flow and incomplete lineage sorting to better understand the systematics of other groups with similarly complex evolutionary histories.

Habitat-linked genetic variation supports microgeographic adaptive divergence in an island-endemic bird species

We investigated the potential mechanisms driving habitat-linked genetic divergence within a bird species endemic to a single 250 km² island. The island scrub-jay (Aphelocoma insularis) exhibits microgeographic divergence in bill morphology across pine-oak ecotones on Santa Cruz Island, California (USA) similar to adaptive differences described in mainland congeners over much larger geographic scales. To test whether individuals exhibit genetic differentiation related to habitat type and divergence in bill length, we genotyped over 3,000 single nucleotide polymorphisms (SNPs) in 123 adult island scrub-jay males from across Santa Cruz Island using restriction site-associated DNA sequencing (RADseq). Neutral landscape genomic analyses revealed that genome-wide genetic differentiation was primarily related to geographic distance and differences in habitat composition. We also found 168 putatively adaptive loci associated with habitat type using multivariate redundancy analysis (RDA) while controlling for spatial effects. Finally, two genome-wide association analyses revealed a polygenic basis to variation in bill length with multiple loci detected in or near genes known to affect bill morphology in other birds. Our findings support the hypothesis that divergent selection at microgeographic scales can cause adaptive divergence in the presence of ongoing gene flow.

Phylogenomics of alligator lizards elucidate diversification patterns across the Mexican Transition Zone and support the recognition of a new genus

Genomic data continue to advance our understanding of species limits and biogeographic patterns. However, there is still no consensus regarding appropriate methods of phylogenomic analysis that make the best use of these heterogeneous data sets. In this study, we used thousands of ultraconserved element (UCE) loci from alligator lizards in the genus Gerrhonotus to compare and contrast species trees inferred using multiple contemporary methods and provide a time frame for biological diversification across the Mexican Transition Zone (MTZ). Concatenated maximum likelihood (ML) and Bayesian analyses provided highly congruent results, with differences limited to poorly supported nodes. Similar topologies were inferred from coalescent analyses in Bayesian Phylogenetics and Phylogeography and SVDquartets, albeit with lower support for some nodes. All divergence times fell within the Miocene, linking speciation to local Neogene vicariance and/or global cooling trends following the mid-Miocene Climatic Optimum. We detected a high level of genomic divergence for a morphologically distinct species restricted to the arid mountains of north-eastern Mexico, and erected a new genus to better reflect evolutionary history. In summary, our results further advocate leveraging the strengths and weaknesses of concatenation and coalescent methods, provide evidence for old divergences for alligator lizards, and indicate that the MTZ continues to harbour substantial unrecognized diversity.

Repeated colonization, adaptive radiation and convergent evolution in the sheet-weaving spiders (Linyphiidae) of the south Pacific Archipelago of Juan Fernandez

We report on the colonization and diversification of linyphiid spiders in the Pacific oceanic archipelago of Juan Fernandez. About 50 spider species occur naturally in these islands, most of them endemic and about half of them are linyphiids. Linyphiidae includes no fewer than 15 species of Laminacauda and three of Neomaso (with several additional undescribed species in the latter genus), all of them single island endemics. There are three additional linyphiid endemic genera, two monotypic and one, Juanfernandezia, with two species. Unlike the rather uniform somatic morphology and small ground sheet webs of the continental Laminacauda and Neomaso species, the Juan Fernandez endemics exhibit morphological features and life history traits that are very rare or unknown in any other linyphiids. A multi-locus phylogenetic analysis confirms at least five independent Juan Fernandez colonizations of Linyphiidae, two within the same genus, and three of which underwent subsequent local diversification. Different calibrations suggest alternative colonization timelines, some at odds with island ages, but all agree on similar diversification timings of the endemic lineages. Rare phenotypic traits (e.g. gigantism, massive chelicerae or elongated legs) evolved multiple times independently within the islands. Based on the remarkable levels of eco-phenotypic differentiation in locally diversified species showing densely packed distributions, we propose that Laminacauda, and probably Neomaso, constitute a case of adaptive radiation.

Integrating Sequence Capture and Restriction Site-Associated DNA Sequencing to Resolve Recent Radiations of Pelagic Seabirds

The diversification of modern birds has been shaped by a number of radiations. Rapid diversification events make reconstructing the evolutionary relationships among taxa challenging due to the convoluted effects of incomplete lineage sorting (ILS) and introgression. Phylogenomic data sets have the potential to detect patterns of phylogenetic incongruence, and to address their causes. However, the footprints of ILS and introgression on sequence data can vary between different phylogenomic markers at different phylogenetic scales depending on factors such as their evolutionary rates or their selection pressures. We show that combining phylogenomic markers that evolve at different rates, such as paired-end double-digest restriction site-associated DNA (PE-ddRAD) and ultraconserved elements (UCEs), allows a comprehensive exploration of the causes of phylogenetic discordance associated with short internodes at different timescales. We used thousands of UCE and PE-ddRAD markers to produce the first well-resolved phylogeny of shearwaters, a group of medium-sized pelagic seabirds that are among the most phylogenetically controversial and endangered bird groups. We found that phylogenomic conflict was mainly derived from high levels of ILS due to rapid speciation events. We also documented a case of introgression, despite the high philopatry of shearwaters to their breeding sites, which typically limits gene flow. We integrated state-of-the-art concatenated and coalescent-based approaches to expand on previous comparisons of UCE and RAD-Seq data sets for phylogenetics, divergence time estimation, and inference of introgression, and we propose a strategy to optimize RAD-Seq data for phylogenetic analyses. Our results highlight the usefulness of combining phylogenomic markers evolving at different rates to understand the causes of phylogenetic discordance at different timescales. [Aves; incomplete lineage sorting; introgression; PE-ddRAD-Seq; phylogenomics; radiations; shearwaters; UCEs.].

Evolutionary history of Sundaland shrews (Eulipotyphla: Soricidae: Crocidura) with a focus on Borneo

The hyperdiverse shrew genus Crocidura is one of few small mammal genera distributed across Sundaland and all of its boundaries. This represents a rare opportunity to study the geological history of this region through the evolutionary history of these shrews. We generate a phylogeny of all recognized species of Sundaland Crocidura and show that most speciation events took place during the Pleistocene, prior to the inundation of the Sunda Shelf around 400 000 years ago. We find east–west differentiation within two separate lineages on Borneo, and that the current taxonomy of its two endemic species does not reflect evolutionary history, but ecophenotypic variation of plastic traits related to elevation. Sulawesi shrews are monophyletic, with a single notable exception: the black-footed shrew (C. nigripes). We show that the black-footed shrew diverged from its relatives on Borneo recently, suggesting a human-assisted breach of Wallace’s line. Overall, the number of Crocidura species, especially on Borneo, probably remains an underestimate.

Fossil-calibrated time tree of Podarcis wall lizards provides limited support for biogeographic calibration models

Podarcis wall lizards are endemic to the Mediterranean Basin where they represent the predominant reptile group. Despite being extensively used as model organisms in evolutionary and ecological studies their phylogeny and historical biogeography are still incompletely resolved. Moreover, molecular clock calibrations used in wall lizard phylogeography are based on the assumption of vicariant speciation triggered by the abrupt Mediterranean Sea level rise at the end of the Messinian salinity crisis (MSC). However, the validity of this biogeographic calibration remains untested. In this study we inferred a robust time tree based on multilocus data and fossil calibrations using both gene concatenation and species-tree approaches and including models with gene-flow. We found five deeply divergent, geographically coherent, and well-supported clades comprising species from i) Iberian Peninsula and North Africa; ii) Western Mediterranean islands, iii) Sicilian and Maltese islands; and iv-v) Balkan region and Aegean islands. The mitochondrial tree shows some inconsistencies with the species tree that warrant future investigation. Diversification of main clades is estimated in a short time frame during the Middle Miocene and might have been associated with a period of global climate cooling with the establishment of a marked climatic zonation in Europe. Cladogenetic events within the main clades are scattered throughout the time tree, from the Late Miocene to the Early Pleistocene, suggesting that speciation events in wall lizards reflect a complex interplay between regional topography, climate and geological history rather than a shared major climatic or paleogeographic event. Our absolute time estimates, as well as a relative dating approach, demonstrate that the assumption of a causal link between sea-level rise at the end of the MSC and the diversification of many island endemics is not justified. This study reinforces the notion that multiple dispersal and vicariant events, at different time frames, are required to explain current allopatric distributions and to account for the historical assembly of Mediterranean biota, and cautions against the use of biogeographic calibrations based on the assumption of vicariance.

Influence of niche breadth and position on the historical biogeography of seafaring scincid lizards

Niche breadth and position can influence diversification among closely related species or populations, yet limited empirical data exist concerning the predictability of the outcomes. We explored the effects of these factors on the evolution of the Emoia atrocostata species group, an insular radiation of lizards in the western Pacific Ocean and Indo-Australasia composed of both endemic and widespread species that differ in niche occupancy. We used molecular data and phylogeographical diffusion models to estimate the timing and patterns of range expansion, and ancestral reconstruction methods to infer shifts in ecology. We show evidence of multidirectional spread from a centre of origin in western Micronesia, and that the phyletic diversity of the group is derived from a putative habitat specialist that survives in the littoral zone. This species is composed of paraphyletic lineages that represent stages or possible endpoints in the continuum toward speciation. Several descendant species have transitioned to either strand or interior forest habitat, but only on remote islands with depauperate terrestrial faunas. Our results suggest that the atrocostata group might be in the early phases of a Wilsonian taxon cycle and that the capacity to tolerate salt stress has promoted dispersal and colonization of remote oceanic islands. Divergence itself, however, is largely driven by geographical isolation rather than shifts in ecology.

Swimming through the sands of the Sahara and Arabian deserts: Phylogeny of sandfish skinks (Scincidae, Scincus) reveals a recent and rapid diversification

Large parts of the Sahara Desert and Arabia are covered by sand seas and sand dunes, which are inhabited by specialized animal communities. For example, many lizards have developed adaptations to life in loose sand, including sand-swimming behavior. The best-known sand swimmers of the Saharo-Arabia are the sandfish skinks (genus Scincus). Although there are currently only four Scincus species recognized, their phylogenetic relationships have not yet been addressed in detail. We use eight genetic markers (three mitochondrial, five nuclear) and a complete sampling of species to infer the relationships within the genus. We employ multiple phylogenetic approaches to reconstruct the evolutionary history of these skinks and to assess the level of reticulation at the onset of their radiation. Our results indicate the presence of five strongly supported species-level lineages, four represented by the currently recognized species and the fifth by S. scincus conirostris, which does not form a clade with S. scincus. Based on these results we elevate the Iranian and northern Arabian S. conirostris to the species level. The two Saharan species, S. albifasciatus and S. scincus, are sister in all analyses. Deeper relationships within the genus, however, remained largely unresolved despite the extensive genetic data set. This basal polytomy, together with the fact that we detected no sign of hybridization in the history of the genus, indicates that the diversification of the five Scincus species was rapid, burst-like, and not followed by secondary hybridization events. Divergence time estimations show a Middle Pliocene crown radiation of the genus (3.3 Mya). We hypothesize that the aridification of the Saharo-Arabia that began in the Late Miocene triggered the initial diversification of Scincus, and that the subsequent expansion of sand deserts enabled their dispersal over the large Saharan and Arabian range. We discuss the evolution of body form in sand swimming lizards and ponder how Scincus retained their fully limbed morphology despite being sand swimmers that are typically limbless.

Competition and hybridization drive interspecific territoriality in birds

Costly interactions between species that arise as a by-product of ancestral similarities in communication signals are expected to persist only under specific evolutionary circumstances. Territorial aggression between species, for instance, is widely assumed to persist only when extrinsic barriers prevent niche divergence or selection in sympatry is too weak to overcome gene flow from allopatry. However, recent theoretical and comparative studies have challenged this view. Here we present a large-scale, phylogenetic analysis of the distribution and determinants of interspecific territoriality. We find that interspecific territoriality is widespread in birds and strongly associated with hybridization and resource overlap during the breeding season. Contrary to the view that territoriality only persists between species that rarely breed in the same areas or where niche divergence is constrained by habitat structure, we find that interspecific territoriality is positively associated with breeding habitat overlap and unrelated to habitat structure. Furthermore, our results provide compelling evidence that ancestral similarities in territorial signals are maintained and reinforced by selection when interspecific territoriality is adaptive. The territorial signals linked to interspecific territoriality in birds depend on the evolutionary age of interacting species, plumage at shallow (within-family) timescales, and song at deeper (between-family) timescales. Evidently, territorial interactions between species have persisted and shaped phenotypic diversity on a macroevolutionary timescale.

Multilocus phylogeny of Paratelmatobiinae (Anura: Leptodactylidae) reveals strong spatial structure and previously unknown diversity in the Atlantic Forest hotspot

The Brazilian Atlantic Forest harbors high levels of anuran diversity and endemism, including several taxa restricted to small geographic ranges. Here, we provide a multilocus phylogeny for Paratelmatobiinae, a leptodactylid subfamily composed of small-ranged species distributed in the Brazilian Atlantic Forest and in the campo rupestre ecosystem. We performed Bayesian inference and maximum likelihood analyses using three mitochondrial and five nuclear markers, and a matrix comprising a broad taxonomic sampling. We then delimitated independently evolving lineages within the group. We recovered Paratelmatobiinae and each of its four genera as monophyletic and robustly supported. Five putatively new species included in our analyses were unambiguously supported in the phylogenetic trees and delimitation analyses. We also recovered other deeply divergent and geographically structured lineages within the four genera of Paratelmatobiinae. Our estimation of divergence times indicates that diversification in the subfamily began in the Eocene and continued until the Pleistocene. We discuss possible scenarios of diversification for the four genera of Paratelmatobiinae, and outline the implications of our findings for taxonomy and conservation.

Cryptic and extensive hybridization between ancient lineages of American crows

Most species and therefore most hybrid zones have historically been defined using phenotypic characters. However, both speciation and hybridization can occur with negligible morphological differentiation. Recently developed genomic tools provide the means to better understand cryptic speciation and hybridization. The Northwestern Crow (Corvus caurinus) and American Crow (Corvus brachyrhynchos) are continuously distributed sister taxa that lack reliable traditional characters for identification. In this first population genomic study of Northwestern and American crows, we use genomic SNPs (nuDNA) and mtDNA to investigate the degree of genetic differentiation between these crows and the extent to which they may hybridize. Our results indicate that American and Northwestern crows have distinct evolutionary histories, supported by two nuDNA ancestry clusters and two 1.1%-divergent mtDNA clades dating to the late Pleistocene, when glacial advances may have isolated crow populations in separate refugia. We document extensive hybridization, with geographic overlap of mtDNA clades and admixture of nuDNA across >900 km of western Washington and western British Columbia. This broad hybrid zone consists of late-generation hybrids and backcrosses, but not recent (e.g., F1) hybrids. Nuclear DNA and mtDNA clines had concordant widths and were both centred in southwestern British Columbia, farther north than previously postulated. Overall, our results suggest a history of reticulate evolution in American and Northwestern crows, perhaps due to recurring neutral expansion(s) from Pleistocene glacial refugia followed by lineage fusion(s). However, we do not rule out a contributing role for more recent potential drivers of hybridization, such as expansion into human-modified habitats.

Eight in one: morphological and molecular analyses reveal cryptic diversity in Amazonian alopoglossid lizards (Squamata: Gymnophthalmoidea)

Among the currently recognized species of Alopoglossus, the Amazonian Alopoglossus angulatus has the widest distribution. We here analyse variation in scutellation and morphometrics of A. angulatus by examining 785 specimens of Alopoglossus. We also analyse intra- and interspecific genetic structure and differentiation using two mitochondrial (Cytb and ND4) and two nuclear (SNCAIP and PRLR) genes from 97 samples. Both morphological and molecular analyses are based on specimens and samples from Brazil, Colombia, Ecuador, French Guiana, Guyana, Peru and Suriname. Our results reveal A. angulatus as a monophyletic group composed of eight independently evolved lineages: A. angulatus s.s. plus three revalidated species plus two newly described species plus two putative species. We provide descriptions of all taxa, except for the putative species, including the first description of the neotype of A. angulatus and redescriptions of resurrected junior synonyms. Illustrations, diagnoses and geographical distribution maps are provided. Gene and species trees are also provided. The two new taxa recognized in this paper, along with the revalidation of three taxa, increase the total number of known species of Alopoglossus from nine to 14.

A precarious future for distinctive peripheral populations of meadow voles (Microtus pennsylvanicus)

Conservation efforts rely on robust taxonomic assessments that should be based on critical assessment of interspecific boundaries, infraspecific variation, and potentially distinctive peripheral populations. The meadow vole (Microtus pennsylvanicus) is widely distributed across North America, including 28 morphologically defined subspecies and numerous isolated populations. Because some subspecies are of high conservation concern, we examined geographic variation across the range of the species to test existing infraspecific taxonomy in terms of local and regional diversification. We sequenced mitochondrial DNA (mtDNA) from 20 subspecies of M. pennsylvanicus and contextualized infraspecific variation through comparison of pairwise genetic distances derived from an extended data set of 63 species of Microtus. We found strong support for at least three divergent clades within M. pennsylvanicus, with observed intraspecific clade divergence exceeding that between several pairwise comparisons of sister species within Microtus. Six nuclear genes were then sequenced to test the validity of mtDNA structure and to further evaluate the possibility of cryptic, species-level diversity using Bayes factor species delimitation (BFD) analyses. BFD consistently and decisively supported multiple species based on the multilocus approach. We propose that taxonomic revision of the meadow vole is required, with the eastern clade now identified as M. pennsylvanicus (Ord 1815), the western clade as M. drummondii (Audubon and Bachman 1853), and the coastal Florida clade as M. dukecampbelli (Woods, Post, and Kilpatrick 1982). We suggest that such an arrangement would more closely reflect evolutionary history and provide critical context for further examination of distinctive southern peripheral populations that harbor novel evolutionary legacies and adaptive potential.

The Genus Brueelia (Phthiraptera: Ischnocera: Philopteridae) of North American Jays and Allies (Aves: Passeriformes: Corvidae), with Descriptions of Five New Species

Five new species of chewing lice in the genus Brueelia Kéler, 1936 , are described from North American jays and allies. They are Brueelia mexicana n. sp. from Aphelocoma woodhouseii cyanotis Ridgway, 1887 ; Brueelia bonnevillensis n. sp. from Aphelocoma woodhouseii nevadae Pitelka, 1945 ; Brueelia diblasiae n. sp. from Cyanocitta stelleri frontalis ( Ridgway, 1873 ); Brueelia tempestwilliamsae n. sp. from Gymnorhinus cyanocephala Wied-Neuwied, 1841 ; Brueelia mcnewae n. sp. from Nucifraga columbiana ( Wilson, 1811 ). An identification key to the Brueelia on corvid hosts is provided.

Mitogenomic phylogenetic analyses of Leptogorgia virgulata and Leptogorgia hebes (Anthozoa: Octocorallia) from the Gulf of Mexico provides insight on Gorgoniidae divergence between Pacific and Atlantic lineages

The use of genetics in recent years has brought to light the need to reevaluate the classification of many gorgonian octocorals. This study focuses on two Leptogorgia species-Leptogorgia virgulata and Leptogorgia hebes-from the northwestern Gulf of Mexico (GOM). We target complete mitochondrial genomes and mtMutS sequences, and integrate this data with previous genetic research of gorgonian corals to resolve phylogenetic relationships and estimate divergence times. This study contributes the first complete mitochondrial genomes for L. ptogorgia virgulata and L. hebes. Our resulting phylogenies stress the need to redefine the taxonomy of the genus Leptogorgia in its entirety. The fossil-calibrated divergence times for Eastern Pacific and Western Atlantic Leptogorgia species based on complete mitochondrial genomes shows that the use of multiple genes results in estimates of more recent speciation events than previous research based on single genes. These more recent divergence times are in agreement with geologic data pertaining to the formation of the Isthmus of Panama.

Pre-quaternary diversification and glacial demographic expansions of Cardiocrinum (Liliaceae) in temperate forest biomes of Sino-Japanese Floristic Region

The Sino-Japanese Floristic Region (SJFR) in East Asia is one of the most diverse temperate floras in the world. However, the relative influence of Neogene palaeogeographical changes and Quaternary climatic fluctuations as causal mechanisms on species diversification remains largely controversial, because most divergence time estimates were inferred from single-locus data and have limited geographic or taxonomic sampling. To evaluate these influences, we use SNP markers from restriction site-associated DNA sequencing (RAD-Seq) loci and expressed sequence tags-simple sequence repeat (EST-SSR) markers to investigate the levels of genetic variation, speciation and demographic history of the temperate-deciduous forest (TDF) endemic Cardiocrinum (Endlicher) Lindley (Liliaceae), a genus comprising three species in China (C. giganteum, C. cathayanum) and Japan (C. cordatum). Phylogenomic and population genomic coalescent-based analyses demonstrated that Late Neogene tectonic/climatic events triggered speciation of Cardiocrinum, and Pleistocene climatic fluctuations had limited influence on its divergence history. Population demographic inference using Approximate Bayesian Computation from EST-SSRs and palaeoclimatic niche models both indicated that all three Cardiocrinum species experienced population expansions during the transition from the LIG to the LGM. We also discussed the implications of these results on the conservation of montane TDF species in the SJFR under ongoing environmental change. Our results improve our understanding of how the constituents of montane TDF across the SJFR responded to previous periods of rapid climate and environmental change in terms of speciation and population demographic processes.

Speciation of North American pygmy shrews (Eulipotyphla: Soricidae) supports spatial but not temporal congruence of diversification among boreal species

Speciation among many animals was rapid through the Pleistocene, impacted by climate and periodic isolation and reconnection. As such, species limits among often morphologically cryptic lineages may remain unresolved despite clear mitogenomic partitioning. Accumulating evidence from phylogeographical studies is revealing congruent regional differentiation of lineages across taxonomic groups that share ecological and evolutionary traits. Here, we analyse multiple DNA loci and morphology to resolve the geography and timeframe associated with evolutionary history of North American pygmy shrews (genus Sorex). We then assess lineage diversification among three co-distributed shrew complexes using phylogenetic and approximate Bayesian computation approaches to test a hypothesis of spatial congruence but temporal incongruence of species formation on a continental scale. Our results indicate consistency in regional lineage distributions, partial congruence of the sequence of divergence, and strong but not definitive support for temporal incongruence, suggesting that successive glacial cycles initiated the process of diversification repeatedly through the Pleistocene. Our results emphasize a continuing need for greater genomic coverage in comparative phylogeography, with persistent challenges. We recognize distinct eastern (Sorex hoyi Baird, 1857) and western (Sorex eximius Osgood, 1901) species of pygmy shrew based on available evidence, but discuss issues with taxonomic designations considering the continuum of speciation throughout the boreal biome.

Complete chloroplast genome sequence and phylogenetic analysis of wasabi (Eutrema japonicum) and its relatives

In Japan, two Eutrema species, wasabi (Eutrema japonicum, the important traditional Japanese condiment) and yuriwasabi (E. tenue), have been recognized as endemic species. We sequenced complete chloroplast (cp) genomes of seven wasabi and yuriwasabi accessions from Japan to study their phylogeny and evolution, using molecular dating of species divergence. Phylogenetic analyses of the complete cp DNA of these two Japanese species and five other Eurasian Eutrema species revealed that wasabi and yuriwasabi did not form a monophyletic group. One yuriwasabi accession (Gifu) formed a clade with E. yunnanense from China, indicating that this accession should be considered as a different species from the other yuriwasabi accessions. We reveal that Japanese Eutrema species diverged from the ‘E. yunnanense–yuriwasabi (Gifu)’ clade approximately 1.3 million years ago (Mya), suggesting that the connection between Japan and the Eurasian continent has existed more recently than the Quaternary period. The abundance of cp sequence data in this study also allowed the detection of genetic differentiation among wasabi cultivars. The two polymorphic sites detected between ‘Fujidaruma’ and ‘Shimane No.3’ were used to develop genotyping markers. The cp genome information provided here will thus inform the evolutionary histories of Japanese Eutrema species and help in genotyping wasabi cultivars.

Historical DNA as a tool to address key questions in avian biology and evolution: A review of methods, challenges, applications, and future directions

Museum specimens play a crucial role in addressing key questions in systematics, evolution, ecology, and conservation. With the advent of high-throughput sequencing technologies, specimens that have long been the foundation of important biological discoveries can inform new perspectives as sources of genomic data. Despite the many possibilities associated with analyzing DNA from historical specimens, several challenges persist. Using avian systems as a model, we review DNA extraction protocols, sequencing technologies, and capture methods that are helping researchers overcome some of these difficulties. We highlight empirical examples in which researchers have used these technologies to address fundamental questions related to avian conservation and evolution. Increasing accessibility to new sequencing technologies will provide researchers with tools to tap into the wealth of information contained within our valuable natural history collections.

Comparison of ultraconserved elements (UCEs) to microsatellite markers for the study of avian hybrid zones: a test in Aphelocoma jays

Objective

Hybrid zones are geographic regions where genetically distinct taxa interbreed, resulting in offspring of mixed ancestry. California Scrub-Jays (Aphelocoma californica) and Woodhouse’s Scrub-Jays (A. woodhouseii) come into secondary contact and hybridize in western Nevada. Although previous work investigated divergence and gene flow between these species using a handful of microsatellite markers, the hybrid zone has not been studied using genome-scale markers, providing an opportunity to assess genome-wide introgression, test for a genetic basis for ecomorphological traits, and compare these estimates to those derived from microsatellites.

Results

Using variant sites flanking ultraconserved elements (UCEs), we performed population assignment and quantified hybrid ancestry for 16 individuals across the zone of secondary contact. Our study included 2468 SNPs distributed throughout the genome, allowing discrimination of genetic affinities of hybrid individuals that were similar to estimates from microsatellites. We show a relationship between bill and wing length and the genetic composition of individuals that was not found in prior work using microsatellites, suggesting a genetic basis for these traits. Our analyses demonstrate the utility of UCEs for the analysis of hybrid zones and provide a basis for future studies to identify the genomic architecture of speciation and phenotypic differences between these incipient species.

Electronic supplementary material

The online version of this article (10.1186/s13104-019-4481-z) contains supplementary material, which is available to authorized users.

Fish diversification at the pace of geomorphological changes: evolutionary history of western Iberian Leuciscinae (Teleostei: Leuciscidae) inferred from multilocus sequence data

The evolutionary history of western Iberian Leuciscinae, obligatory freshwater fish, is directly linked to the evolution of the hydrographic network of the Iberian Peninsula after its isolation from the rest of Europe, which involved dramatic rearrangements such as the transition from endorheic lakes to open basins draining to the Atlantic. Previous phylogenetic research on western Iberian leuciscines, using mainly mitochondrial DNA and more recently one or two nuclear genes, has found contradictory results and there remain many unresolved issues regarding species relationships, taxonomy, and evolutionary history. Moreover, there is a lack of integration between phylogenetic and divergence time estimates and information on the timing of geomorphological changes and paleobasin rearrangements in the Iberian Peninsula. This study presents the first comprehensive fossil-calibrated multilocus coalescent species tree of western Iberian Leuciscinae (including 14 species of Achondrostoma, Iberochondrostoma, Pseudochondrostoma and Squalius endemic to the Iberian Peninsula, seven of which endemic to Portugal) based on seven nuclear genes, and from which we infer their biogeographic history by comparing divergence time estimates to known dated geological events. The phylogenetic pattern suggests slow-paced evolution of leuciscines during the Early-Middle Miocene endorheic phase of the main Iberian river basins, with the shift to exorheism in the late Neogene-Quaternary allowing westward dispersals that resulted in many cladogenetic events and a high rate of endemism in western Iberia. The results of this study also: (i) confirm the paraphyly of S. pyrenaicus with respect to S. carolitertii, and thus the possible presence of a new taxon in the Portuguese Tagus currently assigned to S. pyrenaicus; (ii) support the taxonomic separation of the Guadiana and Sado populations of S. pyrenaicus; (iii) show the need for further population sampling and taxonomic research to clarify the phylogenetic status of A. arcasii from the Minho basin and of the I. lusitanicum populations in the Sado and Tagus basins; and (iv) indicate that A. occidentale, I. olisiponensis and P. duriensis are the most ancient lineages within their respective genera.

Generation and classification of transcriptomes in two Croomia species and molecular evolution of CYC/TB1 genes in Stemonaceae

The genus Croomia (Stemonaceae) is an excellent model for studying the evolution of the Eastern Asia (EA)-Eastern North America (ENA) floristic disjunction and the genetic mechanisms of floral zygomorphy formation. In addition to the presence of both actinomorphic and zygomorphic flowers within the genus, species are disjunctively distributed between EA and ENA. However, due to the limited availability of genomic resources, few studies of Croomia have examined these questions. In this study, we sequenced the floral and leaf transcriptomes of the zygomorphic flowered C roomia heterosepala and the actinomorphic flowered Croomia japonica, and used comparative genomic approaches to investigate the transcriptome evolution of the two closely related species. The sequencing and de novo assembly of transcriptomes from flowers of C. heterosepala (ChFlower), flowers of C. japonica (CjFlower), and leaves of C. japonica (CjLeaf) yielded 57,193, 62,131 and 64,448 unigenes, respectively. In addition, estimation of Ka/Ks ratios for 11,566 potential orthologous groups between ChFlower and CjFlower revealed that only six pairs had Ka/Ks ratios significantly greater than 1 and are likely under positive selection. A total of 429 single copy nuclear genes (SCNGs) and 21,460 expression sequence tags-simple sequence repeats (EST-SSRs) were identified in this study. Specifically, we identified seven CYC/TB1-like genes from Stemonaceae. Phylogenetic and molecular evolution analyses indicated that these CYC/TB1-like genes formed a monophyletic clade (SteTBL1) and were subject to strong purifying selection. The shifts of floral symmetry in Stemonaceae do not appear to be correlated with TBL copy number.

Bayesian Divergence-Time Estimation with Genome-Wide Single-Nucleotide Polymorphism Data of Sea Catfishes (Ariidae) Supports Miocene Closure of the Panamanian Isthmus

The closure of the Isthmus of Panama has long been considered to be one of the best defined biogeographic calibration points for molecular divergence-time estimation. However, geological and biological evidence has recently cast doubt on the presumed timing of the initial isthmus closure around 3 Ma but has instead suggested the existence of temporary land bridges as early as the Middle or Late Miocene. The biological evidence supporting these earlier land bridges was based either on only few molecular markers or on concatenation of genome-wide sequence data, an approach that is known to result in potentially misleading branch lengths and divergence times, which could compromise the reliability of this evidence. To allow divergence-time estimation with genomic data using the more appropriate multispecies coalescent (MSC) model, we here develop a new method combining the single-nucleotide polymorphism-based Bayesian species-tree inference of the software SNAPP with a molecular clock model that can be calibrated with fossil or biogeographic constraints. We validate our approach with simulations and use our method to reanalyze genomic data of Neotropical army ants (Dorylinae) that previously supported divergence times of Central and South American populations before the isthmus closure around 3 Ma. Our reanalysis with the MSC model shifts all of these divergence times to ages younger than 3 Ma, suggesting that the older estimates supporting the earlier existence of temporary land bridges were artifacts resulting at least partially from the use of concatenation. We then apply our method to a new restriction-site associated DNA-sequencing data set of Neotropical sea catfishes (Ariidae) and calibrate their species tree with extensive information from the fossil record. We identify a series of divergences between groups of Caribbean and Pacific sea catfishes around 10 Ma, indicating that processes related to the emergence of the isthmus led to vicariant speciation already in the Late Miocene, millions of years before the final isthmus closure.

Evolutionary insights from comparative transcriptome and transcriptome-wide coalescence analyses in Tetrastigma hemsleyanum

BACKGROUND

Tetrastigma hemsleyanum is of great medicinal importance and used as a model system to address the evolutionary history of warm-temperate evergreen (WTE) forest biomes in East Asia over Neogene time scales. However, further studies on the neutral and adaptive divergence processes of T. hemsleyanum are currently impeded by a lack of genomic resources. In this study, we de novo assembled and annotated a reference transcriptome for two cpDNA lineages (Central-South-East vs. Southwest) of T. hemsleyanum. We further used comparative genomic and multilocus coalescent approaches to investigate the tempo and mode of lineage diversification in T. hemsleyanum.

RESULTS

A total of 52,838 and 65,197 unigenes with an N50 of 1,667 and 1,841 bp for Central-South-East (CSE) and Southwest (SW) lineages, respectively, were recovered, and 6,692 putative orthologs were identified between the two lineages. Estimation of K_a/K_s ratios for these orthologs revealed that ten genes had K_a/K_s values significantly greater than 0.5 (P < 0.05), whereas 2,099 (K_a/K_s < 0.5, P < 0.05) were inferred to be under purifying selection. Based on three bioinformatic strategies, we identified a total of 1,018 single-copy nuclear genes (SCNGs) from the orthologs. We successfully designed eight nuclear gene primer pairs with high intraspecific variation (e.g. h_T = 0.923, π_T = 1.68×10^-3), when surveyed across a subset of T. hemsleyanum individuals. Concordant with the previous cpDNA data, the haplotype networks constructed for most nuclear gene loci clearly identified the two lineages. A multilocus coalescence analysis suggested that the separation between the two lineages appears to have occurred during the mid-Pliocene. Despite their ancient divergence, both lineages experienced expansion at rather localized scales and have continued to exchange genes at a low rate.

CONCLUSIONS

This study demonstrated the utility of transcriptome sequencing as a basis for SCNG development in non-model species and the advantages of integrating multiple nuclear loci for phylogeographic and phylogenetic studies.

Unraveling historical introgression and resolving phylogenetic discord within Catostomus (Osteichthys: Catostomidae)

BACKGROUND

Porous species boundaries can be a source of conflicting hypotheses, particularly when coupled with variable data and/or methodological approaches. Their impacts can often be magnified when non-model organisms with complex histories of reticulation are investigated. One such example is the genus Catostomus (Osteichthys, Catostomidae), a freshwater fish clade with conflicting morphological and mitochondrial phylogenies. The former is hypothesized as reflecting the presence of admixed genotypes within morphologically distinct lineages, whereas the latter is interpreted as the presence of distinct morphologies that emerged multiple times through convergent evolution. We tested these hypotheses using multiple methods, to including multispecies coalescent and concatenated approaches. Patterson's D-statistic was applied to resolve potential discord, examine introgression, and test the putative hybrid origin of two species. We also applied naïve binning to explore potential effects of concatenation.

RESULTS

We employed 14,007 loci generated from ddRAD sequencing of 184 individuals to derive the first highly supported nuclear phylogeny for Catostomus. Our phylogenomic analyses largely agreed with a morphological interpretation,with the exception of the placement of Xyrauchen texanus, which differs from both morphological and mitochondrial phylogenies. Additionally, our evaluation of the putative hybrid species C. columbianus revealed a lack introgression and instead matched the mitochondrial phylogeny. Furthermore, D-statistic tests clarified all discrepancies based solely on mitochondrial data, with agreement among topologies derived from concatenation and multispecies coalescent approaches. Extensive historic introgression was detected across six species-pairs. Potential endemism in the Virgin and Little Colorado Rivers was also apparent, and the former genus Pantosteus was derived as monophyletic, save for C. columbianus.

CONCLUSIONS

Complex reticulated histories detected herein support the hypothesis that introgression was responsible for conflicts that occurred within the mitochondrial phylogeny, and explains discrepancies found between it and previous morphological phylogenies. Additionally, the hybrid origin of C. columbianus was refuted, but with the caveat that more fine-grain sampling is still needed. Our diverse phylogenomic approaches provided largely concordant results, with naïve binning useful in exploring the single conflict. Considerable diversity was found within Catostomus across southwestern North America, with two drainages [Virgin River (UT) and Little Colorado River (AZ)] reflecting unique composition.

A biogeographical profile of the sand cockroach Arenivaga floridensis and its bearing on origin hypotheses for Florida scrub biota

Florida scrub is a xeric ecosystem associated with the peninsula's sand ridges, whose intermittent Pliocene-Pleistocene isolation is considered key to scrub endemism. One scrub origin hypothesis posits endemics were sourced by the Pliocene dispersal of arid-adapted taxa from southwestern North America; a second invokes Pleistocene migration within eastern North America. Only one study to date has explicitly tested these competing hypotheses, supporting an eastern origin for certain scrub angiosperms. For further perspective, we conducted a genetic analysis of an endemic arthropod, the Florida sand cockroach (Arenivaga floridensis), with two aims: (1) to reconstruct the peninsular colonization and residence history of A. floridensis and (2) determine whether its biogeographic profile favors either origin hypothesis. We sequenced the cox2 mitochondrial gene for 237 specimens (65 populations) as well as additional loci (cox1, nuclear H3) for a subset of Florida roaches and congeners. Using Network and Bayesian inference methods, we identified three major lineages whose genetic differentiation and phylogeographical structure correspond with late Pliocene peninsula insularization, indicating Arenivaga was present and broadly distributed in Florida at that time. Stem and crown divergence estimates (6.36 Ma; 2.78 Ma) between A. floridensis and western sister taxa span a period of extensive dispersal by western biota along an arid Gulf Coast corridor. These phylogeographical and phylogenetic results yield a biogeographic profile consistent with the western origin hypothesis. Moreover, age estimates for the roach's peninsular residence complement those of several other endemics, favoring a Pliocene (or earlier) inception of the scrub ecosystem. We argue that eastern versus western hypotheses are not mutually exclusive; rather, a composite history of colonization involving disparate biotas better explains the diverse endemism of Florida scrub.