A large-scale, higher-level, molecular phylogenetic study of the insect order Lepidoptera (moths and butterflies)

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.

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.

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.

Molecular phylogeny, divergence time, biogeography and trends in host plant usage in the agriculturally important tortricid tribe Grapholitini (Lepidoptera: Tortricidae: Olethreutinae)

The leaf-roller moth tribe Grapholitini comprises about 1200 described species and contains numerous notorious pests of fruits and seeds. The phylogeny of the tribe has been little studied using contemporary methods, and the monophyly of several genera remains questionable. In order to provide a more robust phylogenetic framework for the group, we conducted a multiple-gene phylogenetic analysis of 104 species representing 27 genera of Grapholitini and 29 outgroup species. Divergence time, ancestral area, and host plant usage were also inferred to explore evolutionary trends in the tribe. Our analyses indicate that Larisa and Corticivora, traditionally assigned to Grapholitini, are best excluded from the tribe. After removal of these two genera, the tribe is found to be monophyletic, represented by two major lineages-a Dichrorampha clade and a Cydia clade, the latter of which can be divided into seven generic groups. The genus Grapholita was found to be polyphyletic, comprising three different clades, and we propose three genera to accommodate these groups: Grapholita (sensu stricto), Aspila (formerly a subgenus of Grapholita) and Ephippiphora (formerly considered a synonym of Grapholita). We summarize each generic group, including related genera not included in our analysis, providing morphological, pheromone and food plant characters that support particular branches within the molecular hypotheses. Biogeographical analyses indicate that Grapholitini probably originated in the Nearctic, Afrotropical and Neotropical regions in the Lutetian of the middle Eocene (ca. 44.3 Ma). Our results also indicate that most groups in Grapholitini originated from Fabaceae-feeding monophagous or oligophagous ancestors, and that host plant shifts probably promoted species diversification within the tribe.

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.

The genome and sex-dependent responses to temperature in the common yellow butterfly, Eurema hecabe

BACKGROUND

Lepidoptera (butterflies and moths) is one of the most geographically widespread insect orders in the world, and its species play important and diverse ecological and applied roles. Climate change is one of the biggest challenges to biodiversity this century, and lepidopterans are vulnerable to climate change. Temperature-dependent gene expression differences are of relevance under the ongoing climate crisis. However, little is known about how climate affects gene expression in lepidopterans and the ecological consequences of this, particularly with respect to genes with biased expression in one of the sexes. The common yellow butterfly, Eurema hecabe (Family Pieridae), is one of the most geographically widespread lepidopterans that can be found in Asia, Africa, and Australia. Nevertheless, what temperature-dependent effects there may be and whether the effects differ between the sexes remain largely unexplored.

RESULTS

Here, we generated high-quality genomic resources for E. hecabe along with transcriptomes from eight developmental stages. Male and female butterflies were subjected to varying temperatures to assess sex-specific gene expression responses through mRNA and microRNA transcriptomics. We find that there are more temperature-dependent sex-biased genes in females than males, including genes that are involved in a range of biologically important functions, highlighting potential ecological impacts of increased temperatures. Further, by considering available butterfly data on sex-biased gene expression in a comparative genomic framework, we find that the pattern of sex-biased gene expression identified in E. hecabe is highly species-specific, rather than conserved across butterfly species, suggesting that sex-biased gene expression responses to climate change are complex in butterflies.

CONCLUSIONS

Our study lays the foundation for further understanding of differential responses to environmental stress in a widespread lepidopteran model and demonstrates the potential complexity of sex-specific responses of lepidopterans to climate change.

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.

The genome sequence of the Feathered Bright, Incurvaria masculella (Denis & Schiffermüller, 1775)

We present a genome assembly from an individual male Incurvaria masculella (the Feathered Bright; Arthropoda; Insecta; Lepidoptera; Incurvariidae). The genome sequence is 552 megabases in span. Most of the assembly is scaffolded into 26 chromosomal pseudomolecules, including the assembled Z sex chromosome. The mitochondrial genome has also been assembled and is 15.3 kilobases in length.

Comparative genome and phylogenetic analysis revealed the complex mitochondrial genome and phylogenetic position of Conopomorpha sinensis Bradley

Conopomorpha sinensis Bradley is a destructive pest that causes severe economic damage to litchi and longan. Previous C. sinensis research has focused on population life tables, oviposition selectivity, pest population prediction, and control technology. However, there are few studies on its mitogenome and phylogenetic evolution. In this study, we sequenced the whole mitogenome of C. sinensis by the third-generation sequencing, and analyzed the characteristics of its mitogenome by comparative genome. The complete mitogenome of C. sinensis is a typical circular and double-stranded structure. The ENC-plot analyses revealed that natural selection could affect the information of codon bias of the protein-coding genes in the mitogenome of C. sinensis in the evolutionary process. Compared with 12 other Tineoidea species, the trnA-trnF gene cluster of tRNA in the C. sinensis mitogenome appears to have a new arrangement pattern. This new arrangement has not been found in other Tineoidea or other Lepidoptera, which needs further exploration. Meanwhile, a long AT repeated sequence was inserted between trnR and trnA, trnE and trnF, ND1 and trnS in the mitogenome of C. sinensis, and the reason for this sequence remains to be further studied. Furthermore, the results of phylogenetic analysis showed that the litchi fruit borer belonged to Gracillariidae, and Gracillariidae was monophyletic. The results will contribute to an improved understanding of the complex mitogenome and phylogeny of C. sinensis. It also will provide a molecular basis for further research on the genetic diversity and population differentiation of C. sinensis.

The mitochondrial genomes of the Geometroidea (Lepidoptera) and their phylogenetic implications

The Geometroidea is a large superfamily of Lepidoptera in species composition and contains numerous economically important pest species that cause great loss in crop and forest production. However, understanding of mitogenomes remains limited due to relatively fewer mitogenomes previously reported for this megadiverse group. Here, we sequenced and annotated nine mitogenomes for Geometridae and further analyzed the mitogenomic evolution and phylogeny of the whole superfamily. All nine mitogenomes contained 37 mitochondrial genes typical in insects, and gene organization was conserved except for Somatina indicataria. In S. indicataria, the positions of two tRNAs were rearranged. The trnR was located before trnA instead of after trnA typical in Lepidoptera, whereas the trnE was detected rarely on the minority strand (N-strand). This trnR-trnA-trnN-trnS1-trnE-trnF newly recognized in S. indicataria represents the first gene rearrangement reported for Geometroidea and is also unique in Lepidoptera. Besides, nucleotide composition analyses showed little heterogeneity among the four geometrid subfamilies involved herein, and overall, nad6 and atp8 have higher nucleotide diversity and Ka/Ks rate in Geometridae. In addition, the taxonomic assignments of the nine species, historically defined by morphological studies, were confirmed by various phylogenetic analyses based on the hitherto most extensive mitogenomic sampling in Geometroidea.

