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

Phylogeny and divergence time estimation of Io moths and relatives (Lepidoptera: Saturniidae: Automeris)

The saturniid moth genus Automeris includes 145 described species. Their geographic distribution ranges from the eastern half of North America to as far south as Peru. Automeris moths are cryptically colored, with forewings that resemble dead leaves, and conspicuously colored, elaborate eyespots hidden on their hindwings. Despite their charismatic nature, the evolutionary history and relationships within Automeris and between closely related genera, remain poorly understood. In this study, we present the most comprehensive phylogeny of Automeris to date, including 80 of the 145 described species. We also incorporate two morphologically similar hemileucine genera, Pseudautomeris and Leucanella, as well as a morphologically distinct genus, Molippa. We obtained DNA data from both dry-pinned and ethanol-stored museum specimens and conducted Anchored Hybrid Enrichment (AHE) sequencing to assemble a high-quality dataset for phylogenetic analysis. The resulting phylogeny supports Automeris as a paraphyletic genus, with Leucanella and Pseudautomeris nested within, with the most recent common ancestor dating back to 21 mya. This study lays the foundation for future research on various aspects of Automeris biology, including geographical distribution patterns, potential drivers of speciation, and ecological adaptations such as antipredator defense mechanisms.

To and fro in the archipelago: Repeated inter-island dispersal and New Guinea's orogeny affect diversification of Delias, the world's largest butterfly genus

The world's largest butterfly genus Delias, commonly known as Jezebels, comprises ca. 251 species found throughout Asia, Australia, and Melanesia. Most species are endemic to islands in the Indo-Australian Archipelago or to New Guinea and nearby islands in Melanesia, and many species are restricted to montane habitats over 1200 m. We inferred an extensively sampled and well-supported molecular phylogeny of the group to better understand the spatial and temporal dimensions of its diversification. The remarkable diversity of Delias evolved in just ca. 15-16 Myr (crown age). The most recent common ancestor of a clade with most of the species dispersed out of New Guinea ca. 14 Mya, but at least six subsequently diverging lineages dispersed back to the island. Diversification was associated with frequent dispersal of lineages among the islands of the Indo-Australian Archipelago, and the divergence of sister taxa on a single landmass was rare and occurred only on the largest islands, most notably on New Guinea. We conclude that frequent inter-island dispersal during the Neogene-likely facilitated by frequent sea level change-sparked much diversification during that period. Many extant New Guinea lineages started diversifying 5 Mya, suggesting that orogeny facilitated their diversification. Our results largely agree with the most recently proposed species group classification system, and we use our large taxon sample to extend this system to all described species. Finally, we summarize recent insights to speculate how wing pattern evolution, mimicry, and sexual selection might also contribute to these butterflies' rapid speciation and diversification.

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.

A global phylogeny of butterflies reveals their evolutionary history, ancestral hosts and biogeographic origins

Butterflies are a diverse and charismatic insect group that are thought to have evolved with plants and dispersed throughout the world in response to key geological events. However, these hypotheses have not been extensively tested because a comprehensive phylogenetic framework and datasets for butterfly larval hosts and global distributions are lacking. We sequenced 391 genes from nearly 2,300 butterfly species, sampled from 90 countries and 28 specimen collections, to reconstruct a new phylogenomic tree of butterflies representing 92% of all genera. Our phylogeny has strong support for nearly all nodes and demonstrates that at least 36 butterfly tribes require reclassification. Divergence time analyses imply an origin ~100 million years ago for butterflies and indicate that all but one family were present before the K/Pg extinction event. We aggregated larval host datasets and global distribution records and found that butterflies are likely to have first fed on Fabaceae and originated in what is now the Americas. Soon after the Cretaceous Thermal Maximum, butterflies crossed Beringia and diversified in the Palaeotropics. Our results also reveal that most butterfly species are specialists that feed on only one larval host plant family. However, generalist butterflies that consume two or more plant families usually feed on closely related plants.

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.

Profiling, monitoring and conserving caterpillar fungus in the Himalayan region using anchored hybrid enrichment markers

The collection of caterpillar fungus accounts for 50–70% of the household income of thousands of Himalayan communities and has an estimated market value of $5–11 billion across Asia. However, Himalayan collectors are at multiple economic disadvantages compared with collectors on the Tibetan Plateau because their product is not legally recognized. Using a customized hybrid-enrichment probe set and market-grade caterpillar fungus (with samples up to 30 years old) from 94 production zones across Asia, we uncovered clear geography-based signatures of historical dispersal and significant isolation-by-distance among caterpillar fungus hosts. This high-throughput approach can readily distinguish samples from major production zones with definitive geographical resolution, especially for samples from the Himalayan region that form monophyletic clades in our analysis. Based on these results, we propose a two-step procedure to help local communities authenticate their produce and improve this multi-national trade-route without creating opportunities for illegal exports and other forms of economic exploitation. We argue that policymakers and conservation practitioners must encourage the fair trade of caterpillar fungus in addition to sustainable harvesting to support a trans-boundary conservation effort that is much needed for this natural commodity in the Himalayan region.

