Insight from an ultraconserved element bait set designed for hemipteran phylogenetics integrated with genomic resources

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.

Lessons from assembling UCEs: A comparison of common methods and the case of Clavinomia (Halictidae)

Sequence data assembly is a foundational step in high-throughput sequencing, with untold consequences for downstream analyses. Despite this, few studies have interrogated the many methods for assembling phylogenomic UCE data for their comparative efficacy, or for how outputs may be impacted. We study this by comparing the most commonly used assembly methods for UCEs in the under-studied bee lineage Nomiinae and a representative sampling of relatives. Data for 63 UCE-only and 75 mixed taxa were assembled with five methods, including ABySS, HybPiper, SPAdes, Trinity and Velvet, and then benchmarked for their relative performance in terms of locus capture parameters and phylogenetic reconstruction. Unexpectedly, Trinity and Velvet trailed the other methods in terms of locus capture and DNA matrix density, whereas SPAdes performed favourably in most assessed metrics. In comparison with SPAdes, the guided-assembly approach HybPiper generally recovered the highest quality loci but in lower numbers. Based on our results, we formally move Clavinomia to Dieunomiini and render Epinomia once more a subgenus of Dieunomia. We strongly advise that future studies more closely examine the influence of assembly approach on their results, or, minimally, use better-performing assembly methods such as SPAdes or HybPiper. In this way, we can move forward with phylogenomic studies in a more standardized, comparable manner.

A combined phylogenetic strategy illuminates the evolution of Goniodorididae nudibranchs (Mollusca, Gastropoda, Heterobranchia)

Goniodorididae is a family of small dorid nudibranchs distributed worldwide that feed on entoprocts, ascidians, and bryozoans. The evolutionary relationships between its taxa have been uncertain due to the limited taxa available for phylogenetic analyses; some genera being paraphyletic. The family includes a remarkable number of synonymized genera in which the species richness is unequally distributed, while some genera have dozens of species others are monospecific. Some clades are very uniform morphologically while others are considered highly variable. To increase backbone phylogenetic resolution a target enrichment approach of ultra-conserved elements was aimed at representative Goniodorididae species for the first time. Additionally, we increase species representation by including mitochondrial markers cytochrome c oxidase subunit I and ribosomal RNA 16S as well as nuclear Histone 3 and ribosomal RNA 18S from 109 Goniodorididae species, out of approximately 160 currently valid species. Maximum likelihood and Bayesian inference analyses were performed to infer the phylogeny of the family. As a result, two subfamilies and eleven genera were elucidated. The synonymized genera Bermudella, Cargoa, and Ceratodoris are here resurrected and a new genus, Naisdoris gen. nov., is described. The clades included taxa with shared prey preference, showing that trophic behavior could have driven species evolution and morphological uniqueness within the family Goniodorididae.

Parallel Losses of Blue Opsin Correlate with Compensatory Neofunctionalization of UV-Opsin Gene Duplicates in Aphids and Planthoppers

Expanding on previous efforts to survey the visual opsin repertoires of the Hemiptera, this study confirms that homologs of the UV- and LW-opsin subfamilies are conserved in all Hemiptera, while the B-opsin subfamily is missing from the Heteroptera and subgroups of the Sternorrhyncha and Auchenorrhyncha, i.e., aphids (Aphidoidea) and planthoppers (Fulgoroidea), respectively. Unlike in the Heteroptera, which are characterized by multiple independent expansions of the LW-opsin subfamily, the lack of B-opsin correlates with the presence of tandem-duplicated UV-opsins in aphids and planthoppers. Available data on organismal wavelength sensitivities and retinal gene expression patterns lead to the conclusion that, in both groups, one UV-opsin paralog shifted from ancestral UV peak sensitivity to derived blue sensitivity, likely compensating for the lost B-opsin. Two parallel bona fide tuning site substitutions compare to 18 non-corresponding amino acid replacements in the blue-shifted UV-opsin paralogs of aphids and planthoppers. Most notably, while the aphid blue-shifted UV-opsin clade is characterized by a replacement substitution at one of the best-documented UV/blue tuning sites (Rhodopsin site 90), the planthopper blue-shifted UV-opsin paralogs retained the ancestral lysine at this position. Combined, the new findings identify aphid and planthopper UV-opsins as a new valuable data sample for studying adaptive opsin evolution.

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.

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.

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.

New Cretaceous Bugs from Northeastern China Imply the Systematic Position of Pachymeridiidae (Hemiptera: Heteroptera)

