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

Evolution and losses of spines in slug caterpillars (Lepidoptera: Limacodidae)

Larvae of the cosmopolitan family Limacodidae, commonly known as "slug" caterpillars, are well known because of the widespread occurrence of spines with urticating properties, a morpho-chemical adaptive trait that has been demonstrated to protect the larvae from natural enemies. However, while most species are armed with rows of spines ("nettle" caterpillars), slug caterpillars are morphologically diverse with some species lacking spines and thus are nonstinging. It has been demonstrated that the evolution of spines in slug caterpillars may have a single origin and that this trait is possibly derived from nonstinging slug caterpillars, but these conclusions were based on limited sampling of mainly New World taxa; thus, the evolution of spines and other traits within the family remains unresolved. Here, we analyze morphological variation in slug caterpillars within an evolutionary framework to determine character evolution of spines with samples from Asia, Australia, North America, and South America. The phylogeny of the Limacodidae was reconstructed based on a multigene dataset comprising five molecular markers (5.6 Kbp: COI, 28S, 18S, EF-1α, and wingless) representing 45 species from 40 genera and eight outgroups. Based on this phylogeny, we infer that limacodids evolved from a common ancestor in which the larval type possessed spines, and then slug caterpillars without spines evolved independently multiple times in different continents. While larvae with spines are well adapted to avoiding generalist predators, our results imply that larvae without spines may be suited to different ecological niches. Systematic relationships of our dataset indicate six major lineages, several of which have not previously been identified.

A comprehensive molecular phylogeny of Geometridae (Lepidoptera) with a focus on enigmatic small subfamilies

Our study aims to investigate the relationships of the major lineages within the moth family Geometridae, with a focus on the poorly studied Oenochrominae-Desmobathrinae complex, and to translate some of the results into a coherent subfamilial and tribal level classification for the family. We analyzed a molecular dataset of 1,206 Geometroidea terminal taxa from all biogeographical regions comprising up to 11 molecular markers that includes one mitochondrial (COI) and 10 protein-coding nuclear gene regions (wingless, ArgK, MDH, RpS5, GAPDH, IDH, Ca-ATPase, Nex9, EF-1alpha, CAD). The molecular data set was analyzed using maximum likelihood as implemented in IQ-TREE and RAxML. We found high support for the subfamilies Larentiinae, Geometrinae and Ennominae in their traditional scopes. Sterrhinae becomes monophyletic only if Ergavia Walker, Ametris Hübner and Macrotes Westwood, which are currently placed in Oenochrominae, are formally transferred to Sterrhinae. Desmobathrinae and Oenochrominae are found to be polyphyletic. The concepts of Oenochrominae and Desmobathrinae required major revision and, after appropriate rearrangements, these groups also form monophyletic subfamily-level entities. Oenochrominae s.str. as originally conceived by Guenée is phylogenetically distant from Epidesmia and its close relatives. The latter is hereby described as the subfamily Epidesmiinae Murillo-Ramos, Sihvonen & Brehm, subfam. nov. Epidesmiinae are a lineage of "slender-bodied Oenochrominae" that include the genera Ecphyas Turner, Systatica Turner, Adeixis Warren, Dichromodes Guenée, Phrixocomes Turner, Abraxaphantes Warren, Epidesmia Duncan & Westwood and Phrataria Walker. Archiearinae are monophyletic when Dirce and Acalyphes are formally transferred to Ennominae. We also found that many tribes were para- or polyphyletic and therefore propose tens of taxonomic changes at the tribe and subfamily levels. Archaeobalbini stat. rev. Viidalepp (Geometrinae) is raised from synonymy with Pseudoterpnini Warren to tribal rank. Chlorodontoperini Murillo-Ramos, Sihvonen & Brehm, trib. nov. and Drepanogynini Murillo-Ramos, Sihvonen & Brehm, trib. nov. are described as new tribes in Geometrinae and Ennominae, respectively.

The wing pattern of Hydriomena Hübner, [1825] (Lepidoptera: Geometridae: Larentiinae) lacks a predictable relationship with venation

Two simple models have been successfully applied to predict the relationship between wing pattern and venation in various lineages of Lepidoptera. However, neither of these models holds for the geometrid genus Hydriomena Hübner, 1825. Wing patterns in Hydriomena were studied intensively during the 1920s after the description of the nymphalid groundplan, an idealized schematic that outlines the primary elements of butterfly wing patterns; geometrids strongly resemble butterflies and, until recently, were considered to be among their closest relatives. The evolution of wing pattern in Geometridae has been neglected since the 1930s. Here, the relationship between wing pattern and venation is examined for Hydriomena costipunctata Barnes and McDunnough, 1912 and the Hydriomena speciosata (Packard, 1873) group. These two lineages have some of the simplest wing patterns in Hydriomena, consisting of large, well-defined dark and light pattern elements. The relationship between wing pattern and venation varies considerably within and between these lineages and can even vary between the right and left wings on the same individual. Although many different wing patterns were observed among the individuals examined for this study, not one can be reconciled with either of the models that successfully predict the relationship between wing pattern and venation in many other groups of Lepidoptera. This suggests that bands occurring on the wings of Hydriomena are not homologous with those on the wings of butterflies or of Acronictinae (Macroheterocera: Noctuidae), the only other two obtectomeran lineages for which the relationship between wing pattern and venation has been examined in recent years.

Optical Modelling and Phylogenetic Analysis Provide Clues to the Likely Function of Corneal Nipple Arrays in Butterflies and Moths

The lenses in compound eyes of butterflies and moths contain an array of nipple-shaped protuberances, or corneal nipples. Previous work has suggested that these nipples increase light transmittance and reduce the eye glare of moths that are inactive during the day. This work builds on but goes further than earlier analyses suggesting a functional role for these structures including, for the first time, an explanation of why moths are attracted to UV light. Using a phylogenetic approach and 3D optical modelling, we show empirically that these arrays have been independently lost from different groups of moths and butterflies and vary within families. We find differences in the shape of nipples between nocturnal and diurnal species, and that anti-glow reflectance levels are different at different wave-lengths, a result thereby contradicting the currently accepted theory of eye glow for predator avoidance. We find that there is reduced reflectance, and hence greater photon absorption, at UV light, which is probably a reason why moths are attracted to UV. We note that the effective refractive index at the end of the nipples is very close to the refractive index of water, allowing almost all the species with nipples to see without distortion when the eye is partially or completely wet and providing the potential to keep eyes dry. These observations provide a functional explanation for these arrays. Of special interest is the finding that their repeated and independent loss across lepidopteran phylogeny is inconsistent with the explanation that they are being lost in the ‘higher’, more active butterflies.

