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

Molecular evolution of a widely-adopted taxonomic marker (COI) across the animal tree of life

DNA barcodes are widely used for identification and discovery of species. While such use draws on information at the DNA level, the current amassment of ca. 4.7 million COI barcodes also offers a unique resource for exploring functional constraints on DNA evolution. Here, we explore amino acid variation in a crosscut of the entire animal kingdom. Patterns of DNA variation were linked to functional constraints at the level of the amino acid sequence in functionally important parts of the enzyme. Six amino acid sites show variation with possible effects on enzyme function. Overall, patterns of amino acid variation suggest convergent or parallel evolution at the protein level connected to the transition into a parasitic life style. Denser sampling of two diverse insect taxa revealed that the beetles (Coleoptera) show more amino acid variation than the butterflies and moths (Lepidoptera), indicating fundamental difference in patterns of molecular evolution in COI. Several amino acid sites were found to be under notably strong purifying selection in Lepidoptera as compared to Coleoptera. Overall, these findings demonstrate the utility of the global DNA barcode library to extend far beyond identification and taxonomy, and will hopefully be followed by a multitude of work.

Variable wing venation in Agathiphaga (Lepidoptera: Agathiphagidae) is key to understanding the evolution of basal moths

Details of the ancestral groundplan of wing venation in moths remain uncertain, despite approximately a century of study. Here, we describe a 3-branched subcostal vein, a 5-branched medial vein and a 2-branched cubitus posterior vein on the forewing of Agathiphaga vitiensis Dumbleton 1952 from Vanuatu. Such veins had not previously been described in any Lepidoptera. Because wing veins are typically lost during lepidopteran evolutionary history, rarely-if ever-to be regained, the venation of A. vitiensis probably represents the ancestral character state for moths. Wing venation is often used to identify fossil insects as moths, because wing scales are not always preserved; the presence of a supposedly trichopteran 3-branched subcostal vein in crown Lepidoptera may decrease the certainty with which certain amphiesmenopteran fossils from the Mesozoic can be classified. And because plesiomorphic veins can influence the development of lepidopteran wing patterns even if not expressed in the adult wing, the veins described here may determine the location of wing pattern elements in many lepidopteran taxa.

Genomic Mining of Phylogenetically Informative Nuclear Markers in Bark and Ambrosia Beetles

Deep level insect relationships are generally difficult to resolve, especially within taxa of the most diverse and species rich holometabolous orders. In beetles, the major diversity occurs in the Phytophaga, including charismatic groups such as leaf beetles, longhorn beetles and weevils. Bark and ambrosia beetles are wood boring weevils that contribute 12 percent of the diversity encountered in Curculionidae, one of the largest families of beetles with more than 50000 described species. Phylogenetic resolution in groups of Cretaceous age has proven particularly difficult and requires large quantity of data. In this study, we investigated 100 nuclear genes in order to select a number of markers with low evolutionary rates and high phylogenetic signal. A PCR screening using degenerate primers was applied to 26 different weevil species. We obtained sequences from 57 of the 100 targeted genes. Sequences from each nuclear marker were aligned and examined for detecting multiple copies, pseudogenes and introns. Phylogenetic informativeness (PI) and the capacity for reconstruction of previously established phylogenetic relationships were used as proxies for selecting a subset of the 57 amplified genes. Finally, we selected 16 markers suitable for large-scale phylogenetics of Scolytinae and related weevil taxa.

A Comparative Analysis of Sonic Defences in Bombycoidea Caterpillars

Caterpillars have long been used as models for studying animal defence. Their impressive armour, including flamboyant warning colours, poisonous spines, irritating sprays, and mimicry of plant parts, snakes and bird droppings, has been extensively documented. But research has mainly focused on visual and chemical displays. Here we show that some caterpillars also exhibit sonic displays. During simulated attacks, 45% of 38 genera and 33% of 61 species of silk and hawkmoth caterpillars (Bombycoidea) produced sounds. Sonic caterpillars are found in many distantly-related groups of Bombycoidea, and have evolved four distinct sound types- clicks, chirps, whistles and vocalizations. We propose that different sounds convey different messages, with some designed to warn of a chemical defence and others, to startle predators. This research underscores the importance of exploring acoustic communication in juvenile insects, and provides a model system to explore how different signals have evolved to frighten, warn or even trick predators.

Progress, pitfalls and parallel universes: a history of insect phylogenetics

The phylogeny of insects has been both extensively studied and vigorously debated for over a century. A relatively accurate deep phylogeny had been produced by 1904. It was not substantially improved in topology until recently when phylogenomics settled many long-standing controversies. Intervening advances came instead through methodological improvement. Early molecular phylogenetic studies (1985-2005), dominated by a few genes, provided datasets that were too small to resolve controversial phylogenetic problems. Adding to the lack of consensus, this period was characterized by a polarization of philosophies, with individuals belonging to either parsimony or maximum-likelihood camps; each largely ignoring the insights of the other. The result was an unfortunate detour in which the few perceived phylogenetic revolutions published by both sides of the philosophical divide were probably erroneous. The size of datasets has been growing exponentially since the mid-1980s accompanied by a wave of confidence that all relationships will soon be known. However, large datasets create new challenges, and a large number of genes does not guarantee reliable results. If history is a guide, then the quality of conclusions will be determined by an improved understanding of both molecular and morphological evolution, and not simply the number of genes analysed.

Preserving and vouchering butterflies and moths for large-scale museum-based molecular research

Butterflies and moths (Lepidoptera) comprise significant portions of the world’s natural history collections, but a standardized tissue preservation protocol for molecular research is largely lacking. Lepidoptera have traditionally been spread on mounting boards to display wing patterns and colors, which are often important for species identification. Many molecular phylogenetic studies have used legs from pinned specimens as the primary source for DNA in order to preserve a morphological voucher, but the amount of available tissue is often limited. Preserving an entire specimen in a cryogenic freezer is ideal for DNA preservation, but without an easily accessible voucher it can make specimen identification, verification, and morphological work difficult. Here we present a procedure that creates accessible and easily visualized “wing vouchers” of individual Lepidoptera specimens, and preserves the remainder of the insect in a cryogenic freezer for molecular research. Wings are preserved in protective holders so that both dorsal and ventral patterns and colors can be easily viewed without further damage. Our wing vouchering system has been implemented at the University of Maryland (AToL Lep Collection) and the University of Florida (Florida Museum of Natural History, McGuire Center of Lepidoptera and Biodiversity), which are among two of the largest Lepidoptera molecular collections in the world.

Species-Level Para- and Polyphyly in DNA Barcode Gene Trees: Strong Operational Bias in European Lepidoptera

