Next-generation museum genomics: Phylogenetic relationships among palpimanoid spiders using sequence capture techniques (Araneae: Palpimanoidea)

Robust phylogenomics settles controversies of classification and reveals evolution of male embolic complex of the Laufeia clade (Araneae, Salticidae, Euophryini)

The Laufeia clade is a peculiar lineage of euophryine jumping spiders showing rapid divergence of male genital structures, especially the embolic complex that directly interacts with female genitalia during sperm transfer. With the rapid growth of species discovery and the perplexing morphology of male genitalia in the Laufeia clade, the controversy in its classification has become a crucial problem. In this study, we applied a phylogenomic approach using ultra-conserved elements data to infer the phylogeny of the Laufeia clade with extensive taxon sampling. A comparative morphological study was performed to evaluate diagnostic characters and understand the evolution of the male embolic complex within the Laufeia clade. The evolution of microhabitats (foliage, tree trunk, rock and surface litter) was also investigated to uncover the potential link between the microhabitat shifts and male embolic complex divergence. The results provide a strongly supported phylogenetic framework and updated generic concepts for the Laufeia clade. The synapomorphies for the updated genera within the Laufeia clade were identified through character mapping on the phylogeny. Ancestral state reconstruction analyses revealed that the Type I embolic complex (characterized by a disc-like embolic disc with a lamina as its outer edge) was ancestral and gradually evolved into the Type II (without lamina of embolic disc, base of embolic complex often modified into a functional "conductor") and Type III (lacking lamina of embolic disc and base of embolic complex) embolic complex, and that some embolic shapes evolved multiple times independently in different lineages of the Laufeia clade. The shift from foliage-dwelling to tree trunk-dwelling in the common ancestor of the Laufeia clade may have facilitated the divergent evolution of male embolic complex in the Laufeia clade. This study provides a solid foundation for future studies of systematics and evolution of this group.

Species delimitation with limited sampling: An example from rare trapdoor spider genus Cyclocosmia (Mygalomorphae, Halonoproctidae)

The outcome of species delimitation depends on many factors, including conceptual framework, study design, data availability, methodology employed and subjective decision making. Obtaining sufficient taxon sampling in endangered or rare taxa might be difficult, particularly when non-lethal tissue collection cannot be utilized. The need to avoid overexploitation of the natural populations may thus limit methodological framework available for downstream data analyses and bias the results. We test species boundaries in rare North American trapdoor spider genus Cyclocosmia Ausserer (1871) inhabiting the Southern Coastal Plain biodiversity hotspot with the use of genomic data and two multispecies coalescent model methods. We evaluate the performance of each methodology within a limited sampling framework. To mitigate the risk of species over splitting, common in taxa with highly structured populations, we subsequently implement a species validation step via genealogical diversification index (gdi), which accounts for both genetic isolation and gene flow. We delimited eight geographically restricted lineages within sampled North American Cyclocosmia, suggesting that major river drainages in the region are likely barriers to dispersal. Our results suggest that utilizing BPP in the species discovery step might be a good option for datasets comprising hundreds of loci, but fewer individuals, which may be a common scenario for rare taxa. However, we also show that such results should be validated via gdi, in order to avoid over splitting.

Burma Terrane Amber Fauna Shows Connections to Gondwana and Transported Gondwanan Lineages to the Northern Hemisphere (Araneae: Palpimanoidea)

Burmese amber is a significant source of fossils that documents the mid-Cretaceous biota. This deposit was formed around 99 Ma on the Burma Terrane, which broke away from Gondwana and later collided with Asia, although the timing is disputed. Palpimanoidea is a dispersal-limited group that was a dominant element of the Mesozoic spider fauna, and has an extensive fossil record, particularly from Burmese amber. Using morphological and molecular data, evolutionary relationships of living and fossil Palpimanoidea are examined. Divergence dating with fossils as terminal tips shows timing of diversification is contemporaneous with continental breakup.Ancestral range estimations show widespread ancestral ranges that divide into lineages that inherit different Pangean fragments, consistent with vicariance. Our results suggest that the Burmese amber fauna has ties to Gondwana due to a historical connection in the Early Cretaceous, and that the Burma Terrane facilitated biotic exchange by transporting lineages from Gondwana into the Holarctic in the Cretaceous.

Advances in the reconstruction of the spider tree of life: A roadmap for spider systematics and comparative studies

In the last decade and a half, advances in genetic sequencing technologies have revolutionized systematics, transforming the field from studying morphological characters or a few genetic markers, to genomic datasets in the phylogenomic era. A plethora of molecular phylogenetic studies on many taxonomic groups have come about, converging on, or refuting prevailing morphology or legacy-marker-based hypotheses about evolutionary affinities. Spider systematics has been no exception to this transformation and the inter-relationships of several groups have now been studied using genomic data. About 51 500 extant spider species have been described, all with a conservative body plan, but innumerable morphological and behavioural peculiarities. Inferring the spider tree of life using morphological data has been a challenging task. Molecular data have corroborated many hypotheses of higher-level relationships, but also resulted in new groups that refute previous hypotheses. In this review, we discuss recent advances in the reconstruction of the spider tree of life and highlight areas where additional effort is needed with potential solutions. We base this review on the most comprehensive spider phylogeny to date, representing 131 of the 132 spider families. To achieve this sampling, we combined six Sanger-based markers with newly generated and publicly available genome-scale datasets. We find that some inferred relationships between major lineages of spiders (such as Austrochiloidea, Palpimanoidea and Synspermiata) are robust across different classes of data. However, several new hypotheses have emerged with different classes of molecular data. We identify and discuss the robust and controversial hypotheses and compile this blueprint to design future studies targeting systematic revisions of these problematic groups. We offer an evolutionary framework to explore comparative questions such as evolution of venoms, silk, webs, morphological traits and reproductive strategies.

An Opiliones-specific ultraconserved element probe set with a near-complete family-level phylogeny

Sequence capture of ultraconserved elements (UCEs) has transformed molecular systematics across many taxa, with arachnids being no exception. The probe set available for Arachnida has been repeatedly used across multiple arachnid lineages and taxonomic levels, however more specific probe sets for spiders have demonstrated that more UCEs can be recovered with higher probe specificity. In this study, we develop an Opiliones-specific UCE probe set targeting 1915 UCEs using a combination of probes designed from genomes and transcriptomes, as well as the most useful probes from the Arachnida probe set. We demonstrate the effectiveness of this probe set across Opiliones with the most complete family-level phylogeny made to date, including representatives from 61 of 63 currently described families. We also test UCE recovery from historical specimens with degraded DNA, examine population-level data sets, and assess "backwards compatibility" with samples hybridized with the Arachnida probe set. The resulting phylogenies - which include specimens hybridized using both the Opiliones and Arachnida probe sets, historical specimens, and transcriptomes - are largely congruent with previous multi-locus and phylogenomic analyses. The probe set is also "backwards compatible", increasing the number of loci obtained in samples previously hybridized with the Arachnida probe set, and shows high utility down to shallow population-level divergences. This probe set has the potential to further transform Opiliones molecular systematics, resolving many long-standing taxonomic issues plaguing this lineage.

A novel probe set for the phylogenomics and evolution of RTA spiders

Spiders are important models for evolutionary studies of web building, sexual selection and adaptive radiation. The recent development of probes for UCE (ultra-conserved element)-based phylogenomic studies has shed light on the phylogeny and evolution of spiders. However, the two available UCE probe sets for spider phylogenomics (Spider and Arachnida probe sets) have relatively low capture efficiency within spiders, and are not optimized for the retrolateral tibial apophysis (RTA) clade, a hyperdiverse lineage that is key to understanding the evolution and diversification of spiders. In this study, we sequenced 15 genomes of species in the RTA clade, and using eight reference genomes, we developed a new UCE probe set (41 845 probes targeting 3802 loci, labelled as the RTA probe set). The performance of the RTA probes in resolving the phylogeny of the RTA clade was compared with the Spider and Arachnida probes through an in-silico test on 19 genomes. We also tested the new probe set empirically on 28 spider species of major spider lineages. The results showed that the RTA probes recovered twice and four times as many loci as the other two probe sets, and the phylogeny from the RTA UCEs provided higher support for certain relationships. This newly developed UCE probe set shows higher capture efficiency empirically and is particularly advantageous for phylogenomic and evolutionary studies of RTA clade and jumping spiders.

