Phylogenomics reveals the evolutionary timing and pattern of butterflies and moths

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

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

Global biogeographic regions for ants have complex relationships with those for plants and tetrapods

On a global scale, biodiversity is geographically structured into regions of biotic similarity. Delineating these regions has been mostly targeted for tetrapods and plants, but those for hyperdiverse groups such as insects are relatively unknown. Insects may have higher biogeographic congruence with plants than tetrapods due to their tight ecological and evolutionary links with the former, but it remains untested. Here, we develop a global regionalization for a major and widespread insect group, ants, based on the most comprehensive distributional and phylogenetic information to date, and examine its similarity to regionalizations for tetrapods and vascular plants. Our ant regionalization supports the newly proposed Madagascan and Sino-Japanese realms based on tetrapod delineations, and it recovers clusters observed in plants but not in tetrapods, such as the Holarctic and Indo-Pacific realms. Quantitative comparison suggests strong associations among different groups-plants showed a higher congruence with ants than with tetrapods. These results underscore the wide congruence of diverse distribution patterns across the tree of life and the similarities shared by insects and plants that are not captured by tetrapod groups. Our analysis highlights the importance of developing global biogeographic maps for insect groups to obtain a more comprehensive geographic picture of life on Earth.

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.

Molecular phylogeny of Hesperiidae (Lepidoptera) with an emphasis on Asian and African genera

Despite considerable research efforts in recent years, the deeper phylogenetic relationships among skipper butterflies (Hesperiidae) remain unresolved. This is primarily because of limited sampling, especially within Asian and African lineages. In this study, we consolidated previous data and extensively sampled Asian and African taxa to elucidate the phylogenetic relationships within Hesperiidae. The molecular dataset comprised sequences from two mitochondrial and two nuclear gene regions from 563 species that represented 353 genera. Our analyses revealed seven subfamilies within Hesperiidae: Coeliadinae, Euschemoninae, Eudaminae, Pyrginae, Heteropterinae, Trapezitinae, and Hesperiinae. The systematics of most tribes and genera aligned with those of prior studies. However, notable differences were observed in several tribes and genera. Overall, the position of taxa assigned to incertae sedis in Hesperiinae is largely clarified in this study. Our results strongly support the monophyly of the tribe Tagiadini (Pyrginae), and the systematics of some genera are clarified with comprehensive discussion. We recognize 15 tribes within the subfamily Hesperiinae. Of these, nine tribes are discussed in detail: Aeromachini, Astictopterini, Erionotini, Unkanini (new status), Ancistroidini, Ismini (confirmed status), Plastingini (new status), Gretnini (confirmed status), and Eetionini (confirmed status). We propose four subtribes within Astictopterini: Hypoleucina subtrib.n., Aclerosina, Cupithina, and Astictopterina. Furthermore, we describe a new genus (Hyarotoidesgen.n.) and reinstate two genera (Zeareinst.stat. and Separeinst.stat.) as valid. Additionally, we propose several new combinations: Zea mythecacomb.n.,Sepa bononiacomb.n. & reinst.stat., and Sepa umbrosacomb.n. This study, with extensive sampling of Asian and African taxa, greatly enhances the understanding of the knowledge of the skipper tree of life.

Prey can detect predators via electroreception in air

Predators and prey benefit from detecting sensory cues of each other's presence. As they move through their environment, terrestrial animals accumulate electrostatic charge. Because electric charges exert forces at a distance, a prey animal could conceivably sense electrical forces to detect an approaching predator. Here, we report such a case of a terrestrial animal detecting its predators by electroreception. We show that predatory wasps are charged, thus emit electric fields, and that caterpillars respond to such fields with defensive behaviors. Furthermore, the mechanosensory setae of caterpillars are deflected by these electrostatic forces and are tuned to the wingbeat frequency of their insect predators. This ability unveils a dimension of the sensory interactions between prey and predators and is likely widespread among terrestrial animals.

Middle Jurassic insect mines on gymnosperms provide missing links to early mining evolution

We investigated the mining mode of insect feeding, involving larval consumption of a plant's internal tissues, from the Middle Jurassic (165 million years ago) Daohugou locality of Northeastern China. Documentation of mining from the Jurassic Period is virtually unknown, and results from this time interval would address mining evolution during the temporal gap of mine-seed plant diversifications from the previous Late Triassic to the subsequent Early Cretaceous. Plant fossils were examined with standard microscopic procedures for herbivory and used the standard functional feeding group-damage-type system of categorizing damage. All fossil mines were photographed and databased. We examined 2014 plant specimens, of which 27 occurrences on 14 specimens resulted in eight, new, mine damage types (DTs) present on six genera of bennettitalean, ginkgoalean, and pinalean gymnosperms. Three conclusions emerge from this study. First, these mid-Mesozoic mines are morphologically conservative and track plant host anatomical structure rather than plant phylogeny. Second, likely insect fabricators of these mines were three basal lineages of polyphagan beetles, four basal lineages of monotrysian moths, and a basal lineage tenthredinoid sawflies. Third, the nutrition hypothesis, indicating that miners had greater access to nutritious, inner tissues of new plant lineages, best explains mine evolution during the mid-Mesozoic.

Duplications and Losses of the Detoxification Enzyme Glycosyltransferase 1 Are Related to Insect Adaptations to Plant Feeding

Insects have developed sophisticated detoxification systems to protect them from plant secondary metabolites while feeding on plants to obtain necessary nutrients. As an important enzyme in the system, glycosyltransferase 1 (GT1) conjugates toxic compounds to mitigate their harm to insects. However, the evolutionary link between GT1s and insect plant feeding remains elusive. In this study, we explored the evolution of GT1s across different insect orders and feeding niches using publicly available insect genomes. GT1 is widely present in insect species; however, its gene number differs among insect orders. Notably, plant-sap-feeding species have the highest GT1 gene numbers, whereas blood-feeding species display the lowest. GT1s appear to be associated with insect adaptations to different plant substrates in different orders, while the shift to non-plant feeding is related to several losses of GT1s. Most large gene numbers are likely the consequence of tandem duplications showing variations in collinearity among insect orders. These results reveal the potential relationships between the evolution of GT1s and insect adaptation to plant feeding, facilitating our understanding of the molecular mechanisms underlying insect-plant interactions.

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.

Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

Systematics and evolutionary dynamics of insect-fern interactions in the specialized fern-spore feeding Cuprininae (Lepidoptera, Stathmopodidae)

