Identifying conserved genomic elements and designing universal bait sets to enrich them

Biogeography confounds the signal of cospeciation in Batesian mimicry

Since the inception of the field of evolution, mimicry has yielded insights into foundational evolutionary processes, including adaptive peak shifts, speciation, and the emergence and maintenance of phenotypic polymorphisms.1,2,3 In recent years, the coevolutionary processes generating mimicry have gained increasing attention from researchers. Despite significant advances in understanding Batesian and Müllerian mimicry in Lepidopteran systems, few other mimetic systems have received similar detailed research. Here, we present a Batesian mimicry complex involving flightless, armored Pachyrhynchus weevils and their winged Doliops longhorn beetle mimics and examine their coevolutionary patterns within the Philippine archipelagos. Pachyrhynchus weevils are primarily found in the Philippines, where distinct species radiations have occurred on different islands, each with unique color patterns serving as a warning to predators. This defensive trait and mimicry between unrelated species were first described by Wallace in 1889. Notably, the distantly related longhorn beetle Doliops, despite being soft-bodied and ostensibly palatable, mimics the heavily armored, flightless Pachyrhynchus. To address mimicry in this system, we reconstructed the phylogeny of Doliops using a probe set consisting of 38,000 ultraconserved elements. Our study examines the following questions central to understanding the Pachyrhynchus-Doliops mimicry system: (1) to what extent are coevolutionary interactions conserved (i.e., lineage-constrained) and (2) are the codiversification patterns primarily driven by biotic or abiotic factors?4,5,6 To assess color mimicry and cospeciation, we examined the evolution of nanostructure-based warning colors and the effect of island biogeography on cospeciation. Our findings demonstrate the beetle's ability to repeatedly evolve multiple solutions to similar evolutionary challenges, evolving similar color patterns using different types of photonic crystals with varying degrees of order. We revealed that the observed pattern of cospeciation is driven mainly by abiotic factors from their biogeographic history. Unlike the patterns of coevolution seen between angiosperms and insect lineages,7 most ecological interactions do not persist longer than a few million years, leading to patterns of modularity rather than ecological nestedness.4,6,7.

Tenkana, a new genus of jumping spiders (Salticidae, Plexippina) from South Asia, with the new Indian species Tenkanajayamangali

We describe a new plexippine genus, Tenkana gen. nov., supported by phylogenomic data from ultraconserved elements (UCEs), Sanger sequences of four genes, and morphological evidence. The type species, Tenkanamanu (Caleb, Christudhas, Laltanpuii & Chitra, 2014), comb. nov. is transferred from Colopsus, as is Tenkanaarkavathi (Caleb, 2022), comb. nov. The phylogenomic data places Tenkana among the plexippines near Hyllus C.L. Koch, 1846 and Telamonia Thorell, 1887, while the constrained four-gene phylogeny indicates that Tenkana is distinct from Colopsus. Additionally, we describe a new species, Tenkanajayamangali sp. nov.

Ultraconserved elements from transcriptome and genome data provide insight into the phylogenomics of Sternorrhyncha (Insecta: Hemiptera)

Sternorrhyncha, one of the four major suborders of Hemiptera, is a phytophagous taxon inclusive of nearly 18 000 described species. The phylogenetic relationships within the taxon and the earliest-branching lineage of its infraorders remain incompletely understood. This study attempted to illuminate the phylogenetic relationships within Sternorrhyncha through the use of maximum likelihood, Bayesian inference and maximum parsimony analyses, employing ultraconserved element (UCE) data from 39 genomic and 62 transcriptomic datasets and thereby representing most families within the taxon. The probe set Hemiptera 2.7Kv1 was used to recover a total of 2731 UCE loci: from 547 to 1699 (with an average of 1084) across all genomic datasets and from 108 to 849 (with an average of 329) across all transcriptomic datasets. All three types of phylogenetic analyses employed in this study produced robust statistical support for Sternorrhyncha being a monophyletic group. The different methods of phylogenetic analysis produced inconsistent descriptions of topological structure at the infraorder level: while maximum likelihood and Bayesian inference analyses produced strong statistical evidence (100%) indicating the clade Psylloidea + Aleyrodoidea to be a sister of the clade Aphidoidea (Aphidomorpha) + Coccoidea (Coccomorpha), the maximum parsimony analysis failed to recover a similar result. Our results also provide detail on the phylogenetic relationships within each infraorder. This study presents the first use of UCE data to investigate the phylogeny of Sternorrhyncha. It also shows the viability of amalgamating genomic and transcriptomic data in studies of phylogenetic relationships, potentially highlighting a resource-efficient approach for future inquiries into diverse taxa through the integration of varied data sources.

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.

Phylogenomics resolves a 100-year-old debate regarding the evolutionary history of caddisflies (Insecta: Trichoptera)

Trichoptera (caddisfly) phylogeny provides an interesting example of aquatic insect evolution, with rich ecological diversification, especially for underwater architecture. Trichoptera provide numerous critical ecosystem services and are also one of the most important groups of aquatic insects for assessing water quality. The phylogenetic relationships of Trichoptera have been debated for nearly a century. In particular, the phylogenetic position of the "cocoon-makers" within Trichoptera has long been contested. Here, we designed a universal single-copy orthologue and sets of ultraconserved element markers specific for Trichoptera for the first time. Simultaneously, we reconstructed the phylogenetic relationship of Trichoptera based on genome data from 111 species, representing 29 families and 71 genera. Our phylogenetic inference clarifies the probable phylogenetic relationships of "cocoon-makers" within Integripalpia. Hydroptilidae is considered as the basal lineage within Integripalpia, and the families Glossosomatidae, Hydrobiosidae, and Rhyacophilidae formed a monophyletic clade, sister to the integripalpian subterorder Phryganides. The resulting divergence time and ancestral state reconstruction suggest that the most recent common ancestor of Trichoptera appeared in the early Permian and that diversification was strongly correlated with habitat change.

Functional trait mismatch between native and introduced bee pollinators servicing a global fruit crop

BACKGROUND