Predictors of sequence capture in a large-scale anchored phylogenomics project

Next-generation sequencing (NGS) technologies have revolutionized phylogenomics by decreasing the cost and time required to generate sequence data from multiple markers or whole genomes. Further, the fragmented DNA of biological specimens collected decades ago can be sequenced with NGS, reducing the need for collecting fresh specimens. Sequence capture, also known as anchored hybrid enrichment, is a method to produce reduced representation libraries for NGS sequencing. The technique uses single-stranded oligonucleotide probes that hybridize with pre-selected regions of the genome that are sequenced via NGS, culminating in a dataset of numerous orthologous loci from multiple taxa. Phylogenetic analyses using these sequences have the potential to resolve deep and shallow phylogenetic relationships. Identifying the factors that affect sequence capture success could save time, money, and valuable specimens that might be destructively sampled despite low likelihood of sequencing success. We investigated the impacts of specimen age, preservation method, and DNA concentration on sequence capture (number of captured sequences and sequence quality) while accounting for taxonomy and extracted tissue type in a large-scale butterfly phylogenomics project. This project used two probe sets to extract 391 loci or a subset of 13 loci from over 6,000 butterfly specimens. We found that sequence capture is a resilient method capable of amplifying loci in samples of varying age (0–111 years), preservation method (alcohol, papered, pinned), and DNA concentration (0.020 ng/μl - 316 ng/ul). Regression analyses demonstrate that sequence capture is positively correlated with DNA concentration. However, sequence capture and DNA concentration are negatively correlated with sample age and preservation method. Our findings suggest that sequence capture projects should prioritize the use of alcohol-preserved samples younger than 20 years old when available. In the absence of such specimens, dried samples of any age can yield sequence data, albeit with returns that diminish with increasing age.

Sequence Capture From Historical Museum Specimens: Maximizing Value for Population and Phylogenomic Studies

The application of high-throughput, short-read sequencing to degraded DNA has greatly increased the feasibility of generating genomic data from historical museum specimens. While many published studies report successful sequencing results from historical specimens; in reality, success and quality of sequence data can be highly variable. To examine predictors of sequencing quality, and methodological approaches to improving data accuracy, we generated and analyzed genomic sequence data from 115 historically collected museum specimens up to 180 years old. Data span both population genomic and phylogenomic scales, including historically collected specimens from 34 specimens of four species of Australian rock-wallabies (genus Petrogale) and 92 samples from 79 specimens of Australo-Papuan murine rodents (subfamily Murinae). For historical rodent specimens, where the focus was sampling for phylogenomics, we found that regardless of specimen age, DNA sequence libraries prepared from toe pad or bone subsamples performed significantly better than those taken from the skin (in terms of proportion of reads on target, number of loci captured, and data accuracy). In total, 93% of DNA libraries from toe pad or bone subsamples resulted in reliable data for phylogenetic inference, compared to 63% of skin subsamples. For skin subsamples, proportion of reads on target weakly correlated with collection year. Then using population genomic data from rock-wallaby skins as a test case, we found substantial improvement in final data quality by mapping to a high-quality “closest sister” de novo assembly from fresh tissues, compared to mapping to a sample-specific historical de novo assembly. Choice of mapping approach also affected final estimates of the number of segregating sites and Watterson's θ, both important parameters for population genomic inference. The incorporation of accurate and reliable sequence data from historical specimens has important outcomes for evolutionary studies at both population and phylogenomic scales. By assessing the outcomes of different approaches to specimen subsampling, library preparation and bioinformatic processing, our results provide a framework for increasing sequencing success for irreplaceable historical specimens.

Maximizing Molecular Data From Low-Quality Fluid-Preserved Specimens in Natural History Collections

Over the past decade, museum genomics studies have focused on obtaining DNA of sufficient quality and quantity for sequencing from fluid-preserved natural history specimens, primarily to be used in systematic studies. While these studies have opened windows to evolutionary and biodiversity knowledge of many species worldwide, published works often focus on the success of these DNA sequencing efforts, which is undoubtedly less common than obtaining minimal or sometimes no DNA or unusable sequence data from specimens in natural history collections. Here, we attempt to obtain and sequence DNA extracts from 115 fresh and 41 degraded samples of homalopsid snakes, as well as from two degraded samples of a poorly known snake, Hydrablabes periops. Hydrablabes has been suggested to belong to at least two different families (Natricidae and Homalopsidae) and with no fresh tissues known to be available, intractable museum specimens currently provide the only opportunity to determine this snake’s taxonomic affinity. Although our aim was to generate a target-capture dataset for these samples, to be included in a broader phylogenetic study, results were less than ideal due to large amounts of missing data, especially using the same downstream methods as with standard, high-quality samples. However, rather than discount results entirely, we used mapping methods with references and pseudoreferences, along with phylogenetic analyses, to maximize any usable molecular data from our sequencing efforts, identify the taxonomic affinity of H. periops, and compare sequencing success between fresh and degraded tissue samples. This resulted in largely complete mitochondrial genomes for five specimens and hundreds to thousands of nuclear loci (ultra-conserved loci, anchored-hybrid enrichment loci, and a variety of loci frequently used in squamate phylogenetic studies) from fluid-preserved snakes, including a specimen of H. periops from the Field Museum of Natural History collection. We combined our H. periops data with previously published genomic and Sanger-sequenced datasets to confirm the familial designation of this taxon, reject previous taxonomic hypotheses, and make biogeographic inferences for Hydrablabes. A second H. periops specimen, despite being seemingly similar for initial raw sequencing results and after being put through the same protocols, resulted in little usable molecular data. We discuss the successes and failures of using different pipelines and methods to maximize the products from these data and provide expectations for others who are looking to use DNA sequencing efforts on specimens that likely have degraded DNA.

Life Science Identifier (Hydrablabes periops)

urn:lsid:zoobank.org:pub:F2AA44 E2-D2EF-4747-972A-652C34C2C09D.

Metagenomic clustering reveals microbial contamination as an essential consideration in ultraconserved element design for phylogenomics with insect museum specimens

Phylogenomics via ultraconserved elements (UCEs) has led to improved phylogenetic reconstructions across the tree of life. However, inadvertently incorporating non-targeted DNA into the UCE marker design will lead to misinformation being incorporated into subsequent analyses. To date, the effectiveness of basic metagenomic filtering strategies has not been assessed in arthropods. Designing markers from museum specimens requires careful consideration of methods due to the high levels of microbial contamination typically found in such specimens. We investigate if contaminant sequences are carried forward into a UCE marker set we developed from insect museum specimens using a standard bioinformatics pipeline. We find that the methods currently employed by most researchers do not exclude contamination from the final set of targets. Lastly, we highlight several paths forward for reducing contamination in UCE marker design.

How to build phylogenetic species trees with OMA