Fern-spore-feeding (FSF) is rare and found in only four families of Lepidoptera. Stathmopodidae is the most speciose family that contains FSF species, and its subfamily Cuprininae exclusively specializes on FSF. However, three species of Stathmopodinae also specialize on FSF. To better understand the evolutionary history of FSF and, more generally, the significance of specialization on a peculiar host, a phylogenetic and taxonomic revision for this group is necessary. We reconstructed the most comprehensive molecular phylogeny, including one mitochondrial and four nuclear genes, of Stathmopodidae to date, including 137 samples representing 62 species, with a particular focus on the FSF subfamily, Cuprininae, including 33 species (41% of named species) from 6 of the 7 Cuprininae genera. Species from two other subfamilies, Stathmopodinae and Atkinsoniinae, were also included. We found that FSF evolved only once in Stathmopodidae and that the previous hypothesis of multiple origins of FSF was misled by inadequate taxonomy. Moreover, we showed that (1) speciation/extinction rates do not differ significantly between FSF and non-FSF groups and that (2) oligophage is the ancestral character state in Cuprininae. We further revealed that a faster rate of accumulating specialists over time, and thus a higher number of specialists, was achieved by a higher transition rate from oligophagages to specialists compared to the transition rate in the opposite direction. We finish by describing three new genera, Trigonodagen. nov., Petalagen. nov., and Pediformisgen. nov., and revalidating five genera: Cuprina, Calicotis, Thylacosceles, Actinoscelis, Thylacosceloides in Cuprininae, and we provide an updated taxonomic key to genera and a revised global checklist of Cuprininae.

Mitochondrial genomes of four slug moths (Lepidoptera, Limacodidae): Genome description and phylogenetic implications

The family Limacodidae belongs to the superfamily Zygaenoidea, which includes 1672 species commonly referred to as slug moths. Limacodidae larvae are major pests for many economically important plant species and can cause human dermatitis. At present, the structure of the mitochondrial genome (mitogenome), phylogenetic position, and adaptive evolution of slug moths are poorly understood. Herein, the mitogenomes of Parasa lepida, Phlossa conjuncta, Thosea sinensis, and Setora sinensis were sequenced and compared with other available mitogenome sequences to better characterize the mitogenomic diversity and evolution of this moth family. The mitogenomes of P. lepida, P. conjuncta, T. sinensis, and S. sinensis were confirmed to be circular in structure with lengths of 15,575 bp, 15,553 bp, 15,535 bp, and 15,529 bp, respectively. The Limacodidae mitogenomes exhibited similar nucleotide composition, codon usage, RNA structure, and control region patterns, indicating the conservation of the mitogenome in the family Limacodidae. A sliding window, Ka/Ks, and genetic distance analyses revealed that the atp8 and nad6 genes exhibited the highest levels of variability and the most rapid evolutionary rates among the 13 protein-coding genes (PCGs) encoded in these Limacodidae mitogenomes, suggesting that they may offer value as candidate DNA markers. The phylogenetic analysis recovered the overall relationship as Tortricoidea + (Sesiidae + (Zygaenoidea + (Cossoidea/+Choreutoidea + (others)))). Within Zygaenoidea, Limacodidae was recovered as monophyletic, and the phylogenetic relationships were recovered as (Phaudidae + Zyganidae) + Limacodidae in all six phylogenetic trees. The analysis indicated that P. lepida, P. conjuncta, T. sinensis, and S. sinensis are members of the Limacodidae.

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 macroevolutionary role for chromosomal fusion and fission in Erebia butterflies

The impact of large-scale chromosomal rearrangements, such as fusions and fissions, on speciation is a long-standing conundrum. We assessed whether bursts of change in chromosome numbers resulting from chromosomal fusion or fission are related to increased speciation rates in Erebia, one of the most species-rich and karyotypically variable butterfly groups. We established a genome-based phylogeny and used state-dependent birth-death models to infer trajectories of karyotype evolution. We demonstrated that rates of anagenetic chromosomal changes (i.e., along phylogenetic branches) exceed cladogenetic changes (i.e., at speciation events), but, when cladogenetic changes occur, they are mostly associated with chromosomal fissions rather than fusions. We found that the relative importance of fusion and fission differs among Erebia clades of different ages and that especially in younger, more karyotypically diverse clades, speciation is more frequently associated with cladogenetic chromosomal changes. Overall, our results imply that chromosomal fusions and fissions have contrasting macroevolutionary roles and that large-scale chromosomal rearrangements are associated with bursts of species diversification.

Pheromone-mediated command from the female to male clock induces and synchronizes circadian rhythms of the moth Spodoptera littoralis

To extract any adaptive benefit, the circadian clock needs to be synchronized to the 24-h day-night cycles. We have investigated if it is a general property of the brain's circadian clock to recognize social interactions as external time givers. Sociosexual interactions with the opposite sex are universal, prevalent even in the lives of solitary animals. The solitary adult life of the Spodoptera littoralis moth is singularly dedicated to sex, offering an ideal context for exploring the impact of sociosexual cues on circadian timekeeping. We have identified specific olfactory cues responsible for social entrainment, revealing a surprisingly strong influence of pheromone-mediated remote sociosexual interactions on circadian rhythms. Males' free-running rhythms are induced and synchronized by the sex pheromone that the female releases in a rhythmic fashion, highlighting a hierarchical relation between the female and male circadian oscillators. Even a single pulse of the sex pheromone altered clock gene expression in the male brain, surpassing the effect of light on the clock. Our finding of a daytime-dependent, lasting impact of pheromone on male's courtship efficacy indicates that circadian timing in moths is a trait under sexual selection. We have identified specific components of the sex-pheromone blend that lack mate-attractive property but have powerful circadian effects, providing rationale for their continued retention by the female. We show that such volatiles, when shared across sympatric moth species, can trigger communal synchronization. Our results suggest that the sex pheromone released by female moths entrains males' behavioral activity rhythm to ensure synchronized timing of mating.

Comparative genomics reveals the dynamics of chromosome evolution in Lepidoptera

Chromosomes are a central unit of genome organization. One-tenth of all described species on Earth are butterflies and moths, the Lepidoptera, which generally possess 31 chromosomes. However, some species display dramatic variation in chromosome number. Here we analyse 210 chromosomally complete lepidopteran genomes and show that the chromosomes of extant lepidopterans are derived from 32 ancestral linkage groups, which we term Merian elements. Merian elements have remained largely intact through 250 million years of evolution and diversification. Against this stable background, eight lineages have undergone extensive reorganization either through numerous fissions or a combination of fusion and fission events. Outside these lineages, fusions are rare and fissions are rarer still. Fusions often involve small, repeat-rich Merian elements and the sex-linked element. Our results reveal the constraints on genome architecture in Lepidoptera and provide a deeper understanding of chromosomal rearrangements in eukaryotic genome evolution.

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

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

An update and review of arthropod vector sensory systems: Potential targets for behavioural manipulation by parasites and other disease agents

For over a century, vector ecology has been a mainstay of vector-borne disease control. Much of this research has focused on the sensory ecology of blood-feeding arthropods (black flies, mosquitoes, ticks, etc.) with terrestrial vertebrate hosts. Of particular interest are the cues and sensory systems that drive host seeking and host feeding behaviours as they are critical for a vector to locate and feed from a host. An important yet overlooked component of arthropod vector ecology are the phenotypic changes observed in infected vectors that increase disease transmission. While our fundamental understanding of sensory mechanisms in disease vectors has drastically increased due to recent advances in genome engineering, for example, the advent of CRISPR-Cas9, and high-throughput "big data" approaches (genomics, proteomics, transcriptomics, etc.), we still do not know if and how parasites manipulate vector behaviour. Here, we review the latest research on arthropod vector sensory systems and propose key mechanisms that disease agents may alter to increase transmission.