Complete mitochondrial genome of two Thitarodes species (Lepidoptera, Hepialidae), the host moths of Ophiocordyceps sinensis and phylogenetic implications

Thitarodes (Lepidoptera, Hepialidae) is the only genus that hosts to the Ophiocordyceps sinensis, a traditional Chinese medicine considered as a powerful medicinal supplement. In this study, the complete mitochondrial genomes (mitogenomes) of two species, T. damxungensis and T. pui, have been sequenced, which are 15,928 bp and 15,362 bp in size respectively, and both contain 13 protein-coding genes (PCGs), 2 rRNAs, 22 tRNAs and an AT-rich region. Like other hepialoids, the gene arrangement of the mitogenomes of T. damxungensis and T. pui is identical to the ancestral arrangement but differs from those of other lepidopteran species on account of the different arrangements of trnM, trnI, and trnQ. The size of AT-rich region is 545 bp in T. damxungensis and 1030 bp in T. pui. Tandem repetition in the AT-rich region is responsible for the length difference of the A + T-rich region in both species. In Hepialidae, the phylogenetic study based on the dataset of the sequences that combined the protein-coding genes and RNA genes suggested that the species T. yunnanensis should still belong to the genus Thitarodes rather than Ahamns, which is different from the results based on the traditional phylogeny.

Mitochondrial phylogeny and comparative mitogenomics of closely related pine moth pests (Lepidoptera: Dendrolimus)

BACKGROUND

Pine moths, Dendrolimus spp. (Lasiocampidae), are serious economic pests of conifer forests. Six closely related species (Dendrolimus punctatus, D. tabulaeformis, D. spectabilis, D. superans, D. houi, and D. kikuchii) occur in China and cause serious damage to coniferophyte. The complete mito genomes of Dendrolimus genus are significant to resolve the phylogenetic relationship and provide theoretical support in pest control.

METHODS

The complete mitogenomes of three species (D. superans, D. houi, and D. kikuchii) were sequenced based on PCR-amplified with universal primers, which were used to amplify initial fragments. Phylogenetic analyses were carried out with 78 complete mitogenomes of lepidopteran species from 10 superfamilies.

RESULTS

The complete mitochondrial genomes of these three species were 15,417, 15,381, and 15,377 bp in length, separately. The phylogenetic analyses produced consistent results for six Dendrolimus species based on complete mitogenomes, two major clades were formed, one containing D. spectabilis clustered with D. punctatus + D. tabulaeformis, and D. superans as the sister group to this three-taxon clade, the other containing D. kikuchii and D. houi. Comparative analyses of the congeneric mitochondrial genomes were performed, which showed that non-coding regions were more variable than the A+T rich region. The mitochondrial nucleotide diversity was more variable when compared within than among genus, and the concatenated tRNA region was the most conserved and the nd6 genes was the most variable.

The first mitochondrial genome of the family Epicopeiidae and higher-level phylogeny of Macroheterocera (Lepidoptera: Ditrysia)

The Macroheterocera clade contains most of the Lepidoptera species. However, extensive comparative and phylogenetic analyses of this group using complete mitochondrial genomes (mitogenome) are limited particularly in the context of increasing macroheteroceran mitogenomes reported to date. In this study, complete mitogenome of the Epicopeia hainesii is determined as the first Epicopeiidae species with mitogenome available. The whole mitogenome is circular with 15,395 bp long, and is highly biased toward A + T nucleotides (80.6%) in nucleotide composition. Comparative analyses show that gene content and arrangement of macroheteroceran mitogenomes are generally conservative and are typical of Lepidoptera but exceptions exist. In newly sequenced mitogenome, the motif "ATACTAA" is putatively located at the end of gene nad1, rather than in intergenic sequences between trnS2 and nad1 genes routinely observed in Lepidoptera. Interestingly, multiple phylogenetic analyses recover the six macroheteroceran superfamilies as Mimallonoidea + (Drepanoidea + ((Bombycoidea + Lasiocampoidea) + (Noctuoidea + Geometroidea))), providing supports for a large-scale transcriptomic study rather than various mitogenome- and multiple-gene-based investigations. In addition, our analyses consistently place the Epicopeiidae as sister group with Geometroidea, firstly demonstrating that this family is closer with Geometroidea than Drepanoidea based on mitogenome data.

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

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

Mutualistic interactions reshuffle the effects of climate change on plants across the tree of life

Climatically induced local species extinctions may trigger coextinction cascades, thus driving many more species to extinction than originally predicted by species distribution models. Using seven pollination networks across Europe that include the phylogeny and life history traits of plants, we show a substantial variability across networks in climatically predicted plant extinction-and particularly the subsequent coextinction-rates, with much higher values in Mediterranean than Eurosiberian networks. While geographic location best predicts the probability of a plant species to be driven to extinction by climate change, subsequent coextinctions are best predicted by the local network of interactions. These coextinctions not only increase the total number of plant species being driven to extinction but also add a bias in the way the major taxonomic and functional groups are pruned.

The bagworm genome reveals a unique fibroin gene that provides high tensile strength

Arthropod silk is known as a versatile tool, and its variability makes it an attractive biomaterial. Eumeta variegata is a bagworm moth (Lepidoptera, Psychidae) that uses silk throughout all life stages. Notably, the bagworm-specific uses of silk include larval development in a bag coated with silk and plant materials and the use of silk attachments to hang pupae. An understanding at the molecular level of bagworm silk, which enables such unique purposes, is an opportunity to expand the possibilities for artificial biomaterial design. However, very little is known about the bagworm fibroin gene and the mechanical properties of bagworm silk. Here, we report the bagworm genome, including a silk fibroin gene. The genome is approximately 700 Mbp in size, and the newly found fibroin gene has a unique repetitive motif. Furthermore, a mechanical property test demonstrates a phylogenetic relationship between the unique motif and tensile strength of bagworm silk.

Modular structure, sequence diversification and appropriate nomenclature of seroins produced in the silk glands of Lepidoptera