Knowledge of species phylogeny is critical to many fields of biology. In an era of genome data availability, the most common way to make a phylogenetic species tree is by using multiple protein-coding genes, conserved in multiple species. This methodology is composed of several steps: orthology inference, multiple sequence alignment and inference of the phylogeny with dedicated tools. This can be a difficult task, and orthology inference, in particular, is usually computationally intensive and error prone if done ad hoc. This tutorial provides protocols to make use of OMA Orthologous Groups, a set of genes all orthologous to each other, to infer a phylogenetic species tree. It is designed to be user-friendly and computationally inexpensive, by providing two options: (1) Using only precomputed groups with species available on the OMA Browser, or (2) Computing orthologs using OMA Standalone for additional species, with the option of using precomputed orthology relations for those present in OMA. A protocol for downstream analyses is provided as well, including creating a supermatrix, tree inference, and visualization. All protocols use publicly available software, and we provide scripts and code snippets to facilitate data handling. The protocols are accompanied with practical examples.

Phylogeny and secondary sexual trait evolution in Schizocosa wolf spiders (Araneae, Lycosidae) shows evidence for multiple gains and losses of ornamentation and species delimitation uncertainty

Members of the Nearctic spider genus Schizocosa Chamberlin, 1904 have garnered much attention in behavioral studies and over many decades, a number of species have developed as model systems for investigating patterns of sexual selection and multimodal communication. Many of these studies have employed a comparative approach using putative, but not rigorously tested, sister species pairs that have distinctive morphological traits and attendant behaviors. Despite past emphasis on the efficacy of these presumably comparative-based studies of closely related species, generating a robust phylogenetic hypothesis for Schizocosa has been an ongoing challenge. Here, we apply a phylogenomic approach using anchored hybrid enrichment to generate a data set comprising over 400 loci representing a comprehensive taxonomic sample of 23 Nearctic Schizocosa. Our sampling also includes numerous outgroup lycosid genera that allow for a robust evaluation of genus monophyly. Based on analyses using concatenation and coalescent-based methods, we recover a well-supported phylogeny that infers the following: 1) The New World Schizocosa do not form a monophyletic group; 2) Previous hypotheses of North American species require reconsideration along with the composition of species groups; 3) Multiple longstanding model species are not genealogically exclusive and thus are not "good" species; 4) This updated phylogenetic framework establishes a new working paradigm for studying the evolution of characters associated with reproductive communication and mating. Ancestral character state reconstructions show a complex pattern of homoplasy that has likely obfuscated previous attempts to reconstruct relationships and delimit species. Important characters presumably related to sexual selection, such as foreleg pigmentation and dense bristle formation, have undergone repeated gain and loss events, many of which have led to increased morphological divergence between sister-species. Evaluation of these traits in a comparative framework illuminates how sexual selection and natural selection influence character evolution and provides a model for future studies of multimodal communication evolution and function.

Micro-Computed Tomography Reveals a Remarkable Twin Intromittent Organ in Spiders – A Novelty for Arachnids With Direct Sperm Transfer

The modification of male pedipalps into secondary sexual intromittent organs is one of the hallmark characteristics of spiders, yet understanding the development and evolution of male genitalia across the order remains a challenging prospect. The embolus – the sclerite bearing the efferent spermatic duct or spermophor, and used to deliver sperm directly to the female genitalia during copulation – has always been considered the single unambiguously homologous palpal sclerite shared by all spider species, fundamental to the bauplan of the order and to the evolution and functional morphology of spider reproductive systems. Indeed, after two centuries of comparative research on spider reproduction, the presence of a single spermophor and embolus on each of a male spider’s two pedipalps remains a central tenet of evolutionary arachnology. Our findings challenge this premise, and reveal a remarkable twin intromittent organ sperm transfer system in a lineage of Australian palpimanoid spiders, characterized by a bifurcate spermophor and the presence of two efferent ducts leading to a pair of embolic sclerites on each pedipalp. This is the first time such a remarkable conformation has been observed in any group of arachnids with direct sperm transfer, complicating our understanding of palpal sclerite homologies, and challenging ideas about the evolution of spider genitalia.

Phylogenomics and loci dropout patterns of deeply diverged Zodarion ant-eating spiders suggest a high potential of RAD-seq for genus-level spider phylogenetics

RAD sequencing yields large amounts of genome-wide data at a relatively low cost and without requiring previous taxon-specific information, making it ideal for evolutionary studies of highly diversified and neglected organisms. However, concerns about information decay with phylogenetic distance have discouraged its use for assessing supraspecific relationships. Here, using Double Digest Restriction Associated DNA (ddRAD) data, we perform the first deep-level approach to the phylogeny of Zodarion, a highly diversified spider genus. We explore the impact of loci and taxon filtering across concatenated and multispecies coalescent reconstruction methods and investigate the patterns of information dropout in reference to both the time of divergence and the mitochondrial divergence between taxa. We found that relaxed loci-filtering and nested taxon-filtering strategies maximized the amount of molecular information and improved phylogenetic inference. As expected, there was a clear pattern of allele dropout towards deeper time and mitochondrial divergences, but the phylogenetic signal remained strong throughout the phylogeny. Therefore, we inferred topologies that were almost fully resolved, highly supported, and noticeably congruent between setups and inference methods, which highlights overall inconsistency in the taxonomy of Zodarion. Because Zodarion appears to be among the oldest and most mitochondrially diversified spider genera, our results suggest that ddRAD data show high potential for inferring intra-generic relationships across spiders and probably also in other taxonomic groups.

A phylogenomic perspective on the evolutionary history of the stonefly genus Suwallia (Plecoptera: Chloroperlidae) revealed by ultraconserved genomic elements

Evolutionary biologists have long sought to disentangle phylogenetic relationships among taxa spanning the tree of life, an increasingly important task as anthropogenic influences accelerate population declines and species extinctions, particularly in insects. Phylogenetic analyses are commonly used to identify unique evolutionary lineages, to clarify taxonomic designations of the focal taxa, and to inform conservation decisions. Advances in DNA sequencing techniques have increasingly facilitated the ability of researchers to apply genomic methods to phylogenetic analyses, even for non-model organisms. Stoneflies are non-model insects that are important bioindicators of the quality of freshwater habitats and landscape disturbance as they spend the immature stages of their life cycles in fresh water, and the adult stages in terrestrial environments. Phylogenetic relationships within the stonefly genus Suwallia (Insecta: Plecoptera: Chloroperlidae) are poorly understood, and have never been assessed using molecular data. We used DNA sequence data from genome-wide ultraconserved element loci to generate the first molecular phylogeny for the group and assess its monophyly. We found that Palearctic and Nearctic Suwallia do not form reciprocally monophyletic clades, and that a biogeographic history including dispersal, vicariance, and founder event speciation via jump dispersal best explains the geographic distribution of this group. Our results also strongly suggest that Neaviperla forcipata (Neave, 1929) is nested within Suwallia, and the concept of the genus Suwallia should be revised to include it. Thus, we formally propose a new taxonomic combination wherein Neaviperla forcipata (Neave, 1929) is reclassified as Suwallia forcipata (Neave, 1929). Moreover, some Suwallia species (e.g., S. amoenacolens, S. kerzhneri, S. marginata, S. pallidula, and S. starki) exhibit pronounced cryptic diversity that is worthy of further investigation. These findings provide a first glimpse into the evolutionary history of Suwallia, improve our understanding of stonefly diversity in the tribe Suwallini, and highlight areas where additional research is needed.

Unlocking inaccessible historical genomes preserved in formalin