Pervasive mimicry in flight behavior among aposematic butterflies

Flight was a key innovation in the adaptive radiation of insects. However, it is a complex trait influenced by a large number of interacting biotic and abiotic factors, making it difficult to unravel the evolutionary drivers. We investigate flight patterns in neotropical heliconiine butterflies, well known for mimicry of their aposematic wing color patterns. We quantify the flight patterns (wing beat frequency and wing angles) of 351 individuals representing 29 heliconiine and 9 ithomiine species belonging to ten color pattern mimicry groupings. For wing beat frequency and up wing angles, we show that heliconiine species group by color pattern mimicry affiliation. Convergence of down wing angles to mimicry groupings is less pronounced, indicating that distinct components of flight are under different selection pressures and constraints. The flight characteristics of the Tiger mimicry group are particularly divergent due to convergence with distantly related ithomiine species. Predator-driven selection for mimicry also explained variation in flight among subspecies, indicating that this convergence can occur over relatively short evolutionary timescales. Our results suggest that the flight convergence is driven by aposematic signaling rather than shared habitat between comimics. We demonstrate that behavioral mimicry can occur between lineages that have separated over evolutionary timescales ranging from <0.5 to 70 My.

Predator selection on phenotypic variability of cryptic and aposematic moths

Natural selection generally favours phenotypic variability in camouflaged organisms, whereas aposematic organisms are expected to evolve a more uniform warning coloration. However, no comprehensive analysis of the phenotypic consequences of predator selection in aposematic and cryptic species exists. Using state-of-the-art image analysis, we examine 2800 wing images of 82 moth species accessed via three online museum databases. We test whether anti-predator strategy (i.e., camouflage or aposematism) explains intraspecific variation in wing colour and pattern across northern hemisphere moths. In addition, we test two mutually non-exclusive, ecological hypotheses to explain variation in colour pattern: diel-activity or dietary-niche. In this work, taking into account phylogenetic relationships, moth phenotypic variability is best explained by anti-predator strategy with camouflaged moths being more variable in wing patterning than aposematic species.

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.

Integrated phylogenomic approaches in insect systematics

The increased accessibility of genomic and imaging methods, and the improved access to ecological, spatial, and other natural history-related data is allowing for insect systematics to grow and find answers to central evolutionary and taxonomic questions. Today, integrated studies in insect phylogenomics and systematics are combining natural history, behavior, developmental biology, morphology, fossils, geographic range data, and ecological interactions. This integration is contributing to the clarification of evolutionary relationships, and the recognition of the role played by these factors on the evolution of insects. Future work should continue to build on these advances, seeking to further increase open-access databasing and support for natural history research, as well as expand its analytical palettes.

150 Years of Coevolution Research: Evolution and Ecology of Yucca Moths (Prodoxidae) and Their Hosts

Yucca moths (Tegeticula and Parategeticula) are specialized pollinators of yucca plants, possessing unique, tentacle-like mouthparts used to actively collect pollen and deposit it onto the flowers of their hosts. The moths' larvae feed on the developing seeds and fruit tissue. First described in 1873, the yucca-yucca moth pollination system is now considered the archetypical example of a coevolved intimate mutualism. Research conducted over the past three decades has transformed our understanding of yucca moth diversity and host plant interactions. We summarize the current understanding of the diversity, ecology, and evolution of this group, review evidence for coevolution of the insects and their hosts, and describe how the nature of the interaction varies across evolutionary time and ecological contexts. Finally, we identify unresolved questions and areas for future research.

Gregariousness in lepidopteran larvae

The gregarious lifestyle of lepidopteran larvae is diverse and shaped by a complex interplay of ecological and evolutionary factors. Our review showed that the larval-aggregation behavior has been reported in 23 lepidopteran families, indicating multiple evolution of this behavior. Some larvae live in sibling groups throughout all larval instars and even pupation stages, which may result from the kin-selection. In contrast, group fusion may occur among different sibling or foraging groups of larvae and form larger aggregates, and the gregariousness of these species might be driven by the group-selection. While group size and foraging patterns vary greatly across species, it is generally associated with improved larval survivorship and accelerated development. However, the advantages of group living, such as facilitating feeding activities, adjusting the temperature, and defending natural enemies, may diminish along with development, with strong intraspecific competition occurring at later instars, even when food is abundant. Therefore, the group sizes and fission-fusion dynamics of certain gregarious lepidopteran larvae may be a consequence of their cost-benefit balance depending on various biotic and abiotic factors. Trail and aggregation pheromones, silk trails, or body contact contribute to collective movement and group cohesion of gregarious lepidopteran larvae. However, frequent contact among group members may cause the horizontal transmission of pathogens and pesticides, which may bring an integrated pest management strategy controlling gregarious lepidopteran pests.

Novel insight into lepidopteran phylogenetics from the mitochondrial genome of the apple fruit moth of the family Argyresthiidae

BACKGROUND

The order Lepidoptera has an abundance of species, including both agriculturally beneficial and detrimental insects. Molecular data has been used to investigate the phylogenetic relationships of major subdivisions in Lepidoptera, which has enhanced our understanding of the evolutionary relationships at the family and superfamily levels. However, the phylogenetic placement of many superfamilies and/or families in this order is still unknown. In this study, we determine the systematic status of the family Argyresthiidae within Lepidoptera and explore its phylogenetic affinities and implications for the evolution of the order. We describe the first mitochondrial (mt) genome from a member of Argyresthiidae, the apple fruit moth Argyresthia conjugella. The insect is an important pest on apples in Fennoscandia, as it switches hosts when the main host fails to produce crops.

RESULTS

The mt genome of A. conjugella contains 16,044 bp and encodes all 37 genes commonly found in insect mt genomes, including 13 protein-coding genes (PCGs), two ribosomal RNAs, 22 transfer RNAs, and a large control region (1101 bp). The nucleotide composition was extremely AT-rich (82%). All detected PCGs (13) began with an ATN codon and terminated with a TAA stop codon, except the start codon in cox1 is ATT. All 22 tRNAs had cloverleaf secondary structures, except trnS1, where one of the dihydrouridine (DHU) arms is missing, reflecting potential differences in gene expression. When compared to the mt genomes of 507 other Lepidoptera representing 18 superfamilies and 42 families, phylogenomic analyses found that A. conjugella had the closest relationship with the Plutellidae family (Yponomeutoidea-super family). We also detected a sister relationship between Yponomeutoidea and the superfamily Tineidae.

CONCLUSIONS

