Shahan Derkarabetian

Repeated and Time-Correlated Morphological Convergence in Cave-Dwelling Harvestmen (Opiliones, Laniatores) from Montane Western North America

Background Many cave-dwelling animal species display similar morphologies (troglomorphism) that have evolved convergent within and among lineages under the similar selective pressures imposed by cave habitats. Here we study such ecomorphological evolution in cave-dwelling Sclerobuninae harvestmen (Opiliones) from the western United States, providing general insights into morphological homoplasy, rates of morphological change, and the temporal context of cave evolution. Methodology/Principal Findings We gathered DNA sequence data from three independent gene regions, and combined these data with Bayesian hypothesis testing, morphometrics analysis, study of penis morphology, and relaxed molecular clock analyses. Using multivariate morphometric analysis, we find that phylogenetically unrelated taxa have convergently evolved troglomorphism; alternative phylogenetic hypotheses involving less morphological convergence are not supported by Bayesian hypothesis testing. In one instance, this morphology is found in specimens from a high-elevation stony debris habitat, suggesting that troglomorphism can evolve in non-cave habitats. We discovered a strong positive relationship between troglomorphy index and relative divergence time, making it possible to predict taxon age from morphology. Most of our time estimates for the origin of highly-troglomorphic cave forms predate the Pleistocene. Conclusions/Significance While several regions in the eastern and central United States are well-known hotspots for cave evolution, few modern phylogenetic studies have addressed the evolution of cave-obligate species in the western United States. Our integrative studies reveal the recurrent evolution of troglomorphism in a perhaps unexpected geographic region, at surprisingly deep time depths, and in sometimes surprising habitats. Because some newly discovered troglomorphic populations represent undescribed species, our findings stress the need for further biological exploration, integrative systematic research, and conservation efforts in western US cave habitats.

Integrative Taxonomy and Species Delimitation in Harvestmen: A Revision of the Western North American Genus Sclerobunus (Opiliones: Laniatores: Travunioidea)

Alpha taxonomy, and specifically the delimitation of species, is becoming increasingly objective and integrative. The use of coalescent-based methods applied to genetic data is providing new tools for the discovery and delimitation of species. Here, we use an integrative approach via a combination of discovery-based multivariate morphological analyses to detect potential new species. These potential species are then used as a priori species in hypothesis-driven validation analyses with genetic data. This research focuses on the harvestmen genus Sclerobunus found throughout the mountainous regions of western North America. Based on our analyses, we conduct a revision of Sclerobunus resulting in synonymy of Cyptobunus with Sclerobunus including transfer of S. cavicolens comb. nov. and elevation of both subspecies of S. ungulatus: S. ungulatus comb. nov. and S. madhousensis comb. nov., stat. nov. The three subspecies of S. robustus are elevated, S. robustus, S. glorietus stat. nov., and S. idahoensis stat. nov. Additionally, five new species of Sclerobunus are described from New Mexico and Colorado, including S. jemez sp. nov., S. klomax sp. nov., S. skywalkeri sp. nov., S. speoventus sp. nov., and S. steinmanni sp. nov. Several of the newly described species are single-cave endemics, and our findings suggest that further exploration of western North American cave habitats will likely yield additional new species.

High Phylogenetic Utility of an Ultraconserved Element Probe Set Designed for Arachnida

Arachnida is an ancient, diverse and ecologically important animal group that contains a number of species of interest for medical, agricultural and engineering applications. Despite their importance, many aspects of the arachnid tree of life remain unresolved, hindering comparative approaches to arachnid biology. Biologists have made considerable efforts to resolve the arachnid phylogeny; yet, limited and challenging morphological characters, as well as a dearth of genetic resources, have hindered progress. Here, we present a genomic toolkit for arachnids featuring hundreds of conserved DNA regions (ultraconserved elements or UCEs) that allow targeted sequencing of any species in the arachnid tree of life. We used recently developed capture probes designed from conserved regions of available arachnid genomes to enrich a sample of loci from 32 diverse arachnids. Sequence capture returned an average of 487 UCE loci for all species, with a range from 170 to 722. Phylogenetic analysis of these UCEs produced a highly resolved arachnid tree with relationships largely consistent with recent transcriptome‐based phylogenies. We also tested the phylogenetic informativeness of UCE probes within the spider, scorpion and harvestman orders, demonstrating the utility of these markers at shallower taxonomic scales and suggesting that these loci will be useful for species‐level differences. This probe set will open the door to phylogenomic and population genomic studies across the arachnid tree of life, enabling systematics, species delimitation, species discovery and conservation of these diverse arthropods.

Population genomic evidence for multiple Pliocene refugia in a montane-restricted harvestman (Arachnida, Opiliones, Sclerobunus robustus) from the southwestern United States