The proliferation of DNA data is revolutionizing all fields of systematic research. DNA barcode sequences, now available for millions of specimens and several hundred thousand species, are increasingly used in algorithmic species delimitations. This is complicated by occasional incongruences between species and gene genealogies, as indicated by situations where conspecific individuals do not form a monophyletic cluster in a gene tree. In two previous reviews, non-monophyly has been reported as being common in mitochondrial DNA gene trees. We developed a novel web service “Monophylizer” to detect non-monophyly in phylogenetic trees and used it to ascertain the incidence of species non-monophyly in COI (a.k.a. cox1) barcode sequence data from 4977 species and 41,583 specimens of European Lepidoptera, the largest data set of DNA barcodes analyzed from this regard. Particular attention was paid to accurate species identification to ensure data integrity. We investigated the effects of tree-building method, sampling effort, and other methodological issues, all of which can influence estimates of non-monophyly. We found a 12% incidence of non-monophyly, a value significantly lower than that observed in previous studies. Neighbor joining (NJ) and maximum likelihood (ML) methods yielded almost equal numbers of non-monophyletic species, but 24.1% of these cases of non-monophyly were only found by one of these methods. Non-monophyletic species tend to show either low genetic distances to their nearest neighbors or exceptionally high levels of intraspecific variability. Cases of polyphyly in COI trees arising as a result of deep intraspecific divergence are negligible, as the detected cases reflected misidentifications or methodological errors. Taking into consideration variation in sampling effort, we estimate that the true incidence of non-monophyly is ∼23%, but with operational factors still being included. Within the operational factors, we separately assessed the frequency of taxonomic limitations (presence of overlooked cryptic and oversplit species) and identification uncertainties. We observed that operational factors are potentially present in more than half (58.6%) of the detected cases of non-monophyly. Furthermore, we observed that in about 20% of non-monophyletic species and entangled species, the lineages involved are either allopatric or parapatric—conditions where species delimitation is inherently subjective and particularly dependent on the species concept that has been adopted. These observations suggest that species-level non-monophyly in COI gene trees is less common than previously supposed, with many cases reflecting misidentifications, the subjectivity of species delimitation or other operational factors.

PCR primers for 30 novel gene regions in the nuclear genomes of Lepidoptera

We report primer pairs for 30 new gene regions in the nuclear genomes of Lepidoptera that can be amplified using a standard PCR protocol. The new primers were tested across diverse Lepidoptera, including nonditrysians and a wide selection of ditrysians. These new gene regions give a total of 11,043 bp of DNA sequence data and they show similar variability to traditionally used nuclear gene regions in studies of Lepidoptera. We feel that a PCR-based approach still has its place in molecular systematic studies of Lepidoptera, particularly at the intrafamilial level, and our new set of primers now provides a route to generating phylogenomic datasets using traditional methods.

Evidence for common horizontal transmission of Wolbachia among butterflies and moths

Background

Wolbachia is one of the most widespread bacteria on Earth. Previous research on Wolbachia-host interactions indicates that the bacterium is typically transferred vertically, from mother to offspring, through the egg cytoplasm. Although horizontal transmission of Wolbachia from one species to another is reported to be common in arthropods, limited direct ecological evidence is available. In this study, we examine horizontal transmission of Wolbachia using a multilocus sequence typing (MLST) strains dataset and used Wolbachia and Lepidoptera genomes to search for evidence for lateral gene transfer (LGT) in Lepidoptera, one of the most diverse cosmopolitan insect orders. We constructed a phylogeny of arthropod-associated MLST Wolbachia strains and calibrated the age of Wolbachia strains associated with lepidopteran species.

Results

Our results reveal inter-specific, inter-generic, inter-familial, and inter-ordinal horizontal transmission of Wolbachia strains, without discernible geographic patterns. We found at least seven probable cases of horizontal transmission among 31 species within Lepidoptera and between Lepidoptera and other arthropod hosts. The divergence time analysis revealed that Wolbachia is recently (22.6–4.7 mya, 95 % HPD) introduced in Lepidoptera. Analysis of nine Lepidoptera genomes (Bombyx mori, Danaus plexippus, Heliconius melpomene, Manduca sexta, Melitaea cinxia, Papilio glaucus, P. polytes, P. xuthus and Plutella xylostella) yielded one possible instance of Wolbachia LGT.

Conclusions

Our results provide evidence of high incidence of identical and multiple strains of Wolbachia among butterflies and moths, adding Lepidoptera to the growing body of evidence for common horizontal transmission of Wolbachia. This study demonstrates interesting dynamics of this remarkable and influential microorganism.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0660-x) contains supplementary material, which is available to authorized users.

Forewing color pattern in Micropterigidae (Insecta: Lepidoptera): homologies between contrast boundaries, and a revised hypothesis for the origin of symmetry systems

Background

Despite the great importance of lepidopteran wing patterns in various biological disciplines, homologies between wing pattern elements in different moth and butterfly lineages are still not understood. Among other reasons, this may be due to an incomplete understanding of the relationship between color pattern and wing venation; many individual wing pattern elements have a known relationship with venation, but a framework to unite all wing pattern elements with venation is lacking. Though plesiomorphic wing veins are known to influence color patterning even when not expressed in the adult wing, most studies of wing pattern evolution have focused on derived taxa with a reduced suite of wing veins.

Results

The present study aims to address this gap through an examination of Micropterigidae, a very early-diverged moth family in which all known plesiomorphic lepidopteran veins are expressed in the adult wing. The relationship between wing pattern and venation was examined in 66 species belonging to 9 genera. The relationship between venation and pattern element location, predicted based on moths in the family Tortricidae, holds for Sabatinca just as it does for Micropterix. However, the pattern elements that are lightly colored in Micropterix are dark in Sabatinca, and vice-versa. When plotted onto a hypothetical nymphalid wing in accordance with the relationship between pattern and venation discussed here, the wing pattern of Sabatinca doroxena very closely resembles the nymphalid groundplan.

Conclusions

The color difference in pattern elements between Micropterix and Sabatinca indicates that homologies exist among the contrast boundaries that divide wing pattern elements, and that color itself is not a reliable indicator of homology. The similarity between the wing pattern of Sabatinca doroxena and the nymphalid groundplan suggests that the nymphalid groundplan may have originated from a Sabatinca-like wing pattern subjected to changes in wing shape and reduced expression of venation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0687-z) contains supplementary material, which is available to authorized users.

Power, resolution and bias: recent advances in insect phylogeny driven by the genomic revolution

