Reassessing the phylogeny of Cyphophthalmi with phylogenomics: A UCE-based phylogeny of mite harvesters (Opiliones)

Cyphophthalmi (the mite harvesters) are a group of Opiliones with broad interest due to their species being classic examples of short-range endemics and displaying model biogeographical patterns for poor dispersers. Cyphophthalmi phylogeny has received attention using morphology, Sanger-based sequencing data, or transcriptomics. Here we turn to a new type of data, ultraconserved elements (UCEs) and provide a first phylogeny for the entire suborder Cyphophthalmi using such data and including representatives from 36 of the 46 currently recognized genera. Phylogenetic analysis of four occupancy matrices (50%, 75%, 90% and 95%), for a total of 840, 567, 129, and 23 loci, respectively, yielded a well resolved phylogeny with monophyly of Pettalidae, Parasironidae, Stylocellidae and Troglosironidae. However, Neogoveidae appeared paraphyletic with respect to Ogoveidae in all datasets and to Troglosironidae in some, and the traditional Sironidae, which was monophyletic, now appeared paraphyletic with respect to the recently erected family Parasironidae. Our phylogenomic results using UCE data resolve the position of several problematic genera (e.g., Pettalus) and add support to other parts of the tree that received low support in Sanger-based phylogenies. Our work also stresses the possibility to add museum samples to phylogenies although methods for optimizing DNA yield from such small-bodied specimens need further improvement. Finally, this backbone phylogeny demonstrates the feasibility of an all-species phylogeny using UCEs for Cyphophthalmi, and by extension, for all Opiliones.

Testing ultraconserved elements (UCEs) for phylogenetic inference across bivalves (Mollusca: Bivalvia)

Bivalves constitute an important resource for fisheries and as cultural objects. Bivalve phylogenetics has had a long tradition using both morphological and molecular characters, and genomic resources are available for a good number of commercially important species. However, relationships among bivalve families have been unstable and major conflicting results exist between mitogenomics and results based on Sanger-based amplicon sequencing or phylotranscriptomics. Here we design and test an ultraconserved elements probe set for the class Bivalvia with the aim to use hundreds of loci without the need to sequence full genomes or transcriptomes, which are expensive and complex to analyze, and to open bivalve phylogenetics to museum specimens. Our probe set successfully captured 1,513 UCEs for a total of 263,800 bp with an average length of 174.59 ± 3.44 per UCE (ranging from 28 to 842 bp). Phylogenetic testing of this UCE probe set across Bivalvia and within the family Donacidae using different data matrices and methods for phylogenetic inference shows promising results at multiple taxonomic levels. In addition, our probe set was able to capture large numbers of UCEs for museum specimens collected before 1900 and from DNAs properly stored, of which many museums and laboratories are well stocked. Overall, this constitutes a novel and useful resource for bivalve phylogenetics.

An ultraconserved element probe set for velvet worms (Onychophora)

Onychophora are cryptic, soil-dwelling invertebrates known for their biogeographic affinities, diversity of reproductive modes, close phylogenetic relationship to arthropods, and peculiar prey capture mechanism. The 216 valid species of Onychophora are grouped into two families - Peripatopsidae and Peripatidae - and apart from a few relationships among major lineages within these two families, a stable phylogenetic backbone for the phylum has yet to be resolved. This has hindered our understanding of onychophoran biogeographic patterns, evolutionary history, and systematics. Neopatida, the Neotropical clade of peripatids, has proved particularly difficult, with recalcitrant nodes and low resolution, potentially due to rapid radiation of the group during the Cretaceous. Previous studies have had to compromise between number of loci and number of taxa due to limitations of Sanger sequencing and phylotranscriptomics, respectively. Additionally, aspects of their genome size and structure have made molecular phylogenetics difficult and data matrices have been affected by missing data. To address these issues, we leveraged recent, published transcriptomes and the first high quality genome for the phylum and designed a high affinity ultraconserved element (UCE) probe set for Onychophora. This new probe set, consisting of ∼ 20,000 probes that target 1,465 loci across both families, has high locus recovery and phylogenetic utility. Phylogenetic analyses recovered the monophyly of major clades of Onychophora and revealed a novel lineage from the Neotropics that challenges our current understanding of onychophoran biogeographic endemicity. This new resource could drastically increase the power of molecular datasets and potentially allow access to genomic scale data from archival museum specimens to further tackle the issues exasperating onychophoran systematics.

Machine learning approaches delimit cryptic taxa in a previously intractable species complex

Cryptic species are not diagnosable via morphological criteria, but can be detected through analysis of DNA sequences. A number of methods have been developed for identifying species based on genetic data; however, these methods are prone to over-splitting taxa with extreme population structure, such as dispersal-limited organisms. Machine learning methodologies have the potential to overcome this challenge. Here, we apply such approaches, using a large dataset generated through hybrid target enrichment of ultraconserved elements (UCEs). Our study taxon is the Aoraki denticulata species complex, a lineage of extremely low-dispersal arachnids endemic to the South Island of Aotearoa New Zealand. This group of mite harvesters has been the subject of previous species delimitation studies using smaller datasets generated through Sanger sequencing and analytical approaches that rely on multispecies coalescent models and barcoding gap discovery. Those analyses yielded a number of putative cryptic species that seems unrealistic and extreme, based on what we know about species' geographic ranges and genetic diversity in non-cryptic mite harvesters. We find that machine learning approaches, on the other hand, identify cryptic species with geographic ranges that are similar to those seen in other morphologically diagnosable mite harvesters in Aotearoa New Zealand's South Island. We performed both unsupervised and supervised machine learning analyses, the latter with training data drawn either from animals broadly (vagile and non-vagile) or from a custom training dataset from dispersal-limited harvesters. We conclude that applying machine learning approaches to the analysis of UCE-derived genetic data is an effective method for delimiting species in complexes of low-vagility cryptic species, and that the incorporation of training data from biologically relevant analogues can be critically informative.