The genome sequence of the White-barred Gold, Micropterix aruncella (Scopoli, 1763)

We present a genome assembly from an individual female Micropterix aruncella (the White-barred Gold; Arthropoda, Insecta, Lepidoptera; Micropterigidae). The genome sequence is 1,079 megabases in span. Most of the assembly is scaffolded into 31 chromosomal pseudomolecules, including the assembled Z sex chromosome. The mitochondrial genome has also been assembled and is 15.0 kilobases in length.

Characterization of four mitochondrial genomes of Crambidae (Lepidoptera, Pyraloidea) and phylogenetic implications

Loxostege turbidalis, Loxostege aeruginalis, Pyrausta despicata, and Crambus perlellus belong to Crambidae, Pyraloidea. Their mitochondrial genomes (mitogenomes) were successfully sequenced. The mitogenomes of L. turbidalis, L. aeruginalis, P. despicata, and C. perlellus are 15 240 bp, 15 339 bp, 15 389 bp, and 15 440 bp. The four mitogenomes all have a typical insect mitochondrial gene order, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and one A + T rich region (control region). The PCGs are initiated by the typical ATN codons, except CGA for the cox1 gene. Most PCGs terminate with common codon TAA or TAG, the incomplete codon T is found as the stop codon for cox2, nad4, and nad5. Most tRNA genes exhibit typical cloverleaf structure, except trnS1 (AGN) lacking the dihydrouridine (DHU) arm. The secondary structure of rRNA of four mitogenomes were predicted. Poly-T structure and micro-satellite regions are conserved in control regions. The phylogenetic analyses based on 13 PCGs showed the relationships of subfamilies in Pyraloidea. Pyralidae, and Crambidae are monophyletic, respectively. Pyralidae comprises four subfamilies, which form the following topology with high support values: (Galleriinae + ((Pyralinae + Epipaschiinae)+ Phycitinae)). Crambidae includes seven subfamilies and is divided into two lineages. Pyraustinae and Spilomelinae are sister groups of each other, and form the "PS clade." Other five subfamilies (Crambinae, Acentropinae, Scopariinae, Schoenobiinae, and Glaphyriinae) form the "non-PS clade" in the Bayesian inference tree. However, Schoenobiinae is not grouped with the other four subfamilies and located at the base of Crambidae in two maximum likelihood trees.

Characterization of Seventeen Complete Mitochondrial Genomes: Structural Features and Phylogenetic Implications of the Lepidopteran Insects

Lepidoptera (moths and butterflies) are widely distributed in the world, but high-level phylogeny in Lepidoptera remains uncertain. More mitochondrial genome (mitogenome) data can help to conduct comprehensive analysis and construct a robust phylogenetic tree. Here, we sequenced and annotated 17 complete moth mitogenomes and made comparative analysis with other moths. The gene order of trnM-trnI-trnQ in 17 moths was different from trnI-trnQ-trnM of ancestral insects. The number, type, and order of genes were consistent with reported moths. The length of newly sequenced complete mitogenomes ranged from 14,231 bp of Rhagastis albomarginatus to 15,756 bp of Numenes albofascia. These moth mitogenomes were typically with high A+T contents varied from 76.0% to 81.7% and exhibited negative GC skews. Among 13 protein coding genes (PCGs), some unusual initiations and terminations were found in part of newly sequenced moth mitogenomes. Three conserved gene-overlapping regions and one conserved intergenic region were detected among 17 mitogenomes. The phylogenetic relationship of major superfamilies in Macroheterocera was as follows: (Bombycoidea + Lasiocampoidea) + ((Drepanoidea + Geometroidea) + Noctuoidea)), which was different from previous studies. Moreover, the topology of Noctuoidea as (Notodontidae + (Erebidae + Noctuidae)) was supported by high Bayesian posterior probabilities (BPP = 1.0) and bootstrapping values (BSV = 100). This study greatly enriched the mitogenome database of moth and strengthened the high-level phylogenetic relationships of Lepidoptera.

Enhanced Resolution of Evolution and Phylogeny of the Moths Inferred from Nineteen Mitochondrial Genomes

The vast majority (approximately 90%) of Lepidoptera species belong to moths whose phylogeny has been widely discussed and highly controversial. For the further understanding of phylogenetic relationships of moths, nineteen nearly complete mitochondrial genomes (mitogenomes) of moths involved in six major lineages were sequenced and characterized. These mitogenomes ranged from 15,177 bp (Cyclidia fractifasciata) to 15,749 bp (Ophthalmitis albosignaria) in length, comprising of the core 37 mitochondrial genes (13 protein-coding genes (PCGs) + 22 tRNAs + two rRNAs) and an incomplete control region. The order and orientation of genes showed the same pattern and the gene order of trnM-trnI-trnQ showed a typical rearrangement of Lepidoptera compared with the ancestral order of trnI-trnQ-trnM. Among these 13 PCGs, ATP8 exhibited the fastest evolutionary rate, and Drepanidae showed the highest average evolutionary rate among six families involved in 66 species. The phylogenetic analyses based on the dataset of 13 PCGs suggested the relationship of (Notodontidae + (Noctuidae + Erebidae)) + (Geometridae + (Sphingidae + Drepanidae)), which suggested a slightly different pattern from previous studies. Most groups were well defined in the subfamily level except Erebidae, which was not fully consistent across bayesian and maximum likelihood methods. Several formerly unassigned tribes of Geometridae were suggested based on mitogenome sequences despite a not very strong support in partial nodes. The study of mitogenomes of these moths can provide fundamental information of mitogenome architecture, and the phylogenetic position of moths, and contributes to further phylogeographical studies and the biological control of pests.

Movement and olfactory signals: Sexually dimorphic antennae and female flightlessness in moths

Darwin argued a role for sexual selection in the evolution of male sensory structures, including insect antennae, the strength of which will depend upon the importance of early arrival at receptive females. There is remarkable variation in the nature and degree of sexual dimorphism in moth antennae, with males of some species having spectacular, feathery antennae. Although it is widely assumed that these elaborate structures provide greater sensitivity to chemical signals (sex pheromones), the factors underlying the interspecific diversity in male antennal structure and size are poorly understood. Because male antennal morphology may be affected by several female life–history traits, including flight ability, we conducted a phylogenetic comparative analysis to test how these traits are linked, using data from 93 species of moths across 11 superfamilies. Our results reveal that elaborate antennae in males have evolved more frequently in species where females are monandrous. Further, female loss of flight ability evolved more frequently in species where males have elaborate antennae. These results suggest that elaborate antennae have evolved in response to more intense male competition, arising from female monandry, and that the evolution of elaborate antennae in males has, in turn, shaped the evolution of female flightlessness.