Understanding connections between biodiversity and ecosystem services can be enhanced by shifting focus from species richness to functional trait-based approaches, that when paired with comparative phylogenetic methods can provide even deeper insights. We investigated the functional ecology and phylogenetic diversity of pollination services provided by hymenopteran insects visiting apple flowers in orchards surrounded by either 'natural' or 'disturbed' landscapes in New South Wales, Australia. We assessed whether morphological and behavioural traits (hairiness, body size, glossa length, pollen load purity, and probability of loose pollen) exhibited non-random phylogenetic patterns. Then, explored whether bees, the primary pollinators in this system, filled unique or overlapping functional entities (FEs). For each landscape, we calculated phylogenetic diversity and used FEs to assess functional richness, evenness, and diversion.

RESULTS

A phylogenomic matrix based on ultraconserved elements (UCEs; 1,382,620 bp from 1,969 loci) was used to infer a fully-resolved and well-supported maximum likelihood phylogeny for 48 hymenopteran morphospecies. There was no significant difference in species richness between landscape categories. Pollinator communities at natural sites had higher phylogenetic complexity (X = 2.37) and functional divergence (x̄ = 0.74 ± 0.02 s.e.) than disturbed sites (X = 1.65 and x̄ = 0.6 ± 0.01 s.e.). Hairiness showed significant phylogenetic clustering (K = 0.94), whereas body size, glossa length, and loose pollen showed weaker non-random phylogenetic patterns (K between 0.3-0.5). Pollen load purity showed no association with phylogeny. The assemblage of 17 bee morphospecies comprised nine FEs: eight FEs consisted of native bees with three containing 65% of all native bee taxa. The introduced honey bee (Apis mellifera) occupied a unique FE, likely due to its different evolutionary history. Both landscape types supported six FEs each with three overlapping: two native bee FEs and the honey bee FE.

CONCLUSIONS

Bee hairiness was the only functional trait to exhibit demonstrable phylogenetic signal. Despite differences in species richness, and functional and phylogenetic diversity between orchard landscape types, both maintained equal bee FE numbers. While no native bee taxon was analogous to the honey bee FE, four native bee FEs shared the same hairiness level as honey bees. Health threats to honey bee populations in Australia will likely disrupt pollination services to apple, and other pollination-dependent food crops, given the low level of functional redundancy within the investigated pollinator assemblages.

Phylogenomic data exploration with increased sampling provides new insights into the higher-level relationships of butterflies and moths (Lepidoptera)

A robust and stable phylogenetic framework is a fundamental goal of evolutionary biology. As the third largest insect order in the world following Coleoptera and Diptera, Lepidoptera (butterflies and moths) play a central role in almost every terrestrial ecosystem as indicators of environmental change and serve as important models for biologists exploring questions related to ecology and evolutionary biology. However, for such a charismatic insect group, the higher-level phylogenetic relationships among its superfamilies are still poorly resolved. Compared to earlier phylogenomic studies, we increased taxon sampling among Lepidoptera (37 superfamilies and 68 families containing 263 taxa) and acquired a series of large amino-acid datasets from 69,680 to 400,330 for phylogenomic reconstructions. Using these datasets, we explored the effect of different taxon sampling with significant increases in the number of included genes on tree topology by considering a series of systematic errors using maximum-likelihood (ML) and Bayesian inference (BI) methods. Moreover, we also tested the effectiveness in topology robustness among the three ML-based models. The results showed that taxon sampling is an important determinant in tree robustness of accurate lepidopteran phylogenetic estimation. Long-branch attraction (LBA) caused by site-wise heterogeneity is a significant source of bias giving rise to unstable positions of ditrysian groups in phylogenomic reconstruction. Phylogenetic inference showed the most comprehensive framework to reveal the relationships among lepidopteran superfamilies, and presented some newly relationships with strong supports (Papilionoidea was sister to Gelechioidea and Immoidea was sister to Galacticoidea, respectively), but limited by taxon sampling, the relationships within the species-rich and relatively rapid radiation Ditrysia and especially Apoditrysia remain poorly resolved, which need to increase taxon sampling for further phylogenomic reconstruction. The present study demonstrates that taxon sampling is an important determinant for an accurate lepidopteran tree of life and provides some essential insights for future lepidopteran phylogenomic studies.

Phylogenomics of the geometrid tribe Palyadini (Lepidoptera: Geometridae) reveals contrasting patterns of phylogenetic signal in wing colour characters

Next generation sequencing techniques currently represent a practical and efficient way to infer robust evolutionary hypotheses. Palyadini is a small Neotropical tribe of geometrid moths composed of six genera that feature strikingly colourful wings. Here, we investigated patterns of evolution and amount of phylogenetic signal contained in various colour characters featured in the wings of members of this tribe by (i) inferring a robust phylogenetic hypothesis using ultraconserved elements (UCEs), and afterwards, (ii) mapping the morphological characters onto the molecular topology under a parsimonious ancestral character optimization. Our matrix, obtained with 60% completeness, includes 754 UCE loci and 73 taxa (64 ingroup, nine outgroup). Maximum likelihood and parsimony generated largely identical topologies with strongly supported nodes, except for one node inside the genus Opisthoxia. According to our topology, most wing colour characters are reconstructed as homoplastic, particularly at the tribe level, but five of the seven provide evidence supporting common ancestry at the genus level. Our results emphasize, once again, that no character system is infallible, and that more research is necessary to take our understanding of the evolution of wing colour in moths to a level comparable with the knowledge we have for butterflies.

Phylogenomic analysis and molecular identification of true fruit flies

The family Tephritidae in the order Diptera, known as true fruit flies, are agriculturally important insect pests. However, the phylogenetic relationships of true fruit flies, remain controversial. Moreover, rapid identification of important invasive true fruit flies is essential for plant quarantine but is still challenging. To this end, we sequenced the genome of 16 true fruit fly species at coverage of 47-228×. Together with the previously reported genomes of nine species, we reconstructed phylogenetic trees of the Tephritidae using benchmarking universal single-copy ortholog (BUSCO), ultraconserved element (UCE) and anchored hybrid enrichment (AHE) gene sets, respectively. The resulting trees of 50% taxon-occupancy dataset for each marker type were generally congruent at 88% nodes for both concatenation and coalescent analyses. At the subfamily level, both Dacinae and Trypetinae are monophyletic. At the species level, Bactrocera dorsalis is more closely related to Bactrocera latifrons than Bactrocera tryoni. This is inconsistent with previous conclusions based on mitochondrial genes but consistent with recent studies based on nuclear data. By analyzing these genome data, we screened ten pairs of species-specific primers for molecular identification of ten invasive fruit flies, which PCR validated. In summary, our work provides draft genome data of 16 true fruit fly species, addressing the long-standing taxonomic controversies and providing species-specific primers for molecular identification of invasive fruit flies.

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.

