Intrasexual cuticular hydrocarbon dimorphism in a wasp sheds light on hydrocarbon biosynthesis genes in Hymenoptera.

Cuticular hydrocarbons (CHCs) cover the cuticle of insects and serve as desiccation barrier and as semiochemicals. While the main enzymatic steps of CHC biosynthesis are well understood, few of the underlying genes have been identified. Here we show how exploitation of intrasexual CHC dimorphism in a mason wasp, Odynerus spinipes, in combination with whole-genome sequencing and comparative transcriptomics facilitated identification of such genes. RNAi-mediated knockdown of twelve candidate gene orthologs in the honey bee, Apis mellifera, confirmed nine genes impacting CHC profile composition. Most of them have predicted functions consistent with current knowledge of CHC metabolism. However, we found first-time evidence for a fatty acid amide hydrolase also influencing CHC profile composition. In situ hybridization experiments furthermore suggest trophocytes participating in CHC biosynthesis. Our results set the base for experimental CHC profile manipulation in Hymenoptera and imply that the evolutionary origin of CHC biosynthesis predates the arthropods' colonization of land.

Transcriptomics provides a robust framework for the relationships of the major clades of cladobranch sea slugs (Mollusca, Gastropoda, Heterobranchia), but fails to resolve the position of the enigmatic genus Embletonia.

BACKGROUND: The soft-bodied cladobranch sea slugs represent roughly half of the biodiversity of marine nudibranch molluscs on the planet. Despite their global distribution from shallow waters to the deep sea, from tropical into polar seas, and their important role in marine ecosystems and for humans (as targets for drug discovery), the evolutionary history of cladobranch sea slugs is not yet fully understood. RESULTS: To enlarge the current knowledge on the phylogenetic relationships, we generated new transcriptome data for 19 species of cladobranch sea slugs and two additional outgroup taxa (Berthella plumula and Polycera quadrilineata). We complemented our taxon sampling with previously published transcriptome data, resulting in a final data set covering 56 species from all but one accepted cladobranch superfamilies. We assembled all transcriptomes using six different assemblers, selecting those assemblies that provided the largest amount of potentially phylogenetically informative sites. Quality-driven compilation of data sets resulted in four different supermatrices: two with full coverage of genes per species (446 and 335 single-copy protein-coding genes, respectively) and two with a less stringent coverage (667 genes with 98.9% partition coverage and 1767 genes with 86% partition coverage, respectively). We used these supermatrices to infer statistically robust maximum-likelihood trees. All analyses, irrespective of the data set, indicate maximal statistical support for all major splits and phylogenetic relationships at the family level. Besides the questionable position of Noumeaella rubrofasciata, rendering the Facelinidae as polyphyletic, the only notable discordance between the inferred trees is the position of Embletonia pulchra. Extensive testing using Four-cluster Likelihood Mapping, Approximately Unbiased tests, and Quartet Scores revealed that its position is not due to any informative phylogenetic signal, but caused by confounding signal. CONCLUSIONS: Our data matrices and the inferred trees can serve as a solid foundation for future work on the taxonomy and evolutionary history of Cladobranchia. The placement of E. pulchra, however, proves challenging, even with large data sets and various optimization strategies. Moreover, quartet mapping results show that confounding signal present in the data is sufficient to explain the inferred position of E. pulchra, again leaving its phylogenetic position as an enigma.

Evolutionary history and divergence times of Odonata (dragonflies and damselflies) revealed through transcriptomics.

Dragonflies and damselflies are among the earliest flying insects with extant representatives. However, unraveling details of their long evolutionary history, such as egg laying (oviposition) strategies, is impeded by unresolved phylogenetic relationships, particularly in damselflies. Here we present a transcriptome-based phylogenetic reconstruction of Odonata, analyzing 2,980 protein-coding genes in 105 species representing nearly all the order's families. All damselfly and most dragonfly families are recovered as monophyletic. Our data suggest a sister relationship between dragonfly families of Gomphidae and Petaluridae. According to our divergence time estimates, both crown-Zygoptera and -Anisoptera arose during the late Triassic. Egg-laying with a reduced ovipositor apparently evolved in dragonflies during the late Jurassic/early Cretaceous. Lastly, we also test the impact of fossil choice and placement, particularly, of the extinct fossil species, †Triassolestodes asiaticus, and †Proterogomphus renateae on divergence time estimates. We find placement of †Proterogomphus renateae to be much more impactful than †Triassolestodes asiaticus.

