Ultraconserved element (UCE) probe set design: Base genome and initial design parameters critical for optimization

Targeted capture and enrichment approaches have proven effective for phylogenetic study. Ultraconserved elements (UCEs) in particular have exhibited great utility for phylogenomic analyses, with the software package phyluce being among the most utilized pipelines for UCE phylogenomics, including probe design. Despite the success of UCEs, it is becoming increasing apparent that diverse lineages require probe sets tailored to focal taxa in order to improve locus recovery. However, factors affecting probe design and methods for optimizing probe sets to focal taxa remain underexplored. Here, we use newly available beetle (Coleoptera) genomic resources to investigate factors affecting UCE probe set design using phyluce . In particular, we explore the effects of stringency during initial design steps, as well as base genome choice on resulting probe sets and locus recovery. We found that both base genome choice and initial bait design stringency parameters greatly alter the number of resultant probes included in final probe sets and strongly affect the number of loci detected and recovered during in silico testing of these probe sets. In addition, we identify attributes of base genomes that correlated with high performance in probe design. Ultimately, we provide a recommended workflow for using Phyluce to design an optimized UCE probe set that will work across a targeted lineage, and use our findings to develop a new, open‐source UCE probe set for beetles of the suborder Adephaga.     

Advancing mite phylogenomics: Designing ultraconserved elements for Acari phylogeny

Mites (Acari) are one of the most diverse groups of life on Earth; yet, their evolutionary relationships are poorly understood. Also, the resolution of broader arachnid phylogeny has been hindered by an underrepresentation of mite diversity in phylogenomic analyses. To further our understanding of Acari evolution, we design targeted ultraconserved genomic elements (UCEs) probes, intended for resolving the complex relationships between mite lineages and closely related arachnids. We then test our Acari UCE baits in‐silico by constructing a phylogeny using 13 existing Acari genomes, as well as 6 additional taxa from a variety of genomic sources. Our Acari‐specific probe kit improves the recovery of loci within mites over an existing general arachnid UCE probe set. Our initial phylogeny recovers the major mite lineages, yet finds mites to be non‐monophyletic overall, with Opiliones (harvestmen) and Ricinuleidae (hooded tickspiders) rendering Parasitiformes paraphyletic.     

Streamlining universal single‐copy orthologue and ultraconserved element design: A case study in Collembola

Genomic data sets are increasingly central to ecological and evolutionary biology, but far fewer resources are available for invertebrates. Powerful new computational tools and the rapidly decreasing cost of Illumina sequencing are beginning to change this, enabling rapid genome assembly and reference marker extraction. We have developed and tested a practical workflow for developing genomic resources in nonmodel groups with real‐world data on Collembola (springtails), one of the most dominant soil animals on Earth. We designed universal molecular marker sets, single‐copy orthologues (BUSCO s) and ultraconserved elements (UCEs), using three existing and 11 newly generated genomes. Both marker types were tested in silico via marker capture success and phylogenetic performance. The new genomes were assembled with Illumina short reads and 9,585?14,743 protein‐coding genes were predicted with ab initio and protein homology evidence. We identified 1,997 benchmarking universal single‐copy orthologues (BUSCO s) across 14 genomes and created and assessed a custom BUSCO data set for extracting single‐copy genes. We also developed a new UCE probe set containing 46,087 baits targeting 1,885 loci. We successfully captured 1,437?1,865 BUSCO s and 975?1,186 UCEs across 14 genomes. Phylogenomic reconstructions using these markers proved robust, giving new insight on deep‐time collembolan relationships. Our study demonstrates the feasibility of generating thousands of universal markers from highly efficient whole‐genome sequencing, providing a valuable resource for genome‐scale investigations in evolutionary biology and ecology.     

Ultraconserved elements show utility in phylogenetic inference of Adephaga (Coleoptera) and suggest paraphyly of ‘Hydradephaga’