Under the hood: Phylogenomics of hooded tick spiders (Arachnida, Ricinulei) uncovers discordance between morphology and molecules

Ricinulei or hooded tick-spiders are a cryptic and ancient group of arachnids. The order consists of around 100 highly endemic extant species restricted to the Afrotropics and the Neotropics along with 22 fossil species. Their antiquity and low vagility make them an excellent group with which to interrogate biogeographic questions. To date, only four molecular analyses have been conducted on the group and they failed to resolve the relationships of the main lineages and even recovering the non-monophyly of the three genera. These studies were limited to a few Sanger loci or phylogenomic analyses with at most seven ingroup samples. To increase phylogenetic resolution in this little-understood and poorly studied group, we present the most comprehensive phylogenomic study of Ricinulei to date leveraging the Arachnida ultra-conserved element probe set. With a data set of 473 loci across 96 ingroup samples, analyses resolved a monophyletic Neotropical clade consisting of four main lineages. Two of them correspond to the current genera Cryptocellus and Pseudocellus while topology testing revealed one lineage to likely be a phylogenetic reconstruction artefact. The fourth lineage, restricted to Northwestern, Andean South America, is consistent with the Cryptocellus magnus group, likely corresponding to the historical genus Heteroricinoides. Since we did not sample the type species for this old genus, we do not formally re-erect Heteroricinoides but our data suggest the need for a thorough morphological re-examination of Neotropical Ricinulei.

300 million years apart: the extreme case of macromorphological skeletal convergence between deltocyathids and a turbinoliid coral (Anthozoa, Scleractinia)

The integration of morphological and molecular lines of evidence has enabled the family Deltocyathidae to be erected to accommodate Deltocyathus species that were previously ascribed to the family Caryophylliidae. However, although displaying the same morphological characteristics as other species of Deltocyathus , molecular data suggested that D. magnificus was phylogenetically distant from Deltocyathidae, falling within the family Turbinoliidae instead. To elucidate the enigmatic evolutionary history of this species and skeletal microstructural features, the phylogenetic relationships of Deltocyathidae and Turbinoliidae were investigated using nuclear ultraconserved and exon loci and complete mitochondrial genomes. Both nuclear and mitochondrial phylogenomic reconstructions confirmed the position of D. magnificus within turbinolids. Furthermore, a novel mitochondrial gene order was uncovered for Deltocyathidae species. This gene order was not present in Turbinoliidae or in D. magnificus that both have the scleractinian canonical gene order, further indicating the taxonomic utility of mitochondrial gene order. D. magnificus is therefore formally moved to the family Turbinoliidae and accommodated in a new genus (Dennantotrochus Kitahara, Vaga & Stolarski, gen. nov.). Surprisingly, turbinolids and deltocyathids do not differ in microstructural organisation of the skeleton that consists of densely packed, individualised rapid accretion deposits and thickening deposits composed of fibres perpendicular to the skeleton surface. Therefore, although both families are clearly evolutionarily divergent, macromorphological features indicate a case of skeletal convergence while these may still share conservative biomineralisation mechanisms. ZooBank: urn:lsid:zoobank.org:pub:5F1C0E25-3CC6-4D1F-B1F0-CD9D0014678E.

The effect of missing data on evolutionary analysis of sequence capture bycatch, with application to an agricultural pest

Sequence capture is a genomic technique that selectively enriches target sequences before high throughput next-generation sequencing, to generate specific sequences of interest. Off-target or 'bycatch' data are often discarded from capture experiments, but can be leveraged to address evolutionary questions under some circumstances. Here, we investigated the effects of missing data on a variety of evolutionary analyses using bycatch from an exon capture experiment on the global pest moth, Helicoverpa armigera. We added > 200 new samples from across Australia in the form of mitogenomes obtained as bycatch from targeted sequence capture, and combined these into an additional larger dataset to total > 1000 mitochondrial cytochrome c oxidase subunit I (COI) sequences across the species' global distribution. Using discriminant analysis of principal components and Bayesian coalescent analyses, we showed that mitogenomes assembled from bycatch with up to 75% missing data were able to return evolutionary inferences consistent with higher coverage datasets and the broader literature surrounding H. armigera. For example, low-coverage sequences broadly supported the delineation of two H. armigera subspecies and also provided new insights into the potential for geographic turnover among these subspecies. However, we also identified key effects of dataset coverage and composition on our results. Thus, low-coverage bycatch data can offer valuable information for population genetic and phylodynamic analyses, but caution is required to ensure the reduced information does not introduce confounding factors, such as sampling biases, that drive inference. We encourage more researchers to consider maximizing the potential of the targeted sequence approach by examining evolutionary questions with their off-target bycatch where possible-especially in cases where no previous mitochondrial data exists-but recommend stratifying data at different genome coverage thresholds to separate sampling effects from genuine genomic signals, and to understand their implications for evolutionary research.

Congruence between ultraconserved element-based matrices and phylotranscriptomic datasets in the scorpion Tree of Life

Scorpions are ancient and historically renowned for their potent venom. Traditionally, the systematics of this group of arthropods was supported by morphological characters, until recent phylogenomic analyses (using RNAseq data) revealed most of the higher-level taxa to be non-monophyletic. While these phylogenomic hypotheses are stable for almost all lineages, some nodes have been hard to resolve due to minimal taxonomic sampling (e.g. family Chactidae). In the same line, it has been shown that some nodes in the Arachnid Tree of Life show disagreement between hypotheses generated using transcritptomes and other genomic sources such as the ultraconserved elements (UCEs). Here, we compared the phylogenetic signal of transcriptomes vs. UCEs by retrieving UCEs from new and previously published scorpion transcriptomes and genomes, and reconstructed phylogenies using both datasets independently. We reexamined the monophyly and phylogenetic placement of Chactidae, sampling an additional chactid species using both datasets. Our results showed that both sets of genome-scale datasets recovered highly similar topologies, with Chactidae rendered paraphyletic owing to the placement of Nullibrotheas allenii. As a first step toward redressing the systematics of Chactidae, we establish the family Anuroctonidae (new family) to accommodate the genus Anuroctonus.