The complete mitochondrial genome of the photosymbiotic sea slug Berghia stephanieae (Valdés, 2005) (Gastropoda, Nudibranchia).

Berghia stephanieae (Nudibranchia, Cladobranchia) is a photosymbiotic sea slug that feeds exclusively on sea anemones from the genus Exaiptasia. It then specifically incorporates dinoflagellates belonging to the Symbiodiniaceae obtained from their prey. Here, we present the complete mitochondrial genome sequence of B. stephanieae combining Oxford Nanopore long read and Illumina short-read sequencing data. The mitochondrial genome has a total length of 14,786 bp, it contains the 13 protein-encoding genes, 23 tRNAs, and two rRNAs and is similar to other nudibranchs except for the presence of a duplicated tRNA-Ser 1.

Midgut transcriptome assessment of the cockroach-hunting wasp Ampulex compressa (Apoidea: Ampulicidae).

The emerald jewel wasp Ampulex compressa (Hymenoptera: Ampulicidae) is a solitary wasp that is widely known for its specialized hunting of cockroaches as larvae provision. Adult wasps mainly feed on pollen and nectar, while their larvae feed on the cockroachs' body, first as ecto- and later as endoparsitoids. Little is known about the expression of digestive, detoxification and stress-response-related genes in the midgut of A. compressa, or about its transcriptional versatility between life stages. To identify gut-biased genes related to digestion, detoxification, and stress response, we explored the midgut transcriptome of lab-reared A. compressa, for both adults and larvae, by focusing on the top 100 significantly up- and down-regulated genes. From the top 100 significantly differentially expressed genes (DEGs), we identified 39 and 36 DEGs putatively related to digestion and detoxification in the adult wasps and larvae, respectively. The two carbohydrases alpha-glucosidase (containing an alpha-amylase domain) and glycosyl hydrolase family 31, as well as the two proteinases chymotrypsin and trypsin, revealed the highest gene diversity. We identified six significant DEGs related to detoxification, which comprise glutathione S-transferase, cytochrome P450s and UDP-glucuronosyltransferase. The gene expression levels that were significantly expressed in both life stages vary strongly between life stages, as found in genes encoding for chymotrypsin and trypsin or glycosyl hydrolases family 31. The number of genes related to alpha-glucosidase, glycosyl hydrolase family 31, and cytochrome P450s was found to be similar across nine reference hymenopteran species, except for the identified glycosyl hydrolase family 31 gene, which was absent in all reference bee species. Phylogenetic analyses of the latter candidate genes revealed that they cluster together with their homologous genes found in the reference hymenopteran species. These identified candidate genes provide a basis for future comparative genomic and proteomic studies on (ontogenetic) dietary transitions in Hymenoptera.

Viromics of extant insect orders unveil the evolution of the flavi-like superfamily.

Insects are the most diversified and species-rich group of animals and harbor an immense diversity of viruses. Several taxa in the flavi-like superfamily, such as the genus Flavivirus, are associated with insects; however, systematic studies on insect virus genetic diversity are lacking, limiting our understanding of the evolution of the flavi-like superfamily. Here, we examined the diversity of flavi-like viruses within the most complete and up-to-date insect transcriptome collection comprising 1,243 insect species by employing a Flaviviridae RdRp profile hidden Markov model search. We identified seventy-six viral sequences in sixty-one species belonging to seventeen insect, one entognathan, and one arachnidan orders. Phylogenetic analyses revealed that twenty-seven sequences fell within the Flaviviridae phylogeny but did not group with established genera. Despite the large diversity of insect hosts studied, we only detected one virus in a blood-feeding insect, which branched within the genus Flavivirus, indicating that this genus likely diversified only in hematophagous arthropods. Nine new jingmenviruses with novel host associations were identified. One of the jingmenviruses established a deep rooting lineage additional to the insect- and tick-associated clades. Segment co-segregation phylogenies support the separation of tick- and insect-associated groups within jingmenviruses, with evidence for segment reassortment. In addition, fourteen viruses grouped with unclassified flaviviruses encompassing genome length of up to 20 kb. Species-specific clades for Hymenopteran- and Orthopteran-associated viruses were identified. Forty-nine viruses populated three highly diversified clades in distant relationship to Tombusviridae, a plant-infecting virus family, suggesting the detection of three previously unknown insect-associated families that contributed to tombusvirus evolution.