Seroins are small lepidopteran silk proteins known to possess antimicrobial activities. Several seroin paralogs and isoforms were identified in studied lepidopteran species and their classification required detailed phylogenetic analysis based on complete and verified cDNA sequences. We sequenced silk gland-specific cDNA libraries from ten species and identified 52 novel seroin cDNAs. The results of this targeted research, combined with data retrieved from available databases, form a dataset representing the major clades of Lepidoptera. The analysis of deduced seroin proteins distinguished three seroin classes (sn1-sn3), which are composed of modules: A (includes the signal peptide), B (rich in charged amino acids) and C (highly variable linker containing proline). The similarities within and between the classes were 31–50% and 22.5–25%, respectively. All species express one, and in exceptional cases two, genes per class, and alternative splicing further enhances seroin diversity. Seroins occur in long versions with the full set of modules (AB₁C₁B₂C₂B₃) and/or in short versions that lack parts or the entire B and C modules. The classes and the modular structure of seroins probably evolved prior to the split between Trichoptera and Lepidoptera. The diversity of seroins is reflected in proposed nomenclature.

Genome sequence of the Japanese oak silk moth, Antheraea yamamai: the first draft genome in the family Saturniidae

Background

Antheraea yamamai, also known as the Japanese oak silk moth, is a wild species of silk moth. Silk produced by A. yamamai, referred to as tensan silk, shows different characteristics such as thickness, compressive elasticity, and chemical resistance compared with common silk produced from the domesticated silkworm, Bombyx mori. Its unique characteristics have led to its use in many research fields including biotechnology and medical science, and the scientific as well as economic importance of the wild silk moth continues to gradually increase. However, no genomic information for the wild silk moth, including A. yamamai, is currently available.

Findings

In order to construct the A. yamamai genome, a total of 147G base pairs using Illumina and Pacbio sequencing platforms were generated, providing 210-fold coverage based on the 700-Mb estimated genome size of A. yamamai. The assembled genome of A. yamamai was 656 Mb (>2 kb) with 3675 scaffolds, and the N50 length of assembly was 739 Kb with a 34.07% GC ratio. Identified repeat elements covered 37.33% of the total genome, and the completeness of the constructed genome assembly was estimated to be 96.7% by Benchmarking Universal Single-Copy Orthologs v2 analysis. A total of 15 481 genes were identified using Evidence Modeler based on the gene prediction results obtained from 3 different methods (ab initio, RNA-seq-based, known-gene-based) and manual curation.

Conclusions

Here we present the genome sequence of A. yamamai, the first genome sequence of the wild silk moth. These results provide valuable genomic information, which will help enrich our understanding of the molecular mechanisms relating to not only specific phenotypes such as wild silk itself but also the genomic evolution of Saturniidae.

Sexual dimorphism in immunity across animals: a meta-analysis

In animals, sex differences in immunity are proposed to shape variation in infection prevalence and intensity among individuals in a population, with females typically expected to exhibit superior immunity due to life-history trade-offs. We performed a systematic meta-analysis to investigate the magnitude and direction of sex differences in immunity and to identify factors that shape sex-biased immunocompetence. In addition to considering taxonomic and methodological effects as moderators, we assessed age-related effects, which are predicted to occur if sex differences in immunity are due to sex-specific resource allocation trade-offs with reproduction. In a meta-analysis of 584 effects from 124 studies, we found that females exhibit a significantly stronger immune response than do males, but the effect size is relatively small, and became non-significant after controlling for phylogeny. Female-biased immunity was more pronounced in adult than immature animals. More recently published studies did not report significantly smaller effect sizes. Among taxonomic and methodological subsets of the data, some of the largest effect sizes were in insects, further supporting previous suggestions that testosterone is not the only potential driver of sex differences in immunity. Our findings challenge the notion of pervasive biases towards female-biased immunity and the role of testosterone in driving these differences.

Biodiversity, Evolution and Ecological Specialization of Baculoviruses: A Treasure Trove for Future Applied Research

The Baculoviridae, a family of insect-specific large DNA viruses, is widely used in both biotechnology and biological control. Its applied value stems from millions of years of evolution influenced by interactions with their hosts and the environment. To understand how ecological interactions have shaped baculovirus diversification, we reconstructed a robust molecular phylogeny using 217 complete genomes and ~580 isolates for which at least one of four lepidopteran core genes was available. We then used a phylogenetic-concept-based approach (mPTP) to delimit 165 baculovirus species, including 38 species derived from new genetic data. Phylogenetic optimization of ecological characters revealed a general pattern of host conservatism punctuated by occasional shifts between closely related hosts and major shifts between lepidopteran superfamilies. Moreover, we found significant phylogenetic conservatism between baculoviruses and the type of plant growth (woody or herbaceous) associated with their insect hosts. In addition, we found that colonization of new ecological niches sometimes led to viral radiation. These macroevolutionary patterns show that besides selection during the infection process, baculovirus diversification was influenced by tritrophic interactions, explained by their persistence on plants and interactions in the midgut during horizontal transmission. This complete eco-evolutionary framework highlights the potential innovations that could still be harnessed from the diversity of baculoviruses.

The complete mitochondrial genome of Clostera anastomosis (Lepidoptera: Notodontidae) and implication for the phylogenetic relationships of Noctuoidea species

In the present study, the complete mitochondrial genome (mitogenome) of Clostera anastomosis (C. anastomosis) has been determined for the first time. The mitogenome is 15,390 base pairs (bp) in length, comprised of 13 protein-coding genes (PCGs), 2 ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs) and one non-coding control region (CR). The gene order shows a typical trnM rearrangement (trnM-trnI-trnQ) compared to ancestral insects (trnI-trnQ-trnM). Almost all the PCGs have the same start codon (ATN) except for cox1 (CGA), and almost all tRNAs have a typical cloverleaf secondary structure except for trnS1. At the beginning of the CR, we found a conserved motif "ATAGA + poly-T" as found in other lepidopteran insects. There are 20 intergenic regions and 11 overlapping regions, ranging from 1 to 53 bp and 1 to 9 bp, respectively. The A + T content is relatively high across the whole mitogenome. The optimal tree topologies of Noctuoidea were given by the dataset consisting of all 13 PCGs from five families (exclude Oenosandridae). Our trees suggested a topology of (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))) and identified that C. anastomosis belongs to Notodontidae.

Origins and structural properties of novel and de novo protein domains during insect evolution