Varicapitatus sinuolatus gen. et sp. nov. is described from the Yixian Formation of Northeastern China. Based on the new specimens, this study discusses the morphological characteristics and taxonomic position of Pachymeridiidae: Pachymeridiidae belongs to Pentatomomorpha and is more closely related to Lygaeoidea. In Heteroptera, the costal fracture of the forewing is a homoplastic characteristic, already evolved independently several times among most taxa before the Early Cretaceous. The pulvillus under the claw in Pentatomomorpha is also a homoplastic characteristic. In Pachymeridiidae, the forewing costal fracture and absence of pulvillus can be recognized as an independent evolution or convergence, implying that pachymeridiids may have different flight and crawling behaviors distinct from other Pentatomomorpha groups.

Phylogeny of Coreoidea based on mitochondrial genomes show the paraphyly of Coreidae and Alydidae

Coreoidea (Insecta: Hemiptera: Heteroptera) is a widely distributed and agriculturally important bugs. However, the phylogeny of Coreoidea lacked consensus on higher-level relationships and several studies by comparative morphological characters and molecular data suggested the non-monophyly of two families: Coreidae and Alydidae. The mitochondrial genome (mitogenome) has long been thought to be a significant marker to understand phylogenetic relationships, but the mitogenome in Alydidae is scarce to date. In the present study, we gathered the mitogenomes of 28 species from four families of Coreoidea excluding Hyocephalidae (Alydidae, Coreidae, Rhopalidae, and Stenocephalidae), including four newly sequenced mitogenomes of Alydidae, and conducted mitogenomic organization and phylogenetic studies. We used maximum likelihood and Bayesian inference methods to infer the higher-level phylogeny from the perspective of mitogenomes, primarily to investigate the phylogenetic relationship betweeen Coreidae and Alydidae. We add evidence that neither Alydidae nor Coreidae are monophyletic based on mitogenomes. Newly sequenced mitogenomes of Alydidae have traditional gene structure and gene rearrangement was not found. Alydinae was always recovered as closely related to Pseudophloeinae of the coreid subfamily with high support. The placement of the coreid subfamily Hydarinae and alydid subfamily Micrelytrinae are unstable depending on approach used. In terms of the length and nucleotide composition of the protein coding genes in mitogenomes, Pseudophloeinae and Hydarinae of coreid were more similar to Alydidae. The unsettled classification issues of Coreidae and Alydidae by mitogenomes were demonstrated in this work, indicating that further study is needed.

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.

Diversification of the phytophagous lineages of true bugs (Insecta: Hemiptera: Heteroptera) shortly after that of the flowering plants

More than 95% of phytophagous true bug (Hemiptera: Heteroptera) species belong to four superfamilies: Miroidea (Cimicomorpha), Pentatomoidea, Coreoidea, and Lygaeoidea (all Pentatomomorpha). These iconic groups of highly diverse, overwhelmingly phytophagous insects include several economically prominent agricultural and silvicultural pest species, though their evolutionary history has not yet been well resolved. In particular, superfamily- and family-level phylogenetic relationships of these four lineages have remained controversial, and the divergence times of some crucial nodes for phytophagous true bugs have hitherto been little known, which hampers a better understanding of the evolutionary processes and patterns of phytophagous insects. In the present study, we used 150 species and concatenated nuclear and mitochondrial protein-coding genes and rRNA genes to infer the phylogenetic relationships within the Terheteroptera (Cimicomorpha + Pentatomomorpha) and estimated their divergence times. Our results support the monophyly of Cimicomorpha, Pentatomomorpha, Miroidea, Pentatomoidea, Pyrrhocoroidea, Coreoidea, and Lygaeoidea. The phylogenetic relationships across phytophagous lineages are largely congruent at deep nodes across the analyses based on different datasets and tree-reconstructing methods with just a few exceptions. Estimated divergence times and ancestral state reconstructions for feeding habit indicate that phytophagous true bugs explosively radiated in the Early Cretaceous-shortly after the angiosperm radiation-with the subsequent diversification of the most speciose clades (Mirinae, Pentatomidae, Coreinae, and Rhyparochromidae) in the Late Cretaceous.

Divergence time estimation of Galliformes based on the best gene shopping scheme of ultraconserved elements

BACKGROUND

Divergence time estimation is fundamental to understanding many aspects of the evolution of organisms, such as character evolution, diversification, and biogeography. With the development of sequence technology, improved analytical methods, and knowledge of fossils for calibration, it is possible to obtain robust molecular dating results. However, while phylogenomic datasets show great promise in phylogenetic estimation, the best ways to leverage the large amounts of data for divergence time estimation has not been well explored. A potential solution is to focus on a subset of data for divergence time estimation, which can significantly reduce the computational burdens and avoid problems with data heterogeneity that may bias results.