Museum specimens represent an unparalleled record of historical genomic data. However, the widespread practice of formalin preservation has thus far impeded genomic analysis of a large proportion of specimens. Limited DNA sequencing from formalin-preserved specimens has yielded low genomic coverage with unpredictable success. We set out to refine sample processing methods and to identify specimen characteristics predictive of sequencing success. With a set of taxonomically diverse specimens collected between 1962 and 2006 and ranging in preservation quality, we compared the efficacy of several end-to-end whole genome sequencing workflows alongside a k-mer-based trimming-free read alignment approach to maximize mapping of endogenous sequence. We recovered complete mitochondrial genomes and up to 3× nuclear genome coverage from formalin-preserved tissues. Hot alkaline lysis coupled with phenol-chloroform extraction out-performed proteinase K digestion in recovering DNA, while library preparation method had little impact on sequencing success. The strongest predictor of DNA yield was overall specimen condition, which additively interacts with preservation conditions to accelerate DNA degradation. Here, we demonstrate a significant advance in capability beyond limited recovery of a small number of loci via PCR or target-capture sequencing. To facilitate strategic selection of suitable specimens for genomic sequencing, we present a decision-making framework that utilizes independent and nondestructive assessment criteria. Sequencing of formalin-preserved specimens will contribute to a greater understanding of temporal trends in genetic adaptation, including those associated with a changing climate. Our work enhances the value of museum collections worldwide by unlocking genomes of specimens that have been disregarded as a valid molecular resource.

Development and evaluation of a custom bait design based on 469 single-copy protein-coding genes for exon capture of isopods (Philosciidae: Haloniscus)

Transcriptome-based exon capture approaches, along with next-generation sequencing, are allowing for the rapid and cost-effective production of extensive and informative phylogenomic datasets from non-model organisms for phylogenetics and population genetics research. These approaches generally employ a reference genome to infer the intron-exon structure of targeted loci and preferentially select longer exons. However, in the absence of an existing and well-annotated genome, we applied this exon capture method directly, without initially identifying intron-exon boundaries for bait design, to a group of highly diverse Haloniscus (Philosciidae), paraplatyarthrid and armadillid isopods, and examined the performance of our methods and bait design for phylogenetic inference. Here, we identified an isopod-specific set of single-copy protein-coding loci, and a custom bait design to capture targeted regions from 469 genes, and analysed the resulting sequence data with a mapping approach and newly-created post-processing scripts. We effectively recovered a large and informative dataset comprising both short (<100 bp) and longer (>300 bp) exons, with high uniformity in sequencing depth. We were also able to successfully capture exon data from up to 16-year-old museum specimens along with more distantly related outgroup taxa, and efficiently pool multiple samples prior to capture. Our well-resolved phylogenies highlight the overall utility of this methodological approach and custom bait design, which offer enormous potential for application to future isopod, as well as broader crustacean, molecular studies.

Converging on the orb: denser taxon sampling elucidates spider phylogeny and new analytical methods support repeated evolution of the orb web

High throughput sequencing and phylogenomic analyses focusing on relationships among spiders have both reinforced and upturned long-standing hypotheses. Likewise, the evolution of spider webs-perhaps their most emblematic attribute-is being understood in new ways. With a matrix including 272 spider species and close arachnid relatives, we analyze and evaluate the relationships among these lineages using a variety of orthology assessment methods, occupancy thresholds, tree inference methods and support metrics. Our analyses include families not previously sampled in transcriptomic analyses, such as Symphytognathidae, the only araneoid family absent in such prior works. We find support for the major established spider lineages, including Mygalomorphae, Araneomorphae, Synspermiata, Palpimanoidea, Araneoidea and the Retrolateral Tibial Apophysis Clade, as well as the uloborids, deinopids, oecobiids and hersiliids Grade. Resulting trees are evaluated using bootstrapping, Shimodaira-Hasegawa approximate likelihood ratio test, local posterior probabilities and concordance factors. Using structured Markov models to assess the evolution of spider webs while accounting for hierarchically nested traits, we find multiple convergent occurrences of the orb web across the spider tree-of-life. Overall, we provide the most comprehensive spider tree-of-life to date using transcriptomic data and use new methods to explore controversial issues of web evolution, including the origins and multiple losses of the orb web.

Mining museums for historical DNA: advances and challenges in museomics

Historical DNA (hDNA), obtained from museum and herbarium specimens, has yielded spectacular new insights into the history of organisms. This includes documenting historical genetic erosion and extinction, discovering species new to science, resolving evolutionary relationships, investigating epigenetic effects, and determining origins of infectious diseases. However, the development of best-practices in isolating, processing, and analyzing hDNA remain under-explored, due to the substantial diversity of specimen preparation types, tissue sources, archival ages, and collecting histories. Thus, for hDNA to reach its full potential, and justify the destructive sampling of the rarest specimens, more experimental work using time-series collections, and the development of improved methods to correct for data asymmetries and biases due to DNA degradation are required.

Genomic library preparation and hybridization capture of formalin-fixed tissues and allozyme supernatant for population genomics and considerations for combining capture- and RADseq-based single nucleotide polymorphism data sets

Until recently many historical museum specimens were largely inaccessible to genomic inquiry, but high-throughput sequencing (HTS) approaches have allowed researchers to successfully sequence genomic DNA from dried and fluid-preserved museum specimens. In addition to preserved specimens, many museums contain large series of allozyme supernatant samples, but the amenability of these samples to HTS has not yet been assessed. Here, we compared the performance of a target-capture approach using alternative sources of genomic DNA from 10 specimens of spring salamanders (Plethodontidae: Gyrinophilus porphyriticus) collected between 1985 and 1990: allozyme supernatants, allozyme homogenate pellets and formalin-fixed tissues. We designed capture probes based on double-digest restriction-site associated sequencing (RADseq) derived loci from frozen blood samples available for seven of the specimens and assessed the success and consistency of capture and RADseq approaches. This study design enabled direct comparisons of data quality and potential biases among the different data sets for phylogenomic and population genomic analyses. We found that in phylogenetic analyses, all enrichment types for a given specimen clustered together. In principal component space all capture-based samples clustered together, but RADseq samples did not cluster with corresponding capture-based samples. Single nucleotide polymorphism calls were on average 18.3% different between enrichment types for a given individual, but these discrepancies were primarily due to differences in heterozygous/homozygous single nucleotide polymorphism calls. We demonstrate that both allozyme supernatant and formalin-fixed samples can be successfully used for population genomic analyses and we discuss ways to identify and reduce biases associated with combining capture and RADseq data.

Museomics dissects the genetic basis for adaptive seasonal colouration in the least weasel

Dissecting the link between genetic variation and adaptive phenotypes provides outstanding opportunities to understand fundamental evolutionary processes. Here, we use a museomics approach to investigate the genetic basis and evolution of winter coat colouration morphs in least weasels (Mustela nivalis), a repeated adaptation for camouflage in mammals with seasonal pelage colour moults across regions with varying winter snow. Whole-genome sequence data was obtained from biological collections and mapped onto a newly assembled reference genome for the species. Sampling represented two replicate transition zones between nivalis and vulgaris colouration morphs in Europe, which typically develop white or brown winter coats, respectively. Population analyses showed that the morph distribution across transition zones is not a by-product of historical structure. Association scans linked a 200 kb genomic region to colouration morph, which was validated by genotyping museum specimens from inter-morph experimental crosses. Genotyping the wild populations narrowed down the association to pigmentation gene MC1R and pinpointed a candidate amino acid change co-segregating with colouration morph. This polymorphism replaces an ancestral leucine residue by lysine at the start of the first extracellular loop of the protein in the vulgaris morph. A selective sweep signature overlapped the association region in vulgaris, suggesting that past adaptation favoured winter-brown morphs and can anchor future adaptive responses to decreasing winter snow. Using biological collections as valuable resources to study natural adaptations, our study showed a new evolutionary route generating winter colour variation in mammals and that seasonal camouflage can be modulated by changes at single key genes.

New specimens from Mid-Cretaceous Myanmar amber illuminate the phylogenetic placement of Lagonomegopidae (Arachnida: Araneae)