Our results underline the potential importance of mt genomes in comparative genomic analyses of Lepidoptera species and provide valuable evolutionary insight across the tree of Lepidoptera species.

Evolutionary genomics of three agricultural pest moths reveals rapid evolution of host adaptation and immune-related genes

BACKGROUND

Understanding the genotype of pest species provides an important baseline for designing integrated pest management (IPM) strategies. Recently developed long-read sequence technologies make it possible to compare genomic features of nonmodel pest species to disclose the evolutionary path underlying the pest species profiles. Here we sequenced and assembled genomes for 3 agricultural pest gelechiid moths: Phthorimaea absoluta (tomato leafminer), Keiferia lycopersicella (tomato pinworm), and Scrobipalpa atriplicella (goosefoot groundling moth). We also compared genomes of tomato leafminer and tomato pinworm with published genomes of Phthorimaea operculella and Pectinophora gossypiella to investigate the gene family evolution related to the pest species profiles.

RESULTS

We found that the 3 solanaceous feeding species, P. absoluta, K. lycopersicella, and P. operculella, are clustered together. Gene family evolution analyses with the 4 species show clear gene family expansions on host plant-associated genes for the 3 solanaceous feeding species. These genes are involved in host compound sensing (e.g., gustatory receptors), detoxification (e.g., ABC transporter C family, cytochrome P450, glucose-methanol-choline oxidoreductase, insect cuticle proteins, and UDP-glucuronosyl), and digestion (e.g., serine proteases and peptidase family S1). A gene ontology enrichment analysis of rapid evolving genes also suggests enriched functions in host sensing and immunity.

CONCLUSIONS

Our results of family evolution analyses indicate that host plant adaptation and pathogen defense could be important drivers in species diversification among gelechiid moths.

A major endogenous glycoside hydrolase mediating quercetin uptake in Bombyx mori

Quercetin is a common plant flavonoid which is involved in herbivore-plant interactions. Mulberry silkworms (domestic silkworm, Bombyx mori, and wild silkworm, Bombyx mandarina) take up quercetin from mulberry leaves and accumulate the metabolites in the cocoon, thereby improving its protective properties. Here we identified a glycoside hydrolase, named glycoside hydrolase family 1 group G 5 (GH1G5), which is expressed in the midgut and is involved in quercetin metabolism in the domestic silkworm. Our results suggest that this enzyme mediates quercetin uptake by deglycosylating the three primary quercetin glycosides present in mulberry leaf: rutin, quercetin-3-O-malonylglucoside, and quercetin-3-O-glucoside. Despite being located in an unstable genomic region that has undergone frequent structural changes in the evolution of Lepidoptera, GH1G5 has retained its hydrolytic activity, suggesting quercetin uptake has adaptive significance for mulberry silkworms. GH1G5 is also important in breeding: defective mutations which result in discoloration of the cocoon and increased silk yield are homozygously conserved in 27 of the 32 Japanese white-cocoon domestic silkworm strains and 12 of the 30 Chinese ones we investigated.

Dispersal from the Qinghai-Tibet plateau by a high-altitude butterfly is associated with rapid expansion and reorganization of its genome

Parnassius glacialis is a typical "Out of the QTP" alpine butterfly that originated on the Qinghai-Tibet Plateau (QTP) and dispersed into relatively low-altitude mountainous. Here we assemble a chromosome-level genome of P. glacialis and resequence 9 populations in order to explore the genome evolution and local adaptation of this species. These results indicated that the rapid accumulation and slow unequal recombination of transposable elements (TEs) contributed to the formation of its large genome. Several ribosomal gene families showed extensive expansion and selective evolution through transposon-mediated processed pseudogenes. Additionally, massive structural variations (SVs) of TEs affected the genetic differentiation of low-altitude populations. These low-altitude populations might have experienced a genetic bottleneck in the past and harbor genes with selective signatures which may be responsible for the potential adaptation to low-altitude environments. These results provide a foundation for understanding genome evolution and local adaptation for "Out of the QTP" of P. glacialis.

Genomics of the relict species Baronia brevicornis sheds light on its demographic history and genome size evolution across swallowtail butterflies

Relict species, like coelacanth, gingko, tuatara, are the remnants of formerly more ecologically and taxonomically diverse lineages. It raises the questions of why they are currently species-poor, have restrained ecology, and are often vulnerable to extinction. Estimating heterozygosity level and demographic history can guide our understanding of the evolutionary history and conservation status of relict species. However, few studies have focused on relict invertebrates compared to vertebrates. We sequenced the genome of Baronia brevicornis (Lepidoptera: Papilionidae), which is an endangered species, the sister species of all swallowtail butterflies, and is the oldest lineage of all extant butterflies. From a dried specimen, we were able to generate both long-read and short-read data and assembled a genome of 406 Mb for Baronia. We found a fairly high level of heterozygosity (0.58%) compared to other swallowtail butterflies, which contrasts with its endangered and relict status. Taking into account the high ratio of recombination over mutation, demographic analyses indicated a sharp decline of the effective population size initiated in the last million years. Moreover, the Baronia genome was used to study genome size variation in Papilionidae. Genome sizes are mostly explained by transposable elements activities, suggesting that large genomes appear to be a derived feature in swallowtail butterflies as transposable elements activity is recent and involves different transposable elements classes among species. This first Baronia genome provides a resource for assisting conservation in a flagship and relict insect species as well as for understanding swallowtail genome evolution.

Stealth echolocation in aerial hawking bats reflects a substrate gleaning ancestry

Predator-prey co-evolution can escalate into an evolutionary arms race.1 Examples of insect countermeasures to bat echolocation are well-known,2 but presumptive direct counter strategies in bats to insect anti-bat tactics are rare. The emission of very low-intensity calls by the hawking Barbastella barbastellus to circumvent high-frequency moth hearing is the most convincing countermeasure known.2,3 However, we demonstrate that stealth echolocation did not evolve through a high-intensity aerial hawking ancestor becoming quiet as previously hypothesized2,3,4 but from a gleaning ancestor transitioning into an obligate aerial hawker. Our ancestral state reconstructions show that the Plecotini ancestor likely gleaned prey using low-intensity calls typical of gleaning bats and that this ability-and associated traits-was subsequently lost in the barbastelle lineage. Barbastelles did not, however, revert to the oral, high-intensity call emission that other hawking bats use but retained the low-intensity nasal emission of closely related gleaning plecotines despite an extremely limited echolocation range. We further show that barbastelles continue to emit low-intensity calls even under adverse noise conditions and do not broaden the echolocation beam during the terminal buzz, unlike other vespertilionids attacking airborne prey.5,6 Together, our results suggest that barbastelles' echolocation is subject to morphological constraints prohibiting higher call amplitudes and beam broadening in the terminal buzz. We suggest that an abundance of eared prey allowed the co-opting and maintenance of low-intensity, nasal echolocation in today's obligate hawking barbastelle and that this unique foraging behavior7 persists because barbastelles remain a rare, acoustically inconspicuous predator to eared moths. VIDEO ABSTRACT.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