Beyond Drosophila: resolving the rapid radiation of schizophoran flies with phylotranscriptomics.

BACKGROUND: The most species-rich radiation of animal life in the 66 million years following the Cretaceous extinction event is that of schizophoran flies: a third of fly diversity including Drosophila fruit fly model organisms, house flies, forensic blow flies, agricultural pest flies, and many other well and poorly known true flies. Rapid diversification has hindered previous attempts to elucidate the phylogenetic relationships among major schizophoran clades. A robust phylogenetic hypothesis for the major lineages containing these 55,000 described species would be critical to understand the processes that contributed to the diversity of these flies. We use protein encoding sequence data from transcriptomes, including 3145 genes from 70 species, representing all superfamilies, to improve the resolution of this previously intractable phylogenetic challenge. RESULTS: Our results support a paraphyletic acalyptrate grade including a monophyletic Calyptratae and the monophyly of half of the acalyptrate superfamilies. The primary branching framework of Schizophora is well supported for the first time, revealing the primarily parasitic Pipunculidae and Sciomyzoidea stat. rev. as successive sister groups to the remaining Schizophora. Ephydroidea, Drosophila's superfamily, is the sister group of Calyptratae. Sphaeroceroidea has modest support as the sister to all non-sciomyzoid Schizophora. We define two novel lineages corroborated by morphological traits, the 'Modified Oviscapt Clade' containing Tephritoidea, Nerioidea, and other families, and the 'Cleft Pedicel Clade' containing Calyptratae, Ephydroidea, and other families. Support values remain low among a challenging subset of lineages, including Diopsidae. The placement of these families remained uncertain in both concatenated maximum likelihood and multispecies coalescent approaches. Rogue taxon removal was effective in increasing support values compared with strategies that maximise gene coverage or minimise missing data. CONCLUSIONS: Dividing most acalyptrate fly groups into four major lineages is supported consistently across analyses. Understanding the fundamental branching patterns of schizophoran flies provides a foundation for future comparative research on the genetics, ecology, and biocontrol.

Identification of scavenger receptors and thrombospondin-type-1 repeat proteins potentially relevant for plastid recognition in Sacoglossa.

Functional kleptoplasty is a photosymbiotic relationship, in which photosynthetically active chloroplasts serve as an intracellular symbiont for a heterotrophic host. Among Metazoa, functional kleptoplasty is only found in marine sea slugs belonging to the Sacoglossa and recently described in Rhabdocoela worms. Although functional kleptoplasty has been intensively studied in Sacoglossa, the fundamentals of the specific recognition of the chloroplasts and their subsequent incorporation are unknown. The key to ensure the initiation of any symbiosis is the ability to specifically recognize the symbiont and to differentiate a symbiont from a pathogen. For instance, in photosymbiotic cnidarians, several studies have shown that the host innate immune system, in particular scavenger receptors (SRs) and thrombospondin-type-1 repeat (TSR) protein superfamily, is playing a major role in the process of recognizing and differentiating symbionts from pathogens. In the present study, SRs and TSRs of three Sacoglossa sea slugs, Elysia cornigera, Elysia timida, and Elysia chlorotica, were identified by translating available transcriptomes into potential proteins and searching for receptor specific protein and/or transmembrane domains. Both receptors classes are highly diverse in the slugs, and many new domain arrangements for each receptor class were found. The analyses of the gene expression of these three species provided a set of species-specific candidate genes, that is, SR-Bs, SR-Es, C-type lectins, and TSRs, that are potentially relevant for the recognition of kleptoplasts. The results set the base for future experimental studies to understand if and how these candidate receptors are indeed involved in chloroplast recognition.