Over long time scales, protein evolution is characterized by modular rearrangements of protein domains. Such rearrangements are mainly caused by gene duplication, fusion and terminal losses. To better understand domain emergence mechanisms we investigated 32 insect genomes covering a speciation gradient ranging from ~ 2 to ~ 390 mya. We use established domain models and foldable domains delineated by hydrophobic cluster analysis (HCA), which does not require homologous sequences, to also identify domains which have likely arisen de novo, that is, from previously noncoding DNA. Our results indicate that most novel domains emerge terminally as they originate from ORF extensions while fewer arise in middle arrangements, resulting from exonization of intronic or intergenic regions. Many novel domains rapidly migrate between terminal or middle positions and single- and multidomain arrangements. Young domains, such as most HCA-defined domains, are under strong selection pressure as they show signals of purifying selection. De novo domains, linked to ancient domains or defined by HCA, have higher degrees of intrinsic disorder and disorder-to-order transition upon binding than ancient domains. However, the corresponding DNA sequences of the novel domains of de novo origins could only rarely be found in sister genomes. We conclude that novel domains are often recruited by other proteins and undergo important structural modifications shortly after their emergence, but evolve too fast to be characterized by cross-species comparisons alone.

The giant butterfly-moth Paysandisia archon has spectrally rich apposition eyes with unique light-dependent photoreceptor dynamics

The palm borer moth Paysandisia archon (Burmeister, 1880) (fam. Castniidae) is a large, diurnally active palm pest. Its compound eyes consist of ~ 20,000 ommatidia and have apposition optics with interommatidial angles below 1°. The ommatidia contain nine photoreceptor cells and appear structurally similar to those in nymphalid butterflies. Two morphological ommatidial types were identified. Using the butterfly numbering scheme, in type I ommatidia, the distal rhabdom consists exclusively of the rhabdomeres of photoreceptors R1–2; the medial rhabdom has contributions from R1–8. The rhabdom in type II ommatidia is distally split into two sub-rhabdoms, with contributions from photoreceptors R2, R3, R5, R6 and R1, R4, R7, R8, respectively; medially, only R3–8 and not R1–2 contribute to the fused rhabdom. In both types, the pigmented bilobed photoreceptors R9 contribute to the rhabdom basally. Their nuclei reside in one of the lobes. Upon light adaptation, in both ommatidial types, the rhabdoms secede from the crystalline cones and pigment granules invade the gap. Intracellular recordings identified four photoreceptor classes with peak sensitivities in the ultraviolet, blue, green and orange wavelength regions (at 360, 465, 550, 580 nm, respectively). We discuss the eye morphology and optics, the photoreceptor spectral sensitivities, and the adaptation to daytime activity from a phylogenetic perspective.

Sex Chromosomes of the Iconic Moth Abraxas grossulariata (Lepidoptera, Geometridae) and Its Congener A. sylvata

The magpie moth, Abraxas grossulariata, is an iconic species in which female heterogamety was discovered at the beginning of the 20th century. However, the sex chromosomes of this species have not yet been cytologically identified. We describe the sex chromosomes of A. grossulariata and its congener, A. sylvata. Although these species split only around 9.5 million years ago, and both species have the expected WZ/ZZ chromosomal system of sex determination and their sex chromosomes share the major ribosomal DNA (rDNA) representing the nucleolar organizer region (NOR), we found major differences between their karyotypes, including between their sex chromosomes. The species differ in chromosome number, which is 2n = 56 in A. grossularita and 2n = 58 in A. sylvata. In addition, A. grossularita autosomes exhibit massive autosomal blocks of heterochromatin, which is a very rare phenomenon in Lepidoptera, whereas the autosomes of A. sylvata are completely devoid of distinct heterochromatin. Their W chromosomes differ greatly. Although they are largely composed of female-specific DNA sequences, as shown by comparative genomic hybridization, cross-species W-chromosome painting revealed considerable sequence differences between them. The results suggest a relatively rapid molecular divergence of Abraxas W chromosomes by the independent spreading of female-specific repetitive sequences.

Resolving Relationships among the Megadiverse Butterflies and Moths with a Novel Pipeline for Anchored Phylogenomics

The advent of next-generation sequencing technology has allowed for thecollection of large portions of the genome for phylogenetic analysis. Hybrid enrichment and transcriptomics are two techniques that leverage next-generation sequencing and have shown much promise. However, methods for processing hybrid enrichment data are still limited. We developed a pipeline for anchored hybrid enrichment (AHE) read assembly, orthology determination, contamination screening, and data processing for sequences flanking the target "probe" region. We apply this approach to study the phylogeny of butterflies and moths (Lepidoptera), a megadiverse group of more than 157,000 described species with poorly understood deep-level phylogenetic relationships. We introduce a new, 855 locus AHE kit for Lepidoptera phylogenetics and compare resulting trees to those from transcriptomes. The enrichment kit was designed from existing genomes, transcriptomes, and expressed sequence tags and was used to capture sequence data from 54 species from 23 lepidopteran families. Phylogenies estimated from AHE data were largely congruent with trees generated from transcriptomes, with strong support for relationships at all but the deepest taxonomic levels. We combine AHE and transcriptomic data to generate a new Lepidoptera phylogeny, representing 76 exemplar species in 42 families. The tree provides robust support for many relationships, including those among the seven butterfly families. The addition of AHE data to an existing transcriptomic dataset lowers node support along the Lepidoptera backbone, but firmly places taxa with AHE data on the phylogeny. Combining taxa sequenced for AHE with existing transcriptomes and genomes resulted in a tree with strong support for (Calliduloidea $+$ Gelechioidea $+$ Thyridoidea) $+$ (Papilionoidea $+$ Pyraloidea $+$ Macroheterocera). To examine the efficacy of AHE at a shallow taxonomic level, phylogenetic analyses were also conducted on a sister group representing a more recent divergence, the Saturniidae and Sphingidae. These analyses utilized sequences from the probe region and data flanking it, nearly doubled the size of the dataset; resulting trees supported new phylogenetics relationships, especially within the Saturniidae and Sphingidae (e.g., Hemarina derived in the latter). We hope that our data processing pipeline, hybrid enrichment gene set, and approach of combining AHE data with transcriptomes will be useful for the broader systematics community.

A global checklist of the Bombycoidea (Insecta: Lepidoptera)

Background

Bombycoidea is an ecologically diverse and speciose superfamily of Lepidoptera. The superfamily includes many model organisms, but the taxonomy and classification of the superfamily has remained largely in disarray. Here we present a global checklist of Bombycoidea. Following Zwick (2008) and Zwick et al. (2011), ten families are recognized: Anthelidae, Apatelodidae, Bombycidae, Brahmaeidae, Carthaeidae, Endromidae, Eupterotidae, Phiditiidae, Saturniidae and Sphingidae. The former families Lemoniidae and Mirinidae are included within Brahmaeidae and Endromidae respectively. The former bombycid subfamilies Oberthueriinae and Prismostictinae are also treated as synonyms of Endromidae, and the former bombycine subfamilies Apatelodinae and Phitditiinae are treated as families.