Capture by hybridization for full-length barcode-based eukaryotic and prokaryotic biodiversity inventories of deep sea ecosystems

Biodiversity inventory of marine systems remains limited due to unbalanced access to the three ocean dimensions. The use of environmental DNA (eDNA) for metabarcoding allows fast and effective biodiversity inventory and is forecast as a future biodiversity research and biomonitoring tool. However, in poorly understood ecosystems, eDNA results remain difficult to interpret due to large gaps in reference databases and PCR bias limiting the detection of some major phyla. Here, we aimed to circumvent these limitations by avoiding PCR and recollecting larger DNA fragments to improve assignment of detected taxa through phylogenetic reconstruction. We applied capture by hybridization (CBH) to enrich DNA from deep-sea sediment samples and compared the results with those obtained through an up-to-date metabarcoding PCR-based approach (MTB). Originally developed for bacterial communities and targeting 16S rDNA, the CBH approach was applied to 18S rDNA to improve the detection of species forming benthic communities of eukaryotes, with a particular focus on metazoans. The results confirmed the possibility of extending CBH to metazoans with two major advantages: (i) CBH revealed a broader spectrum of prokaryotic, eukaryotic, and particularly metazoan diversity, and (ii) CBH allowed much more robust phylogenetic reconstructions of full-length barcodes with up to 1900 base pairs. This is particularly important for taxa whose assignment is hampered by gaps in reference databases. This study provides a database and probes to apply 18S CBH to diverse marine systems, confirming this promising new tool to improve biodiversity assessments in data-poor ecosystems such as those in the deep sea.

De novo genome assemblies of butterflies

BACKGROUND

The availability of thousands of genomes has enabled new advancements in biology. However, many genomes have not been investigated for their quality. Here we examine quality trends in a taxonomically diverse and well-known group, butterflies (Papilionoidea), and provide draft, de novo assemblies for all available butterfly genomes. Owing to massive genome sequencing investment and taxonomic curation, this is an excellent group to explore genome quality.

FINDINGS

We provide de novo assemblies for all 822 available butterfly genomes and interpret their quality in terms of completeness and continuity. We identify the 50 highest quality genomes across butterflies and conclude that the ringlet, Aphantopus hyperantus, has the highest quality genome. Our post-processing of draft genome assemblies identified 118 butterfly genomes that should not be reused owing to contamination or extremely low quality. However, many draft genomes are of high utility, especially because permissibility of low-quality genomes is dependent on the objective of the study. Our assemblies will serve as a key resource for papilionid genomics, especially for researchers without computational resources.

CONCLUSIONS

Quality metrics and assemblies are typically presented with annotated genome accessions but rarely with de novo genomes. We recommend that studies presenting genome sequences provide the assembly and some metrics of quality because quality will significantly affect downstream results. Transparency in quality metrics is needed to improve the field of genome science and encourage data reuse.

Evolutionary trade-offs between male secondary sexual traits revealed by a phylogeny of the hyperdiverse tribe Eumaeini (Lepidoptera: Lycaenidae)

Male butterflies in the hyperdiverse tribe Eumaeini possess an unusually complex and diverse repertoire of secondary sexual characteristics involved in pheromone production and dissemination. Maintaining multiple sexually selected traits is likely to be metabolically costly, potentially resulting in trade-offs in the evolution of male signals. However, a phylogenetic framework to test hypotheses regarding the evolution and maintenance of male sexual traits in Eumaeini has been lacking. Here, we infer a comprehensive, time-calibrated phylogeny from 379 loci for 187 species representing 91% of the 87 described genera. Eumaeini is a monophyletic group that originated in the late Oligocene and underwent rapid radiation in the Neotropics. We examined specimens of 818 of the 1096 described species (75%) and found that secondary sexual traits are present in males of 91% of the surveyed species. Scent pads and scent patches on the wings and brush organs associated with the genitalia were probably present in the common ancestor of Eumaeini and are widespread throughout the tribe. Brush organs and scent pads are negatively correlated across the phylogeny, exhibiting a trade-off in which lineages with brush organs are unlikely to regain scent pads and vice versa. In contrast, scent patches seem to facilitate the evolution of scent pads, although they are readily lost once scent pads have evolved. Our results illustrate the complex interplay between natural and sexual selection in the origin and maintenance of multiple male secondary sexual characteristics and highlight the potential role of sexual selection spurring diversification in this lineage.

Spatial phylogenetics of butterflies in relation to environmental drivers and angiosperm diversity across North America

Broad-scale, quantitative assessments of insect biodiversity and the factors shaping it remain particularly poorly explored. Here we undertook a spatial phylogenetic analysis of North American butterflies to test whether climate stability and temperature gradients have shaped their diversity and endemism. We also performed the first quantitative comparisons of spatial phylogenetic patterns between butterflies and flowering plants. We expected concordance between the two groups based on shared historical environmental drivers and presumed strong butterfly-host plant specializations. We instead found that biodiversity patterns in butterflies are strikingly different from flowering plants, especially warm deserts. In particular, butterflies show different patterns of phylogenetic clustering compared with flowering plants, suggesting differences in habitat conservation between the two groups. These results suggest that shared biogeographic histories and trophic associations do not necessarily assure similar diversity outcomes. The work has applied value in conservation planning, documenting warm deserts as a North American butterfly biodiversity hotspot.