Comparative mitochondrial genome analysis and phylogenetic relationship among lepidopteran species

Lepidoptera has rich species including many agricultural pests and economical insects around the world. The mitochondrial genomes (mitogenomes) were utilized to explore the phylogenetic relationships between difference taxonomic levels in Lepidoptera. However, the knowledge of mitogenomic characteristics and phylogenetic position about superfamily-level in this order is unresolved. In this study, we integrated 794 mitogenomes consisting of 37 genes and a noncoding control region, which covered 26 lepidopteran superfamilies from newly sequenced and publicly available genomes for comparative genomic and phylogenetic analysis. In primitive taxon, putative start codon of cox1 gene was ATA or ATT instead of CGA, but stop codon of that showed four types, namely TAA, TAG, TA and T. The 7-bp overlap between atp8 and atp6 presented as "ATGATAA". Moreover, the most frequently utilized amino acids were leucine (UUA) in 13 PCGs. Phylogenetic analysis showed that the main backbone relationship in Lepidoptera was (Hepialoidea + (Nepticuloidea + (Adeloidea + (Tischerioidea + (Tineoidea + (Yponomeutoidea + (Gracillarioidea + (Papilionoidea + ((Zygaenoidea + Tortricoidea) + (Gelechioidea + (Pyraloidea + ((Geometroidea + Noctuoidea) + (Lasiocampoidea + Bombycoidea))))))))))))).

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

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

The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture

Simple Summary

In this study, we produced a genomic resource for the light brown apple moth, Epiphyas postvittana, to understand the biological basis of adaptation to a high number of hosts (polyphagy) and the invasive nature of this and other lepidopteran pests. The light brown apple moth is an invasive pest of horticultural plants, with over 500 recorded plant hosts. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, causing significant economic losses for fruit producers. Comparative genomic analyses with other lepidopteran genomes indicate that a high proportion of the genome is made up of repetitive sequences, with the majority of the known elements being DNA transposable elements and retrotransposons. Twenty gene families show significant expansions, including some likely to have a role in its pest status. Finally, population genomics, investigated by a RAD-tag approach, indicated likely patterns of invasion and admixture, with Californian moths most probably being derived from Australia.

Abstract

The light brown apple moth, Epiphyas postvittana is an invasive, polyphagous pest of horticultural systems around the world. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, where it has been found on over 500 plants, including many horticultural crops. We have produced a genomic resource, to understand the biological basis of the polyphagous and invasive nature of this and other lepidopteran pests. The assembled genome sequence encompassed 598 Mb and has an N50 of 301.17 kb, with a BUSCO completion rate of 97.9%. Epiphyas postvittana has 34% of its assembled genome represented as repetitive sequences, with the majority of the known elements made up of longer DNA transposable elements (14.07 Mb) and retrotransposons (LINE 17.83 Mb). Of the 31,389 predicted genes, 28,714 (91.5%) were assigned to 11,438 orthogroups across the Lepidoptera, of which 945 were specific to E. postvittana. Twenty gene families showed significant expansions in E. postvittana, including some likely to have a role in its pest status, such as cytochrome p450s, glutathione-S-transferases and UDP-glucuronosyltransferases. Finally, using a RAD-tag approach, we investigated the population genomics of this pest, looking at its likely patterns of invasion.

A supertree of Northern European macromoths

Ecological and life-history data on the Northern European macromoth (Lepidoptera: Macroheterocera) fauna is widely available and ideal for use in answering phylogeny-based research questions: for example, in comparative biology. However, phylogenetic information for such studies lags behind. Here, as a synthesis of all currently available phylogenetic information on the group, we produce a supertree of 114 Northern European macromoth genera (in four superfamilies, with Geometroidea considered separately), providing the most complete phylogenetic picture of this fauna available to date. In doing so, we assess those parts of the phylogeny that are well resolved and those that are uncertain. Furthermore, we identify those genera for which phylogenetic information is currently too poor to include in such a supertree, or entirely absent, as targets for future work. As an aid to studies involving these genera, we provide information on their likely positions within the macromoth tree. With phylogenies playing an ever more important role in the field, this supertree should be useful in informing future ecological and evolutionary studies.

Global Nepticulidae, Opostegidae, and Tischeriidae (Lepidoptera): temporal dynamics of species descriptions and their authors

This study identifies the number of named and described species of three monotrysian, plant-mining lepidopteran families worldwide: Nepticulidae and Opostegidae (Nepticuloidea), and Tischeriidae (Tischerioidea). At the end of 2021, we estimated that a total of 1000 Nepticulidae species, 197 Opostegidae species, and 170 Tischeriidae species have been described since the taxonomic practice of describing species began in the 18th century. We examine and discuss the history of descriptions and authorship of species worldwide for each of the three families. We found that the total (accumulative) number of species described increased with each time period delineated. About five new species were described per year on average, or about 22 new species were described per year in the 21st century. We recognize researchers with the most number of described species in these three taxa.  

?Description of a new species of the genus Neopseustis Meyrick, 1909 from China, with a new classification of the genus (Lepidoptera, Neopseustoidea, Neopseustidae)

A new species of the genus Neopseustis Meyrick, 1909, Neopseustischentangensis S.Y. Huang & Chen sp. nov., which was confirmed by both morphological and molecular methods, is described from Xizang, China. This is currently the westernmost species in Asia of the primitive lepidopteran family Neopseustidae. The new species is externally reminiscent of N.moxiensis Chen & Owada, 2009; however, it can be easily distinguished from the latter by comparison of the male genitalia and is further distinguished by the large genetic distance in DNA barcodes (COI). The adult and genitalia of the new and similar species have been illustrated. Utilizing our new data, a new classification of the genus is provided, with its members subdivided into four species groups: the meyricki-group, the moxiensis-group, the bicornuta-group, and the chentangensis-group, which are supported by both molecular and morphological evidence. A checklist of the genus and a key to the species groups are also provided.

Identifying two moth species (Lepidoptera: Ditrysia) from Saudi Arabia using mitochondrial 16S rRNA sequences