The beetle suborder Adephaga has been the subject of many phylogenetic reconstructions utilizing a variety of data sources and inference methods. However, no strong consensus has yet emerged on the relationships among major adephagan lineages. Ultraconserved elements (UCEs) have proved useful for inferring difficult or unresolved phylogenies at varying timescales in vertebrates, arachnids and Hymenoptera. Recently, a UCE bait set was developed for Coleoptera using polyphagan genomes and a member of the order Strepsiptera as an outgroup. Here, we examine the utility of UCEs for reconstructing the phylogeny of adephagan families, in the first in vitro application a UCE bait set in Coleoptera. Our final dataset included 305 UCE loci for 18 representatives of all adephagan families except Aspidytidae, and two polyphagan outgroups, with a total concatenated length of 83 547 bp. We inferred trees using maximum likelihood analyses of the concatenated UCE alignment and coalescent species tree methods (astral ii , ASTRID, svdquartets ). Although the coalescent species tree methods had poor resolution and weak support, concatenated analyses produced well‐resolved, highly supported trees. Hydradephaga was recovered as paraphyletic, with Gyrinidae sister to Geadephaga and all other adephagans. Haliplidae was recovered as sister to Dytiscoidea, with Hygrobiidae and Amphizoidae successive sisters to Dytiscidae. Finally, Noteridae was recovered as monophyletic and sister to Meruidae. Given the success of UCE data for resolving phylogenetic relationships within Adephaga, we suggest the potential for further resolution of relationships within Adephaga using UCEs with improved taxon sampling, and by developing Adephaga‐specific probes.     

Phylogenomic species delimitation and host‐symbiont coevolution in the fungus‐farming ant genus Sericomyrmex Mayr (Hymenoptera: Formicidae): ultraconserved elements (UCEs) resolve a recent radiation

Ants in the Neotropical genus Sericomyrmex Mayr cultivate fungi for food. Both ants and fungi are obligate, coevolved symbionts. The taxonomy of Sericomyrmex is problematic because the morphology of the worker caste is generally homogeneous across all of the species within the genus, species limits are vague, and the relationships between them are unknown. We used ultraconserved elements (UCEs) as genome‐scale markers to reconstruct evolutionary history and to infer species boundaries in Sericomyrmex. We recovered an average of ～990 UCE loci for 88 Sericomyrmex samples from across the geographical range of the genus as well as for five outgroup taxa. Using maximum likelihood and species‐tree approaches, we recovered nearly identical topologies across datasets with 50–95% matrix completeness. We identify nine species‐level lineages in Sericomyrmex, including two new species. This is less than the previously described 19 species, even accounting for two species for which we had no UCE samples, which brings the total number of Sericomyrmex species to 11. Divergence‐dating analyses recovered 4.3 Ma as the crown‐group age estimates for Sericomyrmex, indicating a recent, rapid radiation. We also sequenced mitochondrial cytochrome oxidase subunit I (COI) for 125 specimens. Resolution and support for clades in our COI phylogeny are weak, indicating that COI is not an appropriate species‐delimitation tool. However, taxa within species consistently cluster together, suggesting that COI is useful as a species identification (‘DNA barcoding’) tool. We also sequenced internal transcribed spacer (ITS) and large subunit (LSU) for 32 Sericomyrmex fungal cultivars. The fungal phylogeny confirms that Sericomyrmex fungi are generalized higher‐attine cultivars, interspersed with Trachymyrmex‐associated fungal species, indicating cultivar sharing and horizontal transfer between these two genera. Our results indicate that UCEs offer immense potential for delimiting and resolving relationships of problematic, recently diverged species.     

Squamate Conserved Loci (SqCL): A unified set of conserved loci for phylogenomics and population genetics of squamate reptiles

The identification of conserved loci across genomes, along with advances in target capture methods and high‐throughput sequencing, has helped spur a phylogenomics revolution by enabling researchers to gather large numbers of homologous loci across clades of interest with minimal upfront investment in locus design. Target capture for vertebrate animals is currently dominated by two approaches—anchored hybrid enrichment (AHE) and ultraconserved elements (UCE)—and both approaches have proven useful for addressing questions in phylogenomics, phylogeography and population genomics. However, these two sets of loci have minimal overlap with each other; moreover, they do not include many traditional loci that that have been used for phylogenetics. Here, we combine across UCE, AHE and traditional phylogenetic gene locus sets to generate the Squamate Conserved Loci set, a single integrated probe set that can generate high‐quality and highly complete data across all three loci types. We use these probes to generate data for 44 phylogenetically disparate taxa that collectively span approximately 33% of terrestrial vertebrate diversity. Our results generated an average of 4.29 Mb across 4709 loci per individual, of which an average of 2.99 Mb was sequenced to high enough coverage (≥10×) to use for population genetic analyses. We validate the utility of these loci for both phylogenomic and population genomic questions, provide a comparison among these locus sets of their relative usefulness and suggest areas for future improvement.     