Our understanding on the phylogenetic relationships of insects has been revolutionised in the last decade by the proliferation of next generation sequencing technologies (NGS). NGS has allowed insect systematists to assemble very large molecular datasets that include both model and non-model organisms. Such datasets often include a large proportion of the total number of protein coding sequences available for phylogenetic comparison. We review some early entomological phylogenomic studies that employ a range of different data sampling protocols and analyses strategies, illustrating a fundamental renaissance in our understanding of insect evolution all driven by the genomic revolution. The analysis of phylogenomic datasets is challenging because of their size and complexity, and it is obvious that the increasing size alone does not ensure that phylogenetic signal overcomes systematic biases in the data. Biases can be due to various factors such as the method of data generation and assembly, or intrinsic biological feature of the data per se, such as similarities due to saturation or compositional heterogeneity. Such biases often cause violations in the underlying assumptions of phylogenetic models. We review some of the bioinformatics tools available and being developed to detect and minimise systematic biases in phylogenomic datasets. Phylogenomic-scale data coupled with sophisticated analyses will revolutionise our understanding of insect functional genomics. This will illuminate the relationship between the vast range of insect phenotypic diversity and underlying genetic diversity. In combination with rapidly developing methods to estimate divergence times, these analyses will also provide a compelling view of the rates and patterns of lineagenesis (birth of lineages) over the half billion years of insect evolution.

The Complete Mitochondrial Genome of Brachmia macroscopa (Lepidoptera: Gelechiidae) and Its Related Phylogenetic Analysis

The sweet potato leaf folder, Brachmia macroscopa, is an important pest in China. The complete mitogenome, which consists of 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and an A + T-rich region, was sequenced and found to be 15,394 bp in length (GeneBank no. KT354968). The gene order and orientation of the B. macroscopa mitogenome were similar to those of other sequenced lepidopteran species. All of the PCGs started with ATN as the canonical start codon except for cox1, which started with CGA. In regard to stop codons, most PCGs stopped at TAA except for cox2, which stopped at TA, and nad4, which stopped at a single T. Thirteen PCGs of the available species (33 taxa) were used to demonstrate phylogenetic relationships. The ditrysian cluster was supported as a monophyletic clade at high levels by using maximum likelihood and Bayesian methods. The apoditrysian group, covering the Gelechioidea, formed a monophyletic clade with a bootstrap value of 88% and a posterior probability of 1.00. The superfamily Gelechioidea was supported as a monophyletic lineage by a posterior probability of 1.00.

Phylogenetic analyses suggest that diversification and body size evolution are independent in insects

BACKGROUND

Skewed body size distributions and the high relative richness of small-bodied taxa are a fundamental property of a wide range of animal clades. The evolutionary processes responsible for generating these distributions are well described in vertebrate model systems but have yet to be explored in detail for other major terrestrial clades. In this study, we explore the macro-evolutionary patterns of body size variation across families of Hexapoda (insects and their close relatives), using recent advances in phylogenetic understanding, with an aim to investigate the link between size and diversity within this ancient and highly diverse lineage.

RESULTS

The maximum, minimum and mean-log body lengths of hexapod families are all approximately log-normally distributed, consistent with previous studies at lower taxonomic levels, and contrasting with skewed distributions typical of vertebrate groups. After taking phylogeny and within-tip variation into account, we find no evidence for a negative relationship between diversification rate and body size, suggesting decoupling of the forces controlling these two traits. Likelihood-based modeling of the log-mean body size identifies distinct processes operating within Holometabola and Diptera compared with other hexapod groups, consistent with accelerating rates of size evolution within these clades, while as a whole, hexapod body size evolution is found to be dominated by neutral processes including significant phylogenetic conservatism.

CONCLUSIONS

Based on our findings we suggest that the use of models derived from well-studied but atypical clades, such as vertebrates may lead to misleading conclusions when applied to other major terrestrial lineages. Our results indicate that within hexapods, and within the limits of current systematic and phylogenetic knowledge, insect diversification is generally unfettered by size-biased macro-evolutionary processes, and that these processes over large timescales tend to converge on apparently neutral evolutionary processes. We also identify limitations on available data within the clade and modeling approaches for the resolution of trees of higher taxa, the resolution of which may collectively enhance our understanding of this key component of terrestrial ecosystems.

Elusive ditrysian phylogeny: an account of combining systematized morphology with molecular data (Lepidoptera)

BACKGROUND

Ditrysia comprise close to 99 % of all butterflies and moths. The evolutionary relationships among the ditrysian superfamilies have received considerable attention in phylogenetic studies based on DNA and transcriptomic data, but the deepest divergences remain for large parts unresolved or contradictory. To obtain complementary insight into the evolutionary history of the clade, and to test previous hypotheses on the subdivision of Ditrysia based on morphology, we examine the morphology of larvae, pupae and adult males and females of 318 taxa representing nearly all ditrysian superfamilies and families. We present the most comprehensive morphological dataset on Ditrysia to date, consisting of over 500 morphological characters. The data are analyzed alone and combined with sequence data (one mitochondrial and seven nuclear protein-coding gene regions, sequenced from 422 taxa). The full dataset consists of 473 exemplar species. Analyses are performed using maximum likelihood methods, and parsimony methods for the morphological dataset. We explore whether combining morphological data and DNA-data can stabilize taxa that are unstable in phylogenetic studies based on genetic data only.

RESULTS

Morphological characters are found phylogenetically informative in resolving apical nodes (superfamilies and families), but characters serving as evidence of relatedness of larger assemblages are few. Results include the recovery of a monophyletic Tineoidea, Sesioidea and Cossoidea, and a stable position for some unstable taxa (e.g. Epipyropidae, Cyclotornidae, Urodoidea + Schreckensteinioidea). Several such taxa, however, remain unstable even though morphological characters indicate a position in the tree (e.g. Immidae). Evidence supporting affinities between clades are suggested, e.g. a novel larval synapomorphy for Tineidae. We also propose the synonymy of Tineodidae with Alucitidae, syn. nov.

CONCLUSIONS

The large morphological dataset provides information on the diversity and distribution of morphological traits in Ditrysia, and can be used in future research on the evolution of these traits, in identification keys and in identification of fossil Lepidoptera. The "backbone" of the phylogeny for Ditrysia remains largely unresolved. As previously proposed as an explanation for the scarcity of molecular signal in resolving the deeper nodes, this may be due to the rapid radiation of Ditrysia in the Cretaceous.

Misdirection of dosage compensation underlies bidirectional sex-specific death in Wolbachia-infected Ostrinia scapulalis