Hivanua, a new genus of harmochirine jumping spiders from the Marquesas Islands (Araneae, Salticidae, Harmochirina)

The genus Hivanuagen. nov. is established for the harmochirine jumping spiders of the Marquesas Islands, formerly placed in Habronattus F.O. Pickard-Cambridge, 1901 and Havaika Prószyński, 2002. The type species, Hivanuatekaosp. nov. is described, and five species described by Berland are re-illustrated and moved into the genus: Hivanuaflavipes (Berland, 1933), comb. nov., Hivanuanigrescens (Berland, 1933), comb. nov., Hivanuanigrolineata (Berland, 1933), comb. nov., Hivanuarufescens (Berland, 1934), comb. nov., and Hivanuatriangulifera (Berland, 1933), comb. nov. The female epigyne is much like that of Habronattus, Bianor Peckham & Peckham, 1896, and other harmochirines, with a centrally placed coupling pocket and two atria with crescent-shaped edges. The terminal apophysis of the male palp, which is variable throughout the pellenine subgroup of the Harmochirina, is absent in H.rufescens but present in H.tekaosp. nov., in which it is elbowed much as in Habronattus. These Pacific Island harmochirines, like the Havaika of Hawaii, appear to be largely foliage dwellers, unlike most of their continental relatives.

Ghatippuspaschima, a new species and genus of plexippine jumping spider from the Western Ghats of India (Salticidae, Plexippini, Plexippina)

We propose a new genus of plexippine jumping spiders from the Western Ghats of India based on the new species Ghatippuspaschimagen. et sp. nov. While it bears a superficial resemblance to Pancorius in body form and Hyllus in membrane-bearing embolus, our UCE phylogenomic data-the first to resolve broad relationships within the Plexippina-as well as morphological features justify its status as a new genus. In addition to the molecular data and morphological descriptions, we provide photographs of living specimens of Ghatippuspaschimagen. et sp. nov. and information on their natural history.

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

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

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.

Phylogenomics reveals extreme gene tree discordance in a lineage of dominant trees: hybridization, introgression, and incomplete lineage sorting blur deep evolutionary relationships despite clear species groupings in Eucalyptus subgenus Eudesmia

Eucalypts are a large and ecologically important group of plants on the Australian continent, and understanding their evolution is important in understanding evolution of the unique Australian flora. Previous phylogenies using plastome DNA, nuclear-ribosomal DNA, or random genome-wide SNPs, have been confounded by limited genetic sampling or by idiosyncratic biological features of the eucalypts, including widespread plastome introgression. Here we present phylogenetic analyses of Eucalyptus subgenus Eudesmia (22 species from western, northern, central and eastern Australia), in the first study to apply a target-capture sequencing approach using custom, eucalypt-specific baits (of 568 genes) to a lineage of Eucalyptus. Multiple accessions of all species were included, and target-capture data were supplemented by separate analyses of plastome genes (average of 63 genes per sample). Analyses revealed a complex evolutionary history likely shaped by incomplete lineage sorting and hybridization. Gene tree discordance generally increased with phylogenetic depth. Species, or groups of species, toward the tips of the tree are mostly supported, and three major clades are identified, but the branching order of these clades cannot be confirmed with confidence. Multiple approaches to filtering the nuclear dataset, by removing genes or samples, could not reduce gene tree conflict or resolve these relationships. Despite inherent complexities in eucalypt evolution, the custom bait kit devised for this research will be a powerful tool for investigating the evolutionary history of eucalypts more broadly.

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.

Prospecting for Zoonotic Pathogens by Using Targeted DNA Enrichment

More than 60 zoonoses are linked to small mammals, including some of the most devastating pathogens in human history. Millions of museum-archived tissues are available to understand natural history of those pathogens. Our goal was to maximize the value of museum collections for pathogen-based research by using targeted sequence capture. We generated a probe panel that includes 39,916 80-bp RNA probes targeting 32 pathogen groups, including bacteria, helminths, fungi, and protozoans. Laboratory-generated, mock-control samples showed that we are capable of enriching targeted loci from pathogen DNA 2,882‒6,746-fold. We identified bacterial species in museum-archived samples, including Bartonella, a known human zoonosis. These results showed that probe-based enrichment of pathogens is a highly customizable and efficient method for identifying pathogens from museum-archived tissues.

Low hybridization temperatures improve target capture success of invertebrate loci: a case study of leaf-footed bugs (Hemiptera: Coreoidea)

Target capture is widely used in phylogenomic, ecological and functional genomic studies. Bait sets that allow capture from a diversity of species can be advantageous, but high-sequence divergence from baits can limit yields. Currently, only four experimental comparisons of a critical target capture parameter, hybridization temperature, have been published. These have been in vertebrates, where bait divergences are typically low, and none include invertebrates where bait-target divergences may be higher. Most invertebrate capture studies use a fixed, high hybridization temperature to maximize the proportion of on-target data, but many report low locus recovery. Using leaf-footed bugs (Hemiptera: Coreoidea), we investigate the effect of hybridization temperature on capture success of ultraconserved elements targeted by (i) baits developed from divergent hemipteran genomes and (ii) baits developed from less divergent coreoid transcriptomes. Lower temperatures generally resulted in more contigs and improved recovery of targets despite a lower proportion of on-target reads, lower read depth and more putative paralogues. Hybridization temperatures had less of an effect when using transcriptome-derived baits, which is probably due to lower bait-target divergences and greater bait tiling density. Thus, accommodating low hybridization temperatures during target capture can provide a cost-effective, widely applicable solution to improve invertebrate locus recovery.

Museomics and phylogenomics of lovebirds (Psittaciformes, Psittaculidae, Agapornis) using low-coverage whole-genome sequencing