The mitochondrial genetic markers are considered useful tools for discrimination between more closely related lepidopteran taxa. Therefore, the present study aimed to investigate the role of mitochondrial (mt) 16 s rRNA gene in the determination of the taxonomic position for two moth species within Ditrysia clade. Maximum likelihood analysis has indicated a well-supported dendrogram based on the Tamura-Nei model for the recovered lepidopterans. The mt 16 s rRNA query sequences from 24 species within seven families were analyzed. This analysis and bootstrap confidence revealed two major clades representing Glossata suborder within Lepidoptera, with a close relationship of Noctuoidea + (Pyraloidea (Hesperioidea + Papilionoidea)). The subfamily Heliothinae forming a sister group with Risobinae (Noctinae + Hadeninae). In addition, there is a clear observation about the close relation between Phycitinae + Galleriinae within Pyraloidea and Cyrestinae + Limenitidinae within Papilionoidea. The present study supported that the Helicoverpa and Meroptera species are the first accounts of these genera inhabiting Saudi Arabia.

Phylogeny of gracillariid leaf-mining moths: evolution of larval behaviour inferred from phylogenomic and Sanger data

Gracillariidae is the most taxonomically diverse cosmopolitan leaf-mining moth family, consisting of nearly 2000 named species in 105 described genera, classified into eight extant subfamilies. The majority of gracillariid species are internal plant feeders as larvae, creating mines and galls in plant tissue. Despite their diversity and ecological adaptations, their phylogenetic relationships, especially among subfamilies, remain uncertain. Genomic data (83 taxa, 589 loci) were integrated with Sanger data (130 taxa, 22 loci), to reconstruct a phylogeny of Gracillariidae. Based on analyses of both datasets combined and analyzed separately, monophyly of Gracillariidae and all its subfamilies, monophyly of the clade "LAMPO" (subfamilies: Lithocolletinae, Acrocercopinae, Marmarinae, Phyllocnistinae, and Oecophyllembiinae) and relationships of its subclade "AMO" (subfamilies: Acrocercopinae, Marmarinae, and Oecophyllembiinae) were strongly supported. A sister-group relationship of Ornixolinae to the remainder of the family, and a monophyletic leaf roller lineage (Callicercops Vári + Parornichinae) + Gracillariinae, as sister to the "LAMPO" clade were supported by the most likely tree. Dating analyses indicate a mid-Cretaceous (105.3 Ma) origin of the family, followed by a rapid diversification into the nine subfamilies predating the Cretaceous-Palaeogene extinction. We hypothesize that advanced larval behaviours, such as making keeled or tentiform blotch mines, rolling leaves and galling, allowed gracillariids to better avoid larval parasitoids allowing them to further diversify. Finally, we stabilize the classification by formally re-establishing the subfamily ranks of Marmarinae stat.rev., Oecophyllembiinae stat.rev. and Parornichinae stat.rev., and erect a new subfamily, Callicercopinae Li, Ohshima and Kawahara to accommodate the enigmatic genus Callicercops.

Morphological and behavioural differences facilitate tropical butterfly persistence in variable environments

The thermal biology of ectotherms largely determines their abundance and distributions. In general, tropical species inhabiting warm and stable thermal environments tend to have low tolerance to cold and variable environments, which may restrict their expansion into temperate climates. However, the distribution of some tropical species does extend into cooler areas such as tropical borders and high elevation tropical mountains. Behavioural and morphological differences may therefore play important roles in facilitating tropical species to cope with cold and variable climates at tropical edges. We used field-validated biophysical models to estimate body temperatures of butterflies across elevational gradients at three sites in southern China and assessed the contribution of behavioural and morphological differences in facilitating their persistence in tropical and temperate climates. We investigated the effects of temperature on the activity of 4,844 individuals of 144 butterfly species along thermal gradients and tested whether species of different climatic affinities-tropical and widespread (distributed in both temperate and tropical regions)-differed in their thermoregulatory strategies (i.e. basking). In addition, we tested whether thermally related morphology or the strength of solar radiation (when butterflies were recorded) was related to such differences. We found that activities of tropical species were restricted (low abundance) at low air temperatures compared to widespread species. Active tropical species were also more likely to bask at cooler body temperatures than widespread species. Heat gain from behavioural thermoregulation was higher for tropical species (when accounting for species abundance), and heat gain correlated with larger thorax widths but not with measured solar radiation. Our results indicate that physiological intolerance to cold temperatures in tropical species may be compensated through behavioural and morphological responses in thermoregulation in variable subtropical environments. Increasing climatic variability with climate change may render tropical species more vulnerable to cold weather extremes compared to widespread species that are more physiologically suited to variable environments.

Mitochondrial Genomes of Hestina persimilis and Hestinalis nama (Lepidoptera, Nymphalidae): Genome Description and Phylogenetic Implications

Simple Summary

In this study, the mitogenomes of Hestina persimilis and Hestinalis nama were obtained via sanger sequencing. Compared with other mitogenomes of Apaturinae butterflies, conclusions can be made that the mitogenomes of Hestina persimilis and Hestinalis nama are highly conservative. The phylogenetic trees build upon mitogenomic data showing that the relationships among Nymphalidae are similar to previous studies. Hestinalisnama is apart from Hestina, and closely related to Apatura, forming a monophyletic clade.

Abstract

In this study, the complete mitochondrial genomes (mitogenomes) of Hestina persimilis and Hestinalis nama (Nymphalidae: Apaturinae) were acquired. The mitogenomes of H. persimilis and H. nama are 15,252 bp and 15,208 bp in length, respectively. These two mitogenomes have the typical composition, including 37 genes and a control region. The start codons of the protein-coding genes (PCGs) in the two mitogenomes are the typical codon pattern ATN, except CGA in the cox1 gene. Twenty-one tRNA genes show a typical clover leaf structure, however, trnS1(AGN) lacks the dihydrouridine (DHU) stem. The secondary structures of rrnL and rrnS of two species were predicted, and there are several new stem loops near the 5′ of rrnL secondary structure. Based on comparative genomic analysis, four similar conservative structures can be found in the control regions of these two mitogenomes. The phylogenetic analyses were performed on mitogenomes of Nymphalidae. The phylogenetic trees show that the relationships among Nymphalidae are generally identical to previous studies, as follows: Libytheinae\Danainae + ((Calinaginae + Satyrinae) + Danainae\Libytheinae + ((Heliconiinae + Limenitidinae) + (Nymphalinae + (Apaturinae + Biblidinae)))). Hestinalisnama is apart from Hestina, and closely related to Apatura, forming monophyly.

Molecular phylogenetic analysis of the Paguristes tortugae Schmitt, 1933 complex and selected other Paguroidea (Crustacea: Decapoda: Anomura)