Dynamic evolution of locomotor performance independent of changes in extended phenotype use in spiders

Many animals use self-built structures (extended phenotypes) to enhance body functions, such as thermoregulation, prey capture or defence. Yet, it is unclear whether the evolution of animal constructions supplements or substitutes body functions-with disparate feedbacks on trait evolution. Here, using brown spiders (Araneae: marronoid clade), we explored if the evolutionary loss and gain of silken webs as extended prey capture devices correlates with alterations in traits known to play an important role in predatory strikes-locomotor performance (sprint speed) and leg spination (expression of capture spines on front legs). We found that in this group high locomotor performance, with running speeds of over 100 body lengths per second, evolved repeatedly-both in web-building and cursorial spiders. There was no correlation with running speed, and leg spination only poorly correlated, relative to the use of extended phenotypes, indicating that web use does not reduce selective pressures on body functions involved in prey capture and defence per se. Consequently, extended prey capture devices serve as supplements rather than substitutions to body traits and may only be beneficial in conjunction with certain life-history traits, possibly explaining the rare evolution and repeated loss of trapping strategies in predatory animals.

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.

Shotgun metagenomics captures more microbial diversity than targeted 16S rRNA gene sequencing for field specimens and preserved museum specimens

The use of museum specimens for research in microbial evolutionary ecology remains an under-utilized investigative dimension with important potential. Despite this potential, there remain barriers in methodology and analysis to the wide-spread adoption of museum specimens for such studies. Here, we hypothesized that there would be significant differences in taxonomic prediction and related diversity among sample type (museum or fresh) and sequencing strategy (medium-depth shotgun metagenomic or 16S rRNA gene). We found dramatically higher predicted diversity from shotgun metagenomics when compared to 16S rRNA gene sequencing in museum and fresh samples, with this differential being larger in museum specimens. Broadly confirming these hypotheses, the highest diversity found in fresh samples was with shotgun sequencing using the Rep200 reference inclusive of viruses and microeukaryotes, followed by the WoL reference database. In museum-specimens, community diversity metrics also differed significantly between sequencing strategies, with the alpha-diversity ACE differential being significantly greater than the same comparisons made for fresh specimens. Beta diversity results were more variable, with significance dependent on reference databases used. Taken together, these findings demonstrate important differences in diversity results and prompt important considerations for future experiments and downstream analyses aiming to incorporate microbiome datasets from museum specimens.

Neglected no longer: Phylogenomic resolution of higher-level relationships in Solifugae

Advanced sequencing technologies have expedited resolution of higher-level arthropod relationships. Yet, dark branches persist, principally among groups occurring in cryptic habitats. Among chelicerates, Solifugae ("camel spiders") is the last order lacking a higher-level phylogeny and have thus been historically characterized as "neglected [arachnid] cousins". Though renowned for aggression, remarkable running speed, and xeric adaptation, inferring solifuge relationships has been hindered by inaccessibility of diagnostic morphological characters, whereas molecular investigations have been limited to one of 12 recognized families. Our phylogenomic dataset via capture of ultraconserved elements sampling all extant families recovered a well-resolved phylogeny, with two distinct groups of New World taxa nested within a broader Paleotropical radiation. Divergence times using fossil calibrations inferred that Solifugae radiated by the Permian, and most families diverged prior to the Paleogene-Cretaceous extinction, likely driven by continental breakup. We establish Boreosolifugae new suborder uniting five Laurasian families, and Australosolifugae new suborder uniting seven Gondwanan families using morphological and biogeographic signal.

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.

Too early for the ferry: The biogeographic history of the Assamiidae of southeast Asia (Chelicerata: Opiliones, Laniatores)

Opiliones (harvestmen) have come to be regarded as an abundant source of model groups for study of historical biogeography, due to their ancient age, poor dispersal capability, and high fidelity to biogeographic terranes. One of the least understood harvestman groups is the Paleotropical Assamiidae, one of the more diverse families of Opiliones. Due to a labyrinthine taxonomy, poorly established generic and subfamilial boundaries, and the lack of taxonomic keys for the group, few efforts have been undertaken to decipher relationships within this arachnid lineage. Neither the monophyly of the family, nor its exact placement in the harvestman phylogeny, have been established. Here, we assessed the internal phylogeny of Assamiidae using a ten-locus Sanger dataset, sampling key lineages putatively ascribed to this family for five of the ten markers. Our analyses recovered Assamiidae as a monophyletic group, in a clade with the primarily Afrotropical Pyramidopidae and the southeast Asian Beloniscidae. Internal relationships of assamiids disfavored the systematic validity of subfamilies, with biogeography reflecting much better phylogenetic structure than the existing higher-level taxonomy. To assess whether the Asian assamiids came to occupy Indo-Pacific terranes via rafting on the Indian subcontinent, we performed divergence dating to infer the age of the family. Our results show that Indo-Pacific clades are ancient, originating well before the Cretaceous and therefore predate a vicariant mechanism commonly encountered for Paleotropical taxa.

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.

Life-history traits drive spatial genetic structuring in Dinaric cave spiders

The subterranean ecosystem exerts strong selection pressures on the organisms that thrive in it. In response, obligate cave-dwellers have developed a series of morphological, physiological, and behavioral adaptations, such as eye reduction, appendage elongation, low metabolic rates or intermittent activity patterns, collectively referred to as troglomorphism. Traditionally, studies on cave organisms have been hampered by the difficulty of sampling (i.e., small population sizes, temporal heterogeneity in specimen occurrence, challenges imposed by the difficult-to-access nature of caves). Here, we circumvent this limitation by implementing a museomics approach. Specifically, we aim at comparing the genetic population structures of five cave spider species demonstrating contrasting life histories and levels of troglomorphism across different caves in the northern Dinarides (Balkans, Europe). We applied a genome-wide hybridization-capture approach (i.e., HyRAD) to capture DNA from 117 historical samples. By comparing the population genetic structures among five species and by studying isolation by distance, we identified deeper population structuring and more pronounced patterns of isolation by distance in the highly troglomorphic Parastalita stygia and Stalita pretneri ground dwellers, while the three web-building Troglohyphantes species, two of which can occasionally be found in surface habitats, showed less structured populations compatible with higher dispersal ability. The spatial distribution of genetic groups revealed common phylogeographic breaks among lineages across the studied species, which hint at the importance of environmental features in driving dispersal potential and shaping underground diversity.

