Phylogenomics recovers multiple origins of portable case making in caddisflies (Insecta: Trichoptera), nature's underwater architects.

Caddisflies (Trichoptera) are among the most diverse groups of freshwater animals with more than 16 000 described species. They play a fundamental role in freshwater ecology and environmental engineering in streams, rivers and lakes. Because of this, they are frequently used as indicator organisms in biomonitoring programmes. Despite their importance, key questions concerning the evolutionary history of caddisflies, such as the timing and origin of larval case making, remain unanswered owing to the lack of a well-resolved phylogeny. Here, we estimated a phylogenetic tree using a combination of transcriptomes and targeted enrichment data for 207 species, representing 48 of 52 extant families and 174 genera. We calibrated and dated the tree with 33 carefully selected fossils. The first caddisflies originated approximately 295 million years ago in the Permian, and major suborders began to diversify in the Triassic. Furthermore, we show that portable case making evolved in three separate lineages, and shifts in diversification occurred in concert with key evolutionary innovations beyond case making.

Comprehensive phylogeny of Pieridae butterflies reveals strong correlation between diversification and temperature.

Temperature is thought to be a key factor influencing global species richness patterns. We investigate the link between temperature and diversification in the butterfly family Pieridae by combining next generation DNA sequences and published molecular data with fine-grained distribution data. We sampled nearly 600 pierid butterfly species to infer the most comprehensive molecular phylogeny of the family and curated a distribution dataset of more than 800,000 occurrences. We found strong evidence that species in environments with more stable daily temperatures or cooler maximum temperatures in the warm seasons have higher speciation rates. Furthermore, speciation and extinction rates decreased in tandem with global temperatures through geological time, resulting in a constant net diversification.

Large-scale genomic data reveal the phylogeny and evolution of owlet moths (Noctuoidea).

The owlet moths (Noctuoidea; ~43-45K described species) are one of the most ecologically diverse and speciose superfamilies of animals. Moreover, they comprise some of the world's most notorious pests of agriculture and forestry. Despite their contributions to terrestrial biodiversity and impacts on ecosystems and economies, the evolutionary history of Noctuoidea remains unclear because the superfamily lacks a statistically robust phylogenetic and temporal framework. We reconstructed the phylogeny of Noctuoidea using data from 1234 genes (946.4 kb nucleotides) obtained from the genome and transcriptome sequences of 76 species. The relationships among the six families of Noctuoidea were well resolved and consistently recovered based on both concatenation and gene coalescence approaches, supporting the following relationships: Oenosandridae + (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))). A Yule tree prior with three unlinked molecular clocks was identified as the preferred BEAST analysis using marginal-likelihood estimations. The crown age of Noctuoidea was estimated at 74.5 Ma, with most families originating before the end of the Paleogene (23 Ma). Our study provides the first statistically robust phylogenetic and temporal framework for Noctuoidea, including all families of owlet moths, based on large-scale genomic data.

HerediGene Population Study IT infrastructure: A model to support genomic research recruitment and precision public health.

The HerediGene Population Study is a large research study focused on identifying new genetic biomarkers for disease prevention, diagnosis, prognosis, and development of new therapeutics. A substantial IT infrastructure evolved to reach enrollment targets and return results to participants. More than 170,000 participants have been enrolled in the study to date, with 5.87% of those whole genome sequenced and 0.46% of those genotyped harboring pathogenic variants. Among other purposes, this infrastructure supports: (1) identifying candidates from clinical criteria, (2) monitoring for qualifying clinical events (e.g., blood draw), (3) contacting candidates, (4) obtaining consent electronically, (5) initiating lab orders, (6) integrating consent and lab orders into clinical workflow, (7) de-identifying samples and clinical data, (8) shipping/transmitting samples and clinical data, (9) genotyping/sequencing samples, (10) and re-identifying and returning results for participants where applicable. This study may serve as a model for similar genomic research and precision public health initiatives.

The Development of an Infrastructure to Facilitate the Use of Whole Genome Sequencing for Population Health.