Endosymbiotic bacteria of the genus Wolbachia often manipulate the reproductive system of their hosts to propagate themselves in host populations. Ostrinia scapulalis moths infected with Wolbachia (wSca) produce female-only progeny (sex chromosomes: ZW), whereas females cured of the infection by antibiotic treatment produce male-only progeny (ZZ). The occurrence of female- and male-only progeny has been attributed to the specific death of the opposite sex during embryonic and larval development. In this bidirectional sex-specific lethality, embryos destined to die express a phenotypic sex opposite to their genotypic sex. On the basis of these findings, we suggested that wSca carries a genetic factor that feminizes the male host, the W chromosome of the host has lost its feminizing function, and discordance between the genotypic and phenotypic sexes underlies this sex-specific death. In the present study, we examined whether the failure of dosage compensation was responsible for this sex-specific mortality. Quantitative PCRs showed that Z-linked gene expression levels in embryos destined to die were not properly dosage compensated; they were approximately two-fold higher in the male progeny of wSca-infected females and approximately two-fold lower in the female progeny of infected-and-cured females. These results support our hypothesis that misdirection of dosage compensation underlies the sex-specific death.

Color Pattern on the Forewing of Micropterix (Lepidoptera: Micropterigidae): Insights into the Evolution of Wing Pattern and Wing Venation in Moths

Wing patterns are key taxonomic characters that have long been used in descriptions of Lepidoptera; however, wing pattern homologies are not understood among different moth lineages. Here, we examine the relationship between wing venation and wing pattern in the genus Micropterix, among the most basal extant Lepidoptera, in order to evaluate the two existing predictive models that have the potential to establish wing pattern element homologies for the order. The location of wing pattern elements along the costal margin of the wing in Micropterix is consistent with the predictions of the model proposed for Tortricidae by Brown and Powell in 1991, later modified by Baixeras in 2002. The predictive power of this model for such distantly related taxa suggests that the model may hold across various superfamilies within Lepidoptera, and supports the long-held notion that fasciae, not spots, are the most likely primitive wing pattern elements for the order. In addition, the location of wing pattern elements suggests that the wing vein commonly termed Sc1 may in fact be a different vein, which Comstock identified in Trichoptera and referred to as "a."

Sex Chromosome Dosage Compensation in Heliconius Butterflies: Global yet Still Incomplete?

The evolution of heterogametic sex chromosomes is often—but not always—accompanied by the evolution of dosage compensating mechanisms that mitigate the impact of sex-specific gene dosage on levels of gene expression. One emerging view of this process is that such mechanisms may only evolve in male-heterogametic (XY) species but not in female-heterogametic (ZW) species, which will consequently exhibit “incomplete” sex chromosome dosage compensation. However, recent results suggest that at least some Lepidoptera (moths and butterflies) may prove to be an exception to this prediction. Studies in bombycoid moths indicate the presence of a chromosome-wide epigenetic mechanism that effectively balances Z chromosome gene expression between the sexes by reducing Z-linked expression in males. In contrast, strong sex chromosome dosage effects without any reduction in male Z-linked expression were previously reported in a pyralid moth, suggesting a lack of any such dosage compensating mechanism. Here we report an analysis of sex chromosome dosage compensation in Heliconius butterflies, sampling multiple individuals for several different adult tissues (head, abdomen, leg, mouth, and antennae). Methodologically, we introduce a novel application of linear mixed-effects models to assess dosage compensation, offering a unified statistical framework that can estimate effects specific to chromosome, to sex, and their interactions (i.e., a dosage effect). Our results show substantially reduced Z-linked expression relative to autosomes in both sexes, as previously observed in bombycoid moths. This observation is consistent with an increasing body of evidence that some lepidopteran species possess an epigenetic dosage compensating mechanism that reduces Z chromosome expression in males to levels comparable with females. However, this mechanism appears to be imperfect in Heliconius, resulting in a modest dosage effect that produces an average 5–20% increase in male expression relative to females on the Z chromosome, depending on the tissue. Thus our results in Heliconius reflect a mixture of previous patterns reported for Lepidoptera. In Heliconius, a moderate pattern of incomplete dosage compensation persists apparently despite the presence of an epigenetic dosage compensating mechanism. The chromosomal distributions of sex-biased genes show an excess of male-biased and a dearth of female-biased genes on the Z chromosome relative to autosomes, consistent with predictions of sexually antagonistic evolution.

Relationships within Cladobranchia (Gastropoda: Nudibranchia) based on RNA-Seq data: an initial investigation

Cladobranchia (Gastropoda: Nudibranchia) is a diverse (approx. 1000 species) but understudied group of sea slug molluscs. In order to fully comprehend the diversity of nudibranchs and the evolution of character traits within Cladobranchia, a solid understanding of evolutionary relationships is necessary. To date, only two direct attempts have been made to understand the evolutionary relationships within Cladobranchia, neither of which resulted in well-supported phylogenetic hypotheses. In addition to these studies, several others have addressed some of the relationships within this clade while investigating the evolutionary history of more inclusive groups (Nudibranchia and Euthyneura). However, all of the resulting phylogenetic hypotheses contain conflicting topologies within Cladobranchia. In this study, we address some of these long-standing issues regarding the evolutionary history of Cladobranchia using RNA-Seq data (transcriptomes). We sequenced 16 transcriptomes and combined these with four transcriptomes from the NCBI Sequence Read Archive. Transcript assembly using Trinity and orthology determination using HaMStR yielded 839 orthologous groups for analysis. These data provide a well-supported and almost fully resolved phylogenetic hypothesis for Cladobranchia. Our results support the monophyly of Cladobranchia and the sub-clade Aeolidida, but reject the monophyly of Dendronotida.

A FISH-based chromosome map for the European corn borer yields insights into ancient chromosomal fusions in the silkworm

A significant feature of the genomes of Lepidoptera, butterflies and moths, is the high conservation of chromosome organization. Recent remarkable progress in genome sequencing of Lepidoptera has revealed that syntenic gene order is extensively conserved across phylogenetically distant species. The ancestral karyotype of Lepidoptera is thought to be n=31; however, that of the most well-studied moth, Bombyx mori, is n=28, and diverse studies suggest that three chromosomal fusion events occurred in this lineage. To identify the boundaries between predicted ancient fusions involving B. mori chromosomes 11, 23 and 24, we constructed fluorescence in situ hybridization (FISH)-based chromosome maps of the European corn borer, Ostrinia nubilalis (n=31). We first determined a 511 Mb genomic sequence of the Asian corn borer, O. furnacalis, a congener of O. nubilalis, and isolated bacterial artificial chromosomes and fosmid clones that were expected to localize in candidate regions for the boundaries using these sequences. Combined with FISH and genetic analysis, we narrowed down the candidate regions to 40 kb-1.5 Mb, in strong agreement with a previous estimate based on the genome of a butterfly, Melitaea cinxia. The significant difference in the lengths of the candidate regions where no functional genes were observed may reflect the evolutionary time after fusion events.