Historic and Prehistoric Epidemics: An Overview of Sources Available for the Study of Ancient Pathogens

Since life on earth developed, parasitic microbes have thrived. Increases in host numbers, or the conquest of a new species, provide an opportunity for such a pathogen to enjoy, before host defense systems kick in, a similar upsurge in reproduction. Outbreaks, caused by "endemic" pathogens, and epidemics, caused by "novel" pathogens, have thus been creating chaos and destruction since prehistorical times. To study such (pre)historic epidemics, recent advances in the ancient DNA field, applied to both archeological and historical remains, have helped tremendously to elucidate the evolutionary trajectory of pathogens. These studies have offered new and unexpected insights into the evolution of, for instance, smallpox virus, hepatitis B virus, and the plague-causing bacterium Yersinia pestis. Furthermore, burial patterns and historical publications can help in tracking down ancient pathogens. Another source of information is our genome, where selective sweeps in immune-related genes relate to past pathogen attacks, while multiple viruses have left their genomes behind for us to study. This review will discuss the sources available to investigate (pre)historic diseases, as molecular knowledge of historic and prehistoric pathogens may help us understand the past and the present, and prepare us for future epidemics.

Assessing the performance of historical skins and bones for museomics using wolf specimens as a case study

Advances in the field of museomics have promoted a high sampling demand for natural history collections (NHCs), eventually resulting in damage to invaluable resources to understand historical biodiversity. It is thus essential to achieve a consensus about which historical tissues present the best sources of DNA. In this study, we evaluated the performance of different historical tissues from Iberian wolf NHCs in genome-wide assessments. We targeted three tissues—bone (jaw and femur), maxilloturbinal bone, and skin—that have been favored by traditional taxidermy practices for mammalian carnivores. Specifically, we performed shotgun sequencing and target capture enrichment for 100,000 single nucleotide polymorphisms (SNPs) selected from the commercial Canine HD BeadChip across 103 specimens from 1912 to 2005. The performance of the different tissues was assessed using metrics based on endogenous DNA content, uniquely high-quality mapped reads after capture, and enrichment proportions. All samples succeeded as DNA sources, regardless of their collection year or sample type. Skin samples yielded significantly higher amounts of endogenous DNA compared to both bone types, which yielded equivalent amounts. There was no evidence for a direct effect of tissue type on capture efficiency; however, the number of genotyped SNPs was strictly associated with the starting amount of endogenous DNA. Evaluation of genotyping accuracy for distinct minimum read depths across tissue types showed a consistent overall low genotyping error rate (<7%), even at low (3x) coverage. We recommend the use of skins as reliable and minimally destructive sources of endogenous DNA for whole-genome and target enrichment approaches in mammalian carnivores. In addition, we provide a new 100,000 SNP capture array validated for historical DNA (hDNA) compatible to the Canine HD BeadChip for high-quality DNA. The increasing demand for NHCs as DNA sources should encourage the generation of genomic datasets comparable among studies.

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.

Old specimens for old branches: Assessing effects of sample age in resolving a rapid Neotropical radiation of squirrels

Ultraconserved Elements (UCEs) have been useful to resolve challenging phylogenies of non-model clades, unpuzzling long-conflicted relationships in key branches of the Tree of Life at both deep and shallow levels. UCEs are often reliably recovered from historical samples, unlocking a vast number of preserved natural history specimens for analysis. However, the extent to which sample age and preservation method impact UCE recovery as well as downstream inferences remains unclear. Furthermore, there is an ongoing debate on how to curate, filter, and properly analyze UCE data when locus recovery is uneven across sample age and quality. In the present study we address these questions with an empirical dataset composed of over 3800 UCE loci from 219 historical and modern samples of Sciuridae, a globally distributed and ecologically important family of rodents. We provide a genome-scale phylogeny of two squirrel subfamilies (Sciurillinae and Sciurinae: Sciurini) and investigate their placement within Sciuridae. For historical specimens, recovery of UCE loci and mean length per locus were inversely related to sample age; deeper sequencing improved the number of UCE loci recovered but not locus length. Most of our phylogenetic inferences-performed on six datasets with alternative data-filtering strategies, and using three distinct optimality criteria-resulted in distinct topologies. Datasets containing more loci (40% and 50% taxa representativeness matrices) yielded more concordant topologies and higher support values than strictly filtered datasets (60% matrices) particularly with IQ-Tree and SVDquartets, while filtering based on information content provided better topological resolution for inferences with the coalescent gene-tree based approach in ASTRAL-III. We resolved deep relationships in Sciuridae (including among the five currently recognized subfamilies) and relationships among the deepest branches of Sciurini, but conflicting relationships remain at both genus- and species-levels for the rapid Neotropical tree squirrel radiation. Our results suggest that phylogenomic consensus can be difficult and heavily influenced by the age of available samples and the filtering steps used to optimize dataset properties.

Phylogenomics and Diversification of the Schistosomatidae Based on Targeted Sequence Capture of Ultra-Conserved Elements