Author Correction: Phylogenomic analysis sheds light on the evolutionary pathways towards acoustic communication in Orthoptera.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Four myriapod relatives - but who are sisters? No end to debates on relationships among the four major myriapod subgroups.

BACKGROUND: Phylogenetic relationships among the myriapod subgroups Chilopoda, Diplopoda, Symphyla and Pauropoda are still not robustly resolved. The first phylogenomic study covering all subgroups resolved phylogenetic relationships congruently to morphological evidence but is in conflict with most previously published phylogenetic trees based on diverse molecular data. Outgroup choice and long-branch attraction effects were stated as possible explanations for these incongruencies. In this study, we addressed these issues by extending the myriapod and outgroup taxon sampling using transcriptome data. RESULTS: We generated new transcriptome data of 42 panarthropod species, including all four myriapod subgroups and additional outgroup taxa. Our taxon sampling was complemented by published transcriptome and genome data resulting in a supermatrix covering 59 species. We compiled two data sets, the first with a full coverage of genes per species (292 single-copy protein-coding genes), the second with a less stringent coverage (988 genes). We inferred phylogenetic relationships among myriapods using different data types, tree inference, and quartet computation approaches. Our results unambiguously support monophyletic Mandibulata and Myriapoda. Our analyses clearly showed that there is strong signal for a single unrooted topology, but a sensitivity of the position of the internal root on the choice of outgroups. However, we observe strong evidence for a clade Pauropoda+Symphyla, as well as for a clade Chilopoda+Diplopoda. CONCLUSIONS: Our best quartet topology is incongruent with current morphological phylogenies which were supported in another phylogenomic study. AU tests and quartet mapping reject the quartet topology congruent to trees inferred with morphological characters. Moreover, quartet mapping shows that confounding signal present in the data set is sufficient to explain the weak signal for the quartet topology derived from morphological characters. Although outgroup choice affects results, our study could narrow possible trees to derivatives of a single quartet topology. For highly disputed relationships, we propose to apply a series of tests (AU and quartet mapping), since results of such tests allow to narrow down possible relationships and to rule out confounding signal.

Phylogenomic analysis sheds light on the evolutionary pathways towards acoustic communication in Orthoptera.

Acoustic communication is enabled by the evolution of specialised hearing and sound producing organs. In this study, we performed a large-scale macroevolutionary study to understand how both hearing and sound production evolved and affected diversification in the insect order Orthoptera, which includes many familiar singing insects, such as crickets, katydids, and grasshoppers. Using phylogenomic data, we firmly establish phylogenetic relationships among the major lineages and divergence time estimates within Orthoptera, as well as the lineage-specific and dynamic patterns of evolution for hearing and sound producing organs. In the suborder Ensifera, we infer that forewing-based stridulation and tibial tympanal ears co-evolved, but in the suborder Caelifera, abdominal tympanal ears first evolved in a non-sexual context, and later co-opted for sexual signalling when sound producing organs evolved. However, we find little evidence that the evolution of hearing and sound producing organs increased diversification rates in those lineages with known acoustic communication.

Correction to: An integrative phylogenomic approach to elucidate the evolutionary history and divergence times of Neuropterida (Insecta: Holometabola).

An amendment to this paper has been published and can be accessed via the original article.

ABCB transporters in a leaf beetle respond to sequestered plant toxins.