Natural history collections contain specimens that provide important insights into studies of ecology and evolution. With the advancement of high-throughput sequencing, historical DNA (hDNA) from museum specimens has become a valuable source of genomic data to study the evolutionary history of organisms. Low-coverage whole genome sequencing (WGS) has been increasingly applied to museum specimens for analyzing organelle genomes, but is still uncommon for genotyping the nuclear DNA fraction. In this study, we applied low-coverage WGS to phylogenomic analyses of parrots in the genus Agapornis by including both modern samples and historical specimens of ∼100-year-old. Agapornis are small-sized African and Malagasy parrots with diverse characters. Earlier phylogenetic studies failed to resolve the positions of some key lineages, prohibiting a robust interpretation of the biogeography and evolution of these African parrots. Here, we demonstrated the use of low-coverage WGS for generating both mitochondrial and nuclear genomic data, and evaluated data quality differences between modern and historical samples. Our resolved Agapornis phylogeny indicates the ancestor of Agapornis likely colonized Madagascar from Australasia by trans-oceanic dispersal events before dispersing to the African continent. Genome-wide SNPs also allowed us to identify the parental origins of hybrid Agapornis individuals. This study demonstrates the potential of applying low-coverage WGS to phylogenomics and population genomics analyses and illustrates how including historical museum specimens can address outstanding questions regarding the evolutionary history of contemporary lineages.

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.

Target enrichment of long open reading frames and ultraconserved elements to link microevolution and macroevolution in non-model organisms

Despite the increasing accessibility of high-throughput sequencing, obtaining high-quality genomic data on non-model organisms without proximate well-assembled and annotated genomes remains challenging. Here, we describe a workflow that takes advantage of distant genomic resources and ingroup transcriptomes to select and jointly enrich long open reading frames (ORFs) and ultraconserved elements (UCEs) from genomic samples for integrative studies of microevolutionary and macroevolutionary dynamics. This workflow is applied to samples of the African unionid bivalve tribe Coelaturini (Parreysiinae) at basin and continent-wide scales. Our results indicate that ORFs are efficiently captured without prior identification of intron-exon boundaries. The enrichment of UCEs was less successful, but nevertheless produced substantial data sets. Exploratory continent-wide phylogenetic analyses with ORF supercontigs (>515,000 parsimony informative sites) resulted in a fully resolved phylogeny, the backbone of which was also retrieved with UCEs (>11,000 informative sites). Variant calling on ORFs and UCEs of Coelaturini from the Malawi Basin produced ~2000 SNPs per population pair. Estimates of nucleotide diversity and population differentiation were similar for ORFs and UCEs. They were low compared to previous estimates in molluscs, but comparable to those in recently diversifying Malawi cichlids and other taxa at an early stage of speciation. Skimming off-target sequence data from the same enriched libraries of Coelaturini from the Malawi Basin, we reconstructed the maternally-inherited mitogenome, which displays the gene order inferred for the most recent common ancestor of Unionidae. Overall, our workflow and results provide exciting perspectives for integrative genomic studies of microevolutionary and macroevolutionary dynamics in non-model organisms.

Evolution of stridulatory mechanisms: vibroacoustic communication may be common in leaf-footed bugs and allies (Heteroptera: Coreoidea)

Intra- and interspecific communication is crucial to fitness via its role in facilitating mating, territoriality and defence. Yet, the evolution of animal communication systems is puzzling-how do they originate and change over time? Studying stridulatory morphology provides a tractable opportunity to deduce the origin and diversification of a communication mechanism. Stridulation occurs when two sclerotized structures rub together to produce vibratory and acoustic (vibroacoustic) signals, such as a cricket 'chirp'. We investigated the evolution of stridulatory mechanisms in the superfamily Coreoidea (Hemiptera: Heteroptera), a group of insects known for elaborate male fighting behaviours and enlarged hindlegs. We surveyed a large sampling of taxa and used a phylogenomic dataset to investigate the evolution of stridulatory mechanisms. We identified four mechanisms, with at least five evolutionary gains. One mechanism, occurring only in male Harmostini (Rhopalidae), is described for the first time. Some stridulatory mechanisms appear to be non-homoplastic apomorphies within Rhopalidae, while others are homoplastic or potentially homoplastic within Coreidae and Alydidae, respectively. We detected no losses of these mechanisms once evolved, suggesting they are adaptive. Our work sets the stage for further behavioural, evolutionary and ecological studies to better understand the context in which these traits evolve and change.

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.

Taxon-specific ultraconserved element probe design for phylogenetic analyses of scale insects (Hemiptera: Sternorrhyncha: Coccoidea)

Scale insects (Coccoidea) are morphologically specialized members of the order Hemiptera, with 56 families recognized to date. However, the phylogenetic relationships within and among families are poorly resolved. In this study, to further characterize the phylogenetic relationships among scale insects, an ultraconserved element (UCE) probe set was designed specifically for Coccoidea based on three low-coverage whole genome sequences along with three publicly available genomes. An in silico test including eight additional genomes was performed to evaluate the effectiveness of the probe set. Most scale insect lineages were recovered by the phylogenetic analysis. This study recovered the monophyly of neococcoids. The newly developed UCE probe set has the potential to reshape and improve our understanding of the phylogenetic relationships within and among families of scale insects at the genome level.

Geographic isolation drives speciation in Nearctic aphids

Across herbivorous insect clades, species richness and host-use diversity tend to positively covary. This could be because host-use divergence drives speciation, or because it raises the ecological limits on species richness. To evaluate these hypotheses, we performed phylogenetic path model analyses of the species diversity of Nearctic aphids. Here, we show that variation in the species richness of aphid clades is caused mainly by host-use divergence, whereas variation in speciation rates is caused more by divergence in non-host-related niche variables. Aphid speciation is affected by both the evolution of host and non-host-related niche components, but the former is largely caused by the latter. Thus, our analyses suggest that host-use divergence can both raise the ecological limits on species richness and drive speciation, although in the latter case, host-use divergence tends to be a step along the causal path leading from non-host-related niche evolution to speciation.

Using phylogenetic path modelling, geographic isolation is revealed to be the main driver of Nearctic aphid speciation, with subsequent host-use divergence raising the limits on clade species richness.

Evolutionary history of the Australasian Scirtinae (Scirtidae; Coleoptera) inferred from ultraconserved elements