RESULTS

In this study, we obtained thousands of ultraconserved elements (UCEs) from 130 extant galliform taxa, including representatives of all genera, to determine the divergence times throughout galliform history. We tested the effects of different "gene shopping" schemes on divergence time estimation using a carefully, and previously validated, set of fossils. Our results found commonly used clock-like schemes may not be suitable for UCE dating (or other data types) where some loci have little information. We suggest use of partitioning (e.g., PartitionFinder) and selection of tree-like partitions may be good strategies to select a subset of data for divergence time estimation from UCEs. Our galliform time tree is largely consistent with other molecular clock studies of mitochondrial and nuclear loci. With our increased taxon sampling, a well-resolved topology, carefully vetted fossil calibrations, and suitable molecular dating methods, we obtained a high quality galliform time tree.

CONCLUSIONS

We provide a robust galliform backbone time tree that can be combined with more fossil records to further facilitate our understanding of the evolution of Galliformes and can be used as a resource for comparative and biogeographic studies in this group.

A novel exome probe set captures phototransduction genes across birds (Aves) enabling efficient analysis of vision evolution

The diversity of avian visual phenotypes provides a framework for studying mechanisms of trait diversification generally, and the evolution of vertebrate vision, specifically. Previous research has focused on opsins, but to fully understand visual adaptation, we must study the complete phototransduction cascade (PTC). Here, we developed a probe set that captures exonic regions of 46 genes representing the PTC and other light responses. For a subset of species, we directly compared gene capture between our probe set and low-coverage whole genome sequencing (WGS), and we discuss considerations for choosing between these methods. Finally, we developed a unique strategy to avoid chimeric assembly by using "decoy" reference sequences. We successfully captured an average of 64% of our targeted exome in 46 species across 14 orders using the probe set and had similar recovery using the WGS data. Compared to WGS or transcriptomes, our probe set: (1) reduces sequencing requirements by efficiently capturing vision genes, (2) employs a simpler bioinformatic pipeline by limiting required assembly and negating annotation, and (3) eliminates the need for fresh tissues, enabling researchers to leverage existing museum collections. We then utilized our vision exome data to identify positively selected genes in two evolutionary scenarios-evolution of night vision in nocturnal birds and evolution of high-speed vision specific to manakins (Pipridae). We found parallel positive selection of SLC24A1 in both scenarios, implicating the alteration of rod response kinetics, which could improve color discrimination in dim light conditions and/or facilitate higher temporal resolution.

Ultra-Conserved Elements and morphology reciprocally illuminate conflicting phylogenetic hypotheses in Chalcididae (Hymenoptera, Chalcidoidea)

Recent technical advances combined with novel computational approaches have promised the acceleration of our understanding of the tree of life. However, when it comes to hyperdiverse and poorly known groups of invertebrates, studies are still scarce. As published phylogenies will be rarely challenged by future taxonomists, careful attention must be paid to potential analytical bias. We present the first molecular phylogenetic hypothesis for the family Chalcididae, a group of parasitoid wasps, with a representative sampling (144 ingroups and seven outgroups) that covers all described subfamilies and tribes, and 82% of the known genera. Analyses of 538 Ultra-Conserved Elements (UCEs) with supermatrix (RAxML and IQTREE) and gene tree reconciliation approaches (ASTRAL, ASTRID) resulted in highly supported topologies in overall agreement with morphology but reveal conflicting topologies for some of the deepest nodes. To resolve these conflicts, we explored the phylogenetic tree space with clustering and gene genealogy interrogation methods, analyzed marker and taxon properties that could bias inferences and performed a thorough morphological analysis (130 characters encoded for 40 taxa representative of the diversity). This joint analysis reveals that UCEs enable attainment of resolution between ancestry and convergent/divergent evolution when morphology is not informative enough, but also shows that a systematic exploration of bias with different analytical methods and a careful analysis of morphological features is required to prevent publication of artifactual results. We highlight a GC content bias for maximum-likelihood approaches, an artifactual mid-point rooting of the ASTRAL tree and a deleterious effect of high percentage of missing data (>85% missing UCEs) on gene tree reconciliation methods. Based on the results we propose a new classification of the family into eight subfamilies and ten tribes that lay the foundation for future studies on the evolutionary history of Chalcididae.

Drawing the Excalibur bug from the stone: adding credibility to the double-edged sword hypothesis of coreid evolution (Hemiptera, Coreidae)