Schistosomatidae Stiles and Hassall 1898 is a medically significant family of digenetic trematodes (Trematoda: Digenea), members of which infect mammals or birds as definitive hosts and aquatic or amphibious gastropods as intermediate hosts. Currently, there are 17 named genera, for many of which evolutionary interrelationships remain unresolved. The lack of a resolved phylogeny has encumbered our understanding of schistosomatid evolution, specifically patterns of host-use and the role of host-switching in diversification. Here, we used targeted sequence capture of ultra-conserved elements (UCEs) from representatives of 13 of the 17 named genera and 11 undescribed lineages that are presumed to represent either novel genera or species to generate a phylogenomic dataset for the estimation of schistosomatid interrelationships. This study represents the largest phylogenetic effort within the Schistosomatidae in both the number of loci and breadth of taxon sampling. We present a near-comprehensive family-level phylogeny providing resolution to several clades of long-standing uncertainty within Schistosomatidae, including resolution for the placement of the North American mammalian schistosomes, implying a second separate capture of mammalian hosts. Additionally, we present evidence for the placement of Macrobilharzia at the base of the Schistosoma + Bivitellobilharzia radiation. Patterns of definitive and intermediate host use and a strong role for intermediate host-switching are discussed relative to schistosomatid diversification.

Historical specimens and the limits of subspecies phylogenomics in the New World quails (Odontophoridae)

As phylogenomics focuses on comprehensive taxon sampling at the species and population/subspecies levels, incorporating genomic data from historical specimens has become increasingly common. While historical samples can fill critical gaps in our understanding of the evolutionary history of diverse groups, they also introduce additional sources of phylogenomic uncertainty, making it difficult to discern novel evolutionary relationships from artifacts caused by sample quality issues. These problems highlight the need for improved strategies to disentangle artifactual patterns from true biological signal as historical specimens become more prevalent in phylogenomic datasets. Here, we tested the limits of historical specimen-driven phylogenomics to resolve subspecies-level relationships within a highly polytypic family, the New World quails (Odontophoridae), using thousands of ultraconserved elements (UCEs). We found that relationships at and above the species-level were well-resolved and highly supported across all analyses, with the exception of discordant relationships within the two most polytypic genera which included many historical specimens. We examined the causes of discordance and found that inferring phylogenies from subsets of taxa resolved the disagreements, suggesting that analyzing subclades can help remove artifactual causes of discordance in datasets that include historical samples. At the subspecies-level, we found well-resolved geographic structure within the two most polytypic genera, including the most polytypic species in this family, Northern Bobwhites (Colinus virginianus), demonstrating that variable sites within UCEs are capable of resolving phylogenetic structure below the species level. Our results highlight the importance of complete taxonomic sampling for resolving relationships among polytypic species, often through the inclusion of historical specimens, and we propose an integrative strategy for understanding and addressing the uncertainty that historical samples sometimes introduce to phylogenetic analyses.

Poor hDNA-Derived NGS Data May Provide Sufficient Phylogenetic Information of Potentially Extinct Taxa

Museum material is an important source of metadata for past and recent biological events. With current sequencing technologies, it is possible to obtain historical DNA (hDNA) from older material and/or endangered species to answer taxonomic, systematic, and biogeographical questions. However, hDNA from museum collections is often highly degraded, making it difficult to assess relationships at or above the species level. We therefore studied two probably extinct gastropod species of the genus Laevicaspia, which were collected ∼140 years ago in the Caspian Sea, to map “standard” mitochondrial and nuclear markers and assess both the sequencing depth and the proportion of ambiguous sites as an indicator for the phylogenetic quality of the NGS data. Our study resulted in the first phylogenetically informative mitochondrial and nuclear markers for L. caspia. Assessment of both sequencing depth (mean coverage) and proportion of ambiguous sites suggests that our assembled consensus sequences are reliable for this species. In contrast, no informative gastropod-specific DNA was obtained for L. conus, likely due to a high degree of tissue digestion and contamination with non-gastropod DNA. Nevertheless, our results show that hDNA may in principle yield high-quality sequences for species-level phylogenetic analyses, which underlines the importance of museum collections as valuable archives of the biological past.

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.

Using natural history to guide supervised machine learning for cryptic species delimitation with genetic data

The diversity of biological and ecological characteristics of organisms, and the underlying genetic patterns and processes of speciation, makes the development of universally applicable genetic species delimitation methods challenging. Many approaches, like those incorporating the multispecies coalescent, sometimes delimit populations and overestimate species numbers. This issue is exacerbated in taxa with inherently high population structure due to low dispersal ability, and in cryptic species resulting from nonecological speciation. These taxa present a conundrum when delimiting species: analyses rely heavily, if not entirely, on genetic data which over split species, while other lines of evidence lump. We showcase this conundrum in the harvester Theromaster brunneus, a low dispersal taxon with a wide geographic distribution and high potential for cryptic species. Integrating morphology, mitochondrial, and sub-genomic (double-digest RADSeq and ultraconserved elements) data, we find high discordance across analyses and data types in the number of inferred species, with further evidence that multispecies coalescent approaches over split. We demonstrate the power of a supervised machine learning approach in effectively delimiting cryptic species by creating a "custom" training data set derived from a well-studied lineage with similar biological characteristics as Theromaster. This novel approach uses known taxa with particular biological characteristics to inform unknown taxa with similar characteristics, using modern computational tools ideally suited for species delimitation. The approach also considers the natural history of organisms to make more biologically informed species delimitation decisions, and in principle is broadly applicable for taxa across the tree of life.

Museum Genomics

Natural history collections are invaluable repositories of biological information that provide an unrivaled record of Earth's biodiversity. Museum genomics-genomics research using traditional museum and cryogenic collections and the infrastructure supporting these investigations-has particularly enhanced research in ecology and evolutionary biology, the study of extinct organisms, and the impact of anthropogenic activity on biodiversity. However, leveraging genomics in biological collections has exposed challenges, such as digitizing, integrating, and sharing collections data; updating practices to ensure broadly optimal data extraction from existing and new collections; and modernizing collections practices, infrastructure, and policies to ensure fair, sustainable, and genomically manifold uses of museum collections by increasingly diverse stakeholders. Museum genomics collections are poised to address these challenges and, with increasingly sensitive genomics approaches, will catalyze a future era of reproducibility, innovation, and insight made possible through integrating museum and genome sciences.