The Scirtidae Fleming, 1821 has been identified as one of the earliest diverging groups of Polyphagan beetles and is particularly speciose in Australia. However, very little is known about the origin of the Australian scirtids and there is a need for a robust, well-supported phylogeny to guide the genus and species descriptions and understand the relationships among taxa. In this study we carried out a phylogenetic analysis of the Australian Scirtinae Fleming, 1821, using DNA sequence data from ultraconserved elements (UCEs) and included representative taxa from New Zealand, New Caledonia, South America, South Africa and Eurasia in the analysis. Bayesian analyses of a concatenated dataset from 79 taxa recovered four major Southern Hemisphere groupings and two Australian–Eurasian groupings. The Veronatus group mainly consisted of genera from New Zealand, with the three Australian representatives only distantly related to each other. Relaxed molecular clock analyses, using the estimated age of the crown node of the Polyphaga for calibration, support a Gondwanan history for four of the groups of Australian Scirtinae and a northern origin for two groups. Our results highlight the value of commercially available UCEs for resolving the phylogenetic history of ancient groups of Coleoptera.

Detection of F1 hybrids from single-genome data reveals frequent hybridization in Hymenoptera and particularly ants

Hybridization occupies a central role in many fundamental evolutionary processes, such as speciation or adaptation. Yet, despite its pivotal importance in evolution, little is known about the actual prevalence and distribution of current hybridization across the tree of life. Here we develop and implement a new statistical method enabling the detection of F1 hybrids from single-individual genome sequencing data. Using simulations and sequencing data from known hybrid systems, we first demonstrate the specificity of the method, and identify its statistical limits. Next, we showcase the method by applying it to available sequencing data from more than 1,500 species of Arthropods, including Hymenoptera, Hemiptera, Coleoptera, Diptera, and Archnida. Among these taxa, we find Hymenoptera, and especially ants, to display the highest number of candidate F1 hybrids, suggesting higher rates of recent hybridization between previously isolated gene pools in these groups. The prevalence of F1 hybrids was heterogeneously distributed across ants, with taxa including many candidates tending to harbor specific ecological and life-history traits. This work shows how large-scale genomic comparative studies of recent hybridization can be implemented, uncovering the determinants of first-generation hybridization across whole taxa.

Phylogenomics and historical biogeography of the cleptoparasitic bee genus Nomada (Hymenoptera: Apidae) using ultraconserved elements

The genus Nomada Scopoli (Hymenoptera: Apidae) is the largest genus of brood parasitic bees with nearly 800 species found across the globe and in nearly all biogeographic realms except Antarctica. There is no previous molecular phylogeny focused on Nomada despite their high species abundance nor is there an existing comprehensive biogeography for the genus. Using ultraconserved element (UCE) phylogenomic data, we constructed the first molecular phylogeny for the genus Nomada and tested the monophyly of 16 morphologically established species groups. We also estimated divergence dates using fossil calibration points and inferred the origin and of this genus around the globe. Our phylogeny recovered 14 of the 16 previously established species groups as monophyletic. The superba and ruficornis groups, however, were recovered as non-monophyletic and need to be re-evaluated using morphology. Divergence dating and historic biogeographic analyses performed on the phylogenetic reconstruction indicates that Nomada most likely originated in the Holarctic ∼65 Mya. Geodispersal into the southern hemisphere occurred three times; once during the Eocene into the Afrotropics, once during the Oligocene into the Neotropics, and once during the Miocene into Australasia. Geodispersal across the Holarctic was most frequent and occurred repeatedly throughout the Cenozoic era, using the De Geer, Thulean, and the Bering Land Bridges. This is the first instance of a bee using both the Thulean and De Geer land bridges and has implications of how early bee species dispersed throughout the Palearctic in the late Cretaceous and early Paleogene.

Genome and cuticular hydrocarbon-based species delimitation shed light on potential drivers of speciation in a Neotropical ant species complex

Geographic separation that leads to the evolution of reproductive isolation between populations generally is considered the most common form of speciation. However, speciation may also occur in the absence of geographic barriers due to phenotypic and genotypic factors such as chemical cue divergence, mating signal divergence, and mitonuclear conflict. Here, we performed an integrative study based on two genome-wide techniques (3RAD and ultraconserved elements) coupled with cuticular hydrocarbon (CHC) and mitochondrial (mt) DNA sequence data, to assess the species limits within the Ectatomma ruidum species complex, a widespread and conspicuous group of Neotropical ants for which heteroplasmy (i.e., presence of multiple mtDNA variants in an individual) has been recently discovered in some populations from southeast Mexico. Our analyses indicate the existence of at least five distinct species in this complex: two widely distributed across the Neotropics, and three that are restricted to southeast Mexico and that apparently have high levels of heteroplasmy. We found that species boundaries in the complex did not coincide with geographic barriers. We therefore consider possible roles of alternative drivers that may have promoted the observed patterns of speciation, including mitonuclear incompatibility, CHC differentiation, and colony structure. Our study highlights the importance of simultaneously assessing different sources of evidence to disentangle the species limits of taxa with complicated evolutionary histories.

Phylogenomics of paleoendemic lampshade spiders (Araneae, Hypochilidae, Hypochilus), with the description of a new species from montane California

Hypochilus is a relictual lineage of Nearctic spiders distributed disjunctly across the United States in three montane regions (California, southern Rocky Mountains, southern Appalachia). Phylogenetic resolution of species relationships in Hypochilus has been challenging, and conserved morphology coupled with extreme genetic divergence has led to uncertain species limits in some complexes. Here, Hypochilus interspecies relationships have been reconstructed and cryptic speciation more critically evaluated using a combination of ultraconserved elements, mitochondrial CO1 by-catch, and morphology. Phylogenomic data strongly support the monophyly of regional clades and support a ((California, Appalachia), southern Rocky Mountains) topology. In Appalachia, five species are resolved as four lineages (H.thorelli Marx, 1888 and H.coylei Platnick, 1987 are clearly sister taxa), but the interrelationships of these four lineages remain unresolved. The Appalachian species H.pococki Platnick, 1987 is recovered as monophyletic but is highly genetically structured at the nuclear level. While algorithmic analyses of nuclear data indicate many species (e.g., all H.pococki populations as species), male morphology instead reveals striking stasis. Within the California clade, nuclear and mitochondrial lineages of H.petrunkevitchi Gertsch, 1958 correspond directly to drainage basins of the southern Sierra Nevada, with H.bernardino Catley, 1994 nested within H.petrunkevitchi and sister to the southernmost basin populations. Combining nuclear, mitochondrial, geographical, and morphological evidence a new species from the Tule River and Cedar Creek drainages is described, Hypochilusxomotesp. nov. We also emphasize the conservation issues that face several microendemic, habitat-specialized species in this remarkable genus.

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.

Taming the beast: a revised classification of Cortinariaceae based on genomic data