A new genus and species of exaggerated antennae Coreidae is described from Myanmar amber of the Late Cretaceous (Cenomanian stage). Ferriantennaexcaliburgen. et sp. nov. appears related to another Cretaceous coreid with exaggerated antennae, Magnusantenna Du & Chen, 2021, but can be differentiated by the fourth antennal segment which is short and paddle-like, the undulating shape of the pronotum and mesonotum, and the shorter and thicker legs. The new coreid, with elaborately formed antennae and simple hind legs instead of the typical extant coreid morphology with simple antennae and elaborately formed hind legs, begs the question: why were the elaborate features of the antennae lost in favor of ornate hind legs? Features that are large and showy are at higher risk of being attacked by predators or stuck in a poor molt and subjected to autotomy and are therefore lost at a higher rate than simple appendages. We hypothesize that because elaborate antennae play an additional significant sensory role compared to elaborate hind legs, that evolutionarily it is more costly to have elaborate antennae versus elaborate hind legs. Thus, through the millenia, as coreid evolution experimented with elaborate/ornate features, those on the antennae were likely selected against in favor of ornate hind legs.

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.

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.

Fossils reshape the Sternorrhyncha evolutionary tree (Insecta, Hemiptera)

The Sternorrhyncha, which comprise about 18,700 described recent species, is a suborder of the Hemiptera, one of big five most diverse insect orders. In the modern fauna, these tiny phytophages comprise insects of great ecological and economic importance, like aphids (Aphidomorpha), scale insects (Coccidomorpha), whiteflies (Aleyrodomorpha) and psyllids (Psylloidea). Their evolutionary history can be traced back to the Late Carboniferous, but the early stages of their evolution and diversification is poorly understood, with two known extinct groups-Pincombeomorpha and Naibiomorpha variously placed in classifications and relationships hypotheses. Most of the recent Sternorrhyncha groups radiated rapidly during the Cretaceous. Here we report the new finding of very specialised sternorrhynchans found as inclusions in mid-Cretaceous amber from Kachin state (northern Myanmar), which represent another extinct lineage within this hemipteran suborder. These fossils, proposed to be placed in a new infraorder, are revealed to be related to whiteflies and psyllids. We present, also for the first time, the results of phylogenetic analyses covering extinct and extant lineages of the Sternorrhyncha.

Population genetics of two chromatic morphs of the Chagas disease vector Rhodnius pallescens Barber, 1932 in Panamá

Rhodnius pallescens is the principal vector of Chagas disease in Panama. Recently a dark chromatic morph has been discovered in the highlands of Veraguas Province. Limited genetic studies have been conducted with regards to the population structure and dispersal potential of Triatominae vectors, particularly in R. pallescens. Next generation sequencing methods such as RADseq and complete mitochondrial DNA (mtDNA) genome sequencing have great potential for examining vector biology across space and time. Here we utilize a RADseq method (3RAD), along with complete mtDNA sequencing, to examine the population structure of the two chromatic morpho types of R. pallescens in Panama. We sequenced 105 R. pallescens samples from five localities in Panama. We generated a 2216 SNP dataset and 6 complete mtDNA genomes. RADseq showed significant differentiation among the five localities (F_CT = 0.695; P = .004), but most of this was between localities with the dark vs. light chromatic morphs (Veraguas vs. Panama Oeste). The mtDNA genomes showed a 97-98% similarity between dark and light chromatic morphs across all genes and a 502 bp insert in light morphs. Thus, both the RADseq and mtDNA data showed highly differentiated clades with essentially no gene flow between the dark and light chromatic morphs from Veraguas and central Panama respectively. We discuss the growing evidence showing clear distinctions between these two morpho types with the possibility that these are separate species, an area of research that requires further investigation. Finally, we discuss the cost-effectiveness of 3RAD which is a third of the cost compared to other RADseq methods used recently in Chagas disease vector research.

The evolution of autotomy in leaf-footed bugs

Sacrificing body parts is one of many behaviors that animals use to escape predation. This trait, termed autotomy, is classically associated with lizards. However, several other taxa also autotomize, and this trait has independently evolved multiple times throughout Animalia. Despite having multiple origins and being an iconic antipredatory trait, much remains unknown about the evolution of autotomy. Here, we combine morphological, behavioral, and genomic data to investigate the evolution of autotomy within leaf-footed bugs and allies (Insecta: Hemiptera: Coreidae + Alydidae). We found that the ancestor of leaf-footed bugs autotomized and did so slowly; rapid autotomy (<2 min) then arose multiple times. The ancestor likely used slow autotomy to reduce the cost of injury or to escape nonpredatory entrapment but could not use autotomy to escape predation. This result suggests that autotomy to escape predation is a co-opted benefit (i.e., exaptation), revealing one way that sacrificing a limb to escape predation may arise. In addition to identifying the origins of rapid autotomy, we also show that across species variation in the rates of autotomy can be explained by body size, distance from the equator, and enlargement of the autotomizable appendage.