The clinical use of genomic analysis has expanded rapidly resulting in an increased availability and utility of genomic information in clinical care. We have developed an infrastructure utilizing informatics tools and clinical processes to facilitate the use of whole genome sequencing data for population health management across the healthcare system. Our resulting framework scaled well to multiple clinical domains in both pediatric and adult care, although there were domain specific challenges that arose. Our infrastructure was complementary to existing clinical processes and well-received by care providers and patients. Informatics solutions were critical to the successful deployment and scaling of this program. Implementation of genomics at the scale of population health utilizes complicated technologies and processes that for many health systems are not supported by current information systems or in existing clinical workflows. To scale such a system requires a substantial clinical framework backed by informatics tools to facilitate the flow and management of data. Our work represents an early model that has been successful in scaling to 29 different genes with associated genetic conditions in four clinical domains. Work is ongoing to optimize informatics tools; and to identify best practices for translation to smaller healthcare systems.

Genomic selection reveals hidden relatedness and increased breeding efficiency in western redcedar polycross breeding.

Western redcedar (WRC) is an ecologically and economically important forest tree species characterized by low genetic diversity with high self-compatibility and high heartwood durability. Using sequence capture genotyping of target genic and non-genic regions, we genotyped 44 parent trees and 1520 offspring trees representing 26 polycross (PX) families collected from three progeny test sites using 45,378 SNPs. Trees were phenotyped for eight traits related to growth, heartwood and foliar chemistry associated with wood durability and deer browse resistance. We used the genomic realized relationship matrix for paternity assignment, maternal pedigree correction, and to estimate genetic parameters. We compared genomics-based (GBLUP) and two pedigree-based (ABLUP: polycross and reconstructed full-sib [FS] pedigrees) models. Models were extended to estimate dominance genetic effects. Pedigree reconstruction revealed significant unequal male contribution and separated the 26 PX families into 438 FS families. Traditional maternal PX pedigree analysis resulted in up to 51% overestimation in genetic gain and 44% in diversity. Genomic analysis resulted in up to 22% improvement in offspring breeding value (BV) theoretical accuracy, 35% increase in expected genetic gain for forward selection, and doubled selection intensity for backward selection. Overall, all traits showed low to moderate heritability (0.09-0.28), moderate genotype by environment interaction (type-B genetic correlation: 0.51-0.80), low to high expected genetic gain (6.01%-55%), and no significant negative genetic correlation reflecting no large trade-offs for multi-trait selection. Only three traits showed a significant dominance effect. GBLUP resulted in smaller but more accurate heritability estimates for five traits, but larger estimates for the wood traits. Comparison between all, genic-coding, genic-non-coding and intergenic SNPs showed little difference in genetic estimates. In summary, we show that GBLUP overcomes the PX limitations, successfully captures expected historical and hidden relatedness as well as linkage disequilibrium (LD), and results in increased breeding efficiency in WRC.

The western redcedar genome reveals low genetic diversity in a self-compatible conifer.

We assembled the 9.8-Gbp genome of western redcedar (WRC; Thuja plicata), an ecologically and economically important conifer species of the Cupressaceae. The genome assembly, derived from a uniquely inbred tree produced through five generations of self-fertilization (selfing), was determined to be 86% complete by BUSCO analysis, one of the most complete genome assemblies for a conifer. Population genomic analysis revealed WRC to be one of the most genetically depauperate wild plant species, with an effective population size of approximately 300 and no significant genetic differentiation across its geographic range. Nucleotide diversity, π, is low for a continuous tree species, with many loci showing zero diversity, and the ratio of π at zero- to fourfold degenerate sites is relatively high (approximately 0.33), suggestive of weak purifying selection. Using an array of genetic lines derived from up to five generations of selfing, we explored the relationship between genetic diversity and mating system. Although overall heterozygosity was found to decline faster than expected during selfing, heterozygosity persisted at many loci, and nearly 100 loci were found to deviate from expectations of genetic drift, suggestive of associative overdominance. Nonreference alleles at such loci often harbor deleterious mutations and are rare in natural populations, implying that balanced polymorphisms are maintained by linkage to dominant beneficial alleles. This may account for how WRC remains responsive to natural and artificial selection, despite low genetic diversity.