The integration of ecological niche modelling into phylogeographic analyses has allowed for the identification and testing of potential refugia under a hypothesis‐based framework, where the expected patterns of higher genetic diversity in refugial populations and evidence of range expansion of nonrefugial populations are corroborated with empirical data. In this study, we focus on a montane‐restricted cryophilic harvestman, Sclerobunus robustus, distributed throughout the heterogeneous Southern Rocky Mountains and Intermontane Plateau of southwestern North America. We identified hypothetical refugia using ecological niche models (ENMs) across three time periods, corroborated these refugia with population genetic methods using double‐digest RAD‐seq data and conducted population‐level phylogenetic and divergence dating analyses. ENMs identify two large temporally persistent regions in the mid‐latitude highlands. Genetic patterns support these two hypothesized refugia with higher genetic diversity within refugial populations and evidence for range expansion in populations found outside hypothesized refugia. Phylogenetic analyses identify five to six genetically divergent, geographically cohesive clades of S. robustus. Divergence dating analyses suggest that these separate refugia date to the Pliocene and that divergence between clades pre‐dates the late Pleistocene glacial cycles, while diversification within clades was likely driven by these cycles. Population genetic analyses reveal effects of both isolation by distance (IBD) and isolation by environment (IBE), with IBD more important in the continuous mountainous portion of the distribution, while IBE was stronger in the populations inhabiting the isolated sky islands of the south. Using model‐based coalescent approaches, we find support for postdivergence migration between clades from separate refugia.

Phylogenomic reclassification of the world’s most venomous spiders (Mygalomorphae, Atracinae), with implications for venom evolution

Here we show that the most venomous spiders in the world are phylogenetically misplaced. Australian atracine spiders (family Hexathelidae), including the notorious Sydney funnel-web spider Atrax robustus, produce venom peptides that can kill people. Intriguingly, eastern Australian mouse spiders (family Actinopodidae) are also medically dangerous, possessing venom peptides strikingly similar to Atrax hexatoxins. Based on the standing morphology-based classification, mouse spiders are hypothesized distant relatives of atracines, having diverged over 200 million years ago. Using sequence-capture phylogenomics, we instead show convincingly that hexathelids are non-monophyletic, and that atracines are sister to actinopodids. Three new mygalomorph lineages are elevated to the family level, and a revised circumscription of Hexathelidae is presented. Re-writing this phylogenetic story has major implications for how we study venom evolution in these spiders, and potentially genuine consequences for antivenom development and bite treatment research. More generally, our research provides a textbook example of the applied importance of modern phylogenomic research.

Comparative performance of double-digest RAD sequencing across divergent arachnid lineages

Next‐generation sequencing technologies now allow researchers of non‐model systems to perform genome‐based studies without the requirement of a (often unavailable) closely related genomic reference. We evaluated the role of restriction endonuclease (RE) selection in double‐digest restriction‐site‐associated DNA sequencing (ddRADseq) by generating reduced representation genome‐wide data using four different RE combinations. Our expectation was that RE selections targeting longer, more complex restriction sites would recover fewer loci than RE with shorter, less complex sites. We sequenced a diverse sample of non‐model arachnids, including five congeneric pairs of harvestmen (Opiliones) and four pairs of spiders (Araneae). Sample pairs consisted of either conspecifics or closely related congeneric taxa, and in total 26 sample pair analyses were tested. Sequence demultiplexing, read clustering and variant calling were performed in the pyRAD program. The 6‐base pair cutter EcoRI combined with methylated site‐specific 4‐base pair cutter MspI produced, on average, the greatest numbers of intra‐individual loci and shared loci per sample pair. As expected, the number of shared loci recovered for a sample pair covaried with the degree of genetic divergence, estimated with cytochrome oxidase I sequences, although this relationship was non‐linear. Our comparative results will prove useful in guiding protocol selection for ddRADseq experiments on many arachnid taxa where reference genomes, even from closely related species, are unavailable.

A new monster from southwest Oregon forests: Cryptomaster behemoth sp. n. (Opiliones, Laniatores, Travunioidea)

Abstract The monotypic genus Cryptomaster Briggs, 1969 was described based on individuals from a single locality in southwestern Oregon. The described species Cryptomaster leviathan Briggs, 1969 was named for its large body size compared to most travunioid Laniatores. However, as the generic name suggests, Cryptomaster are notoriously difficult to find, and few subsequent collections have been recorded for this genus. Here, we increase sampling of Cryptomaster to 15 localities, extending their known range from the Coast Range northeast to the western Cascade Mountains of southern Oregon. Phylogenetic analyses of mitochondrial and nuclear DNA sequence data reveal deep phylogenetic breaks consistent with independently evolving lineages. We use discovery and validation species delimitation approaches to generate and test species hypotheses, including a coalescent species delimitation method to test multi-species hypotheses. For delimited species, we use light microscopy and SEM to discover diagnostic morphological characters. Although Cryptomaster has a small geographic distribution, this taxon is consistent with other short-range endemics in having deep phylogenetic breaks indicative of species level divergences. Herein we describe Cryptomaster behemoth sp. n., and provide morphological diagnostic characters for identifying Cryptomaster leviathan and Cryptomaster behemoth.

Cryptic elevational zonation in trapdoor spiders (Araneae, Antrodiaetidae, Aliatypus janus complex) from the California southern Sierra Nevada