New approaches to species delimitation and population structure of anthozoans: Two case studies of octocorals using ultraconserved elements and exons

As coral populations decline worldwide in the face of ongoing environmental change, documenting their distribution, diversity and conservation status is now more imperative than ever. Accurate delimitation and identification of species is a critical first step. This task, however, is not trivial as morphological variation and slowly evolving molecular markers confound species identification. New approaches to species delimitation in corals are needed to overcome these challenges. Here, we test whether target enrichment of ultraconserved elements (UCEs) and exons can be used for delimiting species boundaries and population structure within species of corals by focusing on two octocoral genera, Alcyonium and Sinularia, as exemplary case studies. We designed an updated bait set (29,181 baits) to target‐capture 3,023 UCE and exon loci, recovering a mean of 1,910 ± 168 SD per sample with a mean length of 1,055 ± 208 bp. Similar numbers of loci were recovered from Sinularia (1,946 ± 227 SD) and Alcyonium (1,863 ± 177 SD). Species‐level phylogenies were highly supported for both genera. Clustering methods based on filtered single nucleotide polymorphisms delimited species and populations that are congruent with previous allozyme, DNA barcoding, reproductive and ecological data for Alcyonium, and offered further evidence of hybridization among species. For Sinularia, results were congruent with those obtained from a previous study using restriction site associated DNA sequencing. Both case studies demonstrate the utility of target‐enrichment of UCEs and exons to address a wide range of evolutionary and taxonomic questions across deep to shallow timescales in corals.     

Sequence capture phylogenomics of historical ethanol‐preserved museum specimens: Unlocking the rest of the vault

Natural history collections play a crucial role in biodiversity research, and museum specimens are increasingly being incorporated into modern genetics‐based studies. Sequence capture methods have proven incredibly useful for phylogenomics, providing the additional ability to sequence historical museum specimens with highly degraded DNA, which until recently have been deemed less valuable for genetic work. The successful sequencing of ultraconserved elements (UCEs) from historical museum specimens has been demonstrated on multiple tissue types including dried bird skins, formalin‐fixed squamates and pinned insects. However, no study has thoroughly demonstrated this approach for historical ethanol‐preserved museum specimens. Alongside sequencing of “fresh” specimens preserved in >95% ethanol and stored at ?80°C, we used extraction techniques specifically designed for degraded DNA coupled with sequence capture protocols to sequence UCEs from historical museum specimens preserved in 70%–80% ethanol and stored at room temperature, the standard for such ethanol‐preserved museum collections. Across 35 fresh and 15 historical museum samples of the arachnid order Opiliones, an average of 345 UCE loci were included in phylogenomic matrices, with museum samples ranging from six to 495 loci. We successfully demonstrate the inclusion of historical ethanol‐preserved museum specimens in modern sequence capture phylogenomic studies, show a high frequency of variant bases at the species and population levels, and from off‐target reads successfully recover multiple loci traditionally sequenced in multilocus studies including mitochondrial loci and nuclear rRNA loci. The methods detailed in this study will allow researchers to potentially acquire genetic data from millions of ethanol‐preserved museum specimens held in collections worldwide.     

Phylogenomic analysis of the beetle suborder Adephaga with comparison of tailored and generalized ultraconserved element probe performance