Anti-bat ultrasound production in moths is globally and phylogenetically widespread.

Warning signals are well known in the visual system, but rare in other modalities. Some moths produce ultrasonic sounds to warn bats of noxious taste or to mimic unpalatable models. Here, we report results from a long-term study across the globe, assaying moth response to playback of bat echolocation. We tested 252 genera, spanning most families of large-bodied moths, and document anti-bat ultrasound production in 52 genera, with eight subfamily origins described. Based on acoustic analysis of ultrasonic emissions and palatability experiments with bats, it seems that acoustic warning and mimicry are the raison d'être for sound production in most moths. However, some moths use high-duty-cycle ultrasound capable of jamming bat sonar. In fact, we find preliminary evidence of independent origins of sonar jamming in at least six subfamilies. Palatability data indicate that jamming and warning are not mutually exclusive strategies. To explore the possible organization of anti-bat warning sounds into acoustic mimicry rings, we intensively studied a community of moths in Ecuador and, using machine-learning approaches, found five distinct acoustic clusters. While these data represent an early understanding of acoustic aposematism and mimicry across this megadiverse insect order, it is likely that ultrasonically signaling moths comprise one of the largest mimicry complexes on earth.

Hidden Phylogenomic Signal Helps Elucidate Arsenurine Silkmoth Phylogeny and the Evolution of Body Size and Wing Shape Trade-Offs.

One of the key objectives in biological research is understanding how evolutionary processes have produced Earth's diversity. A critical step toward revealing these processes is an investigation of evolutionary tradeoffs-that is, the opposing pressures of multiple selective forces. For millennia, nocturnal moths have had to balance successful flight, as they search for mates or host plants, with evading bat predators. However, the potential for evolutionary trade-offs between wing shape and body size are poorly understood. In this study, we used phylogenomics and geometric morphometrics to examine the evolution of wing shape in the wild silkmoth subfamily Arsenurinae (Saturniidae) and evaluate potential evolutionary relationships between body size and wing shape. The phylogeny was inferred based on 782 loci from target capture data of 42 arsenurine species representing all 10 recognized genera. After detecting in our data one of the most vexing problems in phylogenetic inference-a region of a tree that possesses short branches and no "support" for relationships (i.e., a polytomy), we looked for hidden phylogenomic signal (i.e., inspecting differing phylogenetic inferences, alternative support values, quartets, and phylogenetic networks) to better illuminate the most probable generic relationships within the subfamily. We found there are putative evolutionary trade-offs between wing shape, body size, and the interaction of fore- and hindwing (HW) shape. Namely, body size tends to decrease with increasing HW length but increases as forewing (FW) shape becomes more complex. Additionally, the type of HW (i.e., tail or no tail) a lineage possesses has a significant effect on the complexity of FW shape. We outline possible selective forces driving the complex HW shapes that make Arsenurinae, and silkmoths as a whole, so charismatic. [Anchored hybrid enrichment; Arsenurinae; geometric morphometrics; Lepidoptera; phylogenomics; Saturniidae.].

A target enrichment probe set for resolving the flagellate land plant tree of life.

PREMISE: New sequencing technologies facilitate the generation of large-scale molecular data sets for constructing the plant tree of life. We describe a new probe set for target enrichment sequencing to generate nuclear sequence data to build phylogenetic trees with any flagellate land plants, including hornworts, liverworts, mosses, lycophytes, ferns, and all gymnosperms. METHODS: We leveraged existing transcriptome and genome sequence data to design the GoFlag 451 probes, a set of 56,989 probes for target enrichment sequencing of 451 exons that are found in 248 single-copy or low-copy nuclear genes across flagellate plant lineages. RESULTS: Our results indicate that target enrichment using the GoFlag451 probe set can provide large nuclear data sets that can be used to resolve relationships among both distantly and closely related taxa across the flagellate land plants. We also describe the GoFlag 408 probes, an optimized probe set covering 408 of the 451 exons from the GoFlag 451 probe set that is commercialized by RAPiD Genomics. CONCLUSIONS: A target enrichment approach using the new probe set provides a relatively low-cost solution to obtain large-scale nuclear sequence data for inferring phylogenetic relationships across flagellate land plants.

