Progress, pitfalls and parallel universes: a history of insect phylogenetics
Abstract
The phylogeny of insects has been both extensively studied and vigorously debated for over a century. A relatively accurate deep phylogeny had been produced by 1904. It was not substantially improved in topology until recently when phylogenomics settled many long-standing controversies. Intervening advances came instead through methodological improvement. Early molecular phylogenetic studies (1985–2005), dominated by a few genes, provided datasets that were too small to resolve controversial phylogenetic problems. Adding to the lack of consensus, this period was characterized by a polarization of philosophies, with individuals belonging to either parsimony or maximum-likelihood camps; each largely ignoring the insights of the other. The result was an unfortunate detour in which the few perceived phylogenetic revolutions published by both sides of the philosophical divide were probably erroneous. The size of datasets has been growing exponentially since the mid-1980s accompanied by a wave of confidence that all relationships will soon be known. However, large datasets create new challenges, and a large number of genes does not guarantee reliable results. If history is a guide, then the quality of conclusions will be determined by an improved understanding of both molecular and morphological evolution, and not simply the number of genes analysed.
1. Introduction
We like to think of scientific research as insulated from human bias and personality. Like other fields of science, phylogenetics follows trends as ideas are rejected or accepted, influenced by new information. However, collective consensus comes not just from a series of technological advances and discoveries, but also from human interactions. New ideas are often rejected for years, even if they are supported by strong evidence. These are exciting times for evolutionary biologists as new technologies give us hope that the resolution of the tree of life is within sight. However, times have been exciting for decades and this optimistic sentiment has arisen with every new technology. It was only 25 years ago that phylogenetic trees (box 1) generated with a few hundred nucleotides were considered revolutionary, just as the application of cladistic (box 1) principles with a defined methodology was revolutionary a decade before that. With the large datasets we have today, some previously intractable questions now appear solved. The authors of this work have witnessed many of these changes, and we present our insights on this history, recognizing that others may remember things differently. We focus this review on the relationships among insect orders, missing many fine works on arthropod phylogeny, and intra-ordinal studies. We attempt to maintain a rough chronological order, considering three main periods: morphological phylogenetics, when morphology was the only source of data (roughly before 1990); the Sanger (box 2) sequencing period, where a few genes dominated most studies (roughly 1990–2005); and the current state of the art with datasets so large that traditional ways of analysing them are no longer feasible. New challenges will doubtless arise in the age of big data, but at least we can look back at previous trends with hindsight in order to learn from history.
Phylogenetic terms: cladistics.
Phylogenetic trees. Graphical representations of evolutionary relationships. Synonyms: evolutionary trees, phylogenies, genealogies.
Monophyletic group. A group of organisms (taxon) that is defined by a most recent common ancestor, and all of its descendants. Also known as a clade.
Cladistics. An approach in systematics that bases all classification on ‘clades’ (i.e. monophyletic groups). Cladistics was developed by Hennig and insists that all named groups (taxa) be monophyletic, as evidenced by shared derived characters (‘synapomorphies’). After Hennig's death, a group of cladists started using the term to refer to a set of numerical analytical procedures which aim to reconstruct a phylogeny based on character state matrices and parsimony. ‘Cladistics’ can mean two distinctly very different things: Hennig's, focusing on monophyly and synapomorphy, or the cladists', based on parsimony methods.
Parsimony. A broad scientific principle that prefers simple over complex explanations. In a phylogenetic context, parsimony refers to preferring a tree with the fewest possible character state transformations. Thus, whenever possible, transformations are assumed to be shared among taxa and thus placed on internodes as synapomorphies, rather than as homoplasies.
Sister group. The most closely related taxon to a group of interest.
Ingroup. A taxon under investigation. For this review, the ingroup is Hexapoda.
Outgroup. A taxon outside the group under study. For this review, the outgroups could be any non-hexapod, but the best would be other arthropods.
Character polarization. Determination of the evolutionary direction of a character, which determines whether a character state is ancestral (plesiomorphic) or newly derived (apomorphic).
Phylogenetic terms: analytical.
Synapomorphy. A shared, derived character (feature) that can be used as an argument for a group being monophyletic (box 1).
Homoplasy. Character state evolving more than once on a tree or changing back to its original state (redundant evolution). Parsimony (box 1) attempts to minimize homoplasy. Homoplasy creates phylogenetic noise (misleading signals).
Distance analysis. Methods that reduce all character differences between pairs of taxa to a single value, their pair-wise distance. Trees are then constructed by grouping the most similar taxa. Distance methods are criticized by cladists as being phenetic.
Phenetics. Organisms are grouped or classified based on overall similarity in their phenotype or appearance, rather than on derived character states only.
Likelihood analysis. A statistical method of selecting among possible trees based on the probability of the data under a model of evolution.
Long branch attraction. A phenomenon that misleads phylogenetic reconstruction. On long branches, shared phylogenetic noise (homoplasy) accumulates and overrides the true phylogenetic signal on short internal branches of a phylogenetic tree.
Bootstraps. A subsampling of phylogenetic data which creates a number of pseudoreplicate datasets. These pseudo-replicates are then analysed individually, and their results are summarized on a consensus tree in order to estimate conflicting signal and provide an assessment of support for individual clades. (Jack-knifing is similar, but the new pseudo-replicate datasets are generated by random deletions of columns of characters.)
Branch support. Quantitative measures to assess confidence for particular clades in a phylogeny. Examples include bootstraps, jack-knifing, posterior probabilities and Bremer support. Congruence among independent datasets could also be considered as branch support, but is seldom quantitatively expressed.
Root. A hypothetical taxon assigned as the most recent common ancestor of all the taxa in a phylogeny. A root is used to assess the polarity of a phylogeny. Outgroups (box 1) can be used to help estimate the position of the root.
Node. The point at which an ancestral lineage splits into two lineages in a phylogeny.
Internode. The lines in a branching diagram between nodes (internal branches on a tree). In a phylogeny, an internode represents an ancestral lineage. Synapomorphies occur on internodes. The longer an internode exists, the more chance for synapomorphies (either molecular or morphological) to accumulate. Short internodes are generally the source of controversy, because they have a lower probability of accumulating informative substitutions.
Substitution. An observed change in a character. For molecular data, substitutions are related to mutations, but because lethal mutations are seldom observable, substitutions are mutations that have survived the filter of selection.
Sanger sequencing. The dominant method of DNA sequencing during the 1980s to 2005.
Restriction sites. Short unique motifs scattered throughout the genome which can be cut by certain restriction enzymes, yielding fragments that can be visualized on a gel providing snippets of the DNA sequence information.
The numerous names of orders and other higher-level taxa for a group as diverse as insects pose a significant challenge to the non-entomologist reader. Common names like ‘angel insect’ or ‘gladiators’ are often as obscure as the scientific ones in Latin or Greek. For this review, we direct the reader to the figures for the common names of the orders and to appendix A for a translation of super-ordinal names. We focus especially on four controversial deep-branching taxa: Entognatha, Palaeoptera, Polyneoptera, and Holometabola. The controversy arises from persistent conflicting evidence that suggests contradictory groups. The entognathous hexapods with internalized mouthparts include mostly tiny, wingless, litter-dwelling species that appear very early in the fossil record. The palaeopteran insects comprise mayflies, dragonflies and damselflies, characterized by wings that cannot be folded. Polyneoptera is the name given to a diverse group of insects, such as grasshoppers and close relatives, walking sticks, roaches, mantids, earwigs, stoneflies and some other groups, usually but not always characterized by leathery forewings. The holometabolous orders exhibit complete metamorphosis where the larva undergoes an amazing reorganization of the body during the pupal stage before it changes into the winged adult form. A confusing convention for the non-entomologist is the inconsistent use of the names Hexapoda and Insecta. Hexapoda (insects in the widest sense) comprise all six-legged arthropods, including the three entognathous orders, whereas Insecta excludes the entognathous orders (appendix A). Phylogenetic terms that many readers might not be familiar with are defined in numbered boxes, which are referenced at the first usages of the term.
2. Pre-Hennigian concepts in insect taxonomy and phylogeny
The roots of insect systematics go back to the sixteenth, seventeenth and eighteenth centuries. Important pioneers of entomology were the Italian naturalist Ulisse Aldrovandi (1522–1605), the Dutch doctor and microscopist Jan Swammerdam (1637–1680) and the German naturalist August Johann Rösel von Rosenhof (1705–1759) [1,2]. In the middle of the eighteenth century, the Swedish botanist Carolus Linnaeus (1707–1778) described more than 10 000 species in his Systema Naturae [3], including over 2000 insects. His ordinal names refer to the characteristics of the wings, e.g. Heteroptera (heterogeneous forewing), Hymenoptera (membranous wings) and Coleoptera (sheath-like forewing). Although his views evolved, Linnaeus was an essentialist in his early works, embracing the—at that time—commonly held belief that organisms were given an ‘essence’ by the Creator, which could be slightly modified but never fundamentally changed. Linnaeus' system remains useful to this day because it was based on characters that, unknown to him, are heritable and hierarchically organized through evolution. The Danish entomologist Johann Christian Fabricius (1745–1808) described 9776 insect species. Unlike his mentor Linnaeus, he emphasized the importance of mouthparts and the potential usefulness of genitalia [4]. Another prominent entomologist of the era was Pierre André Latreille (1762–1833). In his major work [5] he outlined insect families for the first time and used a broad spectrum of characters [2]. Together with explicit criteria for homology, this was an important step towards an evolutionary concept of classification.
The evolutionary theory developed by Charles Darwin and Alfred Russell Wallace [6,7] laid a new foundation for classifying organisms, but had limited immediate impact on insect systematics [1]. Ernst Haeckel (1834–1919), an energetic promoter of Darwin's ideas in Germany, dealt with insects among many other groups. His classification included five ‘legions’ based on how insects feed [8]; we see today that it only partly reflected phylogenetic relationships. However, Haeckel presented the first explicit phylogenetic tree of insects [8, p. 710]. In 1904 [9], a remarkable study covering the entire Hexapoda was published by Carl Börner (1880–1953). Börner was a specialist on grape phylloxera (Daktulosphaira vitifoliae), an almost microscopic aphid-like insect that is a major pest of grapes. He was also a collector of springtails (Collembola), small hexapods that are common in leaf litter. As a young scientific assistant, he discussed cephalic structures in great detail. He focused on the hypopharynx, a central element of insect mouthparts and one of the most difficult character systems to explore. Even though his approach lacked a repeatable methodology, his phylogenetic tree (figure 1) comes close to concepts developed decades later. Naturally, since our current cladistic (box 1) concept of reserving names for monophyletic (box 1) groups [10,11] was not developed until the 1950s and 1960s, Börner's classification is partly inconsistent with the branching pattern shown in the tree. For instance, he placed the phenotypically similar Archaeognatha and Zygentoma in the Order Thysanura. (Figure 1; see appendix A here, and throughout, for a definition of taxon names.)
A highly productive North American entomologist of the early-twentiethth century was G.C. Crampton [12,13], whose phylogenetic tree from 1938 [14] was another hypothesis that came remarkably close to modern concepts (see fig. 1 in Engel & Kristensen [2]). Important works were published by Imms [15], Snodgrass [16], Weber [17,18] and also by Handlirsch, who was frequently cited in Hennig's later work [11] (see [19]). Handlirsch attempted a classification reflecting phylogeny, but believed that a purely phylogenetic system was not possible [11]. Even in studies published posthumously in 1937 [20] and 1939 [21] Handlirsch considered the extinct winged Palaeodictyoptera as the ancestors not only of Pterygota but of all other insects including the wingless (apterygote) orders [11].
3. Hennig's breakthrough
Willi Hennig (1913–1976) revolutionized systematics and classification [22] in the last century with his theoretical work, offering clear and repeatable methodology. Works published prior to Hennig are often referred to as ‘intuitive’. This is perhaps unfairly pejorative when you consider that their remarkably accurate phylogenetic insights were often based on expertise gained through meticulous observation, rather than intuitive hunches. However, before Hennig's methods were widely adopted, systematists would postulate relationships based on shared characters that they deemed particularly important. In this respect, phylogenies could be considered as imparted wisdom, rather than science. Hennig's method involved distinguishing ancestral (plesiomorphic) and derived (apomorphic) features. He also developed a more precise concept of monophyly (box 1), under which no descendants of the most recent common ancestor could be excluded from a named group (clade). Hennig reconstructed phylogenies with an iterative, stepwise approach. Using putative shared-derived character states (synapomorphies), he successively established sistergroup (box 1) relationships. Distinguishing ancestral from derived character states required the definition of a taxon outside the group of interest for comparison (an outgroup (box 1)). The outgroup method was introduced as a formal procedure in the early 1980s [23,24] even though it had already implicitly been used by Hennig [22]. Hennig's phylogeny [11], published in 1969, is widely considered to be the starting point of modern insect phylogenetics (figure 2). Despite his precise methodology, his hypotheses were quite similar to earlier trees. They changed with time, as can be seen by comparing the phylogenetic concept presented in Hennig's 1969 phylogeny (figure 2) [11] with his earlier work [10].
Hennig's ‘Phylogenetische Systematik’ [22], was not a completely new concept when it was published in 1950. The botanist Zimmermann developed similar ideas in the 1930s, and Sturtevant used a very similar approach in his taxonomic studies of fruit flies (Drosophilidae; [25]). Moreover, it is apparent that ideas similar to Hennig's were implicitly used before his methods were formalized. It is impossible to consider Börner's phylogeny without recognizing an approach that went beyond intuition. Aside from the primacy of synapomorphies, a major point of Hennig's concept was that classification should be strictly linked to phylogeny. The requirement that named taxa be monophyletic originated with Hennig, but the unique value of synapomorphies was loosely recognized by systematists earlier. Herbert Ross (1908–1978) [26], for instance, was polarizing (box 1) characters relative to a hypothetical ancestor in 1937, and he indicated derived states with marks on the internodes (box 2) of his insect phylogeny in 1955 [27]. The phylogeny in his 1965 textbook [28] is almost as close to current concepts as morphology has ever been. However, as advocated in general by Ernst Mayr [29] (see also Nelson's reply [30]), Ross gave names to paraphyletic groups (groups that do not include all descendants of the deepest ancestor). If your concept of ‘dinosaur’ does not include birds, then you accept paraphyletic taxa too. Systematists today consider, for example, birds to be a subgroup of Sauropsida, a clade that also contains dinosaurs and extant reptiles such as turtles, lizards and crocodiles. Mayr and followers understood that birds had been derived from a paraphyletic assemblage of reptiles, but still found ‘reptilia’ to be a useful term representing a different evolutionary level, just as we sometimes use ‘apterygotes’ as a name for the ancestrally wingless hexapods, even though we understand that they are not a monophyletic group. Generally, when systematists put a name in quotes, it is to indicate that they understand it to be a paraphyletic group, and are waiting for the term to fade into disuse.
A remarkable study was published by the Argentinian entomologist Alvaro Wille [31] in 1960. Although he distinguished ‘primitive’ from ‘specialized’ or ‘unusual’ features, he also characterized groups by a mixture of plesiomorphies and apomorphies. The major clades on his tree, however, were characterized by evolutionary innovations (figure 3). Another important work of the time was Hinton's 1958 review [32]. Hinton made some bold statements that appear untenable today, such as ‘the polyphyletic nature of the old groups Myriapoda and Hexapoda’, but his evaluation of morphological characters, taken largely from the head, including a detailed scrutiny of larval muscles, helped elucidate the evolution of Holometabola.
Gerhard Mickoleit, who graduated under the insect morphologist Hermann Weber at the University of Tübingen and attended seminars given by Hennig in the early 1970s, investigated several groups of insects, with a focus on genital structures, especially the ovipositor. This included thrips [33], Neuropterida (lacewings and close relatives), beetles [34], and fleas, flies and scorpion flies [35–38]. In 1973 [34], he provided specific evidence for the first time for a close relationship between neuropteroids and beetles.
4. Geographical isolation and ‘parallel universes’
The importance of the contributions made by Russian entomologists and insect palaeontologists is reflected by numerous citations by Hennig [11]. Formal names for important higher ranking taxa such as Palaeoptera, Neoptera and Polyneoptera were introduced by Russian scientists [39,40]. Moreover, Russian palaeontologists, notably A. V. Martynov, B. B. Rohdendorf, V. V. Zherikin and A. G. Ponomarenko [41,42] (reviewed in 2002 [43] and 2009 [1]) made immense contributions to the knowledge of fossil insects that provided a critical window for observing past morphology.
Through the entire twentieth century, most Russian entomologists maintained a conservative approach, with traditional descriptions based on morphology, without formal cladistic phylogenetic character evaluations. International collaboration was partly impeded by linguistic barriers, but also by the isolation of the Soviet Union and the pseudo-scientific Lysenkoism, an antigenetic view that was politically favoured [44]. The limited exchange and cooperation is still reflected by the strikingly different nomenclature for high-ranking taxa, such as Scarabaeona for Pterygota, Scarabaeiformes for Holometabola and Scarabaeoidea for Coleoptera [4,43].
A prominent and highly efficient Russian entomologist of the nineteenth century was Victor I. Motschulsky, who published numerous works on biogeographic, faunistic, or systematic aspects of entomology, most of them on beetles [45]. Georgij G. Jacobson had a crucial impact on the development of Russian systematics in the early twentieth century. He is best known as the author of the 1905 magisterial ‘Beetles of Russia, Western Europe and neighbouring countries’ [46], including an impressive catalogue with keys for the identification of all known Eurasian genera. More recently, the palaeoentomologist Alexandr P. Rasnitsyn described approximately 250 new genera and over 800 new species of fossil insects. He suggested a sistergroup relationship between Hymenoptera (sawflies, bees, wasps and ants) and the remaining Holometabola [43] before this was established with formal analyses of extensive morphological or molecular datasets [47–50] (see Ronquist's and others' reanalysis [51]). Rasnitsyn [43] suggested that insect flight originated from gliding [52]. Phylistics, his alternative approach to cladistics, as discussed by Brothers [53] explicitly accepts paraphyletic groups (e.g. †Caloneurida [43])
Before the Internet, geography and language also played a role in isolating phylogenetic communities from Europe, America and East Asia. For example, a profound treatment of insect morphology was presented by René Jeannel in 1949 [54], but has rarely been used outside of the French community. Hennig's work was unknown to most Americans until it was translated into English in 1966 [55]. Although systematists were aware of work in other countries, the Meetings on Insect Phylogeny in Dresden played a major role in fostering collaborations between workers from different parts of the world, although they have not seen significant Russian or Latin American participation. These meetings were organized for the first time in 2003 by Klaus-Dieter Klass and Niels Peder Kristensen [56]. Most members of the ‘1000 insect transcriptome evolution’ (www.1KITE.org) initiative first became acquainted at these meetings. The 1KITE team created our current best estimate of insect ordinal phylogeny (figure 5) with the largest dataset assembled to date. Europe was an ideal meeting place, because the particular brand of cladistic fervour in America that was characterized by name-calling and personal insults was less pronounced there. The Europeans absorbed new ideas quickly, and went about their business in developing new centres of insect phylogenetics based on emerging techniques in morphology, and the refinement of model-based molecular phylogenetics (reviewed in [57]).
5. Post-Hennigian approaches
The classical tradition of insect morphology and phylogeny was upheld on a high level by Niels P. Kristensen (1943–2014) of the Zoologisk Museum in Copenhagen. He published outstanding morphological treatments of lepidopteran key taxa [58–61], profound reviews of insect phylogeny [62–64] and landmark volumes on systematics and morphology of Lepidoptera in the Handbook of Zoology series [65–68]. Even though he never performed computer-assisted analyses, his critical contributions helped refine character interpretations and pointed out problematic phylogenetic issues. A characteristic feature of Kristensen's approach was a deep-rooted scepticism, reflected by largely or completely unresolved parts of his phylogenetic trees. His display of polyneopteran relationships [63] became known as ‘Kristensen's comb’, and if polytomy is preferable to error, then Kristensen's phylogeny was not bested until genomic resources were brought to bear. However, Kristensen remained sceptical even after the publication of large transcriptomic (box 4) works [49,50] (NP Kristensen 2015, personal communication to R.G.B.).
Controversy was common in morphology-based insect phylogenetics, and results were strongly affected by the selection of characters, before very large and well documented datasets emerged in the twenty-first century. Boudreaux's 1979 book [69] on arthropod phylogeny, though criticized by some [43,63,70,71], was cited frequently by others [72–77]. Even though Boudreaux adopted the methods of phylogenetic systematics, according to Kristensen [78] his interpretations often differed from those of Hennig [10,11]. As in the case of the controversial Zoraptera [79], phylogenetic conclusions were often based on characters that were ancestral, ill-defined or homplasious. Jarmila Kukalová-Peck published a summarizing account of insect palaeontology [80], numerous specific studies on extant and extinct insects [81–83] and comprehensive analyses of characters of the wing base and wing venation [84,85]. She advocated the origin of wings from gill-like appendages [86,87] and proposed a clade Cercophora (Diplura + Insecta) for the first time [87]. Her groundplan approach challenged standard cladistic procedures [88] and was criticized by some authors [89]. Her phylogenetic hypotheses, usually based on wing characters, yielded some results inconsistent with earlier [11] and most recent concepts [50].
As in earlier attempts to classify insects (like Haeckel 1896), studies based entirely on wing venation [84] show the weakness of limited character systems, especially when strong functional constraints drive convergent evolution. Nevertheless, in-depth studies of specific body parts, organs or developmental stages can yield important insights. Examples are the circulatory system investigated by Günther Pass [90,91], the female genitalia of polyneopteran groups studied by Klaus Klass [92,93], and embryology, with important contributions made by Ryuichiro Machida and others [94,95]. Throughout the first decade of this century, it was more common in presentations to see these characters mapped onto molecular phylogenies than to have explicit, data-matrix-based phylogenies constructed from these systems. This is understandable, given the recognition that subsets of characters were only part of the whole picture, the general lack of coordination in taxon sampling, and the enormous effort involved in constructing a unified combined data matrix.
Classical Hennigian studies relied on detailed anatomical information obtained for few selected taxa with informal character discussions without data matrices [35–37,62,63,96–98]. These studies treated all taxa within larger groups as a single hypothetical ancestor, reducing characters to reconstructed groundplan states. Modern computer-based analysis is better suited to entering characters of individual representative species into data matrices. Even so, earlier computer-based studies extracted data from the literature [62] and coded entire orders with identical groundplan states [72,74,77]. This was almost inevitable, because thorough anatomical studies using microtome sectioning of a single species [99] could take years. In the early 2000s, new technologies such as micro-computed tomography (µCT) and computer-based three-dimensional reconstructions greatly accelerated the acquisition of high-quality anatomical data [19]. The coordinated efforts of international research teams, using both new and traditional methods [100,101], have yielded matrices of hundreds of characters from different body parts and life stages. For example, a study of Holometabola [48] contained 365 well-documented characters that corroborated current molecular phylogenies.
6. Insect morphology and cladistics
In the late 1970s and 1980s, cladistics ‘evolved’ as a transformed version of Hennigian phylogenetic systematics [102,103], arguably linked with the development of suitable computers and software programs. The Hennigian method of searching for sister taxa required great care in polarizing each character. Polarity in this context refers to the assignment of character states as either ancestral or derived. However, character polarity is automatically determined based on outgroups (box 1; i.e. rooting (box 2)) with computer-based analysis [104]. The first computer program capable of estimating phylogenies was Felsenstein's PHYLIP in 1980. Mickevitch and Farris were developing their program ‘PHYSIS’ near the same time and released it in 1982. It saw limited use, perhaps because of the $5000 price tag. Farris' updated program Hennig86 became available in 1989 [105], and Swofford's PAUP [106] was released free of charge the same year. Along with new molecular data, Whiting et al. [77] presented a morphological matrix-based analysis of most major insect groups in 1997, which was extended to include all hexapod orders by Wheeler et al. in 2001 [74]. That same year, Beutel & Gorb [72] presented a matrix-based morphological analysis of the entire Hexapoda. These morphological phylogenies were largely consistent with earlier hypotheses [10,11,28,31,62,63,64,78]. Wheeler's insect ordinal phylogeny [74] emphasized molecular data, but, without the morphological data, their results were largely unresolved and implausible [107].
7. The dawn of molecular systematics in the early 1990s—molecular work in the Sanger days
A number of studies in the late 1980s explored animal phylogeny, including insects, using direct RNA sequencing of the nuclear small subunit ribosomal RNA gene (18S rRNA) [108,109]. Turbeville et al.'s 1991 work [110] used parsimony (box 1), distance (box 2), and other methods [109]. Their distance analysis grouped the annelids with the molluscs as opposed to the previous assumption that annelids should group with the arthropods based on segmentation. They also recovered Pancrustacea, a group that unites the traditional crustaceans with hexapods. At the time, the Tracheata hypothesis (Myriapoda + Hexapoda) was heavily entrenched, and they suggested that the position of the crustaceans may have been the result of bias introduced by long-branch attraction, and the limited number of characters. Earlier work [108] also recovered Pancrustacea, and suggested that the annelids were distant from the arthropods. These works remind us to be careful of what we dismiss as ‘wrong’, because we now understand Pancrustacea to be strongly supported. Turbeville et al. [110] were aware of branch-length artefacts and alignment (box 3) ambiguity, and made careful, if arbitrary, decisions about data exclusion. Unlike some who followed them, they considered suboptimal trees to be worth discussing. However, given that few were impressed with confirming arthropod monophyly, and still fewer believed that crustaceans should group with Hexapoda, this study, as insightful as it was, did not become a model for future analyses.
Phylogenetic terms from the Sanger days.
Nucleotide compositional bias. When nucleotide frequencies stray significantly from 25% of each DNA base (A, C, T and G). This phenomenon is particularly problematic when bias differs among taxa.
Among-site rate variation (ASRV). When different sites along a sequence vary in their substitution rates. For example, when the substitution rates are higher for third codon positions than for second codon positions, this difference in rates is important to capture in model-based analyses, and argues against equally weighted parsimony. ASRV is also extremely problematic when it varies across lineages in a tree.
Multiple sequence alignment. The process of lining up DNA or amino acid data into columns of presumed homologous positions.
Consensus tree. A graph summarizing a set of trees by showing only clades which are shared among multiple equally favoured solutions or even multiple analyses. A strict consensus tree shows only those relationships found in all trees, whereas a majority-rule consensus depicts the most common resolution.
Sensitivity analysis. A means of exploring the robustness of a conclusion by altering the analytical details that influence it. For example, if one were interested in exploring how alignment parameters (like the penalty for inserting a gap in an alignment) influenced a phylogeny, one could change the input values to create new phylogenies from the new alignments and explore how the resulting trees differ.
Input parameters. Many complex analyses, such as alignment of DNA sequences or phylogeny reconstructions, require a priori specification of a number of parameters. Common input parameters include values for costs or ratios, or parameters of a specific evolutionary model. Input parameters are often derived from empirical data and can drastically alter phylogenetic results.
A 1984 review of insect molecular systematics by Berlocher [111] focused on allozyme gel electrophoresis studies, with discussions of methods used at the time. Before the invention of the polymerase chain reaction (PCR) in 1985, direct sequencing of rRNA was possible, but rare, and only one study of the molecular structure of rRNA, presenting 3 insect 5.8S sequences [112] was mentioned in the review. A 1988 paper by Simon [113] included a table of 30 molecular phylogenetics projects underway at the time, but all of them were as yet unpublished. Thus, the first molecular study that we are aware of that specifically addressed insect phylogeny was published in 1989, when Wheeler [114] discussed separate analyses of insect 18S sequences and restriction sites (box 2). The restriction site analysis (his figure 6) included more taxa than the DNA sequence tree, and supported Metapterygota (damselfly + Neoptera), and Neuropteroidea (beetle + lacewing). It seemed from this work that the 18S rRNA gene was a promising source of characters, especially given that restriction sites alone could result in a reasonable tree. In 1992, Carmean et al. [115] used 18S rRNA to explore relationships among holometabolous insect orders and noticed that flies had an elevated substitution rate, and long regions that had to be excluded from the analysis because they could not be aligned (box 3) with confidence. They surmised that the flies (Diptera) were being drawn to the root (box 2), and, thus, excluded them in most of their analyses. Pashley et al. in 1993 [116] published distance and parsimony analyses of a fragment of 18S rRNA from nine orders of Holometabola. They were able to recover Mecopterida and Amphiesmenoptera, but bootstrap support (box 2) for most groups was very low. Pashley et al. concluded that using different outgroups (box 1) yielded different topologies for poorly supported ingroup taxa. The failure of these analyses to converge on strongly supported results from a few taxa, and fragments of 18S is not surprising. This gene alone has never resolved relationships among all orders of Holometabola, even with many more taxa, although von Reumont's work [117] came very close with combined complete 18S and 28S but without the confounding Strepsiptera (see below).
Mitochondrial data were also explored in the early days of Sanger sequencing (box 2). Liu and Beckenbach [118] explored the mitochondrial cytochrome oxydase II (COII) gene in 10 orders of insects, using a genetic-distance-based analysis [119] and parsimony. Trees from various analyses grouped the cockroach and the termite, and the three species of Hymenoptera (ichneumonid wasp, bee, and ant), and not much else. In a study of arthropod phylogeny, a small fragment of mitochondrial 12S rRNA gene was analysed, and it was proposed that onychophorans (velvet worms), are in fact modified arthropods [120]. As onychophorans are generally considered to be arthropod outgroups, this study was published with fanfare in the journal Science. The statement that ‘These data demonstrate that 12S … can resolve arthropod relationships…’ is strongly contradicted by the highly unusual (and since rejected) phylogeny they recovered.
Through the Sanger sequencing period, molecular phylogenetics focused largely on 18S, 28S and a few mitochondrial genes, mostly 12S rRNA, 16S rRNA and COI. However, rRNAs were difficult to align [121,122] and model [123], and it seemed that the mitochondrial genes were biased and full of misleading signal [118]. Single-copy nuclear genes were seen as a solution, but remained difficult to sequence. The standard markers were easier to amplify, because universal primers were available [123,124,125], and both mitochondrial genes and nuclear rRNAs were present in multiple copies in every cell. A review in 2000 [126] called for coordinated efforts in selecting genes that were compatible across studies, and supported the continued use of 18S rRNA and commonly sequenced mitochondrial markers. In a contrasting opinion, in order to move beyond rRNA and mitochondrial genes, a group of workers at the University of Maryland embarked on a programme to locate and sequence single-copy nuclear protein-coding genes [127]. Of 14 ‘promising candidates’ they identified, several (EF-1α, DDC, POLII and, to a lesser extent, PEPCK) that saw extensive use in insect intraordinal phylogenetics. They would continue to develop useful protocols for amplifying genes such as wingless, CAD rudimentary and others [128–135]. Additional contributions to the arsenal of nuclear genes for insect phylogenetics soon followed [136–138]. Histone H3 and U2 snRNAs were examined [139,140], with the former used extensively despite the fact that neither gene could recover any reasonable higher level groups [141]. Practically all major higher level insect phylogenetic studies in the past decade have relied, at least in part, on single-copy nuclear genes, and these have now become the dominant markers in transcriptome (box 4) analyses. The markers developed by the Maryland workers and others were put to good use across arthropods, among a few orders [142], and within orders such as Lepidoptera, Hymenoptera, Diptera and Coleoptera (see below). However, they were not applied broadly across orders until Wiegmann's 2009 work [47].
Modern terms in molecular phylogenetics.
High-throughput sequencing. A variety of methods that generate huge amounts of DNA sequence data. These methods were called ‘next generation’ during the Sanger sequencing period.
Transcriptome. All the genes that are active in a cell, tissue, or whole organism at the time a specimen is collected. These genes are sequenced by isolating messenger RNA from the cells.
Assembly. High-throughput sequencing typically produces millions of fragments of relatively short sequences. These sequences are aligned, and stitched together into longer fragments called ‘contigs’, which are then sorted into orthologues.
Orthology prediction. An orthologue is a homologous gene. Gene duplication produces non-homologous, similar-looking genes (paralogues) that can disrupt phylogenetic analyses. Orthology prediction is the process of identifying orthologues and distinguishing them from paralogues.
Data masking. Identifying ‘bad’ data, and throwing it out. Problematic data come from a variety of sources. Data can be misaligned, randomized, mostly missing or violate model assumptions.
Partitions. Subsets of the data. Because models are used to analyse data, it is common to subdivide the data into subsets that are assumed to share similar properties so that the models are a closer match to biological reality.
Quartet mapping. A method of branch support that examines a particular node by sampling replicates of four randomly selected taxa that surround the node of interest. The results summarize the percentage of times that these quartets recover alternative topologies.
Pipeline. A sequence of analytical procedures. In computer science, a pipeline is a chain of computer programs that perform analytical steps in series.
Taken as a whole, it is not difficult to see why morphologists would be less than excited by the state of molecular phylogenetics in the early 1990s. In historical context, this was a time when university hiring priorities favoured molecular workers who could sequence a couple of hundred nucleotides from backyard insects, and ‘discover’ relationships that were either already widely accepted or hard to believe. Grant money seemed to be reserved for molecular work. Still, morphological workers seemed to be at an impasse. They agreed on the general outlines of Hennig and Kristensen, and had established Dictyoptera, Mecopterida, Amphiesmenoptera, Antliophora and Neuropterida as monophyletic. However, they were unable to resolve the relationships among the entognathous, palaeopteran, polyneopteran, or holometabolous orders.
The development of PCR made possible the rapid collection of nucleotide sequence data throughout the 1990's and beyond [143,144]. Most molecular workers at this time recognized the limitations of their own data, and were proposing ways to address them. Two early reviews [123,145] discussed strategies for modelling DNA to account for known biases in the way it evolves. Swofford and co-workers influential chapters in the ‘Molecular Systematics' books [146,147] had laid a foundation for understanding the analytical issues, and the programs PAUP* [148] and PAML [149] were available for running model-based (likelihood (box 2)) analyses at a time when computers were finally up to the task of implementing complex substitution models. PAUP is an acronym for ‘Phylogenetic Analysis Using Parsimony’, so the asterisk after the new release referred to ‘and other methods’, such as likelihood and distance. The pull-down menu (graphical user interface or GUI) in PAUP* was an ideal platform for newcomers to learn the complexities of models of DNA evolution and statistically based likelihood phylogeny-building methods. At this time, the need to accommodate biases with either models or differential weighting were becoming obvious [146]. By the mid-1990s the field of insect molecular phylogenetics looked promising, but it had begun to splinter into camps based on analytical methods. There was a brief honeymoon where the logic of cladistic taxonomy was universally adopted, but then co-opted by some who conflated cladistics with parsimony, as they transitioned from intuitive Hennigian methods to computer-based parsimony analyses.
8. The problem with ‘the Strepsiptera problem’: 1995–2010
There was probably no question that occupied the minds of insect systematists more during the late 1990s than the ‘Strepsiptera problem’. This is surprising because strepsipterans are neither diverse nor conspicuous. However, like many parasites with highly modified structural features, these fascinating and unusual insects were difficult to place morphologically. Their ribosomal data would prove to be the battleground over which likelihood and parsimony practitioners would argue, which in turn helped to reveal the problems inherent in parsimony. Most morphological studies placed Strepsiptera (figure 4) as the sister taxon to the beetles (Coleoptera), based on hindwing flight (posteromotorism), and a few other characters [72,150–152], or within a subgroup of Coleoptera [153] (see also [154]). The first of the molecular-based studies addressing this, (but without published data or analytical detail) was a Scientific Correspondence that appeared in Nature in 1994 by Whiting & Wheeler [155]. They proposed that Strepsiptera were the sister taxon of Diptera (flies; the two groups combined named Halteria [77]). Halteria refers to halteres; the gyroscopic reduced hindwing stubs found in flies that are superficially similar to strepsipteran forewings. They surmised that the grouping of flies with strepsipterans was in itself evidence of a homeotic mutational transformation that resulted in halteres flipping from the third thoracic segment in flies to the second in strepsipterans (figure 4). Most morphologists doubted this assertion. Scepticism came from the molecular perspective as well. Carmean & Crespi [156] responded almost immediately that ‘long branches attract flies’, which was unambiguously demonstrated by Huelsenbeck [157], with a likelihood (box 2) analysis of available data [156] (the Whiting data were not yet public). Chalwatzis et al. [158] were the first team to actually publish an analysis of this problem and make available their data. Like Whiting and Wheeler, they used 18S rRNA sequences and recovered Halteria. In a 1996 follow-up work [159], increasing the number of taxa to 26 and including all holometabolous orders, they again recovered Halteria. In addition to parsimony, these authors used the neighbour-joining distance method with a model [160] designed to accommodate nucleotide compositional bias (box 3) among lineages and among-site rate variation (box 3). Strepsipteran 18S was found to be about 1000 nucleotides longer than the next longest 18S and shared extreme and similar AT nucleotide compositional bias with Diptera. All their analyses favoured Halteria, including those designed to correct for the bias they had observed. However, in an analysis they did not show, when site-specific rates were used to correct for among-site rate variation [161,162], the bootstrap (box 2) value for Halteria dropped from 100% to 77%. They cautioned that 18S could be artificially clustering long-branched taxa and looked forward to investigating other genes not linked to rRNA to test their findings.
The largest dataset of the 1990s exploring the phylogeny of Holometabola was presented in 1997 by Whiting et al. [77]. Approximately 1100 18S and 400 28S rRNA positions were aligned with the multiple sequence alignment (box 3) program Malign [163], and analysed with parsimony. Molecular data were then combined with morphological data taken from the literature. Sensitivity to alignment was explored by evaluating trees both with, and without hypervariable regions. Both beetles and neuropteroids were polyphyletic, owing to contamination. The paper is best remembered for its recovery of Halteria, a hypothesis that they would vigorously defend [74,164–167] until morphology [48,168] and additional genes [47,154,169,170] overturned it 15 years after it had been proposed.
9. Long-branch distraction?
Whiting [77,164] rejected the suggestion that Halteria was an artefact of long-branch attraction (box 2), arguing that 18S and 28S corroborated one another. However, these are not independent genes but rather different regions of the same transcript. Countering Whiting's other arguments, Huelsenbeck showed in 1998 [171] that the length of the branches leading to Strepsiptera and Diptera were ‘virtually unparalleled in phylogenetic analysis’. Huelsenbeck's analysis was among the first to apply likelihood to a large number of insect orders. Whiting [77] had argued that the branch leading to the amphiesmenopterans was ‘not far out of range’ of those leading to the Strepsiptera and Diptera. However, this statement missed the central tenet on which long-branch attraction is based. In parsimony analyses, as independent changes are transferred away from the terminal branches leading to both Strepsiptera and Diptera where they occurred, to the internode (box 2) that falsely links them together, the observed terminal branch lengths are underestimates of the real number of independent changes. In other words, parsimony takes two independent changes and assumes they are shared-derived character states. By removing either Strepsiptera or Diptera from the analysis, each remained in the same position relative to the remaining taxa. Their argument was that Strepsiptera could not have been attracted to Diptera given that it ends up in the same place in the tree when Diptera is removed [165]. However, removing either taxon simply caused the remaining long-branch taxon to attract to the next longest branch, the Amphiesmenoptera, a branch that they recognized as almost as long [77] as those leading to Strepsiptera and Diptera (although they had underestimated these branch lengths). Huelsenbeck [171] showed that, given the taxon sampling used by Whiting et al. [77], the branches leading to both Strepsiptera and Diptera were long enough to attract one another with parsimony, and that likelihood analyses could not distinguish among hypotheses. It was recognized that taxa at the end of long branches may actually be sister groups [107,157,171,172], but that the rRNA data in hand could not support any conclusion including Halteria. Friedrich & Tautz [173], confirming the observations of Chalwatzis [158,159], showed in 1998 that there had been an extreme change in both substitution rate and compositional bias (box 3) in the stem dipteran lineage that would pose problems for phylogenetic analyses. These biases were further explored in 2000 by Steel et al. [172] and others [174]. Hwang et al. [189] sequenced additional large subunit rRNA fragments in order to test the Halteria hypothesis. Their parsimony analyses recovered Halteria, which they attributed to long-branch attraction (box 2), whereas their likelihood analyses placed the strepsipteran with the scorpionfly. All these authors recommended that the Halteria hypothesis be dropped. So why did it take phylogenomics to settle the issue? It didn't, and this is not the wisdom of hindsight. For likelihood practitioners, it was settled in 1998. It had been conclusively demonstrated that Halteria was the result of inappropriate methods and obvious predictable bias. Halteria has only been found from the analysis of nuclear rRNA data, and contradicted by every other source of data [47,50,154,168,175]. Genomic [175,176] and transcriptomic (box 4) [50,177] analyses now leave very little room for debate: Strepsiptera belongs as sister to the beetles. The debate was never really about Halteria, but rather, about the philosophical merits of parsimony versus likelihood.
10. Alignment issues
In addition to disagreements over the merits of parsimony, the methods of nucleotide alignment (box 3) played an important role in insect phylogenetics during the Sanger sequencing period [107,163,178–186]. The definition of cladistics had been transformed, and, in the new sense, was characterized by strict and exclusive adherence to parsimony analyses. The first work in insect molecular systematics that covered all orders came from a group of cladists (in the new sense) who were centred at the American Museum of Natural History in New York [74]. They extended the dataset of Whiting with the same rRNA fragments as in the 1997 work [77], but with additional taxa, particularly outside Holometabola. The principal analytical difference was the implementation of simultaneous alignment and tree building [187]. They called this method ‘direct optimization’ when implemented by their program, POY [179]. The molecular data by themselves, presented in their figs 12a 13, and 14, suggested many implausible relationships [107]. However, it seems that the morphological data provided a stabilizing scaffold that mediated the misbehaviour of the molecular data. In order to explore the influence of different analytical assumptions, they presented six combined data trees. The analysis that minimized incongruence among datasets (their fig. 11) would seem to be the favoured hypothesis, although this was not explicitly stated. The summary trees of all assumption sets, shown in their fig. 18a,b, were largely unresolved consensus trees (box 3). However, often their results are cited as their fig. 20, which did not come from an analysis but rather was a ‘discussion tree’ created from nodes the authors favoured from different datasets. The deepest parts of their trees seemed robust to analytical assumptions, whereas relationships among polyneopterans and holometabolans were unstable. Despite published papers that pointed toward branch effects and compositional bias (box 3) for these data [137,156,157,171,172,188,189], they did not take these problems into consideration, favouring parsimony on philosophical grounds.
Many of the differences among phylogenetic hypotheses were the result of differing analytical approaches and ambiguous alignment (box 3) of rRNA data. The history of alignment disputes has been described in detail elsewhere [181,182,185], and some researchers likely turned to nuclear single-copy genes simply to avoid the problems of rRNA alignment altogether, and perhaps the tedious bickering from both sides of this issue [190]. POY has since fallen out of favour with systematists, owing to numerous and diverse criticisms [121,183,184,191], although the possibility of simultaneous alignment and tree building remains [192]. The problem with model-based simultaneous alignment and tree building is that it is difficult to create a biologically reasonable model for gaps.
11. The sensibility of sensitivity
By 2001, the insect systematics community was strongly divided into parsimony, and likelihood camps. It was another set of parallel universes, with different journals (Systematic Biology versus Cladistics), different heroes (Felsenstein versus Farris), different branch support measures (bootstraps versus jack-knifing and Bremer support (box 2)) and even different computer systems (Mac versus PC) brought about by the platforms of their different programs (PAUP versus Winclada, NONA, and TNT). Many in both camps basically dismissed the ideas from the other side as flawed and without merit. Most morphologists found themselves in the parsimony camp, likely because of tradition and the fact that morphological characters are less amenable to modelling than molecular characters. Much of the error in parsimony analysis could have been mediated by differential weighting of characters, upweighting slow sites, and downweighting fast sites [123,193]. Morphologists had always weighted their data, if only by selection of characters that they deemed reliable. However, with DNA there was little interest in differential weighting, as one side rejected it based on the insistence that equal weights were assumption free, and the other side preferred likelihood, because models mimic differential weights with the added benefit of being grounded in statistics [194]. Weights and other parameters (box 3) upon which phylogenetic conclusions depend must be selected by the user. If their selection is arbitrary, then subjectivity remains, but it is transferred from a thinking person to a machine [181,182,185]. Many in the molecular-parsimony camp were dedicated to POY analyses, and believed that they were removing as much subjectivity as possible. In order to deal with the problem of objectively selecting analytical parameters, they developed a brand of sensitivity analyses [74,195] that was based on incongruence length difference (ILD) tests [196]. ILD testing involves comparing subdivisions of the data with combined data, seeking parameters that minimize incongruence. Criticism of ILD testing is beyond the scope of this review, but can be found in many works [182,197–201]. However, even if ILD tests were legitimate, then one must decide which among many parameters should be evaluated, each with an infinite space to explore, with each influencing the behaviour of the others [182]. Grant & Kluge [202], in a particularly radical application of their own view of epistemological consistency, reject the whole idea of sensitivity analyses (box 3), and suggest that all parameters should be equal, and set to 1 on philosophical grounds. Ogden & Whiting [203] applied sensitivity analysis to the ‘Palaeoptera problem’—the phylogenetic positions of dragonflies and damselflies (Odonata) and Mayflies (Ephemeroptera) relative to insects that have the ability to fold their wings (Neoptera). They showed that the results were indeed sensitive to input parameters. In a justification for using a single analytical method, and counter to exploring the influence of the application of model-based methods, they stated that they ‘…do not consider congruence among different methodologies to be a suitable measure of robustness because agreement among inferior methods is nebulous at best’. It seems that this attitude was shared by both sides, as model-based analyses were not explored by the cladistics group until Terry & Whiting in 2005 [204], and parsimony analyses were virtually abandoned by the practitioners of likelihood. We see now that short, ancient internodes (box 2) are always sensitive to assumptions and input parameters. Thus, the nodes unseen by Börner in 1904 collapse with sensitivity analyses (box 3) as it was applied, as seen, for example, in fig. 18 in Wheeler et al. [74].
Concerning the relationships among insect orders, the entire Sanger period provided few if any new insights that were widely agreed upon. Cockroach paraphyly, an apparent exception, had been suggested based on morphology [205]. Part of the lack of resolution came about, because the parsimony and likelihood schools were so far apart, and non-specialists could not choose between them. Even a hypothesis that was supported by practitioners on both sides of the analytical divide—Nonoculata (Protura + Diplura; appendix A)—now seems to be an error (but see [206]). In addition, the common result of finding snow fleas (Mecoptera: Boreidae) closer to the fleas than other mecopterans is now in question. These are very difficult phylogenetic problems. The current prevailing opinion is that model-based analyses outperform parsimony even when parsimony is weighted to be more realistic [123] (despite the editorial in 2016 in Cladistics [207]). Accepting this premise, the philosophically driven parsimony-based insect molecular phylogenies that dominated the literature in the 1990s and 2000s were an unfortunate detour, especially when compounded by the failure to recognize the errors resulting from DNA compositional bias, non-homogeneous substitution rates, alignment error and the inconsistencies of rRNA analysis with POY [107,171,172, 181–184,208].
12. The likelihood camp
The basic principles and performance of likelihood [209–215] were laid out by Felsenstein long before they became standard practice. Likelihood was first introduced into phylogenetic systematics by Cavalli-Sforza and Edwards in 1965 but it was not widely applied because user-friendly programs were not available until PHYLIP was developed in 1980 [216]. Even after likelihood programs were available, it took a while for people to develop an understanding of how models of evolution could lead to an estimate of phylogeny. Parsimony was far easier to grasp. Many were uncomfortable with the number of assumptions required for model-based analyses. However, although it was claimed that the assumptions required for equally weighted parsimony were fewer or non-existent, it is clear that they were simply not defined. If they were defined, then they would be exceedingly complex and unacceptably unrealistic [121,201]. Even if there were fewer assumptions, these few would still lead to error with certainty under common branch length combinations [188]. However, a major obstacle to using likelihood was that it was difficult to analyse more than 10 taxa in a reasonable time frame. Sophisticated models of evolution could not be implemented until computers gradually gained the speed to analyse datasets of typical size, more than 10 years after PHYLIP was introduced. Throughout the 1990s, as computational speed increased, phylogenetic methods based on likelihood as an optimality criterion grew in importance and implementation. In 2000, it could still take weeks on a desktop computer to analyse 500 nucleotides for 50 taxa. Bootstrapping or any kind of branch support (box 2) was difficult if not impossible for likelihood analyses with more than 25 taxa until fast maximum-likelihood programs were developed—PhyML [217], Garli [218] and RAxML [219].
The motivation for implementing likelihood was strong. Felsenstein had demonstrated in 1978 [188] that, under parsimony with some branch length ratios, the addition of data would strengthen support for the wrong tree. He speculated that parsimony would work if rates of evolution were low or sufficiently equal among lineages. Hendy & Penny [220] extended Felsenstein's work to show that neither of these conditions for the success of parsimony would hold once the number of taxa exceeds four. They concluded that rather than unequal rates it was the long branches that were the problem, and introduced the concept of long-branch attraction. The idea that adding more data could not overcome this bias was hard to accept, and we still find the idea expressed that more genes or increased taxon sampling is a panacea. Into the 1980s, most cladists were still basking in the glow of defeating the numerical taxonomists (whom they called pheneticists). Although likelihood is clearly based on individual characters (like parsimony), its statistical underpinnings were incorrectly assumed by some cladists to link it to phenetics (box 2). This is ironic, because, as Tuffley & Steele [221] demonstrated, under certain (unrealistic) models of evolution, likelihood can be equated with parsimony.
Bayesian analysis, which shares many properties with likelihood, was introduced into phylogenetics in 1999 [222], and could be implemented in the user-friendly program MrBayes [223]. By 2001, many in the likelihood school rapidly adopted MrBayes, because it calculated branch support in the form of posterior probabilities at lightning speed [141]. It was later realized that much longer Bayesian runs were necessary to be sure that the program had converged on the optimal answer, especially when data required complex models of evolution. However, the problems with model-based analysis were not entirely because of ignorance or the lack of computing power. The influence of long-branch attraction was debated [164,165,224], but its ubiquity was not fully understood. Models that did not accommodate key elements of reality, such as among-site rate variation (box 3), could be as error prone as parsimony, without parsimony's comfortable philosophical footing based on the perception that it minimized unjustified assumptions. It was uncomfortable to use a method so dependent on models, if one could not justify which model to use. Model selection became a major focus of phylogenetic studies [225–227]. At first models of evolution were tested manually using likelihood ratio tests [228,229], but this became automated in 1998 with ‘Modeltest’ [230]. It seemed that this program almost always suggested the most complex model, which led to the development of decision theory (reviewed in Sullivan & Joyce [231]), including a stronger penalty for increasing the number of parameters (box 3).
In hindsight, if we are to judge by current standards, and our ability to assess accuracy in the light of phylogenomic data, likelihood analyses were both more accurate and philosophically grounded. Most early likelihood practitioners in entomology confined themselves to intraordinal relationships, and their work has become relatively robust, as datasets have become larger. Friedrich & Tautz, in 1995 [232], were among the first to use likelihood to estimate deep arthropod relationships with PHYLIP [216]. Their analysis included three hexapods, and recovered Pancrustacea, and crustacean paraphyly. Likelihood (among other methods) was also used by von Dohlen & Moran [233] in 1995 to demonstrate the paraphyly of ‘Homoptera’. Frati et al. [228] used maximum-likelihood analyses of mitochondrial COII gene data in 1997 to examine relationships among springtails, and demonstrated that including a correction for among-site rate variation (box 3) in the analysis had more of an effect on likelihood scores than the substitution models themselves. Flook & Rowell [234] used likelihood methods to explore the properties of mitochondrial data among orthopterans. Whitfield & Cameron [235] found in 1998 that likelihood outperformed parsimony in their study of hymenopteran 16S rRNA. Lo et al. [236] used likelihood methods in 2000 to demonstrate the paraphyly of cockroaches. In 2001, Kjer et al. were the first entomologists to use Bayesian methods in their study of caddisfly (Trichoptera) phylogeny, and Kjer [107] was the first to include many insect orders with Bayesian methods.
It was not until the mid-2000s that consensus among entomologists in the USA swung towards likelihood analyses for molecular data, but parsimony analyses are still being published because of cultural/historical factors, and is still being actively favoured by the journal Cladistics [207]. It can take a long time for attitudes to shift, and sometimes recollection is subject to ‘retrospective meaning change’ (see discussion by Hull [237]). As with debates over creationism, or climate change, the fact that there are two sides to an issue does not mean that both sides are equally supported. Parsimony for molecular data seems to be supported by faith. Sometimes progress in science comes, not from evidence or flashes of insight, but through strong personalities fading into retirement. (Paraphrasing Max Planck: ‘Science advances one funeral at a time'.)
13. The dominance of ribosomal RNA
Although many papers included small fragments of 28S or histone H3, it was the 18S that dominated results from the Sanger sequencing period, sometimes stabilized by morphological characters [107]. This was partially owing to historical artefact and partially due to the ease of amplifying and sequencing nuclear rRNA. The 18S gene suffers from extensive among-site rate variation [193] and severe alignment problems within some regions, whereas (unlike the 28S) the alignable regions are practically invariant. Kjer [107] explored the properties of the 18S, structurally aligned, using a model-based analysis that accommodated rRNA covariation [238]. His phylogeny was much closer to current consensus than previous parsimony analyses of 18S. In 2005, a large insect phylogeny was presented by Terry & Whiting [207], using histone H3, larger portions of the 18S and 28S, and a modified morphological data matrix from Wheeler et al. [74]. They focused on polyneopterans, including Mantophasmatodea for the first time, and included a Bayesian phylogeny along with their POY-based parsimony [204]. Their Bayesian analysis was a great leap forward. They recovered many of the nodes we now find with larger datasets; many for the first time with molecular data, including Xenonomia, Eukinolabia and Haplocerata (appendix A), which they named, as well as Polyneoptera, Neuropteroidea (Strepsiptera was not included) and Antliophora. In a counterpoint to POY-based analyses, Kjer et al. [141] provided a review of the data of the time, with a commentary on methods. They reported the results of a 15 000 nucleotide multi-gene supermatrix, put together from complete 18S, a large fragment of 28S, EF-1α, histone H3 and mitochondrial 12S, 16S, COI and COII, along with 170 morphological characters from the older sources, such as Hennig, and Kristensen [11,62–64]. Their results came very close to our current consensus, particularly within the polyneopterans. In all Kjer's analyses, Strepsiptera and Zoraptera were excluded, because these taxa exhibited extreme substitution rate accelerations in their rRNA. He was also suspicious of the published Zoraptera sequences. Zoraptera was resequenced and a modified structural alignment [107] was used by Yoshizawa & Johnson [75] in order to place this difficult taxon. They found it to be sister to Dictyoptera, as did Ishiwata [170] with nuclear single-copy genes. They cited morphological support for this relationship from Boudreaux [69] and Kukalová-Peck [82]. Misof's group [239] provided an insect-specific secondary structural model, and re-evaluated Kjer's [107] analysis of 18S, with increased taxon sampling. They found similar results to those from earlier structural alignments although they found Zoraptera grouped with stoneflies (Plecoptera). This analysis included Strepsiptera, which, as in other likelihood analyses of rRNA [157,189] did not group with Diptera, but instead, in this case as sister to an implausible Diptera + ‘Coleoptera’ group, with the long-branch Diptera acting as a second internal root that rendered beetles paraphyletic. For the first time since 1997 [77], the molecular data recovered Hymenoptera as sister to the rest of Holometabola (appendix A; Aparaglossata). The most thorough exploration of rRNA-based insect phylogeny was completed by von Reumont et al. in 2009 [117]. They used an automated alignment algorithm that incorporated rRNA secondary structural information [186], eliminated randomized sites (phylogenetic noise) with the program Aliscore [240], and used more realistic substitution models. Thus, none of the previous criticisms over manual manipulation of alignments and manual data exclusion could be applied to this study, recovering Nonoculata, Ectognatha, Dicondylia, Pterygota, Chiastomyaria, Neoptera, Holometabola, Aparaglossata, Amphiesmenoptera and Mecopterida.
Published phylograms (trees with branch lengths proportional to the number of estimated substitutions) [75,117,141,171,172] illustrate the extreme heterogeneity of rRNA substitution rates among orders, and this property causes problems with standard methods [172], even under likelihood. Protura and Diplura share extreme branch lengths relative to neighbouring Collembola and Archaeognatha, and the rRNA of both has extremely long regions of hypervariability that are difficult to align [241]. Phylograms show that Zoraptera, Strepsiptera and Diptera are also extreme. Odonata evolve more slowly than their neighbours in the tree. Ribosomal RNA analyses frequently recover Nonoculata [74,75,107,117,239,242–245] Chiastomyaria [117,141,239], Dermaptera sister to Plecoptera [107,117,239], and mecopteran paraphyly [74,75,77,107,117,141,246]. The consistency of these results despite the differences in alignment and optimality criteria indicate that rRNA supports these relationships when analysed with existing methods, even though much larger datasets now contradict Nonoculata and mecopteran paraphyly.
14. Other types of data
Besides rRNAs and the few nuclear protein-coding gene studies, there were other novel character systems explored for insect phylogenetics, such as locations of introns, and mitochondrial gene order. Rokas et al. reported in 1999 that an insertion in a homeobox gene [247], shared by Diptera and Lepidoptera, was not found in Strepsiptera, contradicting Halteria. Carapelli et al. [248] noted that Collembola and Diplura shared the loss of an intron within EF-1α. Intron positions in EF-2 were mapped [249], showing that Coleoptera, Lepidoptera and Diptera shared a derived arrangement that was absent in Hymenoptera, predicting our current understanding. A survey of intron positions in EF-1α [250], found a remarkable tree from only six informative characters, but intron positions did show homoplasy (box 2), largely because of independent loss. A study of ecdysone receptors by Bonneton et al. [251] showed a significant rate acceleration that countered Halteria, as did their sequence analysis. Predel & Roth put their analysis of neuropeptides to use in studies of cockroaches, grasshoppers and Mantophasmatodea [252–255]. Xie et al. tabulated the distributions and lengths of 18S hypervariable regions [241], and they reported that some of them could be used as synapomorphies for insect groups. In addition to updating a secondary structural model for insects, they found that Zoraptera and Dermaptera shared the greatest number of hypervariable regions of identical lengths. Boore et al. [256] examined mitochondrial gene order among arthropods in 1995, and they found that Pancrustacea was supported by a mitochondrial gene order character [257]. After this discovery, it was hoped that mitochondrial gene order might help resolve difficult nodes among insect orders, because it was assumed that gene order was highly conserved and unlikely to be homoplastic (box 2). However, the most controversial internodes are likely to be short. We can think of internodes as targets where the size of the target is proportional to the length of time an ancestral lineage exists before it splits. As in archery, small targets are hard to hit. Thus, the probability of hitting extremely short internodes with extremely rare events is extremely low. An understanding of this phenomenon is currently important in genomic studies, where it is hoped that, with an abundance of characters, extremely short internodes may be hit by extremely rare changes in the structure of genomes. Unfortunately, mitochondrial gene order is remarkably conservative in insects, except within Paraneoptera [258–261] and Hymenoptera [262,263], with groups supported by changes in gene order, summarized in a recent review by Cameron [190].
15. Mitochondrial genomes
Mitochondrial data have been the subject of two recent reviews [190,264]. The accumulation of mitochondrial genomes continued through the 2000s [265–269], at a slow pace, but picked up rapidly after 2003 with concerted efforts from Cameron, Song, and Whiting [190]. Currently, whole mtDNA genomes are accumulating very rapidly because they can be efficiently targeted with high-throughput (box 4) methods [270], and are often recoverable as accidental ‘by-catch’ in high-throughput sequencing. It was not until preliminary results from the full-scale efforts to sequence entire mitochondrial genomes were published [271–274] that the extent of the problems with mitochondrial data became clear. Cameron and others concluded that mitochondrial data were promising, but that nucleotide compositional bias among lineages, unequal substitution rates among groups, and other long-branch effects must be carefully considered. Many of the relationships recovered with complete mitochondrial genomes were implausible, and they recommended that mitochondrial genomes be combined with other sources of data. Talavera and Vila found the same problems in 2011 [275], and proposed that deep nodes cannot be reconstructed with the methods of the time (which included Bayesian and likelihood analyses). Simon & Hadrys [274] recommended a similarly cautious view, finding many implausible relationships among orders even when using dense taxon sampling, and careful modelling. Chen et al. [276] found that including a projapygid helped recover dipluran monophyly, but they were still unable to recover hexapod monophyly with extensive taxon sampling among basal hexapods and arthropod outgroups. Cameron's more optimistic review of the phylogenetic implications of insect mitochondrial genomics, summarized model violations and made thorough recommendations for the appropriate treatment of mitochondrial genomes for phylogenetics. The issue of whether mitochondrial data are ‘good’ or ‘bad’ is clearly a gross oversimplification. Many ancient nodes that are accepted and corroborated by other data are recovered from mitochondrial data [274,277–281], and many relationships among polyneopterans are shared between nuclear and mitochondrial analyses [190]. For example, using mtDNA genomes, Wan et al. [281] found many of the nodes that Misof et al. [50] recovered from transcriptomes (box 4) and some of these nodes (figure 5: P,Q,R) have only rarely been seen before. Mitochondrial data consistently recover Mantophasmatodea with Phasmatodea [273,281]. Cameron et al. [282] found Megaloptera sister to Neuroptera, reflecting the results from transcriptomes [50]. The number of nucleotides in any analysis is strongly correlated with branch support. However, as mitochondrial genes are all linked and thus inherited as a unit, once the gene tree is accurately recovered, there is little more in terms of corroboration that the full mitochondrial genomes can add. The motivations and disagreements today, in the era of ‘big data’ phylogenomics, are sometimes centred around those who advocate for more data (in terms of both longer sequences and more taxa), and those who advocate ‘better data’. This disagreement has been with us since the beginning of molecular systematics, and it misses the point that more data, better data and better models are all good things.
16. Work on individual orders
This review has given short shrift to the vast majority of insect phylogenetics papers because of our focus on works addressing higher-level insect phylogeny. Given their almost unimaginable diversity, it is impossible for any individual to be considered an expert for all Hexapoda, and most workers spend their careers exploring particular groups. Here we list a sample of the recent advances from various authors in the phylogeny of Odonata [283–290,292], Ephemeroptera [278,293,294], Plecoptera [295], Dermaptera [296], Embioptera [297], Phasmatodea [298], Dictyoptera [236,279,299–301], Mantodea [302], Orthoptera [234,303–306], Hemiptera [307,308], Psocodea [309], Hymenoptera [310–317], Neuropterida [318–320], Coleoptera [321–323], Diptera [324,325], Lepidoptera [326–334], Trichoptera [335–337], Mecoptera [246,338,339] and Siphonaptera [340,341].
17. Beyond the standard toolbox: multiple genes, transcriptomes and genomes
Although many useful studies are still published with a few genes and morphology, phylogenetics today frequently involves the analysis of very large datasets. Savard et al. [342], in an early use of genomic phylogenetic resources in 2006, analysed 185 nuclear genes from four holometabolous orders, rooted with a grasshopper and an aphid, found results that are consistent with our current best estimates: (Hymenoptera, (Coleoptera, (Lepidoptera, Diptera))). At that time, most studies found Hymenoptera to be weakly supported as sister to Mecopterida; the group including Mecoptera, Siphonaptera, Diptera, Trichoptera and Lepidoptera, so the strong support for their alternative result led the authors to suggest that large datasets could resolve long-standing controversies in insect phylogenies. One of the first studies on Holometabola to break free of the standard rRNA and mitochondrial genes for interordinal analyses reported results from six single-copy nuclear genes [47]. A similar study using nine nuclear genes [154], found nearly identical results, both predicting our current understanding of relationships within Holometabola. Even though the datasets were no larger than previous rRNA-dominated analyses, and significantly smaller than the transcriptomic analyses to come (figure 5), the fact that both papers rejected Halteria independent of rRNA gave them extra impact. Three new nuclear protein-coding genes (DPD1, RPB1 and RPB2) were used in 2011 [170], further rejecting Halteria. Sasaki et al. [343] sequenced over 10 000 nucleotides from these same three genes, and focused their attention on the early splits among hexapods, with significant arthropod outgroups, and polyneopterans. They recovered the unusual result of (Protura, ((Collembola, Diplura), Insecta)), which has not been subsequently corroborated.
18. Data-mining and big, automated phylogeny pipelines
Behaviourists, ecologists and other biologists rely on phyloenetic trees to understand the evolution of complex characteristics. GenBank is data-rich, and the temptation to create pipelines (box 4) to download, combine, filter and analyse these data to produce a tree is strong. Building upon work by Hunt et al. [344] to generate large datasets from public databases, Peters et al. [312] developed a ‘proof-of-concept’ pipeline that mined GenBank for data from Hymenoptera, in order to construct a phylogeny with over 1000 taxa. The concept worked, but the phylogeny suffered from the quality of the original data in GenBank. Bocak et al. [322] constructed a phylogeny with public databases for more than 8000 beetle species. Again, this study proved that such a thing can work, and supported several disputed internal relationships. Zhou et al. [345] produced a phylogeny of over 16 000 barcode haplotypes from Trichoptera, but this study differed in that constraints were used to insulate the phylogeny from predictable errors. These studies provide evidence that producing huge phylogenies from public databases is feasible. However, based on our experience with genomic and morphological data, we caution that without analytical expertise for the specific properties of the data, as well as the insights of taxonomic specialists for a particular group of insects, it is impossible to reconstruct and evaluate the plausibility of phylogenetic relationships. This idea exemplifies the balance between skilled analyses that produce reasonable phylogenies, and the concern that unjustified or capricious decisions could bias phylogenetic conclusions.
18.1. The Palaeoptera problem revisited
An early transcriptomic analysis of seven pterygote orders, rooted with Collembola, grouped the mayflies with Neoptera (the Chiastomyaria hypothesis; appendix A) [346], but the Palaeoptera problem was far from solved, because Regier et al. [347], in a study focusing on arthropods, recovered a contradictory node Palaeoptera. Thomas et al. [348] evaluated the standard Sanger data in 2013, and found support for Palaeoptera. The first very large EST (expressed sequence tags = partial transcriptomes) dataset to evaluate arthropod relationships was published by Meusemann et al. in 2010 [349]. In addition to the size of the dataset, this paper was groundbreaking in terms of filtering the data. Randomized or phylogenetically uninformative sites were algorithmically identified and masked (box 4) with a program called Aliscore [240], and the matrix was optimized with the MARE program [350]. MARE eliminates both genes and taxa that are problematic owing to missing data, resulting in a smaller, but more dense matrix. Meusemann et al. recovered both Palaeoptera and Chiastomyaria using alternative analytical parameters (box 3). In addition, they also found that Hymenoptera was the sister taxon of other Holometabola, and a monophyletic Nonoculata, as in rRNA-dominated analyses.
A consistent pattern with large datasets is that they tend to recover either Palaeoptera [117], Chiastomyaria or both [50,349], but rarely the morphologically favoured Metapterygota (but see [351]). While Misof et al. [50] reported the monophyly of Palaeoptera, the quartet mapping (box 4) analyses reported in their supplementary materials favoured Chiastomyaria and rejected the morphologically favoured Metapterygota. The resolution of Palaeoptera was predicted to be among the most difficult nodes to recover [352], and even now, with millions of nucleotides applied to the question, it must be considered unresolved (figure 5). We continue to see that the problem nodes from morphology still exist, with continued conflict for the placement of Diplura and Zoraptera and the status of Palaeoptera.
18.2. Strepsiptera revisited
Among the first of the truly genomic analyses was Niehuis et al. in 2012 [175], who sequenced the Strepsiptera nuclear genome, and compared it with previously sequenced genomes of two beetles, four hymenopterans, three flies and Bombyx (silkmoth), with two outgroups. They tested four hypotheses regarding the placement of Strepsiptera, and concluded that it belonged with beetles, as originally placed by morphologists. This was further strengthened by McKenna [176], who added a neuropteran genome to the analysis, and by Boussau et al. [177], who added transcriptomes from additional key taxa to genomic sequences, and analysed them with models designed to avoid branch-length effects. They ruled out the possibility of a close relationship between Strepsiptera and either of the beetle families, Rhipiphoridae or Meloidae (among others).
18.3. Insect phylogeny resolved
Many laboratories are now collecting large datasets from hybrid capture techniques, such as anchored hybrid enrichment [353], or ultraconserved elements (UCEs) [354]. Both methods allow for the recovery of data from degraded museum specimens. Anchored hybrid enrichment has the advantage that probes can be designed from transcriptomes, making data from the two sources completely combinable. UCEs have the advantage that probes can be designed without the need for sequences from closely related taxa. The Weirauch (Heteroptera), Johnson/Dietrich (Hemiptera and Psocodea), McKenna (Coleoptera), Wiegmann (Diptera), Kawahara (Lepidoptera), Ward, Borowiec, Schultz and Brady (Formicidae) and Kjer/Frandsen (Trichoptera) laboratories currently have large hybrid-enriched datasets in progress, and, according to their conference presentations, these data are largely resolving long-standing problems with strong bootstrap support.
The Misof et al. insect phylogenomics study based on transcriptomes of 1478 genes [50], published in Science in November 2014, is by far the largest analysis to date of insect relationships and their phylogeny, and provides our current consensus (figure 5). Multiple technical advances occurred between Meusemann et al. in 2010 [349] and the Misof study in 2014 [50]. Both studies involved many of the same authors. Although the Misof study is short in print, one of their strongest contributions is the 200 pages of supplementary materials, which include recommendations for careful assembly (box 4), orthology prediction (box 4), data masking (box 4) and signal optimization. Protein domains were considered as partitions (box 4), and site-specific rate models were developed. Diplura was sister to Insecta, in agreement with Letsch & Simon [355], despite the recovery of Entognatha in other large datasets [347,349]. The 2014 Misof et al. study was the first of the publications from the 1KITE initiative, which has now collected transcriptomes from over 1400 taxa. Subprojects in the works from 1KITE include large datasets targeted at ‘basal hexapods’, Odonata, Polyneoptera, Paraneoptera, Hymenoptera, Coleoptera, Neuropterida, Trichoptera, Lepidoptera and Amphiesmenoptera, along with over 100 side projects dealing with molecular evolution in insects. These side projects put insect phylogenomics at the forefront of the discovery of character systems based on genomic meta-characters, and their impact reaches beyond entomology with development of new phylogenomic approaches. The Misof et al. study has already had a visible impact on phylogenetics of higher order groups by facilitating target enrichment and providing open access data.
19. Integrated phylogenetics
Innovative approaches such as µCT [356] and computer-based reconstruction [357], an optimized combined application of different techniques, and the concept of evolutionary morphology [358] have led to a remarkable renaissance in insect morphology in the last two decades, especially in Europe and Japan. Recent years have been characterized by matrices of increasing size, and a distinctly improved documentation of the characters made possible by the use of a broad array of techniques and an optimized workflow [359]. The Bayesian results of the largest morphological character state matrix used in insect systematics up to that time [48] were fully compatible with transcriptomic studies [49,50], indicating that morphology can still play a role in estimating and corroborating molecular phylogenetics.
Many examples of the value of integrated phylogenetics come from Misof et al. [50]. Perhaps their most unusual result was the grouping of Psocodea with Holometabola, but the possibility that model misspecification had influenced this placement could not be eliminated [50]. Morphological data provide another reason to be sceptical [152]. Zoraptera was not strongly placed in their study either [50], but it was reliably placed in a monophyletic Polyneoptera, as also suggested by recent morphological and embryological studies [79,360]. Despite recent progress, it is obvious that morphology has its limitations. Even large datasets of high quality create partially unsatisfying results, sometimes despite impressive lists of shared-derived character states (synapomorphies) for presumptive clades, as in [360]. Artefactual synapomorphies created by phylogenetic analyses—i.e. ‘cladistic noise’—often suggest results that are in fact insufficiently supported or not supported at all by any convincing features. An example is the ‘clade’ Dictyoptera + (Zoraptera + Plecoptera) supported by a recent parsimony analysis of morphological data by Matsumura et al. [360]. As the authors pointed out, some of the obtained presumptive synapomorphies were obviously the result of misleading redundant evolution. Correlated characters can also cause artefacts in morphology-based phylogenetic reconstructions [361] as addressed in the context of the Palaeoptera problem [362]. Solutions were suggested, based on modified weighting or the exclusion of characters.
19.1. The closest relatives of hexapods
A crucial question was apparently inaccessible to morphological approaches but was largely solved with molecular data: the systematic position of Hexapoda. A monophyletic Tracheata (Myriapoda + Hexapoda) was considered as granted in morphology-based studies, with Hexapoda placed either as the sistergroup of Myriapoda (millipedes and centipedes) or as the sister taxon of a myriapod subgroup [363,364]. Molecular datasets of different size and composition, and analysed with different approaches, consistently yielded a clade Pancrustacea (also called Tetraconata), usually with hexapods placed among paraphyletic crustacean lineages [50,257,346,347,349,366–368]. Even though some morphological arguments for Pancrustacea have been presented [369,370], a formal character analysis is still lacking and the morphological evidence is far from convincing. Possible candidates for the closest relatives of Hexapoda within Pancrustacea include the highly specialized relict group Remipedia, Malacostraca, and possibly the miniaturized Cephalocarida [50,347], although other studies contradict this [117,371]. The tremendous morphological gap between these aquatic groups and the terrestrial hexapods hinders meaningful comparisons of morphological characters and hypotheses of homology. It was pointed out by Klass & Kristensen [73] that the monophyly of Hexapoda is not strongly supported morphologically, with basically only one character complex defining it—the regional specialization of the body into head, thorax and abdomen, with the thorax divided into three-segments and the abdomen into 11. However, as shown in Beutel et al. [372], the Pancrustacea concept has strong implications for the hexapod groundplan. The strongly supported placement of Hexapoda among crustacean groups implies an entire series of additional hexapod autapomorphies, such as terrestrial habits, simplified walking legs, fusion of the second maxillae (labium), the loss of the ventral food rim, the absence of midgut glands and nephridial organs, and others.
20. Confidence and caution
Figure 5b shows an exponential growth in the size of datasets since the late 1980s, and we expect this growth to continue. It is tempting to think that every part of insect phylogeny has now been resolved with large datasets. Almost every node has strong bootstrap support. However, bootstrap support was designed to evaluate stochasticity and, with large datasets, stochasticity is reduced or even eliminated. This would be considered a good thing, if models and assumptions were perfect. However, because of the size of current datasets, small biases in the data, or misspecifications of the model can result in strong bootstrap support for error. We predict that model refinement will be a rich source of discovery in the future. For example, the failure to resolve the Palaeoptera problem, as indicated by quartet mapping (box 4) [50], may point towards a true case of conflicting gene tree histories. However, we are reluctant to assume this biological explanation without appropriate analysis. Such a comfortable explanation for misbehaved data, or inappropriate models, can make analytical failures seem like new discoveries. Every caution available at the time was applied by Misof et al. [50], and we find no reason to doubt their results. Moreover, their cautions are reflected in the expansive supplementary material [50]. We would like to emphasize that phylogenies can never be more than hypotheses, subject to the limitations of models and assumptions. This statement is obvious, but bears repeating in the light of enormous datasets that are now available.
With every advance, from cladistics to Sanger sequencing, and now genomics, we saw a wave of overconfidence that intractable conflicts would be solved, only to learn of new obstacles. We are only beginning to understand the behaviour of large datasets, which is why we cannot write about these innovations from a historical perspective. The discipline will continue to improve. We still look for confirmation and congruence from other sources of data, such as morphology and rare genomic events.
Both morphological and molecular investigations focused on insect systematics have made tremendous progress in the last decade. Similarly, the investigation of extinct insects has increased its pace with advanced morphological techniques allowing a stunningly detailed reconstruction of amber fossils. Although improvements could still be made in communication across different lines of investigation, it is unlikely that we will see many major revolutions in deep insect phylogeny, outside the groups, we have flagged as unresolved. The disagreement over parsimony versus likelihood has been resolved. With genomics, we will likely continue to learn about the function of genes and links between developmental and phylogenetic processes and how these processes change over time and across lineages. Optimized pipelines (box 4) of processing and connecting different sources of evidence is presently a key target for the future. Such pipelines are one of the main aims of 1KITE and associated projects. Continued integration of different disciplines will likely lead to a much better understanding of the complex evolution of insects, revealing why this group of organisms reached unparalleled species diversity and successfully conquered virtually every terrestrial and freshwater environment on the Earth.
taxon | taxa included | evidence | references |
---|---|---|---|
Acercaria | Thysanoptera, Hemiptera, Psocodea | cerci absent, 1 abdominal ganglionic mass, lacinia chisel-like, four Malpighian tubules | [9,11,14,31,39,62,72,74,77]mc[107,141,152,239,351,375] |
Aparaglossata | Holometabola, minus Hymenoptera | loss of paraglossae,ovipositor modified or reduced, max. eight Malpighian tubules | [47–50,77,117,154,170,175,177,239,342,343,349,351,355,375–377] |
Amphiesmenoptera | Lepidoptera, Trichoptera | many (see Kristensen) | generally accepted |
Cercophora | Diplura, Insecta | double claws, 9+9+2 axoneme, cerci | [14,50,72,355] |
Chiastomyaria | Ephemeroptera, Neoptera | indirect flight musculature, copulation with aedeagus | [31,107,117,239,281,355,378] |
Coleopterida | Coleoptera, Strepsiptera | posteromotorism | [9,11,14,31,47–50,62,72,154,170,175,177,343,375,376] |
Condylognatha | Thysanoptera, Hemiptera | mandibular stylet(s) | [9,11,31,50,62,308] |
Dicondylia | Zygentoma, Pterygota | secondary mandibular joint, gonanglum | [9,11,14,31,50,62,72,74]mc [107,117,141,170,276,343,349,347,375] |
Dictyoptera | Blattodea (incl. termites), Mantodea | secondary anterior tentorial bridge, female genital vestibulum, ootheca | generally accepted, although some put roaches with the mantids (based on shared plesiomorphies) |
Dictyoptera+Zoraptera | Dictyoptera, Zoraptera | [69] tentatively [74]mc [251,355,374] | |
Ectognatha (Insecta) | Archaeognatha, Dicondylia | ovipositor, flagellar antenna, Johnston's organ, corpotentorium, terminal filament | generally accepted |
Ellipura | Collembola, Protura | specific entognathy, linea ventralis | [10,11,62,63,72,74]mc [140,343,379–382] |
Entognatha | Collembola, Protura, Diplura | entognathy, eyes partly reduced | [11,31,62,107,117,141,242,248,320,347,349,374,375] |
Halteria | Strepsiptera, Diptera | rRNA (parsimony, NJ), substitution rate similarity, phenetic nucleotide compositional similarity | [74,77,159] |
Haplocerata | Plecoptera, Zoraptera | transcriptomes | [239] |
Hemiptera | Acercaria excl. Psocodea and Thysanoptera | four-segmented labial rostrum, labial endite lobes and palps absent, buccal pump | generally accepted |
Holometabola | Neuropteroidea Mecopterida Hymenoptera | complete Metamorphosis, characters of larvae | generally accepted |
Mecopterida | Amphiesmenoptera, Antliophora | ecdysone receptors, ovipositor absent, telescoping post-abdomen | [10,11,14,31,47–50,62,72,116,117,141,154,170,175,177,204,251,343,349,351,355,374–378] |
Metapterygota | Odonata, Neoptera | secondary mandibular articulation as ball-and-socket joint, modified mandibular muscles, no subimago | [14,62,72,74]mc[77,114,204,351] |
Neoptera | Pterygota, minus Palaeoptera | wings folded back over abdomen, modified wing base, arolium (?) | generally accepted |
Neuropterida | Raphidioptera, Megaloptera, Neuroptera | third valve of ovipositor with intrinsic muscles | generally accepted |
Neuropteroidea | Neuropterida, Coleopterida | modifications of ovipositor (?) | [36,47,49,50,62,72,141,154,159,170,204,343,375] |
Nonoculata | Protura, Diplura | rRNA, EF-1α introns, eyes lost | [107,141,239,242–245,248,249,349,351,375,376,383–385] |
Palaeoptera | Odonata, Ephemeroptera | bristle-like antennae, aquatic larvae | [9,11,50,101,141,159,170,343,347,349,386] |
Pancrustacea | ‘Crustacea’, Hexapoda | rRNA, mitochondrial gene order, four-partite crystalline cone | [50,76,107,108,110,117,141,232,242,244,245,256,276,343,349,347,384,385,368,387–389] |
Paraneoptera | Acercaria, Zoraptera | 6/4 Malpighian tubules, 3/2 tarsomeres, 2/1 abdominal ganglionic complexes | [11,31,62,152,376] |
Paurometabola | Polyneoptera excl. Plecoptera | fan-like folding of hind wing, enlarged euplantulae, terrestrial larvae | [10,375] |
Polyneoptera | Neoptera excl. Acercaria and Holometabola | tegmina, enlarged anal field of hind wing, euplantulae | [14,31,39,50,343,347,355] |
Psocodea | Psocoptera, Phthiraptera | cibarial water-vapour uptake apparatus, antennal rupture facilitating device | [11,14,31,50,72,74]mc[77,107,170,177,239,343,375,390] |
Pterygota | Winged insects (including secondarily wingless orders) | wings, copulation (?) | generally accepted |
Thysanura, s.l. | Archaeognatha, Zygentoma | phenetic similarity (symplesiomorphies) | [9,375] |
Thysanura, s.s. | Zygentoma | sperm coupling, loss of superlinguae | [50,239,373] |
Tracheata | Myriapoda, Hexapoda | tracheal system, Malpighian tubules, spermatophore, loss of second antenna etc. | [140,380] generally accepted before 1990. |
Xenonomia | Grylloblattodea, Mantophasmatodea | rRNA, transcriptomes | [50,100,141,152,204] |
Authors' contributions
K.M.K. and R.G.B. wrote about their respective areas of expertise. C.S. provided comments, fact-checking and additional historical and methodological insights. M.Y. interpreted and summarized Russian works.
Competing interests
We have no competing interests.
Funding
We received no funding for this study.
Acknowledgements
We thank Nicole Tam for preparing the illustrations. Phil Ward, Marek Borowiec, Harald Letsch, Charles Mitter, John Huelsenbeck, Bjorn v. Reumont, Duane McKenna, Günther Pass, Alexander Blanke, Bernhard Misof, and Sabrina Simon provided helpful comments. John Morse provided valuable insights about his advisor, Herbert Ross. K.M.K. thanks the Schlinger endowment for funding.
References
- 1
Beutel RG, Friedrich F, Leschen RA . 2009 Charles Darwin, beetles and phylogenetics. Naturwissenschaften 96, 1293–1312. (doi:10.1007/s00114-009-0601-2) Crossref, PubMed, ISI, Google Scholar - 2
Engel MS, Kristensen NP . 2013 A history of entomological classification. Annu. Rev. Entomol. 58, 585–607. (doi:10.1146/annurev-ento-120811-153536) Crossref, PubMed, ISI, Google Scholar - 3
Linneaus C . 1758 Systema naturæ per regna tria naturæ, secundum classes, ordines, genera, species, cum characteribus, differentiis, synonymis, locis, 10th edn. Stockholm: Laurentius Salvius. Crossref, Google Scholar - 4
- 5
Latreille PA . 1796 Précis des caractères génériques des insectes, disposés dans un ordre naturel. [Internet]. Paris, France: Prèvôt. http://www.biodiversitylibrary.org/bibliography/58411 (accessed 10 Jan 2016). Google Scholar - 6
Darwin CB, Wallace AR . 1858 On the tendency of species to form varieties; and on the perpetuation of varieties and species by natural means of selection. J. Proc. Linn. Soc. Lond. Zool. 3, 45–50. (doi:10.1111/j.1096-3642.1858.tb02500.x) Crossref, Google Scholar - 7
Darwin C . 1859 On the origin of species by means of natural selection, or the preservation of favoured races in the struggle for life. London, UK: John Murray. Crossref, Google Scholar - 8
Haekel E . 1896 Systematische Phylogenie der wirbellosen Thiere. Berlin, Germany: Reimer. Crossref, Google Scholar - 9
- 10
Hennig W . 1953 Kritische Bemerkungen zum phylogenetischen System der Insecten. Beitr. Z. Entomol. Sonderh. 3, 1–85. Google Scholar - 11
Hennig W . 1969 Die Stammesgeschichte der Insekten. Frankfurt am Main: Waldemar Kramer. Google Scholar - 12
Crampton GC . 1918 A phylogenetic study of the terga and wing bases in Embiids, Plecoptera, Dermaptera, and Coleoptera. Psyche (Stuttg). 25, 4–13. (doi:10.1155/1918/96037) Crossref, Google Scholar - 13
Crampton GC . 1928 The eulabium, mentum, submentum and gular region of insects. J. Entomol. Zool. 20, 1–18. Google Scholar - 14
Crampton GC . 1938 The interrelationships and lines of descent of living insects. Psyche (Stuttg). 45, 165–181. (doi:10.1155/1938/18939) Crossref, Google Scholar - 15
- 16
- 17
- 18
- 19
Friedrich F, Matsumura Y, Pohl H, Bai M, Hörnschemeyer T, Beutel RG . 2014 Insect morphology in the age of phylogenomics: innovative techniques and its future role in systematics. Entomol. Sci. 17, 1–24. (doi:10.1111/ens.12053) Crossref, ISI, Google Scholar - 20
Handlirsch A . 1937 Neue Untersuchungen über die fossilen Insekten. Ann. Naturhistorischen. Mus. Wien. 40, 1–40. Google Scholar - 21
Handlirsch A . 1939 ,Neue Untersuchungen über die fossilen Insekten. Ann. Naturhistorischen. Mus. Wien. 49, 1–240. Google Scholar - 22
Hennig W . 1950 Grundzüge einer Theorie der Phylogenetischen Systematik. Berlin: Deutscher Zentralverlag. Google Scholar - 23
Watrous LE, Wheeler QD . 1981 The out-group comparison method of character analysis. Syst. Biol. 30, 1–11. (doi:10.1093/sysbio/30.1.1) Crossref, Google Scholar - 24
Maddison WP, Donoghue MJ, Maddison DR . 1984 Outgroup analysis and parsimony. Syst. Zool. 33, 83–103. (doi:10.2307/2413134) Crossref, Google Scholar - 25
Sturtevant AV . 1939 On the subdivision of the genus Drosophila. Proc. Natl Acad. Sci. USA 25, 137–141. (doi:10.1073/pnas.25.3.137) Crossref, PubMed, Google Scholar - 26
Ross HH . 1937 A generic classification of the Nearctic sawflies (Hymenoptera: Symphyta). Ill. Biol. Monogr. 15, 1–173. (doi:10.5962/bhl.title.50339) Google Scholar - 27
- 28
- 29
Mayr E . 1975 Grundlagen der Zoologischen Systematik. Hamburg, Berlin: Verlaug Paul Parey. Google Scholar - 30
Nelson G . 1974 Darwin-Hennig classification: a reply to Ernst Mayr. Syst. Zool. 23, 452–458. (doi:10.2307/2412551) Crossref, Google Scholar - 31
Wille A . 1960 The phylogeny and relatonships between the insect orders. Rev. Biol. Trop. 8, 93–123. Google Scholar - 32
Hinton HE . 1958 The phylogeny of the panorpoid orders. Annu. Rev. Entomol. 3, 181–206. (doi:10.1146/annurev.en.03.010158.001145) Crossref, ISI, Google Scholar - 33
- 34
Mickoleit G . 1973 Über den Ovipositor der Neuropteroidea und Coleoptera und seine phylogenetische Bedeutung (Insecta, Holometabola). Z. Morphol. Tiere. 74, 37–64. (doi:10.1007/BF00291795) Crossref, Google Scholar - 35
Mickoleit G . 1975 Die Genital- und Postgenitalsegmente der Mecoptera-Weibchen (Insecta, Holometabola). I. Das Exoskelet. Z. Morphol. Tiere. 80, 97–135. (doi:10.1007/BF00281741) Crossref, Google Scholar - 36
Mickoleit G . 1976 Die Genital- und Postgenitalsegmente der Mecoptera-Weibchen (Insecta, Holometabola). II. Das Dach der Genitalkammer. Z. Morphol. Tiere. 88, 133–156. Google Scholar - 37
Mickoleit G . 1978 Die phylogenetischen Beziehungen der Schnabelfliegen-Familien aufgrund morphologischer Ausprägungen der weiblichen Genital- und Postgenitalsegmente (Mecoptera). Entomol. Ger. 4, 258–271. Google Scholar - 38
Mickoleit G . 2008 Die Sperma-Auspreßvorrichtung der Nannochoristidae (Insecta: Mecoptera). Entomol. Gen. 31, 193–226. (doi:10.1127/entom.gen/31/2008/193) Crossref, Google Scholar - 39
Martynov AV . 1925 Über zwei Grundtypen der Flügel bei den Insecten und ihre Evolution. Zoomorphology 4, 465–501. Google Scholar - 40
Martynov AV . 1938 Studies on the geological history and phylogeny of the orders of insects (Pterygota). Tr. Paleont. Inst. Inst. Akad. Nauk SSSR 7, 1–150. Google Scholar - 41
Ponomarenko AG . 1969 The historical development of archostematan beetles. Tr. Paleontol. Inst. Akad. Nauk SSSR 125, 1–238. Google Scholar - 42
Ponomarenko AG . 1977 Adephaga. In Mesozoic Coleoptera, Suborder, etc. (edsArnoldy LV, Jerikin VV, Nikritin LM, Ponomarenko AG ). Tr. Palaeont Inst Akad Nauk SSSR 161, 17–104. Google Scholar - 43
Rasnitsyn AP, Quicke DLJ (eds). 2002 History of insects. Dordrecht: Kluwer Academic Publishers. Crossref, Google Scholar - 44
Hoßfeld U, Olsson L . 2002 From the modern synthesis to Lysenkoism, and back? Science 297, 55. (doi:10.1126/science.1068355) Crossref, PubMed, ISI, Google Scholar - 45
- 46
Jacobson GG . 1905 The beetles of Russia, West Europe and adjacent countries. St Petersburg, Russia: Devrien. Google Scholar - 47
Wiegmann BM, Trautwein MD, Kim J-W, Cassel BK, Bertone MA, Winterton SL, Yeates DK . 2009 Single-copy nuclear genes resolve the phylogeny of the holometabolous insects. BMC Biol. 7, 34. (doi:10.1186/1741-7007-7-34) Crossref, PubMed, ISI, Google Scholar - 48
Beutel RG 2011 Morphological and molecular evidence converge upon a robust phylogeny of the megadiverse Holometabola. Cladistics 27, 341–355. (doi:10.1111/j.1096-0031.2010.00338.x) Crossref, ISI, Google Scholar - 49
Peters RS 2014 The evolutionary history of holometabolous insects inferred from transcriptome-based phylogeny and comprehensive morphological data. BMC Evol. Biol. 14, 52. (doi:10.1186/1471-2148-14-52) Crossref, PubMed, ISI, Google Scholar - 50
Misof B 2014 Phylogenomics resolves the timing and pattern of insect evolution. Science 346, 763–767. (doi:10.1126/science.1257570) Crossref, PubMed, ISI, Google Scholar - 51
Ronquist F, Rasnitsyn AP, Roy A, Eriksson K, Lindgren M . 1999 Phylogeny of the Hymenoptera: a cladistic reanalysis of Rasnitsyn's (1988) data. Zool. Scr. 28, 13–50. (doi:10.1046/j.1463-6409.1999.00023.x) Crossref, ISI, Google Scholar - 52
Rasnitsyn AP . 2003 On the skimming hypothesis of the origin of insect flight. Acta Zool. Cracoviensia. 46, 85–88. Google Scholar - 53
Brothers DJ . 2011 Alexandr Pavlovich Rasnitsyn, (palaeo)entomologist extraordinaire—a personal appreciation. ZooKeys 130, 1–10. (doi:10.3897/zookeys.130.1890) Crossref, ISI, Google Scholar - 54
Jeannel R . 1949 Ordre des Coleopteroides. In Traite de Zoologie (ed.Grasse P-P ), pp. 771–891. Paris, France: Masson et Cie. Google Scholar - 55
Hennig W . 1966 Phylogenetic systematics. Urbana, Illinois, USA: University of Illinois Press. Google Scholar - 56
Klass K-D (ed.) . 2003 Proceedings of the 1st Dresden meeting on insect phylogeny: ‘Phylogenetic Relationships within the Insect Orders’ (Dresden, September 19–21, 2003). Entomol Abh. 61, 119–72. Google Scholar - 57
- 58
Kristensen NP . 1968 The skeletal anatomy of the heads of adult Mnesarchaeidae and Neopseustidae (Lep., Dacnonypha). Ent. Meddr. 36, 137–151. Google Scholar - 59
Kristensen NP . 1968 The anatomy of the head and the alimentary canal of adult Eriocraniidae (Lep., Dacnonypha), Entomol. Meddelelser. 36, 239–315. Google Scholar - 60
Kristensen NP . 1968 The morphological and functional evolution of the mouthparts in adult Lepidoptera (abstract). In Opuscula Entomologica, pp. 69–72. Lund: Entomologiska Sällskapet. Google Scholar - 61
Kristensen NP . 1984 Studies on the morphology and systematics of primitive Lepidoptera. Steenstrupia 10, 141–191. Google Scholar - 62
Kristensen NP . 1975 The phylogeny of hexapod ‘orders’. A critical review of recent accounts. Z. Zool. Syst. Evol.-Forsch. 13, 1–44. (doi:10.1111/j.1439-0469.1975.tb00226.x) Crossref, Google Scholar - 63
Kristensen NP . 1981 Phylogeny of insect orders. Annu. Rev. Entomol. 26, 135–157. (doi:10.1146/annurev.en.26.010181.001031) Crossref, ISI, Google Scholar - 64
Kristensen NP . 1991 Phylogeny of extant hexapods. In The insects of Australia, pp. 126–140. Melbourne: Melbourne University Publishing. Google Scholar - 65
Kristensen NP . 1999 The non-glossatan moths. Lepidoptera: moths and butterflies. 1. Evolution, systematics, and biogeography. In Handbook of zoology, vol. IV, part 35, pp. 41–49. Berlin and New York: De Gruyter. Google Scholar - 66
Kristensen NP . 1999 The homoneurous glossata. Lepidoptera: moths and butterflies. 1. Evolution, systematics, and biogeography. In Handbook of zoology, vol. IV, part 35, pp. 51–63. Berlin: De Gruyter. Google Scholar - 67
Kristensen NP, Skalski AW . 1999 Phylogeny and paleontology. Lepidoptera: moths and butterflies. 1. Evolution, systematics, and biogeography. In Handbook of zoology, vol. IV, part 35, pp. 7–25. Berlin and New York: De Gruyter. Google Scholar - 68
Kristensen NP . 2003 Reproductive organs. Lepidoptera: moths and butterflies. 1. Evolution, systematics, and biogeography. In Handbook of zoology, vol. IV, part 36, pp. 427–447. Berlin: Walter de Gruyter. Google Scholar - 69
Boudreaux HB . 1979 Arthropod phylogeny with special reference to the insects. New York, Chichester, Brisbane, Toronto: Wiley & Sons Inc. Google Scholar - 70
Mashimo Y 2014 100 years Zoraptera—a phantom in insect evolution and the history of its investigation. Insect Syst. Evol. J. Morphol. 45, 371–393. (doi:10.1163/1876312X-45012110) Crossref, ISI, Google Scholar - 71
Wipfler B, Bai M, Schoville S, Dallai R, Uchifune T, Machida R, Cui Y, Beutel RG . 2014 Ice crawlers (Grylloblattodea)–the history of the investigation of a highly unusual group of insects. J. Insect. Biodivers. 2, 1–25. (doi:10.12976/jib/2014.2.2) Crossref, Google Scholar - 72
Beutel RG, Gorb SN . 2001 Ultrastructure of attachment specializations of hexapods (Arthropoda): evolutionary patterns inferred from a revised ordinal phylogeny. J. Zool. Syst. Evol. Res. 39, 177–207. (doi:10.1046/j.1439-0469.2001.00155.x) Crossref, ISI, Google Scholar - 73
Klass K-D, Kristensen NP . 2001 The ground plan and affinities of hexapods: recent progress and open problems. Ann. Soc. Entomol. Fr. 37, 265–298. ISI, Google Scholar - 74
Wheeler WC, Whiting M, Wheeler QD, Carpenter JM . 2001 The phylogeny of the extant hexapod orders. Cladistics 17, 113–169. (doi:10.1111/j.1096-0031.2001.tb00115.x) Crossref, ISI, Google Scholar - 75
Yoshizawa K, Johnson KP . 2005 Aligned 18S for Zoraptera (Insecta): phylogenetic position and molecular evolution. Mol. Phylogenet. Evol. 37, 572–580. (doi:10.1016/j.ympev.2005.05.008) Crossref, PubMed, ISI, Google Scholar - 76
Wheeler WC, Cartwright P, Hayashi CY . 1993 Arthropod phylogeny: a combined approach. Cladistics 9, 1–39. (doi:10.1111/j.1096-0031.1993.tb00207.x) Crossref, ISI, Google Scholar - 77
Whiting MF, Carpenter JC, Wheeler QD, Wheeler WC . 1997 The strepsiptera problem: phylogeny of the holometabolous insect orders inferred from 18S and 28S ribosomal DNA sequences and morphology. Syst. Biol. 46, 1–68. (doi:10.1093/sysbio/46.1.1) PubMed, ISI, Google Scholar - 78
Kristensen NP . 1995 Forty years’ insect phylogenetic systematics. Zool. Beitr. N.F. 36, 83–124. Google Scholar - 79
Mashimo Y, Beutel RG, Dallai R, Lee C-Y, Machida R . 2014 Embryonic development of Zoraptera with special reference to external morphology, and its phylogenetic implications (Insecta). J. Morphol. 275, 295–312. (doi:10.1002/jmor.20215) Crossref, PubMed, ISI, Google Scholar - 80
Kukalova-Peck J . 1991 Fossil history and the evolution of hexapod structures. In The insects of Australia: a textbook for students and research workers (ed. CSIRO), pp. 141–179. Ithaca, NY: Cornell University Press. Google Scholar - 81
Kukalová-Peck J, Brauckmann C . 1990 Wing folding in pterygote insects, and the oldest Diaphanopterodea from the early Late Carboniferous of West Germany. Can. J. Zool. 68, 1104–1111. (doi:10.1139/z90-163) Crossref, ISI, Google Scholar - 82
Kukalova-Peck J, Peck SB . 1993 Zoraptera wing structures: evidence for new genera and relationship with the blattoid orders (Insecta: Blattoneoptera). Syst. Entomol. 18, 333–350. (doi:10.1111/j.1365-3113.1993.tb00670.x) Crossref, ISI, Google Scholar - 83
Kukalová-Peck J, Beutel RG . 2013 Is the carboniferous Adiphlebia lacoana really the ‘oldest beetle’. Eur. J. Entomol. 109, 633–645. (doi:10.14411/eje.2012.075) Crossref, ISI, Google Scholar - 84
Haas F, Kukalova-Peck J . 2001 Dermaptera hindwing structure and folding: new evidence for familial, ordinal and superordinal relationships within Neoptera (Insecta). Eur. J. Entomol. 98, 445–509. (doi:10.14411/eje.2001.065) Crossref, ISI, Google Scholar - 85
Kukalova-Peck J, Lawrence JF . 2004 Relationships among coleopteran suborders and major endoneopteran lineages: evidence from hind wing characters. Eur. J. Entomol. 101, 95–144. (doi:10.14411/eje.2004.018) Crossref, ISI, Google Scholar - 86
Kukalová-Peck J . 1983 Origin of the insect wing and wing articulation from the arthropodan leg. Can. J. Zool. 61, 1618–1669. (doi:10.1139/z83-217) Crossref, ISI, Google Scholar - 87
Kukalová-Peck J . 1987 New carboniferous Diplura, Monura, and Thysanura, the hexapod ground plan, and the role of thoracic side lobes in the origin of wings (Insecta). Can. J. Zool. 65, 2327–2345. (doi:10.1139/z87-352) Crossref, ISI, Google Scholar - 88
Kukalová-Peck J . 2008 Phylogeny of higher taxa in insecta: finding synapomorphies in the extant fauna and separating them from homoplasies. Evol. Biol. 35, 4–51. (doi:10.1007/s11692-007-9013-4) Crossref, ISI, Google Scholar - 89
Béthoux O, Kristensen NP, Engel MS . 2008 Hennigian phylogenetic systematics and the ‘groundplan’ vs. ‘post-groundplan’ approaches: a reply to Kukalová-Peck. Evol. Biol. 35, 317–323. (doi:10.1007/s11692-008-9035-6) Crossref, ISI, Google Scholar - 90
Pass G . 1985 Gross and fine structure of the antennal circulatory organ in cockroaches (Blattodea, Insecta). J. Morph. 185, 255–268. (doi:10.1002/jmor.1051850210) Crossref, PubMed, ISI, Google Scholar - 91
Pass G . 2000 Accessory pulsatile organs: evolutionary innovations in insects. Annu. Rev. Entomol. 45, 495–518. (doi:10.1146/annurev.ento.45.1.495) Crossref, PubMed, ISI, Google Scholar - 92
Klass K-D . 1995 Die Phylogeny der Dictyoptera. Ph.D. thesis. Fakultät fur Biologie, Ludwig Maximilians Universität, Munich, Germany. Google Scholar - 93
Klass K-D. 2003 The female genitalic region in basal earwigs (Insecta: Dermaptera; Pygidicranidae sl). Entomol. Abh. 61, 173–225. Google Scholar - 94
Machida R, Nagashima T, Ando H . 1990 The early embryonic development of the jumping bristletail Pedetontus unimaculatus machida (Hexapoda: Microcoryphia, Machilidae). J. Morphol. 206, 181–195. (doi:10.1002/jmor.1052060205) Crossref, PubMed, ISI, Google Scholar - 95
Machida R, Tojo T, Tsutsumi T, Uchifune T, Klass K-D, Picker MD 2004 Embryonic development of heel-walkers: reference to some prerevolutionary stages (Insecta: Mantophasmatodea). Proc. Arthropod Embryol. Soc. Jpn. 39, 31–39. Google Scholar - 96
Klass K-D . 1998 The proventriculus of the Dicondylia, with comments on evolution and phylogeny of Dictyoptera and Odonata (Insecta). Zool. Anz. 237, 15–42. ISI, Google Scholar - 97
Klass K-D . 1998 Possible homologies in the proventriculi of Dicondylia (Hexapoda) and Malacostraca (Crustacea). Zool. Anz. 237, 43–58. ISI, Google Scholar - 98
Klass K-D . 1998 The ovipositor of Dictyoptera (Insecta): homology and ground-plan of the main elements. Zool. Anz. 236, 69–101. ISI, Google Scholar - 99
Mickoleit G . 1961 Zur Thoraxmorphologie der Thysanoptera. Zool. Jahrb. Abt. Anat. Ontog. Tiere. 79, 1–92. Google Scholar - 100
Wipfler B, Machida R, Müller B, Beutel RG . 2011 On the head morphology of Grylloblattodea (Insecta) and the systematic position of the order, with a new nomenclature for the head muscles of Dicondylia. Syst. Entomol. 36, 241–266. (doi:10.1111/j.1365-3113.2010.00556.x) Crossref, ISI, Google Scholar - 101
Blanke A, Wipfler B, Letsch H, Koch M, Beckmann F, Beutel R, Misof B . 2012 Revival of Palaeoptera—head characters support a monophyletic origin of Odonata and Ephemeroptera (Insecta). Cladistics 28, 560–581. (doi:10.1111/j.1096-0031.2012.00405.x) Crossref, PubMed, ISI, Google Scholar - 102
Patterson C . 1978 Verifiability in systematics. Syst. Zool. 27, 218–222. (doi:10.2307/2412977) Crossref, Google Scholar - 103
- 104
Forey PL, Humphries CJ, Kitching IJ, Scotland RW, Siebert DJ, Williams DM . 1992 Cladistics: a practical course in systematics. The Systematics Association, publication no. 10. Oxford, UK: Clarendon Press. Google Scholar - 105
Farris JS . 1989 Hennig86, a PC-DOS program for phylogenetic analysis. Cladistics 5, 163. (doi:10.1111/j.1096-0031.1989.tb00573.x) Google Scholar - 106
Swofford DL . 1991 PAUP: phylogenetic analysis using parsimony, version 3.1. Champaign, IL: Illinois Natural History Survey. Google Scholar - 107
Kjer KM . 2004 Aligned 18S and insect phylogeny. Syst. Biol. 53, 506–514. (doi:10.1080/10635150490445922) Crossref, PubMed, ISI, Google Scholar - 108
Field KG, Olsen GJ, Lane DJ, Giovannoni SJ, Ghiselin MT, Raff E, Pace N, Raff R . 1988 Molecular phylogeny of the animal kingdom. Science 239, 748–753. (doi:10.1126/science.3277277) Crossref, PubMed, ISI, Google Scholar - 109
Lake JA . 1990 Origin of the Metazoa. Proc. Natl Acad. Sci. USA 87, 763–766. (doi:10.1073/pnas.87.2.763) Crossref, PubMed, ISI, Google Scholar - 110
Turbeville JM, Pfeifer DM, Field KG, Raff RA . 1991 The phylogenetic status of arthropods, as inferred from 18S rRNA sequences. Mol. Biol. Evol. 8, 669–686. PubMed, ISI, Google Scholar - 111
Berlocher SH . 1984 Insect molecular systematics. Annu. Rev. Entomol. 29, 403–433. (doi:10.1146/annurev.en.29.010184.002155) Crossref, ISI, Google Scholar - 112
Fujiwara H, Ishikawa H . 1982 Primary and secondary structures of Tetrahymena and aphid 5.8S rRNAs: structural features of 5.8S rRNA which interacts with the 28S rRNA containing the hiddenbreak. Nucleic Acids Res. 10, 5173–5182. (doi:10.1093/nar/10.17.5173) Crossref, PubMed, ISI, Google Scholar - 113
Simon C . 1988 Evolution of 13- and 17-year periodical Cicadas (Homoptera: Cicadidae: Magicicada). Bull. Entomol. Soc. Am. 34, 163–176. (doi:10.1093/besa/34.4.163) Google Scholar - 114
Wheeler WC . 1989 The systematics of insect ribosomal DNA. In The hierarchy of life molecules and morphology in phylogenetic analysis (edsFernholm KBB, Jörnvall H ), pp. 307–321. Amsterdam, The Netherlands: Elsevier. Google Scholar - 115
Carmean D, Kimsey LS, Berbee ML . 1992 18S rDNA sequences and the holometabolous insects. Mol. Phylogenet. Evol. 1, 270–278. (doi:10.1016/1055-7903(92)90002-X) Crossref, PubMed, ISI, Google Scholar - 116
Pashley DP, McPheron BA, Zimmer EA . 1993 Systematics of the holometabolous insect orders based on 18S ribosomal RNA. Mol. Phylogenet. Evol. 2, 132–142. (doi:10.1006/mpev.1993.1013) Crossref, PubMed, ISI, Google Scholar - 117
von Reumont BM 2009 Can comprehensive background knowledge be incorporated into substitution models to improve phylogenetic analyses? A case study on major arthropod relationships. BMC Evol. Biol. 9, 119. (doi:10.1186/1471-2148-9-119) Crossref, PubMed, ISI, Google Scholar - 118
Liu H, Beckenbach AT . 1992 Evolution of the mitochondrial cytochrome oxidase II gene among 10 orders of insects. Mol. Phylogenet. Evol. 1, 41–52. (doi:10.1016/1055-7903(92)90034-E) Crossref, PubMed, ISI, Google Scholar - 119
Saitou N, Nei M . 1987 The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406–425. PubMed, ISI, Google Scholar - 120
Ballard JW, Olsen GJ, Faith DP, Odgers WA, Rowell DM, Atkinson PW . 1992 Evidence from 12S ribosomal RNA sequences that onychophorans are modified arthropods. Science 258, 1345–1348. (doi:10.1126/science.1455227) Crossref, PubMed, ISI, Google Scholar - 121
Kjer KM . 1995 Use of rRNA secondary structure in phylogenetic studies to identify homologous positions: an example of alignment and data presentation from the frogs. Mol. Phylogenet. Evol. 4, 314–330. (doi:10.1006/mpev.1995.1028) Crossref, PubMed, ISI, Google Scholar - 122
Hickson RE, Simon C, Cooper A, Spicer GS, Sullivan J, Penny D . 1996 Conserved sequence motifs, alignment, and secondary structure for the third domain of animal 12S rRNA. Mol. Biol. Evol. 13, 150–169. (doi:10.1093/oxfordjournals.molbev.a025552) Crossref, PubMed, ISI, Google Scholar - 123
Simon C, Frati F, Beckenbach AT, Crespi B, Liu H, Flook PK . 1994 Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann. Entomol. Soc. Am. 87, 651–701. (doi:10.1093/aesa/87.6.651) Crossref, ISI, Google Scholar - 124
Kocher TD, Thomas WK, Meyer A, Edwards SV, Pääbo S, Villablanca FX, Wilson AC . 1989 Dynamics of mitochondrial DNA evolution in animals: amplification and sequencing with conserved primers. Proc. Natl Acad. Sci. USA 86, 6196–6200. (doi:10.1073/pnas.86.16.6196) Crossref, PubMed, ISI, Google Scholar - 125
Simon C, Pääbo S, Kocher TD, Wilson AC . 1990 Evolution of the mitochondrial ribosomal RNA in insects as shown by the poiymerase chain reaction. In Molecular evolution. UCLA Symposia on Molecular and Cellular Bioiogy, New Series (edsClegg M, O'Brien S ), pp. 235–244. New York, NY: Liss. Google Scholar - 126
Caterino MS, Cho S, Sperling FA . 2000 The current state of insect molecular systematics: a thriving Tower of Babel. Annu. Rev. Entomol. 45, 1–54. (doi:10.1146/annurev.ento.45.1.1) Crossref, PubMed, ISI, Google Scholar - 127
Friedlander TP, Regier JC, Mitter C . 1992 Nuclear gene sequences for higher level phylogenetic analysis: 14 promising candidates. Syst. Biol. 41, 483–490. (doi:10.1093/sysbio/41.4.483) Crossref, ISI, Google Scholar - 128
Friedlander TP, Regier JC, Mitter C . 1994 Phylogenetic information content of five nuclear gene sequences in animals: initial assessment of character sets from concordance and divergence studies. Syst. Biol. 43, 511–525. (doi:10.1093/sysbio/43.4.511) Crossref, ISI, Google Scholar - 129
Cho S, Mitchell A, Regier JC, Mitter C, Poole RW, Friedlander TP 1995 A highly conserved nuclear gene for low-level phylogenetics: elongation factor-1 alpha recovers morphology-based tree for heliothine moths. Mol. Biol. Evol. 12, 650–656. PubMed, ISI, Google Scholar - 130
Friedlander TP, Regier JC, Mitter C, Wagner DL . 1996 A nuclear gene for higher level phylogenetics: phosphoenolpyruvate carboxykinase tracks mesozoic-age divergences within Lepidoptera (Insecta). Mol. Biol. Evol. 13, 594–604. (doi:10.1093/oxfordjournals.molbev.a025619) Crossref, PubMed, ISI, Google Scholar - 131
Friedlander TP, Regier JC, Mitter C . 1997 Initial assessment of character sets from five nuclear gene sequences in animals. In Biodiversity II (edsReaka-Kudla ML, Wilson DE, Wilson EO ), pp. 301–320. Washington, DC: Joseph Henry Press. Google Scholar - 132
Friedlander TP, Regier JC, Mitter C, Wagner DL, Fang QQ . 2000 Evolution of heteroneuran Lepidoptera (Insecta) and the utility of dopa decarboxylase for Cretaceous-age phylogenetics. Zool. J. Linn. Soc. 130, 213–234. (doi:10.1111/j.1096-3642.2000.tb01630.x) Crossref, ISI, Google Scholar - 133
Wiegmann BM, Regier JC, Mitter C, Friedlander TP, Wagner DL, Nielson ES . 2000 Nuclear genes resolve Mesozoic-aged divergences in the insect order Lepidoptera. Mol. Phylogenet. Evol. 15, 242–259. (doi:10.1006/mpev.1999.0746) Crossref, PubMed, ISI, Google Scholar - 134
Fang QQ, Cho S, Regier JC, Mitter C, Matthews M, Poole RW, Friedlander TP, Zhao S . 1997 A new nuclear gene for insect phylogenetics: dopa decarboxylase is informative of relationships within Heliothinae (Lepidoptera: Noctuidae). Syst. Biol. 46, 269–283. (doi:10.1093/sysbio/46.2.269) Crossref, PubMed, ISI, Google Scholar - 135
Fang QQ, Mitchell A, Regier JC, Mitter C, Friedlander TP, Poole RW . 2000 Phylogenetic utility of the nuclear gene dopa decarboxylase in noctuoid moths (Insecta: Lepidoptera: noctuoidea). Mol. Phylogenet. Evol. 15, 473–486. (doi:10.1006/mpev.1999.0758) Crossref, PubMed, ISI, Google Scholar - 136
Burmester T, Massey HC, Zakharkin SO, Benes H . 1998 The evolution of hexamerins and the phylogeny of insects. J. Mol. Evol. 47, 93–108. (doi:10.1007/PL00006366) Crossref, PubMed, ISI, Google Scholar - 137
Campbell DL, Brower AV, Pierce NE . 2000 Molecular evolution of the wingless gene and its implications for the phylogenetic placement of the butterfly family Riodinidae (Lepidoptera: papilionoidea). Mol. Biol. Evol. 17, 684–696. (doi:10.1093/oxfordjournals.molbev.a026347) Crossref, PubMed, ISI, Google Scholar - 138
Ascher JS, Danforth BN, Ji S . 2001 Phylogenetic utility of the major opsin in bees (Hymenoptera: Apoidea): a reassessment. Mol. Phylogenet. Evol. 19, 76–93. (doi:10.1006/mpev.2001.0911) Crossref, PubMed, ISI, Google Scholar - 139
Colgan DJ, McLauchlan A, Wilson GDF, Livingston SP, Edgecombe GD, Macaranas J, Cassis G, Gray MR . 1998 Histone H3 and U2 snRNA DNA sequences and arthropod molecular evolution. Aust. J. Zool. 46, 419. (doi:10.1071/ZO98048) Crossref, ISI, Google Scholar - 140
Edgecombe GD, Wilson GDF, Colgan DJ, Gray MR, Cassis G . 2000 Arthropod cladistics: combined analysis of histone H3 and U2 snRNA sequences and morphology. Cladistics 16, 155–203. (doi:10.1111/j.1096-0031.2000.tb00352.x) Crossref, ISI, Google Scholar - 141
Kjer KM, Carle FL, Litman J, Ware J . 2006 A molecular phylogeny of Hexapoda. Arthropod Syst. Phylogeny 64, 35–44. Google Scholar - 142
Regier JC, Shultz JW, Kambic RE . 2004 Phylogeny of basal hexapod lineages and estimates of divergence times. Ann. Entomol. Soc. Am. 97, 411–419. (doi:10.1603/0013-8746(2004)097[0411:POBHLA]2.0.CO;2) Crossref, ISI, Google Scholar - 143
Saiki RK, Scharf S, Faloona F, Mullis KB, Horn GT, Erlich HA, Arnheim N . 1985 Enzymatic amplification of beta-globin genomic sequences and restriction site analysis for diagnosis of sickle cell anemia. Science 230, 1350–1354. (doi:10.1126/science.2999980) Crossref, PubMed, ISI, Google Scholar - 144
Mullis KB, Faloona FA . 1987 Specific synthesis of DNA in vitro via a polymerase-catalyzed chain reaction. Methods Enzymol. 155, 335–350. (doi:10.1016/0076-6879(87)55023-6) Crossref, PubMed, ISI, Google Scholar - 145
Simon C . 1991 Molecular systematics at the species boundary: exploiting conserved and variable regions of the mitochondrial genome of animals via direct sequencing from enzymatically amplified DNA. In Molecular techniques in taxonomy (edsHewitt GM, Johnston AWB, Young JPW ), pp. 33–71. New York, NY: Springer. Crossref, Google Scholar - 146
Swofford DL, Olsen GJ . 1990 Phylogeny reconstruction. In Molecular systematics (edsHillis DM, Moritz C ), pp. 411–501. Sunderland, MA: Sinauer Associates. Google Scholar - 147
Swofford DL, Olsen GJ, Waddell PJ, Hillis DM . 1996 Phylogenetic inference. In Molecular systematics (edsHillis DM, Moritz C, Mable BK ), pp. 407–514, 2nd edn. Sunderland, MA: Sinauer Assoc. Google Scholar - 148
Swofford DL . PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. Sunderland, MA: Sinauer Associates. Google Scholar - 149
Yang Z . 1996 Phylogenetic analysis by maximum likelihood (PAML). Berkeley, CA: Department of Integrative Biology, University of California at Berkeley. Google Scholar - 150
Kinzelbach R . 1990 The systematic position of Strepsiptera (Insecta). Am. Entomol. 36, 292–303. (doi:10.1093/ae/36.4.292) Crossref, Google Scholar - 151
Kathirithamby J . 1989 Review of the Order Strepsiptera. Syst. Entomol. 14, 41–92. (doi:10.1111/j.1365-3113.1989.tb00265.x) Crossref, ISI, Google Scholar - 152
Beutel RG, Gorb SN . 2006 A revised interpretation of the evolution of attachment structures in Hexapoda with special emphasis on Mantophasmatodea. Arthropod Syst. Phylogeny 64, 3–25. Crossref, Google Scholar - 153
Crowson RA . 1955 The natural classification of the families of Coleoptera. London, UK: Nathaniel Lloyd & Co., Ltd. Google Scholar - 154
McKenna DD, Farrell BD . 2010 9-genes reinforce the phylogeny of Holometabola and yield alternate views on the phylogenetic placement of Strepsiptera. PLoS ONE 5, e11887. (doi:10.1371/journal.pone.0011887) Crossref, PubMed, ISI, Google Scholar - 155
Whiting MF, Wheeler WC . 1994 Insect homeotic transformation. Nature 368, 696. (doi:10.1038/368696a0) Crossref, ISI, Google Scholar - 156
Carmean D, Crespi BJ . 1995 Do long branches attract flies? Nature 373, 666–670. (doi:10.1038/373666b0) Crossref, PubMed, ISI, Google Scholar - 157
Huelsenbeck JP . 1997 Is the Felsenstein zone a fly trap? Syst. Biol. 46, 69–74. (doi:10.1093/sysbio/46.1.69) Crossref, PubMed, ISI, Google Scholar - 158
Chalwatzis N, Bauer A, Stetzer R, Kinzelbach R, Zimmermann RK . 1995 Strongly expanded 18S ribosomal-RNA genes correlated with a peculiar morphology in the insect order of Strepsiptera. Zool. Anal. Complex Syst. 98, 115–126. ISI, Google Scholar - 159
Chalwatzis N, Hauf J, van de Peer Y, Kinzelbach R, Zimmermann RK . 1996 18S ribosomal RNA genes in insects: primary structure of the genes and molecular phylogeny of the Holometabola. Ann. Entomol. Soc. Am. 89, 788–803. (doi:10.1093/aesa/89.6.788) Crossref, ISI, Google Scholar - 160
Lockhart PJ, Steel MA, Hendy MD, Penny D . 1994 Recovering evolutionary trees under a more realistic model of sequence evolution. Mol. Biol. Evol. 11, 605–612. PubMed, ISI, Google Scholar - 161
van de Peer Y, Neefs J-M, De Rijk P, de Wachter R . 1993 Reconstructing evolution from eukaryotic small-ribosomal-subunit RNA sequences: calibration of the molecular clock. J. Mol. Evol. 37, 221–232. (doi:10.1007/BF02407359) Crossref, PubMed, ISI, Google Scholar - 162
Van de Peer Y, Nicolaï S, De Rijk P, De Wachter R . 1996 Database on the structure of small ribosomal subunit RNA. Nucleic Acids Res. 24, 86–91. (doi:10.1093/nar/24.1.86) Crossref, PubMed, ISI, Google Scholar - 163
Wheeler WC, Gladstein DL . A heuristic multiple sequence alignment procedure based on parsimony. New York, NY: American Museum of Natural History. Google Scholar - 164
Whiting MF . 1998 Long-branch distraction and the Strepsiptera. Syst. Biol. 47, 134–138. (doi:10.1080/106351598261076) Crossref, PubMed, ISI, Google Scholar - 165
Siddall ME, Whiting MF . 1999 Long-branch abstractions. Cladistics 15, 9–24. (doi:10.1111/j.1096-0031.1999.tb00391.x) Crossref, ISI, Google Scholar - 166
Whiting MF . 2002 Phylogeny of the holometabolous insect orders: molecular evidence. Zool. Scr. 31, 3–15. (doi:10.1046/j.0300-3256.2001.00093.x) Crossref, ISI, Google Scholar - 167
Whiting M . 2005 Phylogenetic position of Diptera: review of the evidence. In The evolutionary biology of flies (edsYeates DK, Wiegmann BM ), pp. 3–13. New York, NY: Columbia University Press. Google Scholar - 168
Friedrich F, Beutel RG . 2010 Goodbye Halteria? The thoracic morphology of Endopterygota (Insecta) and its phylogenetic implications. Cladistics 26, 579–612. (doi:10.1111/j.1096-0031.2010.00305.x) Crossref, ISI, Google Scholar - 169
Longhorn SJ, Pohl HW, Vogler AP . 2010 Ribosomal protein genes of holometabolan insects reject the Halteria, instead revealing a close affinity of Strepsiptera with Coleoptera. Mol. Phylogenet. Evol. 55, 846–859. (doi:10.1016/j.ympev.2010.03.024) Crossref, PubMed, ISI, Google Scholar - 170
Ishiwata K, Sasaki G, Ogawa J, Miyata T, Su Z-H . 2011 Phylogenetic relationships among insect orders based on three nuclear protein-coding gene sequences. Mol. Phylogenet. Evol. 58, 169–180. (doi:10.1016/j.ympev.2010.11.001) Crossref, PubMed, ISI, Google Scholar - 171
Huelsenbeck JP . 1998 Systematic bias in phylogenetic analysis: is the Strepsiptera problem solved?. Syst. Biol. 47, 519–537. PubMed, ISI, Google Scholar - 172
Steel M, Huson D, Lockhart PJ . 2000 Invariable sites models and their use in phylogeny reconstruction. Syst. Biol. 49, 225–232. (doi:10.1093/sysbio/49.2.225) Crossref, PubMed, ISI, Google Scholar - 173
Friedrich M, Tautz D . 1997 An episodic change of rDNA nucleotide substitution rate has occurred during the emergence of the insect order Diptera. Mol. Biol. Evol. 14, 644–653. (doi:10.1093/oxfordjournals.molbev.a025804) Crossref, PubMed, ISI, Google Scholar - 174
McMahon DP, Hayward A, Kathirithamby J . 2011 The first molecular phylogeny of Strepsiptera (Insecta) reveals an early burst of molecular evolution correlated with the transition to endoparasitism. PLoS ONE 6, e21206. (doi:10.1371/journal.pone.0021206) Crossref, PubMed, ISI, Google Scholar - 175
Niehuis O 2012 Genomic and morphological evidence converge to resolve the enigma of Strepsiptera. Curr. Biol. 22, 1309–1313. (doi:10.1016/j.cub.2012.05.018) Crossref, PubMed, ISI, Google Scholar - 176
McKenna DD . 2014 Molecular phylogenetics and evolution of Coleoptera. In Handbook of zoology. Volume IV, Arthropoda: insecta. Part 38, Coleoptera, beetles. Volume 3, morphology and systematics (Phytophaga) (edsBeutel RG, Leschen RAB ), pp. 1–10. Berlin, Germany: Walter de Gruyter. Google Scholar - 177
Boussau B 2014 Strepsiptera, phylogenomics and the long branch attraction problem. PLoS ONE 9, e107709. (doi:10.1371/journal.pone.0107709) Crossref, PubMed, ISI, Google Scholar - 178
Wheeler WC . 1995 Sequence alignment, parameter sensitivity, and the phylogenetic analysis of molecular data. Syst. Biol. 44, 321–333. (doi:10.1093/sysbio/44.3.321) Crossref, ISI, Google Scholar - 179
Wheeler WC . 1996 Optimization alignment: the end of multiple sequence alignment in phylogenetics? Cladistics 12, 1–9. (doi:10.1111/j.1096-0031.1996.tb00189.x) Crossref, ISI, Google Scholar - 180
Hickson RE, Simon C, Perrey SW . 2000 The performance of several multiple-sequence alignment programs in relation to secondary-structure features for an rRNA sequence. Mol. Biol. Evol. 17, 530–539. (doi:10.1093/oxfordjournals.molbev.a026333) Crossref, PubMed, ISI, Google Scholar - 181
Kjer KM, Gillespie JJ, Ober KA . 2006 Structural homology in ribosomal RNA, and a deliberation on POY. Arthropod Syst. Phylogeny 64, 71–76. Google Scholar - 182
Kjer KM, Gillespie JJ, Ober KA . 2007 Opinions on multiple sequence alignment, and an empirical comparison of repeatability and accuracy between POY and structural alignment. Syst. Biol. 56, 133–146. (doi:10.1080/10635150601156305) Crossref, PubMed, ISI, Google Scholar - 183
Ogden TH, Rosenberg MS . 2006 Multiple sequence alignment accuracy and phylogenetic inference. Syst. Biol. 55, 314–328. (doi:10.1080/10635150500541730) Crossref, PubMed, ISI, Google Scholar - 184
Ogden TH, Rosenberg MS . 2007 Alignment and topological accuracy of the direct optimization approach via POY and traditional phylogenetics via ClustalW+PAUP*. Syst. Biol. 56, 182–193. (doi:10.1080/10635150701281102) Crossref, PubMed, ISI, Google Scholar - 185
Kjer KM, Roshan U, Gillespie JJ . 2009 Structural and evolutionary considerations for multiple sequence alignment of RNA, and the challenges for algorithms that ignore them. See http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.122.8235 (accessed 26 Dec 2015). Google Scholar - 186
Stocsits RR, Letsch H, Hertel J, Misof B, Stadler PF . 2009 Accurate and efficient reconstruction of deep phylogenies from structured RNAs. Nucleic Acids Res. 37, 6184–6193. (doi:10.1093/nar/gkp600) Crossref, PubMed, ISI, Google Scholar - 187
Sankoff D, Cedergren RJ . 1983 Simultaneous comparison of three or more sequences related by a tree. In Time warps, string edits, and macromolecules: the theory and practice of sequence comparison (edsSankotr D, Kruskal JB ), pp. 253–263. Reading, PA: Addison-Wesley. Google Scholar - 188
Felsenstein J . 1978 Cases in which parsimony or compatibility methods will be positively misleading. Syst Zool. 27, 401–410. (doi:10.2307/2412923) Crossref, Google Scholar - 189
Hwang U-W, Kim W, Tautz D, Friedrich M . 1989 Molecular phylogenetics at the Felsenstein zone: approaching the Strepsiptera problem using 5.85 and 28S rDNA sequences. Mol. Phylogenet. Evol. 9, 470–480. (doi:10.1006/mpev.1998.0518) Crossref, ISI, Google Scholar - 190
Cameron SL . 2014 Insect mitochondrial genomics: implications for evolution and phylogeny. Annu. Rev. Entomol. 59, 95–117. (doi:10.1146/annurev-ento-011613-162007) Crossref, PubMed, ISI, Google Scholar - 191
Letsch HO, Kuck P, Stocsits RR, Misof B . 2010 The impact of rRNA secondary structure consideration in alignment and tree reconstruction: simulated data and a case study on the phylogeny of Hexapods. Mol. Biol. Evol. 27, 2507–2521. (doi:10.1093/molbev/msq140) Crossref, PubMed, ISI, Google Scholar - 192
Redelings B, Suchard M . 2005 Joint Bayesian estimation of alignment and phylogeny. Syst. Biol. 54, 401–418. (doi:10.1080/10635150590947041) Crossref, PubMed, ISI, Google Scholar - 193
Simon C, Buckley TR, Frati F, Stewart JB, Beckenbach AT . 2006 Incorporating molecular evolution into phylogenetic analysis, and a new compilation of conserved polymerase chain reaction primers for animal mitochondrial DNA. Annu. Rev. Ecol. Evol. Syst. 37, 547–579. (doi:10.1146/annurev.ecolsys.37.091305.110018) Crossref, ISI, Google Scholar - 194
Kjer KM, Swigonova Z, LaPolla JS, Broughton RE . 2007 Why weight? Mol. Phylogenet. Evol. 43, 999–1004. (doi:10.1016/j.ympev.2007.02.028) Crossref, PubMed, ISI, Google Scholar - 195
Phillips A, Janies D, Wheeler W . 2000 Multiple sequence alignment in phylogenetic analysis. Mol. Phylogenet. Evol. 16, 317–330. (doi:10.1006/mpev.2000.0785) Crossref, PubMed, ISI, Google Scholar - 196
Mickevich MF, Farris JS . 1981 The implications of congruence in Menidia. Syst. Zool. 30, 351–370. (doi:10.2307/2413255) Crossref, Google Scholar - 197
Graham SW, Kohn JR, Morton BR, Eckenwalder E, Barrett CH . 1998 Phylogenetic congruence and discordance among one morphological and three molecular data sets from Pontederiaceae. Syst. Biol. 47, 545–567. (doi:10.1080/106351598260572) Crossref, PubMed, ISI, Google Scholar - 198
Dolphin K, Belshaw R, Orme CDL, Quicke DLJ . 2000 Noise and incongruence: interpreting results of the incongruence length difference test. Mol. Phylogenet. Evol. 17, 401–406. (doi:10.1006/mpev.2000.0845) Crossref, PubMed, ISI, Google Scholar - 199
Dowton M, Austin AD . 2002 Increased incongruence does not necessarily indicate increased phylogenetic accuracy—the behavior of the ILD test in mixed-model analyses. Syst. Biol. 51, 19–31. (doi:10.1080/106351502753475853) Crossref, PubMed, ISI, Google Scholar - 200
Darlu P, Lecointre G . 2002 When does the incongruence length difference test fail? Mol. Biol. Evol. 19, 432–437. (doi:10.1093/oxfordjournals.molbev.a004098) Crossref, PubMed, ISI, Google Scholar - 201
Barker FK, Lutzoni F . 2002 The utility of the incongruence length difference test. Syst. Biol. 51, 625–637. (doi:10.1080/10635150290102302) Crossref, PubMed, ISI, Google Scholar - 202
Grant T, Kluge AG . 2003 Data exploration in phylogenetic inference: scientific, heuristic, or neither. Cladistics 19, 379–418. (doi:10.1111/j.1096-0031.2003.tb00311.x) Crossref, ISI, Google Scholar - 203
Ogden TH, Whiting MF . 2003 The problem with ‘the Palaeoptera problem:’ sense and sensitivity. Cladistics 19, 432–442. (doi:10.1111/j.1096-0031.2003.tb00313.x) Crossref, ISI, Google Scholar - 204
Terry MD, Whiting MF . 2005 Mantophasmatodea and phylogeny of the lower neopterous insects. Cladistics 21, 240–257. (doi:10.1111/j.1096-0031.2005.00062.x) Crossref, ISI, Google Scholar - 205
McKittrick FA . 1965 A contribution to the understanding of cockroach-termite affinities. Ann. Entomol. Soc. Am. 58, 18–22. (doi:10.1093/aesa/58.1.18) Crossref, PubMed, Google Scholar - 206
Dell'Ampio E 2014 Decisive data sets in phylogenomics: lessons from studies on the phylogenetic relationships of primarily wingless insects. Mol. Biol. Evol. 31, 239–249. (doi:10.1093/molbev/mst196) Crossref, PubMed, ISI, Google Scholar - 207
The Editors. 2016 Editorial. Cladistics 32, 1. Google Scholar
- 208
Yoshizawa K . 2010 Direct optimization overly optimizes data. Syst. Entomol. 35, 199–206. (doi:10.1111/j.1365-3113.2010.00526.x) Crossref, ISI, Google Scholar - 209
Felsenstein J . 1985 Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39, 783–791. (doi:10.2307/2408678) Crossref, PubMed, ISI, Google Scholar - 210
Felsenstein J . 1988 Phylogenies and quantitative characters. Annu. Rev. Ecol. Syst. 19, 445–471. (doi:10.1146/annurev.es.19.110188.002305) Crossref, Google Scholar - 211
Felsenstein J . 1988 Phylogenies from molecular sequences: inference and reliability. Annu. Rev. Genet. 22, 521–565. (doi:10.1146/annurev.ge.22.120188.002513) Crossref, PubMed, ISI, Google Scholar - 212
Felsenstein J . 1982 Numerical methods for inferring evolutionary trees. Q. Rev. Biol. 57, 379–404. (doi:10.1086/412935) Crossref, ISI, Google Scholar - 213
Felsenstein J . 1981 A likelihood approach to character weighting and what it tells us about parsimony and compatibility. Biol. J. Linn. Soc. 16, 183–196. (doi:10.1111/j.1095-8312.1981.tb01847.x) Crossref, ISI, Google Scholar - 214
Felsenstein J . 1992 Phylogenies from restriction sites: a maximum-likelihood approach. Evolution 46, 159–173. (doi:10.2307/2409811) Crossref, PubMed, ISI, Google Scholar - 215
Felsenstein J, Sober E . 1986 Parsimony and likelihood: an exchange. Syst. Biol. 35, 617–626. (doi:10.2307/2413121) Google Scholar - 216
Felsenstein J . 1989 PHYLIP—phylogeny interference package (version 3.2). Cladistics 5, 164–166. Google Scholar - 217
Guindon S, Gascuel O . 2003 A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 52, 696–704. (doi:10.1080/10635150390235520) Crossref, PubMed, ISI, Google Scholar - 218
Zwickl DJ . 2006 Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. PhD dissertation, The University of Texas at Austin, Austin, TX. Google Scholar - 219
Stamatakis A . 2006 RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics Oxf. Engl. 22, 2688–2690. (doi:10.1093/bioinformatics/btl446) Crossref, PubMed, ISI, Google Scholar - 220
Hendy MD, Penny D . 1989 A framework for the quantitative study of evolutionary trees. Syst. Zool. 38, 297–309. (doi:10.2307/2992396) Crossref, Google Scholar - 221
Tuffley C, Steel M . 1997 Links between maximum likelihood and maximum parsimony under a simple model of site substitution. Bull. Math. Bio. 59, 581–607. (doi:10.1007/BF02459467) Crossref, Google Scholar - 222
Larget B, Simon DL . 1999 Markov Chain Monte Carlo algorithms for the Bayesian analysis of phylogenetic trees. Mol. Biol. Evol. 16, 750. (doi:10.1093/oxfordjournals.molbev.a026160) Crossref, ISI, Google Scholar - 223
Huelsenbeck JP, Ronquist F . 2001 MrBayes: Bayesian inference of phylogeny. Bioinformatics 17, 754–755. Crossref, PubMed, ISI, Google Scholar - 224
Swofford DL, Waddell PJ, Huelsenbeck JP, Foster PG, Lewis PO, Rogers JS . 2001 Bias in phylogenetic estimation and its relevance to the choice between parsimony and likelihood methods. Syst. Biol. 50, 525–539. (doi:10.1080/106351501750435086) Crossref, PubMed, ISI, Google Scholar - 225
Yang Z . 1994 Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods. J. Mol. Evol. 39, 306–314. (doi:10.1007/BF00160154) Crossref, PubMed, ISI, Google Scholar - 226
Sullivan J, Holsinger KE, Simon C . 1996 The effect of topology on estimates of among-site rate variation. J. Mol. Evol. 42, 308–312. (doi:10.1007/BF02198857) Crossref, PubMed, ISI, Google Scholar - 227
Yang Z . 1996 Among-site rate variation and its impact on phylogenetic analyses. Trends Ecol. Evol. 11, 367–372. (doi:10.1016/0169-5347(96)10041-0) Crossref, PubMed, ISI, Google Scholar - 228
Frati F, Simon C, Sullivan J, Swofford DL . 1997 Evolution of the mitochondrial cytochrome oxidase II gene in Collembola. J. Mol. Evol. 44, 145–158. (doi:10.1007/PL00006131) Crossref, PubMed, ISI, Google Scholar - 229
Sullivan J, Holsinger KE, Simon C . 1995 Among-site rate variation and phylogenetic analysis of 12S rRNA in Sigmodontine rodents. Mol. Biol. Evol. 12, 988–1001. PubMed, ISI, Google Scholar - 230
Posada D, Crandall KA . 1998 MODELTEST: testing the model of DNA substitution. Bioinformatics 14, 817–818. (doi:10.1093/bioinformatics/14.9.817) Crossref, PubMed, ISI, Google Scholar - 231
Sullivan J, Joyce P . 2005 Model selection in phylogenetics. Annu. Rev. Ecol. Evol. Syst. 36, 445–466. (doi:10.1146/annurev.ecolsys.36.102003.152633) Crossref, ISI, Google Scholar - 232
Friedrich M, Tautz D . 1995 Ribosomal DNA phylogeny of the major extant arthropod classes and the evolution of myriapods. Nature 376, 165–167. (doi:10.1038/376165a0) Crossref, PubMed, ISI, Google Scholar - 233
von Dohlen CD, Moran NA . 1995 Molecular phylogeny of the Homoptera: a paraphyletic taxon. J. Mol. Evol. 41, 211–223. (doi:10.1007/BF00170675) Crossref, PubMed, ISI, Google Scholar - 234
Flook PK, Rowell CHF . 1997 The effectiveness of mitochondrial rRNA sequences for the reconstruction of the phylogeny of an insect order (Orthoptera). Mol. Phylogenet. Evol. 8, 177–192. (doi:10.1006/mpev.1997.0425) Crossref, PubMed, ISI, Google Scholar - 235
Whitfield JB, Cameron SA . 1998 Hierarchical analysis of variation in the mitochondrial 16S rRNA gene among Hymenoptera. Mol. Biol. Evol. 15, 1728–1743. (doi:10.1093/oxfordjournals.molbev.a025899) Crossref, PubMed, ISI, Google Scholar - 236
Lo N, Tokuda G, Watanabe H, Rose H, Slaytor M, Maekawa K, Bandi C, Noda H . 2000 Evidence from multiple gene sequences indicates that termites evolved from wood-feeding cockroaches. Curr. Biol. 10, 801–804. (doi:10.1016/S0960-9822(00)00561-3) Crossref, PubMed, ISI, Google Scholar - 237
Hull D . 1988 Science as a process: an evolutionary account of the social and conceptual development of science. Chicago: University of Chicago Press. Google Scholar - 238
Schöniger M, von Haeseler A . 1994 A stochastic model for the evolution of autocorrelated DNA sequences. Mol. Phylogenet. Evol. 3, 240–247. (doi:10.1006/mpev.1994.1026) Crossref, PubMed, ISI, Google Scholar - 239
Misof B, Niehuis O, Bischoff I, Rickert A, Erpenbeck D, Staniczek A . 2007 Towards an 18S phylogeny of hexapods: accounting for group-specific character covariance in optimized mixed nucleotide/doublet models. Zool. Jena Ger. 110, 409–429. (doi:10.1016/j.zool.2007.08.003) PubMed, ISI, Google Scholar - 240
Misof B, Misof K . 2009 A Monte Carlo approach successfully identifies randomness in multiple sequence alignments : a more objective means of data exclusion. Syst. Biol. 58, 21–34. (doi:10.1093/sysbio/syp006) Crossref, PubMed, ISI, Google Scholar - 241
Xie Q, Tian X, Qin Y, Bu W . 2009 Phylogenetic comparison of local length plasticity of the small subunit of nuclear rDNAs among all Hexapoda orders and the impact of hyper-length-variation on alignment. Mol. Phylogenet. Evol. 50, 310–316. (doi:10.1016/j.ympev.2008.10.025) Crossref, PubMed, ISI, Google Scholar - 242
Luan Y, Mallatt JM, Xie R, Yang Y, Yin W . 2005 The phylogenetic positions of three basal-hexapod groups (Protura, Diplura, and Collembola) based on ribosomal RNA gene sequences. Mol. Biol. Evol. 22, 1579–1592. (doi:10.1093/molbev/msi148) Crossref, PubMed, ISI, Google Scholar - 243
Giribet G, Edgecombe GD, Carpenter JM, D'Haese CA, Wheeler WC . 2004 Is Ellipura monophyletic? A combined analysis of basal hexapod relationships with emphasis on the origin of insects. Org. Divers. Evol. 4, 319–340. (doi:10.1016/j.ode.2004.05.001) Crossref, ISI, Google Scholar - 244
Gao Y, Bu Y, Luan Y . 2008 Phylogenetic relationships of basal hexapods reconstructed from nearly complete 18S and 28S rRNA gene sequences. Zool. Sci. 25, 1139–1145. (doi:10.2108/zsj.25.1139) Crossref, PubMed, ISI, Google Scholar - 245
Mallatt J, Craig CW, Yoder MJ . 2010 Nearly complete rRNA genes assembled from across the metazoan animals: effects of more taxa, a structure-based alignment, and paired-sites evolutionary models on phylogeny reconstruction. Mol. Phylogenet. Evol. 55, 1–17. (doi:10.1016/j.ympev.2009.09.028) Crossref, PubMed, ISI, Google Scholar - 246
Whiting MF . 2002 Mecoptera is paraphyletic: multiple genes and phylogeny of Mecoptera and Siphonaptera. Zool. Scr. 31, 93–104. (doi:10.1046/j.0300-3256.2001.00095.x) Crossref, ISI, Google Scholar - 247
Rokas A, Kathirithamby J, Holland PWH . 1999 Intron insertion as a phylogenetic character: the engrailed homeobox of Strepsiptera does not indicate affinity with Diptera. Insect Mol. Biol. 8, 527–530. (doi:10.1046/j.1365-2583.1999.00149.x) Crossref, PubMed, ISI, Google Scholar - 248
Carapelli A, Frati F, Nardi F, Dallai R, Simon C . 2000 Molecular phylogeny of the apterygotan insects based on nuclear and mitochondrial genes. Pedobiologia 44, 361–373. (doi:10.1078/S0031-4056(04)70055-4) Crossref, ISI, Google Scholar - 249
Krauss V . 2004 Phylogenetic mapping of intron positions: a case study of translation initiation factor eIF2. Mol. Biol. Evol. 22, 74–84. (doi:10.1093/molbev/msh255) Crossref, PubMed, ISI, Google Scholar - 250
Djernæs M, Damgaard J . 2006 Exon-intron structure, paralogy and sequenced regions of elongation factor-1 alpha in Hexapoda. Arthropod Syst. Phylogeny 64, 45–52. Google Scholar - 251
Bonneton F, Brunet FG, Kathirithamby J, Laudet V . 2006 The rapid divergence of the ecdysone receptor is a synapomorphy for Mecopterida that clarifies the Strepsiptera problem. Insect Mol. Biol. 15, 351–362. (doi:10.1111/j.1365-2583.2006.00654.x) Crossref, PubMed, ISI, Google Scholar - 252
Predel R, Roth S . 2007 Neuropeptide evolution and the analysis of phylogenetic relationships in Blattaria (Hexapoda). Arthropod Syst. Phylogeny 65, 3–6. Google Scholar - 253
Roth S, Köhler G, Reinhardt K, Predel R . 2007 A discrete neuropeptide difference between two hybridizing grasshopper subspecies. Biol. J. Linn. Soc. 91, 541–548. (doi:10.1111/j.1095-8312.2007.00865.x) Crossref, ISI, Google Scholar - 254
Roth S, Fromm B, Gäde G, Predel R . 2009 A proteomic approach for studying insect phylogeny: CAPA peptides of ancient insect taxa (Dictyoptera, Blattoptera) as a test case. BMC Evol. Biol. 9, 50. (doi:10.1186/1471-2148-9-50) Crossref, PubMed, ISI, Google Scholar - 255
Predel R, Neupert S, Huetteroth W, Kahnt J, Waidelich D, Roth S . 2012 Peptidomics-based phylogeny and biogeography of Mantophasmatodea (Hexapoda). Syst. Biol. 61, 609–629. (doi:10.1093/sysbio/sys003) Crossref, PubMed, ISI, Google Scholar - 256
Boore JL, Collins TM, Stanton D, Daehler LL, Brown WM . 1995 Deducing the pattern of arthropod phytogeny from mitochondrial DNA rearrangements. Nature 376, 163–165. (doi:10.1038/376163a0) Crossref, PubMed, ISI, Google Scholar - 257
Boore JL, Lavrov DV, Brown WM . 1998 Gene translocation links insects and crustaceans. Nature 392, 667–668. (doi:10.1038/33577) Crossref, PubMed, ISI, Google Scholar - 258
Shao R, Campbell NJH, Barker SC . 2001 Numerous gene rearrangements in the mitochondrial genome of the Wallaby louse, Heterodoxus macropus (Phthiraptera). Mol. Biol. Evol. 18, 858–865. (doi:10.1093/oxfordjournals.molbev.a003867) Crossref, PubMed, ISI, Google Scholar - 259
Shao R, Barker SC . 2003 The highly rearranged mitochondrial genome of the plague thrips, Thrips imaginis (Insecta: Thysanoptera): convergence of two novel gene boundaries and an extraordinary arrangement of rRNA genes. Mol. Biol. Evol. 20, 362–370. (doi:10.1093/molbev/msg045) Crossref, PubMed, ISI, Google Scholar - 260
Thao ML, Baumann L, Baumann P . 2004 Organization of the mitochondrial genomes of whiteflies, aphids, and psyllids (Hemiptera, Sternorrhyncha). BMC Evol. Biol. 4, 25. (doi:10.1186/1471-2148-4-25) Crossref, PubMed, ISI, Google Scholar - 261
Wei D-D, Shao R, Yuan M-L, Dou W, Barker SC, Wang J-J . 2012 The multipartite mitochondrial genome of Liposcelis bostrychophila: insights into the evolution of mitochondrial genomes in bilateral animals. PLoS ONE 7, e33973. (doi:10.1371/journal.pone.0033973) Crossref, PubMed, ISI, Google Scholar - 262
Dowton M, Cameron SL, Austin AD, Whiting MF . 2009 Phylogenetic approaches for the analysis of mitochondrial genome sequence data in the Hymenoptera–a lineage with both rapidly and slowly evolving mitochondrial genomes. Mol. Phylogenet. Evol. 52, 512–519. (doi:10.1016/j.ympev.2009.04.001) Crossref, PubMed, ISI, Google Scholar - 263
Dowton M, Cameron SL, Dowavic JI, Austin AD, Whiting MF . 2009 Characterization of 67 mitochondrial tRNA gene rearrangements in the hymenoptera suggests that mitochondrial tRNA gene position is selectively neutral. Mol. Biol. Evol. 26, 1607–1617. (doi:10.1093/molbev/msp072) Crossref, PubMed, ISI, Google Scholar - 264
Mandal SD, Chhakchhuak L, Gurusubramanian G, Kumar NS . 2014 Mitochondrial markers for identification and phylogenetic studies in insects—a review. DNA Barcodes 2, 1. (doi:10.2478/dna-2014-0001) Crossref, Google Scholar - 265
Nardi F, Carapelli A, Fanciulli PP, Dallai R, Frati F . 2001 The complete mitochondrial DNA sequence of the basal hexapod Tetrodontophora bielanensis: evidence for heteroplasmy and tRNA translocations. Mol. Biol. Evol. 18, 1293–1304. (doi:10.1093/oxfordjournals.molbev.a003914) Crossref, PubMed, ISI, Google Scholar - 266
Friedrich M, Muqim N . 2003 Sequence and phylogenetic analysis of the complete mitochondrial genome of the flour beetle Tribolium castanaeum. Mol. Phylogenet. Evol. 26, 502–512. (doi:10.1016/S1055-7903(02)00335-4) Crossref, PubMed, ISI, Google Scholar - 267
Bae JS, Kim I, Sohn HD, Jin BR . 2004 The mitochondrial genome of the firefly, Pyrocoelia rufa: complete DNA sequence, genome organization, and phylogenetic analysis with other insects. Mol. Phylogenet. Evol. 32, 978–985. (doi:10.1016/j.ympev.2004.03.009) Crossref, PubMed, ISI, Google Scholar - 268
Kim I, Cha SY, Yoon MH, Hwang JS, Lee SM, Sohn HD, Jin BR . 2005 The complete nucleotide sequence and gene organization of the mitochondrial genome of the oriental mole cricket, Gryllotalpa orientalis (Orthoptera: Gryllotalpidae). Gene 353, 155–168. (doi:10.1016/j.gene.2005.04.019) Crossref, PubMed, ISI, Google Scholar - 269
Stewart JB, Beckenbach AT . 2005 Insect mitochondrial genomics: the complete mitochondrial genome sequence of the meadow spittlebug Philaenus spumarius (Hemiptera: Auchenorrhyncha: Cercopoidae). Genome 48, 46–54. (doi:10.1139/g04-090) Crossref, PubMed, ISI, Google Scholar - 270
Liu S 2016 Mitochondrial capture enriches mito-DNA 100 fold, enabling PCR-free mitogenomics biodiversity analysis. Mol. Ecol. Resour. 16, 470–479. (doi:10.1111/1755-0998.12472) Crossref, PubMed, ISI, Google Scholar - 271
Cameron SL, Miller KB, D'Haese CA, Whiting MF, Barker SC . 2004 Mitochondrial genome data alone are not enough to unambiguously resolve the relationships of Entognatha, Insecta and Crustacea sensu lato (Arthropoda). Cladistics 20, 534–557. (doi:10.1111/j.1096-0031.2004.00040.x) Crossref, ISI, Google Scholar - 272
Cameron SL, Beckenbach AT, Downton MA, Whiting MF . 2006 Evidence from mitochondrial genomics on interordinal relationships in insects. Arthropod Syst. Phylogeny 64, 27–34. Google Scholar - 273
Cameron SL, Barker SC, Whiting MF . 2006 Mitochondrial genomics and the new insect order Mantophasmatodea. Mol. Phylogenet. Evol. 38, 274–279. (doi:10.1016/j.ympev.2005.09.020) Crossref, PubMed, ISI, Google Scholar - 274
Simon S, Hadrys H . 2013 A comparative analysis of complete mitochondrial genomes among Hexapoda. Mol. Phylogenet. Evol. 69, 393–403. (doi:10.1016/j.ympev.2013.03.033) Crossref, PubMed, ISI, Google Scholar - 275
Talavera G, Vila R . 2011 What is the phylogenetic signal limit from mitogenomes? The reconciliation between mitochondrial and nuclear data in the Insecta class phylogeny. BMC Evol. Biol. 11, 315. (doi:10.1186/1471-2148-11-315) Crossref, PubMed, ISI, Google Scholar - 276
Chen W-J, Koch M, Mallatt JM, Luan Y-X . 2014 comparative analysis of mitochondrial genomes in Diplura (Hexapoda, Arthropoda): taxon sampling is crucial for phylogenetic inferences. Genome Biol. Evol. 6, 105–120. (doi:10.1093/gbe/evt207) Crossref, PubMed, ISI, Google Scholar - 277
Carapelli A, Vannini L, Nardi F, Boore JL, Beani L, Dallai R, Frati F . 2006 The mitochondrial genome of the entomophagous endoparasite Xenos vesparum (Insecta: Strepsiptera). Gene 376, 248–259. (doi:10.1016/j.gene.2006.04.005) Crossref, PubMed, ISI, Google Scholar - 278
Zhang J, Zhou C, Gai Y, Song D, Zhou K . 2008 The complete mitochondrial genome of Parafronurus youi (Insecta: Ephemeroptera) and phylogenetic position of the Ephemeroptera. Gene 424, 18–24. (doi:10.1016/j.gene.2008.07.037) Crossref, PubMed, ISI, Google Scholar - 279
Zhang Y, Xuan W, Zhao J, Zhu C, Jiang G . 2010 The complete mitochondrial genome of the cockroach Eupolyphaga sinensis (Blattaria: Polyphagidae) and the phylogenetic relationships within the Dictyoptera. Mol. Biol. Rep. 37, 3509–3516. (doi:10.1007/s11033-009-9944-1) Crossref, PubMed, ISI, Google Scholar - 280
Plazzi F, Ricci A, Passamonti M . 2011 The mitochondrial genome of Bacillus stick insects (Phasmatodea) and the phylogeny of orthopteroid insects. Mol. Phylogenet. Evol. 58, 304–316. (doi:10.1016/j.ympev.2010.12.005) Crossref, PubMed, ISI, Google Scholar - 281
Wan X, Kim MI, Kim MJ, Kim I . 2012 Complete mitochondrial genome of the free-living earwig, Challia fletcheri (Dermaptera: Pygidicranidae) and phylogeny of Polyneoptera. PLoS ONE 7, e42056. (doi:10.1371/journal.pone.0042056) Crossref, PubMed, ISI, Google Scholar - 282
Cameron SL, Sullivan J, Song H, Miller KB, Whiting MF . 2009 A mitochondrial genome phylogeny of the Neuropterida (lace-wings, alderflies and snakeflies) and their relationship to the other holometabolous insect orders. Zool. Scr. 38, 575–590. (doi:10.1111/j.1463-6409.2009.00392.x) Crossref, ISI, Google Scholar - 283
Rehn AC . 2003 Phylogenetic analysis of higher-level relationships of Odonata. Syst. Entomol. 28, 181–239. (doi:10.1046/j.1365-3113.2003.00210.x) Crossref, ISI, Google Scholar - 284
Ware J, May M, Kjer K . 2007 Phylogeny of the higher Libelluloidea (Anisoptera: Odonata): an exploration of the most speciose superfamily of dragonflies. Mol. Phylogenet. Evol. 45, 289–310. (doi:10.1016/j.ympev.2007.05.027) Crossref, PubMed, ISI, Google Scholar - 285
Bybee SM, Ogden TH, Branham MA, Whiting MF . 2008 Molecules, morphology and fossils: a comprehensive approach to odonate phylogeny and the evolution of the odonate wing. Cladistics 24, 477–514. (doi:10.1111/j.1096-0031.2007.00191.x) Crossref, ISI, Google Scholar - 286
Carle FL, Kjer KM, May ML . 2008 Evolution of Odonata, with special reference to Coenagrionoidea (Zygoptera). Arthropod Syst. Phylogeny 66, 37–44. Google Scholar - 287
Fleck G, Brenk M, Misof B . 2008 Larval and molecular characters help to solve phylogenetic puzzles in the highly diverse dragonfly family Libellulidae (Insecta: Odonata: Anisoptera): the Tetrathemistinae are a polyphyletic group. Org. Divers. Evol. 8, 1–16. (doi:10.1016/j.ode.2006.08.003) Crossref, ISI, Google Scholar - 288
Dumont HJ, Vierstraete A, Vanfleteren JR . 2010 A molecular phylogeny of the Odonata (Insecta). Syst. Entomol. 35, 6–18. (doi:10.1111/j.1365-3113.2009.00489.x) Crossref, ISI, Google Scholar - 289
Davis RB, Nicholson DB, Saunders EL, Mayhew PJ . 2011 Fossil gaps inferred from phylogenies alter the apparent nature of diversification in dragonflies and their relatives. BMC Evol. Biol. 11, 252. (doi:10.1186/1471-2148-11-252) Crossref, PubMed, ISI, Google Scholar - 290
Blanke A, Greve C, Mokso R, Beckman F, Misof B . 2013 An updated phylogeny of Anisoptera including formal convergence analysis of morphological characters. Syst. Entomol. 38, 474–490. (doi:10.1111/syen.12012) Crossref, ISI, Google Scholar - 291
Dijkstra K-DB, Kalkman VJ, Dow RA, Stokvis FR, Van Tol J . 2014 Redefining the damselfly families: a comprehensive molecular phylogeny of Zygoptera (Odonata): molecular phylogeny of damselflies. Syst. Entomol. 39, 68–96. (doi:10.1111/syen.12035) Crossref, ISI, Google Scholar - 292
Carle FL, Kjer KM, May ML . 2015 A molecular phylogeny and classification of Anisoptera (Odonata). Arthropod Syst. Phylogeny 73, 281–301. ISI, Google Scholar - 293
Ogden TH, Gattolliat JL, Sartori M, Staniczek AH, Soldán T, Whiting MF . 2009 Towards a new paradigm in mayfly phylogeny (Ephemeroptera): combined analysis of morphological and molecular data. Syst. Entomol. 34, 616–634. (doi:10.1111/j.1365-3113.2009.00488.x) Crossref, ISI, Google Scholar - 294
Sivaramakrishnan KG, Subramanian KA, Arunachalam M, Kumar CS, Sundar S . 2011 Emerging trends in molecular systematics and molecular phylogeny of mayflies (Insecta: Ephemeroptera). J. Threat Taxa 3, 1975–1980. (doi:10.11609/JoTT.o2661.1975-80) Crossref, Google Scholar - 295
Qian YH, Wu HY, Ji XY, Yu WW, Du YZ . 2014 Mitochondrial genome of the stonefly Kamimuria wangi (Plecoptera: Perlidae) and phylogenetic position of Plecoptera based on mitogenomes. PLoS ONE 9, e86328. (doi:10.1371/journal.pone.0086328) Crossref, PubMed, ISI, Google Scholar - 296
Kocarek P, John V, Hulva P . 2013 When the body hides the ancestry: phylogeny of morphologically modified epizoic earwigs based on molecular evidence. PLoS ONE 8, e66900. (doi:10.1371/journal.pone.0066900) Crossref, PubMed, ISI, Google Scholar - 297
Miller KB, Hayashi C, Whiting MF, Svenson GJ, Edgerly JS . 2012 The phylogeny and classification of Embioptera (Insecta). Syst. Entomol. 37, 550–570. (doi:10.1111/j.1365-3113.2012.00628.x) Crossref, ISI, Google Scholar - 298
Bradler S, Robertson JA, Whiting MF . 2014 A molecular phylogeny of Phasmatodea with emphasis on Necrosciinae, the most species-rich subfamily of stick insects. Syst. Entomol. 39, 205–222. (doi:10.1111/syen.12055) Crossref, ISI, Google Scholar - 299
Ware JL, Litman J, Klass K-D, Spearman LA . 2008 Relationships among the major lineages of Dictyoptera: the effect of outgroup selection on dictyopteran tree topology. Syst. Entomol. 33, 429–450. (doi:10.1111/j.1365-3113.2008.00424.x) Crossref, ISI, Google Scholar - 300
Legendre F, Nel A, Svenson GJ, Robillard T, Pellens R, Grandcolas P . 2015 Phylogeny of Dictyoptera: dating the origin of cockroaches, praying mantises and termites with molecular data and controlled fossil evidence. PLoS ONE 10, e0130127. (doi:10.1371/journal.pone.0130127) Crossref, PubMed, ISI, Google Scholar - 301
Inward D, Beccaloni G, Eggleton P . 2007 Death of an order: a comprehensive molecular phylogenetic study confirms that termites are eusocial cockroaches. Biol. Lett. 3, 331. (doi:10.1098/rsbl.2007.0102) Link, ISI, Google Scholar - 302
Svenson GJ, Whiting MF . 2009 Reconstructing the origins of praying mantises (Dictyoptera, Mantodea): the roles of Gondwanan vicariance and morphological convergence. Cladistics 25, 468–514. (doi:10.1111/j.1096-0031.2009.00263.x) Crossref, ISI, Google Scholar - 303
Flook PK, Klee S, Rowell CHF . 1999 Combined molecular phylogenetic analysis of the Orthoptera (Arthropoda, Insecta) and implications for their higher systematics. Syst. Biol. 48, 233–253. (doi:10.1080/106351599260274) Crossref, PubMed, ISI, Google Scholar - 304
Zhang H-L, Huang Y, Lin L-L, Wang X-Y, Zheng Z-M . 2013 The phylogeny of the Orthoptera (Insecta) as deduced from mitogenomic gene sequences. Zool. Stud. 52, 1–13. (doi:10.1186/1810-522X-52-1) Crossref, ISI, Google Scholar - 305
Song H . 2010 Grasshopper systematics: past, present and future. J. Orthoptera Res. 19, 57–68. (doi:10.1665/034.019.0112) Crossref, Google Scholar - 306
Song H . 2015 300 million years of diversification: elucidating the patterns of orthopteran evolution based on comprehensive taxon and gene sampling. Cladistics 31, 621–651. (doi:10.1111/cla.12116) Crossref, PubMed, ISI, Google Scholar - 307
Song N, Liang A-P, Bu C-P . 2012 A molecular phylogeny of Hemiptera inferred from mitochondrial genome sequences. PLoS ONE 7, e48778. (doi:10.1371/journal.pone.0048778) Crossref, PubMed, ISI, Google Scholar - 308
Cryan JR, Urban JM . 2012 Higher-level phylogeny of the insect order Hemiptera: is Auchenorrhyncha really paraphyletic? Syst. Entomol. 37, 7–21. (doi:10.1111/j.1365-3113.2011.00611.x) Crossref, ISI, Google Scholar - 309
Johnson KP, Yoshizawa K, Smith VS . 2004 Multiple origins of parasitism in lice. Proc. R. Soc. Lond. B 271, 1771–1776. (doi:10.1098/rspb.2004.2798) Link, ISI, Google Scholar - 310
Cameron SA, Mardulyn P . 2001 Multiple molecular data sets suggest independent origins of highly eusocial behavior in bees (Hymenoptera:Apinae). Syst. Biol. 50, 194–214. (doi:10.1080/10635150151125851) Crossref, PubMed, ISI, Google Scholar - 311
Hines HM, Hunt JH, O'Connor TK, Gillespie JJ, Cameron SA . 2007 Multigene phylogeny reveals eusociality evolved twice in vespid wasps. Proc. Natl Acad. Sci. USA 104, 3295–3299. (doi:10.1073/pnas.0610140104) Crossref, PubMed, ISI, Google Scholar - 312
Peters RS, Meyer B, Krogmann L, Borner J, Meusemann K, Schütte K, Niehuis O, Misof B . 2011 The taming of an impossible child: a standardized all-in approach to the phylogeny of Hymenoptera using public database sequences. BMC Biol. 9, 55. (doi:10.1186/1741-7007-9-55) Crossref, PubMed, ISI, Google Scholar - 313
Sharkey MJ 2011 Phylogenetic relationships among superfamilies of Hymenoptera. Cladistics 28, 80–112. (doi:10.1111/j.1096-0031.2011.00366.x) Crossref, ISI, Google Scholar - 314
Klopfstein S, Vilhelmsen L, Heraty JM, Sharkey M, Ronquist F . 2013 The hymenopteran tree of life: evidence from protein-coding genes and objectively aligned ribosomal data. PLoS ONE 8, e69344. (doi:10.1371/journal.pone.0069344) Crossref, PubMed, ISI, Google Scholar - 315
Johnson BR, Borowiec ML, Chiu JC, Lee EK, Atallah J, Ward PS . 2013 Phylogenomics resolves evolutionary relationships among ants, bees, and wasps. Curr. Biol. 23, 2058–2062. (doi:10.1016/j.cub.2013.08.050) Crossref, PubMed, ISI, Google Scholar - 316
Mao M, Gibson T, Dowton M . 2014 Evolutionary dynamics of the mitochondrial genome in the Evaniomorpha (Hymenoptera)—a group with an intermediate rate of gene rearrangement. Genome Biol. Evol. 6, 1862–1874. (doi:10.1093/gbe/evu145) Crossref, PubMed, ISI, Google Scholar - 317
Ward PS . 2014 The phylogeny and evolution of ants. Annu. Rev. Ecol. Evol. Syst. 45, 23–43. (doi:10.1146/annurev-ecolsys-120213-091824) Crossref, ISI, Google Scholar - 318
Winterton SL, Hardy NB, Wiegmann BM . 2010 On wings of lace. Syst. Entomol. 35, 49–378. (doi:10.1111/j.1365-3113.2010.00521.x) Crossref, ISI, Google Scholar - 319
Aspöck U, Haring E, Aspöck H . 2012 The phylogeny of the Neuropterida: long lasting and current controversies and challenges (Insecta: Endopterygota). Arthropod Syst. Phylogeny 70, 119–129. ISI, Google Scholar - 320
Aspöck U . 2002 Phylogeny of the Neuropterida (Insecta: Holometabola). Zool. Scr. 31, 51–55. (doi:10.1046/j.0300-3256.2001.00087.x) Crossref, ISI, Google Scholar - 321
Mckenna DD 2015 The beetle tree of life reveals that Coleoptera survived end-Permian mass extinction to diversify during the Cretaceous terrestrial revolution. Syst. Entomol. 40, 835–880. (doi:10.1111/syen.12132) Crossref, ISI, Google Scholar - 322
Bocak L, Barton C, Crampton-Platt A, Chesters D, Ahrens D, Vogler AP . 2014 Building the Coleoptera tree-of-life for >8000 species: composition of public DNA data and fit with Linnaean classification. Syst. Entomol. 39, 97–110. (doi:10.1111/syen.12037) Crossref, ISI, Google Scholar - 323
Lawrence JF, Ślipiński A, Seago AE, Thayer MK, Newton AF, Marvaldi AE . 2011 Phylogeny of the Coleoptera based on morphological characters of adults and larvae. Ann. Zool. 61, 1–217. (doi:10.3161/000345411X576725) Crossref, ISI, Google Scholar - 324
Wiegmann BM, Yeates DK, Thorne JL, Kishino H . 2003 Time flies, a new molecular time-scale for brachyceran fly evolution without a clock. Syst. Biol. 52, 745–756. (doi:10.1093/sysbio/52.6.745) Crossref, PubMed, ISI, Google Scholar - 325
Wiegmann BM 2011 Episodic radiations in the fly tree of life. Proc. Natl Acad. Sci. USA 108, 5690–5695. (doi:10.1073/pnas.1012675108) Crossref, PubMed, ISI, Google Scholar - 326
Wiegmann BM, Regier JC, Mitter C . 2002 Combined molecular and morphological evidence on the phylogeny of the earliest lepidopteran lineages. Zool. Scr. 31, 67–81. (doi:10.1046/j.0300-3256.2001.00091.x) Crossref, ISI, Google Scholar - 327
Wahlberg N 2005 Synergistic effects of combining morphological and molecular data in resolving the phylogeny of butterflies and skippers. Proc. R. Soc. B. 272, 1577–1586. (doi:10.1098/rspb.2005.3124) Link, ISI, Google Scholar - 328
Wahlberg N, Wheat CW . 2008 Genomic outposts serve the phylogenomic pioneers: designing novel nuclear markers for genomic DNA extractions of Lepidoptera. Syst. Biol. 57, 231–242. (doi:10.1080/10635150802033006) Crossref, PubMed, ISI, Google Scholar - 329
Timmermans MJTN, Lees DC, Simonsen TJ . 2014 Towards a mitogenomic phylogeny of Lepidoptera. Mol. Phylogenet. Evol. 79, 169–178. (doi:10.1016/j.ympev.2014.05.031) Crossref, PubMed, ISI, Google Scholar - 330
Regier JC 2013 A large-scale, higher-level, molecular phylogenetic study of the insect order Lepidoptera (moths and butterflies). PLoS ONE 8, e58568. (doi:10.1371/journal.pone.0058568) Crossref, PubMed, ISI, Google Scholar - 331
Regier JC 2015 A molecular phylogeny for the oldest (nonditrysian) lineages of extant Lepidoptera, with implications for classification, comparative morphology and life-history evolution. Syst. Entomol. 40, 671–704. (doi:10.1111/syen.12129) Crossref, ISI, Google Scholar - 332
Kawahara AY, Breinholt JW . 2014 Phylogenomics provides strong evidence for relationships of butterflies and moths. Proc. R. Soc. B. 281, 20140970. (doi:10.1098/rspb.2014.0970) Link, ISI, Google Scholar - 333
Heikkilä M, Mutanen M, Wahlberg N, Sihvonen P, Kaila L . 2015 Elusive ditrysian phylogeny: an account of combining systematized morphology with molecular data (Lepidoptera). BMC Evol Biol. 15, 260. (doi:10.1186/s12862-015-0520-0) Crossref, PubMed, ISI, Google Scholar - 334
Simonsen TJ, de Jong R, Heikkilä M, Kaila L . 2012 Butterfly morphology in a molecular age—does it still matter in butterfly systematics? Arthropod Struct. Dev. 41, 307–322. (doi:10.1016/j.asd.2012.04.006) Crossref, PubMed, ISI, Google Scholar - 335
Kjer KM, Blahnik RJ, Holzenthal RW . 2001 Phylogeny of Trichoptera (caddisflies): characterization of signal and noise within multiple datasets. Syst. Biol. 50, 781–816. (doi:10.1080/106351501753462812) Crossref, PubMed, ISI, Google Scholar - 336
Holzenthal RW, Blahnik RJ, Kjer KM, Prather AL . 2007 An update on the phylogeny of caddisflies (Trichoptera). In Proc. of the XIIth International Symposium on Trichoptera (edsBueno-Soria J, Barba-Alvearz R, Armitage B ). Columbus, OH: The Caddis Press. Google Scholar - 337
Malm T, Johanson KA, Wahlberg N . 2013 The evolutionary history of Trichoptera (Insecta): a case of successful adaptation to life in freshwater. Syst. Entomol. 38, 459–473. (doi:10.1111/syen.12016) Crossref, ISI, Google Scholar - 338
Hu G-L, Yan G, Xu H, Hua B-Z . 2015 Molecular phylogeny of Panorpidae (Insecta: Mecoptera) based on mitochondrial and nuclear genes. Mol. Phylogenet. Evol. 85, 22–31. (doi:10.1016/j.ympev.2015.01.009) Crossref, PubMed, ISI, Google Scholar - 339
Whiting MF, Whiting AS, Hastriter MW . 2003 A comprehensive phylogeny of Mecoptera and Siphonaptera. Entomol. Abh. 61, 169. Google Scholar - 340
Whiting MF, Whiting AS, Hastriter MW, Dittmar K . 2008 A molecular phylogeny of fleas (Insecta: Siphonaptera): origins and host associations. Cladistics 24, 677–707. (doi:10.1111/j.1096-0031.2008.00211.x) Crossref, ISI, Google Scholar - 341
Zhu Q, Hastriter MW, Whiting MF, Dittmar K . 2015 Fleas (Siphonaptera) are Cretaceous, and evolved with Theria. Mol. Phylogenet. Evol. 90, 129–139. (doi:10.1016/j.ympev.2015.04.027) Crossref, PubMed, ISI, Google Scholar - 342
Savard J, Tautz D, Richards S, Weinstock GM, Gibbs RA, Werren JH, Tettelin H, Lercher MJ . 2006 Phylogenomic analysis reveals bees and wasps (Hymenoptera) at the base of the radiation of Holometabolous insects. Genome Res. 16, 1334–1338. (doi:10.1101/gr.5204306) Crossref, PubMed, ISI, Google Scholar - 343
Sasaki G, Ishiwata K, Machida R, Miyata T, Su Z-H . 2013 Molecular phylogenetic analyses support the monophyly of Hexapoda and suggest the paraphyly of Entognatha. BMC Evol. Biol. 13, 236. (doi:10.1186/1471-2148-13-236) Crossref, PubMed, ISI, Google Scholar - 344
Hunt T 2007 A comprehensive phylogeny of beetles reveals the evolutionary origins of a superradiation. Science 318, 1913–1916. (doi:10.1126/science.1146954) Crossref, PubMed, ISI, Google Scholar - 345
Zhou X et al. 2016 The Trichoptera barcode initiative: a strategy for generating a species-level Tree of Life. Phil. Trans. R. Soc. B 371, 20160025. (doi:10.1098/rstb.2016.0025) Link, ISI, Google Scholar - 346
Simon S, Strauss S, von Haeseler A, Hadrys H . 2009 A phylogenomic approach to resolve the basal pterygote divergence. Mol. Biol. Evol. 26, 2719–2730. (doi:10.1093/molbev/msp191) Crossref, PubMed, ISI, Google Scholar - 347
Regier JC, Shultz JW, Zwick A, Hussey A, Ball B, Wetzer R, Martin JW, Cunningham CW . 2010 Arthropod relationships revealed by phylogenomic analysis of nuclear protein-coding sequences. Nature 463, 1079–1083. (doi:10.1038/nature08742) Crossref, PubMed, ISI, Google Scholar - 348
Thomas JA, Trueman JWH, Rambaut A, Welch JJ . 2013 Relaxed phylogenetics and the palaeoptera problem: resolving deep ancestral splits in the insect phylogeny. Syst. Biol. 62, 285–297. (doi:10.1093/sysbio/sys093) Crossref, PubMed, ISI, Google Scholar - 349
Meusemann K 2010 A phylogenomic approach to resolve the arthropod tree of life. Mol. Biol. Evol. 27, 2451–2464. (doi:10.1093/molbev/msq130) Crossref, PubMed, ISI, Google Scholar - 350
- 351
Simon S, Narechania A, DeSalle R, Hadrys H . 2012 Insect phylogenomics: exploring the source of incongruence using new transcriptomic data. Genome Biol. Evol. 4, 1295–1309. (doi:10.1093/gbe/evs104) Crossref, PubMed, ISI, Google Scholar - 352
Whitfield JB, Kjer KM . 2008 Ancient rapid radiations of insects: challenges for phylogenetic analysis. Annu. Rev. Entomol. 53, 449–472. (doi:10.1146/annurev.ento.53.103106.093304) Crossref, PubMed, ISI, Google Scholar - 353
Lemmon AR, Emme SA, Lemmon EM . 2012 Anchored hybrid enrichment for massively high-throughput phylogenomics. Syst. Biol. 61, 727–744. (doi:10.1093/sysbio/sys049) Crossref, PubMed, ISI, Google Scholar - 354
Faircloth BC, McCormack JE, Crawford NG, Harvey MG, Brumfield RT, Glenn TC . 2012 Ultraconserved elements anchor thousands of genetic markers spanning multiple evolutionary timescales. Syst. Biol. 61, 717–726. (doi:10.1093/sysbio/sys004) Crossref, PubMed, ISI, Google Scholar - 355
Letsch H, Simon S . 2013 Insect phylogenomics: new insights on the relationships of lower neopteran orders (Polyneoptera). Syst. Entomol. 38, 783–793. (doi:10.1111/syen.12028) Crossref, ISI, Google Scholar - 356
Hörnschemeyer T, Beutel RG, Pasop F . 2002 Head structures of Priacma serrata leconte (Coleptera, Archostemata) inferred from X-ray tomography. J. Morphol. 252, 298–314. (doi:10.1002/jmor.1107) Crossref, PubMed, ISI, Google Scholar - 357
Beutel RG, Haas A . 1998 Larval head morphology of Hydroscapha natans LeConte, 1874 (Coleoptera, Myxophaga, Hydroscaphidae) with special reference to miniaturization. Zoomorphology 118, 103–116. (doi:10.1007/s004350050061) Crossref, ISI, Google Scholar - 358
Wirkner CS, Richter S . 2010 Evolutionary morphology of the circulatory system in Peracarida (Malacostraca; Crustacea). Cladistics 26, 143–167. (doi:10.1111/j.1096-0031.2009.00278.x) Crossref, ISI, Google Scholar - 359
Beutel RG, Kristensen NP . 2012 Morphology and insect systematics in the era of phylogenomics. Arthropod Struct. Dev. 41, 303–305. (doi:10.1016/j.asd.2012.05.003) Crossref, PubMed, ISI, Google Scholar - 360
Matsumura Y, Wipfler B, Pohl HW, Dallai R, Machida R, Câmara JT 2015 Cephalic anatomy of Zorotypus weidneri New, 1978: new evidence for a placement of Zoraptera. Arthropod Syst. Phylogeny 3, 85–105. Google Scholar - 361
Holland BR, Spencer HG, Worthy TH, Kennedy M . 2010 Identifying cliques of convergent characters: concerted evolution in the cormorants and shags. Syst. Biol. 59, 433–445. (doi:10.1093/sysbio/syq023) Crossref, PubMed, ISI, Google Scholar - 362
Blanke A, Greve C, Wipfler B, Beutel RG, Holland BR, Misof B . 2013 the identification of concerted convergence in insect heads corroborates Palaeoptera. Syst. Biol. 62, 250–263. (doi:10.1093/sysbio/sys091) Crossref, PubMed, ISI, Google Scholar - 363
Willmann R . 2005 Phylogenese und Systemk der Insecta. In Lehrbuch der Speziellen Zoologie, begründet von A Kaestner, 2 Aufl Bd I: Wirbellose Tiere 5 Teil: Insecta (ed.Dathe HH ), pp. 1–65. Heidelberg, Berlin: Spektrum, Gustav Fischer. Google Scholar - 364
Kraus O, Kraus MK . 1994 Phylogenetic system of the tracheata (Mandibulata): on ‘Myriapoda’—Insecta interrelation-ships, phylogenetic age and primary ecological niches. Verh. Nat. Ver Hambg N.F. 34, 5–31. Google Scholar - 365
Regier JC, Shultz JW . 1997 Molecular phylogeny of the major arthropod groups indicates polyphyly of crustaceans and a new hypothesis for the origin of hexapods. Mol. Biol. Evol. 14, 902–913. (doi:10.1093/oxfordjournals.molbev.a025833) Crossref, PubMed, ISI, Google Scholar - 366
Regier JC, Shultz JW, Kambic RE . 2005 Pancrustacean phylogeny: hexapods are terrestrial crustaceans and maxillopods are not monophyletic. Proc. R. Soc. B 272, 395. (doi:10.1098/rspb.2004.2917) Link, ISI, Google Scholar - 367
Cook CE, Yue Q, Akam M . 2005 Mitochondrial genomes suggest that hexapods and crustaceans are mutually paraphyletic. Proc. R. Soc. B 272, 1295–1304. (doi:10.1098/rspb.2004.3042) Link, ISI, Google Scholar - 368
Cook CE, Smith ML, Telford MJ, Bastianello A, Akam M . 2001 Hox genes and the phylogeny of the arthropods. Curr. Biol. 11, 759–763. (doi:10.1016/S0960-9822(01)00222-6) Crossref, PubMed, ISI, Google Scholar - 369
Richter S . 2002 The Tetraconata concept: hexapod-crustacean relationships and the phylogeny of Crustacea. Org. Divers. Evol. 2, 217–237. (doi:10.1078/1439-6092-00048) Crossref, ISI, Google Scholar - 370
Fanenbruck M, Harzsch S, Wägele JW . 2004 The brain of the Remipedia (Crustacea) and an alternative hypothesis on their phylogenetic relationships. Proc. Natl Acad. Sci. USA 101, 3868–3873. (doi:10.1073/pnas.0306212101) Crossref, PubMed, ISI, Google Scholar - 371
Oakley TH, Wolfe JM, Lindgren AR, Zaharoff AK . 2012 Phylotranscriptomics to bring the understudied into the fold: monophyletic Ostracoda, fossil placement, and pancrustacean phylogeny. Mol. Biol. Evol. 30, 215–233. (doi:10.1093/molbev/mss216) Crossref, PubMed, ISI, Google Scholar - 372
Beutel RG, Friedrich F, Yang X-K, Ge S-Q . 2013 Insect morphology and phylogeny: a textbook for students of entomology. Berlin, Germany: Walter de Gruyter. Crossref, Google Scholar - 373
Blanke A, Koch M, Wipfler B, Wilde F, Misof B . 2014 Head morphology of Tricholepidion gertschi indicates monophyletic Zygentoma. Front. Zool. 11, 16. (doi:10.1186/1742-9994-11-16) Crossref, PubMed, ISI, Google Scholar - 374
Beier M . 1969 5. Klassifikation. In Handbuch der Zoologie IV Band: Arthropoda—2 Hälfte: Insecta, pp. 1–17. Berlin, Germany: Walter De Gruyter. Google Scholar - 375
Hadrys H, Simon S, Kaune B, Schmitt O, Schöner A, Jakob W, Schierwater B . 2012 Isolation of hox cluster genes from insects reveals an accelerated sequence evolution rate. PLoS ONE 7, e34682. (doi:10.1371/journal.pone.0034682) Crossref, PubMed, ISI, Google Scholar - 376
Letsch HO, Meusemann K, Wipfler B, Schütte K, Beutel R, Misof B . 2012 Insect phylogenomics: results, problems and the impact of matrix composition. Proc. R. Soc. B 279, 3282–3290. (doi:10.1098/rspb.2012.0744) Link, ISI, Google Scholar - 377
Savard J, Tautz D, Lercher MJ . 2006 Genome-wide acceleration of protein evolution in flies (Diptera). BMC Evol. Biol. 6, 7. (doi:10.1186/1471-2148-6-7) Crossref, PubMed, ISI, Google Scholar - 378
Staniczek AH . 2000 The mandible of silverfish (Insecta: Zygentoma) and mayflies (Ephemeroptera): its morphology and phylogenetic significance. Zool. Anz. 239, 147–178. ISI, Google Scholar - 379
Bitsch J, Bitsch C, Bourgoin T, D'Haese C . 2004 The phylogenetic position of early hexapod lineages: morphological data contradict molecular data. Syst. Entomol. 29, 433–440. (doi:10.1111/j.0307-6970.2004.00261.x) Crossref, ISI, Google Scholar - 380
Kristensen NP . 1997 The groundplan and basal diversification of the hexapods. In Arthropod Relationships (edsFortey RA, Thomas RH ), pp. 281–293. London, UK: Chapman & Hall. Google Scholar - 381
Bitsch C, Bitsch J . 2000 The phylogenetic interrelationships of the higher taxa of apterygote hexapods. Zool. Scr. 29, 131–156. (doi:10.1046/j.1463-6409.2000.00036.x) Crossref, ISI, Google Scholar - 382
Shao H, Zhang Y, Xie R, Yin W . 1999 Mitochondrial cytochrome b sequences variation of Protura and molecular systematics of Apterygota. Chin. Sci. Bull. 44, 2031–2036. (doi:10.1007/BF02884915) Crossref, Google Scholar - 383
Giribet G, Edgecombe GD, Wheeler WC . 2001 Arthropod phylogeny based on eight molecular loci and morphology. Nature 413, 157–161. (doi:10.1038/35093097) Crossref, PubMed, ISI, Google Scholar - 384
Mallatt JM, Giribet G . 2006 Further use of nearly complete 28S and 18S rRNA genes to classify Ecdysozoa: 37 more arthropods and a kinorhynch. Mol. Phylogenet. Evol. 40, 772–794. (doi:10.1016/j.ympev.2006.04.021) Crossref, PubMed, ISI, Google Scholar - 385
Dell'Ampio E, Szucsich NU, Carapelli A, Frati F, Steiner G, Steinacher A, Pass G . 2009 Testing for misleading effects in the phylogenetic reconstruction of ancient lineages of hexapods: influence of character dependence and character choice in analyses of 28S rRNA sequences. Zool. Scr. 38, 155–170. (doi:10.1111/j.1463-6409.2008.00368.x) Crossref, ISI, Google Scholar - 386
Hovmöller R, Pape T, Kallersjo M . 2002 The Palaeoptera problem: basal pterygote phylogeny inferred from 18S and 28S rDNA sequences. Cladistics 18, 313–323. (doi:10.1006/clad.2002.0199) Crossref, ISI, Google Scholar - 387
Pisani D, Poling LL, Lyons-Weiler M, Hedges B . 2004 The colonization of land by animals: molecular phylogeny and divergence times among arthropods. BMC Biol. 2, 1. (doi:10.1186/1741-7007-2-1) Crossref, PubMed, ISI, Google Scholar - 388
Regier JC, Shultz JW . 2001 Elongation factor-2: a useful gene for arthropod phylogenetics. Mol. Phylogenet. Evol. 20, 136–148. (doi:10.1006/mpev.2001.0956) Crossref, PubMed, ISI, Google Scholar - 389
Rota-Stabelli O, Campbell L, Brinkmann H, Edgecombe GD, Longhorn SJ, Peterson KJ, Pisani D, Philippe H, Telford MJ . 2011 A congruent solution to arthropod phylogeny: phylogenomics, microRNAs and morphology support monophyletic Mandibulata. Proc. R. Soc. B 278, 298–306. (doi:10.1098/rspb.2010.0590) Link, ISI, Google Scholar - 390
Seeger W . 1979 Spezialmerkmale an Eihüllen und Embryonen von Psocoptera im Vergleich zu anderen Paraneoptera (Insecta): Psocoptera als monophyletische Gruppe. Stuttg. Beitr. Naturkunde A 329, 1–57. Google Scholar