Family Cortinariaceae currently includes only one genus, Cortinarius, which is the largest Agaricales genus, with thousands of species worldwide. The species are important ectomycorrhizal fungi and form associations with many vascular plant genera from tropicals to arctic regions. Genus Cortinarius contains a lot of morphological variation, and its complexity has led many taxonomists to specialize in particular on infrageneric groups. The previous attempts to divide Cortinarius have been shown to be unnatural and the phylogenetic studies done to date have not been able to resolve the higher-level classification of the group above section level. Genomic approaches have revolutionized our view on fungal relationships and provide a way to tackle difficult groups. We used both targeted capture sequencing and shallow whole genome sequencing to produce data and to perform phylogenomic analyses of 75 single-copy genes from 19 species. In addition, a wider 5-locus analysis of 245 species, from the Northern and Southern Hemispheres, was also done. Based on our results, a classification of the family Cortinariaceae into ten genera—Cortinarius, Phlegmacium, Thaxterogaster, Calonarius, Aureonarius, Cystinarius, Volvanarius, Hygronarius, Mystinarius, and Austrocortinarius—is proposed. Seven genera, 10 subgenera, and four sections are described as new to science and five subgenera are introduced as new combinations in a new rank. In addition, 41 section names and 514 species names are combined in new genera and four lecto- and epitypes designated. The position of Stephanopus in suborder Agaricineae remains to be studied. Targeted capture sequencing is used for the first time in fungal taxonomy in Basidiomycetes. It provides a cost-efficient way to produce -omics data in species-rich groups. The -omics data was produced from fungarium specimens up to 21 years old, demonstrating the value of museum specimens in the study of the fungal tree of life. This study is the first family revision in Agaricales based on genomics data and hopefully many others will soon follow.

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.

UCEasy: A software package for automating and simplifying the analysis of ultraconserved elements (UCEs)

Background

The use of Ultraconserved Elements (UCEs) as genetic markers in phylogenomics has become popular and has provided promising results. Although UCE data can be easily obtained from targeted enriched sequencing, the protocol for in silico analysis of UCEs consist of the execution of heterogeneous and complex tools, a challenge for scientists without training in bioinformatics. Developing tools with the adoption of best practices in research software can lessen this problem by improving the execution of computational experiments, thus promoting better reproducibility.

New information

We present UCEasy, an easy-to-install and easy-to-use software package with a simple command line interface that facilitates the computational analysis of UCEs from sequencing samples, following the best practices of research software. UCEasy is a wrapper that standardises, automates and simplifies the quality control of raw reads, assembly and extraction and alignment of UCEs, generating at the end a data matrix with different levels of completeness that can be used to infer phylogenetic trees. We demonstrate the functionalities of UCEasy by reproducing the published results of phylogenomic studies of the bird genus Turdus (Aves) and of Adephaga families (Coleoptera) containing genomic datasets to efficiently extract UCEs.

Genome-wide ultraconserved elements resolve phylogenetic relationships and biogeographic history among Neotropical leaf-nosed bats in the genus Anoura (Phyllostomidae)

AnouraGray, 1838 are Neotropical nectarivorous bats and the most speciose genus within the phyllostomid subfamily Glossophaginae. However, Anoura species limits remain debated, and phylogenetic relationships remain poorly known, because previous studies used limited Anoura taxon sampling or focused primarily on higher-level relationships. Here, we conduct the first phylogenomic study of Anoura by analyzing 2039 genome-wide ultraconserved elements (UCEs) sequenced for 42 individuals from 8 Anoura species/lineages plus two outgroups. Overall, our results based on UCEs resolved relationships in the genus and supported (1) the monophyly of small-bodied Anoura species (previously genus Lonchoglossa); (2) monotypic status of A. caudifer; and (3) nested positions of "A. carishina", A. caudifer aequatoris, and A. geoffroyi peruana specimens within A. latidens, A. caudifer and A. geoffroyi, respectively (suggesting that these taxa are not distinct species). Additionally, (4) phylogenetic networks allowing reticulate edges did not explain gene tree discordance better than the species tree (without introgression), indicating that a coalescent model accounting for discordance solely through incomplete lineage sorting fit our data well. Sensitivity analyses indicated that our species tree results were not adversely affected by varying taxon sampling across loci. Tree calibration and Bayesian coalescent analyses dated the onset of diversification within Anoura to around ∼ 6-9 million years ago in the Miocene, with extant species diverging mainly within the past ∼ 4 million years. We inferred a historical biogeographical scenario for Anoura of parapatric speciation fragmenting the range of a wide-ranging ancestral lineage centered in the Central to Northern Andes, along with Pliocene-Pleistocene dispersal or founder event speciation in Amazonia and the Brazilian Atlantic forest during the last ∼ 2.5 million years.

Combining genomic, phenotypic and Sanger sequencing data to elucidate the phylogeny of the two-clawed spiders (Dionycha)

The importance of morphology in the phylogenomic era has recently gained attention, but relatively few studies have combined both types of information when inferring phylogenetic relationships. Sanger sequencing legacy data can also be important for understanding evolutionary relationships. The possibility of combining genomic, morphological and Sanger data in one analysis seems compelling, permitting a more complete sampling and yielding a comprehensive view of the evolution of a group. Here we used these three data types to elucidate the systematics and evolution of the Dionycha, a highly diverse group of spiders relatively underrepresented in phylogenetic studies. The datasets were analyzed separately and combined under different inference methods, including a novel approach for analyzing morphological matrices with commonly used evolutionary models. We tested alternative hypotheses of relationships and performed simulations to investigate the accuracy of our findings. We provide a comprehensive and thorough phylogenetic hypothesis for Dionycha that can serve as a robust framework to test hypotheses about the evolution of key characters. We also show that morphological data might have a phylogenetic impact, even when massively outweighed by molecular data. Our approach to analyze morphological data may serve as an alternative to the proposed practice of arbitrarily partitioning, weighting, and choosing between parsimony and stochastic models. As a result of our findings, we propose Trachycosmidae new rank for a group of Australian genera formerly included in Trochanteriidae and Gallieniellidae, and consider Ammoxenidae as a junior synonym of Gnaphosidae. We restore the family rank for Prodidomidae, but transfer the subfamily Molycriinae to Gnaphosidae. Drassinella is transferred to Liocranidae, Donuea to Corinnidae, and Mahafalytenus to Viridasiidae.

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

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