Phylogeny and classification of Odonata using targeted genomics.

Dragonflies and damselflies are a charismatic, medium-sized insect order (~6300 species) with a unique potential to approach comparative research questions. Their taxonomy and many ecological traits for a large fraction of extant species are relatively well understood. However, until now, the lack of a large-scale phylogeny based on high throughput data with the potential to connect both perspectives has precluded comparative evolutionary questions for these insects. Here, we provide an ordinal hypothesis of classification based on anchored hybrid enrichment using a total of 136 species representing 46 of the 48 families or incertae sedis, and a total of 478 target loci. Our analyses recovered the monophyly for all three suborders: Anisoptera, Anisozygoptera and Zygoptera. Although the backbone of the topology was reinforced and showed the highest support values to date, our genomic data was unable to stronglyresolve portions of the topology. In addition, a quartet sampling approach highlights the potential evolutionary scenarios that may have shaped evolutionary phylogeny (e.g., incomplete lineage sorting and introgression) of this taxon. Finally, in light of our phylogenomic reconstruction and previous morphological and molecular information we proposed an updated odonate classification and define five new families (Amanipodagrionidae fam. nov., Mesagrionidae fam. nov., Mesopodagrionidae fam. nov., Priscagrionidae fam. nov., Protolestidae fam. nov.) and reinstate another two (Rhipidolestidae stat. res., Tatocnemididae stat. res.). Additionally, we feature the problematic taxonomic groupings for examination in future studies to improve our current phylogenetic hypothesis.

Phylogenomics reveals the evolutionary timing and pattern of butterflies and moths.

Butterflies and moths (Lepidoptera) are one of the major superradiations of insects, comprising nearly 160,000 described extant species. As herbivores, pollinators, and prey, Lepidoptera play a fundamental role in almost every terrestrial ecosystem. Lepidoptera are also indicators of environmental change and serve as models for research on mimicry and genetics. They have been central to the development of coevolutionary hypotheses, such as butterflies with flowering plants and moths' evolutionary arms race with echolocating bats. However, these hypotheses have not been rigorously tested, because a robust lepidopteran phylogeny and timing of evolutionary novelties are lacking. To address these issues, we inferred a comprehensive phylogeny of Lepidoptera, using the largest dataset assembled for the order (2,098 orthologous protein-coding genes from transcriptomes of 186 species, representing nearly all superfamilies), and dated it with carefully evaluated synapomorphy-based fossils. The oldest members of the Lepidoptera crown group appeared in the Late Carboniferous (∼300 Ma) and fed on nonvascular land plants. Lepidoptera evolved the tube-like proboscis in the Middle Triassic (∼241 Ma), which allowed them to acquire nectar from flowering plants. This morphological innovation, along with other traits, likely promoted the extraordinary diversification of superfamily-level lepidopteran crown groups. The ancestor of butterflies was likely nocturnal, and our results indicate that butterflies became day-flying in the Late Cretaceous (∼98 Ma). Moth hearing organs arose multiple times before the evolutionary arms race between moths and bats, perhaps initially detecting a wide range of sound frequencies before being co-opted to specifically detect bat sonar. Our study provides an essential framework for future comparative studies on butterfly and moth evolution.

A phylogenomic framework, evolutionary timeline and genomic resources for comparative studies of decapod crustaceans.

Comprising over 15 000 living species, decapods (crabs, shrimp and lobsters) are the most instantly recognizable crustaceans, representing a considerable global food source. Although decapod systematics have received much study, limitations of morphological and Sanger sequence data have yet to produce a consensus for higher-level relationships. Here, we introduce a new anchored hybrid enrichment kit for decapod phylogenetics designed from genomic and transcriptomic sequences that we used to capture new high-throughput sequence data from 94 species, including 58 of 179 extant decapod families, and 11 of 12 major lineages. The enrichment kit yields 410 loci (greater than 86 000 bp) conserved across all lineages of Decapoda, more clade-specific molecular data than any prior study. Phylogenomic analyses recover a robust decapod tree of life strongly supporting the monophyly of all infraorders, and monophyly of each of the reptant, 'lobster' and 'crab' groups, with some results supporting pleocyemate monophyly. We show that crown decapods diverged in the Late Ordovician and most crown lineages diverged in the Triassic-Jurassic, highlighting a cryptic Palaeozoic history, and post-extinction diversification. New insights into decapod relationships provide a phylogenomic window into morphology and behaviour, and a basis to rapidly and cheaply expand sampling in this economically and ecologically significant invertebrate clade.