New information

This checklist represents the first effort to synthesize the current taxonomic treatment of the entire superfamily. It includes 12,159 names and references to their authors, and it accounts for the recent burst in species and subspecies descriptions within family Saturniidae (ca. 1,500 within the past 10 years) and to a lesser extent in Sphingidae (ca. 250 species over the same period). The changes to the higher classification of Saturniidae proposed by Nässig et al. (2015) are rejected as premature and unnecessary. The new tribes, subtribes and genera described by Cooper (2002) are here treated as junior synonyms. We also present a new higher classification of Sphingidae, based on Kawahara et al. (2009), Barber and Kawahara (2013) and a more recent phylogenomic study by Breinholt et al. (2017), as well as a reviewed genus and species level classification, as documented by Kitching (2018).

Conservation and lineage-specific rearrangements in the GOBP/PBP gene complex of distantly related ditrysian Lepidoptera

General odorant binding proteins (GOBPs) and pheromone binding proteins (PBPs) form a monophyletic subfamily of insect odorant binding proteins (OBPs) specific for Lepidoptera, butterflies and moths. The GOBP/PBP genes include six subgroups (GOBP1–2, PBP-A–D) previously reported to form a complex arrayed in a conserved order in representative moths (superfamily Bombycoidea) and butterflies (Nymphalidae). Although our knowledge of lepidopteran genomes has increased greatly recently, the structure of the GOBP/PBP complex has been studied only for species that represent limited lineages of the highly diverged Ditrysia. To understand the evolution of this functionally important gene complex, we determined 69–149 kb genomic sequences that include GOBP2 and five PBP genes in three Ostrinia moths (Pyraloidea), O. nubilalis, O. furnacalis, and O. latipennis, using bacterial artificial chromosome (BAC) and fosmid clones. The structure of the GOBP2/PBP gene cluster was well conserved despite the different sex pheromone composition utilized by the three moths. Five expressed PBP genes in Ostrinia moths were the result of two duplications of PBP-A genes. Surprisingly, an allele containing a fusion gene between tandemly arrayed PBP-A genes was observed in O. nubilalis. We also revealed duplication and intra-chromosomal translocation of the GOBP1 gene in P. xylostella by fluorescence in situ hybridization (FISH) analysis. Additionally, we compared the structure of the GOBP/PBP gene complex of seventeen species covering six superfamilies and twelve families of the lepidopteran clade, Ditrysia, and found the gene order was basically conserved despite the frequent occurrence of lineage-specific gains, losses, inversions and translocations of these genes, compared with their neighboring genes. Our findings support the hypothesis that the structure of the GOBP/PBP gene complex was already established in the common ancestor of Ditrysia.

Lepidoptera genomes: current knowledge, gaps and future directions

Butterflies and moths (Lepidoptera) are one of the most ecologically diverse and speciose insect orders. With recent advances in genomics, new Lepidoptera genomes are regularly being sequenced, and many of them are playing principal roles in genomics studies, particularly in the fields of phylo-genomics and functional genomics. Thus far, assembled genomes are only available for <10 of the 43 Lepidoptera superfamilies. Nearly all are model species, found in the speciose clade Ditrysia. Community support for Lepidoptera genomics is growing with successful management and dissemination of data and analytical tools in centralized databases. With genomic studies quickly becoming integrated with ecological and evolutionary research, the Lepidoptera community will unquestionably benefit from new high-quality reference genomes that are more evenly distributed throughout the order.

Construction of a Species-Level Tree of Life for the Insects and Utility in Taxonomic Profiling

Although comprehensive phylogenies have proven an invaluable tool in ecology and evolution, their construction is made increasingly challenging both by the scale and structure of publically available sequences. The distinct partition between gene-rich (genomic) and species-rich (DNA barcode) data is a feature of data that has been largely overlooked, yet presents a key obstacle to scaling supermatrix analysis. I present a phyloinformatics framework for draft construction of a species-level phylogeny of insects (Class Insecta). Matrix-building requires separately optimized pipelines for nuclear transcriptomic, mitochondrial genomic, and species-rich markers, whereas tree-building requires hierarchical inference in order to capture species-breadth while retaining deep-level resolution. The phylogeny of insects contains 49,358 species, 13,865 genera, 760 families. Deep-level splits largely reflected previous findings for sections of the tree that are data rich or unambiguous, such as inter-ordinal Endopterygota and Dictyoptera, the recently evolved and relatively homogeneous Lepidoptera, Hymenoptera, Brachycera (Diptera), and Cucujiformia (Coleoptera). However, analysis of bias, matrix construction and gene-tree variation suggests confidence in some relationships (such as in Polyneoptera) is less than has been indicated by the matrix bootstrap method. To assess the utility of the insect tree as a tool in query profiling several tree-based taxonomic assignment methods are compared. Using test data sets with existing taxonomic annotations, a tendency is observed for greater accuracy of species-level assignments where using a fixed comprehensive tree of life in contrast to methods generating smaller de novo reference trees. Described herein is a solution to the discrepancy in the way data are fit into supermatrices. The resulting tree facilitates wider studies of insect diversification and application of advanced descriptions of diversity in community studies, among other presumed applications. [Data integration; data mining; insects; phylogenomics; phyloinformatics; tree of life.].

Higher predation risk for insect prey at low latitudes and elevations

Biotic interactions underlie ecosystem structure and function, but predicting interaction outcomes is difficult. We tested the hypothesis that biotic interaction strength increases toward the equator, using a global experiment with model caterpillars to measure predation risk. Across an 11,660-kilometer latitudinal gradient spanning six continents, we found increasing predation toward the equator, with a parallel pattern of increasing predation toward lower elevations. Patterns across both latitude and elevation were driven by arthropod predators, with no systematic trend in attack rates by birds or mammals. These matching gradients at global and regional scales suggest consistent drivers of biotic interaction strength, a finding that needs to be integrated into general theories of herbivory, community organization, and life-history evolution.

The pivotal role of aristaless in development and evolution of diverse antennal morphologies in moths and butterflies

Background