Old Brains in Alcohol: The Usability of Legacy Collection Material to Study the Spider Neuroarchitecture

Natural history collections include rare and significant taxa that might otherwise be unavailable for comparative studies. However, curators must balance the needs of current and long-term research. Methods of data extraction that minimize the impact on specimens are therefore favored. Micro-CT has the potential to expose new character systems based on internal anatomy to taxonomic and phylogenetic analysis without dissection or thin sectioning for histology. However, commonly applied micro-CT protocols involve critical point drying, which permanently changes the specimen. Here, we apply a minimally destructive method of specimen preparation for micro-CT investigation of spider neuroanatomy suitable for application to legacy specimens in natural history collections. We used two groups of female spiders of the common species Araneus diadematus—freshly captured (n = 11) vs. legacy material between 70 and 90 years old (n = 10)—to qualitatively and quantitatively assess the viability of micro-CT scanning and the impact of aging on their neuroarchitecture. We statistically compared the volumes of the supraesophageal ganglion (syncerebrum) and used 2D geometric morphometrics to analyze variations in the gross shape of the brain. We found no significant differences in the brain shape or the brain volume relative to the cephalothorax size. Nonetheless, a significant difference was observed in the spider size. We considered such differences to be explained by environmental factors rather than preservation artifacts. Comparison between legacy and freshly collected specimens indicates that museum specimens do not degrade over time in a way that might bias the study results, as long as the basic preservation conditions are consistently maintained, and where lapses in preservation have occurred, these can be identified. This, together with the relatively low-impact nature of the micro-CT protocol applied here, could facilitate the use of old, rare, and valuable material from collections in studies of internal morphology.

Bird nests as botanical time capsules: DNA barcoding identifies the contents of contemporary and historical nests

Bird nests in natural history collections are an abundant yet vastly underutilized source of genetic information. We sequenced the nuclear ribosomal internal transcribed spacer to identify plant species used as nest material in two contemporary (2003 and 2018) and two historical (both 1915) nest specimens constructed by Song Sparrows (Melospiza melodia) and Savannah Sparrows (Passerculus sandwichensis). A total of 13 (22%) samples yielded single, strong bands that could be identified using GenBank resources: six plants (Angiospermae), six green algae (Chlorophyta), and one ciliate (Ciliophora). Two native plant species identified in the nests included Festuca microstachys, which was introduced to the nest collection site by restoration practitioners, and Rosa californica, identified in a nest collected from a lost habitat that existed about 100 years ago. Successful sequencing was correlated with higher sample mass and DNA quality, suggesting future studies should select larger pieces of contiguous material from nests and materials that appear to have been fresh when incorporated into the nest. This molecular approach was used to distinguish plant species that were not visually identifiable, and did not require disassembling the nest specimens as is a traditional practice with nest material studies. The many thousands of nest specimens in natural history collections hold great promise as sources of genetic information to address myriad ecological questions.

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.

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.

Uncovering the genomic and metagenomic research potential in old ethanol-preserved snakes

Natural history museum collections worldwide represent a tremendous resource of information on past and present biodiversity. Fish, reptiles, amphibians and many invertebrate collections have often been preserved in ethanol for decades or centuries and our knowledge on the genomic and metagenomic research potential of such material is limited. Here, we use ancient DNA protocols, combined with shotgun sequencing to test the molecular preservation in liver, skin and bone tissue from five old (1842 to 1964) museum specimens of the common garter snake (Thamnophis sirtalis). When mapping reads to a T. sirtalis reference genome, we find that the DNA molecules are highly damaged with short average sequence lengths (38-64 bp) and high C-T deamination, ranging from 9% to 21% at the first position. Despite this, the samples displayed relatively high endogenous DNA content, ranging from 26% to 56%, revealing that genome-scale analyses are indeed possible from all specimens and tissues included here. Of the three tested types of tissue, bone shows marginally but significantly higher DNA quality in these metrics. Though at least one of the snakes had been exposed to formalin, neither the concentration nor the quality of the obtained DNA was affected. Lastly, we demonstrate that these specimens display a diverse and tissue-specific microbial genetic profile, thus offering authentic metagenomic data despite being submerged in ethanol for many years. Our results emphasize that historical museum collections continue to offer an invaluable source of information in the era of genomics.

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.

Comparison of sequence-capture and ddRAD approaches in resolving species and populations in hexacorallian anthozoans

Genome-level sequencing is the next step in understanding species-level relationships within Anthozoa (soft corals, anemones, stony corals, and their kin) as morphological and PCR-directed (single-locus) sequencing methods often fall short of differentiating species. The sea anemone genus Metridium is a common northern temperate sea anemone whose species are difficult to differentiate using morphology alone. Here we use Metridium as a case study to confirm the low level of information available in six loci for species differentiation commonly sequenced for Actiniaria and explore and compare the efficacy of ddRAD and sequence-capture methods in species-level systematics and biogeographic studies. We produce phylogenetic trees from concatenated datasets and perform DAPC and STRUCTURE analyses using SNP data. The six conventional loci are not able to consistently differentiate species within Metridium. The sequence-capture dataset resulted in high support and resolution for both current species and relationships between geographic areas. The ddRAD datasets displayed ambiguity among species, and support between major geographic groupings was not as high as the sequence-capture datasets. The level of resolution and support resulting from the sequence-capture data, combined with the ability to add additional individuals and expand beyond the genus Metridium over time, emphasizes the utility of sequence-capture methods for both systematics and future biogeographic studies within anthozoans. We discuss the strengths and weaknesses of the genomic approaches in light of our findings and suggest potential implications for the biogeography of Metridium based on our sampling.

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.

Partitioned Gene-Tree Analyses and Gene-Based Topology Testing Help Resolve Incongruence in a Phylogenomic Study of Host-Specialist Bees (Apidae: Eucerinae)