A highly conserved plant volatile odorant receptor detects a sex pheromone component of the greater wax moth, Galleria mellonella (Lepidoptera: Pyralidae)

Odorant receptors (ORs) are key specialized units for mate and host finding in moths of the Ditrysia clade, to which 98% of the lepidopteran species belong. Moth ORs have evolved to respond to long unsaturated acetates, alcohols, or aldehydes (Type I sex pheromones), falling into conserved clades of pheromone receptors (PRs). These PRs might have evolved from old lineages of non-Ditrysian moths that use plant volatile-like pheromones. However, a Ditrysian moth called the greater wax moth, Galleria mellonella (a worldwide-distributed pest of beehives), uses C9-C11 saturated aldehydes as the main sex pheromone components (i.e., nonanal and undecanal). Thus, these aldehydes represent unusual components compared with the majority of moth species that use, for instance, Type I sex pheromones. Current evidence shows a lack of consensus in the amount of ORs for G. mellonella, although consistent in that the moth does not have conserved PRs. Using genomic data, 62 OR candidates were identified, 16 being new genes. Phylogeny showed no presence of ORs in conserved PR clades. However, an OR with the highest transcript abundance, GmelOR4, appeared in a conserved plant volatile-detecting clade. Functional findings from the HEK system showed the OR as sensitive to nonanal and 2-phenylacetaldehyde, but not to undecanal. It is believed that to date GmelOR4 represents the first, but likely not unique, OR with a stable function in detecting aldehydes that help maintain the life cycle of G. mellonella around honey bee colonies.

Gene family expansion analysis and identification of the histone family in Spodoptera frugiperda

Spodoptera frugiperda, a major invasive pest, causes severe damage to various economically important crops. Previous comparative genomics studies have revealed a close association between the invasiveness of S. frugiperda and its genome. In recent years, a vast amount of genome from lepidopteran species has become available, offering an opportunity for a more detailed and comprehensive understanding of the biological characteristics of S. frugiperda. In this study, we conducted a comprehensive comparative genomics analysis of S. frugiperda using genome from 46 lepidopteran species. We found the highest number of gene family expansion events in S. frugiperda, indicating that gene family expansion is a crucial mechanism in its adaptive evolution. The expanded gene families are enriched in various biological processes, including nutrient metabolism, development, stress response, reproduction, and immune processes, suggesting that the expansion of these gene families likely contributes to the strong environmental adaptability of S. frugiperda. Furthermore, we identified the expansion of histone gene families in S. frugiperda which resulted from chromosome segmental duplications after the divergence from closely related species. Expression analysis of histone genes indicated that certain members might exert an influence on the growth and reproduction processes of S. frugiperda. Overall, our study deepens our understanding of the biological characteristics of S. frugiperda, providing a theoretical basis for the comprehensive management and sustained control of S. frugiperda and other lepidopteran pests in the future.

Acoustic camouflage increases with body size and changes with bat echolocation frequency range in a community of nocturnally active Lepidoptera

Body size is an important trait in predator-prey dynamics as it is often linked to detection, as well as the success of capture or escape. Larger prey, for example, often runs higher risk of detection by their predators, which imposes stronger selection on their anti-predator traits compared to smaller prey. Nocturnal Lepidoptera (moths) vary strongly in body size, which has consequences for their predation risk, as bigger moths return stronger echoes for echolocating bats. To compensate for increased predation risk, larger moths are therefore expected to have improved anti-predator defences. Moths are covered by different types of scales, which for a few species are known to absorb ultrasound, thus providing acoustic camouflage. Here, we assessed whether moths differ in their acoustic camouflage in a size-dependent way by focusing on their body scales and the different frequency ranges used by bats. We used a sonar head to measure 3D echo scans of a total of 111 moth specimens across 58 species, from eight different families of Lepidoptera. We scanned all the specimens and related their echo-acoustic target strength to various body size measurements. Next, we removed the scales covering the thorax and abdomen and scanned a subset of specimens again to assess the sound absorptive properties of these scales. Comparing intact specimens with descaled specimens, we found almost all species to absorb ultrasound, reducing detection risk on average by 8%. Furthermore, the sound absorptive capacities of body scales increased with body size suggesting that larger species benefit more from acoustic camouflage. The size-dependent effect of camouflage was in particular pronounced for the higher frequencies (above 29 kHz), with moth species belonging to large-bodied families consequently demonstrating similar target strengths compared to species from small-bodied families. Finally, we found the families to differ in frequency range that provided the largest reduction in detection risk, which may be related to differences in predation pressure and predator communities of these families. In general, our findings have important implications for predator-prey interactions across eco-evolutionary timescales and may suggest that acoustic camouflage played a role in body size evolution of nocturnally active Lepidoptera.

Elevated rates of dietary generalization in eusocial lineages of the secondarily herbivorous bees

BACKGROUND

Within the Hymenoptera, bees are notable for their relationship with flowering plants, being almost entirely dependent on plant pollen and nectar. Though functionally herbivorous, as a result of their role as pollinators, bees have received comparatively little attention as models for insect herbivory. Bees often display dietary specialization, but quantitative comparison against other herbivorous insects has not previously been conducted.

RESULTS

In the most comprehensive analysis to date for 860 bee species, dietary specialization amounted to 50.1% of studied species collecting pollen from between 1 and 2 botanical families with a relatively long tail of dietary generalists, with 11.1% of species collecting from more than 10 botanical families. This distribution deviated from the truncated Pareto distribution of dietary breadth seen in other herbivorous insect lineages. However, this deviation was predominantly due to eusocial bee lineages, which show a range of dietary breadths that conformed to a normal distribution, while solitary bees show a typical truncated distribution not strongly different from other herbivorous insects. We hypothesize that the relatively low level of dietary specialization in bees as a whole reflects the relaxation of the constraints typically observed in herbivorous insects with a comparatively reduced importance of plant chemistry and comparatively increased importance of phenology and foraging efficiency. The long flight periods of eusocial bees that are necessary to allow overlapping generations both allows and necessitates the use of multiple flowering resources, whereas solitary bees with short flight periods have more limited access to varied resources within a constrained activity period.

CONCLUSIONS

Collectively, solitary bees show slightly lower specialization compared to other herbivorous insects, possibly due to their balanced relationship with plants, rather than direct antagonism such as seen in the direct consumption of plant tissues. An additional factor may be the mediocre diversity of bees at low latitudes combined with low levels of dietary specialization, whereas these areas typically display a high rate of specialization by herbivorous insects in general. Though the most important factors structuring dietary specialization in bees appear to differ from many other herbivorous insects, solitary bees show a surprisingly similar overall pattern of dietary specialization.

Unraveling the complex evolutionary history of lepidopteran chromosomes through ancestral chromosome reconstruction and novel chromosome nomenclature

BACKGROUND