Morphological characters, as presently applied to describe members of the Paguristes tortugae Schmitt, 1933 species complex, appear to be of limited value in inferring phylogenetic relationships within the genus, and may have similarly misinformed understanding of relationships between members of this complex and those presently assigned to the related genera Areopaguristes Rahayu McLaughlin, 2010 and Pseudopaguristes McLaughlin, 2002. Previously undocumented observations of similarities and differences in color patterns among populations additionally suggest genetic divergences within some species, or alternatively seem to support phylogenetic groupings of some species. In the present study, a Maximum Likelihood (ML) phylogenetic analysis was undertaken based on the H3, 12S mtDNA, and 16S mtDNA sequences of 148 individuals, primarily representatives of paguroid species from the western Atlantic. This molecular analysis supported a polyphyletic Diogenidae Ortmann, 1892, although incomplete taxonomic sampling among the genera of Diogenidae limits the utility of this finding for resolving family level relationships. Several hypotheses regarding the evolutionary relationships among hermit crab genera were refuted by the Kishino-Hasegawa (KH). Shimodaira-Hasegawa (SH) and Approximately Unbiased (AU) tree topology tests, among them the hypothesis that Areopaguristes is monophyletic. A lack of support for the monophyly of Areopaguristes calls into question the phylogenetic validity of gill number for the differentiation of Paguristes, Areopaguristes, and Pseudopaguristes. The study was inconclusive with regard to the relationships among these three genera, but previously unknown diversity within both Paguristes and Areopaguristes was demonstrated. Existence of an undescribed species confounded under the name Paguristes tortugae Schmitt, 1933 was supported by genetics, morphology, and coloration. A second undescribed species with remarkable similarity to Areopaguristes hummi Wass, 1955 was discovered based on genetics and coloration.

Convergent Evolution of Wingbeat-Powered Anti-Bat Ultrasound in the Microlepidoptera

Bats and moths provide a textbook example of predator-prey evolutionary arms races, demonstrating adaptations, and counter adaptations on both sides. The evolutionary responses of moths to the biosonar-led hunting strategies of insectivorous bats include convergently evolved hearing structures tuned to detect bat echolocation frequencies. These allow many moths to detect hunting bats and manoeuvre to safety, or in the case of some taxa, respond by emitting sounds which startle bats, jam their biosonar, and/or warn them of distastefulness. Until now, research has focused on the larger macrolepidoptera, but the recent discovery of wingbeat-powered anti-bat sounds in a genus of deaf microlepidoptera (Yponomeuta), suggests that the speciose but understudied microlepidoptera possess further and more widespread anti-bat defences. Here we demonstrate that wingbeat-powered ultrasound production, likely providing an anti-bat function, appears to indeed be spread widely in the microlepidoptera; showing that acoustically active structures (aeroelastic tymbals, ATs) have evolved in at least three, and likely four different regions of the wing. Two of these tymbals are found in multiple microlepidopteran superfamilies, and remarkably, three were found in a single subfamily. We document and characterise sound production from four microlepidopteran taxa previously considered silent. Our findings demonstrate that the microlepidoptera contribute their own unwritten chapters to the textbook bat-moth coevolutionary arms race.

Complete mitochondrial genome of Acanthopsyche nigraplaga (Lepidoptera: Psychidae)

We report the mitochondrial genome (mitogenome) of a case-making moth Acanthopsyche nigraplaga Wileman, 1911 (Lepidoptera: Psychidae). The 15,704 bp long complete mitogenome comprises a typical set of genes [13 protein-coding genes (PCGs), 2 rRNA genes, and 22 tRNA genes] and one major non-coding, A + T-rich region, with an arrangement identical to that observed in most lepidopteran mitogenomes. Twelve of the 13 PCGs of the A. nigraplaga mitogenome initiate with a typical ATN start codon, however COI contains the atypical CGA start codon that is common for lepidopteran COI genes. A phylogenetic analysis using concatenated nucleotide sequences of the 13 PCGs and 2 rRNA genes using the Bayesian inference method fully resolved A. nigraplaga in a monophyletic clade within the Psychidae. Acanthopsyche nigraplaga was situated in a sister position to Eumeta variegata and Mahasena oolona with high nodal support. As more mitogenome sequences are available further scrutinized analysis for the superfamily Tineoidea including Psychidae will be possible.

Unique Members of the Adipokinetic Hormone Family in Butterflies and Moths (Insecta, Lepidoptera)

Lepidoptera is amongst one of the four most speciose insect orders and ecologically very successful because of their ability to fly. Insect flight is always aerobic and exacts a high metabolic demand on the animal. A family of structurally related neuropeptides, generically referred to as adipokinetic hormones (AKHs), play a key role in triggering the release of readily utilizable fuel metabolites into the hemolymph from the storage forms in the fat body. We used mass spectrometry to elucidate AKH sequences from 34 species of Lepidoptera and searched the literature and publicly available databases to compile (in a phylogenetic context) a comprehensive list of all Lepidoptera sequences published/predicted from a total of 76 species. We then used the resulting set of 15 biochemically characterized AKHs in a physiological assay that measures lipid or carbohydrate mobilization in three different lepidopteran species to learn about the functional cross-activity (receptor-ligand interactions) amongst the different butterfly/moth families. Our results include novel peptide structures, demonstrate structural diversity, phylogenetic trends in peptide distribution and order-specificity of Lepidoptera AKHs. There is almost an equal occurrence of octa-, nona-, and decapeptides, with an unparalleled emphasis on nonapeptides than in any insect order. Primitive species make Peram-CAH-II, an octapeptide found also in other orders; the lepidopteran signature peptide is Manse-AKH. Not all of the 15 tested AKHs are active in Pieris brassicae; this provides insight into structure-activity specificity and could be useful for further investigations into possible biorational insecticide development.

Complete mitochondrial genome of Pterodecta felderi (Lepidoptera: Callidulidae)

We report the mitochondrial genome (mitogenome) of Pterodecta felderi (Callidulidae: Lepidoptera), which is the first mitogenome sequences in the family Callidulidae, a monotypic family in the superfamily Calliduloidea. The 15,340-bp long complete mitogenome consists of a typical set of genes (13 protein-coding genes [PCGs], 2 rRNA genes, and 22 tRNA genes) and 1 major non-coding A + T-rich region, which are arranged in a way that is frequently observed in Lepidoptera. Of the 13 PCGs, 12 P. felderi start with ATN, except for COI, which starts with CGA. The P. felderi mitogenome consists of 210-bp long intergenic-spacer sequences and 27-bp long overlaps. Phylogenetic analysis of superfamilial relationships in the lepidopteran clade Obtectomera with concatenated sequences of the 13 PCGs and 2 rRNA genes using the Bayesian inference method showed that Calliduloidea, which is only represented by P. felderi, was placed as the most basal lineage about Macroheterocera (Lasiocampoidea, Bombycoidea, Mimallonoidea, Noctuoidea, and Drepanoidea), Papilionoidea, and Pyraloidea.

The Evolution of Endophagy in Herbivorous Insects