Unioverse: A phylogenomic resource for reconstructing the evolution of freshwater mussels (Bivalvia, Unionoida).

Freshwater mussels (order Unionoida) are a diverse radiation of parasitic bivalves that require temporary larval encystment on vertebrate hosts to complete metamorphosis to free-living juveniles. The freshwater mussel-fish symbiosis represents a useful relationship for understanding eco-evolutionary dynamics in freshwater ecosystems but the practicality of this promising model system is undermined by the absence of a stable freshwater mussel phylogeny. Inadequate character sampling is the primary analytical impediment obfuscating a coherent phylogeny of freshwater mussels, specifically the lack of nuclear molecular markers appropriate for reconstructing supraspecific relationships and testing macroevolutionary hypotheses. The objective of this study is to develop a phylogenomic resource, specifically an anchored hybrid enrichment probe set, capable of capturing hundreds of molecular markers from taxa distributed across the entirety of freshwater mussel biodiversity. Our freshwater mussel specific anchored hybrid enrichment probe set, called Unioverse, successfully captures hundreds of nuclear protein-coding loci from all major lineages of the Unionoida and will facilitate more data-rich and taxonomically inclusive reconstructions of freshwater mussel evolution. We demonstrate the utility of this resource at three disparate evolutionary scales by estimating a backbone phylogeny of the Bivalvia with a focus on the Unionoida, reconstructing the subfamily-level relationships of the Unionidae, and recovering the systematic position of the phylogenetically unstable genus Plectomerus.

Anchored hybrid enrichment phylogenomics resolves the backbone of erebine moths.

The subfamily Erebinae (Lepidoptera, Erebidae) includes approximately 10,000 species with many still undescribed. It is one of the most diverse clades within the moth superfamily Noctuoidea and encompasses a diversity of ecological habits. Erebine caterpillars feed on a broad range of host plants including several economically important crops. Adults possess a unique array of adaptations for predator defense, including some of the most sensitive hearing organs (tympana) across the Lepidoptera and striking wing coloration to startle visual predators. Despite the relevance of these moths to agriculture and ecological research, a robust phylogenetic framework is lacking. Here we used anchored hybrid enrichment, a relatively new approach in phylogenomics, to resolve relationships among the subfamily. Using the recently developed Lep1 anchored hybrid enrichment probe set, 658 gene fragments with an average length of 320 bp were captured from an exemplar set of 75 erebine species, representing 73 genera and 23 tribes. While the total number of erebine tribes is not firmly established, this represents at least 75% of known tribal level diversity. Anchored hybrid enrichment data were partitioned by locus and by codon position for maximum likelihood phylogenetic analysis and coalescent-based species-tree approaches. Results from our study provided strong nodal support (BP ≥ 95) for nearly all nodes in the partitioned ML tree, solidifying many relationships that were previously uncertain or moderately supported based on morphology or a smaller number of gene fragments. Likelihood analyses confidently resolved the placement of Acantholipini as a sister tribe to Sypnini and all other Erebinae. The remaining tribes were placed in a single, strongly supported clade split into two major subclades. Additionally, 25 tropical species that did not have previous tribal assignments are confidently placed on the phylogeny. Statistical comparisons with Shimodaira-Hasegawa (SH) tests found that our maximum likelihood trees were significantly more likely than alternative hypotheses. This study demonstrates the utility of anchored phylogenomics for resolving relationships within subfamilies of Lepidoptera.

Four hundred shades of brown: Higher level phylogeny of the problematic Euptychiina (Lepidoptera, Nymphalidae, Satyrinae) based on hybrid enrichment data.