Lepidoptera is one of the most species-rich animal groups, with substantial karyotype variations among species due to chromosomal rearrangements. Knowledge of the evolutionary patterns of lepidopteran chromosomes still needs to be improved.

RESULTS

Here, we used chromosome-level genome assemblies of 185 lepidopteran insects to reconstruct an ancestral reference genome and proposed a new chromosome nomenclature. Thus, we renamed over 5000 extant chromosomes with this system, revealing the historical events of chromosomal rearrangements and their features. Additionally, our findings indicate that, compared with autosomes, the Z chromosome in Lepidoptera underwent a fast loss of conserved genes, rapid acquisition of lineage-specific genes, and a low rate of gene duplication. Moreover, we presented evidence that all available 67 W chromosomes originated from a common ancestor chromosome, with four neo-W chromosomes identified, including one generated by fusion with an autosome and three derived through horizontal gene transfer. We also detected nearly 4000 inter-chromosomal gene movement events. Notably, Geminin is transferred from the autosome to the Z chromosome. When located on the autosome, Geminin shows female-biased expression, but on the Z chromosome, it exhibits male-biased expression. This contributes to the sexual dimorphism of body size in silkworms.

CONCLUSIONS

Our study sheds light on the complex evolutionary history of lepidopteran chromosomes based on ancestral chromosome reconstruction and novel chromosome nomenclature.

Critical re-examination of known purported fossil Bombycoidea (Lepidoptera)

We critically re-examine 17 records of fossils currently assigned to the lepidopteran superfamily Bombycoidea, which includes the silk moths, emperor moths and hawk moths. These records include subfossils, compression and impression fossils, permineralizations and ichnofossils. We assess whether observable morphological features warrant their confident assignment to the superfamily. None of the examined fossils displays characters that allow unequivocal identification as Sphingidae, but three fossils and a subfossil (Mioclanis shanwangiana Zhang, Sun and Zhang, 1994, two fossil larvae, and a proboscis in asphaltum) have combinations of diagnostic features that support placement in the family. The identification of a fossil pupa as Bunaeini (Saturniidae) is well supported. The other fossils that we evaluate lack definitive bombycoid and, in several cases, even lepidopteran characters. Some of these dubious fossils have been used as calibration points in earlier studies casting doubt on the resulting age estimates. All fossil specimens reliably assigned to Bombycoidea are relatively young, the earliest fossil evidence of the superfamily dating to the middle Miocene.

UPrimer: A Clade-Specific Primer Design Program Based on Nested-PCR Strategy and Its Applications in Amplicon Capture Phylogenomics

Amplicon capture is a promising target sequence capture approach for phylogenomic analyses, and the design of clade-specific nuclear protein-coding locus (NPCL) amplification primers is crucial for its successful application. In this study, we developed a primer design program called UPrimer that can quickly design clade-specific NPCL amplification primers based on genome data, without requiring manual intervention. Unlike other available primer design programs, UPrimer uses a nested-PCR strategy that greatly improves the amplification success rate of the designed primers. We examined all available metazoan genome data deposited in NCBI and developed NPCL primer sets for 21 metazoan groups with UPrimer, covering a wide range of taxa, including arthropods, mollusks, cnidarians, echinoderms, and vertebrates. On average, each clade-specific NPCL primer set comprises ∼1,000 NPCLs. PCR amplification tests were performed in 6 metazoan groups, and the developed primers showed a PCR success rate exceeding 95%. Furthermore, we demonstrated a phylogenetic case study in Lepidoptera, showing how NPCL primers can be used for phylogenomic analyses with amplicon capture. Our results indicated that using 100 NPCL probes recovered robust high-level phylogenetic relationships among butterflies, highlighting the utility of the newly designed NPCL primer sets for phylogenetic studies. We anticipate that the automated tool UPrimer and the developed NPCL primer sets for 21 metazoan groups will enable researchers to obtain phylogenomic data more efficiently and cost-effectively and accelerate the resolution of various parts of the Tree of Life.

Opsin Gene Duplication in Lepidoptera: Retrotransposition, Sex Linkage, and Gene Expression

Color vision in insects is determined by signaling cascades, central to which are opsin proteins, resulting in sensitivity to light at different wavelengths. In certain insect groups, lineage-specific evolution of opsin genes, in terms of copy number, shifts in expression patterns, and functional amino acid substitutions, has resulted in changes in color vision with subsequent behavioral and niche adaptations. Lepidoptera are a fascinating model to address whether evolutionary change in opsin content and sequence evolution are associated with changes in vision phenotype. Until recently, the lack of high-quality genome data representing broad sampling across the lepidopteran phylogeny has greatly limited our ability to accurately address this question. Here, we annotate opsin genes in 219 lepidopteran genomes representing 33 families, reconstruct their evolutionary history, and analyze shifts in selective pressures and expression between genes and species. We discover 44 duplication events in opsin genes across ∼300 million years of lepidopteran evolution. While many duplication events are species or family specific, we find retention of an ancient long-wavelength-sensitive (LW) opsin duplication derived by retrotransposition within the speciose superfamily Noctuoidea (in the families Nolidae, Erebidae, and Noctuidae). This conserved LW retrogene shows life stage-specific expression suggesting visual sensitivities or other sensory functions specific to the early larval stage. This study provides a comprehensive order-wide view of opsin evolution across Lepidoptera, showcasing high rates of opsin duplications and changes in expression patterns.

Contrasting laboratory and field outcomes of bat-moth interactions

Predator-prey interactions are important but difficult to study in the field. Therefore, laboratory studies are often used to examine the outcomes of predator-prey interactions. Previous laboratory studies have shown that moth hearing and ultrasound production can help prey avoid being eaten by bats. We report here that laboratory behavioural outcomes may not accurately reflect the outcomes of field bat-moth interactions. We tested the success rates of two bat species capturing moths with distinct anti-bat tactics using behavioural experiments. We compared the results with the dietary composition of field bats using next-generation DNA sequencing. Rhinolophus episcopus and Rhinolophus osgoodi had a lower rate of capture success when hunting for moths that produce anti-bat clicks than for silent eared moths and earless moths. Unexpectedly, the success rates of the bats capturing silent eared moths and earless moths did not differ significantly from each other. However, the field bats had a higher proportion of silent eared moths than that of earless moths and that of clicking moths in their diets. The difference between the proportions of silent eared moths and earless moths in the bat diets can be explained by the difference between their abundance in bat foraging habitats. These findings suggest that moth defensive tactics, bat countertactics and moth availability collectively shape the diets of insectivorous bats. This study illustrates the importance of using a combination of behavioural experiments and molecular genetic techniques to reveal the complex interactions between predators and prey in nature.

Natural history and morphology of immature stages of Tolype medialis (Jones, 1912) (Lepidoptera: Lasiocampidae: Macromphaliinae)