Putting the genome in insect phylogenomics

Next-generation sequencing technologies provide a substantial increase in the size of molecular phylogenetic datasets that can be obtained for studies of insect systematics. Several new genome reduction approaches are leveraging these technologies to generate large phylogenomic datasets: targeted amplicon sequencing, target capture, and transcriptome sequencing. Although cost effective, these approaches provide limited data for questions outside of phylogenetics. For many groups of insects, sequencing the entire genome at modest coverage is feasible. Using these genomic reads, an automated Target Restricted Assembly Method (aTRAM) can use the results of blast searches to assemble thousands of single copy ortholog genes across a group of interest. These locally assembled genes can then be compiled into very large phylogenomic datasets. These genomic libraries have the advantage in that they also contain reads from the mitochondrial genome and symbiont genomes, as well the entire insect genome, and can be leveraged for additional studies beyond phylogenetics.

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

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

Off-target capture data, endosymbiont genes and morphology reveal a relict lineage that is sister to all other singing cicadas

Phylogenetic asymmetry is common throughout the tree of life and results from contrasting patterns of speciation and extinction in the paired descendant lineages of ancestral nodes. On the depauperate side of a node, we find extant ‘relict’ taxa that sit atop long, unbranched lineages. Here, we show that a tiny, pale green, inconspicuous and poorly known cicada in the genus Derotettix, endemic to degraded salt-plain habitats in arid regions of central Argentina, is a relict lineage that is sister to all other modern cicadas. Nuclear and mitochondrial phylogenies of cicadas inferred from probe-based genomic hybrid capture data of both target and non-target loci and a morphological cladogram support this hypothesis. We strengthen this conclusion with genomic data from one of the cicada nutritional bacterial endosymbionts, Sulcia, an ancient and obligate endosymbiont of the larger plant-sucking bugs (Auchenorrhyncha) and an important source of maternally inherited phylogenetic data. We establish Derotettiginae subfam. nov. as a new, monogeneric, fifth cicada subfamily, and compile existing and new data on the distribution, ecology and diet of Derotettix. Our consideration of the palaeoenvironmental literature and host-plant phylogenetics allows us to predict what might have led to the relict status of Derotettix over 100 Myr of habitat change in South America.

Adapterama I: universal stubs and primers for 384 unique dual-indexed or 147,456 combinatorially-indexed Illumina libraries (iTru & iNext)

Massively parallel DNA sequencing offers many benefits, but major inhibitory cost factors include: (1) start-up (i.e., purchasing initial reagents and equipment); (2) buy-in (i.e., getting the smallest possible amount of data from a run); and (3) sample preparation. Reducing sample preparation costs is commonly addressed, but start-up and buy-in costs are rarely addressed. We present dual-indexing systems to address all three of these issues. By breaking the library construction process into universal, re-usable, combinatorial components, we reduce all costs, while increasing the number of samples and the variety of library types that can be combined within runs. We accomplish this by extending the Illumina TruSeq dual-indexing approach to 768 (384 + 384) indexed primers that produce 384 unique dual-indexes or 147,456 (384 × 384) unique combinations. We maintain eight nucleotide indexes, with many that are compatible with Illumina index sequences. We synthesized these indexing primers, purifying them with only standard desalting and placing small aliquots in replicate plates. In qPCR validation tests, 206 of 208 primers tested passed (99% success). We then created hundreds of libraries in various scenarios. Our approach reduces start-up and per-sample costs by requiring only one universal adapter that works with indexed PCR primers to uniquely identify samples. Our approach reduces buy-in costs because: (1) relatively few oligonucleotides are needed to produce a large number of indexed libraries; and (2) the large number of possible primers allows researchers to use unique primer sets for different projects, which facilitates pooling of samples during sequencing. Our libraries make use of standard Illumina sequencing primers and index sequence length and are demultiplexed with standard Illumina software, thereby minimizing customization headaches. In subsequent Adapterama papers, we use these same primers with different adapter stubs to construct amplicon and restriction-site associated DNA libraries, but their use can be expanded to any type of library sequenced on Illumina platforms.