Herbivorous feeding inside plant tissues, or endophagy, is a common lifestyle across Insecta, and occurs in insect taxa that bore, roll, tie, mine, gall, or otherwise modify plant tissues so that the tissues surround the insects while they are feeding. Some researchers have developed hypotheses to explain the adaptive significance of certain endophytic lifestyles (e.g., miners or gallers), but we are unaware of previous efforts to broadly characterize the adaptive significance of endophagy more generally. To fill this knowledge gap, we characterized the limited set of evolutionary selection pressures that could have encouraged phytophagous insects to feed inside plants, and then consider how these factors align with evidence for endophagy in the evolutionary history of orders of herbivorous insects. Reviewing the occurrence of endophytic taxa of various feeding guilds reveals that the pattern of evolution of endophagy varies strongly among insect orders, in some cases being an ancestral trait (e.g., Coleoptera and Lepidoptera) while being more derived in others (e.g., Diptera). Despite the large diversity of endophagous lifestyles and evolutionary trajectories that have led to endophagy in insects, our consideration of selection pressures leads us to hypothesize that nutritionally based factors may have had a stronger influence on evolution of endophagy than other factors, but that competition, water conservation, and natural enemies may have played significant roles in the development of endophagy.

Ecology and evolution of cycad-feeding Lepidoptera

Cycads are an ancient group of tropical gymnosperms that are toxic to most animals - including humans - though the larvae of many moths and butterflies (order: Lepidoptera) feed on cycads with apparent immunity. These insects belong to distinct lineages with varying degrees of specialisation and diverse feeding ecologies, presenting numerous opportunities for comparative studies of chemically mediated eco-evolutionary dynamics. This review presents the first evolutionary evaluation of cycad-feeding among Lepidoptera along with a comprehensive review of their ecology. Our analysis suggests that multiple lineages have independently colonised cycads from angiosperm hosts, yet only a few clades appear to have radiated following their transitions to cycads. Defensive traits are likely important for diversification, as many cycad specialists are warningly coloured and sequester cycad toxins. The butterfly family Lycaenidae appears to be particularly predisposed to cycad-feeding and several cycadivorous lycaenids are warningly coloured and chemically defended. Cycad-herbivore interactions provide a promising but underutilised study system for investigating plant-insect coevolution, convergent and divergent adaptations, and the multi-trophic significance of defensive traits; therefore the review ends by suggesting specific research gaps that would be fruitfully addressed in Lepidoptera and other cycad-feeding insects.

Ecological Drift and Directional Community Change in an Isolated Mediterranean Forest Reserve-Larger Moth Species Under Higher Threat

Long-term data are important to understand the changes in ecological communities over time but are quite rare for insects. We analyzed such changes using historic museum collections. For our study area, an isolated forest reserve in North-East Italy, data from the past 80 yr were available. We used records of 300 moth species to analyze whether extinction risk was linked to their body size or to their degree of ecological specialization. Specialization was scored 1) by classifying larval food affiliations, habitat preferences, and the northern distributional limit and 2) by analyzing functional dispersion (FDis) within species assemblages over time. Our results show that locally extinct species (mean wingspan: 37.0 mm) were larger than persistent (33.2 mm) or previously unrecorded ones (30.7 mm), leading to a smaller mean wingspan of the moth community over time. Some ecological filters appear to have selected against bigger species. By using coarse specialization categories, we did not observe any relationship with local extinction risk. However, FDis, calculated across 12 species traits, significantly decreased over time. We conclude that simple classification systems might fail in reflecting changes in community-wide specialization. Multivariate approaches such as FDis may provide deeper insight, as they reflect a variety of ecological niche dimensions. With the abandonment of extensive land use practices, natural succession seems to have shifted the moth community toward a preponderance of forest-affiliated species, leading to decreased FDis values. Multivariate analyses of species composition also confirmed that the moth community has significantly changed during the last 80 yr.

Identification and characterization of the vasa gene in the diamondback moth, Plutella xylostella

Vasa is an ATP-dependent RNA helicase, participating in multiple biological processes. It has been widely used as a germ cell marker and its promoter has become a key component of several genetic pest control systems. Here we present the vasa gene structure and its promoter activity in Plutella xylostella, one of the most destructive pests of cruciferous crops. Full length Pxvasa cDNA sequences were obtained, revealing 14 exons and at least 30 alternatively spliced transcripts. Inferred amino acid sequences showed nine conserved DEAD-box family protein motifs with partial exclusion from some isoforms. Real-time quantitative PCR indicated the up-regulation of Pxvasa in both female and male adults compared with other developmental stages, and the expression levels of Pxvasa were found to be much higher in adult gonads, especially ovaries, than in other tissues. The putative promoter region of Pxvasa was sequenced and several ecdysone-induced transcription factor (TF) binding sites were predicted in silico. To further analyze the promoter region, two upstream regulatory fragments of different lengths were tested as putative promoters in transient cell and embryo expression assays, one of which was subsequently utilized to drive Cas9 expression in vivo. A transgenic line was recovered and the expression patterns of Cas9 and native Pxvasa were profiled in adult tissues and eggs with RT-PCR. This work provides the foundation for further studies on the gene functions of Pxvasa as well as the potential application of its promoter in genetic manipulation of P. xylostella.

Evolutionary constraints on flicker fusion frequency in Lepidoptera

Flying insects occupy both diurnal and nocturnal niches, and their visual systems encounter distinct challenges in both conditions. Visual adaptations, such as superposition eyes of moths, enhance sensitivity to low light levels but trade off with spatial and temporal resolution. Conversely, apposition eyes of butterflies enable high spatial resolution but are poorly sensitive in dim light. Although diel activity patterns of insects influence visual processing, their role in evolution of visual systems is relatively unexplored. Lepidopteran insects present an excellent system to study how diel activity patterns and phylogenetic position influence the visual transduction system. We addressed this question by comparing electroretinography measurements of temporal response profiles of diverse Lepidoptera to light stimuli that were flickering at different frequencies. Our data show that the eyes of diurnal butterflies are sensitive to visual stimuli of higher temporal frequencies than nocturnal moths. Hesperiid skippers, which are typically diurnal or crepuscular, exhibit intermediate phenotypes with peak sensitivity across broader frequency range. Across all groups, species within families exhibited similar phenotypes irrespective of diel activity. Thus, Lepidopteran photoreceptors may have diversified under phylogenetic constraints, and shifts in their sensitivity to higher temporal frequencies occurred concomitantly with the evolution of diurnal lifestyles.

Complete mitochondrial genome of Amorophaga japonica Robinson, 1986 (Lepidoptera: Tineidae)