New lagonomegopid spiders are described from Mid-Cretaceous Myanmar (Burmese) amber. Two new genera and species based on single specimens, Scopomegops fax gen. & sp. nov. and Hiatomegops spinalis gen. & sp. nov. are described. Two specimens belonging to Lineaburmops beigeli are further described. Additionally, after re-examining the holotype of Odontomegops titan, a detailed description of its basal ventral abdomen is added here. A phylogenetic analysis was performed to investigate the phylogenetic placement of Lagonomegopidae. A matrix of 79 morphological characters, scored for six lagonomegopid taxa and 26 non-lagonomegopid taxa, was analysed through parsimony and Bayesian phylogenetic inference. Our results recover extant Palpimanoidea as a monophyletic group and partly suggest that Lagonomegopidae is the sister-group to extant Palpimanoidea. The external sexual organs, retrolateral tibial apophysis on the male palp and tracheal spiracle in lagonomegopids are discussed.

HyRAD-X Exome Capture Museomics Unravels Giant Ground Beetle Evolution

Advances in phylogenomics contribute toward resolving long-standing evolutionary questions. Notwithstanding, genetic diversity contained within more than a billion biological specimens deposited in natural history museums remains recalcitrant to analysis owing to challenges posed by its intrinsically degraded nature. Yet that tantalizing resource could be critical in overcoming taxon sampling constraints hindering our ability to address major evolutionary questions. We addressed this impediment by developing phyloHyRAD, a new bioinformatic pipeline enabling locus recovery at a broad evolutionary scale from HyRAD-X exome capture of museum specimens of low DNA integrity using a benchtop RAD-derived exome-complexity-reduction probe set developed from high DNA integrity specimens. Our new pipeline can also successfully align raw RNAseq transcriptomic and ultraconserved element reads with the RAD-derived probe catalog. Using this method, we generated a robust timetree for Carabinae beetles, the lack of which had precluded study of macroevolutionary trends pertaining to their biogeography and wing-morphology evolution. We successfully recovered up to 2,945 loci with a mean of 1,788 loci across the exome of specimens of varying age. Coverage was not significantly linked to specimen age, demonstrating the wide exploitability of museum specimens. We also recovered fragmentary mitogenomes compatible with Sanger-sequenced mtDNA. Our phylogenomic timetree revealed a Lower Cretaceous origin for crown group Carabinae, with the extinct Aplothorax Waterhouse, 1841 nested within the genus Calosoma Weber, 1801 demonstrating the junior synonymy of Aplothorax syn. nov., resulting in the new combination Calosoma  burchellii (Waterhouse, 1841) comb. nov. This study compellingly illustrates that HyRAD-X and phyloHyRAD efficiently provide genomic-level data sets informative at deep evolutionary scales.

Museomics and phylogenomics with protein-encoding ultraconserved elements illuminate the evolution of life history and phallic morphology of flesh flies (Diptera: Sarcophagidae)

BACKGROUND

The common name of the Flesh flies (Sarcophagidae) usually relates them with organisms feeding on decomposing organic matter, although the biology of one of the largest radiations among insects also includes predation, coprophagy, and even kleptoparasitism. The question of whether the ancestor of all sarcophagids was a predator or a decomposer, or in association to which host have sarcophagids evolved, has thus always piqued the curiosity of flesh fly specialists. Such curiosity has often been hindered by both the impossibility of having a well-supported phylogeny of Sarcophagidae and its sister group to trace live habits and the scarcity of information on the biology of the group. Using a phylogenomic dataset of protein-encoding ultraconserved elements from representatives of all three subfamilies of Sarcophagidae as ingroup and a large Calyptratae outgroup, a robust phylogenetic framework and timescale are generated to understand flesh fly systematics and the evolution of their life histories.

RESULTS

The evolutionary history for Sarcophagidae reconstructed here differs considerably from previous hypotheses. Within subfamily Sarcophaginae, a group of predatory flies, including genera Lepidodexia and Boettcheria, emerged as sister-group to the rest of Sarcophaginae. The genera Oxysarcodexia, Ravinia, and Tricharaea, long considered archaic and early-branching coprophagous and sarcosaprophagous lineages, were found nested well within the Sarcophaginae as sister-group to the sarcosaprophagous Microcerella. Predation on invertebrates is suggested as the ancestral and dominant strategy throughout the early evolution of flesh flies. Several transitions from predation to sarcosaprophagy and coprophagy occur across the sarcophagid phylogenetic tree, in contrast with almost no transitions from sarcosaprophagy or coprophagy to predatory habits. Regarding the morphological evolution of flesh flies, there might be a concerted evolution of male genitalia traits, such as the phallotrema position and the juxta, or the vesica and the folding of the phallotrema. One diversification rate shift was inferred in the evolution of sarcophagids, which is related to the origin of genus Sarcophaga.

CONCLUSIONS

This study has a significant impact on understanding sarcophagid evolution and highlights the importance of having a robust phylogenetic framework to reconstruct the ancestral character state of biological and morphological characters. I discuss the evolution of life histories of the family in relation to their hosts or substrates and outline how sarcosaprophagy, coprophagy, and kleptoparasitism behavior on various hosts may have evolved from predation on invertebrates. This study provides a phylogenetic framework for further physiological and comparative genomic work between predatory, sarcosaprophagous, coprophagous, and kleptoparasitic lineages, which could also have significant implications for the evolution of diverse life histories in other Diptera.

Interrogating Genomic-Scale Data to Resolve Recalcitrant Nodes in the Spider Tree of Life

Genome-scale data sets are converging on robust, stable phylogenetic hypotheses for many lineages; however, some nodes have shown disagreement across classes of data. We use spiders (Araneae) as a system to identify the causes of incongruence in phylogenetic signal between three classes of data: exons (as in phylotranscriptomics), noncoding regions (included in ultraconserved elements [UCE] analyses), and a combination of both (as in UCE analyses). Gene orthologs, coded as amino acids and nucleotides (with and without third codon positions), were generated by querying published transcriptomes for UCEs, recovering 1,931 UCE loci (codingUCEs). We expected that congeners represented in the codingUCE and UCEs data would form clades in the presence of phylogenetic signal. Noncoding regions derived from UCE sequences were recovered to test the stability of relationships. Phylogenetic relationships resulting from all analyses were largely congruent. All nucleotide data sets from transcriptomes, UCEs, or a combination of both recovered similar topologies in contrast with results from transcriptomes analyzed as amino acids. Most relationships inferred from low-occupancy data sets, containing several hundreds of loci, were congruent across Araneae, as opposed to high occupancy data matrices with fewer loci, which showed more variation. Furthermore, we found that low-occupancy data sets analyzed as nucleotides (as is typical of UCE data sets) can result in more congruent relationships than high occupancy data sets analyzed as amino acids (as in phylotranscriptomics). Thus, omitting data, through amino acid translation or via retention of only high occupancy loci, may have a deleterious effect in phylogenetic reconstruction.

Using historical genome-wide DNA to unravel the confused taxonomy in a songbird lineage that is extinct in the wild

Urgent conservation action for terminally endangered species is sometimes hampered by taxonomic uncertainty, especially in illegally traded animals that are often cross-bred in captivity. To overcome these problems, we used a genomic approach to analyze historical DNA from museum samples across the Asian Pied Starling (Gracupica contra) complex in tropical Asia, a popular victim of the ongoing songbird crisis whose distinct Javan population ("Javan Pied Starling") is extinct in the wild and subject to admixture in captivity. Comparing genomic profiles across the entire distribution, we detected three deeply diverged lineages at the species level characterized by a lack of genomic intermediacy near areas of contact. Our study demonstrates that the use of historical DNA can be instrumental in delimiting species in situations of taxonomic uncertainty, especially when modern admixture may obfuscate species boundaries. Results of our research will enable conservationists to commence a dedicated ex situ breeding program for the Javan Pied Starling, and serve as a blueprint for similar conservation problems involving terminally endangered species subject to allelic infiltration from close congeners.

Phylogenomic Resolution of Sea Spider Diversification through Integration of Multiple Data Classes