Identification of the Female-Produced Sex Pheromone of the Leafminer Holocacista capensis Infesting Grapevine in South Africa

We report the first identification of a sex pheromone in a heliozelid moth, Holocacista capensis van Nieukerken & Geertsema. This leafminer recently infested grapevine in South Africa. Compared to solvent extraction of pheromone glands, solid phase microextraction (SPME) proved to be highly effective for collection of the pheromone from calling females. The volatiles collected by SPME were analyzed by gas chromatography with electroantennographic detection (GC/EAD). Three compounds eliciting electrophysiological activity from the male antenna were identified as (Z)-5-tetradecenal, (Z)-7-tetradecenal, and (Z)-9-hexadecenal by coupled gas chromatography-mass spectrometry (GC/MS). GC/MS analysis of dimethyldisulphide (DMDS) derivatives of fatty acyl moieties in pheromone gland extracts confirmed the presence of the corresponding putative pheromone precursors with double bonds in the same position and with Z geometry. Field trapping experiments in a South African vineyard confirmed that both (Z)-5-tetradecenal and (Z)-7-tetradecenal are essential for the attraction of male H. capensis, whereas addition of (Z)-9-hexadecenal to the blend did not affect the attractiveness. The composition of the pheromone is discussed in relation to the phylogeny of this family of moths.

Wolbachia in butterflies and moths: geographic structure in infection frequency

INTRODUCTION

Butterflies and moths (Lepidoptera) constitute one of the most diverse insect orders, and play an important role in ecosystem function. However, little is known in terms of their bacterial communities. Wolbachia, perhaps the most common and widespread intracellular bacterium on Earth, can manipulate the physiology and reproduction of its hosts, and is transmitted vertically from mother to offspring, or sometimes horizontally between species. While its role in some hosts has been studied extensively, its incidence across Lepidoptera is poorly understood. A recent analysis using a beta-binomial model to infer the between-species distribution of prevalence estimated that approximately 40 % of arthropod species are infected with Wolbachia, but particular taxonomic groups and ecological niches seem to display substantially higher or lower incidences. In this study, we took an initial step and applied a similar, maximum likelihood approach to 300 species of Lepidoptera (7604 individuals from 660 populations) belonging to 17 families and 10 superfamilies, and sampled from 36 countries, representing all continents excluding Antarctica.

RESULTS

Approximately a quarter to a third of individuals appear to be infected with Wolbachia, and around 80 % of Lepidoptera species are infected at a non-negligible frequency. This incidence estimate is very high compared to arthropods in general. Wolbachia infection in Lepidoptera is shown to vary between families, but there is no evidence for closely related groups to show similar infection levels. True butterflies (Papilionoidea) are overrepresented in our data, however, our estimates show this group can be taken as a representative for the other major lepidopteran superfamilies. We also show substantial variation in infection level according to geography - closer locations tend to show similar infection levels. We further show that variation in geography is due to a latitudinal gradient in Wolbachia infection, with lower frequencies towards higher latitudes.

CONCLUSIONS

Our comprehensive survey of Wolbachia infection in Lepidoptera suggests that infection incidence is very high, and provides evidence that climate and geography are strong predictors of infection frequency.

The complete mitochondrial genome of Papilio glaucus and its phylogenetic implications

Due to the intriguing morphology, lifecycle, and diversity of butterflies and moths, Lepidoptera are emerging as model organisms for the study of genetics, evolution and speciation. The progress of these studies relies on decoding Lepidoptera genomes, both nuclear and mitochondrial. Here we describe a protocol to obtain mitogenomes from Next Generation Sequencing reads performed for whole-genome sequencing and report the complete mitogenome of Papilio (Pterourus) glaucus. The circular mitogenome is 15,306 bp in length and rich in A and T. It contains 13 protein-coding genes (PCGs), 22 transfer-RNA-coding genes (tRNA), and 2 ribosomal-RNA-coding genes (rRNA), with a gene order typical for mitogenomes of Lepidoptera. We performed phylogenetic analyses based on PCG and RNA-coding genes or protein sequences using Bayesian Inference and Maximum Likelihood methods. The phylogenetic trees consistently show that among species with available mitogenomes Papilio glaucus is the closest to Papilio (Agehana) maraho from Asia.

Phylogenetic relationships of subfamilies in the family Hesperiidae (Lepidoptera: Hesperioidea) from China

Hesperiidae is one of the largest families of butterflies. Our knowledge of the higher systematics on hesperiids from China is still very limited. We infer the phylogenetic relationships of the subfamilies of Chinese skippers based on three mitochondrial genes (cytochrome b (Cytb), the NADH dehydrogenase subunit 1 (ND1) and cytochrome oxidase I (COI)). In this study, 30 species in 23 genera were included in the Bayesian and maximum likelihood analyses. The subfamily Coeliadinae, Eudaminae, Pyrginae and Heteropterinae were recovered as a monophyletic clade with strong support. The subfamily Hesperiinae formed a clade, but support for monophyly was weak. Our results imply that the five subfamilies of Chinese Hesperiidae should be divided into: Coeliadinae, Eudaminae, Pyrginae, Heteropterinae and Hesperiinae. The relationships of the five subfamilies should be as follows: Coeliadinae + (Eudaminae + (Pyrginae + (Heteropterinae + Hesperiinae))).

Tempo and mode of antibat ultrasound production and sonar jamming in the diverse hawkmoth radiation

The bat-moth arms race has existed for over 60 million y, with moths evolving ultrasonically sensitive ears and ultrasound-producing organs to combat bat predation. The evolution of these defenses has never been thoroughly examined because of limitations in simultaneously conducting behavioral and phylogenetic analyses across an entire group. Hawkmoths include >1,500 species worldwide, some of which produce ultrasound using genital stridulatory structures. However, the function and evolution of this behavior remain largely unknown. We built a comprehensive behavioral dataset of hawkmoth hearing and ultrasonic reply to sonar attack using high-throughput field assays. Nearly half of the species tested (57 of 124 species) produced ultrasound to tactile stimulation or playback of bat echolocation attack. To test the function of ultrasound, we pitted big brown bats (Eptesicus fuscus) against hawkmoths over multiple nights and show that hawkmoths jam bat sonar. Ultrasound production was immediately and consistently effective at thwarting attack and bats regularly performed catching behavior without capturing moths. We also constructed a fossil-calibrated, multigene phylogeny to study the evolutionary history and divergence times of these antibat strategies across the entire family. We show that ultrasound production arose in multiple groups, starting in the late Oligocene (∼ 26 Ma) after the emergence of insectivorous bats. Sonar jamming and bat-detecting ears arose twice, independently, in the Miocene (18-14 Ma) either from earless hawkmoths that produced ultrasound in response to physical contact only, or from species that did not respond to touch or bat echolocation attack.