Relationships within satyrine butterflies have been notoriously difficult to resolve using both morphology and Sanger sequencing methods, and this is particularly true for the mainly Neotropical subtribe Euptychiina, which contains about 400 described species. Known larvae of Euptychiina feed on grasses and sedges, with the exception of the genus Euptychia, which feed on mosses and lycopsids, and the butterflies occur widely in rainforest, cloudforest and grassland habitats, where they are often abundant. Several previous molecular and morphological studies have made significant progress in tackling the systematics of the group, but many relationships remain unresolved, with long-branch-attraction artifacts being a major problem. Additionally, the monophyly of the clade remains uncertain, with Euptychia possibly not being closely related to the remainder of the clade. Here we present a backbone phylogeny of the subtribe based on 106 taxa, 368 nuclear loci, and over 180,000 bps obtained through hybrid enrichment. Using both concatenation and species tree approaches (IQ-TREE, EXABAYES, ASTRAL), we can for the first time strongly confirm the monophyly of Euptychiina with Euptychia being the sister group to the remainder of the clade. The Euptychiina is divided into nine well supported clades, but the placement of a few genera such as Hermeuptychia, Pindis and the Chloreuptychia catharina group still remain uncertain. As partially indicated in previous studies, the genera Cissia, Chloreuptychia, Magneuptychia, Megisto, Splendeuptychia and Euptychoides, among others, were found to be highly polyphyletic and revisions are in preparation. The phylogeny will provide a strong backbone for the analysis of datasets in development that are much more taxonomically comprehensive but have orders of magnitude fewer loci. This study therefore represents a critical step towards resolving the higher classification and studying the evolution of this highly diverse lineage.

Origin and macroevolution of micro-moths on sunken Hawaiian Islands.

The origins and evolution of Hawaiian biodiversity are a matter of controversy, and the mechanisms of lineage diversification for many organisms on this remote archipelago remain unclear. Here we focus on the poorly known endemic leaf-mining moth genus Philodoria (Lepidoptera, Gracillariidae), whose species feed on a diversity of Hawaiian plant lineages, many of which are critically endangered. We use anchored hybrid enrichment to assemble the first phylogenomic dataset (507 loci) for any Hawaiian animal taxon. To uncover the timing and pattern of diversification of these moths, we apply two frequently used dating calibration strategies, biogeographic calibrations and secondary calibrations. Island calibrations on their own resulted in much younger and unrealistic dates compared to strategies that relied on secondary calibrations. Philodoria probably originated on the now partially sunken islands of Laysan or Lisianski, approximately 21 Ma, and were associated with host plants in the families Ebenaceae, Malvaceae or Primulaceae. Major feeding groups associated with specific host-plant families originated soon after the plants colonized the islands. Allopatric isolation and host shifts, in concert and independently, probably play major roles in the diversification of Philodoria Our dating results indicate that Philodoria is among the oldest known Hawaiian arthropod lineages, and that island calibrations alone can lead to unrealistically young dates.

Anchored phylogenomics illuminates the skipper butterfly tree of life.

BACKGROUND: Butterflies (Papilionoidea) are perhaps the most charismatic insect lineage, yet phylogenetic relationships among them remain incompletely studied and controversial. This is especially true for skippers (Hesperiidae), one of the most species-rich and poorly studied butterfly families. METHODS: To infer a robust phylogenomic hypothesis for Hesperiidae, we sequenced nearly 400 loci using Anchored Hybrid Enrichment and sampled all tribes and more than 120 genera of skippers. Molecular datasets were analyzed using maximum-likelihood, parsimony and coalescent multi-species phylogenetic methods. RESULTS: All analyses converged on a novel, robust phylogenetic hypothesis for skippers. Different optimality criteria and methodologies recovered almost identical phylogenetic trees with strong nodal support at nearly all nodes and all taxonomic levels. Our results support Coeliadinae as the sister group to the remaining skippers, the monotypic Euschemoninae as the sister group to all other subfamilies but Coeliadinae, and the monophyly of Eudaminae plus Pyrginae. Within Pyrginae, Celaenorrhinini and Tagiadini are sister groups, the Neotropical firetips, Pyrrhopygini, are sister to all other tribes but Celaenorrhinini and Tagiadini. Achlyodini is recovered as the sister group to Carcharodini, and Erynnini as sister group to Pyrgini. Within the grass skippers (Hesperiinae), there is strong support for the monophyly of Aeromachini plus remaining Hesperiinae. The giant skippers (Agathymus and Megathymus) once classified as a subfamily, are recovered as monophyletic with strong support, but are deeply nested within Hesperiinae. CONCLUSIONS: Anchored Hybrid Enrichment sequencing resulted in a large amount of data that built the foundation for a new, robust evolutionary tree of skippers. The newly inferred phylogenetic tree resolves long-standing systematic issues and changes our understanding of the skipper tree of life. These resultsenhance understanding of the evolution of one of the most species-rich butterfly families.