It's a moth! It's a butterfly! It's the complete mitochondrial genome of the American moth-butterfly Macrosoma conifera (Warren, 1897) (Insecta: Lepidoptera: Hedylidae)!

The taxonomic placement of the moth-butterfly, Macrosoma conifera (Warren 1897) (Lepidoptera: Hedylidae), has been controversial. The 15,344 bp complete M. conifera circular mitogenome, assembled by genome skimming, consists of 81.7% AT nucleotides, 22 tRNAs, 13 protein-coding genes, 2 rRNAs and a control region in the typical butterfly gene order. Macrosoma conifera COX1 features an atypical CGA start codon while ATP6, COX1, COX2, and ND5 exhibit incomplete stop codons completed by the post-transcriptional addition of 3' A residues. Phylogenetic reconstruction places M. conifera as sister to the skippers (Hesperiidae), which is consistent with several recent phylogenetic analyses.

Is Sexual Conflict a Driver of Speciation? A Case Study With a Tribe of Brush-footed Butterflies

Understanding the evolutionary mechanisms governing the uneven distribution of species richness across the tree of life is a great challenge in biology. Scientists have long argued that sexual conflict is a key driver of speciation. This hypothesis, however, has been highly debated in light of empirical evidence. Recent advances in the study of macroevolution make it possible to test this hypothesis with more data and increased accuracy. In the present study, we use phylogenomics combined with four different diversification rate analytical approaches to test whether sexual conflict is a driver of speciation in brush-footed butterflies of the tribe Acraeini. The presence of a sphragis, an external mating plug found in most species among Acraeini, was used as a proxy for sexual conflict. Diversification analyses statistically rejected the hypothesis that sexual conflict is associated with shifts in diversification rates in Acraeini. This result contrasts with earlier studies and suggests that the underlying mechanisms driving diversification are more complex than previously considered. In the case of butterflies, natural history traits acting in concert with abiotic factors possibly play a stronger role in triggering speciation than does sexual conflict. [Acraeini butterflies; arms race; exon capture phylogenomics; Lepidoptera macroevolution; sexual selection; sphragis.].

A Guide to Carrying Out a Phylogenomic Target Sequence Capture Project

High-throughput DNA sequencing techniques enable time- and cost-effective sequencing of large portions of the genome. Instead of sequencing and annotating whole genomes, many phylogenetic studies focus sequencing effort on large sets of pre-selected loci, which further reduces costs and bioinformatic challenges while increasing coverage. One common approach that enriches loci before sequencing is often referred to as target sequence capture. This technique has been shown to be applicable to phylogenetic studies of greatly varying evolutionary depth. Moreover, it has proven to produce powerful, large multi-locus DNA sequence datasets suitable for phylogenetic analyses. However, target capture requires careful considerations, which may greatly affect the success of experiments. Here we provide a simple flowchart for designing phylogenomic target capture experiments. We discuss necessary decisions from the identification of target loci to the final bioinformatic processing of sequence data. We outline challenges and solutions related to the taxonomic scope, sample quality, and available genomic resources of target capture projects. We hope this review will serve as a useful roadmap for designing and carrying out successful phylogenetic target capture studies.

Sequence capture across large phylogenetic scales by using pooled PCR-generated baits: A case study of Lepidoptera

Sequence capture across large phylogenetic scales is not easy because hybridization capture is only effective when the genetic distance between the bait and target is small. Here, we propose a simple but effective strategy to tackle this issue: pooling DNA from a number of selected representative species of different clades to prepare PCR-generated baits to minimize the genetic distance between the bait and target. To demonstrate the utility of this strategy, we newly developed a set of universal nuclear markers (including 94 nuclear protein-coding genes) for Lepidoptera, a superdiverse insect group. We used a DNA pool from six lepidopteran species (representing six superfamilies) to prepare PCR baits for the 94 markers. These homemade PCR baits were used to capture sequence data from 43 species of 17 lepidopteran families, and 94% of the target loci were recovered. We constructed two data sets from the obtained data (one containing ~90 kb target coding sequences and the other containing ~120 kb target + flanking coding sequences). Both data sets yielded highly similar and well-resolved trees with 90% of nodes having >95% bootstrap support. Our capture experiment indicated that using DNA mixtures pooled from different clade-representative species of Lepidoptera to prepare PCR baits can reliably capture a large number of targeted nuclear markers across different Lepidoptera lineages. We hope that this newly developed nuclear marker set will serve as a new phylogenetic tool for Lepidoptera phylogenetics, and the PCR bait preparation strategy can facilitate the application of sequence capture techniques by researchers to accelerate data collection.

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.

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.