Brain composition in Godyris zavaleta, a diurnal butterfly, Reflects an increased reliance on olfactory information

Interspecific comparisons of brain structure can inform our functional understanding of brain regions, identify adaptations to species-specific ecologies, and explore what constrains adaptive changes in brain structure, and coevolution between functionally related structures. The value of such comparisons is enhanced when the species considered have known ecological differences. The Lepidoptera have long been a favored model in evolutionary biology, but to date descriptions of brain anatomy have largely focused on a few commonly used neurobiological model species. We describe the brain of Godyris zavaleta (Ithomiinae), a member of a subfamily of Neotropical butterflies with enhanced reliance on olfactory information. We demonstrate for the first time the presence of sexually dimorphic glomeruli within a distinct macroglomerular complex (MGC) in the antennal lobe of a diurnal butterfly. This presents a striking convergence with the well-known moth MGC, prompting a discussion of the potential mechanisms behind the independent evolution of specialized glomeruli. Interspecific analyses across four Lepidoptera further show that the relative size of sensory neuropils closely mirror interspecific variation in sensory ecology, with G. zavaleta displaying levels of sensory investment intermediate between the diurnal monarch butterfly (Danaus plexippus), which invests heavily in visual neuropil, and night-flying moths, which invest more in olfactory neuropil. We identify several traits that distinguish butterflies from moths, and several that distinguish D. plexippus and G. zavaleta. Our results illustrate that ecological selection pressures mold the structure of invertebrate brains, and exemplify how comparative analyses across ecologically divergent species can illuminate the functional significance of variation in brain structure.

Explosive radiation or uninformative genes? Origin and early diversification of tachinid flies (Diptera: Tachinidae)

Molecular phylogenetic studies at all taxonomic levels often infer rapid radiation events based on short, poorly resolved internodes. While such rapid episodes of diversification are an important and widespread evolutionary phenomenon, much of this poor phylogenetic resolution may be attributed to the continuing widespread use of "traditional" markers (mitochondrial, ribosomal, and some nuclear protein-coding genes) that are often poorly suited to resolve difficult, higher-level phylogenetic problems. Here we reconstruct phylogenetic relationships among a representative set of taxa of the parasitoid fly family Tachinidae and related outgroups of the superfamily Oestroidea. The Tachinidae are one of the most species rich, yet evolutionarily recent families of Diptera, providing an ideal case study for examining the differential performance of loci in resolving phylogenetic relationships and the benefits of adding more loci to phylogenetic analyses. We assess the phylogenetic utility of nine genes including both traditional genes (e.g., CO1 mtDNA, 28S rDNA) and nuclear protein-coding genes newly developed for phylogenetic analysis. Our phylogenetic findings, based on a limited set of taxa, include: a close relationship between Tachinidae and the calliphorid subfamily Polleninae, monophyly of Tachinidae and the subfamilies Exoristinae and Dexiinae, subfamily groupings of Dexiinae+Phasiinae and Tachininae+Exoristinae, and robust phylogenetic placement of the somewhat enigmatic genera Strongygaster, Euthera, and Ceracia. In contrast to poor resolution and phylogenetic incongruence of "traditional genes," we find that a more selective set of highly informative genes is able to more precisely identify regions of the phylogeny that experienced rapid radiation of lineages, while more accurately depicting their phylogenetic context. Although much expanded taxon sampling is necessary to effectively assess the monophyly of and relationships among major tachinid lineages and their relatives, we show that a small number of well-chosen nuclear protein-coding genes can successfully resolve even difficult phylogenetic problems.

Four new species of the genus Diolcogaster Ashmead, 1900 (Hymenoptera: Braconidae: Microgastrinae) from South East Asia with a key to the Indian species

Four new species of the genus Diolcogaster Ashmead, 1900 (Hymenoptera: Braconidae: Microgastrinae) are described and illustrated: Diolcogaster andamanensis n. sp. from the Andaman Islands, and D. duocolor n. sp., D. longistria n. sp. and D. solitarium n. sp. from mainland India. The solitary larval parasitoid D. solitarium was reared from Gatesclarkeana sp. (Lepidoptera: Tortricidae). A new combination, Diolcogaster tomentosae (Wilkinson, 1930) n. comb., is proposed for the Indian species Protomicroplitis tomentosae (Wilkinson, 1930) along with its redescription and documentation of the gregarious cocoons associated with the pyralid (Epipaschiinae) host feeding on Terminalia cattappa L.

The fossil record and taphonomy of butterflies and moths (Insecta, Lepidoptera): implications for evolutionary diversity and divergence-time estimates

Background

It is conventionally accepted that the lepidopteran fossil record is significantly incomplete when compared to the fossil records of other, very diverse, extant insect orders. Such an assumption, however, has been based on cumulative diversity data rather than using alternative statistical approaches from actual specimen counts.

Results

We reviewed documented specimens of the lepidopteran fossil record, currently consisting of 4,593 known specimens that are comprised of 4,262 body fossils and 331 trace fossils. The temporal distribution of the lepidopteran fossil record shows significant bias towards the late Paleocene to middle Eocene time interval. Lepidopteran fossils also record major shifts in preservational style and number of represented localities at the Mesozoic stage and Cenozoic epoch level of temporal resolution. Only 985 of the total known fossil specimens (21.4%) were assigned to 23 of the 40 extant lepidopteran superfamilies. Absolute numbers and proportions of preservation types for identified fossils varied significantly across superfamilies. The secular increase of lepidopteran family-level diversity through geologic time significantly deviates from the general pattern of other hyperdiverse, ordinal-level lineages.

Conclusion

Our statistical analyses of the lepidopteran fossil record show extreme biases in preservation type, age, and taxonomic composition. We highlight the scarcity of identified lepidopteran fossils and provide a correspondence between the latest lepidopteran divergence-time estimates and relevant fossil occurrences at the superfamily level. These findings provide caution in interpreting the lepidopteran fossil record through the modeling of evolutionary diversification and in determination of divergence time estimates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0290-8) contains supplementary material, which is available to authorized users.

A Burst of miRNA Innovation in the Early Evolution of Butterflies and Moths