New targets acquired: Improving locus recovery from the Angiosperms353 probe set

PREMISE

Universal target enrichment kits maximize utility across wide evolutionary breadth while minimizing the number of baits required to create a cost-efficient kit. The Angiosperms353 kit has been successfully used to capture loci throughout the angiosperms, but the default target reference file includes sequence information from only 6-18 taxa per locus. Consequently, reads sequenced from on-target DNA molecules may fail to map to references, resulting in fewer on-target reads for assembly, and reducing locus recovery.

METHODS

We expanded the Angiosperms353 target file, incorporating sequences from 566 transcriptomes to produce a 'mega353' target file, with each locus represented by 17-373 taxa. This mega353 file is a drop-in replacement for the original Angiosperms353 file in HybPiper analyses. We provide tools to subsample the file based on user-selected taxon groups, and to incorporate other transcriptome or protein-coding gene data sets.

RESULTS

Compared to the default Angiosperms353 file, the mega353 file increased the percentage of on-target reads by an average of 32%, increased locus recovery at 75% length by 49%, and increased the total length of the concatenated loci by 29%.

DISCUSSION

Increasing the phylogenetic density of the target reference file results in improved recovery of target capture loci. The mega353 file and associated scripts are available at: https://github.com/chrisjackson-pellicle/NewTargets.

Evolution of host use in fungivorous ciid beetles (Coleoptera: Ciidae): Molecular phylogeny focusing on Japanese taxa

Consumer-resource interactions between trophic levels are ubiquitous and important factors in shaping the diversity of insects. However, dietary patterns such as host specificity and conservatism have been insufficiently examined in fungivorous insects. Here we reconstructed the evolutionary history of host use in fungivorous ciid beetles (Coleoptera: Ciidae) and tested for host conservatism. Phylogenetic relationships among 49 species from Japan were inferred by using a large sequence data set from ultraconserved elements (UCEs). In addition, sequences of three genes (COI, 28S rRNA, 18S rRNA) were analyzed to reconstruct the phylogeny for 130 OTUs from a broader range of taxa and geographic regions using the UCE tree as a backbone topology. We found that Ciini and Orophiini are not recovered as reciprocally monophyletic groups. As previously suggested, the largest genus Cis Latreille was also not monophyletic. Ancestral-state reconstruction of host use in both datasets showed that Ciidae species were clustered by host-use group across the tree. This pattern was confirmed by the significantly lower transition rate compared with expectations under the random shift hypothesis. The observed conservatism in host use implied these beetles possess unique adaptations to specific fungal taxa, just as herbivorous insects are adapted to specific plant taxa.

New alignment-based sequence extraction software (ALiBaSeq) and its utility for deep level phylogenetics

Despite many bioinformatic solutions for analyzing sequencing data, few options exist for targeted sequence retrieval from whole genomic sequencing (WGS) data with the ultimate goal of generating a phylogeny. Available tools especially struggle at deep phylogenetic levels and necessitate amino-acid space searches, which may increase rates of false positive results. Many tools are also difficult to install and may lack adequate user resources. Here, we describe a program that uses freely available similarity search tools to find homologs in assembled WGS data with unparalleled freedom to modify parameters. We evaluate its performance compared to other commonly used bioinformatics tools on two divergent insect species (>200 My) for which annotated genomes exist, and on one large set each of highly conserved and more variable loci. Our software is capable of retrieving orthologs from well-curated or unannotated, low or high depth shotgun, and target capture assemblies as well or better than other software as assessed by recovering the most genes with maximal coverage and with a low rate of false positives throughout all datasets. When assessing this combination of criteria, ALiBaSeq is frequently the best evaluated tool for gathering the most comprehensive and accurate phylogenetic alignments on all types of data tested. The software (implemented in Python), tutorials, and manual are freely available at https://github.com/AlexKnyshov/alibaseq.

A phylogenetic and taxonomic review of baviine jumping spiders (Araneae, Salticidae, Baviini)

The systematics and taxonomy of the tropical Asian jumping spiders of the tribe Baviini is reviewed, with a molecular phylogenetic study (UCE sequence capture, traditional Sanger sequencing) guiding a reclassification of the group’s genera. The well-studied members of the group are placed into six genera: Bavia Simon, 1877, Indopadilla Caleb & Sankaran, 2019, Padillothorax Simon, 1901, Piranthus Thorell, 1895, Stagetillus Simon, 1885, and one new genus, Maripanthus Maddison, gen. nov. The identity of Padillothorax is clarified, and Bavirecta Kanesharatnam & Benjamin, 2018 synonymized with it. Hyctiota Strand, 1911 is synonymized with Stagetillus. The molecular phylogeny divides the baviines into three clades, the Piranthus clade with a long embolus (Piranthus, Maripanthus), the genus Padillothorax with a flat body and short embolus, and the Bavia clade with a higher body and (usually) short embolus (remaining genera). In general, morphological synapomorphies support or extend the molecularly delimited groups. Eighteen new species are described: Bavianessagyna, Indopadillabamilin, I.kodagura, I.nesinor, I.redunca, I.redynis, I.sabivia, I.vimedaba, Maripanthusdraconis (type species of Maripanthus), M.jubatus, M.reinholdae, Padillothoraxbadut, P.mulu, Piranthusapi, P.bakau, P.kohi, P.mandai, and Stagetillusirri, all sp. nov., with taxonomic authority W. Maddison. The distinctions between baviines and the astioid Nungia Żabka, 1985 are reviewed, leading to four species being moved into Nungia from Bavia and other genera. Fifteen new combinations are established: Baviamaurerae (Freudenschuss & Seiter, 2016), Indopadillaannamita (Simon, 1903), I.kahariana (Prószyński & Deeleman-Reinhold, 2013), I.sonsorol (Berry, Beatty & Prószyński, 1997), I.suhartoi (Prószyński & Deeleman-Reinhold, 2013), Maripanthusmenghaiensis (Cao & Li, 2016), M.smedleyi (Reimoser, 1929), Nungiahatamensis (Thorell, 1881), N.modesta (Keyserling, 1883), N.papakula (Strand, 1911), N.xiaolonghaensis (Cao & Li, 2016), Padillothoraxcasteti (Simon, 1900), P.exilis (Cao & Li, 2016), P.flavopunctus (Kanesharatnam & Benjamin, 2018), Stagetillusbanda (Strand, 1911), all comb. nov. One combination is restored, Baviacapistrata (C. L. Koch, 1846). Five of these new or restored combinations correct previous errors of placing species in genera that have superficially similar palps but extremely different body forms, in fact belonging in distantly related tribes, emphasizing that the general shape of male palps should be used with caution in determining relationships. A little-studied genus, Padillothorus Prószyński, 2018, is tentatively assigned to the Baviini. Ligdus Thorell, 1895 is assigned to the Ballini.

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 & Derkarabetian, fam. nov., consisting of Lomanella Pocock, 1903, and a newly described genus Abaddon Derkarabetian & 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 & Giribet, gen. nov., with one species, T. fadriquei Derkarabetian, Prieto & 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

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.