Incongruence among phylogenetic results has become a common occurrence in analyses of genome-scale data sets. Incongruence originates from uncertainty in underlying evolutionary processes (e.g., incomplete lineage sorting) and from difficulties in determining the best analytical approaches for each situation. To overcome these difficulties, more studies are needed that identify incongruences and demonstrate practical ways to confidently resolve them. Here, we present results of a phylogenomic study based on the analysis 197 taxa and 2,526 ultraconserved element (UCE) loci. We investigate evolutionary relationships of Eucerinae, a diverse subfamily of apid bees (relatives of honey bees and bumble bees) with >1,200 species. We sampled representatives of all tribes within the group and >80% of genera, including two mysterious South American genera, Chilimalopsis and Teratognatha. Initial analysis of the UCE data revealed two conflicting hypotheses for relationships among tribes. To resolve the incongruence, we tested concatenation and species tree approaches and used a variety of additional strategies including locus filtering, partitioned gene-trees searches, and gene-based topological tests. We show that within-locus partitioning improves gene tree and subsequent species-tree estimation, and that this approach, confidently resolves the incongruence observed in our data set. After exploring our proposed analytical strategy on eucerine bees, we validated its efficacy to resolve hard phylogenetic problems by implementing it on a published UCE data set of Adephaga (Insecta: Coleoptera). Our results provide a robust phylogenetic hypothesis for Eucerinae and demonstrate a practical strategy for resolving incongruence in other phylogenomic data sets.

Enhancing CT imaging: A safe protocol to stain and de-stain rare fetal museum specimens using diffusible iodine-based staining (diceCT)

Computed tomography (CT) scanning is being increasingly employed in the study of natural history, particularly to investigate the internal anatomy of unique specimens in museum collections. Different techniques to enhance the contrast between tissues have been developed to improve the quality of the scans while preserving the integrity of these rare specimens. Diffusible iodine-based contrast enhanced computed tomography (diceCT) was found to be particularly effective and reversible for staining tissues in formalin preserved specimens. While it can also be effectively employed to stain ethanol-preserved specimens of small size, the reversibility of this process and the applicability to large-bodied animals has never been thoroughly tested. Here, we describe a novel diceCT protocol developed to stain and de-stain ethanol-preserved prenatal specimens of baleen whales (Mysticeti, Cetacea). These large (10-90 cm in length only considering early fetal stages) specimens present unique challenges as they are rare in collections and irreplaceable, therefore it is imperative to not damage them with the staining process. Before trying this protocol on baleen whales' specimens, we conducted a pilot study on commercially available fetal pigs using the same parameters. This allowed us to optimize the staining time to obtain the best results in CT scanning and to test first-hand the effect of de-staining on ethanol-based specimens. External coloration of the specimens is also a concern with iodine-staining, as stained specimens assume a bright red color that needs to be removed from both internal and external tissues before they can be stored. To test the reversibility of the stain in ethanol-preserved specimens with fur, we also conducted a small experiment using commercially acquired domestic mice. After these trials were successful, we applied the staining and de-staining protocol to multiple fetal specimens of mysticetes up to 90 cm in length, both ethanol- and formalin-preserved. Specimens were soaked in a solution of 1% pure iodine in 70% ethanol for 1-28 days, according to their size. After scanning, specimens are soaked in a solution of 3% sodium thiosulfate in 70% ethanol that is able to completely wash out the iodine from the tissues in a shorter time frame, between a few hours and 14 days. The same concentrations were used for formalin-preserved specimens, but DI water was used as solvent instead of ethanol. The staining technique proved particularly useful to enhance the contrast difference between cartilage, mineralized bone, teeth, and the surrounding soft tissues even when the specimens where scanned in medical-grade CT scans. The specimens did not suffer any visible damage or shrinkage due to the procedure, and in both the fetal samples and in the mice the color of the stain was completely removed by the de-staining process. We conclude therefore that this protocol can be safely applied to a variety of ethanol-preserved museum specimens to enhance the quality of the CT scanning and highlight internal morphological features without recurring to dissection or other irreversible procedures. We also provide tips to best apply this protocol, from how to mix the solutions to how to minimize the staining time.

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 re-evaluation of Triaenonychoidea (Opiliones : Laniatores), and systematics of Triaenonychidae, including new families, genera and species

The Opiliones superfamily Triaenonychoidea currently includes two families, the monogeneric New Zealand–endemic Synthetonychiidae Forster, 1954 and Triaenonychidae Sørensen, 1886, a diverse family distributed mostly throughout the temperate Gondwanan terranes, with ~110 genera and ~500 species and subspecies currently described. Traditionally, Triaenonychidae has been divided into subfamilies diagnosed by very few morphological characters largely derived from the troublesome ‘Roewerian system’ of morphology, and classifications based on this system led to many complications. Recent research within Triaenonychoidea using morphology and traditional multilocus data has shown multiple deeply divergent lineages, non-monophyly of Triaenonychidae, and non-monophyly of subfamilies, necessitating a revision based on phylogenomic data. We used sequence capture of ultraconserved elements across 164 samples to create a 50% taxon occupancy matrix with 704 loci. Using phylogenomic and morphological examinations, we explored family-level relationships within Triaenonychoidea, including describing two new families: (1) Lomanellidae Mendes &amp; Derkarabetian, fam. nov., consisting of Lomanella Pocock, 1903, and a newly described genus Abaddon Derkarabetian &amp; Baker, gen. nov. with one species, A. despoliator Derkarabetian, sp. nov.; and (2) the elevation to family of Buemarinoidae Karaman, 2019, consisting of Buemarinoa Roewer, 1956, Fumontana Shear, 1977, Flavonuncia Lawrence, 1959, and a newly described genus Turonychus Derkarabetian, Prieto &amp; Giribet, gen. nov., with one species, T. fadriquei Derkarabetian, Prieto &amp; Giribet, sp. nov. With our dataset we also explored phylogenomic relationships within Triaenonychidae with an extensive taxon set including samples representing ~80% of the genus-level diversity. Based on our results we (1) discuss systematics of this family including the historical use of subfamilies, (2) reassess morphology in the context of our phylogeny, (3) hypothesise placement for all unsampled genera, (4) highlight lineages most in need of taxonomic revision, and (5) provide an updated species-level checklist. Aside from describing new taxa, our study provides the phylogenomic context necessary for future evolutionary and systematic research across this diverse lineage. ZooBank Registration: urn:lsid:zoobank.org:pub:81683834-98AB-43AA-B25A-C28C6A404F41

