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 errors. 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, the quality of conclusions will be determined by an improved understanding of both molecular and morphological evolution, and not simply the number of genes analyzed.

Keywords

Insects, morphology, parsimony, cladistics, likelihood, phylogenomics

Authors’ contributions

Kjer and Beutel wrote about their respective areas of expertise. Simon provided comments, fact-checking, and additional historical and methodological insights. Yavorskaya interpreted and summarized Russian works.

 

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 from human interactions as well. 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. But 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 generated with a few hundred nucleotides were considered revolutionary, just as the application of cladistic1 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 Sanger2 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.

The numerous names of orders and other higher-level taxa1 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 Table 1 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 3 entognathous orders, whereas Insecta excludes the entognathous orders (Table 1). Phylogenetic terms that many readers might not be familiar with are defined in numbered boxes, which are referenced by superscripts at the first several usages of the term, referring to the box containing the definition.

 

Pre-Hennigian concepts in insect taxonomy and phylogeny

The roots of insect systematics go back to the 16th, 17th and 18th 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 18th century, the Swedish botanist Carolus Linnaeus (1707-1778) described more than 10,000 species in his “Systema Naturae” (3), including over 2,000 insects. His ordinal names refer to the characteristics of the wings, e.g., Heteroptera (heterogeneous forewing), Hymenoptera (membranous forewing), 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 for classification based on phylogeny (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 today see that it only partly reflected phylogenetic relationships. However, Haeckel presented the first explicit phylogenetic tree of insects (Haeckel 1896 (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 (Fig. 1) comes close to concepts developed decades later. Naturally, since our current cladistic1 concept of reserving names for monophyletic1 groups (10)(11) was not developed until the 1950s and 60s, 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. (Fig. 1; See table 1 here, and throughout, for a definition of taxon1 names).

Terms 1

Monophyletic group: A group of organisms 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 parsimony2. “Cladistics” can mean two very different things.

Taxon (pl. taxa): A named group of organisms. Under cladistic principles, all taxa are monophyletic.

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).

 

A highly productive North American entomologist of the early 20th century was G.C. Crampton (12)(13), who’s 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 Handlirsch, who was frequently cited in Hennig’s later work (11) (see Friedrich et al. (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).

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 Henning 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 monophyly1, 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 sistergroup1 relationships. Distinguishing ancestral from derived character states required the definition of a taxon1 outside the group of interest for comparison  (an outgroup1). 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 (2). Hennig’s phylogeny (11), published in 1969, is widely considered to be the starting point of modern insect phylogenetics (Fig. 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 (Fig. 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; e.g., (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 monophyletic1 originated with Hennig, but the unique value of synapomorphies was loosely recognized by systematists earlier. Herbert Ross (1908-1978)(26), for instance, was polarizing1 characters relative to a hypothetical ancestor in 1937, and he indicated derived states with marks on the internodes2 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 gotten. 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 (Fig. 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 main focus on genital structures, especially the ovipositor. This included thrips (33), Neuropterida (lace wings and close relatives), beetles (34), and fleas, flies and scorpion flies (35)(36)(37)(38). In 1973 (34), he provided specific evidence for the first time for a close relationship between neuropteroids and beetles.

Geographic isolation and “parallel universes”

The importance of the contributions made by Russian entomologists and insect paleontologists 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 paleontologists, 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.

Through the entire 20th century, most Russian entomologists maintained a conservative approach, with traditional descriptions based on morphology, without formal cladistic1 phylogenetic character evaluations. International collaboration was partly impeded by linguistic barriers, but also by the isolation of the Soviet Union and the pseudoscientific Lysenkoism, an anti-genetic 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 19th 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 20th century. He is best known as the author of the 1905 magisterial “Beetles of Russia, Western Europe and neighboring 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 sistergroup1 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 data sets (47)(48)(49)(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 (44))

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 (Fig. 5) with the largest dataset assembled to date. Europe was an ideal meeting place because the particular brand of cladistic fervor 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 centers of insect phylogenetics based on emerging techniques in morphology, and the refinement of model-based molecular phylogenetics (reviewed in (57)).

 

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)(59)(60)(61), profound reviews of insect phylogeny (62)(63)(64), and landmark volumes on systematics and morphology of Lepidoptera in the Handbook of Zoology series (65)(66)(67)(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 skepticism, 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 skeptical even after the publication of large transcriptomic4 works (49)(50) (pers. comm. to RGB).

Controversy was common in morphology-based insect phylogenetics and results were strongly affected by the selection of characters, before very large and well documented data sets emerged in the 21st century. Boudreaux’s 1979 book (69) on arthropod phylogeny, though criticized by some (63)(70)(71)(43), was cited frequently by others (72)(73)(74)(75)(76)(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 paleontology (80), numerous specific studies on extant and extinct insects (81)(82)(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 (e.g., 82). 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)(36)(37)(62)(63)(96)(97)(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 (e.g. (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 (e.g. (99)) could take years. In the early 2000s, new technologies such as micro-computed tomography (µCT) and computer-based 3D reconstructions greatly accelerated the acquisition of high quality anatomical data (19). The coordinated efforts of international research teams, using both new and traditional methods (e.g. (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.

 

Insect morphology and cladistics

In the late 1970s and 1980s cladistics “evolved” as a transformed version of Hennigian phylogenetic systematics (e.g. (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 outgroups1 (i.e., rooting2) with computer-based analysis (e.g. (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 Henning86 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 an analysis of a morphological matrix-based analysis of most major 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 (e.g.,(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).

 

The dawn of molecular systematics in the early ’90s – 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 parsimony2, distance2, and other methods (109). Their distance analysis grouped the annelids with the mollusks 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 artifacts and alignment 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.

 

Terms 2

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.

Synapomorphy: A shared, derived character (feature) that can used as an argument for a group being monophyletic.

Homoplasy: Character state evolving more than once on a tree or changing back to its original state. Parsimony attempts to minimize homoplasy. Parsimony attempts to minimize homoplasy. Homoplasy creates phylogenetic noise (conflicting 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 classified based on overall similarity in their phenotype or appearance, rather than focussing 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 pseudoreplicates are then analyzed 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. (Jackknifing is similar, but in jackknifing, the new pseudoreplicate 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, jackknifing, 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. Defining an outgroup also establishes a 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, since they have a very low probability of accumulating informative substitions.

Substitution: An observed change in a character. For molecular data, substitutions are related to mutations, but since 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-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.

 

A 1984 review of insect molecular systematics by Stewart 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 Chris 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 18S sequences and restriction sites2. The restriction site analysis (his Fig. 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 substitution2 rate, and long regions that had to be excluded from the analysis because they could not be aligned with confidence. They surmised that the flies (Diptera) were being drawn to the root2, and thus, excluded them in most of their analyses Pashley et al. in 1993 (116) published distance2 and parsimony2 analyses of a fragment of 18S rRNA from nine orders of Holometabola. They were able to recover Mecopterida, and Amphiesmenoptera, but bootstrap2 support for most groups was very low. Pashley et al. concluded that using different outgroups1 yielded different topologies for poorly supported ingroup1 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.

Mitochondrial data were also explored in the early days of Sanger2 sequencing. Liu and Beckenbach (118) explored the mitochondrial cytochrome oxydase II (COII) gene in 10 orders of insects, using a genetic-distance-based2 analysis (119) and parsimony2. 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 analyzed, and it was proposed that onychophorans (velvet worms), are in fact modified arthropods (120). Since 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 Sanger2 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 program to locate and sequence single copy nuclear protein coding genes (127). Of 14 “promising candidates” they identified, several (EF-1a, DDC, POLII, and to a lesser extent, PEPCK) saw extensive use in insect intra-ordinal phylogenetics. They would continue to develop useful protocols for amplifying genes such as wingless, CAD and others (128)(129)(130)(131)(132)(133)(134)(135). Additional contributions to the arsenal of nuclear genes for insect phylogenetics soon followed (136)(137)(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 transcriptome4 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).

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 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 modeling DNA to account for known biases in the way it evolves. Swofford and colleague’s 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 (likelihood2) analyses at a time when computers were finally up to the task of implementing complex substitution2 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 (109). 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.

 

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 likelihood2 and parsimony practitioners would argue, which in turn helped to reveal the problems inherent in parsimony. Most morphological studies placed Strepsiptera (Fig. 4) as the sister taxon to the beetles (Coleoptera), based on hind-wing flight (posteromotorism), and a few other characters (72)(150)(151)(152), or within a subgroup of Coleoptera (153) (see also McKenna and Farrell (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 and 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 hind wing stubs found in flies that are superficially similar to strepsipteran forewings. They surmised that the grouping of flies with strepsiperans was in itself evidence of a homeotic mutational transformation that resulted in halteres flipping from the 3rd thoracic segment in flies to the second in strepsipterans (Fig. 4). Most morphologists doubted this assertion. Skepticism came from the molecular perspective as well. Carmean and Crespi (156) responded almost immediately that “long branches attract flies”, which was unambiguously demonstrated by Huelsenbeck (157), with a likelihood2 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 neighbor-joining distance method with a model (160) designed to accommodate nucleotide compositional bias3 among lineages and among-site rate variation3. 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 bootstrap2 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.

 

Terms 3.

Nucleotide compositional bias: When nucleotide frequencies stray significantly from 25% each. 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 3rd codon positions than for 2nd 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 only showing 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, while 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 specification of a number of parameters. Common input parameters include numbers that indicate costs or ratios, or parameters of a specific evolutionary model. Input parameters are often derived from empirical data and can drastically alter phylogenetic results.

 

 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 alignment3 program Malign (163), and analyzed 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, due to contamination. The paper is best remembered for its recovery of Halteria, a hypothesis that they would vigorously defend (74)(164)(165)(166)(167) until morphology (48)(168) and additional genes (154)(47)(154)(169)(170) overturned it 15 years after it had been proposed.

 

 

Long-Branch distraction?

Whiting (77)(164) rejected the suggestion that Halteria was an artifact of long branch attraction2, 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 internode2 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 and 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 bias3 in the stem dipteran lineage that would pose problems for phylogenetic analyses. These biases were further explored in 2000 by Steel and colleagues (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 attraction2, while 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 transcriptomic4 (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.

 

Alignment Issues

In addition to disagreements over the merits of parsimony, the methods of nucleotide alignment3 played an important role in insect phylogenetics during the Sanger sequencing period (107)(163)(178)(179)(180)(181)(182)(183)(184)(185)(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 centered 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 principle 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 itself, presented in their Figs 12a 13, and 14, suggested many implausible relationships (107). However, it seems that the morphological data as shown in their Fig. 10, provided a reasonable scaffold on which to structure 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 Figs. 18a&b, were largely unresolved consensus trees3. 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 nodes2 the authors favoured from different datasets. The deepest parts of their trees seemed robust to analytical assumptions, while relationships among polyneopterans and holometabolans were unstable. Despite published papers that pointed toward branch effects and compositional bias3 for these data (137)(156)(157)(171)(172)(188)(189), they did not take these problems into consideration, favoring parsimony on philosophical grounds.

Many of the differences among phylogenetic hypotheses were the result of differing analytical approaches and ambiguous alignment3 of rRNA data. The history of alignment disputes have 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 favor with systematists, due 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.

 

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 vs. Cladistics), different heroes (Felsenstein vs. Farris), different branch support2 measures (bootstraps vs. jack-knifing and Bremer support) and even different computer systems (Mac vs. PC) brought about by the platforms of their different programs (PAUP vs. 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 due to tradition and the fact that morphological characters are less amenable to modeling 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). Morphologist 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 parameters3 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 to 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)(198)(199)(200)(201). However, even if ILD tests were legitimate, one must decide which among many parameters should be evaluated, each with an infinite space to explore, with each influencing the behavior of the others (182). Grant and Kluge (202), in a particularly radical application of their own view of epistemological consistency, reject the whole idea of sensitivity analyses3, and suggest that all parameters should be equal, and set to 1 on philosophical grounds. Ogden and 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 and Whiting in 2005 (204), and parsimony analyses were virtually abandoned by the practitioners of likelihood. We see now that short, ancient internodes2 are always sensitive to assumptions and input parameters. Thus, the nodes2 unseen by Börner in 1904 collapse with sensitivity analyses3 as it was applied, as seen, for example in fig. 18 in Wheeler et al. 2001 (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; Table 1)–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) (but see 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 compositional bias, non-homogenous substitution rates, alignment error and the inconsistencies of rRNA analysis with POY (107)(171)(172)(181)(182)(183)(184)(208).

 

The likelihood camp

The basic principles and performance of likelihood (209)(210)(211)(212)(213)(214)(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, 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 analyze more than 10 taxa in a reasonable time frame. Sophisticated models of evolution could not be implemented until computers gradually gained the speed to analyze datasets of typical size, more than ten years after PHYLIP was introduced. Throughout the 1990’s phylogenetic methods based on likelihood as an optimality criterion increased in importance and implementation.

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 and Penny (218) 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 phenetics2. In 2000, it could still take weeks on a desktop computer to analyze 500 nucleotides for 50 taxa. Bootstrapping or any kind of branch support2 was difficult if not impossible for likelihood analyses with more than 25 taxa until fast maximum likelihood programs were developed–PhyML (219), Garli (220) and RAxML (221)(222).

Bayesian analysis, which shares many properties with likelihood, was introduced into phylogenetics in 1999 (223), and could be implemented in the user friendly program MrBayes (224). 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 due to 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 variation3, 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)(226)(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 and Joyce (231)), including a stronger penalty for increasing the number of parameters3.

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 gotten larger. Friedrich and Tautz, in 1995 (232), were among the first to use likelihood to estimate deep arthropod relationships with Phylip (216). Their analysis included 3 hexapods, and recovered Pancrustacea, and crustacean paraphyly. Likelihood (among other methods) was also used by von Dohlen and 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 variation3 in the analysis had more of an effect on likelihood scores than the substitution models themselves. Flook and Rowell (234) used likelihood methods to explore the properties of mitochondrial data among orthopterans. Whitfield and 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 (238) were the first entomologists to use Bayesian methods in their study of Caddisfly (Trichoptera) phylogeny, and Kjer (101) was the first to include many insect orders with Bayesian methods.

It was not until the mid 2000s that consensus among entomologists in the U.S. swung toward likelihood analyses for molecular data, but parsimony analyses are still being published due to 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”).

 

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 due to historical artifact 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, while (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 and 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 nodes2 we now find with larger datasets; many for the first time with molecular data, including Xenonomia, Eukinolabia, and Haplocerata (Table 1), 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-1a, Histone H3, and mitochondrial 12S, 16S, COI and COII, along with 170 morphological characters from the older sources, such as Hennig, and Kristensen. 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 and 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 root2 that rendered beetles paraphyletic. For the first time since 1997 (77) the molecular data recovered Hymenoptera as sister to the rest of Holometabola (Table 1. 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 neighboring 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 neighbors in the tree. Ribosomal RNA analyses frequently recover Nonoculata (74)(75)(107)(117)(239)(242)(243)(244)(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 analyzed with existing methods, even though much larger datasets now contradict Nonoculata and mecopteran paraphyly.

 

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 and colleagues 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-1a. Intron positions in EF2 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-1a (250), found a remarkable tree from only 6 informative characters, but intron positions did show homoplasy2, largely due to 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 and Roth have put their analysis of neuropeptides to use in studies of cockroaches, grasshoppers and Mantophasmatodea (252)(253)(254)(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 nodes2 among insect orders, because it was assumed that gene order was highly conserved and unlikely to be homoplastic2. However, the most controversial internodes2 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)(259)(260)(261) and Hymenoptera (262)(263), with groups supported by changes in gene order, summarized in a recent review by Cameron (190).

Mitochondrial genomes

Mitochondrial data have been the subject of two recent reviews (190)(264). The accumulation of mitochondrial genomes continued through the 2000s (265)(266)(267)(268)(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 throughput4 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)(272)(273)(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 and Hadrys (274) recommended a similarly cautious view, finding many implausible relationships among orders even when using dense taxon sampling, and careful modeling. 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)(278)(279)(280)(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 transcriptomes4 and some of these nodes (Fig. 5: P,Q,R) have only rarely been seen before. Mitochondrial data consistently recover Mantophasmatodea with Phasmatodea (281)(273). Cameron et al. (282) found Megaloptera sister to Neuroptera, reflecting results from transcriptomes (50). The number of nucleotides in any analysis is strongly correlated with branch support2. However, since 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 centered around those who advocate for more data (both in terms of 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.

 

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)(284)(285)(286)(287)(288)(289)(290) (286)(292), Ephemeroptera (293)(294)(278), Plecoptera (295), Dermaptera (296), Embioptera (297), Phasmatodea (298), Dictyoptera (236)(299)(300)(279)(301), Mantodea (302), Orthoptera (234)(303)(304)(305)(306), Hemiptera (307)(308), Psocodea (309), Hymenoptera (310)(311)(312)(313)(314)(315)(316)(317), Neuropterida (318)(319)(320), Coleoptera (321)(322)(322)(323), Diptera (324)(325), Lepidoptera (326)(327)(328)(329)(330)(331) (332)(333)(334), Trichoptera (335)(336)(337), Mecoptera (246)(338)(339) and Siphonaptera (340)(341).

 

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, analyzed 185 nuclear genes from four holometabolous orders, rooted2 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, 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 (Fig. 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. (338) sequenced over 10K nucleotides from these same 3 genes, and focused their attention to 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.

      Data-mining and big, automated phylogeny pipelines – Behaviorists, ecologists, and other biologists rely on phylogenetic trees to understand the evolution of complex characteristics. GenBank is data rich and the temptation to create pipelines4 to download, combine, filter, and analyze these data to produce a tree is strong. Building upon work by Hunt and Vogler (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 a 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 phylogenenetic 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.

The Palaeoptera Problem Revisited. An early transcriptomic analysis of seven pterygote orders, rooted with Collembola, grouped the mayflies with Neoptera (the Chiastomyaria hypothesis; Table 1) (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 masked4 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 due to missing data, resulting in a smaller, but more dense matrix. Meusemann et al. recovered both Palaeoptera and Chiastomyaria using alternative analytical parameters3. In addition, they also found that Hymenoptera were 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 (349) (50), but rarely the morphologically favoured Metapterygota (but see (351)). While Misof et al. (50) reported the monophyly of Palaeoptera, the quartet mapping4 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 (Fig. 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.

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 to previously sequenced genomes of 2 beetles, 4 hymenopterans, 3 flies, and Bombyx (silkmoth), with 2 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 analyzed 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 famiIies, Rhipiphoridae or Meloidae (among others).

Insect Phylogeny Resolved. Many labs are now collecting large datasets from hybrid capture techniques, such as anchored hybrid enrichment (353), or ultraconserved elements (354) (UCEs). 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, Schultz, and Brady (Formicidae) and Kjer/Frandsen (Trichoptera) labs currently have large hybrid enriched datasets in progress, and according to their conference presentations, these data are largely resolving long-standing problems with strong bootstrap2 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 (Fig. 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 assembly4, orthology prediction4, data masking4, and signal optimization. Protein domains were considered as partitions4, and site-specific rate models were developed. Diplura was sister to Insecta, in agreement with Letsch and S. 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.

Terms 4.

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 orthologs.

Orthology prediction: An ortholog is a homologous gene. Gene duplication produces non-homologous, similar-looking genes (paralogs) that can disrupt phylogenetic analyses. Orthology prediction is the process of identifying orthologs, and distinguishing them from paralogs.

Data masking: Identifying “bad” data, and throwing it out. Problematic data comes from a variety of sources. Data can be misaligned, randomized, mostly missing, or it can violate model assumptions.

Partitions: Subsets of the data. Since models are used to analyze 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 4 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.

 

Integrated Phylogenetics

Innovative approaches such as micro-computed tomography (356) and computer-based reconstruction (e.g., (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 the most unusual result from Misof et al. 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 skeptical (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 (360)(79). Despite recent progress, it is obvious that morphology has its limitations. Even large data sets of high quality create partially unsatisfying results, sometimes despite impressive lists of shared derived character states (synapomorphies) for presumptive clades, as in (360). Artifactual 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 and colleagues (2015). As the authors pointed out, some of the obtained presumptive synapomorphies were obviously the result of misleading redundant evolution. Correlated characters can also cause artifacts in morphology-based phylogenetic reconstructions as addressed (361) in the context of the Palaeoptera problem (362)(342). Solutions were suggested, based on modified weighting or the exclusion of characters.

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 either placed as sistergroup of Myriapoda (millipedes and centipedes), or as the sister taxon of a myriapod subgroup (363)(364). Molecular data sets of different size and composition and analyzed with different approaches, consistently yielded a clade Pancrustacea (also called Tetraconata), usually with hexapods placed among paraphyletic crustacean lineages (50)(257)(347)(349)(346)(366)(367)(368). Even though some morphological arguments for Pancrustacea have been presented (369)(370), a formal character analysis is still lacking and the morphological evidence far from convincing. Possible candidates for 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 and 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 3-segments and an 11-segmented abdomen. 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 2nd maxillae (labium), the loss of the ventral food rim, the absence of midgut glands and nephridial organs, and others.

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. But 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 mapping4 (50), may point toward a true case of incomplete lineage sorting. Many of our difficult to resolve nodes may have complex and 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 materials (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 behavior 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 vs. 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 pipelines4 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 Earth.

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, Alex Blanke, Bernhard Misof, and Sabrina Simon, provided helpful comments. John Morse provided valuable insights about his advisor, Herbert Ross. KMK thanks the Schlinger endowment for funding.

Competing interests

We have no competing interests.

 

Literature cited

  1. Beutel RG, Friedrich F, Leschen RA. Charles Darwin, beetles and phylogenetics. Naturwissenschaften. 2009;96(11):1293–312.
  2. Engel MS, Kristensen NP. A History of Entomological Classification. Annu Rev Entomol. 2013; 58(1):585–607.
  3. Linneaus C. Systema naturæ per regna tria naturæ, secundum classes, ordines, genera, species, cum characteribus, differentiis, synonymis, locis. 10th ed. Stockholm: Laurentius Salvius; 1758.
  4. Grimaldi DA, Engel MS. Evolution of the Insects. Cambridge University Press;
  5. Latreille PA. Précis des caractères génériques des insectes, disposés dans un ordre naturel. [Internet]. Paris,: Prèvôt; 1796 [cited 2016 Jan 10]. Available from: http://www.biodiversitylibrary.org/bibliography/58411
  6. Darwin CB, Wallace AR. On the tendency of species to form varieties; and on the perpetuation of varieties and species by natural means of selection. [Read 1 July]. J Proc Linn Soc Lond Zool. 1858;3:45–50.
  7. Darwin C. On the origin of species by means of natural selection, or the preservation of favored races in the struggle for life. London: John Murray;
  8. Haekel E. Phylogeny. 1896. Systematische Phylogenie der wirbellosen Thiere. Reimer, Berlin.
  9. Börner C. Zur Systematik der Hexapoda. Zool Anz. 1904;27:511–33.
  10. Hennig W. Kritische Bemerkungen zum phylogenetischen System der Insecten. Beitr Zur Entomol Sonderh. 1953;3:1–85.
  11. Hennig W. Die Stammesgeschichte der Insekten. Frankfurt am Main: Waldemar Kramer; 1969.
  12. Crampton GC. A phylogenetic study of the terga and wing bases in Embiids, Plecoptera, Dermaptera, and Coleoptera. Psyche (Stuttg). 1918;25:4–13.
  13. Crampton GC. The eulabium, mentum, submentum and gular region of insects. J Entomol Zool. 1928;20:1–18.
  14. Crampton GC. The interrelationships and lines of descent of living insects. Psyche (Stuttg). 1938;45(4):165–81.
  15. Imms AD. The ancestry of insects. Trans Soc Br Entomol. 3:1–32.
  16. Snodgrass RE. The Principles of Insect Morphology. New York, USA: McGraw-Hill Book Company; 1935.
  17. Weber H. Lehrbuch der Entomologie. Jena, Germany: Gustav Fischer; 1933. 1-726 p.
  18. Weber H. Grundriß der Insektenkunde. Jena: Gustav Fischer; 1938. 1225 p.
  19. Friedrich F, Matsumura Y, Pohl H, Bai M, Hörnschemeyer T, Beutel RG. Insect morphology in the age of phylogenomics: innovative techniques and its future role in systematics. Entomol Sci. 2014;17(1):1–24.
  20. Handlirsch A. Neue Untersuchungen über die fossilen Insekten. Ann Naturhistorischen Mus Wien. 1937;40:1–40.
  21. Handlirsch A. ,Neue Untersuchungen über die fossilen Insekten. Ann Naturhistorischen Mus Wien. 1939;49:1–240.
  22. Hennig W. Grundzüge einer Theorie der Phylogenetischen Systematik. Berlin: Deutscher Zentralverlag; 1950.
  23. Watrous LE, Wheeler QD. The Out-Group Comparison Method of Character Analysis. Syst Biol. 1981;30(1):1–11.
  24. Maddison WP, Donoghue MJ, Maddison DR. Outgroup analysis and parsimony. Syst Zool. 1984;33:83–103.
  25. Sturtevant AV. On the subdivision of the genus Drosophila. Proc Natl Acad Sci. 1939;25:137–41.
  26. Ross HH. A generic classification of the Nearctic sawflies (Hymenoptera: Symphyta). Ill Biol Monogr. 1937;15(2):1–173.
  27. Ross HH. The evolution of the insect orders. Entomol News. 1955;66(8):197–208.
  28. Ross HH. A Textbook of Entomology. 3rd ed. New York: J. Wiley; 1965.
  29. Mayr E. Grundlagen der Zoologischen Systematik. Hamburg, Berlin: Verlaug Paul Parey; 1975.
  30. Nelson G. Darwin-Hennig classification: a reply to Ernst Mayr. Syst Zool. 23(3):452–8.
  31. Wille A. The phylogeny and relatonships between the insect orders. Rev Biol Trop. 1960;8:93–123.
  32. Hinton HE. The phylogeny of the panorpoid orders. Annu Rev Entomol. 1958;3(1):181–206.
  33. Mickoleit G. Zur Thoraxmorphologie der Thysanoptera. (Thesis dissertation). Zool Jahrb Abt Anat Ontog Tiere. 1960;79:1–92.
  34. Mickoleit G. Über den Ovipositor der Neuropteroidea und Coleoptera und seine phylogenetische Bedeutung (Insecta, Holometabola). Z Für Morphol Tiere. 1973;74(1):37–64.
  35. Mickoleit G. Die Genital- und Postgenitalsegmente der Mecoptera-Weibchen (Insecta, Holometabola) I Das Exoskelet. Z Für Morphol Tiere. 1975;80:97–135.
  36. Mickoleit G. Die Genital- und Postgenitalsegmente der Mecoptera-Weibchen (Insecta, Holometabola) II Das Dach der Genitalkammer. Zeitschrift für Morphologie der Tiere. Z Für Morphol Tiere. 1976;88:133–56.
  37. Mickoleit G. Die phylogenetischen Beziehungen der Schnabelfliegen-Familien aufgrund morphologischer Ausprägungen der weiblichen Genital- und Postgenitalsegmente (Mecoptera). Entomol Ger. 1978;4:258–71.
  38. Mickoleit G. Die Sperma-Auspreßvorrichtung der Nannochoristidae (Insecta: Mecoptera). Entomol Gen. 2008;31:193–226.
  39. Martynov AV. Über zwei Grundtypen der Flügel bei den Insecten und ihre Evolution. Zoomorphology. 1925;4(3):465–501.
  40. Martynov AV. Studies on the geological history and phylogeny of the orders of insects (Pterygota). Tr Paleont Inst Inst Akad Nauk SSSR. 1938;7:1–150.
  41. Ponomarenko AG. The historical development of archoste­matan beetles. Tr Paleont Inst Akad Nauk SSSR. 1969;125:1–238.
  42. Ponomarenko AG. Mesozoic Coleo­ptera Arnoldy LV, Jerikin VV, Nikri­tin LM, Ponomarenko AG (eds) Suborder Adephaga, etc. pp. 3-104. Tr Palaeont Inst Akad Nauk SSSR. 1977;161:1–183.
  43. Rasnitsyn AP, Quicke DLJ. Rasnitsyn, A.P., Quicke, D.L.J. (eds.) History of Insects. Dordrecht: Kluwer Academic Publishers; 2002.
  44. Hoßfeld U, Olsson L. From the Modern Synthesis to Lysenkoism, and Back? Science. 2002;297(5578):55.
  45. Motschulsky VI. Die Käfer Russlands. I. Insecta Carabica. Moscow: Gautier; 1850.
  46. Jacobson GG. The beetles of Russia, West Europe and adjacent countries. Devrien; 1905.
  47. Wiegmann BM, Trautwein MD, Kim J-W, Cassel BK, Bertone MA, Winterton SL, et al. Single-copy nuclear genes resolve the phylogeny of the holometabolous insects. BMC Biol. 2009;7:34.
  48. Beutel RG, Friedrich F, Hörnschemeyer T, Pohl H, Hünefeld F, Beckmann F, et al. Morphological and molecular evidence converge upon a robust phylogeny of the megadiverse Holometabola. Cladistics. 2011;27(4):341–55.
  49. Peters RS, Meusemann K, Petersen M, Mayer C, Wilbrandt J, Ziesmann T, et al. The evolutionary history of holometabolous insects inferred from transcriptome-based phylogeny and comprehensive morphological data. BMC Evol Biol. 2014;14:52.
  50. Misof B, Liu S, Meusemann K, Peters RS, Donath A, Mayer C, et al. Phylogenomics resolves the timing and pattern of insect evolution. Science. 2014;346(6210):763–7.
  51. Ronquist F, Rasnitsyn AP, Roy A, Eriksson K, Lindgren M. Phylogeny of the Hymenoptera: A cladistic reanalysis of Rasnitsyn’s (1988) data. Zool Scr. 1999;28(1-2):13–50.
  52. Rasnitsyn AP. On the skimming hypothesis of the origin of insect flight. Acta Zool Cracoviensia. 2003;46:85–8.
  53. Brothers DJ. Alexandr Pavlovich Rasnitsyn, (palaeo)entomologist extraordinaire – a personal appreciation. ZooKeys. 2011;(130):1–10.
  54. Jeannel R. Ordre des Coleopteroides. In: Traite de Zoologie. P.-P.Grasse; 1949. p. 771–891.
  55. Hennig W. Phylogenetic Systematics. Urbana, Illinois, USA: University of Illinois Press; 1966.
  56. Klass K-D, editor. Proceedings of the 1st Dresden meeting on insect phylogeny: “Phylogenetic Relationships within the Insect Orders” (Dresden, September 19–21, 2003). Entomol Abh. 2003;61:119–72.
  57. Felsenstein J. Inferring Phylogenies. Sunderland, MA.: Sinauer Associates Inc.; 2004.
  58. Kristensen NP. The skeletal anatomy of the heads of adult Mnesarchaeidae and Neopseustidae (Lep., Dacnonypha). Ent Meddr. 1968;36:137–51.
  59. Kristensen NP. The Anatomy of the Head and the Alimentary Canal of Adult Eriocraniidae (Lep., Dacnonypha), in Entomologiske Meddelelser. In Copenaghen: Entomologisk Forening; 1968. p. 239–315.
  60. Kristensen NP. The Morphological and Functional Evolution of the Mouthparts in Adult Lepidoptera (abstract). In: Opuscula Entomologica. Lund: Entomologiska sällskapet; 1968. p. 69–72.
  61. Kristensen NP. Studies on the morphology and systematics of primitive Lepidoptera. Steenstrupia. 1984;10:141–91.
  62. Kristensen NP. The phylogeny of hexapod “orders”. A critical review of recent accounts. Z Zool Syst Evol-Forsch. 1975;13:1–44.
  63. Kristensen NP. Phylogeny of Insect Orders. Annu Rev Entomol. 1981;26(1):135–57.
  64. Kristensen NP. Phylogeny of extant hexapods. In: The Insects of Australia. Melbourne: Melbourne University Publishing; 1991. p. 126–40.
  65. Kristensen NP. The Non-Glossatan Moths. Lepidoptera: Moths and Butterflies. 1. Evolution, Systematics, and Biogeography. In: Handbook of Zoology Vol IV, Part 35. Berlin and New York: De Gruyter; 1999. p. 41–9.
  66. Kristensen NP. The Homoneurous Glossata. In: Lepidoptera: Moths and Butterflies 1 Evolution, Systematics, and Biogeography Handbook of Zoology Vol IV, Part 35. Berlin and New York: De Gruyter; 1999. p. 51–63.
  67. Kristensen NP, Skalski AW. Phylogeny and paleontology. In: Lepidoptera: Moths and Butterflies 1 Evolution, Systematics, and Biogeography Handbook of Zoology Vol IV, Part 35. Berlin and New York: De Gruyter; 1999. p. 7–25.
  68. Kristensen NP. Reproductive organs, Lepidoptera:Moths and butterflies 2. In: Handbook of Zoology IV/36. 2003. p. 427–47.
  69. Boudreaux HB. Arthropod Phylogeny with Special Reference to the Insects. New York, Chichester, Brisbane, Toronto: Wiley & Sons Inc; 1979.
  70. Mashimo Y, Matsumura Y, Machida R, Dallai R, Gottardo M, Yoshizawa K, et al. 100 years Zoraptera – a phantom in insect evolution and the history of its investigation. Insect Systematics & Evolution. J Morphol. 2014;45(4):371–93.
  71. Wipfler B, Bai M, Schoville S, Dallai R, Uchifune T, Machida R, et al. Ice Crawlers (Grylloblattodea) – the history of the investigation of a highly unusual group of insects. J Insect Biodivers. 2014;2(2):1–25.
  72. Beutel RG, Gorb SN. Ultrastructure of attachment specializations of hexapods (Arthropoda): evolutionary patterns inferred from a revised ordinal phylogeny. J Zool Syst Evol Res. 2001;39(4):177–207.
  73. Klass K-D, Kristensen NP. The ground plan and afÞnities of hexapods: recent progress and open problems. Ann Soc Entomol Fr. 2001;37:265–98.
  74. Wheeler WC, Whiting M, Wheeler QD, Carpenter JM. The Phylogeny of the Extant Hexapod Orders. Cladistics. 2001;17(2):113–69.
  75. Kazunori Yoshizawa KPJ. Yoshizawa, K. and K. P. Johnson. Aligned 18S for Zoraptera (Insecta): phylogenetic position and molecular evolution. Molecular Phylogenetics and Evolution 37:2. Mol Phylogenet Evol. 2005;37(2):572–80.
  76. Wheeler WC, Cartwright P, Hayashi CY. Arthropod Phylogeny: A combined approach. Cladistics. 1993;9(1):1–39.
  77. Whiting MF, Carpenter JC, Wheeler QD, Wheeler WC. The Strepsiptera Problem: Phylogeny of the Holometabolous Insect Orders Inferred from 18S and 28S Ribosomal DNA Sequences and Morphology. Syst Biol. 1997;46(1):1–68.
  78. Kristensen NP. Forty years‘ insect phylogenetic systematics. Zool. Beitr. N.F. Zool Beitre NF. 1995;36(1):83–124.
  79. Mashimo Y, Beutel RG, Dallai R, Lee C-Y, Machida R. Embryonic development of Zoraptera with special reference to external morphology, and its phylogenetic implications (Insecta). J Morphol. 2014;275(3):295–312.
  80. Kukalova-Peck J. Fossil history and the evolution of hexapod structures. Ch. 6. In: The Insects of Australia A textbook for students and research workers Ed CSIRO. Ithaca, New York, USA: Cornell University Press; 1991. p. 141–79.
  81. Kukalová-Peck J, Brauckmann C. Wing folding in pterygote insects, and the oldest Diaphanopterodea from the early Late Carboniferous of West Germany. Can J Zool. 1990;68:1104–11.
  82. Kukalova‐Peck J, Peck SB. Zoraptera wing structures: evidence for new genera and relationship with the blattoid orders (Insecta: Blattoneoptera). Syst Entomol. 1993 Oct 1;18(4):333–50.
  83. Kukalová-Peck J, Beutel RG. Is the Carboniferous †Adiphlebia lacoana really the. EJE. 2013;109(4):633–45.
  84. Haas F, Kukalova-Peck J. Dermaptera hindwing structure and folding: New evidence for familial, ordinal and superordinal relationships within Neoptera (Insecta). Eur J Entomol. 2001;98(4):445–509.
  85. Kukalova-Peck J, Lawrence JF. Relationships among coleopteran suborders and major endoneopteran lineages: Evidence from hind wing characters. Eur J Entomol. 2004;101(1):95–144.
  86. Kukalová-Peck J. Origin of the insect wing and wing articulation from the arthropodan leg. Can J Zool. 1983;61:1618–69.
  87. Kukalová-Peck J. 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. 1987 Oct 1;65(10):2327–45.
  88. Kukalová-Peck J. Phylogeny of Higher Taxa in Insecta: Finding Synapomorphies in the Extant Fauna and Separating Them from Homoplasies. Evol Biol. 2008;35(1):4–51.
  89. Béthoux O, Kristensen NP, Engel M S. Hennigian phylogenetic systematics and the “groundplan” vs. “post-groundplan” approaches: A reply to Kukalová-Peck. Evol Biol. 35(4):317–23.
  90. Pass G. Gross and fine structure of the antennal circulatory organ in cockroaches (Blattodea, Insecta). J Morph. 1985;185:255–68.
  91. Pass G. Accessory Pulsatile Organs: Evolutionary Innovations in Insects. Annu Rev Entomol. 2000;45(1):495–518.
  92. Klass K-D. Die Phylogeny der Dictyoptera. Ph.D. thesis. Munich: Fakulta’t fur Biologie. Ludwig Maximilians Universita’t.; 1995.
  93. Klass K-D, et al. Place saver. 2003;
  94. Machida R, Nagashima T, Ando H. The early embryonic development of the jumping bristletail Pedetontus unimaculatus machida (Hexapoda: Microcoryphia, Machilidae). J Morphol. 1990;206(2):181–95.
  95. Machida R, Tojo T, Tsutsumi T, Uchifune T, Klass K-D, Picker MD, et al. Embryonic development of heel-walkers: Reference to some prerevolutionary stages (Insecta: Mantophasmatodea). Proc Arthropod Embryol Soc Jpn. 2004;39:31–9.
  96. Klass K-D. The proventriculus of the Dicondylia, with comments on evolution and phylogeny of Dictyoptera and Odonata (Insecta). Zool Anz. 1998;(237):15–42.
  97. Klass K-D. Possible homologies in the proventriculi of Dicondylia (Hexapoda) and Malacostraca (Crustacea). Zool Anz. 1998;(237):43–58.
  98. Klass K-D. The ovipositor of Dictyoptera (Insecta): homology and ground-plan of the main elements. Zool Anz. 1998;(236):69–101.
  99. Mickoleit G. Zur Thoraxmorphologie der Thysanoptera. Zool Jahrb Abt Anat Ontog Tiere. 1961;79:1–92.
  100. Wipfler B, Machida R, Müller B, Beutel RG. 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. 2011;36(2):241–66.
  101. Blanke A, Wipfler B, Letsch H, Koch M, Beckmann F, Beutel R, et al. Revival of Palaeoptera—head characters support a monophyletic origin of Odonata and Ephemeroptera (Insecta). Cladistics. 2012;28(6):560–81.
  102. Patterson C. Verifiability in systematics. Syst Zool. 1978;27:218–22.
  103. Patterson C. Cladistics. Biologist. 1980;27(5):234–40.
  104. Forey PL, Humphries CJ, Kitching IJ, Scotland RW, Siebert DJ, Williams DM. Cladistics: a Practical Course in Systematics. The Systematics Assn. Pubn. No. 10. Oxford, UK: Clarendon Press; 1992.
  105. Farris JS. Hennig86: a PC-DOS program for phylogenetic analysis. Cladistics. 1989;5:163.
  106. Swofford DL. PAUP: Phylogenetic Analysis Using Parsimony, Version 3.1 Computer program distributed by the Illinois Natural History Survey, Champaign, Illinois. 1991.
  107. Kjer KM. Aligned 18S and insect phylogeny. Syst Biol. 2004;53(3):506–14.
  108. Field KG, Olsen GJ, Lane DJ, Giovannoni SJ, Ghiselin MT, Raff EC, et al. Molecular phylogeny of the animal kingdom. Science. 1988;239(4841 Pt 1):748–53.
  109. Lake JA. Origin of the Metazoa. Proc Natl Acad Sci. 1990;87:763–6.
  110. Turbeville JM, Pfeifer DM, Field KG, Raff RA. The phylogenetic status of arthropods, as inferred from 18S rRNA sequences. Mol Biol Evol. 1991;8(5):669–86.
  111. Berlocher SH. Insect Molecular Systematics. Annu Rev Entomol. 1984;29(1):403–33.
  112. Fujiwara H, Ishikawa H. 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. 1982;10.
  113. Simon C. Evolution of 13- and 17-Year Periodical Cicadas (Homoptera: Cicadidae: Magicicada). Bull Entomol Soc Am. 1988;34(4):163–76.
  114. Wheeler WC. The systematics of insect ribosomal DNA. In: The hierarchy of life Molecules and morphology in phylogenetic analysis (K B B Fernholm, and H Jörnvall, ed). Amsterdam: Elsevier; 1989. p. 307–21.
  115. Carmean D, Kimsey LS, Berbee ML. 18S rDNA sequences and the holometabolous insects. Mol Phylogenet Evol. 1992;1:270–8.
  116. Pashley DP, McPheron BA, Zimmer EA. Systematics of the holometabolous insect orders based on 18S ribosomal RNA. Mol Phylogenet Evol. 1993;2:132–42.
  117. Reumont BM von, Meusemann K, Szucsich NU, Dell’Ampio E, Gowri-Shankar V, Bartel D, et al. Can comprehensive background knowledge be incorporated into substitution models to improve phylogenetic analyses? A case study on major arthropod relationships. BMC Evol Biol. 2009;9:119.
  118. Liu H, Beckenbach AT. Evolution of the mitochondrial cytochrome oxidase II gene among 10 orders of insects. Mol Phylogenet Evol. 1992;1(1):41–52.
  119. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25.
  120. Ballard JW, Olsen GJ, Faith DP, Odgers WA, Rowell DM, Atkinson PW. Evidence from 12S ribosomal RNA sequences that onychophorans are modified arthropods. Science. 1992;258(5086):1345–8.
  121. Kjer KM. 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. 1995;4(3):314–30.
  122. Hickson RE, Simon C, Cooper A, Spicer GS, Sullivan J, Penny D. Conserved sequence motifs, alignment, and secondary structure for the third domain of animal 12S rRNA. Mol Biol Evol. 1996;13:150–69.
  123. Simon C, Frati F, Beckenbach AT, Crespi, Liu H, Flook PK. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am. 1994;87:651–701.
  124. Kocher TD, Thomas WK, Meyer A, Edwards SV, Pääbo S, Villablanca FX, et al. Dynamics of mitochondrial DNA evolution in animals: amplification and sequencing with conserved primers. Proc Natl Acad Sci U S A. 1989;86(16):6196–200.
  125. Simon C, Pääbo S, Kocher TD, Wilson AC. Evolution of the mitochondrial ribosomal RNA in insects as shown by the poiymerase chain reaction, In M. Clegg & S. O’Brien [eds.]. In: Molecular evolution UCLA Symposia on Molecular and Cellular Bioiogy, New Series. Liss, New York; 1990. p. 235–44.
  126. Caterino MS, Cho S, Sperling FA. The current state of insect molecular systematics: a thriving Tower of Babel. Annu Rev Entomol. 2000;45:1–54.
  127. Friedlander TP, Regier JC, Mitter C. Nuclear gene sequences for higher level phylogenetic analysis: 14 promising candidates. Syst Biol. 1992;41:483–90.
  128. Friedlander TP, Regier JC, Mitter C. Phylogenetic information content of five nuclear gene sequences in animals: initial assessment of character sets from concordance and divergence studies. Syst Biol. 1994;43(4):511–25.
  129. Cho S, Mitchell A, Regier JC, Mitter C, Poole RW, Friedlander TP, et al. A highly conserved nuclear gene for low-level phylogenetics: elongation factor-1 alpha recovers morphology-based tree for heliothine moths. Mol Biol Evol. 1995;12(4):650–6.
  130. Friedlander TP, Regier JC, Mitter C, Wagner DL. A nuclear gene for higher level phylogenetics: phosphoenolpyruvate carboxykinase tracks mesozoic-age divergences within Lepidoptera (Insecta). Mol Biol Evol. 1996;13(4):594–604.
  131. Friedlander TP, Regier JC, Mitter C. Initial assessment of character sets from five nuclear gene sequences in animals. Pp. 301-320. In: Biodiversity II, ML Reaka-Kudla, DE Wilson & EO Wilson Eds. Washington, DC.: Joseph Henry Press; 1997. p. 301–20.
  132. Friedlander TP, Regier JC, Mitter C, Wagner DL, Fang QQ. Evolution of heteroneuran Lepidoptera (Insecta) and the utility of dopa decarboxylase for Cretaceous-age phylogenetics. Zool J Linn Soc. 2000;130:213–34.
  133. Wiegmann BM, Regier JC, Mitter C, Friedlander TP, Wagner DL, Nielson ES. Nuclear genes resolve Mesozoic-aged divergences in the insect order Lepidoptera. Mol Phylogenet Evol. 2000;15:242–59.
  134. Fang QQ, Cho S, Regier JC, Mitter C, Matthews M, Poole RW, et al. A new nuclear gene for insect phylogenetics: dopa decarboxylase is informative of relationships within Heliothinae (Lepidoptera: Noctuidae). Syst Biol. 1997;46(2):269–83.
  135. Fang QQ, Mitchell A, Regier JC, Mitter C, Friedlander TP, Poole RW. Phylogenetic utility of the nuclear gene dopa decarboxylase in noctuoid moths (Insecta: Lepidoptera: noctuoidea). Mol Phylogenet Evol. 2000;15(3):473–86.
  136. Burmester T, Massey HC, Zakharkin SO, Benes H. The evolution of hexamerins and the phylogeny of insects. J Mol Evol. 1998;47(1):93–108.
  137. Campbell DL, Brower AV, Pierce NE. Molecular evolution of the wingless gene and its implications for the phylogenetic placement of the butterfly family Riodinidae (Lepidoptera: papilionoidea). Mol Biol Evol. 2000;17(5):684–96.
  138. Ascher JS, Danforth BN, Ji S. Phylogenetic utility of the major opsin in bees (Hymenoptera: Apoidea): a reassessment. Mol Phylogenet Evol. 2001;19(1):76–93.
  139. Colgan DJ, McLauchlan A, Wilson GDF, Livingston SP, Edgecombe GD, Macaranas J, et al. Histone H3 and U2 snRNA DNA sequences and arthropod molecular evolution. Aust J Zool. 1998;46(5):419.
  140. Edgecombe GD, Wilson GDF, Colgan DJ, Gray MR, Cassis G. Arthropod cladistics: Combined analysis of histone H3 and U2 snRNA sequences and morphology. Cladistics. 2000;16(2):155–203.
  141. Kjer KM, Carle FL, Litman J, Ware J. A molecular phylogeny of Hexapoda. Arthropod Syst Phylogeny. 2006;64(1):35–44.
  142. Regier JC, Shultz JW, Kambic RE. Phylogeny of basal hexapod lineages and estimates of divergence times. Ann Entomol Soc Am. 2004;97:411–9.
  143. Saiki RK, Scharf S, Faloona F, Mullis KB, Horn GT, Erlich HA, et al. Enzymatic amplification of beta-globin genomic sequences and restriction site analysis for diagnosis of sickle cell anemia. Science. 1985;230(4732):1350–4.
  144. Mullis KB, Faloona FA. Specific synthesis of DNA in vitro via a polymerase-catalyzed chain reaction. Methods Enzymol. 1987;155:335–50.
  145. Simon C. 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, ed GM Hewitt, AWB Johnston, JPW Young,. New York: Springer Verlag, NATO Advanced Studies Institute; 1991. p. H57:33–71.
  146. Swofford DL, Olsen GJ. Phylogeny reconstruction. In “Molecular Systematics”.(Eds DM Hillis and C. Moritz.) pp. 411–501. Sinauer Associates: Sunderland, Massachusetts; 1990.
  147. Swofford DL, Olsen GJ, Waddell PJ, Hillis DM. Phylogenetic inference. In: Molecular Systematics, ed DM Hillis, C Moritz, BK Mable. 2nd ed. Sunderland, MA: Sinauer Assoc; 1996. p. 407–514.
  148. Swofford DL. PAUP*. Phylogenetic Analysis Using Parsimony (*and Other
 Methods). Version 4. Sunderland, MA.: Sinauer Associates;
  149. Yang Z. Phylogenetic Analysis by Maximum Likelihood (PAML). Department of Integrative Biology, University of California at Berkeley.; 1996.
  150. Kinzelbach R. The systematic position of Strepsiptera (Insecta). Am Entomol. 1990;36:292–303.
  151. Kathirithamby J. Review of the Order Strepsiptera. Syst Entomol. 1989;14(1):41–92.
  152. Beutel RG, Gorb SN. A revised interpretation of the evolution of attachment structures in Hexapoda with special emphasis on Mantophasmatodea. Arthropod Syst Phylogeny. 2006;64(1):3–25.
  153. Crowson RA. The natural classification of the families of Coleoptera. London: Nathaniel Lloyd & Co., Ltd.; 1955. 187 p.
  154. McKenna DD, Farrell BD. 9-Genes Reinforce the Phylogeny of Holometabola and Yield Alternate Views on the Phylogenetic Placement of Strepsiptera. PLoS ONE. 2010;5(7):e11887.
  155. Whiting MF, Wheeler WC. Insect homeotic transformation. Nature. 1994;(368):696.
  156. Carmean D, Crespi BJ. Do long branches attract flies? Nature. 1995;373:666–70.
  157. Huelsenbeck JP. Is the Felsenstein zone a fly trap? Syst Biol. 1997;46(1):69–74.
  158. Chalwatzis N, Bauer A, Stetzer R, Kinzelbach R, Zimmermann RK. Strongly expanded 18S ribosomal-RNA genes correlated with a peculiar morphology in the insect order of Strepsiptera. Zool Anal Complex Syst. 1995;98:115–26.
  159. Chalwatzis N, Hauf J, van de Peer Y, Kinzelbach R, Zimmermann RK. 18S ribosomal RNA genes in insects: primary structure of the genes and molecular phylogeny of the Holometabola. Ann Entomol Soc Am. 1996;89(6):788–803.
  160. Lockhart PJ, Steel MA, Hendy MD, Penny D. Recovering evolutionary trees under a more realistic model of sequence evolution. Mol Biol Evol. 1994;11(4):605–12.
  161. van de Peer Y, Neefs J-M, De Rijk P, de Wachter R. Reconstructing evolution from eukaryotic small-ribosomal-subunit RNA sequences: Calibration of the molecular clock. J Mol Evol. 1993;37:221–32.
  162. Van de Peer Y, Nicolaï S, De Rijk P, De Wachter R. Database on the structure of small ribosomal subunit RNA. Nucleic Acids Res. 1996;24(1):86–91.
  163. Wheeler WC, Gladstein DL. A heuristic multiple sequence alignment procedure based on parsimony. New York, USA: American Museum of Natural History;
  164. Whiting MF. Long-branch distraction and the Strepsiptera. Syst Biol. 1998;47(1):134–8.
  165. Siddall ME, Whiting MF. Long‐Branch Abstractions. Cladistics. 1999;15(1):9–24.
  166. Whiting MF. Phylogeny of the holometabolous insect orders: molecular evidence. Zool Scr. 2002;31(1):3–15.
  167. Whiting M. Phylogenetic position of Diptera: Review of the evidence. In: The Evolutionary Biology of Flies Edited by Yeates DK, Wiegmann BM. New York, USA: Columbia University Press; 2005. p. 3–13.
  168. Friedrich F, Beutel RG. Goodbye Halteria? The thoracic morphology of Endopterygota (Insecta) and its phylogenetic implications. Cladistics. 2010;26(6):579–612.
  169. Longhorn SJ, Pohl HW, Vogler AP. Ribosomal protein genes of holometabolan insects reject the Halteria, instead revealing a close affinity of Strepsiptera with Coleoptera. Mol Phylogenet Evol. 2010;55(3):846–59.
  170. Ishiwata K, Sasaki G, Ogawa J, Miyata T, Su Z-H. Phylogenetic relationships among insect orders based on three nuclear protein-coding gene sequences. Mol. Phylogenet. Mol Phylogenet Evol. 2011;58:169–80.
  171. Huelsenbeck JP. Systematic bias in phylogenetic analysis: is the Strepsiptera problem solved? Syst Biol. 1998;47:519–37.
  172. Steel M, Huson D, Lockhart PJ. Invariable sites models and their use in phylogeny reconstruction. Syst Biol. 2000;49(2):225–32.
  173. Friedrich M, Tautz D. An episodic change of rDNA nucleotide substitution rate has occurred during the emergence of the insect order Diptera. Mol Biol Evol. 1997;14(6):644–53.
  174. McMahon DP, Hayward A, Kathirithamby J. The First Molecular Phylogeny of Strepsiptera (Insecta) Reveals an Early Burst of Molecular Evolution Correlated with the Transition to Endoparasitism. Escriva H, editor. PLoS ONE. 2011;6(6):e21206.
  175. Niehuis O, Hartig G, Grath S, Pohl H, Lehmann J, Tafer H, et al. Genomic and morphological evidence converge to resolve the enigma of Strepsiptera. Curr Biol CB. 2012;22(14):1309–13.
  176. McKenna DD. Molecular Phylogenetics and Evolution of Coleoptera. In: Handbook of Zoology Volume IV Arthropoda: Insecta Part 38 Coleoptera, Beetles Volume 3: Morphology and Systematics (Phytophaga) Beutel, R G and R A B Leschen (Eds). Berlin: Walter de Gruyter; 2014. p. 1–10.
  177. Boussau B, Walton Z, Delgado JA, Collantes F, Beani L, Stewart IJ, et al. Strepsiptera, Phylogenomics and the Long Branch Attraction Problem. PLoS ONE. 2014 Oct 1;9(10):e107709.
  178. Wheeler WC. Sequence alignment, parameter sensitivity, and the phylogenetic analysis of molecular data. Syst Biol. 1995;44(3):321–33.
  179. Wheeler WC. Optimization alignment: The end of multiple sequence alignment in phylogenetics? Cladistics. 1996;12:1–9.
  180. Hickson RE, Simon C, Perrey SW. The Performance of Several Multiple-Sequence Alignment Programs in Relation to Secondary-Structure Features for an rRNA Sequence. Mol Biol Evol. 2000;17(4):530–9.
  181. Kjer KM, Gillespie JJ, Ober KA. Structural homology in ribosomal RNA, and a deliberation on POY. Arthropod Syst Phylogeny. 2006;64:71–6.
  182. Kjer KM, Gillespie JJ, Ober KA. Opinions on multiple sequence alignment, and an empirical comparison of repeatability and accuracy between POY and structural alignment. Syst Biol. 2007;56(1):133–46.
  183. Ogden TH, Rosenberg MS. Multiple Sequence Alignment Accuracy and Phylogenetic Inference. Syst Biol. 2006;55(2):314–28.
  184. Ogden TH, Rosenberg MS. Alignment and topological accuracy of the direct optimization approach via POY and traditional phylogenetics via ClustalW + PAUP*. Syst Biol. 2007;56:182–93.
  185. Kjer KM, Roshan U, Gillespie JJ. Structural and evolutionary considerations for multiple sequence alignment of RNA, and the challenges for algorithms that ignore them. 2009 [cited 2015 Dec 26]; Available from: http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.122.8235
  186. Stocsits RR, Letsch H, Hertel J, Misof B, Stadler PF. Accurate and efficient reconstruction of deep phylogenies from structured RNAs. Nucleic Acids Res. 2009 Oct 1;37(18):6184–93.
  187. Sankoff D, Cedergren RJ. Simultaneous Comparison of three or more sequences related by a tree In D. Sankotr & J. B. Kruskal [eds.]. In: Time warps, string edits, and macromolecules: the theory and practice of sequence comparison. Reading, Pennsylvania: Addison-Wesley; 1983. p. 253–63.
  188. Felsenstein J. Cases in which Parsimony or Compatibility Methods Will be Positively Misleading. Syst Zool. 1978;27(4):401–10.
  189. Hwang U-W, Kim W, Tautz D, Friedrich M. Molecular phylogenetics at the Felsenstein Zone: Approaching the Strepsiptera problem using 5.85 and 28S rDNA sequences. Mol Phylogenet Evol. 1998;9:470–80.
  190. Cameron SL. Insect Mitochondrial Genomics: Implications for Evolution and Phylogeny. Annu Rev Entomol. 2014;59(1):95–117.
  191. Letsch HO, Kuck P, Stocsits RR, Misof B. 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. 2010;27(11):2507–21.
  192. Redelings B, Suchard M. Joint Bayesian estimation of alignment and phylogeny. Syst Biol. 2005;54(3):401–18.
  193. Simon C, Buckley TR, Frati F, Stewart JB, Beckenbach AT. 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. 2006;37:547–79.
  194. Kjer KM, Swigonova Z, LaPolla JS, Broughton RE. Why weight? Mol Phylogenet Evol. 2007;43(3):999–1004.
  195. Phillips A, Janies D, Wheeler W. Multiple sequence alignment in phylogenetic analysis. Mol Phylogenet Evol. 2000;16(3):317–30.
  196. Mickevich MF, Farris JS. The implications of congruence in Menidia. Syst Zool. 1981;30:351–70.
  197. Graham SW, Kohn JR, Morton BR, Eckenwalder E, Barrett CH. Phylogenetic congruence and discordance among one morphological and three molecular data sets from Pontederiaceae. Syst Biol. 1998;47:545–67.
  198. Dolphin K, Belshaw R, Orme CDL, Quicke DLJ. Noise and Incongruence: Interpreting Results of the Incongruence Length Difference Test. Mol Phylogenet Evol. 2000;17(3):401–6.
  199. Dowton M, Austin AD. Increased incongruence does not necessarily indicate increased phylogenetic accuracy-The behavior of the ILD test in mixed-model analyses. Syst Biol. 2002;51:19–31.
  200. Darlu P, Lecointre G. When Does the Incongruence Length Difference Test Fail? Mol Biol Evol. 2002;19(4):432–7.
  201. Barker FK, Lutzoni F. The utility of the incongruence length difference test. Syst Biol. 51:625–37.
  202. Grant T, Kluge AG. Data exploration in phylogenetic inference: scientific, heuristic, or neither. Cladistics. 2003 Oct 1;19(5):379–418.
  203. Ogden TH, Whiting MF. The problem with “the Palaeoptera Problem,” sense and sensitivity. Cladistics. 2003;19:432–42.
  204. Terry MD, Whiting MF. Mantophasmatodea and phylogeny of the lower neopterous insects. Cladistics. 2005;21(3):240–57.
  205. McKittrick FA. A contribution to the understanding of cockroach-termite affinities. Ann Entomol Soc Am. 1965;58(1):18–22.
  206. Dell’Ampio E, Meusemann K, Szucsich NU, Peters RS, Meyer B, Borner J, et al. Decisive Data Sets in Phylogenomics: Lessons from Studies on the Phylogenetic Relationships of Primarily Wingless Insects. Mol Biol Evol. 2014;31(1):239–49.
  207. Editors the. Editorial. Cladistics. 2016;32(1):1.
  208. Yoshizawa K. Direct optimization overly optimizes data. Syst Entomol. 2010;35:199–206.
  209. Felsenstein J. Confidence Limits on Phylogenies: An Approach Using the Bootstrap. Evolution. 1985;39(4):783.
  210. Felsenstein J. Phylogenies and quantitative characters. Annu Rev Ecol Syst. 1988;445–71.
  211. Felsenstein J. Phylogenies from molecular sequences: inference and reliability. Annu Rev Genet. 1988;22(1):521–65.
  212. Felsenstein J. Numerical methods for inferring evolutionary trees. Q Rev Biol. 1982;379–404.
  213. Felsenstein J. A likelihood approach to character weighting and what it tells us about parsimony and compatibility. Biol J Linn Soc. 1981;16(3):183–96.
  214. Felsenstein J. Phylogenies from restriction sites: A maximum-likelihood approach. Evolution. 1992;159–73.
  215. Felsenstein J, Sober E. Parsimony and likelihood: an exchange. Syst Biol. 1986;35(4):617–26.
  216. Felsenstein J. PHYLIP-phylogeny interference package (version 3.2). Cladistics. 1989;5:164–6.
  217. Penny D, Steel MA. Parsimony, likelihood, and the role of models in molecular phylogenetics. Mol Biol Evol. 2000;17:839–50.
  218. Hendy MD, Penny D. A framework for the quantitative study of evolutionary trees. Syst Zool. 1989;38:297–309.
  219. Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003 Oct;52(5):696–704.
  220. Zwickl DJ. Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. Ph.D. dissertation,. Austin, TX: The University of Texas at Austin;
  221. Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinforma Oxf Engl. 2006;22(21):2688–90.
  222. Stamatakis A, Hoover P, Rougemont J. A rapid bootstrap algorithm for the RAxML Web servers. Syst Biol. 2008;57:758–71.
  223. Larget B, Simon DL. Markov Chasin Monte Carlo Algorithms for the Bayesian Analysis of Phylogenetic Trees. Mol Biol Evol. 1999;16(6):750.
  224. Swofford DL, Waddell PJ, Huelsenbeck JP, Foster PG, Lewis PO, Rogers JS. Bias in phylogenetic estimation and its relevance to the choice between parsimony and likelihood methods. Syst Biol. 2001;50(4):525–39.
  225. Yang Z. Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: Approximate methods. J Mol Evol. 1994;39(3):306–14.
  226. Sullivan J, Holsinger KE, Simon C. The effect of topology on estimates of among-site rate variation. J Mol Evol. 1996;42(2):308–12.
  227. Yang Z. Among-site rate variation and its impact on phylogenetic analyses. Trends Ecol Evol. 1996;11(9):367–72.
  228. Frati F, Simon C, Sullivan J, Swofford DL. Evolution of the mitochondrial cytochrome oxidase II gene in Collembola. J Mol Evol. 1997;44(2):145–58.
  229. Sullivan J, Holsinger KE, Simon C. Among-site rate variation and phylogenetic analysis of 12S rRNA in Sigmodontine rodents. Mol Biol Evol. 1995;12:988–1001.
  230. Posada D, Crandall KA. MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998;14(9):817–8.
  231. Sullivan J, Joyce P. Model Selection in Phylogenetics. Annu Rev Ecol Evol Syst. 2005;36(1):445–66.
  232. Friedrich M, Tautz D. Ribosomal DNA phylogeny of the major extant arthropod classes and the evolution of myriapods. Nature. 1995;376(6536):165–7.
  233. von Dohlen CD, Moran NA. Molecular phylogeny of the Homoptera: a paraphyletic taxon. J Mol Evol. 1995;41(2):211–23.
  234. Flook PK, Rowell CHF. The effectiveness of mitochondrial rRNA sequences for the reconstruction of the phylogeny of an insect order (Orthoptera). Mol Phylogenet Evol. 1997;8:177–92.
  235. Whitfield JB, Cameron SA. Hierarchical analysis of variation in the mitochondrial 16S rRNA gene among Hymenoptera. Mol Biol Evol. 1998;15(12):1728–43.
  236. Lo N, Tokuda G, Watanabe H, Rose H, Slaytor M, Maekawa K, et al. Evidence from multiple gene sequences indicates that termites evolved from wood-feeding cockroaches. Curr Biol. 2000;10(13):801–4.
  237. Hull D. Science as a Process: An Evolutionary Account of the Social and Conceptual Development of Science. In: Appleman P, editor. Darwin: A Norton Critical Edition, 3rd edition. Chicago: University of Chicago Press; 1988. p. 361–3.
  238. Schöniger M, Von Haeseler A. A stochastic model for the evolution of autocorrelated DNA sequences. Mol Phylogenet Evol. 1994;3(3):240–7.
  239. Misof B, Niehuis O, Bischoff I, Rickert A, Erpenbeck D, Staniczek A. Towards an 18S phylogeny of hexapods: accounting for group-specific character covariance in optimized mixed nucleotide/doublet models. Zool Jena Ger. 2007;110(5):409–29.
  240. Misof B, Misof K. A Monte Carlo Approach Successfully Identifies Randomness in Multiple Sequence Alignments : A More Objective Means of Data Exclusion. Syst Biol. 2009;58(1):21–34.
  241. Xie Q, Tian X, Qin Y, Bu W. 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. 2009;50(2):310–6.
  242. Luan Y, Mallatt JM, Xie R, Yang Y, Yin W. The Phylogenetic Positions of Three Basal-Hexapod Groups (Protura, Diplura, and Collembola) Based on Ribosomal RNA Gene Sequences. Mol Biol Evol. 2005;22(7):1579–92.
  243. Giribet G, Edgecombe GD, Carpenter JM, D’Haese CA, Wheeler WC. Is Ellipura monophyletic? A combined analysis of basal hexapod relationships with emphasis on the origin of insects. Org Divers Evol. 2004;4(4):319–40.
  244. Gao Y, Bu Y, Luan Y. Phylogenetic relationships of basal hexapods reconstructed from nearly complete 18S and 28S rRNA gene sequences. Zool Sci. 2008;25:1139–45.
  245. Mallatt J, Craig CW, Yoder MJ. 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. Mol Phylogenet Evol. 2009;55(1):1–17.
  246. Whiting MF. Mecoptera is paraphyletic: multiple genes and phylogeny of Mecoptera and Siphonaptera. Zool Scr. 2002;31(1):93–104.
  247. Rokas A, Kathirithamby J, Holland PWH. Intron insertion as a phylogenetic character: the engrailed homeobox of Strepsiptera does not indicate affinity with Diptera. Insect Mol Biol. 1999;8(4):527–30.
  248. Carapelli A, Frati F, Nardi F, Dallai R, Simon C. Molecular phylogeny of the apterygotan insects based on nuclear and mitochondrial genes. Pedobiologia. 2000;44(3–4):361–73.
  249. Krauss V. Phylogenetic Mapping of Intron Positions: A Case Study of Translation Initiation Factor eIF2. Mol Biol Evol. 2004;22(1):74–84.
  250. Djernæs M, Damgaard J. Exon-intron structure, paralogy and sequenced regions of elongation factor-1 alpha in Hexapoda. Arthropod Syst Phylogeny. 2006;64:45–52.
  251. Bonneton F, Brunet FG, Kathirithamby J, Laudet V. The rapid divergence of the ecdysone receptor is a synapomorphy for Mecopterida that clarifies the Strepsiptera problem. Insect Mol Biol. 2006;15(3):351–62.
  252. Predel R, Roth S. Neuropeptide evolution and the analysis of phylogenetic relationships in Blattaria (Hexapoda). Arthropod Syst Phylogeny. 2007;65:3–6.
  253. Roth S, Köhler G, Reinhardt K, Predel R. A discrete neuropeptide difference between two hybridizing grasshopper subspecies. Biol J Linn Soc. 2007;91(4):541–8.
  254. Roth S, Fromm B, Gäde G, Predel R. A proteomic approach for studying insect phylogeny: CAPA peptides of ancient insect taxa (Dictyoptera, Blattoptera) as a test case. BMC Evol Biol. 2009;9:50.
  255. Predel R, Neupert S, Huetteroth W, Kahnt J, Waidelich D, Roth S. Peptidomics-based phylogeny and biogeography of Mantophasmatodea (Hexapoda). Syst Biol. 2012;61(4):609–29.
  256. Boore JL, Collins TM, Stanton D, Daehler LL, Brown WM. Deducing the pattern of arthropod phytogeny from mitochondrial DNA rearrangements. Nature. 1995;376(6536):163–5.
  257. Boore JL, Lavrov DV, Brown WM. Gene translocation links insects and crustaceans. Nature. 1998;392(6677):667–8.
  258. Shao R, Campbell NJH, Barker SC. Numerous Gene Rearrangements in the Mitochondrial Genome of the Wallaby Louse, Heterodoxus macropus (Phthiraptera). Mol Biol Evol. 2001;18(5):858–65.
  259. Shao R, Barker SC. 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. 2003;20(3):362–70.
  260. Thao ML, Baumann L, Baumann P. Organization of the mitochondrial genomes of whiteflies, aphids, and psyllids (Hemiptera, Sternorrhyncha). BMC Evol Biol. 2004;4:25.
  261. Wei D-D, Shao R, Yuan M-L, Dou W, Barker SC, Wang J-J. The Multipartite Mitochondrial Genome of Liposcelis bostrychophila : Insights into the Evolution of Mitochondrial Genomes in Bilateral Animals. PLOS ONE. 2012;7(3):e33973.
  262. Dowton M, Cameron SL, Austin AD, Whiting MF. 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. 2009;52(2):512–9.
  263. Dowton M, Cameron SL, Dowavic JI, Austin AD, Whiting MF. Characterization of 67 Mitochondrial tRNA Gene Rearrangements in the Hymenoptera Suggests That Mitochondrial tRNA Gene Position Is Selectively Neutral. Mol Biol Evol. 2009;26(7):1607–17.
  264. Mandal SD, Chhakchhuak L, Gurusubramanian G, Kumar NS. Mitochondrial markers for identification and phylogenetic studies in insects – A Review. DNA Barcodes. 2014;2:1.
  265. Nardi F, Carapelli A, Fanciulli PP, Dallai R, Frati F. The complete mitochondrial DNA sequence of the basal hexapod Tetrodontophora bielanensis: evidence for heteroplasmy and tRNA translocations. Mol Biol Evol. 2001;18(7):1293–304.
  266. Friedrich M, Muqim N. Sequence and phylogenetic analysis of the complete mitochondrial genome of the flour beetle Tribolium castanaeum. Mol Phylogenet Evol. 2003;26(3):502–12.
  267. Bae JS, Kim I, Sohn HD, Jin BR. The mitochondrial genome of the firefly, Pyrocoelia rufa: complete DNA sequence, genome organization, and phylogenetic analysis with other insects. Mol Phylogenet Evol. 2004;32(3):978–85.
  268. Kim I, Cha SY, Yoon MH, Hwang JS, Lee SM, Sohn HD, et al. The complete nucleotide sequence and gene organization of the mitochondrial genome of the oriental mole cricket, Gryllotalpa orientalis (Orthoptera: Gryllotalpidae). Gene. 2005;353(2):155–68.
  269. Stewart JB, Beckenbach AT. Insect mitochondrial genomics: the complete mitochondrial genome sequence of the meadow spittlebug Philaenus spumarius (Hemiptera: Auchenorrhyncha: Cercopoidae). Genome Natl Res Counc Can Génome Cons Natl Rech Can. 2005;48(1):46–54.
  270. Liu S, Wang X, Xie L, Tan M, Li Z, Su X, et al. Mitochondrial capture enriches mito‐DNA 100 fold, enabling PCR‐free mitogenomics biodiversity analysis. Mol Ecol Resour. 2016;16(2):470–9.
  271. Cameron SL, Miller KB, D’Haese CA, Whiting MF, Barker SC. Mitochondrial genome data alone are not enough to unambiguously resolve the relationships of Entognatha, Insecta and Crustacea sensu lato (Arthropoda). Cladistics. 2004;20(6):534–57.
  272. Cameron SL, Beckenbach AT, Downton MA, Whiting MF. Evidence from Mitochondrial Genomics on Interordinal Relationships in Insects. Arthropod Syst Phylogeny. 2006;64(1):27–34.
  273. Cameron SL, Barker SC, Whiting MF. Mitochondrial genomics and the new insect order Mantophasmatodea. Mol Phylogenet Evol. 2006;38(1):274–9.
  274. Simon S, Hadrys H. A comparative analysis of complete mitochondrial genomes among Hexapoda. Mol Phylogenet Evol. 2013;69(2):393–403.
  275. Talavera G, Vila R. What is the phylogenetic signal limit from mitogenomes? The reconciliation between mitochondrial and nuclear data in the Insecta class phylogeny. BMC Evol Biol. 2011;11:315.
  276. Chen W-J, Koch M, Mallatt JM, Luan Y-X. Comparative Analysis of Mitochondrial Genomes in Diplura (Hexapoda, Arthropoda): Taxon Sampling Is Crucial for Phylogenetic Inferences. Genome Biol Evol. 2014;6(1):105–20.
  277. Carapelli A, Vannini L, Nardi F, Boore JL, Beani L, Dallai R, et al. The mitochondrial genome of the entomophagous endoparasite Xenos vesparum (Insecta: Strepsiptera). Gene. 2006;376:248–59.
  278. Zhang J, Zhou C, Gai Y, Song D, Zhou K. The complete mitochondrial genome of Parafronurus youi (Insecta: Ephemeroptera) and phylogenetic position of the Ephemeroptera. Gene. 2008;424(1–2):18–24.
  279. Zhang Y, Xuan W, Zhao J, Zhu C, Jiang G. The complete mitochondrial genome of the cockroach Eupolyphaga sinensis (Blattaria: Polyphagidae) and the phylogenetic relationships within the Dictyoptera. Mol Biol Rep. 2010;37(7):3509–16.
  280. Plazzi F, Ricci A, Passamonti M. The mitochondrial genome of Bacillus stick insects (Phasmatodea) and the phylogeny of orthopteroid insects. Mol Phylogenet Evol. 2011;58(2):304–16.
  281. Wan X, Kim MI, Kim MJ, Kim I. Complete mitochondrial genome of the free-living earwig, Challia fletcheri (Dermaptera: Pygidicranidae) and phylogeny of Polyneoptera. PloS One. 2012;7(8):e42056.
  282. Cameron SL, Sullivan J, Song H, Miller KB, Whiting MF. A mitochondrial genome phylogeny of the Neuropterida (lace-wings, alderflies and snakeflies) and their relationship to the other holometabolous insect orders. Zool Scr. 2009;38(6):575–90.
  283. Rehn AC. Phylogenetic analysis of higher-level relationships of Odonata. Syst Entomol. 2003;28:181–239.
  284. Ware J, May M, Kjer K. Phylogeny of the higher Libelluloidea (Anisoptera: Odonata): An exploration of the most speciose superfamily of dragonflies. Mol Phylogenet Evol. 2007;45(1):289–310.
  285. Bybee SM, Ogden TH, Branham MA, Whiting MF. Molecules, morphology and fossils: a comprehensive approach to odonate phylogeny and the evolution of the odonate wing. Cladistics. 2008;24(4):477–514.
  286. Carle FL, Kjer KM, May ML. Evolution of Odonata, with special reference to Coenagrionoidea (Zygoptera). Arthropod Syst Phylogeny. 2008;66(1):37–44.
  287. Fleck G, Brenk M, Misof B. 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. 2008;8(1):1–16.
  288. Dumont HJ, Vierstraete A, Vanfleteren JR. A molecular phylogeny of the Odonata (Insecta). Syst Entomol. 2010;35(1):6–18.
  289. Davis RB, Nicholson DB, Saunders EL, Mayhew PJ. Fossil gaps inferred from phylogenies alter the apparent nature of diversification in dragonflies and their relatives. BMC Evol Biol. 2011;11(1):252.
  290. Blanke A, Greve C, Mokso R, Beckman F, Misof B. An updated phylogeny of Anisoptera including formal convergence analysis of morphological characters. Syst Entomol. 38:474–90.
  291. Dijkstra K-DB, Kalkman VJ, Dow RA, Stokvis FR, Van Tol J. Redefining the damselfly families: a comprehensive molecular phylogeny of Zygoptera (Odonata): Molecular phylogeny of damselflies. Syst Entomol. 2014;39(1):68–96.
  292. Carle FL, Kjer KM, May ML. A molecular phylogeny and classification of Anisoptera (Odonata). Arthropod Syst Phylogeny. 2015;73(2):281–301.
  293. Ogden TH, Gattolliat JL, Sartori M, Staniczek AH, Soldán T, Whiting MF. Towards a new paradigm in mayfly phylogeny (Ephemeroptera): combined analysis of morphological and molecular data. Syst Entomol. 2009;34(4):616–34.
  294. Sivaramakrishnan KG, Subramanian KA, Arunachalam M, Kumar CS, Sundar S. Emerging trends in molecular systematics and molecular phylogeny of mayflies (Insecta: Ephemeroptera). J Threat Taxa. 2011;3(8):1975–80.
  295. Yu-Han Q, Hai-Yan W, Xiao-Yu J, Wei-Wei Y, Yu-Zhou D. Mitochondrial Genome of the Stonefly Kamimuria wangi (Plecoptera: Perlidae) and Phylogenetic Position of Plecoptera Based on Mitogenomes. PLoS ONE. 2014;9(1):e86328.
  296. Kocarek P, John V, Hulva P. When the Body Hides the Ancestry: Phylogeny of Morphologically Modified Epizoic Earwigs Based on Molecular Evidence. PLoS ONE. 2013;8(6):e66900.
  297. Miller KB, Hayashi C, Whiting MF, Svenson GJ, Edgerly JS. The phylogeny and classification of Embioptera (Insecta). Syst Entomol. 2012;37(3):550–70.
  298. Bradler S, Robertson JA, Whiting MF. A molecular phylogeny of Phasmatodea with emphasis on Necrosciinae, the most species-rich subfamily of stick insects. Syst Entomol. 2014;39(2):205–22.
  299. Ware JL, Litman J, Klass K-D, Spearman LA. Relationships among the major lineages of Dictyoptera: the effect of outgroup selection on dictyopteran tree topology. Syst Entomol. 2008;33(3):429–50.
  300. Legendre F, Nel A, Svenson GJ, Robillard T, Pellens R, Grandcolas P. Phylogeny of Dictyoptera: Dating the Origin of Cockroaches, Praying Mantises and Termites with Molecular Data and Controlled Fossil Evidence. PLoS One. 2015;10(7):e0130127.
  301. Inward D, Beccaloni G, Eggleton P. Death of an order: a comprehensive molecular phylogenetic study confirms that termites are eusocial cockroaches. Biol Lett. 2007;3(3):331.
  302. Svenson GJ, Whiting MF. Reconstructing the origins of praying mantises (Dictyoptera, Mantodea): the roles of Gondwanan vicariance and morphological convergence. Cladistics. 2009;25(5):468–514.
  303. Flook PK, Klee S, Rowell CHF. Combined molecular phylogenetic analysis of the Orthoptera (Arthropoda, Insecta) and implications for their higher systematics. Syst Biol. 1999;48:233–53.
  304. Zhang H-L, Huang Y, Lin L-L, Wang X-Y, Zheng Z-M. The phylogeny of the Orthoptera (Insecta) as deduced from mitogenomic gene sequences. Zool Stud. 2013;52(1):1–13.
  305. Song H. Grasshopper Systematics: Past, Present and Future. J Orthoptera Res. 2010;19(1):57–68.
  306. Song H. 300 million years of diversification: elucidating the patterns of orthopteran evolution based on comprehensive taxon and gene sampling. Cladistics. 2015;
  307. Song N, Liang A-P, Bu C-P. A Molecular Phylogeny of Hemiptera Inferred from Mitochondrial Genome Sequences. PLoS ONE. 2012;7(11):e48778.
  308. Cryan JR, Urban JM. Higher-level phylogeny of the insect order Hemiptera: is Auchenorrhyncha really paraphyletic? Syst Entomol. 2012;37(1):7–21.
  309. Johnson KP, Yoshizawa K, Smith VS. Multiple origins of parasitism in lice. Proc R Soc Lond. 2004;271:1771–6.
  310. Cameron SA, Mardulyn P. Multiple Molecular Data Sets Suggest Independent Origins of Highly Eusocial Behavior in Bees (Hymenoptera:Apinae). Syst Biol. 2001;50(2):194–214.
  311. Hines HM, Hunt JH, O’Connor TK, Gillespie JJ, Cameron SA. Multigene phylogeny reveals eusociality evolved twice in vespid wasps. Proc Natl Acad Sci U S A. 2007;104(9):3295–9.
  312. Peters RS, Meyer B, Krogmann L, Borner J, Meusemann K, Schütte K, et al. The taming of an impossible child: a standardized all-in approach to the phylogeny of Hymenoptera using public database sequences. BMC Biol. 2011;9:55.
  313. Sharkey MJ, Carpenter JM, Vilhelmsen L, Heraty J, Liljeblad J, Dowling APG, et al. Sharkey, M.J., et al. Phylogenetic relationships among superfamilies of Hymenoptera. Cladistics. 2011;28(1):80–112.
  314. Klopfstein S, Vilhelmsen L, Heraty JM, Sharkey M, Ronquist F. The Hymenopteran Tree of Life: Evidence from Protein-Coding Genes and Objectively Aligned Ribosomal Data. PLoS ONE. 2013;8(8):e69344.
  315. Johnson BR, Borowiec ML, Chiu JC, Lee EK, Atallah J, Ward PS. Phylogenomics resolves evolutionary relationships among ants, bees, and wasps. Curr Biol. 2013;23:2058–62.
  316. Mao M, Gibson T, Dowton M. Evolutionary Dynamics of the Mitochondrial Genome in the Evaniomorpha (Hymenoptera)—A Group with an Intermediate Rate of Gene Rearrangement. Genome Biol Evol. 2014;6(7):1862–74.
  317. Ward PS. The Phylogeny and Evolution of Ants. Annu Rev Ecol Evol Syst. 2014;45(1):23–43.
  318. Winterton SL, Hardy NB, Wiegmann BM. On wings of lace. Syst Entomol. 2010;35(3):49–378.
  319. Aspöck U, Haring E, Aspöck H. The phylogeny of the Neuropterida: long lasting and current controversies and challenges (Insecta: Endopterygota). Arthropod Syst Phylogeny. 2012;70:119–29.
  320. Aspöck U. Phylogeny of the Neuropterida (Insecta: Holometabola). Zool Scr. 2002;31(1):51–5.
  321. Mckenna DD, Wild AL, Kanda K, Bellamy CL, Beutel RG, Caterino MS, et al. The beetle tree of life reveals that Coleoptera survived end‐Permian mass extinction to diversify during the Cretaceous terrestrial revolution. Syst Entomol. 2015;40(4):835–80.
  322. Bocak L, Barton C, Crampton-Platt A, Chesters D, Ahrens D, Vogler AP. Building the Coleoptera tree-of-life for >8000 species: composition of public DNA data and fit with Linnaean classification. Syst Entomol. 2014;39(1):97–110.
  323. Lawrence JF, Ślipiński A, Seago AE, Thayer MK, Newton AF, Marvaldi AE. Phylogeny of the Coleoptera Based on Morphological Characters of Adults and Larvae. Ann Zool. 2011;61(1):1–217.
  324. Wiegmann BM, Yeates DK, Thorne JL, Kishino H. Time flies, a new molecular time-scale for brachyceran fly evolution without a clock. Syst Biol. 2003;52(6):745–56.
  325. Wiegmann BM, Trautwein MD, Winkler IS, Barr NB, Kim J-W, Lambkin C, et al. Episodic radiations in the fly tree of life. Proc Natl Acad Sci. 2011;108(14):5690–5.
  326. Wiegmann BM, Regier JC, Mitter C. Combined molecular and morphological evidence on the phylogeny of the earliest lepidopteran lineages. Zool Scr. 2002;31(1):67–81.
  327. Wahlberg N, Braby MF, Brower AV., de Jong R, Lee M-M, Nylin S, et al. Synergistic effects of combining morphological and molecular data in resolving the phylogeny of butterflies and skippers. Proc R Soc B Biol Sci. 2005;272(1572):1577–86.
  328. Wahlberg N, Wheat CW. Genomic Outposts Serve the Phylogenomic Pioneers: Designing Novel Nuclear Markers for Genomic DNA Extractions of Lepidoptera. Syst Biol. 2008;57(2):231–42.
  329. Timmermans MJTN, Lees DC, Simonsen TJ. Towards a mitogenomic phylogeny of Lepidoptera. Mol Phylogenet Evol. 2014;79:169–78.
  330. Regier JC, Mitter C, Zwick A, Bazinet AL, Cummings MP, Kawahara AY, et al. A large-scale, higher-level, molecular phylogenetic study of the insect order Lepidoptera (moths and butterflies). PLoS One. 2013;8(3):e58568.
  331. Regier JC, Mitter C, Kristensen NP, Davis DR, Van Nieukerken EJ, Rota J, et al. A molecular phylogeny for the oldest (nonditrysian) lineages of extant Lepidoptera, with implications for classification, comparative morphology and life-history evolution. Syst Entomol. 2015;40(4):671–704.
  332. Kawahara AY, Breinholt JW. Phylogenomics provides strong evidence for relationships of butterflies and moths. Proc R Soc Lond B Biol Sci. 2014;281(1788):20140970.
  333. Heikkilä M, Mutanen M, Wahlberg N, Sihvonen P, Kaila L. Elusive ditrysian phylogeny: an account of combining systematized morphology with molecular data (Lepidoptera). BMC Evol Biol. 2015;15(1):260.
  334. Simonsen TJ, de Jong R, Heikkilä M, Kaila L. Butterfly morphology in a molecular age – Does it still matter in butterfly systematics? Arthropod Struct Dev. 2012;41(4):307–22.
  335. Kjer KM, Blahnik RJ, Holzenthal RW. Phylogeny of Trichoptera (caddisflies): characterization of signal and noise within multiple datasets. Syst Biol. 2001;50(6):781–816.
  336. Holzenthal RW, Blahnik RJ, Kjer KM, Prather AL. An update on the phylogeny of Caddisflies (Trichoptera). Proceedings of the XIIth International Symposium on Trichoptera. In: Proceedings of the XIIth International Symposium on Trichoptera. The Caddis Press, Bueno-Soria, R. Barba-Alvearz and B. Armitage (Eds).; 2007.
  337. Malm T, Johanson KA, Wahlberg N. The evolutionary history of Trichoptera (Insecta): A case of successful adaptation to life in freshwater. Syst Entomol. 2013;38(3):459–73.
  338. Hu G-L, Yan G, Xu H, Hua B-Z. Molecular phylogeny of Panorpidae (Insecta: Mecoptera) based on mitochondrial and nuclear genes. Mol Phylogenet Evol. 2015;85:22–31.
  339. Whiting MF, Whiting AS, Hastriter MW. A comprehensive phylogeny of Mecoptera and Siphonaptera. Entomol Abh. 2003;61:169.
  340. Whiting MF, Whiting AS, Hastriter MW, Dittmar K. A molecular phylogeny of fleas (Insecta: Siphonaptera): origins and host associations. Cladistics. 2008;24(5):677–707.
  341. Zhu Q, Hastriter MW, Whiting MF, Dittmar K. Fleas (Siphonaptera) are Cretaceous, and evolved with Theria. Mol Phylogenet Evol. 2015;90:129–39.
  342. Savard J, Tautz D, Richards S, Weinstock GM, Gibbs RA, Werren JH, et al. Phylogenomic analysis reveals bees and wasps (Hymenoptera) at the base of the radiation of Holometabolous insects. Genome Res. 2006;16(11):1334–8.
  343. Sasaki G, Ishiwata K, Machida R, Miyata T, Su Z-H. Molecular phylogenetic analyses support the monophyly of Hexapoda and suggest the paraphyly of Entognatha. BMC Evol Biol. 2013;13:236.
  344. Hunt T, Bergsten J, Levkanicova Z, Papadopoulou A, John OS, Wild R, et al. A Comprehensive Phylogeny of Beetles Reveals the Evolutionary Origins of a Superradiation. Science. 2008;318(5858):1913–6.
  345. Zhou X, Frandsen PB, Holzenthal RW, Kjer KM. The Trichoptera barcode initiative: A strategy for generating a species-level tree of life. Trans R Soc -B. 2016; in press
  346. Simon S, Strauss S, von Haeseler A, Hadrys H. A Phylogenomic Approach to Resolve the Basal Pterygote Divergence. Mol Biol Evol. 2009;26(12):2719–30.
  347. Regier JC, Shultz JW, Zwick A, Hussey A, Ball B, Wetzer R, et al. Arthropod relationships revealed by phylogenomic analysis of nuclear protein-coding sequences. Nature. 201;463(7284):1079–83.
  348. Thomas JA, Trueman JWH, Rambaut A, Welch JJ. Relaxed phylogenetics and the palaeoptera problem: resolving deep ancestral splits in the insect phylogeny. Syst Biol. 2013;62(2):285–97.
  349. Meusemann K, Reumont BM von, Simon S, Roeding F, Strauss S, Kück P, et al. A Phylogenomic Approach to Resolve the Arthropod Tree of Life. Mol Biol Evol. 2010;27(11):2451–64.
  350. Meyer B, Misof B. MARE v0.1-rc. 2010.
  351. Simon S, Narechania A, DeSalle R, Hadrys H. Insect Phylogenomics: Exploring the Source of Incongruence Using New Transcriptomic Data. Genome Biol Evol. 2012;4(12):1295–309.
  352. Whitfield JB, Kjer KM. Ancient rapid radiations of insects: challenges for phylogenetic analysis. Annu Rev Entomol. 2008;53:449–72.
  353. Lemmon AR, Emme SA, Lemmon EM. Anchored hybrid enrichment for massively high-throughput phylogenomics. Syst Biol. 2012;61(5):727–44.
  354. Faircloth BC, McCormack JE, Crawford NG, Harvey MG, Brumfield RT, Glenn TC. Ultraconserved elements anchor thousands of genetic markers spanning multiple evolutionary timescales. Syst Biol. 2012;61(5):717–26.
  355. Letsch H, Simon S. Insect phylogenomics: new insights on the relationships of lower neopteran orders (Polyneoptera). Syst Entomol. 2013;38(4):783–93.
  356. Hörnschemeyer T, Beutel RG, Pasop F. Head structures of Priacma serrata leconte (coleptera, archostemata) inferred from X-ray tomography. J Morphol. 2002;252(3):298–314.
  357. Beutel RG, Haas A. Larval head morphology of Hydroscapha natans LeConte, 1874 (Coleoptera, Myxophaga, Hydroscaphidae) with special reference to miniaturization. Zoomorphology. 1998;118:103–16.
  358. Wirkner CS, Richter S. Evolutionary morphology of the circulatory system in Peracarida (Malacostraca; Crustacea). Cladistics. 2010;26(2):143–67.
  359. Beutel RG, Kristensen NP. Morphology and insect systematics in the era of phylogenomics. Arthropod Struct Dev. 2012;41(4):303–5.
  360. Matsumura Y, Wipfler B, Pohl HW, Dallai R, Machida R, Câmara JT, et al. Cephalic anatomy of Zorotypus weidneri New, 1978: new evidence for a placement of Zoraptera. Arthropod Syst Phylogeny. 2015;3(1):85–105.
  361. Holland BR, Spencer HG, Worthy TH, Kennedy M. Identifying Cliques of Convergent Characters: Concerted Evolution in the Cormorants and Shags. Syst Biol. 2010;59(4):433–45.
  362. Blanke A, Greve C, Wipfler B, Beutel RG, Holland BR, Misof B. The Identification of Concerted Convergence in Insect Heads Corroborates Palaeoptera. Syst Biol. 2013;62(2):250–63.
  363. Willmann R. Phylogenese und Systemk der Insecta. In: Dathe, HH (ed) Lehrbuch der Speziellen Zoologie, begründet von A Kaestner, 2 Aufl Bd I: Wirbellose Tiere 5 Teil: Insecta. Heidelberg, Berlin: Spektrum, Gustav Fischer; 2005. p. 1–65.
  364. Kraus O, Kraus MK. Phylogenetic system of the tracheata (Mandibulata): on “Myriapoda” – Insecta interrelation- ships, phylogenetic age and primary ecological niches. Verhandlungen Naturwissenschaftlichen Ver Hambg NF. 1994;34:5–31.
  365. Regier JC, Shultz JW. Molecular phylogeny of the major arthropod groups indicates polyphyly of crustaceans and a new hypothesis for the origin of hexapods. Mol Biol Evol. 1997;14(9):902–13.
  366. Regier JC, Shultz JW, Kambic RE. Pancrustacean phylogeny: hexapods are terrestrial crustaceans and maxillopods are not monophyletic. Proc R Soc B Biol Sci. 2005;272(1561):395.
  367. Cook CE, Yue Q, Akam M. Mitochondrial genomes suggest that hexapods and crustaceans are mutually paraphyletic. Proc R Soc B Biol Sci. 2005;272(1569):1295–304.
  368. Cook CE, Smith ML, Telford MJ, Bastianello A, Akam M. Hox genes and the phylogeny of the arthropods. Curr Biol. 2001;11(10):759–63.
  369. Richter S. The Tetraconata concept: hexapod-crustacean relationships and the phylogeny of Crustacea. Org Divers Evol. 2002;2(3):217–37.
  370. Fanenbruck M, Harzsch S, Wägele JW. The brain of the Remipedia (Crustacea) and an alternative hypothesis on their phylogenetic relationships. Proc Natl Acad Sci U S A. 2004;101(11):3868–73.
  371. Oakley TH, Wolfe JM, Lindgren AR, Zaharoff AK. Phylotranscriptomics to Bring the Understudied into the Fold: Monophyletic Ostracoda, Fossil Placement, and Pancrustacean Phylogeny. Mol Biol Evol. 2012;30(1):215–33.
  372. Beutel RG, Friedrich F, Yang X-K, Ge S-Q. Insect Morphology and Phylogeny: A Textbook for Students of Entomology [Internet]. Walter de Gruyter; 2013
  373. Blanke A, Koch M, Wipfler B, Wilde F, Misof B. Head morphology of Tricholepidion gertschi indicates monophyletic Zygentoma. Front Zool. 2014;11:16.
  374. Beier M. 5. Klassifikation. In: Handbuch der Zoologie IV Band: Arthropoda – 2 Hälfte: Insecta. Berlin: Walter De Gruyter; 1969. p. 1–17.
  375. Hadrys H, Simon S, Kaune B, Schmitt O, Schöner A, Jakob W, et al. Isolation of Hox Cluster Genes from Insects Reveals an Accelerated Sequence Evolution Rate. PLoS ONE. 2012;7(6):e34682.
  376. Letsch HO, Meusemann K, Wipfler B, Schütte K, Beutel R, Misof B. Insect phylogenomics: results, problems and the impact of matrix composition. Proc R Soc Lond B Biol Sci. 2012;279(1741):3282–90.
  377. Savard J, Tautz D, Lercher MJ. Genome-wide acceleration of protein evolution in flies (Diptera). BMC Evol Biol. 2006;6:7.
  378. Staniczek AH. The mandible of silverfish (Insecta: Zygentoma) and Mayflies (Ephemeroptera): its morphology and phylogenetic significance. Zool Anz. 2000;239(2):147–78.
  379. Bitsch J, Bitsch C, Bourgoin T, D’Haese C. The phylogenetic position of early hexapod lineages: morphological data contradict molecular data. Syst Entomol. 2004;29(4):433–40.
  380. Kristensen NP. The groundplan and basal diversification of the hexapods. Fortey, R. A., Thomas, R. H. (eds.). In: Arthropod Relationships. London: Chapman & Hall; 1997. p. 281–93.
  381. Bitsch C, Bitsch J. The phylogenetic interrelationships of the higher taxa of apterygote hexapods. Zool Scr. 2000;29:131–56.
  382. Shao H, Zhang Y, Xie R, Yin W. Mitochondria1 cytochromeb sequences variation of Protura and molecular systematics of Apterygota. Chin Sci Bull. 1999;44(22):2031–6.
  383. Giribet G, Edgecombe GD, Wheeler WC. Arthropod phylogeny based on eight molecular loci and morphology. Nature. 2001;413:157–61.
  384. Mallatt JM, Giribet G. Further use of nearly complete 28S and 18S rRNA genes to classify Ecdysozoa: 37 more arthropods and a kinorhynch. Mol Phylogenet Evol. 2006;40:772–94.
  385. Dell’Ampio E, Szucsich NU, Carapelli A, Frati F, Steiner G, Steinacher A, et al. 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. 2009;38:155–70.
  386. Hovmoller R, Pape T, Kallersjo M. The Palaeoptera Problem: Basal Pterygote Phylogeny Inferred from 18S and 28S rDNA Sequences. Cladistics. 2002;18(3):313–23.
  387. Pisani D, Poling LL, Lyons-Weiler M, Hedges B. The colonization of land by animals: molecular phylogeny and divergence times among arthropods. BMC Biol. 2004;2:1.
  388. Regier JC, Shultz JW. Elongation Factor-2: A Useful Gene for Arthropod Phylogenetics. Mol Phylogenet Evol. 2001;20(1):136–48.
  389. Rota-Stabelli O, Campbell L, Brinkmann H, Edgecombe GD, Longhorn SJ, Peterson KJ, et al. A congruent solution to arthropod phylogeny: phylogenomics, microRNAs and morphology support monophyletic Mandibulata. Proc Biol Sci. 2011;278(1703):298–306.
  390. Seeger W. Spezialmerkmale an Eihüllen und Embryonen von Psocoptera im Vergleich zu anderen Paraneoptera (Insecta): Psocoptera als monophyletische Gruppe. Stuttg Beitr Zur Naturkunde A. 1979;329:1–57.
  391. Blanke A, Machida R, Szucsich NU, Wilde F, Misof B. Mandibles with two joints evolved much earlier in the history of insects: dicondyly is a synapomorphy of bristletails, silverfish and winged insects. Syst Entomol. 2015;40(2):357–64.

 

Figure Legends

Fig. 1. Phylogeny modified from Börner 1904. Taxa are named by modern convention.

 

Fig. 2. Hennig’s 1969 phylogeny(11), combined and modified from the original figures. Numerals indicate fossils as Hennig listed in his figures: 1. Rhyniella; 2. Eopterum (no longer considered an insect); 3. Rhyniognatha; 4. Monura; 5. Triassomachilis; 6. Triplosoba pulchella ; 7. Permoplecoptera; 8. Alleged subgroups of Ephemeroptera; 9. Erasipteron; 10. Protodonata (Meganisoptera); 11. Protanisoptera; 12. Protozygoptera; 13. Stemgroup of Anisozygoptera+Anisoptera; 14. Sheimia sojanensis; 15. Protoelytroptera; 16. Mesoforficula and others; 17. Puknoblattina; 18. Paleozoic “Problattoidea” and Blattodea; 19. Oedischia; 20. Glosselytodea; 21. Sthenaropodidae; 22. Oedishiidae, Elcanidae; 23. Tettavus; 24. Triassolocusta; 25. Tcholmanvissia ; 26. “Paraplecoptera” sensu Sharov (now Eoblattida Handlirsch 1906); 27. Protoperlaria (now Prothorthoptera Handlirsch 1906); 28. Perlopsis and other definitive Plecoptera; 29. Permopsocodea; 30. Procicadellopsis; 31. Archipsyllidae; 32. Permothrips longipennis; 33. Permaphidopsis; 34. Mesococcus asiaticus; 35. Archescytinidae; 36. Cicadopsyllidae; 37. Permaleurodes rotundatus; 38. Auchenorrhycha; 39. Paraknightia; 40. Boreocixius; 41. Permosialis; 42. Palaeohemerobiidae and Permithonidae, sensu Carpenter; 43. Tshekardocoleus and other branches; 44. Archezyela; 45. Mecoptera from Australia; 46., 47. Paratrichoptera; 48. Microptysma; 49. Microptysmodes

 

Fig. 3. Modified from Wille 1960 (31). Taxa are named by modern convention.

 

Fig. 4. Main image: Electron micrograph of a male Stylops ovinae (Strepsiptera). All insects have a 3 segmented thorax, each with a pair of legs. Wings, when present, are found on the second and third segments. Strepsipterans have reduced forewings modified as sense organs (arrows) attached to the small middle thoracic segment. Their anterior thoracic segment is greatly reduced. Flies have similarly reduced hindwings, attached to the 3rd thoracic segment. The 3rd thoracic segments of Strepsipterans and beetles are highly expanded, containing the functional wings and associated flight muscles. Image copyright Hans Pohl, used with permission. Insert: Wikipedia creative commons.

 

Fig. 5. Current consensus, modified from Misof et al.(50). Previous studies mentioned in this review are numbered and color coded on the left, with nodes they supported on the right. Red; morphology, without formalized data matrices. Orange: Morphology, with computer analysis. Blue: Sanger sequenced data in which rRNA played a predominant role. Black: Sanger sequenced multiple nuclear protein-coding genes. Green: Large genomic or transcriptomic data. Panel B. The sizes of datasets, plotted through time. The Y axis is on a log scale. Colors as in Panel A. Data size is calculated by multiplying the number of taxa by the number of characters. For works where amino acids were used as characters, we multiplied the number of characters by 3, so that these datasets were comparable to nucleotide datasets. Transcriptome work often has many missing data, so that character numbers were multiplied by the proportion of data present.

 

Table 1. Names of higher taxa used in the text, and the groups they define. Other groups are indicated on Fig. 5. Taxa in bold are supported in the current consensus. Taxa in italics have been strongly rejected. Those without indication are (in the opinion of these authors) targets for additional attention Citations are numbered as in the literature cited section, with the name of the first author and the last two digits of the publication year, for quick reference: (9)=Börner 04; (10)= Hennig 53; (11)=Hennig 69; (14)=Crampton 38; (31)=Wille 60; (34)=Mickoleit 73; (39)=Martynov 25; (47)=Wiegmann 09; (48)=Beutel 11; (49)=Peters 14; (50)=Misof 14; (62)=Kristensen 75; (63)=Kristensen 81; (69)=Boudreaux 79; (72)=Beutel 01; (74)=Wheeler 01; (76)=Wheeler 93; (77)=Whiting 97; (100)=Wipfler 11; (101)=Blanke 12; (107)=Kjer 04; (108)=Field 88; (110)=Turbeville 91; (114)=Wheeler 89; (116)=Pashley 93; (117)=Reumont 09; (140)=Edgecombe 00; (141)=Kjer 06; (152)=Beutel 06; (154)=McKenna 10; (159)=Chalwatzis 96; (170)=Ishiwata 11; (171)=Huelsenbeck 98; (175)=Niehuis 12; (177)=Boussau 14; (204)=Terry 05; (232)=Friedrich 95; (239)=Misof 07; (242)=Luan 05; (243)=Giribet 04; (244)=Gao 08; (245)=Mallatt 09; (248)=Carapelli 00; (249)=Krauss 04; (251)=Bonneton 06; (256)=Boore 95; (281)=Wan 12; (276)=Chen 14; (308)=Cryan 12; (319)=Aspöck 02; (434)=Sasaki 13; (347)=Regier 10; (349)=Meusemann 10; (351)=Simon 12; (355)=Letsch 13; (367)=Cook 01; (373)=Blanke 14; (374)=Beier 69; (375)=Hadrys 12; (376)=Letsch 12; (377)=Savard 06; (378)=Staniczek 00; (379)=Bitsch 04; (380)=Kristensen 97; (381)=Bitsch 00; (382)=Shao 99; (383)=Giribet 01; (384)=Mallatt 06; (385)=Dell’Ampio 09; (386)=Hovmoller 02; (387)=Pisani 04; (388)=Regier 01; (389)= Rota-Stabelli (390)= Seeger 79 (391)= Blanke 15

Leave a Reply

Contact Us At

Subscribe

Join our mailing list today and benefit from our free ebooks, daily deals, and discount