Despite significant advances in invertebrate phylogenomics over the past decade, the higher-level phylogeny of Pycnogonida (sea spiders) remains elusive. Due to the inaccessibility of some small-bodied lineages, few phylogenetic studies have sampled all sea spider families. Previous efforts based on a handful of genes have yielded unstable tree topologies. Here, we inferred the relationships of 89 sea spider species using targeted capture of the mitochondrial genome, 56 conserved exons, 101 ultraconserved elements, and 3 nuclear ribosomal genes. We inferred molecular divergence times by integrating morphological data for fossil species to calibrate 15 nodes in the arthropod tree of life. This integration of data classes resolved the basal topology of sea spiders with high support. The enigmatic family Austrodecidae was resolved as the sister group to the remaining Pycnogonida and the small-bodied family Rhynchothoracidae as the sister group of the robust-bodied family Pycnogonidae. Molecular divergence time estimation recovered a basal divergence of crown group sea spiders in the Ordovician. Comparison of diversification dynamics with other marine invertebrate taxa that originated in the Paleozoic suggests that sea spiders and some crustacean groups exhibit resilience to mass extinction episodes, relative to mollusk and echinoderm lineages.

Phylogenetic Systematics and Evolution of the Spider Infraorder Mygalomorphae Using Genomic Scale Data

The infraorder Mygalomorphae is one of the three main lineages of spiders comprising over 3000 nominal species. This ancient group has a worldwide distribution that includes among its ranks large and charismatic taxa such as tarantulas, trapdoor spiders, and highly venomous funnel-web spiders. Based on past molecular studies using Sanger-sequencing approaches, numerous mygalomorph families (e.g., Hexathelidae, Ctenizidae, Cyrtaucheniidae, Dipluridae, and Nemesiidae) have been identified as non-monophyletic. However, these data were unable to sufficiently resolve the higher-level (intra- and interfamilial) relationships such that the necessary changes in classification could be made with confidence. Here, we present a comprehensive phylogenomic treatment of the spider infraorder Mygalomorphae. We employ 472 loci obtained through anchored hybrid enrichment to reconstruct relationships among all the mygalomorph spider families and estimate the timeframe of their diversification. We sampled nearly all currently recognized families, which has allowed us to assess their status, and as a result, propose a new classification scheme. Our generic-level sampling has also provided an evolutionary framework for revisiting questions regarding silk use in mygalomorph spiders. The first such analysis for the group within a strict phylogenetic framework shows that a sheet web is likely the plesiomorphic condition for mygalomorphs, as well as providing insights to the ancestral foraging behavior for all spiders. Our divergence time estimates, concomitant with detailed biogeographic analysis, suggest that both ancient continental-level vicariance and more recent dispersal events have played an important role in shaping modern day distributional patterns. Based on our results, we relimit the generic composition of the Ctenizidae, Cyrtaucheniidae, Dipluridae, and Nemesiidae. We also elevate five subfamilies to family rank: Anamidae (NEW RANK), Euagridae (NEW RANK), Ischnothelidae (NEW RANK), Pycnothelidae (NEW RANK), and Bemmeridae (NEW RANK). Three families Entypesidae (NEW FAMILY), Microhexuridae (NEW FAMILY), and Stasimopidae (NEW FAMILY), and one subfamily Australothelinae (NEW SUBFAMILY) are newly proposed. Such a major rearrangement in classification, recognizing nine newly established family-level rank taxa, is the largest the group has seen in over three decades. [Biogeography; molecular clocks; phylogenomics; spider web foraging; taxonomy.].

Phylogenomics and biogeography of leptonetid spiders (Araneae:Leptonetidae)

Leptonetidae are rarely encountered spiders, usually associated with caves and mesic habitats, and are disjunctly distributed across the Holarctic. Data from ultraconserved elements (UCEs) were used in concatenated and coalescent-based analyses to estimate the phylogenetic history of the family. Our taxon sample included close outgroups, and 90% of described leptonetid genera, with denser sampling in North America and Mediterranean Europe. Two data matrices were assembled and analysed; the first ‘relaxed’ matrix includes the maximum number of loci and the second ‘strict’ matrix is limited to the same set of core orthologs but with flanking introns mostly removed. A molecular dating analysis incorporating fossil and geological calibration points was used to estimate divergence times, and dispersal–extinction–cladogenesis analysis (DEC) was used to infer ancestral distributions. Analysis of both data matrices using maximum likelihood and coalescent-based methods supports the monophyly of Archoleptonetinae and Leptonetinae. However, relationships among Archoleptonetinae, Leptonetinae, and Austrochiloidea are poorly supported and remain unresolved. Archoleptonetinae is elevated to family rank Archoleptonetidae (new rank) and Leptonetidae (new status) is restricted to include only members of the subfamily Leptonetinae; a taxonomic review with morphological diagnoses is provided for both families. Four well supported lineages within Leptonetidae (new status) are recovered: (1) the Calileptoneta group, (2) the Leptoneta group, (3) the Paraleptoneta group, and (4) the Protoleptoneta group. Most genera within Leptonetidae are monophyletic, although Barusia, Cataleptoneta, and Leptoneta include misplaced species and require taxonomic revision. The origin of Archoleptonetidae (new rank), Leptonetidae, and the four main lineages within Leptonetidae date to the Cretaceous. DEC analysis infers the Leptoneta and Paraleptoneta groups to have ancestral distributions restricted to Mediterranean Europe, whereas the Calileptoneta and Protoleptoneta groups include genera with ancestral distributions spanning eastern and western North America, Mediterranean Europe, and east Asia. Based on a combination of biology, estimated divergence times, and inferred ancestral distributions we hypothesise that Leptonetidae was once widespread across the Holarctic and their present distributions are largely the result of vicariance. Given the wide disjunctions between taxa, we broadly interpret the family as a Holarctic relict fauna and hypothesise that they were once part of the Boreotropical forest ecosystem.

New Light on Historical Specimens Reveals a New Species of Ladybird (Coleoptera: Coccinellidae): Morphological, Museomic, and Phylogenetic Analyses

Natural history collections house an important source of genetic data from yet unexplored biological diversity. Molecular data from museum specimens remain underexploited, which is mainly due to the degradation of DNA from specimens over time. However, Next-Generation Sequencing (NGS) technology can now be used to sequence "old" specimens. Indeed, many of these specimens are unique samples of nomenclatural types and can be crucial for resolving systematic or biogeographic scientific questions. Two ladybird beetle specimens from Patagonia corresponding to a new species of the genus Eriopis Mulsant were found in the collections of the Muséum national d'Histoire naturelle (MNHN), Paris. Here, we describe Eriopis patagonia Salazar, sp. nov. Total DNA of one of the two specimens was sequenced by NGS using a paired-end Illumina approach. We reconstruct and characterize the mitochondrial genome of this species (16,194 bp). Then, the protein-coding genes (PCGs) and ribosomal RNAs (rRNAs) were used to infer by maximum likelihood and Bayesian Inference the phylogenetic position of E. patagonia among 27 representatives of Coccinellidae. Phylogenetic analysis confirmed the position of Eriopis as sister group to Cycloneda Crotch. Hence, we highlight the high potential of sequencing technology for extracting molecular information from old specimens, which are used here for the systematic study of a genus, while demonstrating the importance of preserving biological collections.

Morphology and performance of the 'trap-jaw' cheliceral strikes in spiders (Araneae, Mecysmaucheniidae)