Phytophagous insects can tolerate and detoxify toxic compounds present in their host plants and have evolved intricate adaptations to this end. Some insects even sequester the toxins for their defence. This necessitates specific mechanisms, especially carrier proteins that regulate uptake and transport to specific storage sites or protect sensitive tissues from noxious compounds. We identified three ATP-binding cassette subfamily B (ABCB) transporters from the transcriptome of the cardenolide-sequestering leaf beetle Chrysochus auratus and analysed their functional role in the sequestration process. These were heterologously expressed and tested for their ability to interact with various potential substrates: verapamil (standard ABCB substrate), the cardenolides digoxin (commonly used), cymarin (present in the species's host plant) and calotropin (present in the ancestral host plants). Verapamil stimulated all three ABCBs and each was activated by at least one cardenolide, however, they differed as to which they were activated by. While the expression of the most versatile transporter fits with a protective role in the blood-brain barrier, the one specific for cymarin shows an extreme abundance in the elytra, coinciding with the location of the defensive glands. Our data thus suggest a key role of ABCBs in the transport network needed for cardenolide sequestration.

An integrative phylogenomic approach to elucidate the evolutionary history and divergence times of Neuropterida (Insecta: Holometabola).

BACKGROUND: The latest advancements in DNA sequencing technologies have facilitated the resolution of the phylogeny of insects, yet parts of the tree of Holometabola remain unresolved. The phylogeny of Neuropterida has been extensively studied, but no strong consensus exists concerning the phylogenetic relationships within the order Neuroptera. Here, we assembled a novel transcriptomic dataset to address previously unresolved issues in the phylogeny of Neuropterida and to infer divergence times within the group. We tested the robustness of our phylogenetic estimates by comparing summary coalescent and concatenation-based phylogenetic approaches and by employing different quartet-based measures of phylogenomic incongruence, combined with data permutations. RESULTS: Our results suggest that the order Raphidioptera is sister to Neuroptera + Megaloptera. Coniopterygidae is inferred as sister to all remaining neuropteran families suggesting that larval cryptonephry could be a ground plan feature of Neuroptera. A clade that includes Nevrorthidae, Osmylidae, and Sisyridae (i.e. Osmyloidea) is inferred as sister to all other Neuroptera except Coniopterygidae, and Dilaridae is placed as sister to all remaining neuropteran families. Ithonidae is inferred as the sister group of monophyletic Myrmeleontiformia. The phylogenetic affinities of Chrysopidae and Hemerobiidae were dependent on the data type analyzed, and quartet-based analyses showed only weak support for the placement of Hemerobiidae as sister to Ithonidae + Myrmeleontiformia. Our molecular dating analyses suggest that most families of Neuropterida started to diversify in the Jurassic and our ancestral character state reconstructions suggest a primarily terrestrial environment of the larvae of Neuropterida and Neuroptera. CONCLUSION: Our extensive phylogenomic analyses consolidate several key aspects in the backbone phylogeny of Neuropterida, such as the basal placement of Coniopterygidae within Neuroptera and the monophyly of Osmyloidea. Furthermore, they provide new insights into the timing of diversification of Neuropterida. Despite the vast amount of analyzed molecular data, we found that certain nodes in the tree of Neuroptera are not robustly resolved. Therefore, we emphasize the importance of integrating the results of morphological analyses with those of sequence-based phylogenomics. We also suggest that comparative analyses of genomic meta-characters should be incorporated into future phylogenomic studies of Neuropterida.

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera.

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO2 detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

Re-assessing the diversity of negative strand RNA viruses in insects.

The spectrum of viruses in insects is important for subjects as diverse as public health, veterinary medicine, food production, and biodiversity conservation. The traditional interest in vector-borne diseases of humans and livestock has drawn the attention of virus studies to hematophagous insect species. However, these represent only a tiny fraction of the broad diversity of Hexapoda, the most speciose group of animals. Here, we systematically probed the diversity of negative strand RNA viruses in the largest and most representative collection of insect transcriptomes from samples representing all 34 extant orders of Hexapoda and 3 orders of Entognatha, as well as outgroups, altogether representing 1243 species. Based on profile hidden Markov models we detected 488 viral RNA-directed RNA polymerase (RdRp) sequences with similarity to negative strand RNA viruses. These were identified in members of 324 arthropod species. Selection for length, quality, and uniqueness left 234 sequences for analyses, showing similarity to genomes of viruses classified in Bunyavirales (n = 86), Articulavirales (n = 54), and several orders within Haploviricotina (n = 94). Coding-complete genomes or nearly-complete subgenomic assemblies were obtained in 61 cases. Based on phylogenetic topology and the availability of coding-complete genomes we estimate that at least 20 novel viral genera in seven families need to be defined, only two of them monospecific. Seven additional viral clades emerge when adding sequences from the present study to formerly monospecific lineages, potentially requiring up to seven additional genera. One long sequence may indicate a novel family. For segmented viruses, cophylogenies between genome segments were generally improved by the inclusion of viruses from the present study, suggesting that in silico misassembly of segmented genomes is rare or absent. Contrary to previous assessments, significant virus-host codivergence was identified in major phylogenetic lineages based on two different approaches of codivergence analysis in a hypotheses testing framework. In spite of these additions to the known spectrum of viruses in insects, we caution that basing taxonomic decisions on genome information alone is challenging due to technical uncertainties, such as the inability to prove integrity of complete genome assemblies of segmented viruses.

The evolution and genomic basis of beetle diversity.

The order Coleoptera (beetles) is arguably the most speciose group of animals, but the evolutionary history of beetles, including the impacts of plant feeding (herbivory) on beetle diversification, remain poorly understood. We inferred the phylogeny of beetles using 4,818 genes for 146 species, estimated timing and rates of beetle diversification using 89 genes for 521 species representing all major lineages and traced the evolution of beetle genes enabling symbiont-independent digestion of lignocellulose using 154 genomes or transcriptomes. Phylogenomic analyses of these uniquely comprehensive datasets resolved previously controversial beetle relationships, dated the origin of Coleoptera to the Carboniferous, and supported the codiversification of beetles and angiosperms. Moreover, plant cell wall-degrading enzymes (PCWDEs) obtained from bacteria and fungi via horizontal gene transfers may have been key to the Mesozoic diversification of herbivorous beetles-remarkably, both major independent origins of specialized herbivory in beetles coincide with the first appearances of an arsenal of PCWDEs encoded in their genomes. Furthermore, corresponding (Jurassic) diversification rate increases suggest that these novel genes triggered adaptive radiations that resulted in nearly half of all living beetle species. We propose that PCWDEs enabled efficient digestion of plant tissues, including lignocellulose in cell walls, facilitating the evolution of uniquely specialized plant-feeding habits, such as leaf mining and stem and wood boring. Beetle diversity thus appears to have resulted from multiple factors, including low extinction rates over a long evolutionary history, codiversification with angiosperms, and adaptive radiations of specialized herbivorous beetles following convergent horizontal transfers of microbial genes encoding PCWDEs.

Improved annotation of protein-coding genes boundaries in metazoan mitochondrial genomes.

With the rapid increase of sequenced metazoan mitochondrial genomes, a detailed manual annotation is becoming more and more infeasible. While it is easy to identify the approximate location of protein-coding genes within mitogenomes, the peculiar processing of mitochondrial transcripts, however, makes the determination of precise gene boundaries a surprisingly difficult problem. We have analyzed the properties of annotated start and stop codon positions in detail, and use the inferred patterns to devise a new method for predicting gene boundaries in de novo annotations. Our method benefits from empirically observed prevalances of start/stop codons and gene lengths, and considers the dependence of these features on variations of genetic codes. Albeit not being perfect, our new approach yields a drastic improvement in the accuracy of gene boundaries and upgrades the mitochondrial genome annotation server MITOS to an even more sophisticated tool for fully automatic annotation of metazoan mitochondrial genomes.

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

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

Robust reference gene design and validation for expression studies in the large milkweed bug, Oncopeltus fasciatus, upon cardiac glycoside stress.

Despite its history as a developmental and evolutionary model organism, gene expression analysis in the large milkweed bug, Oncopeltus fasciatus, has rarely been explored using quantitative real-time PCR. The strength of this method depends greatly on the endogenous controls used for normalization, which are lacking for the milkweed bug system. Here, to fill in this gap in our knowledge, we validated the stability of a set of ten candidate reference genes identified from the O. fasciatus transcriptome, and did so upon exposure to a dietary toxin, a cardiac glycoside, and across four different exposure periods. To increase robustness against gDNA contaminants, genome resources were used to design intron-bridging primers. A comprehensive stability validation by the Bestkeeper, Normfinder, geNorm and comparative ΔCt methods identified ef1a and tubulin as the most stable genes across treatments and time points, whereas 18S rRNA was the most unstable. However, accounting for the temporal scale indicated that time point confined normalizers might enable higher quantification accuracy for treatment comparison. Overall this study demonstrates: (i) a robust RT-qPCR primer design approach is possible for non-model organisms where genome annotation is often incomplete, and (ii) the importance of detailed reference gene stability exploration in multifactorial experimental designs.