BACKGROUND
five years ago, the backbone branching pattern of the evolutionary tree for one of the largest groups of insects, moths and butterflies , was still for the most part unresolved.
our state of knowledge of phylogenetic affinities within this insect order, with over  <dig>  <dig> described species  <cit> , was well depicted by the comb-like summary tree composed by kristensen and skalski  <cit> .
well-supported hypotheses existed only for the relationships among the most ancient superfamilies.
only very cautious and speculative hypotheses had been proposed for the relationships among the superfamilies of the enormous clade of more advanced lepidopterans, known as the ditrysia  <cit> .
this ditrysian clade encompasses nearly  <dig> % of all butterflies and moths  <cit> .
ditrysia are currently divided into  <dig> superfamilies  .
the size of the superfamilies varies considerably in the number of described species, from several small monotypic superfamilies based on a single genus with only a single or a few species to noctuoidea with over  <dig>  <dig> species.
some of these monotypic superfamilies are characterized by features that do not correspond to those diagnosing any of the other superfamilies.
the current circumscriptions of the ditrysian superfamilies are based on morphological features and have been generally accepted.
recent studies based on genetic data have corroborated the monophyly of most of the superfamilies .
however, the interrelationships among most of these superfamilies lack any convincing supporting evidence.
this lack of knowledge of evolutionary affinities has hindered further research on this enormous suborder  <cit> .fig.
<dig> current classification of ditrysia, modified from van nieukerken et al.
<cit>
the initial challenges to inferring evolutionary relationships among the ditrysian superfamilies have been attributed to their morphological homogeneity, and thus to the difficulty of finding distinct, unique morphological characters uniting superfamilies and supporting descent from a common ancestor  <cit> .
untold numbers of detailed studies on the morphology of lepidoptera, especially of genital structures and wing venation, have furthered our understanding of the relatedness of species, genera and families, but the relationships among larger assemblages have remained obscure on morphological grounds.
the most important studies in comparative morphology aimed at resolving the relationships among ditrysian lepidoptera were by brock  <cit>  and minet  <cit> .
the works by minet have been influential and the basis for the division of ditrysia into three nested clades, apoditrysia, obtectomera and macrolepidoptera.
these hypotheses are, however, based on a limited number of characters; two for apoditrysia, two for obtectomera, and one for macrolepidoptera.
furthermore, they have not been verified across all ditrysian superfamilies, families or subfamilies.
in recent years, several studies based on molecular methods have brought new information to bear on the evolutionary history among and within the ditrysian superfamilies.
these studies range from wide-taxon sampling multigene studies to studies with fewer taxa but with large-scale gene sampling or genomic data .
the impact of these studies on our knowledge of the relatedness of the lepidopteran superfamilies is beyond question.
they have supported some of the previous hypotheses of relationships of ditrysian lepidoptera, but have challenged or refuted others, e.g.
the monophyly of tineoidea or the position of papilionoidea.
the main division of ditrysia into the nested clades apoditrysia, obtectomera and macrolepidoptera has been backed to some extent, but with changes.
for example, genetic evidence has supported the morphology-based suspicion of rammert  <cit>  and kaila  <cit>  that gelechioidea are a superfamily in obtectomera to which they were not thought to belong .
also, there is mounting evidence that butterflies, skippers and hedylids  are more closely related to ‘microlepidopteran’ groups of obtectomeran moths than to the core macrolepidoptera.
the traditional concept of the superfamily papilionoidea has also changed.
the monophyly of tineoidea has repeatedly been found to be unfounded, and the monophyly of zygaenoidea has also been challenged as cyclotornidae and epipyropidae have not been unambiguously included.
despite these advances, the evolutionary history of the ditrysian lepidoptera is still full of questions.
the low support values obtained in many studies, especially along the apoditrysian “backbone”, have made the suggested relationships uncertain.
a proposed reason for the difficulties in inferring evolutionary relationships among the ditrysian superfamilies is their postulated explosive diversification in the mid- to late cretaceous into the early tertiary, contemporaneous with the radiation of flowering plants  <cit> .
when diversification is rapid, little molecular or morphological change has time to accumulate as evidence of common ancestry.
this weak signal may also be obscured by subsequent changes in the dna and morphology.
with a limited spectrum of genes it may not be possible to extract this weak signal  <cit> .
to overcome this, phylogeneticists, in lepidopterology as elsewhere, have turned to comparative genomics.
by comparing large portions of or entire genomes, clarification of the deeper evolutionary relationships is hoped to be achieved.
recent phylogenomic studies have indeed shown much stronger support values for relationships among the superfamilies than the first studies based on individual gene sequences.
however, low taxon sampling still is an issue in these studies.
with low taxon sampling, distantly related terminal taxa may group together with strong support values in the absence of more closely related taxa.
a wider taxon sampling and more refined models of molecular evolution are needed to overcome this problem  <cit> .
conflicting hypotheses with strong support values  are hard to evaluate if alternative evidence is not available.
such evidence is traditionally sought in the morphology of the organisms.
synapomorphies, i.e.
shared derived characters, are evidence of descent from a most recent common ancestor and inform us of a common evolutionary history.
however, as morphology has never been systematically studied across all ditrysian groups, and the data subjected to a rigorous phylogenetic analysis, we lack morphological information of their affinities.
the alleged morphological homogeneity of ditrysia, the laboriousness of comparative morphology, and also the immense number of species of the order have without doubt held back their classification.
several earlier hypotheses based on morphological features that claim to diagnose certain groups also remain untested across larger groups.
a frequent problem in phylogenetic analyses is taxa that for some reason  do not find a stable phylogenetic position and may even have a nearly random placement.
such taxa are often called rogue taxa  <cit> .
the presence of such taxa may effectively remove resolution of phylogenetic trees .
they have also been shown to affect the topology of the entire tree, thus distorting the phylogenetic relationships among other clades  <cit> .
rogue taxa exert a great influence on the statistical support values of the nodes in phylogenetic trees, obscuring the meaning of these values.
increasing genetic data may not necessarily solve the problem, as fast-evolving taxa may display deviating patterns across large parts of their genomes.
exclusion of rogues is indeed a common practice to try to salvage reliable-looking topologies for the remaining ‘well-behaved’ taxa  <cit> .
moreover, it is likely that a lot of such pruning goes unreported as it is done during analysis, long before manuscript preparation.
exclusion of rogues, however, means discarding data, and is therefore not an optimal way to treat these taxa, which may be crucial for the study in question.
taxa may behave as rogues in molecular or morphological analyses, or both.
however, there is evidence that by combining morphological and molecular data, the negative effect of rogues is reduced and these taxa may find a stable position  <cit> .
thorough examination of the morphology across all ditrysian groups is also needed to obtain information on character evolution.
in a phylogenetic context, morphological structures may show an order of evolution – a polarity.
some characters may also explain the success of certain groups of lepidoptera.
such characters are also called key innovations as they enable, for example, the exploitation of resources, help avoid predators or survive adverse conditions  <cit> .
the time of emergence of such characters can also be investigated to assess whether certain factors in the geologic history of the earth have favored lepidoptera with particular attributes.
fossil lepidoptera are extremely hard to identify and assign to particular taxa, because in adults the scale covering hides diagnostic characters, and fossils of soft-bodied larvae are rare  <cit> .
however, a better understanding of the morphological diversity of lepidoptera and of characteristic features of different groups could prove useful for reliable identification of fossil lepidoptera.
the time of appearance of a certain character is also of interest.
dated and securely identified fossils are used as calibration points in dating phylogenies to obtain estimates of times of divergence  <cit> .
to address these needs, we examine the morphology of ditrysian lepidoptera across  <dig> of the  <dig> currently recognized superfamilies with the aim of providing phylogenetic information on the relationships among them.
we examine  <dig> morphological characters from larvae, pupae and adult males and females of over  <dig> exemplar taxa.
we aimed to obtain samples as comprehensively as possible across lepidopteran subfamilies, and added a number of taxa with unknown or unsupported placement.
special emphasis was placed in the comprehensive inclusion of non-macrolepidoptera groups and character sampling of all life stages.
we have aimed to examine real samples of as many exemplars as possible ourselves.
these data have been supplemented and cross-checked using literature sources.
we analyze the morphological data alone and combined with the molecular data set published by mutanen et al.
<cit> , supplemented with additional taxa and sequences for some taxa.
we report the results of the phylogenetic analyses and compare these to the results of other phylogenetic studies on ditrysia.
we re-examine previous hypotheses based on morphology and present new morphological evidence supporting affinities among clades.
we also discuss the effects of homoplasy and methodological issues on the results.
methods
taxa
the exemplar taxa for the study were chosen to represent as many ditrysian subfamilies as possible.
we estimate that about  <dig> % of the  <dig> subfamilies recognized by kristensen et al.
<cit>  are included in the present study.
an exact number is impossible to give as the concepts of many subfamilies vary or are ambiguously defined.
the full dataset consists of  <dig> exemplar species .
for  <dig> taxa there are only morphological data, for  <dig> taxa only dna data, and for  <dig> terminal taxa there are both morphological and molecular data.
in some cases either the life stages or sequence data were not available from the same species, and so were supplemented by data obtained from a closely related species in the same genus.
these cases are indicated in additional file  <dig>  in addition,  <dig> non-ditrysian taxa were included in the analyses.
all trees were rooted using micropterix calthella.
the authors examined the morphology of eight of the  <dig> non-ditrysian outgroup taxa: eriocrania semipurpurella, hepialus humuli, andesiana lamellata, lampronia capitella, incurvaria pectinea, azaleodes micronipha, ptyssoptera sp.
and tischeria ekebladella.
the remaining  <dig> taxa were represented by only molecular data.
the molecular data are largely the same as in mutanen et al.
<cit> , but excluding polypogon strigilatus, and supplemented by  <dig> additional taxa and  <dig> sequences .
these additional taxa are taxa for which fresh material had become available since the publication by mutanen et al.
or other relevant taxa and data published in previous phylogenetic studies .
examination of morphological characters
morphological data were collected from  <dig> species .
larval characters were available and coded for  <dig> of these species, pupae for  <dig> species, and adult characters for  <dig> species.
when all life stages of same species were not available, character coding was done from a closely related species in the same genus for which material was available.
there are seven such cases in our morphological data, indicated in additional file  <dig>
larvae examined were preserved either as dry inflated or in alcohol.
pupae or pupal exuviae were also examined from either dry or alcohol preserved samples, or whenever possible, both.
adult characters were coded from mounted specimens.
the wings of the adults were removed and the body treated in  <dig> % koh solution to clear it of scales, lipids and muscle tissue.
the exoskeleton of the adult was kept in alcohol.
the specimens were obtained mostly from the collections of the finnish museum of natural history, but also from several other museums and a large network of international collaborators .
larvae and pupae were examined for external characters.
due to the paucity of material of the potentially most informative first instar larvae, the larval characters were coded from final instars.
adults were coded mostly for characters of the exoskeleton.
characters known to have much intraspecific variation and/or considerable problems with assessment of homology among more distantly related groups, such as genital structures and wing venation, were not included in the present matrix.
wing venation characters, except for the distinction between homoneurous and heteroneurous venation, were left out from the analyses.
although homologous characters of wing venation are fairly easy to establish for many of the larger taxa, even they posed considerable problems with homologies and continuous nature of variation.
these problems were even more pronounced with many groups of smaller moths  for which homologies, usually losses of veins and considerable intraspecific variation  <cit> , were impossible to establish.
inclusion of venation characters would have considerably increased the amount of uncertain homologies, thus more likely decreasing rather than increasing the reliability of our morphological data.
the head was not opened.
the basal wing sclerites were not studied because the observation of the sclerites in the small-sized exemplars would have needed very time-consuming and delicate preparation of special slides.
preliminary exploration of these structures indicated that features of these structures would largely prove continuous, thus a thorough study of them was deemed implausible in the time-frame allocated to the study.
as the focus of the study was to seek morphological evidence to infer phylogenetic relationships among the superfamilies, less attention was paid to intra-family level relationships.
some apparent autapomorphies were coded when they were known to be present in other members of the taxon in question.
literature from which characters for some of the taxa were coded is listed below the character list in additional file  <dig>  many of the characters used in this study are not original but were initially taken from literature, and examined across the whole taxon sampling.
character observations were made solely using light microscopy.
a scanning electron microscope would have greatly expanded the number and probably also accuracy of observable characters, but was not available.
the observations were made with leica mz <dig>  mz <dig> , leitz diaplan phase contrast microscope and wild m <dig>  character data management and storage was implemented in morphobank  <cit> .
the morphological character matrix comprises  <dig> binary or multistate characters .
of these  <dig> are larval characters,  <dig> are pupal characters, and  <dig> are adult characters.
both inapplicable and unclear character states were coded with question marks.
continuous characters were not included in the study.
those typically include characters showing high plasticity, characters based on relative variation in size, shape, or degree of sclerotization, all these being factors that complicate assigning characters into discrete character states.
our focus was to find unambiguous characters suitable as potential diagnostic characters in identification keys.
molecular data
genetic data, totaling  <dig> base pairs, were sequenced from  <dig> taxa.
dna samples were obtained mostly from museum specimens and additional specimens collected during field work in  <dig> in new zealand  and tasmania .
genetic data includes cytochrome oxidase subunit i gene  from the mitochondrial genome and seven genes from the nuclear genome; elongation factor-1a , ribosomal protein s <dig> , carbamoyl phosphate synthetase domain protein , cytosolic malate dehydrogenase , isocitrate dehydrogenase , wingless, and glyceraldehyde-3-phosphate dehydrogenase .
sequence lengths, gene summaries and genbank accession numbers are presented in additional files  <dig>   <dig> and  <dig>  respectively.
dna was extracted from legs, sometimes from other body parts using qiagen’s dneasy extraction kit.
dna amplification and sequencing largely followed the protocol explained in wahlberg and wheat  <cit> , but with slight modification and optimization applied to clean-ups of pcr reactions and a different sequencing reaction purification kit  applied later.
sequencing was performed mainly with an abi  <dig> capillary sequencer , but a smaller part was performed with an abi prismr 3130l capillary sequencer .
details on dna extraction are explained in mutanen et al.
<cit> .
the program voseq  <cit>  was used to generate sequence data files and gene summaries .
matrices were prepared for concatenation in mesquite  <cit> .
analyses
the parsimony analysis of morphological data  was done with tnt tree search using new technology search  <cit>  with the following settings; all four search options  selected; find minimum length  <dig> times; initial addseq 15; random seed ten.
morphological characters were treated as unordered and equally weighted.
inapplicable, unclear and missing characters were coded with question marks.
for examination of trees and character distributions the maximum parsimony trees were exported to winclada  <cit> , where they were treated as follows: select all trees; keep best trees only.
data sets of pupal, larval characters and adult characters were examined separately and combined.
morphological data were also analyzed with maximum likelihood methods carried out with raxml  <cit>  in cipres  <cit>  using the mk model  <cit>  for morphological data.
dna data and combined datasets were analyzed using maximum likelihood methods carried out with raxml  <cit>  in cipres  <cit> .
analyses of combined data based on three different data sets were run:data set of  <dig> taxa including only taxa for which both dna and morphological data were available;
data set of  <dig> taxa including all the taxa for which dna data were available, for  <dig> of these taxa also morphological data were available and were included, and;
data set including all  <dig> taxa:  <dig> of these having both molecular and morphological data,  <dig> taxa with only dna data and  <dig> taxa with only morphological data .
analyses were run with all codon positions retained and with third positions excluded from all genes except ef-1a, which evolves more slowly than any of the other genes  <cit> .
this methodology was also followed in mutanen et al.
<cit> .
in another set of analyses, both molecular and morphological data were categorized and partitioned according to the rate of evolution of the characters with the programs tiger  <cit>  and ratepartitions .
molecular data were divided into seven partitions, morphological data into five.
the analyses were run with all partitions present but also with the partitions with the fastest evolving molecular characters removed .
omission of further partitions began to break up well-established clades and produce spurious groupings of taxa.
the fastest evolving morphological character partition  was also omitted in some trial analyses.
however, it became clear that the fastest evolving morphological character partition also contained phylogenetically informative characters and its removal caused false disruptions of even closely related taxa, and was therefore retained in subsequent analyses.
the gtr + g model was used for the dna data and the mk model  <cit>  for morphological data.
the ascertainment bias correction option  was selected for partitions of morphological data, which include only variable characters.
the option to conduct bootstrapping and search for the best-scoring ml-tree in a single program run was selected.
supports for nodes were evaluated with 100– <dig> bootstrap replicates.
the resulting trees were examined in figtree .
mesquite  <cit>  was used to explore character distribution on the trees based on combined data.
we compared our result to other recent studies of ditrysian phylogeny based on molecular data.
comparison of the results with these studies was not straightforward as many of the other studies have much lower taxon sampling and several superfamilies included in the present study are not represented.
comparison of support values obtained for certain clades in different studies is also problematic.
we have deliberately retained unstable taxa in our analyses although their presence likely lowers support values, and, especially in consensus trees obtained by parsimony analyses, their presence significantly collapses resolution.
to overcome this, and more effectively filter information from the mp trees, we explored them individually.
the character and taxon sampling of the present study were not specifically designed for researching the relationships among subfamilies and families, but we do report on these if they significantly deviate from current concepts.
relationships among non-ditrysian lepidoptera were not explored.
many of the morphological characters supporting subfamily, family or superfamily level clades have been found by previous authors and we do not repeat all of them.
more detailed information on characters and corresponding references are found in additional file  <dig>
RESULTS
we report the results of the maximum likelihood  and parsimony  analyses of the morphological dataset alone, and compare these results to those of the ml analysis of combined morphological and molecular data.
we consider the total evidence ml analysis   our main result.
this decision as to the priority of results is based on the principle of inclusion of the maximum amount of data, even though we recognize that this may also bias the results, as may, for example, the inclusion of unstable taxa.
moreover, as källersjö et al.
<cit>  have shown, third codon positions, even if possibly often saturated, also contain phylogenetic signal that may be informative when combined with other data.
we also compare the results with those obtained with third codon positions excluded  and the “noisiest” molecular character partitions i.e.
those with the fastest evolving characters  removed .
the fastest evolving sites were identified with the programs tiger, which identifies rapidly-evolving characters  and ratepartitions, which partitions data according to their rate of evolution .
with the fastest evolving molecular character partitions removed, the results were nearly identical to those obtained with all data retained but produced some spurious groupings, which contradict results of other recent studies .
in the results of analyses based on data with the noisiest partitions removed, bootstrap support values for the clades were in general lower or more or less equal to those obtained for the same clades in the analyses based on data with noisy characters retained.
the results of the analyses based on dna data alone  are not discussed separately as they are very similar to those in mutanen et al.
<cit> , discussed therein.fig.
<dig> phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data;  <dig> taxa;  <dig> characters .
***both dna and morphological data; **only morphological data; *only dna data
fig.
<dig> phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data;  <dig> taxa;  <dig> characters .
***both dna and morphological data; **only morphological data; *only dna data
fig.
<dig> phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data;  <dig> taxa;  <dig> characters .
***both dna and morphological data; **only morphological data; *only dna data
analyses based on morphological data only
the parsimony analysis yielded  <dig> equally most parsimonious trees .
the resolution of the resulting trees was generally uniform, but interrelationships of certain superfamilies varied among the trees.
the strict consensus tree is very poorly resolved due to the adverse effect of some taxa assuming two alternative positions , as well as the paraphyly of gelechioidea observed in some of the mp trees.
an example of one of the most parsimonious trees is shown in additional file  <dig>  most superfamilies are recovered as monophyletic.
however, some results appear spurious, such as the consistent grouping of the non-ditrysian hepialus in tineoidea.
bombycoidea, and in some mp trees also drepanoidea and geometroidea, were nested within noctuoidea.
the unstable position of ochsenheimeria is likely due to the internal feeding mode of its larva having caused morphological modifications that remove it from the otherwise constantly recovered association with ypsolopha , the larvae of which feed externally on leaves of various trees and shrubs.
other differences in the topology were minor rearrangements of apical genera.
the tree resulting from the maximum likelihood analysis of the morphological data set is shown in additional file  <dig>  there is no support for the deeper nodes, but most superfamilies are recovered as monophyletic except for those in macroheterocera, which are intermixed.
separate analyses on data sets of larval, pupal and adult morphology were also run, but are not discussed, because the resulting trees were nearly unresolved due to the small number of characters in comparison to the number of taxa included in these analyses.
combined analyses
phylogenetic relationships
the results are in many respects similar to those obtained in the study by mutanen et al.
<cit>  as the molecular data are largely the same, especially when third codon positions were excluded.
considering that the morphological data in the present study consisted of only  <dig> % of the total amount of data , the phylogenetic signal from these characters appears disproportionately strong, and their inclusion significantly affects the topology of several clades, while supporting the molecular results elsewhere.
as will be discussed in more detail below, morphological attributes stabilized some taxa that had been unstable in dna-only analyses and thus propose a hypothesis about their phylogenetic affinities.
likewise, molecular data anchored taxa that were unstable in our analyses of morphological data alone.
clades supported by both morphological and dna-data in general obtained higher support values in the combined analyses than when either type of data was analyzed alone.
however, we also obtained some intriguing patterns in the topology that strongly conflicted with the results of several recent studies based solely on genetic data.
these patterns show the strength of morphological data over the eight gene regions and the usefulness of different types of data in evaluating different phylogenetic hypotheses.
the overall topology of the phylogenetic trees based on the three different data sets  is very similar and the weakly supported backbone nodes are weak in all trees .
although the support values are the highest for the smallest data set including  <dig> taxa , and weakest in the trees including all  <dig> taxa , we choose to present the latter .
the topology of the trees is essentially the same and the tree based on the larger taxon sampling also includes interesting information on the placement and performance of taxa with only either morphological or molecular data.
we report support values obtained from all three analyses .
phylogeny of main ditrysian clades
ditrysia
the monophyly of ditrysia is recovered in the maximum likelihood analyses based on dna, morphological and combined data sets .
the support values for the ditrysian clade in the ml analyses of combined data are as follows:  <dig> taxa:  <dig> %;  <dig> taxa:  <dig> %;  <dig> taxa:  <dig> %.
in the parsimony analysis of morphology alone, hepialus is placed among ditrysia .
probable reasons for this are that in immature stages there appears to be no morphological characterization that would distinguish ditrysia from monotrysians, and the methodological decision to exclude morphological characters only applicable for resolving interrelationships within monotrysia.
the adult characters supporting the ditrysian clade are the name-giving ditrysian female reproductive system i.e.
one opening for mating and a separate one for laying eggs  <cit> , and the absence of structures found in some but not all the outgroup taxa such as the division of sternum  <dig>  an elongated metathoracic trochantin, and the tergosternal connection of tergal origin.
division of ditrysia
our results support previous assertions on the scarcity of clear morphological characters that would define larger assemblages of ditrysian superfamilies  <cit> , and no new such characters were found.
analysis of the morphological data alone or combined with the molecular data set did not improve the resolution of the deeper nodes obtained by mutanen et al.
<cit> .
support values for several of the deeper nodes , but also of many other nodes, are low.
this is in part due to our deliberate decision to not remove unstable taxa from the analyses.
the unstable taxa varied depending on the data set and method used.
in the maximum likelihood analyses of combined data, the following taxa were unstable: imma, millieria, douglasiidae, tanaoctena, cadmogenes and several small superfamilies.
there were also other taxa with unstable positions, but only within superfamilies, e.g.
lyonetia in yponomeutoidea.
characters previously proposed by minet  <cit>  to circumscribe subclades of ditrysia were tested across the 300+ taxa representing  <dig> out of  <dig> superfamilies.
the results show that although several obviously contain phylogenetic signal, there is ambiguity in the interpretation and distribution of details of these characters, as already implied by minet  <cit> .
apoditrysia
characters of the adult sternum  <dig> have been proposed to define the apoditrysia, i.e.
all ditrysia except tineoidea, gracillarioidea, yponomeutoidea and gelechioidea  <cit> .
recent studies have, however, repeatedly suggested that gelechioidea belong to obtectomera, a clade nested in apoditrysia.
apodemes of sternite  <dig> of apoditrysia were defined by minet  <cit>  to be short, posteriorly enlarged and with anterior corners of sternite  <dig> distinctly produced laterally.
the division between non-apoditrysian and apoditrysian sternum types is also described as the tortricid type versus the tineid type of thoraco-abdominal articulation  <cit> .
the apoditrysian type of sternum  <dig> is indeed never present in the tineoidea, gracillarioidea and yponomeutoidea, and present in most apoditrysians.
however, a condition more similar to the non-apoditrysians is found in some cossoidea  <cit> , zygaenoidea  <cit> , millieriidae  <cit> , gelechioidea  <cit> , and in the present study, it is observed that clear lateral extensions are occasionally absent even in some macroheterocera, a clade nested in obtectomera.
in some taxa, sternite  <dig> is modified due to the presence of hearing organs  making the presence or absence of the lateral extensions difficult to assess.
in mimallonoidea, the presence of the character was also ambiguous.
sternum  <dig> has an additional modification that, according to minet  <cit> , is an autapomorphy for the earlier concept of drepanoidea : the lateral extension of sternum  <dig> is long and reaches in front of the spiracle all the way to the anterior corner of sternite  <dig>  in addition to drepanidae and epicopeiidae, this modification is present in pterophoroidea, schreckensteinia, coleophora in gelechioidea, carposinoidea, crambidae, calliduloidea and papilionoidea  .
obtectomera
a pupal structure – body appendages  firmly glued to the body  – has been considered a character of considerable phylogenetic significance.
a further specialization, i.e.
pupa with immobile intersegments between abdominal segments  <dig> and  <dig> , was used by minet  <cit> , building upon the work of chapman  <cit> , to define the clade obtectomera comprising pyraloids, butterflies and other macrolepidoptera.
in some groups the abdominal segments 1– <dig> are variably fused, but the mobility of 3– <dig>  is clear-cut.
therefore the latter was taken into account in character coding.
presence of this pupal type in yponomeutoidea, gelechioidea, epermenioidea and alucitoidea was interpreted to represent convergence  <cit> .
now it appears that carposinoidea , epermenioidea, alucitoidea and gelechioidea belong to obtectomera sensu minet  <cit> , so convergence needs only to be assumed for some yponomeutoidea and possibly for the unstable immoidea.
in the present study, this interpretation is supported when third codon positions or partitions  <dig> and  <dig> are removed.
with third codon positions included, the cossoidea + sesioidea + zygaenoidea assemblage is nested in obtectomera, which, on a morphological basis, seems less plausible than the conventional view of them being true non-obtectomerans.
the case of yponomeutoidea is not quite straightforward.
the majority of yponomeutoids studied was observed to have the pupa with intersegment 3– <dig> not entirely either mobile or immobile, but rather having a limited ability for movement.
this condition was considered its own character state, supported by experiments where pupal exuviae were softened so that the potential for mobility could be evaluated.
several yponomeutoids with their pupal intersegment 3– <dig> found entirely immobile are taxa whose thoracic structures largely cover the ventral surface of the abdomen, which alone may effectively prohibit the ability to move abdominal segments and lead to the fusion of the thoracic segments.
immoidea, with an obtect pupa, was associated with choreutidae, i.e.
non-obtectomerans, although with negligible support.
the pupal skin in immoidea is thin, and the level of intersegmental mobility is ambiguous.
minet  <cit>  considered their pupa obtectomeran while common  <cit>  interpreted it to be of the non-obtectomeran type.
another character later proposed by minet  <cit>  to redefine the obtectomera is the presence of a dorsal lobe or protrusion on the adult pretarsus.
as can be noted also from the pictures in minet  <cit> , the presence of the setose lobe on the pulvillus is not always obvious.
the lobe can be small and not easily distinguishable from the non-obtectomeran type of pulvillus, especially if the latter has long setae on the dorsal side.
this character was also examined across all taxa in the present study, and the different character states applied here were: pulvillus simple, pulvillus with setose outgrowth, or pulvillus bifid.
setose outgrowths were clearly visible at the base of the dorsal side of the pulvillus in most taxa assigned to obtectomera.
in gelechioidea, however, the pulvilli were simple.
in a few bombycoids and in pieridae, the pulvillus was also more of the simple type without projecting setae.
in other papilionoids, the pulvillus was distinctly bifid and differed from those with a protrusion with setae.
imma and copromorpha also had distinct setose lobes on their pulvilli.
the resolution limits imposed by the light microscopes used in the present study lowered confidence in coding this character for some of the smaller taxa.
in conclusion, our result suggests that even though the obtectomera is probably largely a valid clade, there are several superfamilies whose assignment to this clade is uncertain, and the degree of homoplasy appears higher than generally presumed.
fänger  <cit>  also examined this character across over  <dig> ditrysian families with a scanning electron microscope and confirmed its validity with some exceptions.
macrolepidoptera and macroheterocera
the apomorphy proposed to delimit macrolepidoptera  <cit> , the first axillary sclerite of the forewing with an elongate angle, was not included in the present matrix due to the continuous nature of the character, which compromises its usefulness.
the traditional concept of macrolepidoptera has recently changed as molecular analyses have repeatedly suggested that pyraloidea appears to be more close to the core macrolepidoptera than papilionoidea and calliduloidea.
the clade including the rest of the traditional macrolepidoptera  is referred to as macroheterocera  <cit> .
in several recent molecular analyses, mimallonoidea has been placed as the sister-group of macroheterocera.
this is also the result obtained in our analysis with the three codon positions included .
in analyses without the third codon positions or fast evolving partitions , the pattern changes and the arrangement of clades varies, with all alternatives having low support.
although the macroheteroceran clade is recovered in the combined analysis, no unique morphological character was found to unite them.
interrelationships and composition of ditrysian superfamilies
the most basal clade of ditrysia recovered in the combined analysis  is a monophyletic tineoidea with the following support values: 268: 88; 422:  <dig> and 473:  <dig> %.
they are also monophyletic in the result based on dna with all three codon positions retained and those of combined data with the fastest evolving partitions removed .
tineoidea are divided into two clades.
one clade includes the monophyletic families eriocottidae and psychidae , the other tineidae .
in recent phylogenetic studies based on dna data alone, a monophyletic tineoidea has only been obtained in the studies by regier et al.
<cit>  and  <cit>  in analyses of unaltered nucleotides from all three nucleotide positions .
in the present study, when the third codon positions are removed , tineoidea are paraphyletic, a result similar to that obtained in other recent studies based on sequence data  <cit> .
tineoidea are then paraphyletic, with eriocottidae, psychidae, and tineidae as separate clades.
the present study did not include exemplar species of meessiidae , included in the phylogenetic study on tineoidea by regier et al.
<cit> .
in the results of their study, representatives of meessiidae were excluded from all previously recognized tineoid families and placed as sister to all other tineoidea + ditrysia, thus also contributing to the paraphyly of the superfamily.
had meessiidae been included in the present study, the monophyly of the superfamily could have been challenged.
tineoidea were also recovered as monophyletic in the maximum likelihood analysis of morphological data , but not in the parsimony analysis .
adults in all three tineoid families included in the present taxon sampling are characterized by the presence of a distinctly shortened proboscis with loose galeae  and pseudapophyses found in most females  <cit> .
shortening or absence of galeae is, however, found throughout the ditrysian superfamilies.
the looseness of galeae could not be coded because most of the short galeae in other groups also loosened when treated in koh solution, and on non-denuded specimens the scale coverage effectively prohibited the reliable observation of this character.
light microscopes were insufficient for comparison of differences in the microstructures holding galeae together.
no support for the monophyly of tineoidea was found from larval or pupal characters.
within tineoidea, monophyly of tineidae is supported by a unique synapomorphy, a nearly triangular cap dorsally covering the base of the larval antenna .
dryadaula, recently proposed to be a new family, dryadaulidae regier et al.
<cit> , is nested in tineidae, possibly because only morphological data from the adult male and female were available and the larval synapomorphy was not studied.
psychidae larvae share the median fusion in all thoracic legs; the members of the subfamily oiketicinae are characterized by the horizontal position of the larval spiracles of t <dig> fig.
<dig> unique tineid apomorphy: more or less triangular cap dorsally covering base of larval antenna.
morophaga choragella
metapherna robinson & nielsen  <cit>  is an australian genus that has been placed in myrmecozelinae in tineidae, although robinson & nielsen  <cit>  noted that it has, e.g., an unusually long proboscis for a tineoid.
in the analyses of dna data-only and of combined dna and morphological data , metapherna salsa  is placed between tineoidea and the rest of ditrysia .
in the morphology-only ml analysis metapherna has an unresolved position within non-apoditrysians .
in the mp analysis it is placed either as a basal ditrysian or in yponomeutoidea, but never in the tineoidea .
based on genetic evidence and their divergent morphology, species of the metapherna group  obviously may eventually require their own family.
in the analyses of combined data , metapherna is basal to a large clade  divided into a rather weakly supported clade including taxa assigned to yponomeutoidea and gracillarioidea , and the apoditrysian clade .
yponomeutoidea plus gracillarioidea, with douglasiidae excluded, have also been found to be monophyletic in several other recent phylogenetic studies on the phylogeny of ditrysia , and also in two studies focusing on the phylogenetic affinities within these superfamilies  <cit> .
in analyses based on morphology only, douglasiidae was included in gracillarioidea, as has also been formerly suggested  <cit>  .
however, in our preferred result  it is associated, with uncertain position, among apoditrysia, as in those recent molecular studies in which it has been included.
in the present study, we did not find any unique characters common to all members of gracillarioidea and yponomeutoidea.
this grouping was supported by a combination of homoplasious characters such as characters of antennal scaling and the tendency to have the proximal pair of tibial spurs of the hindleg on the proximal half of the leg.
however, the transverse costa at the base of sternum a <dig> found in most yponomeutoidea and proposed to be a possible groundplan autapomorphy for the superfamily  <cit>  was also found to be present in roeslerstammiidae : roeslerstammia , vanicela, amphithera and the female of thereutis.
a similar costa, although longer and often reaching the sternal apodemes, was present on sternum  <dig> of choreutoidea .
some tortricids also exhibited a sclerotized bar or patch in a similar position.fig.
<dig> transverse costa behind anterior margin of sternum ii.
a
yponomeuta evonymellus ; b
roeslerstammia erxlebella
yponomeutoidea were found to be monophyletic  and the presence of the pleural lobes on sternum  <dig> of the males was verified.
pleural lobes were found to be present in nearly all yponomeutoids.
in some cases, the expansion of the lobe was not strong and it was difficult to determine if lobes were present.
they were absent in cadmogenes literata , which is currently placed in plutellidae, but may eventually need to be placed in its own family  <cit> .
in our study, only morphological data were available for cadmogenes.
its position was unstable, and it grouped most often with some other unstable taxa, such as millieria or tinagma.
as in previous studies, the position of lyonetia was unstable within yponomeutoidea.
intriguingly, it was rarely associated with leucoptera , which has been considered a subfamily of lyonetiidae.
this result is in line with that obtained by sohn et al.
<cit> .
in our result heliodinidae  was grouped together with ypsolopha and ochsenheimeria, which deviates from all former suggestions, where it most often is considered close to bedelliidae  <cit> .
gracillarioidea  had no distinct immature or adult characters common to all taxa in the character sampling of the present study.
several small superfamilies, families and taxa, the position of which has greatly varied in analyses based on genetic data — millieriidae, tinagma, tanaoctena, epermenioidea, immoidea, carposinoidea, pterophoroidea, and choreutoidea — do not find a stable position in the combined analysis either.
they tend to form one or two agglomerations, or join other groups with extremely low support values.
a repeatedly recovered pattern is that alucitoidea, epermenioidea and carposinoidea form a monophylum.
however, interesting observations were made regarding these taxa.
characters listed below could prove useful in diagnosing families even though they may not have phylogenetic signal at higher levels.
for example, the sculpture of the larval mentum is dentate in two groups, viz.
choreutidae and alucitoidea s. l., respectively.
the larvae of carposinidae have a characteristic bilobed extension between the pre- and submentum.
in the analyses based on molecular data  only, schreckensteinia is an unstable taxon, taking multiple alternative positions in the resulting phylograms.
in the combined analysis, schreckensteinia consistently groups with the urodidae .
they share several, yet to some extent homoplastic, synapomorphies in their larvae, notably the closely set stemmata, granulose sculpture of the mentum, medially narrowed prolegs with mesoseries on a3- <dig>  and the closely set prolegs on a <dig>  the adult urodus  has pigmented patches on the cuticle of its mesothorax, which are also found in roeslerstammiidae .
immidae have some obtectomeran features such as the setose outgrowth on the dorsal side of the pulvilli  <cit> , and a possibly obtectomeran pupa .
immidae have previously been placed in sesioidea  <cit>  and as an unplaced obtectomeran superfamily of their own  <cit> .
alucitoidea are monophyletic and are joined by isonomeutis amauropa meyrick,  <dig> in the combined analysis .
isonomeutis has previously been placed in copromorphidae.
in our analysis, both with and without the third codon positions or fast evolving partitions removed, alucitoidea, carposinoidea and epermenioidea are grouped together .
the family alucitidae is consistently nested in tineodidae.
a tongue-shaped lobe laterad of the larval submentum is an autapomorphy for isonomeutis + alucitidae + tineodidae.
on the basis of this evidence , we here synonymize tineodidae with alucitidae, syn.
nov.
the adult character common to pterophoroidea, here only including pterophoridae, is the acute angle of the metepimeron  <cit> .
however, the metepimeron of tanycnema anomala  is also rather acute.
pterophoroidea also have the prespiracular extension of sternum  <dig> discussed above in the section on apoditrysia.
galacticoidea are consistently recovered as the sister group to tortricoidea , which is also found in the results of the ml analyses of morphology .
galacticidae also group with tortricidae in regier et al.
<cit> .
tortricoidea are monophyletic in all analyses.
tortricoidea have been characterized by the flat papillae anales of the females  <cit> .
however, similarly flat papillae were occasionally also present in several other groups, e.g.
zygaenoidea, pyraloidea, prodidactis, calliduloidea and noctuoidea.
these taxa, however, have attributes that suggest the character has evolved repeatedly.
the heliocosma group was removed from tortricidae by horak and common  <cit> , and left without family assignment.
in the study by mutanen et al.
<cit>  it was, together with another unassigned taxon, piestoceros, grouped with brachodidae.
regier et al.
<cit>  recovered heliocosma together with tortricoidea.
in our result of the analysis with three codon positions included, heliocosma, piestoceros and choristis form a monophylum that, with very weak support, forms the sister-group of the cossoidea + sesioidea + zygaenoidea assemblage .
with the third codon positions removed, they are associated with immoidea .
when the fast evolving partitions are removed, they group with millieria .
this group is thus one of the most unstable ones in lepidoptera.
no hypothesis of its position can be preferred over others at present and morphology does not offer clear clues about its affinities either.
zygaenoidea , sesioidea   and cossoidea  are all, respectively, found to be monophyletic in our results and are grouped together in the same clade  .
in the smallest data set  synechodes coniophora , however, is placed within cossidae.
no other brachodids were included in the analyses of the  <dig> and  <dig> taxa data sets due to the unavailability of dna data.
in the results based on dna data only zygaenoidea are monophyletic , but sesiidae and cossidae are intermixed .
in other recent studies based on genetic data, sesiidae and cossidae have not been found to be monophyletic either .
sesiidae and cossidae are also considered to belong to the same superfamily according to the most recent classification of lepidoptera  <cit> .
epipyropidae and/or cyclotornidae were not included in zygaenoidea in mutanen et al.
<cit> , cho et al.
<cit> , regier et al.
<cit>  or bazinet et al.
<cit> , but in the present study belong there with moderate support .
the support seems largely due to synapomorphies of immatures shared with at least some other zygaenoidea: prolegs of larval a <dig> are approximate and arise on a common protuberance ; the spiracle of abdominal segment  <dig> of pupae is exposed , and there is a sculpted, laterally extended flange on the pupal eyepiece .
species in both families are fully or partially parasitic at the larval stage.
there are indications that molecular evolution of parasitic species may often be accelerated  <cit> , possibly explaining the long branches and unstable phylogenetic affinities of epipyropidae and cyclotornidae in molecular data sets.
this is also the first evidence that epipyropids and cyclotornids might be sister groups .fig.
<dig>
cyclotorna larval segment a <dig>  ventral view, showing swelling on which both prolegs are positioned
minet  <cit>  proposed the presence of two asymmetrical pits on the lower frontoclypeus of adults as a possible autapomorphy for cossidae.
such pits were present in several, but not all cossids.
at times, it was difficult to tell the difference between such pits and large scars of scales i.e.
the round sockets in which scales were attached, on the frontoclypeus.
better magnification would be necessary to observe the details.
in the analyses based on molecular data only, sesioidea are not monophyletic, a result also found in mutanen et al.
<cit> ; regier et al.
<cit> , cho et al.
<cit>  and bazinet et al.
<cit> .
in the combined analysis, however, they are  though with low bootstrap support .
pennisetia , which in our dna-only analyses was not with sesiidae, joins other sesiids because it shares several adult characteristics with them, such as the sesiid wing-locking mechanism and a y-shaped suture/depression on the mesepimeron  <cit> .
adult characters common to brachodidae, sesiidae and castniidae are the elongated posterior tendons of the metafurca  and the extension on the metafurca  <cit> .
several females in these families also have a telescopic ovipositor.
the immature stages of sesioids and cossoids are generally uniform, suggesting no obvious support for the monophyly of sesioidea.
the anteriorly more heavily pigmented ocular diaphragm was suggested as an autapomorphy of sesioidea  <cit> .
we were not able to examine the pigmentation of the ocular diaphragm, because of the varying success in treating the samples in koh solution so that the diaphragm would be sufficiently clearly visible.
the size of the patagia, which has been mentioned as a synapomorphy of sesiidae, castniidae and brachodidae, is a continuous character with much variation among lepidoptera and therefore was not included in the present study.
in addition, paysandisia  shares the character considered autapomorphic for sesiinae, the bag-like protuberance concealing the pleural suture  <cit> .
in our result, cossoidea + zygaenoidea + sesioidea are in obtectomera.
the support for this finding is, however, negligible and disappears when third codon positions are removed .
several recent studies present evidence for the exclusion of these from obtectomera.
our result is nevertheless intriguing as a similar pattern with the cossoidea + zygaenoidea + sesioidea assemblage nested within obtectomera has also been recovered by regier et al.
<cit>  in some of the analyses based on five protein-coding genes.
if this pattern indeed proves the most likely, it emphasizes the level of convergence in morphology .
in the morphology-only analyses, the superfamilies are not in obtectomera.
choreutidae is a family without known phylogenetic position.
choreutidae had previously been placed in a number of different families until minet  <cit>  moved them to a superfamily of their own.
in our results, choreutoidea are monophyletic, and are placed in a very weakly supported clade together with immoidea.
in the result obtained without third codon positions , the placement of choreutoidea shifts to a monophylum together with galacticidae and tortricoidea, a finding also recovered by the mp and ml of morphological data .
the morphological data collected from the choreutids in the present study contain, however, intriguing features.
the cervical sclerites of the adult prothorax are very similar to those of brachodidae as being strongly turned inwards in both groups .
in brachodids, the lower arm of the sclerites is also shortened.
similar, but not as strongly turned cervical sclerites were present in castniidae and some cossidae.
other lepidoptera do not have such wide and inwardly turned cervical sclerites .
brock  <cit>  and heppner  <cit>  have proposed a close affiliation of these families but based on characters not included in our study .
the dentate sculpture of the larval mentum characterizes choreutidae, but also the isonomeutis-tineodidae-alucitidae clade, which does not seem closely related to choreutidae.fig.
<dig> cervical sclerites.
a
anthophila fabriciana ;  brachodes appendiculata ;  cydia nigricana
millieriidae are a small microlepidopteran family, which was formerly considered a subfamily in choreutidae  <cit> .
heppner considered millieriidae to be “somewhat of a missing link between other choreutidae and the remainder of the sesioidea”.
based on evidence from morphological and molecular studies  millieriinae were excluded from choreutidae and elevated to family level by rota  <cit> .
rota and kristensen  <cit>  compared the thoraco-abdominal articulation of millieriidae to those of choreutidae and found that unlike choreutids, millieriidae do not have lateral extensions on the apodemes on sternum  <dig>  in the study by mutanen et al.
<cit>  millieria was placed together with the unstable douglasiidae.
in the present study, millieria is grouped, with weak support, in a clade that contains cadmogenes, schreckensteinioidea and urodoidea .
in morphology-only analyses, millieria is associated either with choreutidae or roeslerstammiidae .
the former position may be explained by the external resemblance of the adult of millieria to choreutidae in having a scaled base of the proboscis and scales between the eye and antenna.
however, both features are also found in other groups.
the immature stages of millieria give little information.
the larva is an internal feeder, and devoid of obvious traits that could help group it with other taxa.
its non-obtectomeran pupa is also generalized in appearance, thus not providing data for phylogenetic insights, either.
the position of millieriidae with roeslerstammiidae is based on very homoplasious characters, such as the external invisibility of sutures between pupal abdominal segments 8– <dig>  thus without convincing morphological support.
the affinities of papilionoidea have varied among recent molecule-based studies.
common to all of them is that the superfamily appears to be closer to the ‘microlepidopteran’ grade instead of being among the core macrolepidoptera as has been the conventional view.
another pattern repeatedly recovered in studies based on sequence data  and a study combining molecular and morphological data  <cit> , has been that papilionidae is the sister of the other papilionoid families, including hedylidae and hesperiidae.
in the present study, a surprising result is the emergence of a topology with the conventional concept of papilionoidea, i.e.
hedylidae as the sister group of the rest of papilionoidea .
the position of hesperiidae depends on the data set analyzed and on the inclusion or exclusion of third codon positions.
if these are excluded, hedylidae are still the most basally arising group of papilionoids, but papilionidae and hesperiidae shift places .
this result is in conflict with the recent studies based on genetic data in which both hedylidae and hesperiidae are nested within papilionoidea.
the closest relatives of papilionoidea are not recovered with unequivocal support in recent molecular studies.
the present result suggests that callidulidae is a sister-group to papilionoidea , and if third codon positions are removed, callidulidae are sister to papilionoidea + macroheterocera.
a monophyletic gelechioidea is recovered in the combined  and morphology-only ml analyses.
the adult character supporting the grouping is the scaled haustellum.
the character is not unique to gelechioidea, being also found in tischerioidea, pyraloidea, choreutoidea and millieriidae.
a monophyletic gelechioidea has been found in all recent molecular studies on ditrysian phylogeny and targeted studies on gelechioidea .
although gelechioidea are found to be monophyletic, support values for the clade are never very high when more than just a few exemplar species are included in the analyses.
the large phylogenetic studies based on genetic data have also shown that gelechioidea belong to the obtectomeran clade, contradicting the earlier assumptions of their affinities with the non-apoditrysian yponomeutoidea  <cit> .
the closest relatives to this massive radiation of small moths are not known.
the results of the present study do not give unambiguous clues to resolve this.
several recent studies have recovered a sister relationship between gelechioidea and thyridoidea  <cit>  and in kawahara and breinholt  <cit>  they are sister to the clade callidulidae + thyrididae.
our study suggests, though with negligible support, a sister group relationship between gelechioidea and a clade comprising alucitoidea, epermenioidea and carposinoidea , or if third codon positions are excluded, a sister group position to all other obtectomera .
in our results, hyblaeoidea  are the sister-group of pyraloidea  .
in the results based on the 268- and 422- taxon data sets, copromorpha lichenitis joins hyblaeoidea  whereas in the 478-taxon set it is associated with  <dig> % support with phycomorpha metachrysa, another copromorphoid, and together they weakly form a clade with carposinidae, currently also placed in copromorphoidea.
the grouping of hyblaea and prodidactis is supported by a membranous projection present on the coxae of the male and a modification of the apex of the larval spinneret, as has been explained in kaila et al.
<cit> .
such a modified larval spinneret is also present in thyridoidea and several papilionoids, including all hesperiidae examined.
coxal processes were also found in several species in geometridae, in cyclidiinae , and in anchinia daphnella .
the affinity of prodidactis to several other families has previously been investigated by epstein and brown  <cit>  with inconclusive results, and the genus was assigned a family of its own.
based on the two above-mentioned characters in common with hyblaea and congruent evidence from five genes, kaila et al.
<cit>  placed prodidactidae in hyblaeoidea.
in the present study, the results of the analysis based on combined data and those based on dna data only, the clade formed by hyblaea and prodidactis tended to group basal to pyraloidea .
morphological characters supporting this position are not clear.
the papillae anales of the female of prodidactis are flat, as they are in several pyraloid species.
in the maximum likelihood analysis of morphological data, hyblaea and prodidactis also group together, their position being close to mimallonoidea and thyridoidea.
apart from the present study and mutanen et al.
<cit> , hyblaeidae and prodidactidae have not been present in phylogenetic studies based on genetic data.
thyridoidea are recovered as monophyletic with strong support .
thyrididae are a somewhat mobile group in our analyses, but often group as sister to callidulidae and papilionoidea.
we found the adult character proposed by minet  <cit>  to support the monophyly of the superfamily, the lateral lobe on a <dig>  very difficult to observe.
out of our five exemplars of thyridoidea, the character was more or less clearly present only in thyris fenestrella.
due to the difficulties in identifying this character, it was not included in our matrix.
however, a number of thyridoid species showed on abdominal segment a <dig> a structure similar to the postspiracular bar found in papilionoidea.
the presence of this bar could associate thyridoids in the same clade with papilionoidea, a result also observed in the molecular analysis by mutanen et al.
<cit> .
in most samples of papilionoidea, the postspiracular bars are strongly sclerotized and extend diagonally from the posterior edges of tergum  <dig> towards the anterior edge of sternum  <dig>  the bar may be long and reach behind the spiracle, but sometimes short and harder to define.
in addition, the bar is not always easy to distinguish from other sclerotized, but more perpendicularly directed, extensions of the posterior edge of tergum  <dig> or the anterior edge of tergum  <dig> reaching toward the sternum.
several of the female thyridoidea have fine hair-like bristles by the hindwing acanthi in addition to the usual stouter bristles.
in other recent studies the closest relatives of thyridoidea have varied, but they seem most often to group either close to papilionoidea, gelechioidea, calliduloidea, pterophoroidea or other small apoditrysian families  .
pyraloidea are monophyletic with strong support .
the most obvious adult characters supporting the clade are the scaled haustellum and the characteristically-shaped abdominal hearing organs.
differences in the hearing organs separate pyralidae and crambidae  <cit> .
the sclerotized structures of the tympanal organs have been regarded as homologous in pyralidae, crambidae, dudgeonidae, female uraniidae and geometridae because they are modified abdominal apodemes .
a monophyletic pyraloidea, with monophyletic pyralidae and crambidae, have also been found in all studies based on genetic data in which representatives of both families have been present .
mimallonoidea  are placed as the sister-group of the macroheterocera  .
this position has also been obtained in several other studies based on genetic data .
in regier et al.
<cit> , mimallonidae are in a clade with axia  and doa  and together form the sister-group of drepanidae.
no particular character in the adult morphology explains this position as the sister-group of macroheterocera.
in the morphology-only analyses mimallonoidea are monophyletic but are associated with thyridoidea and hyblaeoidea .
in the results based on the 268- and 422-taxon sets, mimallo amilia, as sole representative of the superfamily in these analyses, is placed as sister to bombycoidea or lasiocampoidea .
the macroheteroceran clade ; 422:  <dig> % ; 473:  <dig> % ) is one of the most strongly supported clades, but the support values for the relationships among the superfamilies are not strong.
the obtained topology  is as follows:  +  + noctuoidea)).
the macroheteroceran superfamilies are recovered monophyletic in analyses based on or including genetic data  lasiocampoidea are nested in bombycoidea), and the monophyly of most is supported by morphology .
however, in morphology-only analyses the phylogenetic signal of these morphological traits does not suffice to recover the clades and they are intermixed .
the proposed relationships among the superfamilies cannot either be unequivocally explained on morphological grounds.
in recent molecular studies, the proposed relationships among the macroheteroceran superfamilies vary .
in the result of the analysis of combined data geometroidea s. l. is obtained .
sematuridae and epicopeiidae are sister groups and together form the sister group of uraniidae + geometridae .
the obtained monophyly and within superfamily relationships are identical to those recovered in a very recent study based on molecular data on the phylogeny of geometroidea  <cit> , which, however, are based on subsample of the present molecular dataset.
rajaei et al.
<cit>  also characterized a new geometroid family, pseudobistonidae, based on pseudobiston pinratanai inoue,  <dig>  not included in the present taxon sampling.
in their study, pseudobistonidae is placed as sister-group of epicopeiidae.
rajaei et al.
<cit>  proposed five apomorphies supporting the monophyly of geometroidea.
two of these are wing characters not included in the present study.
the authors acknowledge that the other three characters are not clear-cut features.
only one of these characters, the reduction of the fenestrae laterales, was included in our study.
although the character is not always easy to define and may also be found in other ditrysia, our observations are generally in agreement with those by rajaei et al.
a list of characters regarded as synapomorphies supporting the relationships among families within geometroidea is provided by rajaei et al.
<cit> .
several of these characters are continuous in nature and are not included in the present study.
in the present study, the monophyly of geometroidea is, however, occasionally broken in the analyses of the 268- and 422-taxon sets and when the third codon positions or the two fastest evolving partitions are removed as sematuridae + epicopeiidae are placed in a position distant from the rest of geometroidea .
epicopeiidae, traditionally placed in drepanoidea  <cit> , also group with sematuridae in regier et al.
<cit> , and are sister to the uraniidae + geometridae clade.
in cho et al.
<cit>  epicopeiidae and sematuridae are also sister-groups in the nt <dig> partitioned analyses, but not in the “all non-synonymous” analyses.
bazinet et al.
<cit>  find epicopeiidae as sister to uraniidae + geometridae .
in mutanen et al.
<cit> , epicopeiidae group with lasiocampidae, and sematuridae is found as the sister-group of uraniidae.
together sematuridae and uraniidae form the sister-group of geometridae.
minet  <cit>  suggested the pointed and elongate ventral arm of the tegula was an apomorphy for geometroidea  and considered the presence of the sharp tegula in epicopeiidae and mimallonidae to be the result of parallel evolution.
the ventral process of the tegula is indeed sharp in sematuridae, epicopeiidae, geometridae, but also in bombycoidea, lasiocampoidea, mimallonoidea and noctuoidea, thus in most macroheterocera.
in drepanidae , the ventral arm of the tegula is blunt, .
the tegula of uraniidae is elongate but spatulate.
minet and scoble  <cit>  suggested the most reliable autapomorphy to unite geometroidea and drepanoidea sensu minet  <cit>   to be the strong anterior extension of pupal forelegs, and the presence of a transverse dorsal groove on pupal segment a <dig>  yet with several secondary losses for both.
the latter character was indeed found only among some representatives of these groups, while the former was highly homoplastic.
apart from these families, it was observed also in bucculatricidae , attevidae and praydidae , urodoidea, schreckensteinioidea, pterophoroidea, thyridoidea, lasiocampoidea, bombycoidea, cimelioidea, most papilionoidea and many noctuoidea.
a morphological character in favor of the relatedness of sematurids with uraniids could be the knob behind the antennae, situated medially of the chaetosemata, present in both sematura and urania .
similar structures, although with more variation in shape, from flat pores to larger triangular protrusions, were found widely in taxa belonging to scopulini/sterrhinae and singular representatives of larentiinae and archiearinae in a study by sihvonen and kaila  <cit> .
such protrusions could thus be found more widely in taxa belonging to geometridae/geometroidea.fig.
<dig> protrusions behind antennae.
sematura lunus.
character is also present in urania leilus and several geometrids
drepanoidea are monophyletic and associate with axia, epicimelia  and doa    as has also been suggested by recent molecular studies .
in several of the recent studies based on molecular data, drepanidae are the sister-group of the other macroheteroceran superfamilies.
in the analysis with the fastest evolving partitions removed and of the 268-taxon set, the clade comprising drepanoidea, cimeliidae and doidae are sister to other macroheterocerans.
however, when the partitions are not removed, the clade is nested in macroheterocera .
minet and scoble  <cit>  proposed a modification of the larval mandibles, a large, flat area that is ventrally delimited by a carina, to unite drepanidae and epicopeiidae.
the coding of this character was found to be continuous when other lepidoptera, especially macroheterocera, were examined in detail, and was thus impossible to code with any certainty, compromising its significance.
it was excluded from the present analysis.
the sister-group relationship of bombycoidea and lasiocampoidea  is recovered in the present result  as it has been in several studies based on genetic data .
in mutanen et al.
<cit>  and the ml analysis based on dna-data only of the present study , the sister-group relationship between the two superfamilies was not recovered suggesting that it is the addition of morphological data that brings about this relationship in the present result.
in the reduced taxon sets , with less taxa from both superfamilies, the sister relationship is not recovered either.
instead, lasiocampoidea are nested in bombycoidea.
in the ml tree of morphology only, lasiocampoidea, although not monophyletic, are in a clade with bombycoidea .
the adult character supporting the affinity of these superfamilies, also found in some other groups, is the closeness of the prescutal clefts in several taxa  <cit> .
the taxon sampling for bombycoidea in the present study is relatively low compared to the sampling of other large superfamilies.
this is because the inter-familial relationships and splitting into subfamilies within bombycoidea has been treated in more detail in recent phylogenetic studies  <cit> .
brahmaeidae swinhoe,  <dig> and lemoniidae neumoegen & dyar,  <dig> were synonymized by zwick  <cit>  based on evidence from dna and morphology.
in the present result, brahmaea does not show a close relationship to lemonia and dactyloceras, probably due to the lack of dna and adult character data.
the monophyly of noctuoidea is strongly supported  and constituent families are also recovered as monophyletic in analyses based on combined  and dna data .
in the analyses based on morphological data alone , noctuoidea are not monophyletic although they share the well-known synapomorphy, i.e.
the metathoracic tympanal organs.
the other established synapomorphy, the presence of two md setae on larval metathorax  <cit> , was excluded in the present analysis due to the impossibility to code it because many of ‘macrolepidopteran’ larvae are covered with microsetae that hide the possible presence of this character in this assemblage.
the monophyly of noctuoidea has convincingly been recovered in all recent molecular studies, and its family systematics with both morphological and molecular insight, has been elucidated by zahiri .
therefore, we do not discuss this superfamily further here.
discussion
our objective was to rigorously explore the morphology of ditrysian lepidoptera and to analyze the extensively coded data to infer a hypothesis for the phylogenetic relationships within ditrysia.
the aim was to obtain supplementary insight from morphology to that derived from dna-based data, and compare this information to the recently published dna-based phylogenetic hypotheses .
we consider this a way to supplement, evaluate and weigh these partially conflicting hypotheses.
a feature in common to all these dna-based studies is the weak or non-existing support for groups among the most basal and the most “derived” taxa, i.e.
for the so-called ‘backbone’ of the lepidopteran phylogeny.
for this reason, we put special emphasis on acquiring as comprehensive data as possible on the groups that supposedly belong to this area of lepidopteran phylogeny.
somewhat lesser attention was paid to inter-family relationships within superfamilies that have recently received specific attention  <cit> .
we combined morphological data obtained from larvae, pupae and adult specimens, coded from nearly the same sampling of exemplar taxa as in the study by mutanen et al.
<cit> .
the sampling covers nearly all recognized ditrysian families with a comprehensive representation of their subfamilies, and is supplemented with a number of taxa with unclear affinities.
as has been the result in dna-based studies, even comprehensive morphological data were not able to bring forth strong support for any particular hypothesis on the interrelationships of the backbone ditrysia.
possible reasons for this are the amount of homoplasy, general lack of informative characters except for apical nodes, conflicts between morphological and molecular data, and our decision to include and retain unstable taxa in the analyses.
these matters are discussed below.
lack of informative characters
problems in inferring deeper level relationships in ditrysia with molecular methods have been explained by ancient rapid radiation of most major lepidopteran lineages .
in ancient rapid radiation, lineage splitting occurs with close temporal spacing and very few character patterns are left as evidence of the relatedness of lineages that diverged from one another.
the little phylogenetic signal that remains may be obscured by biases in sequence data , and/or by homoplastic changes on the apical branches  <cit> .
rapid radiation and difficulties in inferring relationships among lineages united by short nodes is a problem in phylogenetic analyses of many insect groups that have undergone spectacular diversification in a short period of time in the cretaceous-tertiary  <cit> .
recent divergence time estimates also suggest a significant increase in the diversification rate in ditrysia around this time period,  <dig> million years ago  <cit> .
our examination of the morphology across the ditrysian superfamilies and families also confirms the lack of clear morphological characters for making inferences about these deeper evolutionary relationships.
a possible explanation for this is that the first taxa of each monophylum supposedly had, and their extant progeny may still retain, predominantly plesiomorphic traits compared to the more derived members of the clade.
paucity of apomorphies among the basal representatives thus leads a predictable pattern that evidence for inferring interrelationships of monophyla is scarce, if such ‘ancient’-looking taxa still exist and are included among the studied taxa.
the collected morphological data were found useful in reconstructing evolutionary relationships at the apical nodes.
our study affirms that most superfamilies are diagnosed by at least one fairly clear apomorphy, although diagnostic characters may not be found in all life stages.
however, the morphology-only analyses were unable to recover some of the well-established superfamilies.
an intriguing feature in our results is that even though no clear synapomorphies have been proposed to diagnose certain superfamilies nor were such characters found in the course of the present study or those suggested in earlier studies were refuted, such superfamilies were nevertheless also recovered in morphology-only analyses.
an example is the superfamily gracillarioidea.
although the exemplar taxa in the present study represent most of the ditrysian superfamilies, it has to be kept in mind that an unknown number of lineages have likely gone extinct in the millions of years that have followed the advent of ditrysia.
it is possible that along with the extinction of these lineages informative characters  have disappeared making the estimation of phylogenetic relationships among the ditrysian superfamilies even more difficult, just as in cases of insufficient taxon sampling where important phylogenetic information is missing.
at present, it is not known how heavily, for example, the cretaceous-paleogene mass extinction event affected the diversity of lepidoptera  <cit> .
the observed morphological homogeneity of ditrysia at the deeper level can safely be considered to be due to real biological and historical causes, but limitations on character coding are also set by methods used or material available.
some characters of potentially great influence may have been coded incorrectly for some taxa due to difficulties in making reliable observations.
many such characters were probably overlooked, and others had to be excluded from the final analysis.
examples include the presence or absence of the pupal frontoclypeal and epicranial sulci, which are easier to observe in exuviae than in pupae with uneclosed adults.
many setal characters of larvae were excluded from the majority of the macrolepidoptera, as their final-instar larvae are usually covered by small setae that obscure the pattern of primary setae.
this problem could have possibly been overcome if first-instar larvae could have been examined instead of final-instar larvae.
this would have, however, dramatically reduced sampling due to the lack of available material.
another option, using final instar larvae, supplemented with first instar larvae, would have been potentially misleading, as it would have been impossible to avoid duplication of some features in the data because of problems of homology.
a problem in coding adult characters is the high plasticity of structures and thus the difficulty of dividing characters into discrete states.
many characters are easily delimited when only a subset of the taxa is examined, i.e.
only extreme states are observed, but the boundaries become exceedingly blurred as more taxa are examined and more variation and intermediate forms are encountered.
many characters had to be discarded from the final analyses because of uncertain or artificial boundaries of states.
typically, e.g., the degree of sclerotization in various body parts, and variation in strength of sutures and sulci, proved impossible to code unambiguously across the taxa.
the use of morphometric tools could arguably increase the amount of character data as continuous characters such as length, size and details of shape could possibly be objectively divided into character states.
uncertainty over the homology of certain structures  across all groups was also a problem.
certain structures may be related to the size of the insect, but these are not easy to detect and interpret.
this was also a reason for our decision to omit wing venation and male genital characters entirely.
both morphological and molecular characters contain information about the phylogeny of organisms but the importance of morphological data in the future of phylogenetics has been disputed see e.g.
.
the problem with morphological data is less with quality than with quantity.
giribet  <cit>  discusses how morphology could be made “scalable” by utilizing new techniques that allow larger acquisition of morphological and anatomical data.
such new techniques include immunostaining and confocal laser microscopy, and also 3-d imaging techniques that allow reconstructions of the internal structures of organism without having to destroy them .
this enables the non-destructive study of rare museum specimens, ultrastructure of insects in amber fossils  <cit>  and in some cases live specimens  <cit> .
other new techniques to add to the list are magnetic resonance imaging   <cit> , advanced digital microscopy, computer-assisted tomography and image recognition  <cit> .
the tools to share and work on morphological information have also undergone a revolution, e.g.
morphbank : : biological imaging , morphobank  <cit>  and various discussion groups and web forums .
homoplasy
the frequency of homoplasy  was found to be high in both larval and pupal characters.
therefore, few characters of the immature stages, according to our analysis, are useful in diagnosing monophyletic groups.
not even the non-ditrysian vs. ditrysian grades could be distinguished by data coded from immature stage morphology alone.
also, the presumably strong pupal structures that according to previous hypotheses  <cit>  distinguish basal groups of ditrysia and obtectomera failed to recover these taxa unequivocally.
one possible explanation is the repeated shifts to different modes of life in ditrysia, as has been shown to be the pattern at lower phylogenetic levels too  <cit> .
extreme examples of taxa with significant amounts of convergence are leaf-mining larvae.
these larvae often have considerable convergent modifications in their morphology and reduction of structures.
parsimony analyses of morphology may fail to group taxa correctly if the level of homoplasy is very high  <cit> .
examples of such taxa in our data are lyonetiidae, douglasiidae, ochsenheimeria and millieria.
another case of obvious convergence is the tendency of diurnal lepidopterans in a number of clades to have clubbed antennae  and papilionoidea.
these examples demonstrate the potential effect of convergence on the overall topology, especially when there are several traits that by adaptation have led to appearance similar enough to make their recognition as distinct features hard.
in many cases more thorough study of the morphology of allegedly homologous characters, e.g.
“clubbed” antennae, is needed as new information could reveal them to be quite distinct despite their superficial similarity.
a morphological trait in lepidoptera regarded to be phylogenetically very significant is the larval proleg shape.
it has often been used as a distinguishing feature to separate ‘macrolepidoptera’ from ‘microlepidoptera’, but appears to be convergent in several lineages, viz.
schreckensteinioidea – urodoidea, zygaenidae, papilionoidea and macroheterocera, but excluding pyraloidea.
these groups usually have the following basic structure of the proleg: the base of the proleg is elongate and forms the greater part of the proleg; the proleg planta is asymmetric, laterally bulbous, and crochets are usually arranged as longitudinal mesoseries.
the usual ‘microlepidopteran’ type of proleg consists of a small proleg base that is little more than a ring encircling the planta.
the planta is either reduced or cylindrical, but not asymmetric, and the crochets most often are formed as a full or nearly full circle.
even though different combinations of these features were observed in some taxa, the “macrolepidopteran” combination is similar in papilionoidea s. lato and macroheterocera.
these groups have, as a rule, a similar larval mode of life as they are external leaf-feeders, usually gnawing on the edge of the leaf of host plants.
this favors adaptations of the proleg shape to optimize the grip on the edge of the leaf.
most pyraloid larvae live more or less concealed and do not need such specialized proleg modifications.
the structure of larval prolegs may play a crucial role in the topology obtained in analyses with morphological data, as the shape is not one feature, but a combination of several traits.
when a set of such non-independently evolved characters are coded from non-related taxa, they may have enough weight to pull the taxa together.
the high level of homoplasy in the morphological traits is in itself interesting.
just as dna can evolve fast enough to reach phylogenetic saturation, it has been shown that morphological characters can also show saturation  <cit> .
constraints exerted by intrinsic causes such as development, function or persistent selective trends do not allow for the continuous rise of new characters.
clades seem to exhaust their available character space and this leads to homoplasy.
however, convergence does not only concern the morphology of organisms.
genes can also be under selective pressure for a particular function and converge at the nucleotide level  <cit> .
although we aimed to avoid size-related characters in our morphological matrix, and tested the effect of removing characters we suspected could be correlated with size , characters or character combinations less obviously correlated with size may have remained.
as the signal for the deep nodes is weak, it is possible that such characters led to artificial clustering of groups by species’ overall size.
combining and comparing morphological and molecular evidence
even though morphological characters comprised only  <dig> to  <dig> % of the total amount of data, depending on the analysis, the incorporation of morphological data had, in places, significant effects on the results.
adding morphological data established more stable positions for some taxa, such as cyclotornidae and epipyropidae, which in molecular analyses have defied well-established positioning within their putative superfamily.
likewise, the monophyly of sesioidea  and cossoidea were now recovered, even though with entirely different morphological traits to those suggested earlier.
tineoidea excluding metapherna were recovered as monophyletic.
in dna-based studies metapherna, an alleged tineid not included in previous analyses, was never found to associate with tineoidea.
rather, it was recovered to form its own clade linking tineoidea and the more advanced lepidoptera – a result already hinted at by robinson and nielsen  <cit> .
it is not uncommon that morphology and genes tell different stories about the evolutionary history of organisms.
this is caused by the differences in the rates of change and the evolutionary mechanisms acting on these rates.
conflict between datasets and resulting phylogenies may be caused by convergent evolution, methodological reasons , or uncertain homology assessment.
cases of conflict between datasets are intriguing as they may allow the identification of the underlying biological processes causing the conflict  <cit> .
when conflicting clades arise in analyses, it is important to consider the effect of the methods chosen to analyze the data, but also to seek alternative evidence and information about the ecology of the organisms.
combined with the molecular dataset consisting of eight gene regions, the morphological data were able to distort not only weakly supported relationships, but also fairly strongly supported ones, such as the interrelationships in papilionoidea  <cit> .
without alternative sources of evidence, conflicting results would be difficult to evaluate.
recent studies based on different sets of genes to those in the present study provide some insight into the matter.
the taxon sampling of those studies has been, however, often lower and not permitting simple comparison.
also, the high support values for proposed relationships in low taxon sampling studies based on transcriptomic, genomic or proteomic data may be misleading.
in the absence of more closely related taxa, more distantly related taxa may be drawn together but still obtain relatively high support values  <cit> .
this phenomenon is illuminated by the conflicting hypotheses with strong support values in large molecular data sets of .
on the other hand, many patterns with low support observed in the study of mutanen et al.
<cit>  with dense taxon sampling, but lower amount of molecular data, have subsequently been obtained in other studies with much larger gene sampling.
conflicting patterns of relationship between molecular and morphological data, each with good support, is found, for example, in the position of douglasiidae.
the morphological result always places this family in gracillarioidea, as in the superfamily concept of davis and robinson  <cit> , also based on morphology.
however, molecular evidence is strongly in conflict with this result, unequivocally placing them in apoditrysia, which is also the case for our combined data analysis.
unstable taxa
in this study, we attempted to cover the diversity of ditrysian lepidoptera as comprehensively as possible.
therefore, we paid special attention to the inclusion of taxa with unknown phylogenetic positions, or which have been unstable in recent studies based on dna-data.
our aim was to examine whether the combination of extensive morphological data could provide information to infer their likely relationships.
we presumed that taxa appearing unstable in genomic data are not necessarily the same as those appearing as such in morphological data.
we further expected that in analyses with both types of data combined at least some of the rogues would become more stable.
combined data indeed gave moderate support at least in some analyses  for the inclusion of epipyropidae and cyclotornidae in zygaenoidea, a result that in molecular studies had either not been obtained, or only with a weak indication for epipyropidae.
the combination of data also brought out the sister-group relationship between bombycoidea and lasiocampoidea, not previously recovered with dna data similar to those used in this study.
although bootstrap support is not strong , the evidence supporting the monophyly of schreckensteinioidea + urodoidea seems a plausible hypothesis when morphological data are also considered.
similarly, when compared with results of studies based on molecular data only, the association of galacticoidea with tortricoidea also gains moderate to strong support in our analyses  for the first time.
a monophylum comprising alucitoidea, epermenioidea and carposinoidea is repeatedly recovered in the analyses of combined data, although with low support values.
adding morphological information also increased support for the placement of prodidactidae in hyblaeoidea.
on the other hand, there are cases where sequence data contradicts morphological information.
for instance, genetic data strongly suggests that douglasiidae are not gracillarioids, but apoditrysians.
nevertheless, combining data did not stabilize the position of many taxa that have been unstable in previous studies.
the unassigned heliocosma group, choreutidae and immoidea are typical examples of such taxa.
as the “backbone” of ditrysian lepidoptera did not gain much further support from combining morphology and molecular data, one can argue that several of the larger assemblages behave in an unstable manner, with presumably close groups displaying alternative arrangements among each other, e.g.
calliduloidea, papilionoidea and thyridoidea.
the cossoidea + sesioidea + zygaenoidea assemblage is particularly intriguing as it displays different positions and compositions among data sets.
in the results of morphology-only analyses, zygaenoidea were not unambiguously associated with cossoidea + sesioidea, which, in turn, often form a sister-group relationship with tineoid group families.
in the combined data with all third codon positions included, this three-superfamily assemblage was placed within obtectomera.
this relationship was broken when the third codon position nucleotides were removed.
an intriguing pattern in our data is that several of the unstable taxa are internal feeders.
these include cemiostominae, lyonetiinae, bedelliidae, heliodinidae , millieriidae and douglasiidae.
similar trends have been found in gelechioidea, where a wide array of life habits is found  <cit> .
although not tested using quantitative methods, these taxa are often observed to display disproportionally long branches and more often display multiple positions in phylograms.
this raises the question whether shifting to leaf-mining or other demanding life modes requires not only considerable morphological adaptation, but also rapid genetic evolution in many loci.
this would support the notion that taxa with long branches, i.e.
those with a disproportionally fast evolutionary rate, would have a higher likelihood to behave as rogues.
in such cases, adding more genetic data might not be able to solve the problem.
another example of taxa with an unusual mode of life are epipyropidae and cyclotornidae, whose larvae have a two-phased development, and are ectoparasites of other insects; a life-history trait that is possibly associated with rapid molecular evolution  <cit> .
as discussed above, stable placement of these families has not been reached with dna data.
CONCLUSIONS
our results support previous assertions of the paucity of clear morphological characters diagnosing larger assemblages of ditrysian superfamilies.
previously proposed characters are tenable to some extent but alone do not allow the unequivocal phylogenetic placement of taxa.
the degree of homoplasy is high and likely hid phylogenetic signal in the analyses.
nevertheless, most currently recognized superfamilies were recovered, and were found to have at least one previously recorded diagnostic character.
in some cases, the suggested apomorphies appeared to be supported only weakly, or remained entirely unsupported .
examination of the distribution of over  <dig> morphological traits across ditrysia revealed characters bearing evidence on the evolutionary relatedness among smaller groups, such as families  and in some cases superfamilies .
these characters also have potential use in identification keys and in identification of fossil samples.
the inclusion of morphological data also allowed to “pin down” several taxa that for some reason appear unstable in analyses of molecular data.
in the near future, we will begin to see the results of the first large-scale studies based on genome-wide data aimed at solving the deep divergences in the evolutionary tree of ditrysia.
the morphological data collected in the course of this study will certainly prove useful in evaluating and validating the results.
likewise, the results of future phylogenomic studies will help explain the distribution of some of the morphological characters and identify possible misinterpretations.
comparing structures across the ditrysian subfamilies will also gain us insight into the evolution and diversity of morphological characters.
the distribution of these characters and their reversals across the different clades can be used in future studies on the evolution of these traits.
availability of supporting data
the data sets supporting the results of this article are available as additional files  <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig> and  <dig>  sequence data are available in genbank .
morphological character data are available in morphobank, project  <dig>  http://www.morphobank.org/.
additional files
additional file 1:
list of exemplar taxa in study.
additional file 2:
sequence lengths.
additional file 3:
gene summary.
additional file 4:
genbank accession numbers.
additional file 5:
character list.
additional file 6:
character matrix.
additional file 7:
phylogenetic tree from maximum likelihood analysis of combined data with third codon position removed from molecular data except from ef1a.
<dig> taxa,  <dig> characters .
additional file 8:
phylogenetic tree from maximum likelihood analysis of dna data with third codon position retained,  <dig> taxa,  <dig> bp.
additional file 9:
phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data,  <dig> taxa.
molecular data were categorized according to the rate of evolution with the program tiger and partitioned with ratepartitions.
two fastest evolving partitions  removed.
additional file 10:
one of the  <dig> equally most parsimonious trees  from the parsimony analysis of morphological data,  <dig> taxa,  <dig> characters.
additional file  <dig> lists morphological data for  <dig> taxa, but in the present analysis, data for two species of thereutis were concatenated into one terminal taxon as were those of two species of cyclotorna making the total number of taxa in the analysis  <dig>
additional file 11:
phylogenetic tree from maximum likelihood analysis of morphology,  <dig> taxa,  <dig> characters.
additional file  <dig> lists morphological data for  <dig> taxa, but in the present analysis, data for two species of thereutis were concatenated into one terminal taxon as were those of two species of cyclotorna making the total number of taxa  <dig>
additional file 12:
phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data,  <dig> taxa.
the data set includes all the taxa for which both morphological and molecular data were available.
additional file 13:
phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data,  <dig> taxa.
the data set includes all the taxa for which dna data were available.
morphological data were available for  <dig> of these taxa.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
lk and ps had the original idea for the study.
the study was designed by all authors.
mh, lk and ps carried out the preparation and morphological examination of the samples.
molecular data were prepared by mm and nw.
mh, nw and lk carried out the phylogenetic analyses.
all authors contributed to writing the manuscript.
all authors read and approved the final manuscript.
we would like to express our deepest gratitude to the following people and collections for loan of material: leif aarvik, john w. brown & m. alma solis , marc epstein, nico elfferich, marianne horak & ted edwards , ole karsholt , jari kaitila, joël minet, kari nupponen, teemu klemetti, leptree project and people, texas a&m university and many other members of the international community of lepidopterists.
laura törmälä is thanked for laboratory work.
we are also thankful to jurate de prins, joaquin baixeras, and toomas tammaru for their comments on an earlier version of this manuscript.
library secretary of the zoological unit  marjatta mikkonen was of great help in the search for literature.
marianna teräväinen corrected the language.
we would also like to thank charles mitter and an anonymous reviewer for comments that greatly improved the manuscript.
this research was funded by the finnish cultural foundation, oskar Öflund foundation, ella and georg ehrnrooth foundation, and the university of helsinki ; academy of finland, grant # <dig> .

 five years ago, the backbone branching pattern of the evolutionary tree for one of the largest groups of insects, moths and butterflies, was still for the most part unresolved. our state of knowledge of phylogenetic affinities within this insect order, with over described species, was well depicted by the comb-like summary tree composed by kristensen and skalski. well-supported hypotheses existed only for the relationships among the most ancient superfamilies. only very cautious and speculative hypotheses had been proposed for the relationships among the superfamilies of the enormous clade of more advanced lepidopterans, known as the ditrysia. this ditrysian clade encompasses nearly % of all butterflies and moths. ditrysia are currently divided into superfamilies. the size of the superfamilies varies considerably in the number of described species, from several small monotypic superfamilies based on a single genus with only a single or a few species to noctuoidea with over species. some of these monotypic superfamilies are characterized by features that do not correspond to those diagnosing any of the other superfamilies. the current circumscriptions of the ditrysian superfamilies are based on morphological features and have been generally accepted. recent studies based on genetic data have corroborated the monophyly of most of the superfamilies. however, the interrelationships among most of these superfamilies lack any convincing supporting evidence. this lack of knowledge of evolutionary affinities has hindered further research on this enormous suborder.fig. current classification of ditrysia, modified from van nieukerken et al. the initial challenges to inferring evolutionary relationships among the ditrysian superfamilies have been attributed to their morphological homogeneity, and thus to the difficulty of finding distinct, unique morphological characters uniting superfamilies and supporting descent from a common ancestor. untold numbers of detailed studies on the morphology of lepidoptera, especially of genital structures and wing venation, have furthered our understanding of the relatedness of species, genera and families, but the relationships among larger assemblages have remained obscure on morphological grounds. the most important studies in comparative morphology aimed at resolving the relationships among ditrysian lepidoptera were by brock and minet. the works by minet have been influential and the basis for the division of ditrysia into three nested clades, apoditrysia, obtectomera and macrolepidoptera. these hypotheses are, however, based on a limited number of characters; two for apoditrysia, two for obtectomera, and one for macrolepidoptera. furthermore, they have not been verified across all ditrysian superfamilies, families or subfamilies. in recent years, several studies based on molecular methods have brought new information to bear on the evolutionary history among and within the ditrysian superfamilies. these studies range from wide-taxon sampling multigene studies to studies with fewer taxa but with large-scale gene sampling or genomic data. the impact of these studies on our knowledge of the relatedness of the lepidopteran superfamilies is beyond question. they have supported some of the previous hypotheses of relationships of ditrysian lepidoptera, but have challenged or refuted others, e.g. the monophyly of tineoidea or the position of papilionoidea. the main division of ditrysia into the nested clades apoditrysia, obtectomera and macrolepidoptera has been backed to some extent, but with changes. for example, genetic evidence has supported the morphology-based suspicion of rammert and kaila that gelechioidea are a superfamily in obtectomera to which they were not thought to belong. also, there is mounting evidence that butterflies, skippers and hedylids are more closely related to ‘microlepidopteran’ groups of obtectomeran moths than to the core macrolepidoptera. the traditional concept of the superfamily papilionoidea has also changed. the monophyly of tineoidea has repeatedly been found to be unfounded, and the monophyly of zygaenoidea has also been challenged as cyclotornidae and epipyropidae have not been unambiguously included. despite these advances, the evolutionary history of the ditrysian lepidoptera is still full of questions. the low support values obtained in many studies, especially along the apoditrysian “backbone”, have made the suggested relationships uncertain. a proposed reason for the difficulties in inferring evolutionary relationships among the ditrysian superfamilies is their postulated explosive diversification in the mid- to late cretaceous into the early tertiary, contemporaneous with the radiation of flowering plants. when diversification is rapid, little molecular or morphological change has time to accumulate as evidence of common ancestry. this weak signal may also be obscured by subsequent changes in the dna and morphology. with a limited spectrum of genes it may not be possible to extract this weak signal. to overcome this, phylogeneticists, in lepidopterology as elsewhere, have turned to comparative genomics. by comparing large portions of or entire genomes, clarification of the deeper evolutionary relationships is hoped to be achieved. recent phylogenomic studies have indeed shown much stronger support values for relationships among the superfamilies than the first studies based on individual gene sequences. however, low taxon sampling still is an issue in these studies. with low taxon sampling, distantly related terminal taxa may group together with strong support values in the absence of more closely related taxa. a wider taxon sampling and more refined models of molecular evolution are needed to overcome this problem. conflicting hypotheses with strong support values are hard to evaluate if alternative evidence is not available. such evidence is traditionally sought in the morphology of the organisms. synapomorphies, i.e. shared derived characters, are evidence of descent from a most recent common ancestor and inform us of a common evolutionary history. however, as morphology has never been systematically studied across all ditrysian groups, and the data subjected to a rigorous phylogenetic analysis, we lack morphological information of their affinities. the alleged morphological homogeneity of ditrysia, the laboriousness of comparative morphology, and also the immense number of species of the order have without doubt held back their classification. several earlier hypotheses based on morphological features that claim to diagnose certain groups also remain untested across larger groups. a frequent problem in phylogenetic analyses is taxa that for some reason do not find a stable phylogenetic position and may even have a nearly random placement. such taxa are often called rogue taxa. the presence of such taxa may effectively remove resolution of phylogenetic trees. they have also been shown to affect the topology of the entire tree, thus distorting the phylogenetic relationships among other clades. rogue taxa exert a great influence on the statistical support values of the nodes in phylogenetic trees, obscuring the meaning of these values. increasing genetic data may not necessarily solve the problem, as fast-evolving taxa may display deviating patterns across large parts of their genomes. exclusion of rogues is indeed a common practice to try to salvage reliable-looking topologies for the remaining ‘well-behaved’ taxa. moreover, it is likely that a lot of such pruning goes unreported as it is done during analysis, long before manuscript preparation. exclusion of rogues, however, means discarding data, and is therefore not an optimal way to treat these taxa, which may be crucial for the study in question. taxa may behave as rogues in molecular or morphological analyses, or both. however, there is evidence that by combining morphological and molecular data, the negative effect of rogues is reduced and these taxa may find a stable position. thorough examination of the morphology across all ditrysian groups is also needed to obtain information on character evolution. in a phylogenetic context, morphological structures may show an order of evolution – a polarity. some characters may also explain the success of certain groups of lepidoptera. such characters are also called key innovations as they enable, for example, the exploitation of resources, help avoid predators or survive adverse conditions. the time of emergence of such characters can also be investigated to assess whether certain factors in the geologic history of the earth have favored lepidoptera with particular attributes. fossil lepidoptera are extremely hard to identify and assign to particular taxa, because in adults the scale covering hides diagnostic characters, and fossils of soft-bodied larvae are rare. however, a better understanding of the morphological diversity of lepidoptera and of characteristic features of different groups could prove useful for reliable identification of fossil lepidoptera. the time of appearance of a certain character is also of interest. dated and securely identified fossils are used as calibration points in dating phylogenies to obtain estimates of times of divergence. to address these needs, we examine the morphology of ditrysian lepidoptera across of the currently recognized superfamilies with the aim of providing phylogenetic information on the relationships among them. we examine morphological characters from larvae, pupae and adult males and females of over exemplar taxa. we aimed to obtain samples as comprehensively as possible across lepidopteran subfamilies, and added a number of taxa with unknown or unsupported placement. special emphasis was placed in the comprehensive inclusion of non-macrolepidoptera groups and character sampling of all life stages. we have aimed to examine real samples of as many exemplars as possible ourselves. these data have been supplemented and cross-checked using literature sources. we analyze the morphological data alone and combined with the molecular data set published by mutanen et al., supplemented with additional taxa and sequences for some taxa. we report the results of the phylogenetic analyses and compare these to the results of other phylogenetic studies on ditrysia. we re-examine previous hypotheses based on morphology and present new morphological evidence supporting affinities among clades. we also discuss the effects of homoplasy and methodological issues on the results. methods taxa the exemplar taxa for the study were chosen to represent as many ditrysian subfamilies as possible. we estimate that about % of the subfamilies recognized by kristensen et al. are included in the present study. an exact number is impossible to give as the concepts of many subfamilies vary or are ambiguously defined. the full dataset consists of exemplar species. for taxa there are only morphological data, for taxa only dna data, and for terminal taxa there are both morphological and molecular data. in some cases either the life stages or sequence data were not available from the same species, and so were supplemented by data obtained from a closely related species in the same genus. these cases are indicated in additional file in addition, non-ditrysian taxa were included in the analyses. all trees were rooted using micropterix calthella. the authors examined the morphology of eight of the non-ditrysian outgroup taxa: eriocrania semipurpurella, hepialus humuli, andesiana lamellata, lampronia capitella, incurvaria pectinea, azaleodes micronipha, ptyssoptera sp. and tischeria ekebladella. the remaining taxa were represented by only molecular data. the molecular data are largely the same as in mutanen et al., but excluding polypogon strigilatus, and supplemented by additional taxa and sequences. these additional taxa are taxa for which fresh material had become available since the publication by mutanen et al. or other relevant taxa and data published in previous phylogenetic studies. examination of morphological characters morphological data were collected from species. larval characters were available and coded for of these species, pupae for species, and adult characters for species. when all life stages of same species were not available, character coding was done from a closely related species in the same genus for which material was available. there are seven such cases in our morphological data, indicated in additional file larvae examined were preserved either as dry inflated or in alcohol. pupae or pupal exuviae were also examined from either dry or alcohol preserved samples, or whenever possible, both. adult characters were coded from mounted specimens. the wings of the adults were removed and the body treated in % koh solution to clear it of scales, lipids and muscle tissue. the exoskeleton of the adult was kept in alcohol. the specimens were obtained mostly from the collections of the finnish museum of natural history, but also from several other museums and a large network of international collaborators. larvae and pupae were examined for external characters. due to the paucity of material of the potentially most informative first instar larvae, the larval characters were coded from final instars. adults were coded mostly for characters of the exoskeleton. characters known to have much intraspecific variation and/or considerable problems with assessment of homology among more distantly related groups, such as genital structures and wing venation, were not included in the present matrix. wing venation characters, except for the distinction between homoneurous and heteroneurous venation, were left out from the analyses. although homologous characters of wing venation are fairly easy to establish for many of the larger taxa, even they posed considerable problems with homologies and continuous nature of variation. these problems were even more pronounced with many groups of smaller moths for which homologies, usually losses of veins and considerable intraspecific variation, were impossible to establish. inclusion of venation characters would have considerably increased the amount of uncertain homologies, thus more likely decreasing rather than increasing the reliability of our morphological data. the head was not opened. the basal wing sclerites were not studied because the observation of the sclerites in the small-sized exemplars would have needed very time-consuming and delicate preparation of special slides. preliminary exploration of these structures indicated that features of these structures would largely prove continuous, thus a thorough study of them was deemed implausible in the time-frame allocated to the study. as the focus of the study was to seek morphological evidence to infer phylogenetic relationships among the superfamilies, less attention was paid to intra-family level relationships. some apparent autapomorphies were coded when they were known to be present in other members of the taxon in question. literature from which characters for some of the taxa were coded is listed below the character list in additional file many of the characters used in this study are not original but were initially taken from literature, and examined across the whole taxon sampling. character observations were made solely using light microscopy. a scanning electron microscope would have greatly expanded the number and probably also accuracy of observable characters, but was not available. the observations were made with leica mz mz, leitz diaplan phase contrast microscope and wild m character data management and storage was implemented in morphobank. the morphological character matrix comprises binary or multistate characters. of these are larval characters, are pupal characters, and are adult characters. both inapplicable and unclear character states were coded with question marks. continuous characters were not included in the study. those typically include characters showing high plasticity, characters based on relative variation in size, shape, or degree of sclerotization, all these being factors that complicate assigning characters into discrete character states. our focus was to find unambiguous characters suitable as potential diagnostic characters in identification keys. molecular data genetic data, totaling base pairs, were sequenced from taxa. dna samples were obtained mostly from museum specimens and additional specimens collected during field work in in new zealand and tasmania. genetic data includes cytochrome oxidase subunit i gene from the mitochondrial genome and seven genes from the nuclear genome; elongation factor-1a, ribosomal protein s, carbamoyl phosphate synthetase domain protein, cytosolic malate dehydrogenase, isocitrate dehydrogenase, wingless, and glyceraldehyde-3-phosphate dehydrogenase. sequence lengths, gene summaries and genbank accession numbers are presented in additional files and respectively. dna was extracted from legs, sometimes from other body parts using qiagen’s dneasy extraction kit. dna amplification and sequencing largely followed the protocol explained in wahlberg and wheat, but with slight modification and optimization applied to clean-ups of pcr reactions and a different sequencing reaction purification kit applied later. sequencing was performed mainly with an abi capillary sequencer, but a smaller part was performed with an abi prismr 3130l capillary sequencer. details on dna extraction are explained in mutanen et al.. the program voseq was used to generate sequence data files and gene summaries. matrices were prepared for concatenation in mesquite. analyses the parsimony analysis of morphological data was done with tnt tree search using new technology search with the following settings; all four search options selected; find minimum length times; initial addseq 15; random seed ten. morphological characters were treated as unordered and equally weighted. inapplicable, unclear and missing characters were coded with question marks. for examination of trees and character distributions the maximum parsimony trees were exported to winclada, where they were treated as follows: select all trees; keep best trees only. data sets of pupal, larval characters and adult characters were examined separately and combined. morphological data were also analyzed with maximum likelihood methods carried out with raxml in cipres using the mk model for morphological data. dna data and combined datasets were analyzed using maximum likelihood methods carried out with raxml in cipres. analyses of combined data based on three different data sets were run:data set of taxa including only taxa for which both dna and morphological data were available; data set of taxa including all the taxa for which dna data were available, for of these taxa also morphological data were available and were included, and; data set including all taxa: of these having both molecular and morphological data, taxa with only dna data and taxa with only morphological data. analyses were run with all codon positions retained and with third positions excluded from all genes except ef-1a, which evolves more slowly than any of the other genes. this methodology was also followed in mutanen et al.. in another set of analyses, both molecular and morphological data were categorized and partitioned according to the rate of evolution of the characters with the programs tiger and ratepartitions. molecular data were divided into seven partitions, morphological data into five. the analyses were run with all partitions present but also with the partitions with the fastest evolving molecular characters removed. omission of further partitions began to break up well-established clades and produce spurious groupings of taxa. the fastest evolving morphological character partition was also omitted in some trial analyses. however, it became clear that the fastest evolving morphological character partition also contained phylogenetically informative characters and its removal caused false disruptions of even closely related taxa, and was therefore retained in subsequent analyses. the gtr + g model was used for the dna data and the mk model for morphological data. the ascertainment bias correction option was selected for partitions of morphological data, which include only variable characters. the option to conduct bootstrapping and search for the best-scoring ml-tree in a single program run was selected. supports for nodes were evaluated with 100– bootstrap replicates. the resulting trees were examined in figtree. mesquite was used to explore character distribution on the trees based on combined data. we compared our result to other recent studies of ditrysian phylogeny based on molecular data. comparison of the results with these studies was not straightforward as many of the other studies have much lower taxon sampling and several superfamilies included in the present study are not represented. comparison of support values obtained for certain clades in different studies is also problematic. we have deliberately retained unstable taxa in our analyses although their presence likely lowers support values, and, especially in consensus trees obtained by parsimony analyses, their presence significantly collapses resolution. to overcome this, and more effectively filter information from the mp trees, we explored them individually. the character and taxon sampling of the present study were not specifically designed for researching the relationships among subfamilies and families, but we do report on these if they significantly deviate from current concepts. relationships among non-ditrysian lepidoptera were not explored. many of the morphological characters supporting subfamily, family or superfamily level clades have been found by previous authors and we do not repeat all of them. more detailed information on characters and corresponding references are found in additional file  we report the results of the maximum likelihood and parsimony analyses of the morphological dataset alone, and compare these results to those of the ml analysis of combined morphological and molecular data. we consider the total evidence ml analysis our main result. this decision as to the priority of results is based on the principle of inclusion of the maximum amount of data, even though we recognize that this may also bias the results, as may, for example, the inclusion of unstable taxa. moreover, as källersjö et al. have shown, third codon positions, even if possibly often saturated, also contain phylogenetic signal that may be informative when combined with other data. we also compare the results with those obtained with third codon positions excluded and the “noisiest” molecular character partitions i.e. those with the fastest evolving characters removed. the fastest evolving sites were identified with the programs tiger, which identifies rapidly-evolving characters and ratepartitions, which partitions data according to their rate of evolution. with the fastest evolving molecular character partitions removed, the results were nearly identical to those obtained with all data retained but produced some spurious groupings, which contradict results of other recent studies. in the results of analyses based on data with the noisiest partitions removed, bootstrap support values for the clades were in general lower or more or less equal to those obtained for the same clades in the analyses based on data with noisy characters retained. the results of the analyses based on dna data alone are not discussed separately as they are very similar to those in mutanen et al., discussed therein.fig. phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data; taxa; characters. ***both dna and morphological data; **only morphological data; *only dna data fig. phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data; taxa; characters. ***both dna and morphological data; **only morphological data; *only dna data fig. phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data; taxa; characters. ***both dna and morphological data; **only morphological data; *only dna data analyses based on morphological data only the parsimony analysis yielded equally most parsimonious trees. the resolution of the resulting trees was generally uniform, but interrelationships of certain superfamilies varied among the trees. the strict consensus tree is very poorly resolved due to the adverse effect of some taxa assuming two alternative positions, as well as the paraphyly of gelechioidea observed in some of the mp trees. an example of one of the most parsimonious trees is shown in additional file most superfamilies are recovered as monophyletic. however, some results appear spurious, such as the consistent grouping of the non-ditrysian hepialus in tineoidea. bombycoidea, and in some mp trees also drepanoidea and geometroidea, were nested within noctuoidea. the unstable position of ochsenheimeria is likely due to the internal feeding mode of its larva having caused morphological modifications that remove it from the otherwise constantly recovered association with ypsolopha, the larvae of which feed externally on leaves of various trees and shrubs. other differences in the topology were minor rearrangements of apical genera. the tree resulting from the maximum likelihood analysis of the morphological data set is shown in additional file there is no support for the deeper nodes, but most superfamilies are recovered as monophyletic except for those in macroheterocera, which are intermixed. separate analyses on data sets of larval, pupal and adult morphology were also run, but are not discussed, because the resulting trees were nearly unresolved due to the small number of characters in comparison to the number of taxa included in these analyses. combined analyses phylogenetic relationships the results are in many respects similar to those obtained in the study by mutanen et al. as the molecular data are largely the same, especially when third codon positions were excluded. considering that the morphological data in the present study consisted of only % of the total amount of data, the phylogenetic signal from these characters appears disproportionately strong, and their inclusion significantly affects the topology of several clades, while supporting the molecular results elsewhere. as will be discussed in more detail below, morphological attributes stabilized some taxa that had been unstable in dna-only analyses and thus propose a hypothesis about their phylogenetic affinities. likewise, molecular data anchored taxa that were unstable in our analyses of morphological data alone. clades supported by both morphological and dna-data in general obtained higher support values in the combined analyses than when either type of data was analyzed alone. however, we also obtained some intriguing patterns in the topology that strongly conflicted with the results of several recent studies based solely on genetic data. these patterns show the strength of morphological data over the eight gene regions and the usefulness of different types of data in evaluating different phylogenetic hypotheses. the overall topology of the phylogenetic trees based on the three different data sets is very similar and the weakly supported backbone nodes are weak in all trees. although the support values are the highest for the smallest data set including taxa, and weakest in the trees including all taxa, we choose to present the latter. the topology of the trees is essentially the same and the tree based on the larger taxon sampling also includes interesting information on the placement and performance of taxa with only either morphological or molecular data. we report support values obtained from all three analyses. phylogeny of main ditrysian clades ditrysia the monophyly of ditrysia is recovered in the maximum likelihood analyses based on dna, morphological and combined data sets. the support values for the ditrysian clade in the ml analyses of combined data are as follows: taxa: %; taxa: %; taxa: %. in the parsimony analysis of morphology alone, hepialus is placed among ditrysia. probable reasons for this are that in immature stages there appears to be no morphological characterization that would distinguish ditrysia from monotrysians, and the methodological decision to exclude morphological characters only applicable for resolving interrelationships within monotrysia. the adult characters supporting the ditrysian clade are the name-giving ditrysian female reproductive system i.e. one opening for mating and a separate one for laying eggs, and the absence of structures found in some but not all the outgroup taxa such as the division of sternum an elongated metathoracic trochantin, and the tergosternal connection of tergal origin. division of ditrysia our results support previous assertions on the scarcity of clear morphological characters that would define larger assemblages of ditrysian superfamilies, and no new such characters were found. analysis of the morphological data alone or combined with the molecular data set did not improve the resolution of the deeper nodes obtained by mutanen et al.. support values for several of the deeper nodes, but also of many other nodes, are low. this is in part due to our deliberate decision to not remove unstable taxa from the analyses. the unstable taxa varied depending on the data set and method used. in the maximum likelihood analyses of combined data, the following taxa were unstable: imma, millieria, douglasiidae, tanaoctena, cadmogenes and several small superfamilies. there were also other taxa with unstable positions, but only within superfamilies, e.g. lyonetia in yponomeutoidea. characters previously proposed by minet to circumscribe subclades of ditrysia were tested across the 300+ taxa representing out of superfamilies. the results show that although several obviously contain phylogenetic signal, there is ambiguity in the interpretation and distribution of details of these characters, as already implied by minet. apoditrysia characters of the adult sternum have been proposed to define the apoditrysia, i.e. all ditrysia except tineoidea, gracillarioidea, yponomeutoidea and gelechioidea. recent studies have, however, repeatedly suggested that gelechioidea belong to obtectomera, a clade nested in apoditrysia. apodemes of sternite of apoditrysia were defined by minet to be short, posteriorly enlarged and with anterior corners of sternite distinctly produced laterally. the division between non-apoditrysian and apoditrysian sternum types is also described as the tortricid type versus the tineid type of thoraco-abdominal articulation. the apoditrysian type of sternum is indeed never present in the tineoidea, gracillarioidea and yponomeutoidea, and present in most apoditrysians. however, a condition more similar to the non-apoditrysians is found in some cossoidea, zygaenoidea, millieriidae, gelechioidea, and in the present study, it is observed that clear lateral extensions are occasionally absent even in some macroheterocera, a clade nested in obtectomera. in some taxa, sternite is modified due to the presence of hearing organs making the presence or absence of the lateral extensions difficult to assess. in mimallonoidea, the presence of the character was also ambiguous. sternum has an additional modification that, according to minet, is an autapomorphy for the earlier concept of drepanoidea: the lateral extension of sternum is long and reaches in front of the spiracle all the way to the anterior corner of sternite in addition to drepanidae and epicopeiidae, this modification is present in pterophoroidea, schreckensteinia, coleophora in gelechioidea, carposinoidea, crambidae, calliduloidea and papilionoidea. obtectomera a pupal structure – body appendages firmly glued to the body – has been considered a character of considerable phylogenetic significance. a further specialization, i.e. pupa with immobile intersegments between abdominal segments and, was used by minet, building upon the work of chapman, to define the clade obtectomera comprising pyraloids, butterflies and other macrolepidoptera. in some groups the abdominal segments 1– are variably fused, but the mobility of 3– is clear-cut. therefore the latter was taken into account in character coding. presence of this pupal type in yponomeutoidea, gelechioidea, epermenioidea and alucitoidea was interpreted to represent convergence. now it appears that carposinoidea, epermenioidea, alucitoidea and gelechioidea belong to obtectomera sensu minet, so convergence needs only to be assumed for some yponomeutoidea and possibly for the unstable immoidea. in the present study, this interpretation is supported when third codon positions or partitions and are removed. with third codon positions included, the cossoidea + sesioidea + zygaenoidea assemblage is nested in obtectomera, which, on a morphological basis, seems less plausible than the conventional view of them being true non-obtectomerans. the case of yponomeutoidea is not quite straightforward. the majority of yponomeutoids studied was observed to have the pupa with intersegment 3– not entirely either mobile or immobile, but rather having a limited ability for movement. this condition was considered its own character state, supported by experiments where pupal exuviae were softened so that the potential for mobility could be evaluated. several yponomeutoids with their pupal intersegment 3– found entirely immobile are taxa whose thoracic structures largely cover the ventral surface of the abdomen, which alone may effectively prohibit the ability to move abdominal segments and lead to the fusion of the thoracic segments. immoidea, with an obtect pupa, was associated with choreutidae, i.e. non-obtectomerans, although with negligible support. the pupal skin in immoidea is thin, and the level of intersegmental mobility is ambiguous. minet considered their pupa obtectomeran while common interpreted it to be of the non-obtectomeran type. another character later proposed by minet to redefine the obtectomera is the presence of a dorsal lobe or protrusion on the adult pretarsus. as can be noted also from the pictures in minet, the presence of the setose lobe on the pulvillus is not always obvious. the lobe can be small and not easily distinguishable from the non-obtectomeran type of pulvillus, especially if the latter has long setae on the dorsal side. this character was also examined across all taxa in the present study, and the different character states applied here were: pulvillus simple, pulvillus with setose outgrowth, or pulvillus bifid. setose outgrowths were clearly visible at the base of the dorsal side of the pulvillus in most taxa assigned to obtectomera. in gelechioidea, however, the pulvilli were simple. in a few bombycoids and in pieridae, the pulvillus was also more of the simple type without projecting setae. in other papilionoids, the pulvillus was distinctly bifid and differed from those with a protrusion with setae. imma and copromorpha also had distinct setose lobes on their pulvilli. the resolution limits imposed by the light microscopes used in the present study lowered confidence in coding this character for some of the smaller taxa. in conclusion, our result suggests that even though the obtectomera is probably largely a valid clade, there are several superfamilies whose assignment to this clade is uncertain, and the degree of homoplasy appears higher than generally presumed. fänger also examined this character across over ditrysian families with a scanning electron microscope and confirmed its validity with some exceptions. macrolepidoptera and macroheterocera the apomorphy proposed to delimit macrolepidoptera, the first axillary sclerite of the forewing with an elongate angle, was not included in the present matrix due to the continuous nature of the character, which compromises its usefulness. the traditional concept of macrolepidoptera has recently changed as molecular analyses have repeatedly suggested that pyraloidea appears to be more close to the core macrolepidoptera than papilionoidea and calliduloidea. the clade including the rest of the traditional macrolepidoptera is referred to as macroheterocera. in several recent molecular analyses, mimallonoidea has been placed as the sister-group of macroheterocera. this is also the result obtained in our analysis with the three codon positions included. in analyses without the third codon positions or fast evolving partitions, the pattern changes and the arrangement of clades varies, with all alternatives having low support. although the macroheteroceran clade is recovered in the combined analysis, no unique morphological character was found to unite them. interrelationships and composition of ditrysian superfamilies the most basal clade of ditrysia recovered in the combined analysis is a monophyletic tineoidea with the following support values: 268: 88; 422: and 473: %. they are also monophyletic in the result based on dna with all three codon positions retained and those of combined data with the fastest evolving partitions removed. tineoidea are divided into two clades. one clade includes the monophyletic families eriocottidae and psychidae, the other tineidae. in recent phylogenetic studies based on dna data alone, a monophyletic tineoidea has only been obtained in the studies by regier et al. and in analyses of unaltered nucleotides from all three nucleotide positions. in the present study, when the third codon positions are removed, tineoidea are paraphyletic, a result similar to that obtained in other recent studies based on sequence data. tineoidea are then paraphyletic, with eriocottidae, psychidae, and tineidae as separate clades. the present study did not include exemplar species of meessiidae, included in the phylogenetic study on tineoidea by regier et al.. in the results of their study, representatives of meessiidae were excluded from all previously recognized tineoid families and placed as sister to all other tineoidea + ditrysia, thus also contributing to the paraphyly of the superfamily. had meessiidae been included in the present study, the monophyly of the superfamily could have been challenged. tineoidea were also recovered as monophyletic in the maximum likelihood analysis of morphological data, but not in the parsimony analysis. adults in all three tineoid families included in the present taxon sampling are characterized by the presence of a distinctly shortened proboscis with loose galeae and pseudapophyses found in most females. shortening or absence of galeae is, however, found throughout the ditrysian superfamilies. the looseness of galeae could not be coded because most of the short galeae in other groups also loosened when treated in koh solution, and on non-denuded specimens the scale coverage effectively prohibited the reliable observation of this character. light microscopes were insufficient for comparison of differences in the microstructures holding galeae together. no support for the monophyly of tineoidea was found from larval or pupal characters. within tineoidea, monophyly of tineidae is supported by a unique synapomorphy, a nearly triangular cap dorsally covering the base of the larval antenna. dryadaula, recently proposed to be a new family, dryadaulidae regier et al., is nested in tineidae, possibly because only morphological data from the adult male and female were available and the larval synapomorphy was not studied. psychidae larvae share the median fusion in all thoracic legs; the members of the subfamily oiketicinae are characterized by the horizontal position of the larval spiracles of t fig. unique tineid apomorphy: more or less triangular cap dorsally covering base of larval antenna. morophaga choragella metapherna robinson & nielsen is an australian genus that has been placed in myrmecozelinae in tineidae, although robinson & nielsen noted that it has, e.g., an unusually long proboscis for a tineoid. in the analyses of dna data-only and of combined dna and morphological data, metapherna salsa is placed between tineoidea and the rest of ditrysia. in the morphology-only ml analysis metapherna has an unresolved position within non-apoditrysians. in the mp analysis it is placed either as a basal ditrysian or in yponomeutoidea, but never in the tineoidea. based on genetic evidence and their divergent morphology, species of the metapherna group obviously may eventually require their own family. in the analyses of combined data, metapherna is basal to a large clade divided into a rather weakly supported clade including taxa assigned to yponomeutoidea and gracillarioidea, and the apoditrysian clade. yponomeutoidea plus gracillarioidea, with douglasiidae excluded, have also been found to be monophyletic in several other recent phylogenetic studies on the phylogeny of ditrysia, and also in two studies focusing on the phylogenetic affinities within these superfamilies. in analyses based on morphology only, douglasiidae was included in gracillarioidea, as has also been formerly suggested. however, in our preferred result it is associated, with uncertain position, among apoditrysia, as in those recent molecular studies in which it has been included. in the present study, we did not find any unique characters common to all members of gracillarioidea and yponomeutoidea. this grouping was supported by a combination of homoplasious characters such as characters of antennal scaling and the tendency to have the proximal pair of tibial spurs of the hindleg on the proximal half of the leg. however, the transverse costa at the base of sternum a found in most yponomeutoidea and proposed to be a possible groundplan autapomorphy for the superfamily was also found to be present in roeslerstammiidae: roeslerstammia, vanicela, amphithera and the female of thereutis. a similar costa, although longer and often reaching the sternal apodemes, was present on sternum of choreutoidea. some tortricids also exhibited a sclerotized bar or patch in a similar position.fig. transverse costa behind anterior margin of sternum ii. a yponomeuta evonymellus; b roeslerstammia erxlebella yponomeutoidea were found to be monophyletic and the presence of the pleural lobes on sternum of the males was verified. pleural lobes were found to be present in nearly all yponomeutoids. in some cases, the expansion of the lobe was not strong and it was difficult to determine if lobes were present. they were absent in cadmogenes literata, which is currently placed in plutellidae, but may eventually need to be placed in its own family. in our study, only morphological data were available for cadmogenes. its position was unstable, and it grouped most often with some other unstable taxa, such as millieria or tinagma. as in previous studies, the position of lyonetia was unstable within yponomeutoidea. intriguingly, it was rarely associated with leucoptera, which has been considered a subfamily of lyonetiidae. this result is in line with that obtained by sohn et al.. in our result heliodinidae was grouped together with ypsolopha and ochsenheimeria, which deviates from all former suggestions, where it most often is considered close to bedelliidae. gracillarioidea had no distinct immature or adult characters common to all taxa in the character sampling of the present study. several small superfamilies, families and taxa, the position of which has greatly varied in analyses based on genetic data — millieriidae, tinagma, tanaoctena, epermenioidea, immoidea, carposinoidea, pterophoroidea, and choreutoidea — do not find a stable position in the combined analysis either. they tend to form one or two agglomerations, or join other groups with extremely low support values. a repeatedly recovered pattern is that alucitoidea, epermenioidea and carposinoidea form a monophylum. however, interesting observations were made regarding these taxa. characters listed below could prove useful in diagnosing families even though they may not have phylogenetic signal at higher levels. for example, the sculpture of the larval mentum is dentate in two groups, viz. choreutidae and alucitoidea s. l., respectively. the larvae of carposinidae have a characteristic bilobed extension between the pre- and submentum. in the analyses based on molecular data only, schreckensteinia is an unstable taxon, taking multiple alternative positions in the resulting phylograms. in the combined analysis, schreckensteinia consistently groups with the urodidae. they share several, yet to some extent homoplastic, synapomorphies in their larvae, notably the closely set stemmata, granulose sculpture of the mentum, medially narrowed prolegs with mesoseries on a3- and the closely set prolegs on a the adult urodus has pigmented patches on the cuticle of its mesothorax, which are also found in roeslerstammiidae. immidae have some obtectomeran features such as the setose outgrowth on the dorsal side of the pulvilli, and a possibly obtectomeran pupa. immidae have previously been placed in sesioidea and as an unplaced obtectomeran superfamily of their own. alucitoidea are monophyletic and are joined by isonomeutis amauropa meyrick, in the combined analysis. isonomeutis has previously been placed in copromorphidae. in our analysis, both with and without the third codon positions or fast evolving partitions removed, alucitoidea, carposinoidea and epermenioidea are grouped together. the family alucitidae is consistently nested in tineodidae. a tongue-shaped lobe laterad of the larval submentum is an autapomorphy for isonomeutis + alucitidae + tineodidae. on the basis of this evidence, we here synonymize tineodidae with alucitidae, syn. nov. the adult character common to pterophoroidea, here only including pterophoridae, is the acute angle of the metepimeron. however, the metepimeron of tanycnema anomala is also rather acute. pterophoroidea also have the prespiracular extension of sternum discussed above in the section on apoditrysia. galacticoidea are consistently recovered as the sister group to tortricoidea, which is also found in the results of the ml analyses of morphology. galacticidae also group with tortricidae in regier et al.. tortricoidea are monophyletic in all analyses. tortricoidea have been characterized by the flat papillae anales of the females. however, similarly flat papillae were occasionally also present in several other groups, e.g. zygaenoidea, pyraloidea, prodidactis, calliduloidea and noctuoidea. these taxa, however, have attributes that suggest the character has evolved repeatedly. the heliocosma group was removed from tortricidae by horak and common, and left without family assignment. in the study by mutanen et al. it was, together with another unassigned taxon, piestoceros, grouped with brachodidae. regier et al. recovered heliocosma together with tortricoidea. in our result of the analysis with three codon positions included, heliocosma, piestoceros and choristis form a monophylum that, with very weak support, forms the sister-group of the cossoidea + sesioidea + zygaenoidea assemblage. with the third codon positions removed, they are associated with immoidea. when the fast evolving partitions are removed, they group with millieria. this group is thus one of the most unstable ones in lepidoptera. no hypothesis of its position can be preferred over others at present and morphology does not offer clear clues about its affinities either. zygaenoidea, sesioidea and cossoidea are all, respectively, found to be monophyletic in our results and are grouped together in the same clade. in the smallest data set synechodes coniophora, however, is placed within cossidae. no other brachodids were included in the analyses of the and taxa data sets due to the unavailability of dna data. in the results based on dna data only zygaenoidea are monophyletic, but sesiidae and cossidae are intermixed. in other recent studies based on genetic data, sesiidae and cossidae have not been found to be monophyletic either. sesiidae and cossidae are also considered to belong to the same superfamily according to the most recent classification of lepidoptera. epipyropidae and/or cyclotornidae were not included in zygaenoidea in mutanen et al., cho et al., regier et al. or bazinet et al., but in the present study belong there with moderate support. the support seems largely due to synapomorphies of immatures shared with at least some other zygaenoidea: prolegs of larval a are approximate and arise on a common protuberance; the spiracle of abdominal segment of pupae is exposed, and there is a sculpted, laterally extended flange on the pupal eyepiece. species in both families are fully or partially parasitic at the larval stage. there are indications that molecular evolution of parasitic species may often be accelerated, possibly explaining the long branches and unstable phylogenetic affinities of epipyropidae and cyclotornidae in molecular data sets. this is also the first evidence that epipyropids and cyclotornids might be sister groups.fig. cyclotorna larval segment a ventral view, showing swelling on which both prolegs are positioned minet proposed the presence of two asymmetrical pits on the lower frontoclypeus of adults as a possible autapomorphy for cossidae. such pits were present in several, but not all cossids. at times, it was difficult to tell the difference between such pits and large scars of scales i.e. the round sockets in which scales were attached, on the frontoclypeus. better magnification would be necessary to observe the details. in the analyses based on molecular data only, sesioidea are not monophyletic, a result also found in mutanen et al.; regier et al., cho et al. and bazinet et al.. in the combined analysis, however, they are though with low bootstrap support. pennisetia, which in our dna-only analyses was not with sesiidae, joins other sesiids because it shares several adult characteristics with them, such as the sesiid wing-locking mechanism and a y-shaped suture/depression on the mesepimeron. adult characters common to brachodidae, sesiidae and castniidae are the elongated posterior tendons of the metafurca and the extension on the metafurca. several females in these families also have a telescopic ovipositor. the immature stages of sesioids and cossoids are generally uniform, suggesting no obvious support for the monophyly of sesioidea. the anteriorly more heavily pigmented ocular diaphragm was suggested as an autapomorphy of sesioidea. we were not able to examine the pigmentation of the ocular diaphragm, because of the varying success in treating the samples in koh solution so that the diaphragm would be sufficiently clearly visible. the size of the patagia, which has been mentioned as a synapomorphy of sesiidae, castniidae and brachodidae, is a continuous character with much variation among lepidoptera and therefore was not included in the present study. in addition, paysandisia shares the character considered autapomorphic for sesiinae, the bag-like protuberance concealing the pleural suture. in our result, cossoidea + zygaenoidea + sesioidea are in obtectomera. the support for this finding is, however, negligible and disappears when third codon positions are removed. several recent studies present evidence for the exclusion of these from obtectomera. our result is nevertheless intriguing as a similar pattern with the cossoidea + zygaenoidea + sesioidea assemblage nested within obtectomera has also been recovered by regier et al. in some of the analyses based on five protein-coding genes. if this pattern indeed proves the most likely, it emphasizes the level of convergence in morphology. in the morphology-only analyses, the superfamilies are not in obtectomera. choreutidae is a family without known phylogenetic position. choreutidae had previously been placed in a number of different families until minet moved them to a superfamily of their own. in our results, choreutoidea are monophyletic, and are placed in a very weakly supported clade together with immoidea. in the result obtained without third codon positions, the placement of choreutoidea shifts to a monophylum together with galacticidae and tortricoidea, a finding also recovered by the mp and ml of morphological data. the morphological data collected from the choreutids in the present study contain, however, intriguing features. the cervical sclerites of the adult prothorax are very similar to those of brachodidae as being strongly turned inwards in both groups. in brachodids, the lower arm of the sclerites is also shortened. similar, but not as strongly turned cervical sclerites were present in castniidae and some cossidae. other lepidoptera do not have such wide and inwardly turned cervical sclerites. brock and heppner have proposed a close affiliation of these families but based on characters not included in our study. the dentate sculpture of the larval mentum characterizes choreutidae, but also the isonomeutis-tineodidae-alucitidae clade, which does not seem closely related to choreutidae.fig. cervical sclerites. a anthophila fabriciana; brachodes appendiculata; cydia nigricana millieriidae are a small microlepidopteran family, which was formerly considered a subfamily in choreutidae. heppner considered millieriidae to be “somewhat of a missing link between other choreutidae and the remainder of the sesioidea”. based on evidence from morphological and molecular studies millieriinae were excluded from choreutidae and elevated to family level by rota. rota and kristensen compared the thoraco-abdominal articulation of millieriidae to those of choreutidae and found that unlike choreutids, millieriidae do not have lateral extensions on the apodemes on sternum in the study by mutanen et al. millieria was placed together with the unstable douglasiidae. in the present study, millieria is grouped, with weak support, in a clade that contains cadmogenes, schreckensteinioidea and urodoidea. in morphology-only analyses, millieria is associated either with choreutidae or roeslerstammiidae. the former position may be explained by the external resemblance of the adult of millieria to choreutidae in having a scaled base of the proboscis and scales between the eye and antenna. however, both features are also found in other groups. the immature stages of millieria give little information. the larva is an internal feeder, and devoid of obvious traits that could help group it with other taxa. its non-obtectomeran pupa is also generalized in appearance, thus not providing data for phylogenetic insights, either. the position of millieriidae with roeslerstammiidae is based on very homoplasious characters, such as the external invisibility of sutures between pupal abdominal segments 8– thus without convincing morphological support. the affinities of papilionoidea have varied among recent molecule-based studies. common to all of them is that the superfamily appears to be closer to the ‘microlepidopteran’ grade instead of being among the core macrolepidoptera as has been the conventional view. another pattern repeatedly recovered in studies based on sequence data and a study combining molecular and morphological data, has been that papilionidae is the sister of the other papilionoid families, including hedylidae and hesperiidae. in the present study, a surprising result is the emergence of a topology with the conventional concept of papilionoidea, i.e. hedylidae as the sister group of the rest of papilionoidea. the position of hesperiidae depends on the data set analyzed and on the inclusion or exclusion of third codon positions. if these are excluded, hedylidae are still the most basally arising group of papilionoids, but papilionidae and hesperiidae shift places. this result is in conflict with the recent studies based on genetic data in which both hedylidae and hesperiidae are nested within papilionoidea. the closest relatives of papilionoidea are not recovered with unequivocal support in recent molecular studies. the present result suggests that callidulidae is a sister-group to papilionoidea, and if third codon positions are removed, callidulidae are sister to papilionoidea + macroheterocera. a monophyletic gelechioidea is recovered in the combined and morphology-only ml analyses. the adult character supporting the grouping is the scaled haustellum. the character is not unique to gelechioidea, being also found in tischerioidea, pyraloidea, choreutoidea and millieriidae. a monophyletic gelechioidea has been found in all recent molecular studies on ditrysian phylogeny and targeted studies on gelechioidea. although gelechioidea are found to be monophyletic, support values for the clade are never very high when more than just a few exemplar species are included in the analyses. the large phylogenetic studies based on genetic data have also shown that gelechioidea belong to the obtectomeran clade, contradicting the earlier assumptions of their affinities with the non-apoditrysian yponomeutoidea. the closest relatives to this massive radiation of small moths are not known. the results of the present study do not give unambiguous clues to resolve this. several recent studies have recovered a sister relationship between gelechioidea and thyridoidea and in kawahara and breinholt they are sister to the clade callidulidae + thyrididae. our study suggests, though with negligible support, a sister group relationship between gelechioidea and a clade comprising alucitoidea, epermenioidea and carposinoidea, or if third codon positions are excluded, a sister group position to all other obtectomera. in our results, hyblaeoidea are the sister-group of pyraloidea. in the results based on the 268- and 422- taxon data sets, copromorpha lichenitis joins hyblaeoidea whereas in the 478-taxon set it is associated with % support with phycomorpha metachrysa, another copromorphoid, and together they weakly form a clade with carposinidae, currently also placed in copromorphoidea. the grouping of hyblaea and prodidactis is supported by a membranous projection present on the coxae of the male and a modification of the apex of the larval spinneret, as has been explained in kaila et al.. such a modified larval spinneret is also present in thyridoidea and several papilionoids, including all hesperiidae examined. coxal processes were also found in several species in geometridae, in cyclidiinae, and in anchinia daphnella. the affinity of prodidactis to several other families has previously been investigated by epstein and brown with inconclusive results, and the genus was assigned a family of its own. based on the two above-mentioned characters in common with hyblaea and congruent evidence from five genes, kaila et al. placed prodidactidae in hyblaeoidea. in the present study, the results of the analysis based on combined data and those based on dna data only, the clade formed by hyblaea and prodidactis tended to group basal to pyraloidea. morphological characters supporting this position are not clear. the papillae anales of the female of prodidactis are flat, as they are in several pyraloid species. in the maximum likelihood analysis of morphological data, hyblaea and prodidactis also group together, their position being close to mimallonoidea and thyridoidea. apart from the present study and mutanen et al., hyblaeidae and prodidactidae have not been present in phylogenetic studies based on genetic data. thyridoidea are recovered as monophyletic with strong support. thyrididae are a somewhat mobile group in our analyses, but often group as sister to callidulidae and papilionoidea. we found the adult character proposed by minet to support the monophyly of the superfamily, the lateral lobe on a very difficult to observe. out of our five exemplars of thyridoidea, the character was more or less clearly present only in thyris fenestrella. due to the difficulties in identifying this character, it was not included in our matrix. however, a number of thyridoid species showed on abdominal segment a a structure similar to the postspiracular bar found in papilionoidea. the presence of this bar could associate thyridoids in the same clade with papilionoidea, a result also observed in the molecular analysis by mutanen et al.. in most samples of papilionoidea, the postspiracular bars are strongly sclerotized and extend diagonally from the posterior edges of tergum towards the anterior edge of sternum the bar may be long and reach behind the spiracle, but sometimes short and harder to define. in addition, the bar is not always easy to distinguish from other sclerotized, but more perpendicularly directed, extensions of the posterior edge of tergum or the anterior edge of tergum reaching toward the sternum. several of the female thyridoidea have fine hair-like bristles by the hindwing acanthi in addition to the usual stouter bristles. in other recent studies the closest relatives of thyridoidea have varied, but they seem most often to group either close to papilionoidea, gelechioidea, calliduloidea, pterophoroidea or other small apoditrysian families. pyraloidea are monophyletic with strong support. the most obvious adult characters supporting the clade are the scaled haustellum and the characteristically-shaped abdominal hearing organs. differences in the hearing organs separate pyralidae and crambidae. the sclerotized structures of the tympanal organs have been regarded as homologous in pyralidae, crambidae, dudgeonidae, female uraniidae and geometridae because they are modified abdominal apodemes. a monophyletic pyraloidea, with monophyletic pyralidae and crambidae, have also been found in all studies based on genetic data in which representatives of both families have been present. mimallonoidea are placed as the sister-group of the macroheterocera. this position has also been obtained in several other studies based on genetic data. in regier et al., mimallonidae are in a clade with axia and doa and together form the sister-group of drepanidae. no particular character in the adult morphology explains this position as the sister-group of macroheterocera. in the morphology-only analyses mimallonoidea are monophyletic but are associated with thyridoidea and hyblaeoidea. in the results based on the 268- and 422-taxon sets, mimallo amilia, as sole representative of the superfamily in these analyses, is placed as sister to bombycoidea or lasiocampoidea. the macroheteroceran clade; 422: %; 473: % ) is one of the most strongly supported clades, but the support values for the relationships among the superfamilies are not strong. the obtained topology is as follows:  +  + noctuoidea)). the macroheteroceran superfamilies are recovered monophyletic in analyses based on or including genetic data lasiocampoidea are nested in bombycoidea), and the monophyly of most is supported by morphology. however, in morphology-only analyses the phylogenetic signal of these morphological traits does not suffice to recover the clades and they are intermixed. the proposed relationships among the superfamilies cannot either be unequivocally explained on morphological grounds. in recent molecular studies, the proposed relationships among the macroheteroceran superfamilies vary. in the result of the analysis of combined data geometroidea s. l. is obtained. sematuridae and epicopeiidae are sister groups and together form the sister group of uraniidae + geometridae. the obtained monophyly and within superfamily relationships are identical to those recovered in a very recent study based on molecular data on the phylogeny of geometroidea, which, however, are based on subsample of the present molecular dataset. rajaei et al. also characterized a new geometroid family, pseudobistonidae, based on pseudobiston pinratanai inoue, not included in the present taxon sampling. in their study, pseudobistonidae is placed as sister-group of epicopeiidae. rajaei et al. proposed five apomorphies supporting the monophyly of geometroidea. two of these are wing characters not included in the present study. the authors acknowledge that the other three characters are not clear-cut features. only one of these characters, the reduction of the fenestrae laterales, was included in our study. although the character is not always easy to define and may also be found in other ditrysia, our observations are generally in agreement with those by rajaei et al. a list of characters regarded as synapomorphies supporting the relationships among families within geometroidea is provided by rajaei et al.. several of these characters are continuous in nature and are not included in the present study. in the present study, the monophyly of geometroidea is, however, occasionally broken in the analyses of the 268- and 422-taxon sets and when the third codon positions or the two fastest evolving partitions are removed as sematuridae + epicopeiidae are placed in a position distant from the rest of geometroidea. epicopeiidae, traditionally placed in drepanoidea, also group with sematuridae in regier et al., and are sister to the uraniidae + geometridae clade. in cho et al. epicopeiidae and sematuridae are also sister-groups in the nt partitioned analyses, but not in the “all non-synonymous” analyses. bazinet et al. find epicopeiidae as sister to uraniidae + geometridae. in mutanen et al., epicopeiidae group with lasiocampidae, and sematuridae is found as the sister-group of uraniidae. together sematuridae and uraniidae form the sister-group of geometridae. minet suggested the pointed and elongate ventral arm of the tegula was an apomorphy for geometroidea and considered the presence of the sharp tegula in epicopeiidae and mimallonidae to be the result of parallel evolution. the ventral process of the tegula is indeed sharp in sematuridae, epicopeiidae, geometridae, but also in bombycoidea, lasiocampoidea, mimallonoidea and noctuoidea, thus in most macroheterocera. in drepanidae, the ventral arm of the tegula is blunt,. the tegula of uraniidae is elongate but spatulate. minet and scoble suggested the most reliable autapomorphy to unite geometroidea and drepanoidea sensu minet to be the strong anterior extension of pupal forelegs, and the presence of a transverse dorsal groove on pupal segment a yet with several secondary losses for both. the latter character was indeed found only among some representatives of these groups, while the former was highly homoplastic. apart from these families, it was observed also in bucculatricidae, attevidae and praydidae, urodoidea, schreckensteinioidea, pterophoroidea, thyridoidea, lasiocampoidea, bombycoidea, cimelioidea, most papilionoidea and many noctuoidea. a morphological character in favor of the relatedness of sematurids with uraniids could be the knob behind the antennae, situated medially of the chaetosemata, present in both sematura and urania. similar structures, although with more variation in shape, from flat pores to larger triangular protrusions, were found widely in taxa belonging to scopulini/sterrhinae and singular representatives of larentiinae and archiearinae in a study by sihvonen and kaila. such protrusions could thus be found more widely in taxa belonging to geometridae/geometroidea.fig. protrusions behind antennae. sematura lunus. character is also present in urania leilus and several geometrids drepanoidea are monophyletic and associate with axia, epicimelia and doa as has also been suggested by recent molecular studies. in several of the recent studies based on molecular data, drepanidae are the sister-group of the other macroheteroceran superfamilies. in the analysis with the fastest evolving partitions removed and of the 268-taxon set, the clade comprising drepanoidea, cimeliidae and doidae are sister to other macroheterocerans. however, when the partitions are not removed, the clade is nested in macroheterocera. minet and scoble proposed a modification of the larval mandibles, a large, flat area that is ventrally delimited by a carina, to unite drepanidae and epicopeiidae. the coding of this character was found to be continuous when other lepidoptera, especially macroheterocera, were examined in detail, and was thus impossible to code with any certainty, compromising its significance. it was excluded from the present analysis. the sister-group relationship of bombycoidea and lasiocampoidea is recovered in the present result as it has been in several studies based on genetic data. in mutanen et al. and the ml analysis based on dna-data only of the present study, the sister-group relationship between the two superfamilies was not recovered suggesting that it is the addition of morphological data that brings about this relationship in the present result. in the reduced taxon sets, with less taxa from both superfamilies, the sister relationship is not recovered either. instead, lasiocampoidea are nested in bombycoidea. in the ml tree of morphology only, lasiocampoidea, although not monophyletic, are in a clade with bombycoidea. the adult character supporting the affinity of these superfamilies, also found in some other groups, is the closeness of the prescutal clefts in several taxa. the taxon sampling for bombycoidea in the present study is relatively low compared to the sampling of other large superfamilies. this is because the inter-familial relationships and splitting into subfamilies within bombycoidea has been treated in more detail in recent phylogenetic studies. brahmaeidae swinhoe, and lemoniidae neumoegen & dyar, were synonymized by zwick based on evidence from dna and morphology. in the present result, brahmaea does not show a close relationship to lemonia and dactyloceras, probably due to the lack of dna and adult character data. the monophyly of noctuoidea is strongly supported and constituent families are also recovered as monophyletic in analyses based on combined and dna data. in the analyses based on morphological data alone, noctuoidea are not monophyletic although they share the well-known synapomorphy, i.e. the metathoracic tympanal organs. the other established synapomorphy, the presence of two md setae on larval metathorax, was excluded in the present analysis due to the impossibility to code it because many of ‘macrolepidopteran’ larvae are covered with microsetae that hide the possible presence of this character in this assemblage. the monophyly of noctuoidea has convincingly been recovered in all recent molecular studies, and its family systematics with both morphological and molecular insight, has been elucidated by zahiri. therefore, we do not discuss this superfamily further here. discussion our objective was to rigorously explore the morphology of ditrysian lepidoptera and to analyze the extensively coded data to infer a hypothesis for the phylogenetic relationships within ditrysia. the aim was to obtain supplementary insight from morphology to that derived from dna-based data, and compare this information to the recently published dna-based phylogenetic hypotheses. we consider this a way to supplement, evaluate and weigh these partially conflicting hypotheses. a feature in common to all these dna-based studies is the weak or non-existing support for groups among the most basal and the most “derived” taxa, i.e. for the so-called ‘backbone’ of the lepidopteran phylogeny. for this reason, we put special emphasis on acquiring as comprehensive data as possible on the groups that supposedly belong to this area of lepidopteran phylogeny. somewhat lesser attention was paid to inter-family relationships within superfamilies that have recently received specific attention. we combined morphological data obtained from larvae, pupae and adult specimens, coded from nearly the same sampling of exemplar taxa as in the study by mutanen et al.. the sampling covers nearly all recognized ditrysian families with a comprehensive representation of their subfamilies, and is supplemented with a number of taxa with unclear affinities. as has been the result in dna-based studies, even comprehensive morphological data were not able to bring forth strong support for any particular hypothesis on the interrelationships of the backbone ditrysia. possible reasons for this are the amount of homoplasy, general lack of informative characters except for apical nodes, conflicts between morphological and molecular data, and our decision to include and retain unstable taxa in the analyses. these matters are discussed below. lack of informative characters problems in inferring deeper level relationships in ditrysia with molecular methods have been explained by ancient rapid radiation of most major lepidopteran lineages. in ancient rapid radiation, lineage splitting occurs with close temporal spacing and very few character patterns are left as evidence of the relatedness of lineages that diverged from one another. the little phylogenetic signal that remains may be obscured by biases in sequence data, and/or by homoplastic changes on the apical branches. rapid radiation and difficulties in inferring relationships among lineages united by short nodes is a problem in phylogenetic analyses of many insect groups that have undergone spectacular diversification in a short period of time in the cretaceous-tertiary. recent divergence time estimates also suggest a significant increase in the diversification rate in ditrysia around this time period, million years ago. our examination of the morphology across the ditrysian superfamilies and families also confirms the lack of clear morphological characters for making inferences about these deeper evolutionary relationships. a possible explanation for this is that the first taxa of each monophylum supposedly had, and their extant progeny may still retain, predominantly plesiomorphic traits compared to the more derived members of the clade. paucity of apomorphies among the basal representatives thus leads a predictable pattern that evidence for inferring interrelationships of monophyla is scarce, if such ‘ancient’-looking taxa still exist and are included among the studied taxa. the collected morphological data were found useful in reconstructing evolutionary relationships at the apical nodes. our study affirms that most superfamilies are diagnosed by at least one fairly clear apomorphy, although diagnostic characters may not be found in all life stages. however, the morphology-only analyses were unable to recover some of the well-established superfamilies. an intriguing feature in our results is that even though no clear synapomorphies have been proposed to diagnose certain superfamilies nor were such characters found in the course of the present study or those suggested in earlier studies were refuted, such superfamilies were nevertheless also recovered in morphology-only analyses. an example is the superfamily gracillarioidea. although the exemplar taxa in the present study represent most of the ditrysian superfamilies, it has to be kept in mind that an unknown number of lineages have likely gone extinct in the millions of years that have followed the advent of ditrysia. it is possible that along with the extinction of these lineages informative characters have disappeared making the estimation of phylogenetic relationships among the ditrysian superfamilies even more difficult, just as in cases of insufficient taxon sampling where important phylogenetic information is missing. at present, it is not known how heavily, for example, the cretaceous-paleogene mass extinction event affected the diversity of lepidoptera. the observed morphological homogeneity of ditrysia at the deeper level can safely be considered to be due to real biological and historical causes, but limitations on character coding are also set by methods used or material available. some characters of potentially great influence may have been coded incorrectly for some taxa due to difficulties in making reliable observations. many such characters were probably overlooked, and others had to be excluded from the final analysis. examples include the presence or absence of the pupal frontoclypeal and epicranial sulci, which are easier to observe in exuviae than in pupae with uneclosed adults. many setal characters of larvae were excluded from the majority of the macrolepidoptera, as their final-instar larvae are usually covered by small setae that obscure the pattern of primary setae. this problem could have possibly been overcome if first-instar larvae could have been examined instead of final-instar larvae. this would have, however, dramatically reduced sampling due to the lack of available material. another option, using final instar larvae, supplemented with first instar larvae, would have been potentially misleading, as it would have been impossible to avoid duplication of some features in the data because of problems of homology. a problem in coding adult characters is the high plasticity of structures and thus the difficulty of dividing characters into discrete states. many characters are easily delimited when only a subset of the taxa is examined, i.e. only extreme states are observed, but the boundaries become exceedingly blurred as more taxa are examined and more variation and intermediate forms are encountered. many characters had to be discarded from the final analyses because of uncertain or artificial boundaries of states. typically, e.g., the degree of sclerotization in various body parts, and variation in strength of sutures and sulci, proved impossible to code unambiguously across the taxa. the use of morphometric tools could arguably increase the amount of character data as continuous characters such as length, size and details of shape could possibly be objectively divided into character states. uncertainty over the homology of certain structures across all groups was also a problem. certain structures may be related to the size of the insect, but these are not easy to detect and interpret. this was also a reason for our decision to omit wing venation and male genital characters entirely. both morphological and molecular characters contain information about the phylogeny of organisms but the importance of morphological data in the future of phylogenetics has been disputed see e.g.. the problem with morphological data is less with quality than with quantity. giribet discusses how morphology could be made “scalable” by utilizing new techniques that allow larger acquisition of morphological and anatomical data. such new techniques include immunostaining and confocal laser microscopy, and also 3-d imaging techniques that allow reconstructions of the internal structures of organism without having to destroy them. this enables the non-destructive study of rare museum specimens, ultrastructure of insects in amber fossils and in some cases live specimens. other new techniques to add to the list are magnetic resonance imaging, advanced digital microscopy, computer-assisted tomography and image recognition. the tools to share and work on morphological information have also undergone a revolution, e.g. morphbank:: biological imaging, morphobank and various discussion groups and web forums. homoplasy the frequency of homoplasy was found to be high in both larval and pupal characters. therefore, few characters of the immature stages, according to our analysis, are useful in diagnosing monophyletic groups. not even the non-ditrysian vs. ditrysian grades could be distinguished by data coded from immature stage morphology alone. also, the presumably strong pupal structures that according to previous hypotheses distinguish basal groups of ditrysia and obtectomera failed to recover these taxa unequivocally. one possible explanation is the repeated shifts to different modes of life in ditrysia, as has been shown to be the pattern at lower phylogenetic levels too. extreme examples of taxa with significant amounts of convergence are leaf-mining larvae. these larvae often have considerable convergent modifications in their morphology and reduction of structures. parsimony analyses of morphology may fail to group taxa correctly if the level of homoplasy is very high. examples of such taxa in our data are lyonetiidae, douglasiidae, ochsenheimeria and millieria. another case of obvious convergence is the tendency of diurnal lepidopterans in a number of clades to have clubbed antennae and papilionoidea. these examples demonstrate the potential effect of convergence on the overall topology, especially when there are several traits that by adaptation have led to appearance similar enough to make their recognition as distinct features hard. in many cases more thorough study of the morphology of allegedly homologous characters, e.g. “clubbed” antennae, is needed as new information could reveal them to be quite distinct despite their superficial similarity. a morphological trait in lepidoptera regarded to be phylogenetically very significant is the larval proleg shape. it has often been used as a distinguishing feature to separate ‘macrolepidoptera’ from ‘microlepidoptera’, but appears to be convergent in several lineages, viz. schreckensteinioidea – urodoidea, zygaenidae, papilionoidea and macroheterocera, but excluding pyraloidea. these groups usually have the following basic structure of the proleg: the base of the proleg is elongate and forms the greater part of the proleg; the proleg planta is asymmetric, laterally bulbous, and crochets are usually arranged as longitudinal mesoseries. the usual ‘microlepidopteran’ type of proleg consists of a small proleg base that is little more than a ring encircling the planta. the planta is either reduced or cylindrical, but not asymmetric, and the crochets most often are formed as a full or nearly full circle. even though different combinations of these features were observed in some taxa, the “macrolepidopteran” combination is similar in papilionoidea s. lato and macroheterocera. these groups have, as a rule, a similar larval mode of life as they are external leaf-feeders, usually gnawing on the edge of the leaf of host plants. this favors adaptations of the proleg shape to optimize the grip on the edge of the leaf. most pyraloid larvae live more or less concealed and do not need such specialized proleg modifications. the structure of larval prolegs may play a crucial role in the topology obtained in analyses with morphological data, as the shape is not one feature, but a combination of several traits. when a set of such non-independently evolved characters are coded from non-related taxa, they may have enough weight to pull the taxa together. the high level of homoplasy in the morphological traits is in itself interesting. just as dna can evolve fast enough to reach phylogenetic saturation, it has been shown that morphological characters can also show saturation. constraints exerted by intrinsic causes such as development, function or persistent selective trends do not allow for the continuous rise of new characters. clades seem to exhaust their available character space and this leads to homoplasy. however, convergence does not only concern the morphology of organisms. genes can also be under selective pressure for a particular function and converge at the nucleotide level. although we aimed to avoid size-related characters in our morphological matrix, and tested the effect of removing characters we suspected could be correlated with size, characters or character combinations less obviously correlated with size may have remained. as the signal for the deep nodes is weak, it is possible that such characters led to artificial clustering of groups by species’ overall size. combining and comparing morphological and molecular evidence even though morphological characters comprised only to % of the total amount of data, depending on the analysis, the incorporation of morphological data had, in places, significant effects on the results. adding morphological data established more stable positions for some taxa, such as cyclotornidae and epipyropidae, which in molecular analyses have defied well-established positioning within their putative superfamily. likewise, the monophyly of sesioidea and cossoidea were now recovered, even though with entirely different morphological traits to those suggested earlier. tineoidea excluding metapherna were recovered as monophyletic. in dna-based studies metapherna, an alleged tineid not included in previous analyses, was never found to associate with tineoidea. rather, it was recovered to form its own clade linking tineoidea and the more advanced lepidoptera – a result already hinted at by robinson and nielsen. it is not uncommon that morphology and genes tell different stories about the evolutionary history of organisms. this is caused by the differences in the rates of change and the evolutionary mechanisms acting on these rates. conflict between datasets and resulting phylogenies may be caused by convergent evolution, methodological reasons, or uncertain homology assessment. cases of conflict between datasets are intriguing as they may allow the identification of the underlying biological processes causing the conflict. when conflicting clades arise in analyses, it is important to consider the effect of the methods chosen to analyze the data, but also to seek alternative evidence and information about the ecology of the organisms. combined with the molecular dataset consisting of eight gene regions, the morphological data were able to distort not only weakly supported relationships, but also fairly strongly supported ones, such as the interrelationships in papilionoidea. without alternative sources of evidence, conflicting results would be difficult to evaluate. recent studies based on different sets of genes to those in the present study provide some insight into the matter. the taxon sampling of those studies has been, however, often lower and not permitting simple comparison. also, the high support values for proposed relationships in low taxon sampling studies based on transcriptomic, genomic or proteomic data may be misleading. in the absence of more closely related taxa, more distantly related taxa may be drawn together but still obtain relatively high support values. this phenomenon is illuminated by the conflicting hypotheses with strong support values in large molecular data sets of. on the other hand, many patterns with low support observed in the study of mutanen et al. with dense taxon sampling, but lower amount of molecular data, have subsequently been obtained in other studies with much larger gene sampling. conflicting patterns of relationship between molecular and morphological data, each with good support, is found, for example, in the position of douglasiidae. the morphological result always places this family in gracillarioidea, as in the superfamily concept of davis and robinson, also based on morphology. however, molecular evidence is strongly in conflict with this result, unequivocally placing them in apoditrysia, which is also the case for our combined data analysis. unstable taxa in this study, we attempted to cover the diversity of ditrysian lepidoptera as comprehensively as possible. therefore, we paid special attention to the inclusion of taxa with unknown phylogenetic positions, or which have been unstable in recent studies based on dna-data. our aim was to examine whether the combination of extensive morphological data could provide information to infer their likely relationships. we presumed that taxa appearing unstable in genomic data are not necessarily the same as those appearing as such in morphological data. we further expected that in analyses with both types of data combined at least some of the rogues would become more stable. combined data indeed gave moderate support at least in some analyses for the inclusion of epipyropidae and cyclotornidae in zygaenoidea, a result that in molecular studies had either not been obtained, or only with a weak indication for epipyropidae. the combination of data also brought out the sister-group relationship between bombycoidea and lasiocampoidea, not previously recovered with dna data similar to those used in this study. although bootstrap support is not strong, the evidence supporting the monophyly of schreckensteinioidea + urodoidea seems a plausible hypothesis when morphological data are also considered. similarly, when compared with results of studies based on molecular data only, the association of galacticoidea with tortricoidea also gains moderate to strong support in our analyses for the first time. a monophylum comprising alucitoidea, epermenioidea and carposinoidea is repeatedly recovered in the analyses of combined data, although with low support values. adding morphological information also increased support for the placement of prodidactidae in hyblaeoidea. on the other hand, there are cases where sequence data contradicts morphological information. for instance, genetic data strongly suggests that douglasiidae are not gracillarioids, but apoditrysians. nevertheless, combining data did not stabilize the position of many taxa that have been unstable in previous studies. the unassigned heliocosma group, choreutidae and immoidea are typical examples of such taxa. as the “backbone” of ditrysian lepidoptera did not gain much further support from combining morphology and molecular data, one can argue that several of the larger assemblages behave in an unstable manner, with presumably close groups displaying alternative arrangements among each other, e.g. calliduloidea, papilionoidea and thyridoidea. the cossoidea + sesioidea + zygaenoidea assemblage is particularly intriguing as it displays different positions and compositions among data sets. in the results of morphology-only analyses, zygaenoidea were not unambiguously associated with cossoidea + sesioidea, which, in turn, often form a sister-group relationship with tineoid group families. in the combined data with all third codon positions included, this three-superfamily assemblage was placed within obtectomera. this relationship was broken when the third codon position nucleotides were removed. an intriguing pattern in our data is that several of the unstable taxa are internal feeders. these include cemiostominae, lyonetiinae, bedelliidae, heliodinidae, millieriidae and douglasiidae. similar trends have been found in gelechioidea, where a wide array of life habits is found. although not tested using quantitative methods, these taxa are often observed to display disproportionally long branches and more often display multiple positions in phylograms. this raises the question whether shifting to leaf-mining or other demanding life modes requires not only considerable morphological adaptation, but also rapid genetic evolution in many loci. this would support the notion that taxa with long branches, i.e. those with a disproportionally fast evolutionary rate, would have a higher likelihood to behave as rogues. in such cases, adding more genetic data might not be able to solve the problem. another example of taxa with an unusual mode of life are epipyropidae and cyclotornidae, whose larvae have a two-phased development, and are ectoparasites of other insects; a life-history trait that is possibly associated with rapid molecular evolution. as discussed above, stable placement of these families has not been reached with dna data.  our results support previous assertions of the paucity of clear morphological characters diagnosing larger assemblages of ditrysian superfamilies. previously proposed characters are tenable to some extent but alone do not allow the unequivocal phylogenetic placement of taxa. the degree of homoplasy is high and likely hid phylogenetic signal in the analyses. nevertheless, most currently recognized superfamilies were recovered, and were found to have at least one previously recorded diagnostic character. in some cases, the suggested apomorphies appeared to be supported only weakly, or remained entirely unsupported. examination of the distribution of over morphological traits across ditrysia revealed characters bearing evidence on the evolutionary relatedness among smaller groups, such as families and in some cases superfamilies. these characters also have potential use in identification keys and in identification of fossil samples. the inclusion of morphological data also allowed to “pin down” several taxa that for some reason appear unstable in analyses of molecular data. in the near future, we will begin to see the results of the first large-scale studies based on genome-wide data aimed at solving the deep divergences in the evolutionary tree of ditrysia. the morphological data collected in the course of this study will certainly prove useful in evaluating and validating the results. likewise, the results of future phylogenomic studies will help explain the distribution of some of the morphological characters and identify possible misinterpretations. comparing structures across the ditrysian subfamilies will also gain us insight into the evolution and diversity of morphological characters. the distribution of these characters and their reversals across the different clades can be used in future studies on the evolution of these traits. availability of supporting data the data sets supporting the results of this article are available as additional files and sequence data are available in genbank. morphological character data are available in morphobank, project http://www.morphobank.org/. additional files additional file 1: list of exemplar taxa in study. additional file 2: sequence lengths. additional file 3: gene summary. additional file 4: genbank accession numbers. additional file 5: character list. additional file 6: character matrix. additional file 7: phylogenetic tree from maximum likelihood analysis of combined data with third codon position removed from molecular data except from ef1a. taxa, characters. additional file 8: phylogenetic tree from maximum likelihood analysis of dna data with third codon position retained, taxa, bp. additional file 9: phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data, taxa. molecular data were categorized according to the rate of evolution with the program tiger and partitioned with ratepartitions. two fastest evolving partitions removed. additional file 10: one of the equally most parsimonious trees from the parsimony analysis of morphological data, taxa, characters. additional file lists morphological data for taxa, but in the present analysis, data for two species of thereutis were concatenated into one terminal taxon as were those of two species of cyclotorna making the total number of taxa in the analysis additional file 11: phylogenetic tree from maximum likelihood analysis of morphology, taxa, characters. additional file lists morphological data for taxa, but in the present analysis, data for two species of thereutis were concatenated into one terminal taxon as were those of two species of cyclotorna making the total number of taxa additional file 12: phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data, taxa. the data set includes all the taxa for which both morphological and molecular data were available. additional file 13: phylogenetic tree from maximum likelihood analysis of combined morphological and molecular data, taxa. the data set includes all the taxa for which dna data were available. morphological data were available for of these taxa. competing interests the authors declare that they have no competing interests. authors’ contributions lk and ps had the original idea for the study. the study was designed by all authors. mh, lk and ps carried out the preparation and morphological examination of the samples. molecular data were prepared by mm and nw. mh, nw and lk carried out the phylogenetic analyses. all authors contributed to writing the manuscript. all authors read and approved the final manuscript. we would like to express our deepest gratitude to the following people and collections for loan of material: leif aarvik, john w. brown & m. alma solis, marc epstein, nico elfferich, marianne horak & ted edwards, ole karsholt, jari kaitila, joël minet, kari nupponen, teemu klemetti, leptree project and people, texas a&m university and many other members of the international community of lepidopterists. laura törmälä is thanked for laboratory work. we are also thankful to jurate de prins, joaquin baixeras, and toomas tammaru for their comments on an earlier version of this manuscript. library secretary of the zoological unit marjatta mikkonen was of great help in the search for literature. marianna teräväinen corrected the language. we would also like to thank charles mitter and an anonymous reviewer for comments that greatly improved the manuscript. this research was funded by the finnish cultural foundation, oskar Öflund foundation, ella and georg ehrnrooth foundation, and the university of helsinki; academy of finland, grant #.
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. the data are analyzed alone and combined with sequence data. results include the recovery of a monophyletic tineoidea, sesioidea and cossoidea, and a stable position for some unstable taxa. we present the most comprehensive morphological dataset on ditrysia to date, consisting of over morphological characters. morphological characters are found phylogenetically informative in resolving apical nodes, but characters serving as evidence of relatedness of larger assemblages are few. several such taxa, however, remain unstable even though morphological characters indicate a position in the tree. 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. analyses are performed using maximum likelihood methods, and parsimony methods for the morphological dataset. 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. 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 taxa representing nearly all ditrysian superfamilies and families. we explore whether combining morphological data and dna-data can stabilize taxa that are unstable in phylogenetic studies based on genetic data only. the full dataset consists of exemplar species. ditrysiamorphologylepidopteraphylogenyunstable taxatotal evidencedeep divergenceshttp://dx.doi.org/ /501100003125suomen kulttuurirahastohttp://dx.doi.org/ /501100005763oskar Öflunds stiftelsehttp://dx.doi.org/ /501100003502ella ja georg ehrnroothin säätiöhttp://dx.doi.org/ /100007797helsingin yliopistohttp://dx.doi.org/ /501100002341suomen akatemia1110906kaila lauri issue-copyright-statement© the author 2015
 ditrysia comprise close to % 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 taxa representing nearly all ditrysian superfamilies and families. we present the most comprehensive morphological dataset on ditrysia to date, consisting of over morphological characters. the data are analyzed alone and combined with sequence data. the full dataset consists of 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.  morphological characters are found phylogenetically informative in resolving apical nodes, 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. several such taxa, however, remain unstable even though morphological characters indicate a position in the tree. 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.  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. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords ditrysiamorphologylepidopteraphylogenyunstable taxatotal evidencedeep divergenceshttp://dx.doi.org/ /501100003125suomen kulttuurirahastohttp://dx.doi.org/ /501100005763oskar Öflunds stiftelsehttp://dx.doi.org/ /501100003502ella ja georg ehrnroothin säätiöhttp://dx.doi.org/ /100007797helsingin yliopistohttp://dx.doi.org/ /501100002341suomen akatemia1110906kaila lauri issue-copyright-statement© the author 2015 