Adephaga is the second largest suborder of beetles (Coleoptera) and they serve as important arthropod predators in both aquatic and terrestrial ecosystems. The suborder is divided into Geadephaga comprising terrestrial families and Hydradephaga for aquatic lineages. Despite numerous studies, phylogenetic relationships among the adephagan families and monophyly of the Hydradephaga itself remain in question. Here we conduct a comprehensive phylogenomic analysis of the suborder using ultraconserved elements (UCEs). This study presents the first in vitro test of a newly developed UCE probe set customized for use within Adephaga that includes both probes tailored specifically for the suborder, alongside generalized Coleoptera probes previously found to work in adephagan taxa. We assess the utility of the entire probe set, as well as comparing the tailored and generalized probes alone for reconstructing evolutionary relationships. Our analyses recovered strong support for the paraphyly of Hydradephaga with whirligig beetles (Gyrinidae) placed as sister to all other adephagan families. Geadephaga was strongly supported as monophyletic and placed sister to a clade composed of Haliplidae + Dytiscoidea. Monophyly of Dytiscoidea was strongly supported with relationships among the dytiscoid families resolved and strongly supported. Relationships among the subfamilies of Dytiscidae were strongly supported but largely incongruent with prior phylogenetic estimates for the family. The results of our UCE probe comparison showed that tailored probes alone outperformed generalized probes alone, as well as the full combined probe set (containing both types of probes), under decreased taxon sampling. When taxon sampling was increased, the full combined probe set outperformed both tailored probes and generalized probes alone. This study provides further evidence that UCE probe sets customized for a focal group result in a greater number of recovered loci and substantially improve phylogenomic analysis.     

Hidden histories of gene flow in highland birds revealed with genomic markers

Genomic studies are revealing that divergence and speciation are marked by gene flow, but it is not clear whether gene flow has played a prominent role during the generation of biodiversity in species‐rich regions of the world where vicariance is assumed to be the principal mode by which new species form. We revisit a well‐studied organismal system in the Mexican Highlands, Aphelocoma jays, to test for gene flow among Mexican sierras. Prior results from mitochondrial DNA (mtDNA) largely conformed to the standard model of allopatric divergence, although there was also evidence for more obscure histories of gene flow in a small sample of nuclear markers. We tested for these ‘hidden histories’ using genomic markers known as ultraconserved elements (UCEs) in concert with phylogenies, clustering algorithms and newer introgression tests specifically designed to detect ancient gene flow (e.g. ABBA/BABA tests). Results based on 4303 UCE loci and 2500 informative SNPs are consistent with varying degrees of gene flow among highland areas. In some cases, gene flow has been extensive and recent (although perhaps not ongoing today), whereas in other cases there is only a trace signature of ancient gene flow among species that diverged as long as 5 million years ago. These results show how a species complex thought to be a model for vicariance can reveal a more reticulate history when a broader portion of the genome is queried. As more organisms are studied with genomic data, we predict that speciation‐with‐bouts‐of‐gene‐flow will turn out to be a common mode of speciation.     

Detection of genome specific monomorphic loci in Bos taurus and Bubalus bubalis with oligodeoxyribonucleotide probe

A synthetic oligodeoxyribonucleotide probe (OAT36) comprising nine repeats of 5'GACA 3′ and several enzymes were used to analyse cow, (Bos taurus) and buffalo (Bubalus bubalis) genomes and a number of monomorphic loci were detected in both the species. Different animals from the same species showed an almost ‘similar’ monomorphic hybridization pattern but animals from two separate species showed a different ‘genome specific’ pattern. The overall hybridization with any enzyme and probe combination was found to be unique to one species. This forms the basis of genome specific hybridization which is substantiated by our zoo-blot hybridization studies. The evolutionary aspect of these loci in the context of sequence polymorphisms is discussed.     

Monophyletic groups within 'higher land birds'– comparison of morphological and molecular data

The relationships within the ‘higher land birds’ and putatively related taxa are analysed in a study using 89 morphological characters and DNA sequences of three nuclear, protein‐coding genes, c‐myc, RAG‐1, and myoglobin intron II. Separate analyses of the different data sets and a ‘total evidence’ analysis in which the data sets of the morphological and molecular analyses were combined are compared. All three analyses support the hitherto disputed sister group relationship between Pici (Ramphastidae, Indicatoridae and Picidae) and Galbulae (Galbulidae and Bucconidae). Previously unrecognized osteological synapomorphies of this clade are presented. All analyses further resulted in monophyly of the taxon [Aegothelidae + (Apodidae/Hemiprocnidae + Trochilidae)]. Analysis of the morphological data and of the combined data set also supported monophyly of the taxon [Strigiformes + (Falconidae + Accipitridae)]. The morphological data further support monophyly of the taxon (Upupidae + Bucerotidae). Other placements in the three analyses received either no or only weak bootstrap support.