MicroRNAs (miRNAs) are involved in posttranscriptional regulation of gene expression. Because several miRNAs are known to affect the stability or translation of developmental regulatory genes, the origin of novel miRNAs may have contributed to the evolution of developmental processes and morphology. Lepidoptera (butterflies and moths) is a species-rich clade with a well-established phylogeny and abundant genomic resources, thereby representing an ideal system in which to study miRNA evolution. We sequenced small RNA libraries from developmental stages of two divergent lepidopterans, Cameraria ohridella (Horse chestnut Leafminer) and Pararge aegeria (Speckled Wood butterfly), discovering 90 and 81 conserved miRNAs, respectively, and many species-specific miRNA sequences. Mapping miRNAs onto the lepidopteran phylogeny reveals rapid miRNA turnover and an episode of miRNA fixation early in lepidopteran evolution, implying that miRNA acquisition accompanied the early radiation of the Lepidoptera. One lepidopteran-specific miRNA gene, miR-2768, is located within an intron of the homeobox gene invected, involved in insect segmental and wing patterning. We identified cubitus interruptus (ci) as a likely direct target of miR-2768, and validated this suppression using a luciferase assay system. We propose a model by which miR-2768 modulates expression of ci in the segmentation pathway and in patterning of lepidopteran wing primordia.

Phylotranscriptomics: saturated third codon positions radically influence the estimation of trees based on next-gen data

Recent advancements in molecular sequencing techniques have led to a surge in the number of phylogenetic studies that incorporate large amounts of genetic data. We test the assumption that analyzing large number of genes will lead to improvements in tree resolution and branch support using moths in the superfamily Bombycoidea, a group with some interfamilial relationships that have been difficult to resolve. Specifically, we use a next-gen data set that included 19 taxa and 938 genes (∼1.2M bp) to examine how codon position and saturation might influence resolution and node support among three key families. Maximum likelihood, parsimony, and species tree analysis using gene tree parsimony, on different nucleotide and amino acid data sets, resulted in largely congruent topologies with high bootstrap support compared with prior studies that included fewer loci. However, for a few shallow nodes, nucleotide and amino acid data provided high support for conflicting relationships. The third codon position was saturated and phylogenetic analysis of this position alone supported a completely different, potentially misleading sister group relationship. We used the program RADICAL to assess the number of genes needed to fix some of these difficult nodes. One such node originally needed a total of 850 genes but only required 250 when synonymous signal was removed. Our study shows that, in order to effectively use next-gen data to correctly resolve difficult phylogenetic relationships, it is necessary to assess the effects of synonymous substitutions and third codon positions.

Environmental determinism of community structure across trophic levels: moth assemblages and substrate type in the rain forests of south-western China

Soil type may drive vegetation structure. In turn, the richness, identity and diversity of arthropod herbivores may be related to plant diversity through specific host plant relationships in a location. We test the hypothesis that the soil type (calcicolous vs alluvial soils) will drive the assemblage structure of a dominant group of arthropod herbivores: the moths. We used sampling sites in rain-forest fragments in south-western China around the Xishuangbanna Tropical Botanical Gardens (21°41′N, 101°25′E) to test this hypothesis. We used Pennsylvania style light traps to take point samples of macromoths and pyraloids from four sampling sites in forest remnants on a limestone geological base and four from alluvial-based forest. A total of 3165 moths (1739 from limestone-based and 1255 from alluvium-based forests) was collected representing 1255 species. The limestone-based sites showed statistically similar levels of species richness and other alpha diversity indices to the four alluvium-based site. Nevertheless the sites were clearly significantly different in terms of species composition. Analysis of contrasting similarity (‘beta’ diversity) indices suggested that there was ‘leakage’ between the two classes of sites when ‘rare’ species were emphasized in the calculations. We used an indicator value procedure to select species that most characterized this separation. We expect that these differences reflect associated changes in plant assemblage structure acting through the herbivorous habits of larval moths. Accordingly, in any assessment of landscape level diversity the nature of the substrate and its associated vegetation is clearly of great importance. This observation also has consequences for the design of conservation programmes.

Phylogenetic distribution of extant richness suggests metamorphosis is a key innovation driving diversification in insects

Insects and their six-legged relatives (Hexapoda) comprise more than half of all described species and dominate terrestrial and freshwater ecosystems. Understanding the macroevolutionary processes generating this richness requires a historical perspective, but the fossil record of hexapods is patchy and incomplete. Dated molecular phylogenies provide an alternative perspective on divergence times and have been combined with birth-death models to infer patterns of diversification across a range of taxonomic groups. Here we generate a dated phylogeny of hexapod families, based on previously published sequence data and literature derived constraints, in order to identify the broad pattern of macroevolutionary changes responsible for the composition of the extant hexapod fauna. The most prominent increase in diversification identified is associated with the origin of complete metamorphosis, confirming this as a key innovation in promoting insect diversity. Subsequent reductions are recovered for several groups previously identified as having a higher fossil diversity during the Mesozoic. In addition, a number of recently derived taxa are found to have radiated following the development of flowering plant (angiosperm) floras during the mid-Cretaceous. These results reveal that the composition of the modern hexapod fauna is a product of a key developmental innovation, combined with multiple and varied evolutionary responses to environmental changes from the mid Cretaceous floral transition onward.

Ancient expansion of the hox cluster in lepidoptera generated four homeobox genes implicated in extra-embryonic tissue formation

Gene duplications within the conserved Hox cluster are rare in animal evolution, but in Lepidoptera an array of divergent Hox-related genes (Shx genes) has been reported between pb and zen. Here, we use genome sequencing of five lepidopteran species (Polygonia c-album, Pararge aegeria, Callimorpha dominula, Cameraria ohridella, Hepialus sylvina) plus a caddisfly outgroup (Glyphotaelius pellucidus) to trace the evolution of the lepidopteran Shx genes. We demonstrate that Shx genes originated by tandem duplication of zen early in the evolution of large clade Ditrysia; Shx are not found in a caddisfly and a member of the basally diverging Hepialidae (swift moths). Four distinct Shx genes were generated early in ditrysian evolution, and were stably retained in all descendent Lepidoptera except the silkmoth which has additional duplications. Despite extensive sequence divergence, molecular modelling indicates that all four Shx genes have the potential to encode stable homeodomains. The four Shx genes have distinct spatiotemporal expression patterns in early development of the Speckled Wood butterfly (Pararge aegeria), with ShxC demarcating the future sites of extraembryonic tissue formation via strikingly localised maternal RNA in the oocyte. All four genes are also expressed in presumptive serosal cells, prior to the onset of zen expression. Lepidopteran Shx genes represent an unusual example of Hox cluster expansion and integration of novel genes into ancient developmental regulatory networks.

Combover/CG10732, a novel PCP effector for Drosophila wing hair formation