Phylogenomic Analysis of Ultraconserved Elements Resolves the Evolutionary and Biogeographic History of Segmented Trapdoor Spiders

The segmented trapdoor spiders (Liphistiidae) are the sole surviving family of the suborder Mesothelae, which forms the sister lineage to all other living spiders. Liphistiids have retained a number of plesiomorphic traits and their present-day distribution is limited to East and Southeast Asia. Studying this group has the potential to shed light on the deep evolutionary history of spiders, but the phylogeny and divergence times of the family have not been resolved with confidence. We performed phylogenomic and molecular dating analyses of 2,765 ultraconserved element loci from 185 liphistiid taxa. Our analyses show that the crown group of Liphistiidae appeared in the mid-Cretaceous at 102 Ma (95% credibility interval 92-113 Ma), but it was not until the Neogene that much of the diversification within the family occurred in mainland Southeast and East Asia. This diversification was coincident with tectonic events such as the extension of the East Asian continental margin, as well as geological upheavals in Indochina induced by the collision between India and Asia. Our study highlights the important role of major tectonic events in shaping the evolutionary history, present-day diversity, and geographical distribution of mesothele and liphistiid spiders.

Comprehensive phylogeny of Myrmecocystus honey ants highlights cryptic diversity and infers evolution during aridification of the American Southwest

The New World ant genus Myrmecocystus Wesmael, 1838 (Formicidae: Formicinae: Lasiini) is endemic to arid and semi-arid habitats of the western United States and Mexico. Several intriguing life history traits have been described for the genus, the best-known of which are replete workers, that store liquified food in their largely expanded crops and are colloquially referred to as "honeypots". Despite their interesting biology and ecological importance for arid ecosystems, the evolutionary history of Myrmecocystus ants is largely unknown and the current taxonomy presents an unsatisfactory systematic framework. We use ultraconserved elements to infer the evolutionary history of Myrmecocystus ants and provide a comprehensive, dated phylogenetic framework that clarifies the molecular systematics within the genus with high statistical support, reveals cryptic diversity, and reconstructs ancestral foraging activity. Using maximum likelihood, Bayesian and species tree approaches on a data set of 134 ingroup specimens (including samples from natural history collections and type material), we recover largely identical topologies that leave the position of only few clades uncertain and cover the intra- and interspecific variation of 28 of the 29 described and six undescribed species. In addition to traditional support values, such as bootstrap and posterior probability, we quantify genealogical concordance to estimate the effects of conflicting evolutionary histories on phylogenetic inference. Our analyses reveal that the current taxonomic classification of the genus is inconsistent with the molecular phylogenetic inference, and we identify cryptic diversity in seven species. Divergence dating suggests that the split between Myrmecocystus and its sister taxon Lasius occurred in the early Miocene. Crown group Myrmecocystus started diversifying about 14.08 Ma ago when the gradual aridification of the southwestern United States and northern Mexico led to formation of the American deserts and to adaptive radiations of many desert taxa.

Ultraconserved element bait set for trypanosomatida target enrichment and phylogenetics

Lack of knowledge of taxonomic biodiversity and reliable genetic markers in Trypanosomatidae limit our understanding of their phylogenetic relationships. Ultraconserved elements (UCEs) have improved phylogenetic analyses and inferences in many vertebrate and invertebrate taxa. However, it is unknown whether protozoans have these markers, their abundance, and if these could be reliably used for phylogenetics. In this study I design a target enrichment bait set for UCE loci for this group. In silico testing showed good loci recovery rates across 63 taxa and produced consistent, highly supported phylogenetic trees. This bait set adds a new resource of useful genetic markers for Trypanosomatidae phylogenetics.

Rediscovering a Forgotten System of Symbiosis: Historical Perspective and Future Potential

While the majority of symbiosis research is focused on bacteria, microbial eukaryotes play important roles in the microbiota and as pathogens, especially the incredibly diverse Fungi kingdom. The recent emergence of widespread pathogens in wildlife (bats, amphibians, snakes) and multidrug-resistant opportunists in human populations (Candida auris) has highlighted the importance of better understanding animal-fungus interactions. Regardless of their prominence there are few animal-fungus symbiosis models, but modern technological advances are allowing researchers to utilize novel organisms and systems. Here, I review a forgotten system of animal-fungus interactions: the beetle-fungus symbioses of Drugstore and Cigarette beetles with their symbiont Symbiotaphrina. As pioneering systems for the study of mutualistic symbioses, they were heavily researched between 1920 and 1970, but have received only sporadic attention in the past 40 years. Several features make them unique research organisms, including (1) the symbiont is both extracellular and intracellular during the life cycle of the host, and (2) both beetle and fungus can be cultured in isolation. Specifically, fungal symbionts intracellularly infect cells in the larval and adult beetle gut, while accessory glands in adult females harbor extracellular fungi. In this way, research on the microbiota, pathogenesis/infection, and mutualism can be performed. Furthermore, these beetles are economically important stored-product pests found worldwide. In addition to providing a historical perspective of the research undertaken and an overview of beetle biology and their symbiosis with Symbiotaphrina, I performed two analyses on publicly available genomic data. First, in a preliminary comparative genomic analysis of the fungal symbionts, I found striking differences in the pathways for the biosynthesis of two B vitamins important for the host beetle, thiamine and biotin. Second, I estimated the most recent common ancestor for Drugstore and Cigarette beetles at 8.8-13.5 Mya using sequence divergence (CO1 gene). Together, these analyses demonstrate that modern methods and data (genomics, transcriptomes, etc.) have great potential to transform these beetle-fungus systems into model systems again.