Mecysmaucheniidae spiders have evolved ultra-fast cheliceral strikes 4 times independently. The mechanism for producing these high-speed strikes is likely due to a latch/spring system that allows for stored energy to be rapidly released. This study examined two different sister lineages: Zearchaea has ultra-fast cheliceral strikes and Aotearoa, based on external morphology of the clypeus, is hypothesized to have slower strikes. Using high-speed videography, I first gathered kinematic data on each taxon. Then, using histology and data from micro-computed tomography scanning, I examined internal cheliceral muscle morphology to test whether shifts in muscle anatomy correspond to performance differences in cheliceral strike. Results from high-speed video analysis revealed that Zearchaea achieves peak angular velocities of 25.0×103±4.8×103 rad s-1 (mean±s.d.) in durations of 0.0843±0.017 ms. The fastest recorded strike had a peak angular and linear velocity of 30.8×103 rad s-1 and 18.2 m s-1, respectively. The slower striking sister species, Aotearoa magna, was three orders of magnitude slower in velocity and longer in duration. Histology revealed sarcomere length differences, with some muscles optimized for force, and other muscles for speed. 3D printed models revealed structural differences that explain how the chelicerae hinge open and close. Combining all of this evidence, I put forth a hypothesis for the ultra-fast trap-jaw mechanism. This research documents the morphological shifts that accompany ultra-fast movements that result in increased rotation in joints and increased muscle specialization.

Advantages of an easy-to-use DNA extraction method for minimal-destructive analysis of collection specimens

Here we present and justify an approach for minimal-destructive DNA extraction from historic insect specimens for next generation sequencing applications. An increasing number of studies use insects from museum collections for biodiversity research. However, the availability of specimens for molecular analyses has been limited by the degraded nature of the DNA gained from century-old museum material and the consumptive nature of most DNA extraction procedures. The method described in this manuscript enabled us to successfully extract DNA from specimens as old as 241 years using a minimal-destructive approach. The direct comparison of the DNeasy extraction Kit and the Monarch® PCR & DNA Clean-up Kit showed a significant increase of 17.3-fold higher DNA yield extracted with the Monarch Oligo protocol on average. By using an extraction protocol originally designed for oligonucleotide clean-up, we were able to combine overcoming the restrictions by target fragment size and strand state, with minimising time consumption and labour-intensity. The type specimens used for the minimal-destructive DNA extraction exhibited no significant external change or post-extraction damage, while sufficient DNA was retrieved for analyses.

Sequence Capture Phylogenomics of True Spiders Reveals Convergent Evolution of Respiratory Systems

The common ancestor of spiders likely used silk to line burrows or make simple webs, with specialized spinning organs and aerial webs originating with the evolution of the megadiverse "true spiders" (Araneomorphae). The base of the araneomorph tree also concentrates the greatest number of changes in respiratory structures, a character system whose evolution is still poorly understood, and that might be related to the evolution of silk glands. Emphasizing a dense sampling of multiple araneomorph lineages where tracheal systems likely originated, we gathered genomic-scale data and reconstructed a phylogeny of true spiders. This robust phylogenomic framework was used to conduct maximum likelihood and Bayesian character evolution analyses for respiratory systems, silk glands, and aerial webs, based on a combination of original and published data. Our results indicate that in true spiders, posterior book lungs were transformed into morphologically similar tracheal systems six times independently, after the evolution of novel silk gland systems and the origin of aerial webs. From these comparative data, we put forth a novel hypothesis that early-diverging web-building spiders were faced with new energetic demands for spinning, which prompted the evolution of similar tracheal systems via convergence; we also propose tests of predictions derived from this hypothesis.[Book lungs; discrete character evolution; respiratory systems; silk; spider web evolution; ultraconserved elements.].

A linked-read approach to museomics: Higher quality de novo genome assemblies from degraded tissues

High-throughput sequencing technologies are a proposed solution for accessing the molecular data in historical specimens. However, degraded DNA combined with the computational demands of short-read assemblies has posed significant laboratory and bioinformatics challenges for de novo genome assembly. Linked-read or "synthetic long-read" sequencing technologies, such as 10× Genomics, may provide a cost-effective alternative solution to assemble higher quality de novo genomes from degraded tissue samples. Here, we compare assembly quality (e.g., genome contiguity and completeness, presence of orthogroups) between four new deer mouse (Peromyscus spp.) genomes assembled using linked-read technology and four published genomes assembled from a single shotgun library. At a similar price-point, these approaches produce vastly different assemblies, with linked-read assemblies having overall higher contiguity and completeness, measured by larger N50 values and greater number of genes assembled, respectively. As a proof-of-concept, we used annotated genes from the four Peromyscus linked-read assemblies and eight additional rodent taxa to generate a phylogeny, which reconstructed the expected relationships among species with 100% support. Although not without caveats, our results suggest that linked-read sequencing approaches are a viable option to build de novo genomes from degraded tissues, which may prove particularly valuable for taxa that are extinct, rare or difficult to collect.

Analyzing drivers of speciation in the Southern Ocean using the sea spider species complex Colossendeis megalonyx as a test case

Colossendeis megalonyx Hoek, 1881 has the broadest distribution of all sea spiders in the Southern Ocean. Previous studies have detected several evolutionarily young lineages within this taxon and interpreted them as a result of allopatric speciation in a few shelf refuges during glacial maxima. However, alternative scenarios such as ecological speciation in sympatry have rarely been considered or tested. Here, we generated the most extensive genomic and morphometric data set on the C. megalonyx species complex to (i) comprehensively describe species diversity, (ii) explore intraspecific connectivity between populations located around Antarctica, and (iii) systematically test for positive selection indicative of adaptive speciation. We successfully applied a target hybrid enrichment approach and recovered all 1607 genes targeted. Phylogenomic analysis was consistent with previous findings and, moreover, increased the resolution of branching within lineages. We found specimens of phylogenetically well-separated lineages occurring in sympatry to be genetically distinct from each other and gene flow between geographically separated populations of the same lineages to be restricted. Evidence for positive selection was found for four genes associated with structural and neuronal functions. Hence, there is an indication for positive selection in the C. megalonyx species complex, yet its specific contribution to the speciation process remains to be explored further. Finally, morphometric analyses revealed multiple significant differences between lineages, but a clear separation proved difficult. Our study highlights the relevance of positive selection as a potential driver for speciation in the Southern Ocean.

Effects of sample age on data quality from targeted sequencing of museum specimens: what are we capturing in time?

Background

Next generation sequencing (NGS) can recover DNA data from valuable extant and extinct museum specimens. However, archived or preserved DNA is difficult to sequence because of its fragmented, damaged nature, such that the most successful NGS methods for preserved specimens remain sub-optimal. Improving wet-lab protocols and comprehensively determining the effects of sample age on NGS library quality are therefore of vital importance. Here, I examine the relationship between sample age and several indicators of library quality following targeted NGS sequencing of ~ 1300 loci using 271 samples of pinned moth specimens (Helicoverpa armigera) ranging in age from 5 to 117 years.

Results

I find that older samples have lower DNA concentrations following extraction and thus require a higher number of indexing PCR cycles during library preparation. When sequenced reads are aligned to a reference genome or to only the targeted region, older samples have a lower number of sequenced and mapped reads, lower mean coverage, and lower estimated library sizes, while the percentage of adapters in sequenced reads increases significantly as samples become older. Older samples also show the poorest capture success, with lower enrichment and a higher improved coverage anticipated from further sequencing.

Conclusions

Sample age has significant, measurable impacts on the quality of NGS data following targeted enrichment. However, incorporating a uracil-removing enzyme into the blunt end-repair step during library preparation could help to repair DNA damage, and using a method that prevents adapter-dimer formation may result in improved data yields.

Spider-specific probe set for ultraconserved elements offers new perspectives on the evolutionary history of spiders (Arachnida, Araneae)