Lasiocampidae belongs to superfamily Lasiocampoidea and contains more than a thousand species nearly distributed worldwide. Despite the great species richness and wide distribution, this group has internal phylogenetic relationships still little explored and with few studies on the morphology and biology of its immatures. This study describes the immature stages of the neotropical species Tolype medialis (Jones, 1912), focusing on the morphology and natural history. The eggs of T. medialis are oviposited freely inside a conical structure, and the larvae showed gregarious behavior in all instars. The seventh and eighth instar bear a pair of abdominal rounded flattened reddish brown glands on the segments A1, A2, A7, and A8 that produce a wax-like secretion that covers the pupae and the internal walls of the cocoon. In order to add information to the Lasiocampidae family, we compare and discuss these and other traits from the morphology and natural history of T. medialis immatures.

Regional Pollination Activity by Moth Migration in Athetis lepigone

Nocturnal moths (Lepidoptera) are important pollinators of a wide range of plant species. Understanding the foraging preferences of these insects is essential for their scientific management. However, this information is lacking for most moth species. The present study was therefore conducted to delineate the host plant feeding and pollination ranges of an agriculturally important nocturnal moth species Athetis lepigone by identifying the pollen species adhering to their bodies during long-distance migration. Pollen grains were dislodged from 1871 A. lepigone migrants captured on Beihuang Island in the Bohai Strait between 2020 and 2021. This region is a key seasonal migration pathway for A. lepigone in northern China. Almost 20% of all moths sampled harbored pollens, providing direct evidences that this moth species may serve as pollinators. Moreover, at least 39 pollen taxa spanning 21 plant families and 31 genera were identified, with a preference for Asteraceae, Amaranthaceae, and Pinaceae. Additionally, the pollen adherence ratios and taxa varied with moth sex, inter-annual changes, and seasonal fluctuations. Most importantly, the pollen taxa were correlated with insect migration stages and indicated that A. lepigone bidirectionally migrates between central China (Shandong, Hebei, and Henan Provinces) and northeastern China (Liaoning Province). Overall, the findings of the present work provide valuable information on the pollination behavior, geographical origins, and pollination regions of A. lepigone moths and could facilitate the design and optimization of efficacious local and regional management strategies for this important insect.

Calibrated microphone array recordings reveal that a gleaning bat emits low-intensity echolocation calls even in open-space habitat

Echolocating bats use ultrasound for orientation and prey capture in darkness. Ultrasound is strongly attenuated in air. Consequently, aerial-hawking bats generally emit very intense echolocation calls to maximize detection range. However, call levels vary more than tenfold (>20 dB) between species and are tightly linked to the foraging strategy. The brown long-eared bat (Plecotus auritus) is a primarily gleaning, low-amplitude species that may occasionally hawk airborne prey. We used state-of-the-art calibrated acoustic 3D-localization and automated call analysis to measure P. auritus' source levels. Plecotus auritus emits echolocation calls of low amplitude (92 dB rmsSPL re. 20 µPa at 10 cm) even while flying in open-space. While P. auritus thus probably benefits from delayed evasive manoeuvres of eared insects, we propose that low-amplitude echolocation did not evolve as an adaptive countermeasure, but is limited by morphological constraints.

Loss of Binding Capabilities in an Ecologically Important Odorant Receptor of the Fall Armyworm, Spodoptera frugiperda, by a Single Point Mutation

Olfaction plays a crucial role in locating food sources, mates, and spawning sites in the fall armyworm (FAW), Spodoptera frugiperda (Lepidoptera: Noctuidae). In the current study, SfruOR14, a highly conserved odorant receptor (OR) in lepidopteran species, was newly uncovered in S. frugiperda. In two-electrode voltage clamp recordings, the SfruOR14/Orco complex was narrowly tuned to six volatile compounds including phenylacetaldehyde (PAA), benzaldehyde, heptaldehyde, (E)-2-hexen-1-al, cinnamaldehyde, and 2-phenylethanol, among which PAA showed the strongest binding affinity. Subsequent homology modeling and molecular docking revealed that Phe79, His83, Tyr149, Pro176, Gln177, Leu202, and Thr348 in SfruOR14 were the key binding residues against the six ligands. Finally, as a result of site-directed mutagenesis, the SfruOR14His83Ala mutant completely lost its binding capabilities toward all ligands. Taken together, our findings provide valuable insights into understanding the interaction between SfruOR14 and the chemical ligands including PAA, which can help to design novel olfactory modulators for pest control.

Chromosomal-level Genome Assembly of the Coffee Bee Hawk Moth Reveals the Evolution of Chromosomes and the Molecular Basis of Distinct Phenotypes

Cephonodes hylas, the coffee bee hawk moth is a hawk moth species with unique characteristics, such as larvae feeding on gardenia, overcoming the toxicity of its iridoid glycosides, diurnal adults, and transparent wings. Although C. hylas is a fascinating model for molecular biological research, genome sequence analysis-based genetic approaches to elucidate these peculiarities have not yet been undertaken. We successfully achieved de novo genome assembly at the chromosome level of C. hylas comparable to the Lepidoptera model organism, silkworm. Additionally, 16,854 protein-coding genes were annotated, and the constructed genome sequence and annotated genes were of the highest quality BUSCO completion compared to closely related species. Comparative genome analysis revealed the process of chromosomal evolution from the Bombycoidea ancestral (n = 31) genome and changes in turnover at the chromosome level associated with chromosomal fusion events, such as the rate of repetitive sequence insertion. These analyses were only possible because the genome was constructed at the chromosome level. Additionally, increased the nonsynonymous/synonymous rate (dN/dS) ratios were observed in multiple photoreceptor-related genes that were strongly associated with the acquisition of diurnal activity. Furthermore, tandemly duplicated expanded genes containing many digestive and other enzymes and larval midgut-specific expression were also confirmed. These genes may be involved in the metabolism of genipin, a toxin found in gardenias. Using the genome sequence of C. hylas determined at the chromosome level, we have successfully identified new insights into the chromosomal evolution of Bombycoidea, as well as the relationship between the genome sequence and its characteristic traits.

Divergence and convergence of gut microbiomes of wild insect pollinators