The complete mitochondrial genome (mitogenome) of Amorophaga japonica Robinson, 1986 (Lepidoptera: Tineidae), comprises 15,027 base pairs (bp) and contains a typical set of genes (13 protein-coding genes [PCGs], 2 rRNA genes, and 22 tRNA genes), and 1 non-coding region. The genome has an arrangement, trnW-trnY-trnC, instead of typical trnW- trnC-trnY at the ND2 and COI junction. This arrangement is unique in lepidopteran mitogenomes. Unlike most lepidopteran insects, which have CGA as the start codon for the COI gene sequence, A. japonica COI had a typical ATT codon. The A + T-rich region was unusually short, with only 199 bp. Phylogenetic analyses with concatenated sequences of the 13 PCGs and two rRNA genes using the Bayesian inference method placed A. japonica in Tineidae as a sister to the cofamilial species, Tineola bisselliella, with high nodal support (Bayesian posterior probability [BPP] = 0.99), presenting the superfamily Tineoidea in a monophyletic group with a BPP of 0.99. Gracillarioidea, represented by three species of Gracillariidae, formed a monophyletic group with the highest BPP, but the Leucoptera malifoliella in Yponomeutoidea was unusually grouped together with the Gracillarioidea with the highest nodal support. As more mitogenome sequences are available, further analysis to infer the relationships among superfamilies of Lepidoptera might be possible.

A multilocus analysis of Epicopeiidae (Lepidoptera, Geometroidea) provides new insights into their relationships and the evolutionary history of mimicry

The family Epicopeiidae is a small group of day-flying moths, known for mimicking many different groups of butterflies and moths. So far, there still lacks a reliable phylogenetic framework of Epicopeiidae that is necessary to our understanding of the evolutionary process of their mimicry. In this study, we sequenced 94 nuclear protein-coding markers for 56 epicopeiid samples and 11 outgroups, covering all ten genera of Epicopeiidae. We used homemade PCR-generated baits to capture target sequences, which allowed us to utilize old and dried specimens that were difficult to handle by conventional PCR + Sanger sequencing. Maximum likelihood and Bayesian analyses of the newly obtained dataset (86,388 bp) at both DNA and protein levels produced identical phylogenies with strong support. The non-mimicry genus Deuveia is the sister group of other epicopeiid genera. Epicopeia and Nossa are not monophyletic, and these two genera nest together to form a clade. We also estimated divergence times of Epicopeiidae and found that their initial diversification happened in Eocene about 41 million years ago. The ancestral state reconstruction of mimicry type for this family suggested that thelast common ancestor of epicopeiid moths is non-mimetic, and the Riodinidae-mimicry type evolved first. In summary, our work provides a comprehensive and robust time-calibrated phylogeny of Epicopeiidae that provides a sound framework for revising their classification and interpreting character evolution.

Distal leg structures of the Aculeata (Hymenoptera): A comparative evolutionary study of Sceliphron (Sphecidae) and Formica (Formicidae)

The distal parts of the legs of Sceliphron caementarium (Sphecidae) and Formica rufa (Formicidae) are documented and discussed with respect to phylogenetic and functional aspects. The prolegs of Hymenoptera offer an array of evolutionary novelties, mainly linked with two functional syndromes, walking efficiently on different substrates and cleaning the body surface. The protibial-probasitarsomeral cleaning device is almost always well-developed. A complex evolutionary innovation is a triple set of tarsal and pretarsal attachment devices, including tarsal plantulae, probasitarsomeral spatulate setae, and an arolium with an internal spring-like arcus, a dorsal manubrium, and a ventral planta. The probasitarsal adhesive sole and a complex arolium are almost always preserved, whereas the plantulae are often missing. Sceliphron has retained most hymenopteran ground plan features of the legs, and also Formica, even though the adhesive apparatus of Formicidae shows some modifications, likely linked to ground-oriented habits of most ants. Plantulae are always absent in extant ants, and the arolium is often reduced in size, and sometimes vestigial. The arolium contains resilin in both examined species. Additionally, resilin enriched regions are also present in the antenna cleaners of both species, although they differ in which of the involved structures is more flexible, the calcar in Sceliphron and the basitarsal comb in Formica. Functionally, the hymenopteran distal leg combines (a) interlocking mechanisms (claws, spine-like setae) and (b) adhesion mechanisms (plantulae, arolium). On rough substrate, claws and spine-like setae interlock with asperities and secure a firm grip, whereas the unfolding arolium generates adhesive contact on smooth surfaces. Differences of the folded arolium of Sceliphron and Formica probably correlate with differences in the mechanism of folding/unfolding.

Disparate expression specificities coded by a shared Hox-C enhancer

Can a single regulatory sequence be shared by two genes undergoing functional divergence? Here we describe a single promiscuous enhancer within the Drosophila Antennapedia Complex, EO053, that directs aspects of the expression of two adjacent genes, pb (a Hox2 ortholog) and zen2 (a divergent Hox3 paralog), with disparate spatial and temporal expression patterns. We were unable to separate the pb-like and zen2-like specificities within EO053, and we identify sequences affecting both expression patterns. Importantly, genomic deletion experiments demonstrate that EO053 cooperates with additional pb- and zen2-specific enhancers to regulate the mRNA expression of both genes. We examine sequence conservation of EO053 within the Schizophora, and show that patterns of synteny between the Hox2 and Hox3 orthologs in Arthropods are consistent with a shared regulatory relationship extending prior to the Hox3/zen divergence. Thus, EO053 represents an example of two genes having evolved disparate outputs while utilizing this shared regulatory region.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Morphological, chemical and electrophysiological investigations of Telchin licus (Lepidoptera: Castniidae)

The giant sugarcane borer Telchin licus (Drury, 1773) (Lepidoptera: Castniidae) is a day-flying moth pest of sugarcane, pineapples and bananas. To better understand the chemical communication in this species, we examined the morphology of its olfactory system and the chemical composition of its body parts. The ventral surface of the clubbed antennae of T. licus has six morphological types of sensilla: sensilla trichodea, basiconica, chaetica, squamiforma, coeloconica, and auricillica. The telescopic ovipositor shows no evidence of a sexual gland, or female-specific compounds. On the other hand, the midleg basitarsus of males releases (E,Z)-2,13-octadecadienol and (Z,E)-2,13-octadecadienol, which are electroantennographically active in both sexes. These compounds are known female sex pheromones in the Sesiidae family and are male-specific compounds in another castniid moth, although further investigations are necessary to elucidate their ecological role in the Castniidae family.

Regional assemblages shaped by historical and contemporary factors: Evidence from a species-rich insect group