Phylogenetics of moth-like butterflies (Papilionoidea: Hedylidae) based on a new 13-locus target capture probe set.

The Neotropical moth-like butterflies (Hedylidae) are perhaps the most unusual butterfly family. In addition to being species-poor, this family is predominantly nocturnal and has anti-bat ultrasound hearing organs. Evolutionary relationships among the 36 described species are largely unexplored. A new, target capture, anchored hybrid enrichment probe set ('BUTTERFLY2.0') was developed to infer relationships of hedylids and some of their butterfly relatives. The probe set includes 13 genes that have historically been used in butterfly phylogenetics. Our dataset comprised of up to 10,898 aligned base pairs from 22 hedylid species and 19 outgroups. Eleven of the thirteen loci were successfully captured from all samples, and the remaining loci were captured from ≥94% of samples. The inferred phylogeny was consistent with recent molecular studies by placing Hedylidae sister to Hesperiidae, and the tree had robust support for 80% of nodes. Our results are also consistent with morphological studies, with Macrosoma tipulata as the sister species to all remaining hedylids, followed by M. semiermis sister to the remaining species in the genus. We tested the hypothesis that nocturnality evolved once from diurnality in Hedylidae, and demonstrate that the ancestral condition was likely diurnal, with a shift to nocturnality early in the diversification of this family. The BUTTERFLY2.0 probe set includes standard butterfly phylogenetics markers, captures sequences from decades-old museum specimens, and is a cost-effective technique to infer phylogenetic relationships of the butterfly tree of life.

A Comprehensive and Dated Phylogenomic Analysis of Butterflies.

Butterflies (Papilionoidea), with over 18,000 described species [1], have captivated naturalists and scientists for centuries. They play a central role in the study of speciation, community ecology, biogeography, climate change, and plant-insect interactions and include many model organisms and pest species [2, 3]. However, a robust higher-level phylogenetic framework is lacking. To fill this gap, we inferred a dated phylogeny by analyzing the first phylogenomic dataset, including 352 loci (> 150,000 bp) from 207 species representing 98% of tribes, a 35-fold increase in gene sampling and 3-fold increase in taxon sampling over previous studies [4]. Most data were generated with a new anchored hybrid enrichment (AHE) [5] gene kit (BUTTERFLY1.0) that includes both new and frequently used (e.g., [6]) informative loci, enabling direct comparison and future dataset merging with previous studies. Butterflies originated around 119 million years ago (mya) in the late Cretaceous, but most extant lineages diverged after the Cretaceous-Paleogene (K-Pg) mass-extinction 65 mya. Our analyses support swallowtails (Papilionidae) as sister to all other butterflies, followed by skippers (Hesperiidae) + the nocturnal butterflies (Hedylidae) as sister to the remainder, indicating a secondary reversal from diurnality to nocturnality. The whites (Pieridae) were strongly supported as sister to brush-footed butterflies (Nymphalidae) and blues + metalmarks (Lycaenidae and Riodinidae). Ant association independently evolved once in Lycaenidae and twice in Riodinidae. This study overturns prior notions of the taxon's evolutionary history, as many long-recognized subfamilies and tribes are para- or polyphyletic. It also provides a much-needed backbone for a revised classification of butterflies and for future comparative studies including genome evolution and ecology.