A molecular phylogeny of historical and contemporary specimens of an under-studied micro-invertebrate group

Resolution of relationships at lower taxonomic levels is crucial for answering many evolutionary questions, and as such, sufficiently varied species representation is vital. This latter goal is not always achievable with relatively fresh samples. To alleviate the difficulties in procuring rarer taxa, we have seen increasing utilization of historical specimens in building molecular phylogenies using high throughput sequencing. This effort, however, has mainly focused on large-bodied or well-studied groups, with small-bodied and under-studied taxa under-prioritized. Here, we utilize both historical and contemporary specimens, to increase the resolution of phylogenetic relationships among a group of under-studied and small-bodied metazoans, namely, cheilostome bryozoans. In this study, we pioneer the sequencing of air-dried cheilostomes, utilizing a recently developed library preparation method for low DNA input. We evaluate a de novo mitogenome assembly and two iterative methods, using the sequenced target specimen as a reference for mapping, for our sequences. In doing so, we present mitochondrial and ribosomal RNA sequences of 43 cheilostomes representing 37 species, including 14 from historical samples ranging from 50 to 149 years old. The inferred phylogenetic relationships of these samples, analyzed together with publicly available sequence data, are shown in a statistically well-supported 65 taxa and 17 genes cheilostome tree, which is also the most broadly sampled and largest to date. The robust phylogenetic placement of historical samples whose contemporary conspecifics and/or congenerics have been sequenced verifies the appropriateness of our workflow and gives confidence in the phylogenetic placement of those historical samples for which there are no close relatives sequenced. The success of our workflow is highlighted by the circularization of a total of 27 mitogenomes, seven from historical cheilostome samples. Our study highlights the potential of utilizing DNA from micro-invertebrate specimens stored in natural history collections for resolving phylogenetic relationships among species.

A molecular phylogeny of the circum-Antarctic Opiliones family Neopilionidae

The Opiliones family Neopilionidae is restricted to the terranes of the former temperate Gondwana: South America, Africa, Australia, New Caledonia and New Zealand. Despite decades of morphological study of this unique fauna, it has been difficult reconciling the classic species of the group (some described over a century ago) with recent cladistic morphological work and previous molecular work. Here we attempted to investigate the pattern and timing of diversification of Neopilionidae by sampling across the distribution range of the family and sequencing three markers commonly used in Sanger-based approaches (18S rRNA, 28S rRNA and cytochrome-c oxidase subunit I). We recovered a well-supported and stable clade including Ballarra (an Australian ballarrine) and the Enantiobuninae from South America, Australia, New Caledonia and New Zealand, but excluding Vibone (a ballarrine from South Africa). We further found a division between West and East Gondwana, with the South American Thrasychirus/Thrasychiroides always being sister group to an Australian–Zealandian (i.e. Australia + New Zealand + New Caledonia) clade. Resolution of the Australian–Zealandian taxa was analysis-dependent, but some analyses found Martensopsalis, from New Caledonia, as the sister group to an Australian–New Zealand clade. Likewise, the species from New Zealand formed a clade in some analyses, but Mangatangi often came out as a separate lineage from the remaining species. However, the Australian taxa never constituted a monophyletic group, with Ballarra always segregating from the remaining Australian species, which in turn constituted 1–3 clades, depending on the analysis. Our results identify several generic inconsistencies, including the possibility of Thrasychiroides nested within Thrasychirus, Forsteropsalis being paraphyletic with respect to Pantopsalis, and multiple lineages of Megalopsalis in Australia. In addition, the New Zealand Megalopsalis need generic reassignment: Megalopsalis triascuta will require its own genus and M. turneri is here transferred to Forsteropsalis, as Forsteropsalis turneri (Marples, 1944), comb. nov.

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.

Natural history collections are critical resources for contemporary and future studies of urban evolution

Urban environments are among the fastest changing habitats on the planet, and this change has evolutionary implications for the organisms inhabiting them. Herein, we demonstrate that natural history collections are critical resources for urban evolution studies. The specimens housed in these collections provide great potential for diverse types of urban evolution research, and strategic deposition of specimens and other materials from contemporary studies will determine the resources and research questions available to future urban evolutionary biologists. As natural history collections are windows into the past, they provide a crucial historical timescale for urban evolution research. While the importance of museum collections for research is generally appreciated, their utility in the study of urban evolution has not been explicitly evaluated. Here, we: (a) demonstrate that museum collections can greatly enhance urban evolution studies, (b) review patterns of specimen use and deposition in the urban evolution literature, (c) analyze how urban versus rural and native versus nonnative vertebrate species are being deposited in museum collections, and (d) make recommendations to researchers, museum professionals, scientific journal editors, funding agencies, permitting agencies, and professional societies to improve archiving policies. Our analyses of recent urban evolution studies reveal that museum specimens can be used for diverse research questions, but they are used infrequently. Further, although nearly all studies we analyzed generated resources that could be deposited in natural history collections (e.g., collected specimens), a minority (12%) of studies actually did so. Depositing such resources in collections is crucial to allow the scientific community to verify, replicate, and/or re-visit prior research. Therefore, to ensure that adequate museum resources are available for future urban evolutionary biology research, the research community-from practicing biologists to funding agencies and professional societies-must make adjustments that prioritize the collection and deposition of urban specimens.