The relative roles of ecological niche conservatism versus niche divergence in promoting montane speciation remains an important topic in biogeography. Here, our aim was to test whether lineage diversification in a species complex of trapdoor spiders corresponds with riverine barriers or with an ecological gradient associated with elevational tiering. Aliatypus janus was sampled from throughout its range, with emphasis on populations in the southern Sierra Nevada Mountains of California. We collected multi-locus genetic data to generate a species tree for A. janus and its close relatives. Coalescent based hypothesis tests were conducted to determine if genetic breaks within A. janus conform to riverine barriers. Ecological niche models (ENM) under current and Last Glacial Maximum (LGM) conditions were generated and hypothesis tests of niche conservatism and divergence were performed. Coalescent analyses reveal deeply divergent genetic lineages within A. janus, likely corresponding to cryptic species. Two primary lineages meet along an elevational gradient on the western slopes of the southern Sierra Nevada Mountains. ENMs under both current and LGM conditions indicate that these groups occupy largely non-overlapping niches. ENM hypothesis testing rejected niche identity between the two groups, and supported a sharp ecological gradient occurring where the groups meet. However, the niche similarity test indicated that the two groups may not inhabit different background niches. The Sierra Nevada Mountains provide a natural laboratory for simultaneously testing ecological niche divergence and conservatism and their role in speciation across a diverse range of taxa. Aliatypus janus represents a species complex with cryptic lineages that may have diverged due to parapatric speciation along an ecological gradient, or been maintained by the evolution of ecological niche differences following allopatric speciation.

Sequence capture phylogenomics of eyeless Cicurina spiders from Texas caves, with emphasis on US federally-endangered species from Bexar County (Araneae, Hahniidae)

Abstract Morphological, mitochondrial, and nuclear phylogenomic data were combined to address phylogenetic and species delimitation questions in cave-limited Cicurina spiders from central Texas. Special effort was focused on specimens and cave locations in the San Antonio region (Bexar County), home to four eyeless species listed as US Federally Endangered. Sequence capture experiments resulted in the recovery of ~200–400 homologous ultra-conserved element (UCE) nuclear loci across taxa, and nearly complete COI mitochondrial DNA sequences from the same set of individuals. Some of these nuclear and mitochondrial sequences were recovered from “standard” museum specimens without special preservation of DNA material, including museum specimens preserved in the 1990s. Multiple phylogenetic analyses of the UCE data agree in the recovery of two major lineages of eyeless Cicurina in Texas. These lineages also differ in mitochondrial clade membership, female genitalic morphology, degree of troglomorphy (as measured by relative leg length), and are mostly allopatric across much of Texas. Rare sympatry was confirmed in Bexar County, where members of the two major clades sometimes co-exist in the same karst feature. Both nuclear phylogenomic and mitochondrial data indicate the existence of undescribed species from the San Antonio region, although further sampling and collection of adult specimens is needed to explicitly test these hypotheses. Our data support the two following species synonymies (Cicurina venii Gertsch, 1992 = Cicurina madla Gertsch, 1992; Cicurina loftini Cokendolpher, 2004 = Cicurina vespera Gertsch, 1992), formally proposed here. Overall, our taxonomy-focused research has many important conservation implications, and again highlights the fundamental importance of robust taxonomy in conservation research.

A stable phylogenomic classification of Travunioidea (Arachnida, Opiliones, Laniatores) based on sequence capture of ultraconserved elements

Abstract Molecular phylogenetics has transitioned into the phylogenomic era, with data derived from next-generation sequencing technologies allowing unprecedented phylogenetic resolution in all animal groups, including understudied invertebrate taxa. Within the most diverse harvestmen suborder, Laniatores, most relationships at all taxonomic levels have yet to be explored from a phylogenomics perspective. Travunioidea is an early-diverging lineage of laniatorean harvestmen with a Laurasian distribution, with species distributed in eastern Asia, eastern and western North America, and south-central Europe. This clade has had a challenging taxonomic history, but the current classification consists of ~77 species in three families, the Travuniidae, Paranonychidae, and Nippononychidae. Travunioidea classification has traditionally been based on structure of the tarsal claws of the hind legs. However, it is now clear that tarsal claw structure is a poor taxonomic character due to homoplasy at all taxonomic levels. Here, we utilize DNA sequences derived from capture of ultraconserved elements (UCEs) to reconstruct travunioid relationships. Data matrices consisting of 317–677 loci were used in maximum likelihood, Bayesian, and species tree analyses. Resulting phylogenies recover four consistent and highly supported clades; the phylogenetic position and taxonomic status of the enigmatic genus Yuria is less certain. Based on the resulting phylogenies, a revision of Travunioidea is proposed, now consisting of the Travuniidae, Cladonychiidae, Paranonychidae (Nippononychidae is synonymized), and the new family Cryptomastridae Derkarabetian & Hedin, fam. n., diagnosed here. The phylogenetic utility and diagnostic features of the intestinal complex and male genitalia are discussed in light of phylogenomic results, and the inappropriateness of the tarsal claw in diagnosing higher-level taxa is further corroborated.