Antennae are multi-segmented appendages and main odor-sensing organs in insects. In Lepidoptera (moths and butterflies), antennal morphologies have diversified according to their ecological requirements. While diurnal butterflies have simple, rod-shaped antennae, nocturnal moths have antennae with protrusions or lateral branches on each antennal segment for high-sensitive pheromone detection. A previous study on the Bombyx mori (silk moth) antenna, forming two lateral branches per segment, during metamorphosis has revealed the dramatic change in expression of antennal patterning genes to segmentally reiterated, branch-associated pattern and abundant proliferation of cells contributing almost all the dorsal half of the lateral branch. Thus, localized cell proliferation possibly controlled by the branch-associated expression of antennal patterning genes is implicated in lateral branch formation. Yet, actual gene function in lateral branch formation in Bombyx mori and evolutionary mechanism of various antennal morphologies in Lepidoptera remain elusive.

Results

We investigated the function of several genes and signaling specifically in lateral branch formation in Bombyx mori by the electroporation-mediated incorporation of siRNAs or morpholino oligomers. Knock down of aristaless, a homeobox gene expressed specifically in the region of abundant cell proliferation within each antennal segment, during metamorphosis resulted in missing or substantial shortening of lateral branches, indicating its importance for lateral branch formation. aristaless expression during metamorphosis was lost by knock down of Distal-less and WNT signaling but derepressed by knock down of Notch signaling, suggesting the strict determination of the aristaless expression domain within each antennal segment by the combinatorial action of them. In addition, analyses of pupal aristaless expression in antennae with various morphologies of several lepidopteran species revealed that the aristaless expression pattern has a striking correlation with antennal shapes, whereas the segmentally reiterated expression pattern was observed irrespective of antennal morphologies.

Conclusions

Our results presented here indicate the significance of aristaless function in lateral branch formation in B. mori and imply that the diversification in the aristaless expression pattern within each antennal segment during metamorphosis is one of the significant determinants of antennal morphologies. According to these findings, we propose a mechanism underlying development and evolution of lepidopteran antennae with various morphologies.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1124-2) contains supplementary material, which is available to authorized users.

A phylogenomic analysis of lichen-feeding tiger moths uncovers evolutionary origins of host chemical sequestration

Host species utilize a variety of defenses to deter feeding, including secondary chemicals. Some phytophagous insects have evolved tolerance to these chemical defenses, and can sequester secondary defense compounds for use against their own predators and parasitoids. While numerous studies have examined plant-insect interactions, little is known about lichen-insect interactions. Our study focused on reconstructing the evolution of lichen phenolic sequestration in the tiger moth tribe Lithosiini (Lepidoptera: Erebidae: Arctiinae), the most diverse lineage of lichen-feeding moths, with 3000 described species. We built an RNA-Seq dataset and examined the adult metabolome for the presence of lichen-derived phenolics. Using the transcriptomic dataset, we recover a well-resolved phylogeny of the Lithosiini, and determine that the metabolomes within species are more similar than those among species. Results from an initial ancestral state reconstruction suggest that the ability to sequester phenolics produced by a single chemical pathway preceded generalist sequestration of phenolics produced by multiple chemical pathways. We conclude that phenolics are consistently and selectively sequestered within Lithosiini. Furthermore, sequestration of compounds from a single chemical pathway may represent a synapomorphy of the tribe, and the ability to sequester phenolics produced by multiple pathways arose later. These findings expand on our understanding of the interactions between Lepidoptera and their lichen hosts.

The deep conservation of the Lepidoptera Z chromosome suggests a non-canonical origin of the W

Moths and butterflies (Lepidoptera) usually have a pair of differentiated WZ sex chromosomes. However, in most lineages outside of the division Ditrysia, as well as in the sister order Trichoptera, females lack a W chromosome. The W is therefore thought to have been acquired secondarily. Here we compare the genomes of three Lepidoptera species (one Dytrisia and two non-Dytrisia) to test three models accounting for the origin of the W: (1) a Z-autosome fusion; (2) a sex chromosome turnover; and (3) a non-canonical mechanism (e.g., through the recruitment of a B chromosome). We show that the gene content of the Z is highly conserved across Lepidoptera (rejecting a sex chromosome turnover) and that very few genes moved onto the Z in the common ancestor of the Ditrysia (arguing against a Z-autosome fusion). Our comparative genomics analysis therefore supports the secondary acquisition of the Lepidoptera W by a non-canonical mechanism, and it confirms the extreme stability of well-differentiated sex chromosomes.

Origins of Aminergic Regulation of Behavior in Complex Insect Social Systems

Neuromodulators are conserved across insect taxa, but how biogenic amines and their receptors in ancestral solitary forms have been co-opted to control behaviors in derived socially complex species is largely unknown. Here we explore patterns associated with the functions of octopamine (OA), serotonin (5-HT) and dopamine (DA) in solitary ancestral insects and their derived functions in eusocial ants, bees, wasps and termites. Synthesizing current findings that reveal potential ancestral roles of monoamines in insects, we identify physiological processes and conserved behaviors under aminergic control, consider how biogenic amines may have evolved to modulate complex social behavior, and present focal research areas that warrant further study.

Cytochrome P450 diversification and hostplant utilization patterns in specialist and generalist moths: Birth, death and adaptation

Across insect genomes, the size of the cytochrome P450 monooxygenase (CYP) gene superfamily varies widely. CYPome size variation has been attributed to reciprocal adaptive radiations in insect detoxification genes in response to plant biosynthetic gene radiations driven by co-evolution between herbivores and their chemically defended hostplants. Alternatively, variation in CYPome size may be due to random "birth-and-death" processes, whereby exponential increase via gene duplications is limited by random decay via gene death or transition via divergence. We examined CYPome diversification in the genomes of seven Lepidoptera species varying in host breadth from monophagous (Bombyx mori) to highly polyphagous (Amyelois transitella). CYPome size largely reflects the size of Clan 3, the clan associated with xenobiotic detoxification, and to some extent phylogenetic age. Consistently across genomes, families CYP6, CYP9 and CYP321 are most diverse and CYP6AB, CYP6AE, CYP6B, CYP9A and CYP9G are most diverse among subfamilies. Higher gene number in subfamilies is due to duplications occurring primarily after speciation and specialization ("P450 blooms"), and the genes are arranged in clusters, indicative of active duplicating loci. In the parsnip webworm, Depressaria pastinacella, gene expression levels in large subfamilies are high relative to smaller subfamilies. Functional and phylogenetic data suggest a correlation between highly dynamic loci (reflective of extensive gene duplication, functionalization and in some cases loss) and the ability of enzymes encoded by these genes to metabolize hostplant defences, consistent with an adaptive, nonrandom process driven by ecological interactions.

Diversification shifts in leafroller moths linked to continental colonization and the rise of angiosperms