The polarization of cells is essential for the proper functioning of most organs. Planar Cell Polarity (PCP), the polarization within the plane of an epithelium, is perpendicular to apical-basal polarity and established by the non-canonical Wnt/Fz-PCP signaling pathway. Within each tissue, downstream PCP effectors link the signal to tissue specific readouts such as stereocilia orientation in the inner ear and hair follicle orientation in vertebrates or the polarization of ommatidia and wing hairs in Drosophila melanogaster. Specific PCP effectors in the wing such as Multiple wing hairs (Mwh) and Rho Kinase (Rok) are required to position the hair at the correct position and to prevent ectopic actin hairs. In a genome-wide screen in vitro, we identified Combover (Cmb)/CG10732 as a novel Rho kinase substrate. Overexpression of Cmb causes the formation of a multiple hair cell phenotype (MHC), similar to loss of rok and mwh. This MHC phenotype is dominantly enhanced by removal of rok or of other members of the PCP effector gene family. Furthermore, we show that Cmb physically interacts with Mwh, and cmb null mutants suppress the MHC phenotype of mwh alleles. Our data indicate that Cmb is a novel PCP effector that promotes to wing hair formation, a function that is antagonized by Mwh.

The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera

Previous studies have reported that chromosome synteny in Lepidoptera has been well conserved, yet the number of haploid chromosomes varies widely from 5 to 223. Here we report the genome (393 Mb) of the Glanville fritillary butterfly (Melitaea cinxia; Nymphalidae), a widely recognized model species in metapopulation biology and eco-evolutionary research, which has the putative ancestral karyotype of n=31. Using a phylogenetic analyses of Nymphalidae and of other Lepidoptera, combined with orthologue-level comparisons of chromosomes, we conclude that the ancestral lepidopteran karyotype has been n=31 for at least 140 My. We show that fusion chromosomes have retained the ancestral chromosome segments and very few rearrangements have occurred across the fusion sites. The same, shortest ancestral chromosomes have independently participated in fusion events in species with smaller karyotypes. The short chromosomes have higher rearrangement rate than long ones. These characteristics highlight distinctive features of the evolutionary dynamics of butterflies and moths.

Butterflies and moths (Lepidoptera) vary in chromosome number. Here, the authors sequence the genome of the Glanville fritillary butterfly, Melitaea cinxia, show it has the ancestral lepidopteran karyotype and provide insight into how chromosomal fusions have shaped karyotype evolution in butterflies and moths.

Phylogenomics provides strong evidence for relationships of butterflies and moths

Butterflies and moths constitute some of the most popular and charismatic insects. Lepidoptera include approximately 160 000 described species, many of which are important model organisms. Previous studies on the evolution of Lepidoptera did not confidently place butterflies, and many relationships among superfamilies in the megadiverse clade Ditrysia remain largely uncertain. We generated a molecular dataset with 46 taxa, combining 33 new transcriptomes with 13 available genomes, transcriptomes and expressed sequence tags (ESTs). Using HaMStR with a Lepidoptera-specific core-orthologue set of single copy loci, we identified 2696 genes for inclusion into the phylogenomic analysis. Nucleotides and amino acids of the all-gene, all-taxon dataset yielded nearly identical, well-supported trees. Monophyly of butterflies (Papilionoidea) was strongly supported, and the group included skippers (Hesperiidae) and the enigmatic butterfly-moths (Hedylidae). Butterflies were placed sister to the remaining obtectomeran Lepidoptera, and the latter was grouped with greater than or equal to 87% bootstrap support. Establishing confident relationships among the four most diverse macroheteroceran superfamilies was previously challenging, but we recovered 100% bootstrap support for the following relationships: ((Geometroidea, Noctuoidea), (Bombycoidea, Lasiocampoidea)). We present the first robust, transcriptome-based tree of Lepidoptera that strongly contradicts historical placement of butterflies, and provide an evolutionary framework for genomic, developmental and ecological studies on this diverse insect order.

A Linnaeus NG (TM) interactive key to the Lithocolletinae of North-West Europe aimed at accelerating the accumulation of reliable biodiversity data (Lepidoptera, Gracillariidae)

We present an interactive key that is available online through any web browser without the need to install any additional software, making it an easily accessible tool for the larger public. The key can be found at http://identify.naturalis.nl/lithocolletinae. The key includes all 86 North-West European Lithocolletinae, a subfamily of smaller moths (“micro-moths”) that is commonly not treated in field guides. The user can input data on several external morphological character systems in addition to distribution, host plant and even characteristics of the larval feeding traces to reach an identification. We expect that this will enable more people to contribute with reliable observation data on this group of moths and alleviate the workload of taxonomic specialists, allowing them to focus on other new keys or taxonomic work.

Complete mitochondrial genomes of five skippers (Lepidoptera: Hesperiidae) and phylogenetic reconstruction of Lepidoptera

We sequenced mitogenomes of five skippers (family Hesperiidae, Lepidoptera) to obtain further insight into the characteristics of butterfly mitogenomes and performed phylogenetic reconstruction using all available gene sequences (PCGs, rRNAs, and tRNAs) from 85 species (20 families in eight superfamilies). The general genomic features found in the butterflies also were found in the five skippers: a high A+T composition (79.3%-80.9%), dominant usage of TAA stop codon, similar skewness pattern in both strands, consistently length intergenic spacer sequence between tRNA(Gln) and ND2 (64-87 bp), conserved ATACTAA motif between tRNA(Ser (UCN)) and ND1, and characteristic features of the A+T-rich region (the ATAGA motif, varying length of poly-T stretch, and poly-A stretch). The start codon for COI was CGA in four skippers as typical, but Lobocla bifasciatus evidently possessed canonical ATG as start codon. All species had the ancestral arrangement tRNA(Asn)/tRNA(Ser (AGN)), instead of the rearrangement tRNA(Ser (AGN))/tRNA(Asn), found in another skipper species (Erynnis). Phylogenetic analyses using all available genes (PCGs, rRNAS, and tRNAs) yielded the consensus superfamilial relationships ((((((Bombycoidea+Noctuoidea+Geometroidea)+Pyraloidea)+Papilionoidea)+Tortricoidea)+Yponomeutoidea)+Hepialoidea), confirming the validity of Macroheterocera (Bombycoidea, Noctuoidea, and Geometroidea in this study) and its sister relationship to Pyraloidea. Within Rhopalocera (butterflies and skippers) the familial relationships (Papilionidae+(Hesperiidae+(Pieridae+((Lycaenidae+Riodinidae)+Nymphalidae)))) were strongly supported in all analyses (0.98-1 by BI and 96-100 by ML methods), rendering invalid the superfamily status for Hesperioidea. On the other hand, current mitogenome-based phylogeny did not find consistent superfamilial relationships among Noctuoidea, Geometroidea, and Bombycoidea and the familial relationships within Bombycoidea between analyses, requiring further taxon sampling in future studies.