Understanding diversity patterns requires accounting for the roles of both historical and contemporary factors in the assembly of communities. Here, we compared diversity patterns of two moth assemblages sampled from Taihang and Yanshan mountains in Northern China and performed ancestral range reconstructions using the Multi-State Speciation and Extinction model, to track the origins of these patterns. Further, we estimated diversification rates of the two moth assemblages and explored the effects of contemporary ecological factors. From 7,788 specimens we identified 835 species belonging to 23 families, using both DNA barcode analysis and morphology. Moths in Yanshan mountains showed higher species diversity than in Taihang mountains. Ancestral range analysis indicated Yanshan as the origin, with significant historical dispersals from Yanshan to Taihang. Asymmetrical diversification, population expansion, along with frequent and considerable gene flow were detected between communities. Moreover, dispersal limitation or the joint effect of environment filtering and dispersal limitation were inferred as main driving forces shaping current diversity patterns. In summary, we demonstrate that a multiscale (community, population and species level) analysis incorporating both historical and contemporary factors can be useful in delineating factors contributing to community assembly and patterning in diversity.

Symmetry systems on the wings of Dichromodes Guenée (Lepidoptera: Geometridae) are unconstrained by venation

The nymphalid groundplan, an idealized schematic illustrating the essential elements of butterfly wing patterns, predicts a consistent relationship between color pattern and wing venation. Moths in the family Geometridae have wing shapes and patterns that often resemble those of butterflies, and until recently, this family was believed to be among butterflies’ closest relatives. However, an examination of the geometrid genus Dichromodes Guenée, 1858 shows no consistent relationship between the central symmetry system and wing venation. Whereas the distal edge of the central symmetry system is predicted to reach the costal margin proximal to the Subcostal vein in butterflies and acronictine moths, it has no consistent relationship with the Subcostal, Radius, or Radial Sector 1 veins in Dichromodes. This finding highlights developmental diversity that was previously overlooked due to the overwhelming preference for butterflies in studies of lepidopteran wing patterns.

Absence of W Chromosome in Psychidae Moths and Implications for the Theory of Sex Chromosome Evolution in Lepidoptera

Moths and butterflies (Lepidoptera) are the largest group with heterogametic females. Although the ancestral sex chromosome system is probably Z0/ZZ, most lepidopteran species have the W chromosome. When and how the W chromosome arose remains elusive. Existing hypotheses place the W origin either at the common ancestor of Ditrysia and Tischeriidae, or prefer independent origins of W chromosomes in these two groups. Due to their phylogenetic position at the base of Ditrysia, bagworms (Psychidae) play an important role in investigating the W chromosome origin. Therefore, we examined the W chromosome status in three Psychidae species, namely Proutia betulina, Taleporia tubulosa, and Diplodoma laichartingella, using both classical and molecular cytogenetic methods such as sex chromatin assay, comparative genomic hybridization (CGH), and male vs. female genome size comparison by flow cytometry. In females of all three species, no sex chromatin was found, no female-specific chromosome regions were revealed by CGH, and a Z-chromosome univalent was observed in pachytene oocytes. In addition, the genome size of females was significantly smaller than males. Overall, our study provides strong evidence for the absence of the W chromosome in Psychidae, thus supporting the hypothesis of two independent W chromosome origins in Tischeriidae and in advanced Ditrysia.

Phylogenomics reveals the evolutionary timing and pattern of butterflies and moths

Lepidoptera play key roles in many biological systems. Butterflies are hypothesized to have evolved contemporaneously with flowering plants, and moths are thought to have gained anti-bat defenses in response to echolocating predatory bats, but these hypotheses have largely gone untested. Using a transcriptomic, dated evolutionary tree of Lepidoptera, we demonstrate that the most recent common ancestor of Lepidoptera is considerably older than previously hypothesized. The oldest moths in crown Lepidoptera were present in the Carboniferous, some 300 million years ago, and began to diversify largely in synchrony with angiosperms. We show that multiple lineages of moths independently evolved hearing organs well before the origin of bats, rejecting the hypothesis that lepidopteran hearing organs arose in response to these predators.

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

First mitochondrial genome from Yponomeutidae (Lepidoptera, Yponomeutoidea) and the phylogenetic analysis for Lepidoptera

The complete mitochondrial genome (mitogenome) of Yponomeuta montanatus is sequenced and compared with other published yponomeutoid mitogenomes. The mitogenome is circular, 15,349 bp long, and includes the typical metazoan mitochondrial genes (13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes) and an A + T-rich region. All 13 protein-coding genes use a typical start codon ATN, the one exception being cox1, which uses CGA across yponomeutoid mitogenomes. Comparative analyses further show that the secondary structures of tRNAs are conserved, including loss of the Dihydorouidine (DHU) arm in trnS1 (AGN), but remarkable nucleotide variation has occurred mainly in the DHU arms and pseudouridine (TψC) loops. A + T-rich regions exhibit substantial length variation among yponomeutoid mitogenomes, and conserved sequence blocks are recognized but some of them are not present in all species. Multiple phylogenetic analyses confirm the position of Y. montanatus in Yponomeutoidea. However, the superfamily-level relationships in the Macroheterocera clade in Lepidoptera recovered herein show considerable difference with that recovered in previous mitogenomic studies, raising the necessity of extensive phylogenetic investigation when more mitogenomes become available for this clade.

Phylogenomics resolves major relationships and reveals significant diversification rate shifts in the evolution of silk moths and relatives

BACKGROUND

Silkmoths and their relatives constitute the ecologically and taxonomically diverse superfamily Bombycoidea, which includes some of the most charismatic species of Lepidoptera. Despite displaying spectacular forms and diverse ecological traits, relatively little attention has been given to understanding their evolution and drivers of their diversity. To begin to address this problem, we created a new Bombycoidea-specific Anchored Hybrid Enrichment (AHE) probe set and sampled up to 571 loci for 117 taxa across all major lineages of the Bombycoidea, with a newly developed DNA extraction protocol that allows Lepidoptera specimens to be readily sequenced from pinned natural history collections.

RESULTS

The well-supported tree was overall consistent with prior morphological and molecular studies, although some taxa were misplaced. The bombycid Arotros Schaus was formally transferred to Apatelodidae. We identified important evolutionary patterns (e.g., morphology, biogeography, and differences in speciation and extinction), and our analysis of diversification rates highlights the stark increases that exist within the Sphingidae (hawkmoths) and Saturniidae (wild silkmoths).

CONCLUSIONS

Our study establishes a backbone for future evolutionary, comparative, and taxonomic studies of Bombycoidea. We postulate that the rate shifts identified are due to the well-documented bat-moth "arms race". Our research highlights the flexibility of AHE to generate genomic data from a wide range of museum specimens, both age and preservation method, and will allow researchers to tap into the wealth of biological data residing in natural history collections around the globe.