Tectonic dynamics and niche availability play intertwined roles in determining patterns of diversification. Such drivers explain the current distribution of many clades, whereas events such as the rise of angiosperms can have more specific impacts, such as on the diversification rates of herbivores. The Tortricidae, a diverse group of phytophagous moths, are ideal for testing the effects of these determinants on the diversification of herbivorous clades. To estimate ancestral areas and diversification patterns in Tortricidae, a complete tribal-level dated tree was inferred using molecular markers (one mitochondrial and five nuclear) and calibrated using fossil constraints. We found that Tortricidae diverged from their sister group c. 120 Myr ago (Ma) and diversified c. 97 Ma, a timeframe synchronous with the rise of angiosperms in the Early-mid Cretaceous. Ancestral areas analysis, based on updated Wallace's biogeographical regions, supports the hypothesis of a Gondwanan origin of Tortricidae in the South American plate. We also detected an increase in speciation rate that coincided with the peak of angiosperm diversification in the Cretaceous. This in turn probably was further heightened by continental colonization of the Palaeotropics when angiosperms became dominant by the end of the Late Cretaceous.

Evolutionary escalation: the bat-moth arms race

Echolocation in bats and high-frequency hearing in their insect prey make bats and insects an ideal system for studying the sensory ecology and neuroethology of predator-prey interactions. Here, we review the evolutionary history of bats and eared insects, focusing on the insect order Lepidoptera, and consider the evidence for antipredator adaptations and predator counter-adaptations. Ears evolved in a remarkable number of body locations across insects, with the original selection pressure for ears differing between groups. Although cause and effect are difficult to determine, correlations between hearing and life history strategies in moths provide evidence for how these two variables influence each other. We consider life history variables such as size, sex, circadian and seasonal activity patterns, geographic range and the composition of sympatric bat communities. We also review hypotheses on the neural basis for anti-predator behaviours (such as evasive flight and sound production) in moths. It is assumed that these prey adaptations would select for counter-adaptations in predatory bats. We suggest two levels of support for classifying bat traits as counter-adaptations: traits that allow bats to eat more eared prey than expected based on their availability in the environment provide a low level of support for counter-adaptations, whereas traits that have no other plausible explanation for their origination and maintenance than capturing defended prey constitute a high level of support. Specific predator counter-adaptations include calling at frequencies outside the sensitivity range of most eared prey, changing the pattern and frequency of echolocation calls during prey pursuit, and quiet, or 'stealth', echolocation.

Conflicting phylogenomic signals reveal a pattern of reticulate evolution in a recent high-Andean diversification (Asteraceae: Astereae: Diplostephium)

High-throughput sequencing is helping biologists to overcome the difficulties of inferring the phylogenies of recently diverged taxa. The present study analyzes the phylogenetic signal of genomic regions with different inheritance patterns using genome skimming and ddRAD-seq in a species-rich Andean genus (Diplostephium) and its allies. We analyzed the complete nuclear ribosomal cistron, the complete chloroplast genome, a partial mitochondrial genome, and a nuclear-ddRAD matrix separately with phylogenetic methods. We applied several approaches to understand the causes of incongruence among datasets, including simulations and the detection of introgression using the D-statistic (ABBA-BABA test). We found significant incongruence among the nuclear, chloroplast, and mitochondrial phylogenies. The strong signal of hybridization found by simulations and the D-statistic among genera and inside the main clades of Diplostephium indicate reticulate evolution as a main cause of phylogenetic incongruence. Our results add evidence for a major role of reticulate evolution in events of rapid diversification. Hybridization and introgression confound chloroplast and mitochondrial phylogenies in relation to the species tree as a result of the uniparental inheritance of these genomic regions. Practical implications regarding the prevalence of hybridization are discussed in relation to the phylogenetic method.

Computational identification of the selenocysteine tRNA (tRNASec) in genomes

Selenocysteine (Sec) is known as the 21st amino acid, a cysteine analogue with selenium replacing sulphur. Sec is inserted co-translationally in a small fraction of proteins called selenoproteins. In selenoprotein genes, the Sec specific tRNA (tRNA^Sec) drives the recoding of highly specific UGA codons from stop signals to Sec. Although found in organisms from the three domains of life, Sec is not universal. Many species are completely devoid of selenoprotein genes and lack the ability to synthesize Sec. Since tRNA^Sec is a key component in selenoprotein biosynthesis, its efficient identification in genomes is instrumental to characterize the utilization of Sec across lineages. Available tRNA prediction methods fail to accurately predict tRNA^Sec, due to its unusual structural fold. Here, we present Secmarker, a method based on manually curated covariance models capturing the specific tRNA^Sec structure in archaea, bacteria and eukaryotes. We exploited the non-universality of Sec to build a proper benchmark set for tRNA^Sec predictions, which is not possible for the predictions of other tRNAs. We show that Secmarker greatly improves the accuracy of previously existing methods constituting a valuable tool to identify tRNA^Sec genes, and to efficiently determine whether a genome contains selenoproteins. We used Secmarker to analyze a large set of fully sequenced genomes, and the results revealed new insights in the biology of tRNA^Sec, led to the discovery of a novel bacterial selenoprotein family, and shed additional light on the phylogenetic distribution of selenoprotein containing genomes. Secmarker is freely accessible for download, or online analysis through a web server at http://secmarker.crg.cat.

Most proteins are made of twenty amino acids. However, there is a small group of proteins that incorporate a 21st amino acid, Selenocysteine (Sec). These proteins are called selenoproteins and are present in some, but not all, species from the three domains of life. Sec is inserted in selenoproteins in response to the UGA codon, normally a stop codon. A Sec specific tRNA (tRNA^Sec), which only exists in the organisms that synthesize selenoproteins recognizes the UGA codon. tRNA^Sec is not only indispensable for Sec incorporation into selenoproteins, but also for Sec synthesis, since Sec is synthesized on its own tRNA. The structure of tRNA^Sec differs from that of canonical tRNAs, and general tRNA detection methods fail to accurately predict it. We developed Secmarker, a tRNA^Sec specific identification tool based on the characteristic structural features of the tRNA^Sec. Our benchmark shows that Secmarker produces nearly flawless tRNA^Sec predictions. We used Secmarker to scan all currently available genome sequences. The analysis of the highly accurate predictions obtained revealed new insights into the biology of tRNA^Sec.

Phylogeny and Evolution of Lepidoptera