Pollination services provided by wild insect pollinators are critical to natural ecosystems and crops around the world. There is an increasing appreciation that the gut microbiota of these insects influences their health and consequently their services. However, pollinator gut microbiota studies have focused on well-described social bees, but rarely include other, more phylogenetically divergent insect pollinators. To expand our understanding, we explored the insect pollinator microbiomes across three insect orders through two DNA sequencing approaches. First, in an exploratory 16S amplicon sequencing analysis of taxonomic community assemblages, we found lineage-specific divergences of dominant microbial genera and microbiota community composition across divergent insect pollinator genera. However, we found no evidence for a strong broad-scale phylogenetic signal, which we see for community relatedness at finer scales. Subsequently, we utilized metagenomic shotgun sequencing to obtain metagenome-assembled genomes and assess the functionality of the microbiota from pollinating flies and social wasps. We uncover a novel gut microbe from pollinating flies in the family Orbaceae that is closely related to Gilliamella spp. from social bees but with divergent functions. We propose this novel species be named Candidatus Gilliamella eristali. Further metagenomes of dominant fly and wasp microbiome members suggest that they are largely not host-insect adapted and instead may be environmentally derived. Overall, this study suggests selective processes involving ecology or physiology, or neutral processes determining microbe colonization may predominate in the turnover of lineages in insect pollinators broadly, while evolution with hosts may occur only under certain circumstances and on smaller phylogenetic scales. IMPORTANCE Wild insect pollinators provide many key ecosystem services, and the microbes associated with these insect pollinators may influence their health. Therefore, understanding the diversity in microbiota structure and function, along with the potential mechanisms shaping the microbiota across diverse insect pollinators, is critical. Our study expands beyond existing knowledge of well-studied social bees, like honey bees, including members from other bee, wasp, butterfly, and fly pollinators. We infer ecological and evolutionary factors that may influence microbiome structure across diverse insect pollinator hosts and the functions that microbiota members may play. We highlight significant differentiation of microbiomes among diverse pollinators. Closer analysis suggests that dominant members may show varying levels of host association and functions, even in a comparison of closely related microbes found in bees and flies. This work suggests varied importance of ecological, physiological, and non-evolutionary filters in determining structure and function across largely divergent wild insect pollinator microbiomes.

A mosaic of endogenous and plant-derived courtship signals in moths

Insects rely on olfaction to guide a wide range of adaptive behaviors, including mate and food localization, mate choice, oviposition site selection, kin recognition, and predator avoidance.1 In nocturnal insects, such as moths2 and cockroaches,3 mate finding is stimulated predominantly by long-range species-specific sex pheromones, typically emitted by females. During courtship, at close range, males in most moth species emit a blend of pheromone compounds from an everted, often large, pheromone gland. While long-distance communication with sex pheromones has been remarkably well characterized in thousands of moth species,2,4 close-range chemosensory sexual communication remains poorly understood. We reveal that in the moth Chloridea virescens, the male pheromone consists of three distinct classes of compounds: de novo biosynthesized alcohols, aldehydes, acetates, and carboxylic acids that resemble the female's emissions; newly identified compounds that are unique to the male pheromone, such as aliphatic polyunsaturated hydrocarbons; and sequestered plant secondary compounds and hormone derivatives, including methyl salicylate (MeSA). Thus, males employ a mosaic pheromone blend of disparate origins that may serve multiple functions during courtship. We show that two olfactory receptors in female antennae are tuned to MeSA, which facilitates female acceptance of the male. Because MeSA is emitted by plants attacked by pathogens and herbivores,5 the chemosensory system of female moths was likely already tuned to this plant volatile, and males appear to exploit the female's preadapted sensory bias. Interestingly, while female moths (largely nocturnal) and butterflies (diurnal) diverged in their use of sensory modalities in sexual communication,6 MeSA is used by males of both lineages.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Lepidopteran Synteny Units reveal deep chromosomal conservation in butterflies and moths

DNA is compacted into individual particles or chromosomes that form the basic units of inheritance. However, different animals and plants have widely different numbers of chromosomes. This means that we cannot readily tell which chromosomes are related to which. Here, we describe a simple technique that looks at the similarity of genes on each chromosome and thus gives us a true picture of their homology or similarity through evolutionary time. We use this new system to look at the chromosomes of butterflies and moths or Lepidoptera. We term the associated synteny units, Lepidopteran Synteny Units (LSUs). Using a sample of butterfly and moth genomes from across evolutionary time, we show that LSUs form a simple and reliable method of tracing chromosomal homology back through time. Surprisingly, this technique reveals that butterfly and moth chromosomes show conserved blocks dating back to their sister group the Trichoptera. As Lepidoptera have holocentric chromosomes, it will be interesting to see if similar levels of synteny are shown in groups of animals with monocentric chromosomes. The ability to define homology via LSU analysis makes it considerably easier to approach many questions in chromosomal evolution.

Scoutknife: A naïve, whole genome informed phylogenetic robusticity metric

Background: The phylogenetic bootstrap, first proposed by Felsenstein in 1985, is a critically important statistical method in assessing the robusticity of phylogenetic datasets. Core to its concept was the use of pseudo sampling - assessing the data by generating new replicates derived from the initial dataset that was used to generate the phylogeny. In this way, phylogenetic support metrics could overcome the lack of perfect, infinite data. With infinite data, however, it is possible to sample smaller replicates directly from the data to obtain both the phylogeny and its statistical robusticity in the same analysis. Due to the growth of whole genome sequencing, the depth and breadth of our datasets have greatly expanded and are set to only expand further. With genome-scale datasets comprising thousands of genes, we can now obtain a proxy for infinite data. Accordingly, we can potentially abandon the notion of pseudo sampling and instead randomly sample small subsets of genes from the thousands of genes in our analyses. Methods: We introduce Scoutknife, a jackknife-style subsampling implementation that generates 100 datasets by randomly sampling a small number of genes from an initial large-gene dataset to jointly establish both a phylogenetic hypothesis and assess its robusticity. We assess its effectiveness by using 18 previously published datasets and 100 simulation studies. Results: We show that Scoutknife is conservative and informative as to conflicts and incongruence across the whole genome, without the need for subsampling based on traditional model selection criteria. Conclusions: Scoutknife reliably achieves comparable results to selecting the best genes on both real and simulation datasets, while being resistant to the potential biases caused by selecting for model fit. As the amount of genome data grows, it becomes an even more exciting option to assess the robusticity of phylogenetic hypotheses.

Base Composition, Codon Usage, and Patterns of Gene Sequence Evolution in Butterflies

Coding sequence evolution is influenced by both natural selection and neutral evolutionary forces. In many species, the effects of mutation bias, codon usage, and GC-biased gene conversion (gBGC) on gene sequence evolution have not been detailed. Quantification of how these forces shape substitution patterns is therefore necessary to understand the strength and direction of natural selection. Here, we used comparative genomics to investigate the association between base composition and codon usage bias on gene sequence evolution in butterflies and moths (Lepidoptera), including an in-depth analysis of underlying patterns and processes in one species, Leptidea sinapis. The data revealed significant G/C to A/T substitution bias at third codon position with some variation in the strength among different butterfly lineages. However, the substitution bias was lower than expected from previously estimated mutation rate ratios, partly due to the influence of gBGC. We found that A/T-ending codons were overrepresented in most species, but there was a positive association between the magnitude of codon usage bias and GC-content in third codon positions. In addition, the tRNA-gene population in L. sinapis showed higher GC-content at third codon positions compared to coding sequences in general and less overrepresentation of A/T-ending codons. There was an inverse relationship between synonymous substitutions and codon usage bias indicating selection on synonymous sites. We conclude that the evolutionary rate in Lepidoptera is affected by a complex interaction between underlying G/C -> A/T mutation bias and partly counteracting fixation biases, predominantly conferred by overall purifying selection, gBGC, and selection on codon usage.