Phylogenomic methods have proven useful for resolving deep nodes and recalcitrant groups in the spider tree of life. Across arachnids, transcriptomic approaches may generate thousands of loci, and target-capture methods, using the previously designed arachnid-specific probe set, can target a maximum of about 1,000 loci. Here, we develop a specialized target-capture probe set for spiders that contains over 2,000 ultraconserved elements (UCEs) and then demonstrate the utility of this probe set through sequencing and phylogenetic analysis. We designed the 'spider-specific' probe set using three spider genomes (Loxosceles, Parasteatoda and Stegodyphus) and ensured that the newly designed probe set includes UCEs from the previously designed Arachnida probe set. The new 'spider-specific' probes were used to sequence UCE loci in 51 specimens. The remaining samples included five spider genomes and taxa that were enriched using Arachnida probe set. The 'spider-specific' probes were also used to gather loci from a total of 84 representative taxa across Araneae. On mapping these 84 taxa to the Arachnida probe set, we captured at most 710 UCE loci, while the spider-specific probe set captured up to 1,547 UCE loci from the same taxon sample. Phylogenetic analyses using maximum likelihood and coalescent methods corroborate most nodes resolved by recent transcriptomic analyses, but not all (e.g. UCE data suggest monophyly of 'symphytognathoids'). Our preferred hypothesis based on topology tests, suggests monophyly of the 'symphytognathoids' (the miniature orb weavers), which in previous studies has only been supported by a combination of morphological and behavioural characters.

More than skin and bones: Comparing extraction methods and alternative sources of DNA from avian museum specimens

Next-generation sequencing has greatly expanded the utility and value of museum collections by revealing specimens as genomic resources. As the field of museum genomics grows, so does the need for extraction methods that maximize DNA yields. For avian museum specimens, the established method of extracting DNA from toe pads works well for most specimens. However, for some specimens, especially those of birds that are very small or very large, toe pads can be a poor source of DNA. In this study, we apply two DNA extraction methods (phenol-chloroform and silica column) to three different sources of DNA (toe pad, skin punch and bone) from 10 historical avian museum specimens. We show that a modified phenol-chloroform protocol yielded significantly more DNA than a silica column protocol (e.g., Qiagen DNeasy Blood & Tissue Kit) across all tissue types. However, extractions using the silica column protocol contained longer fragments on average than those using the phenol-chloroform protocol, probably as a result of loss of small fragments through the silica column. While toe pads yielded more DNA than skin punches and bone fragments, skin punches proved to be a reliable alternative source of DNA and might be especially appealing when toe pad extractions are impractical. Overall, we found that historical bird museum specimens contain substantial amounts of DNA for genomic studies under most extraction scenarios, but that a phenol-chloroform protocol consistently provides the high quantities of DNA required for most current genomic protocols.

Golden Orbweavers Ignore Biological Rules: Phylogenomic and Comparative Analyses Unravel a Complex Evolution of Sexual Size Dimorphism

Instances of sexual size dimorphism (SSD) provide the context for rigorous tests of biological rules of size evolution, such as Cope's rule (phyletic size increase), Rensch's rule (allometric patterns of male and female size), as well as male and female body size optima. In certain spider groups, such as the golden orbweavers (Nephilidae), extreme female-biased SSD (eSSD, female:male body length $\ge$2) is the norm. Nephilid genera construct webs of exaggerated proportions, which can be aerial, arboricolous, or intermediate (hybrid). First, we established the backbone phylogeny of Nephilidae using 367 anchored hybrid enrichment markers, then combined these data with classical markers for a reference species-level phylogeny. Second, we used the phylogeny to test Cope and Rensch's rules, sex specific size optima, and the coevolution of web size, type, and features with female and male body size and their ratio, SSD. Male, but not female, size increases significantly over time, and refutes Cope's rule. Allometric analyses reject the converse, Rensch's rule. Male and female body sizes are uncorrelated. Female size evolution is random, but males evolve toward an optimum size (3.2-4.9 mm). Overall, female body size correlates positively with absolute web size. However, intermediate sized females build the largest webs (of the hybrid type), giant female Nephila and Trichonephila build smaller webs (of the aerial type), and the smallest females build the smallest webs (of the arboricolous type). We propose taxonomic changes based on the criteria of clade age, monophyly and exclusivity, classification information content, and diagnosability. Spider families, as currently defined, tend to be between 37 million years old and 98 million years old, and Nephilidae is estimated at 133 Ma (97-146), thus deserving family status. We, therefore, resurrect the family Nephilidae Simon 1894 that contains Clitaetra Simon 1889, the Cretaceous GeratonephilaPoinar and Buckley (2012), Herennia Thorell 1877, IndoetraKuntner 2006, new rank, Nephila Leach 1815, Nephilengys L. Koch 1872, Nephilingis Kuntner 2013, Palaeonephila Wunderlich 2004 from Tertiary Baltic amber, and TrichonephilaDahl 1911, new rank. We propose the new clade Orbipurae to contain Araneidae Clerck 1757, Phonognathidae Simon 1894, new rank, and Nephilidae. Nephilid female gigantism is a phylogenetically ancient phenotype (over 100 Ma), as is eSSD, though their magnitudes vary by lineage.

Phylogenomic analysis and revised classification of atypoid mygalomorph spiders (Araneae, Mygalomorphae), with notes on arachnid ultraconserved element loci

The atypoid mygalomorphs include spiders from three described families that build a diverse array of entrance web constructs, including funnel-and-sheet webs, purse webs, trapdoors, turrets and silken collars. Molecular phylogenetic analyses have generally supported the monophyly of Atypoidea, but prior studies have not sampled all relevant taxa. Here we generated a dataset of ultraconserved element loci for all described atypoid genera, including taxa (Mecicobothrium and Hexurella) key to understanding familial monophyly, divergence times, and patterns of entrance web evolution. We show that the conserved regions of the arachnid UCE probe set target exons, such that it should be possible to combine UCE and transcriptome datasets in arachnids. We also show that different UCE probes sometimes target the same protein, and under the matching parameters used here show that UCE alignments sometimes include non-orthologs. Using multiple curated phylogenomic matrices we recover a monophyletic Atypoidea, and reveal that the family Mecicobothriidae comprises four separate and divergent lineages. Fossil-calibrated divergence time analyses suggest ancient Triassic (or older) origins for several relictual atypoid lineages, with late Cretaceous/early Tertiary divergences within some genera indicating a high potential for cryptic species diversity. The ancestral entrance web construct for atypoids, and all mygalomorphs, is reconstructed as a funnel-and-sheet web.

Comparative morphology of cheliceral muscles using high-resolution X-ray microcomputed-tomography in palpimanoid spiders (Araneae, Palpimanoidea)

Spiders are important predators in terrestrial ecosystems, yet we know very little about the principal feeding structures of spiders, the chelicerae, which are functionally equivalent to "jaws" or "mandibles" and are an extremely important aspect of spider biology. In particular, members of Palpimanoidea have evolved highly unusual cheliceral morphologies and functions, including high-speed, ballistic movements in mecysmaucheniid spiders, the fastest arachnid movements known thus far, and the elongated, highly maneuverable chelicerae of archaeids that use an attack-at-a-distance strategy. Here, using micro-Computed-Tomography scanning techniques, we perform a comparative study to examine cheliceral muscle morphology in six different spider specimens representing five palpimanoid families. We provide a hypothesis for homology in palpimanoid cheliceral muscles and then compare and contrast these findings with previous studies on other non-palpimanoid spiders. We document and discuss two sets of cheliceral muscles in palpimanoids that have not been previously observed in other spiders or which may represent a position shift compared to other spiders. In the palpimanoids, Palpimanus sp., Huttonia sp., and Colopea sp. showed similar cheliceral muscle anatomy. In Eriauchenius ranavalona, which has highly maneuverable chelicerae, some of the muscles have a more horizontal orientation, and there is a greater degree of cheliceral muscle divergence. In Zearchaea sp. and Aotearoa magna, some muscles have also shifted to a more horizontal orientation, and in Zearchaea sp., a species with a ballistic, high-speed predatory strike, there is a loss of cheliceral muscles. This research is a first step toward understanding cheliceral form and function across spiders.