Until recently, deep-level phylogeny in Lepidoptera, the largest single radiation of plant-feeding insects, was very poorly understood. Over the past two decades, building on a preceding era of morphological cladistic studies, molecular data have yielded robust initial estimates of relationships both within and among the ∼43 superfamilies, with unsolved problems now yielding to much larger data sets from high-throughput sequencing. Here we summarize progress on lepidopteran phylogeny since 1975, emphasizing the superfamily level, and discuss some resulting advances in our understanding of lepidopteran evolution.

Detailed Morphology of All Life Stages of the Agave Red Worm, Comadia redtenbacheri (Hammerschmidt) (Lepidoptera: Cossidae)

The agave red worm, Comadia redtenbacheri (Hammerschmidt), is an important source of food and income in Mexico. Despite its importance, several aspects of its biology, morphology, and behavior remain poorly studied. In this work, we describe and illustrate the morphology of all the life stages that may aid in understanding certain aspects of its biology. To obtain all life stages, last instar larvae were collected from agave plants and allowed to pupate; after the adults emerged, they were allowed to mate and oviposit. The frenulum is longer in males; epiphysis I is longer in females than in males; the abdomen bears two types of tubercles of unknown function. Eggs present a reticulate chorion; primary rosette cells are highly variable in shape; the micropylar formula is (10-14): (12-13). First instar larvae are white, becoming red as they develop; L3 in the prothorax is subprimary; the SV setal group in A1 is comprised of only SV1 on first instar larvae; last instars have several secondary setae. Pupae are adecticous and obtect; there are rows of spines on the dorsum of the abdomen. The biological significance of some of the findings is discussed.

The rapidly changing landscape of insect phylogenetics

Insect phylogenetics is being profoundly changed by many innovations. Although rapid developments in genomics have center stage, key progress has been made in phenomics, field and museum science, digital databases and pipelines, analytical tools, and the culture of science. The importance of these methodological and cultural changes to the pace of inference of the hexapod Tree of Life is discussed. The innovations have the potential, when synthesized and mobilized in ways as yet unforeseen, to shine light on the million or more clades in insects, and infer their composition with confidence. There are many challenges to overcome before insects can enter the 'phylocognisant age', but because of the promise of genomics, phenomics, and informatics, that is now an imaginable future.

Farnesyl biliverdins IXα are novel ligands of biliproteins from moths of the Noctuoidea superfamily: A chemosystematic view of the Lepidoptera

Bilins, derived from biliverdin IXα, are known from animals, plants and microorganisms, where they play vital roles as light-absorbing pigments. Bilins occur also in many insects. Recently, we discovered in insects a novel structural type of bilins with a farnesyl substituent at pyrrole ring A of biliverdin IXα. The first of these unusual bilins with a molecular mass of 852 (C₄₈H₆₀O₁₀N₄) was identified in Cerura vinula, subsequently in Spodoptera littoralis; both species are members of the Noctuoidea superfamily of moths. From an evolutionary point of view, it was of interest to examine other species and families of this monophyletic clade. Here, we show that other moths species in this clade (three Notodontidae species, one Erebidae species, and one Noctuidae species) have farnesylated biliverdins IXα that are present as a mixture of three bilins, differing by the number of oxygen atoms (O_8-10). These bilins are associated with typical hemolymph storage proteins, which were identified by mass spectroscopic sequencing of tryptic peptides as arylphorins (a class of 500-kDa hexamerins) in the Notodontidae and Erebidae families, and as 350-kDa very high-density lipoproteins in the Noctuidae family. Circular dichroism spectroscopy revealed that the bilins adopt opposite conformations in complex with the two different classes of proteins. At present, farnesylated biliverdins and IXα-isomers of bilins in general are known only from species of the Noctuoidea clade; the sister clades of Bombycoidea and Papilionoidea synthesise the IXγ-isomer of biliverdin and derivatives thereof.

Conserved chemosensory proteins in the proboscis and eyes of Lepidoptera

Odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) are endowed with several different functions besides being carriers for pheromones and odorants. Based on a previous report of a CSP acting as surfactant in the proboscis of the moth Helicoverpa armigera, we revealed the presence of orthologue proteins in two other moths Plutella xylostella and Chilo suppressalis, as well as two butterflies Papilio machaon and Pieris rapae, using immunodetection and proteomic analysis. The unusual conservation of these proteins across large phylogenetic distances indicated a common specific function for these CSPs. This fact prompted us to search for other functions of these proteins and discovered that CSPs are abundantly expressed in the eyes of H. armigera and possibly involved as carriers for carotenoids and visual pigments. This hypothesis is supported by ligand-binding experiments and docking simulations with retinol and β-carotene. This last orange pigment, occurring in many fruits and vegetables, is an antioxidant and the precursor of visual pigments. We propose that structurally related CSPs solubilise nutritionally important carotenoids in the proboscis, while they act as carriers of both β-carotene and its derived products 3-hydroxyretinol and 3-hydroxyretinal in the eye. The use of soluble olfactory proteins, such as CSPs, as carriers for visual pigments in insects, here reported for the first time, parallels the function of retinol-binding protein in vertebrates, a lipocalin structurally related to vertebrate odorant-binding proteins.

Cool habitats support darker and bigger butterflies in Australian tropical forests

Morphology mediates the relationship between an organism's body temperature and its environment. Dark organisms, for example, tend to absorb heat more quickly than lighter individuals, which could influence their responses to temperature. Therefore, temperature‐related traits such as morphology may affect patterns of species abundance, richness, and community assembly across a broad range of spatial scales. In this study, we examined variation in color lightness and body size within butterfly communities across hot and cool habitats in the tropical woodland–rainforest ecosystems of northeast Queensland, Australia. Using thermal imaging, we documented the absorption of solar radiation relative to color lightness and wingspan and then built a phylogenetic tree based on available sequences to analyze the effects of habitat on these traits within a phylogenetic framework. In general, darker and larger individuals were more prevalent in cool, closed‐canopy rainforests than in immediately adjacent and hotter open woodlands. In addition, darker and larger butterflies preferred to be active in the shade and during crepuscular hours, while lighter and smaller butterflies were more active in the sun and midday hours—a pattern that held after correcting for phylogeny. Our ex situ experiment supported field observations that dark and large butterflies heated up faster than light and small butterflies under standardized environmental conditions. Our results show a thermal consequence of butterfly morphology across habitats and how environmental factors at a microhabitat scale may affect the distribution of species based on these traits. Furthermore, this study highlights how butterfly species might differentially respond to warming based on ecophysiological traits and how thermal refuges might emerge at microclimatic and habitat scales.