Michael F. Whiting

Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified

Nuclear mitochondrial pseudogenes (numts) are nonfunctional copies of mtDNA in the nucleus that have been found in major clades of eukaryotic organisms. They can be easily coamplified with orthologous mtDNA by using conserved universal primers; however, this is especially problematic for DNA barcoding, which attempts to characterize all living organisms by using a short fragment of the mitochondrial cytochrome c oxidase I (COI) gene. Here, we study the effect of numts on DNA barcoding based on phylogenetic and barcoding analyses of numt and mtDNA sequences in two divergent lineages of arthropods: grasshoppers and crayfish. Single individuals from both organisms have numts of the COI gene, many of which are highly divergent from orthologous mtDNA sequences, and DNA barcoding analysis incorrectly overestimates the number of unique species based on the standard metric of 3% sequence divergence. Removal of numts based on a careful examination of sequence characteristics, including indels, in-frame stop codons, and nucleotide composition, drastically reduces the incorrect inferences of the number of unique species, but even such rigorous quality control measures fail to identify certain numts. We also show that the distribution of numts is lineage-specific and the presence of numts cannot be known a priori. Whereas DNA barcoding strives for rapid and inexpensive generation of molecular species tags, we demonstrate that the presence of COI numts makes this goal difficult to achieve when numts are prevalent and can introduce serious ambiguity into DNA barcoding.

Mecoptera is paraphyletic: multiple genes and phylogeny of Mecoptera and Siphonaptera

Phylogenetic relationships among members of the Mecoptera and Siphonaptera were inferred from DNA sequence data. Four loci (18S and 28S ribosomal DNA, cytochrome oxidase II and elongation factor‐1α) were sequenced for 69 taxa selected to represent major flea and mecopteran lineages. Phylogenetic analyses of these data support a paraphyletic Mecoptera with two major lineages: Nannochoristidae + (Siphonaptera + Boreidae) and Meropidae + ((Choristidae + Apteropanorpidae) (Panorpidae + (Panorpidae + Bittacidae))). The flea family Ctenophthalmidae is paraphyletic, and the Ceratophylloidea is monophyletic. Morphological evidence is discussed which is congruent with the placement of Siphonaptera as sister group to Boreidae.

Loss and recovery of wings in stick insects

The evolution of wings was the central adaptation allowing insects to escape predators, exploit scattered resources, and disperse into new niches, resulting in radiations into vast numbers of species. Despite the presumed evolutionary advantages associated with full-sized wings (macroptery), nearly all pterygote (winged) orders have many partially winged (brachypterous) or wingless (apterous) lineages, and some entire orders are secondarily wingless (for example, fleas, lice, grylloblattids and mantophasmatids), with about 5% of extant pterygote species being flightless. Thousands of independent transitions from a winged form to winglessness have occurred during the course of insect evolution; however, an evolutionary reversal from a flightless to a volant form has never been demonstrated clearly for any pterygote lineage. Such a reversal is considered highly unlikely because complex interactions between nerves, muscles, sclerites and wing foils are required to accommodate flight. Here we show that stick insects (order Phasmatodea) diversified as wingless insects and that wings were derived secondarily, perhaps on many occasions. These results suggest that wing developmental pathways are conserved in wingless phasmids, and that ‘re-evolution’ of wings has had an unrecognized role in insect diversification.

Mantophasmatodea and phylogeny of the lower neopterous insects

Polyneoptera is a name sometimes applied to an assemblage of 11 insect orders comprising the lower neopterous or “orthopteroid” insects. These orders include familiar insects such as Orthoptera (grasshoppers), Blattodea (roaches), Isoptera (termites) (Mantodea) praying mantises, Dermaptera (earwigs), Phasmatodea (stick insects), Plecoptera (stoneflies), as well as the more obscure, Embiidina (web‐spinners), Zoraptera (angel insects) and Grylloblattodea (ice‐crawlers). Many of these insect orders exhibit a high degree of morphological specialization, a condition that has led to multiple phylogenetic hypotheses and little consensus among investigators. We present a phylogenetic analysis of the polyneopteran orders representing a broad range of their phylogenetic diversity and including the recently described Mantophasmatodea. These analyses are based on complete 18S rDNA, 28S rDNA, Histone 3 DNA sequences, and a previously published morphology matrix coded at the ordinal level. Extensive analyses utilizing different alignment methodologies and parameter values across a majority of possible ranges were employed to test for sensitivity of the results to ribosomal alignment and to explore patterns across the theoretical alignment landscape. Multiple methodologies support the paraphyly of Polyneoptera, the monophyly of Dictyoptera, Orthopteroidea (sensu Kukalova‐Peck; i.e. Orthoptera + Phasmatodea + Embiidina), and a group composed of Plecoptera + Dermaptera + Zoraptera. Sister taxon relationships between Embiidina + Phasmatodea in a group called “Eukinolabia”, and Dermaptera + Zoraptera (“Haplocercata”) are also supported by multiple analyses. This analysis also supports a sister taxon relationship between the newly described Mantophasmatodea, which are endemic to arid portions of southern Africa, and Grylloblattodea, a small order of cryophilic insects confined to the north‐western Americas and north‐eastern Asia, in a group termed “Xenonomia”. This placement, coupled with the morphological disparity of the two groups, validates the ordinal status of Mantophasmatodea.

Long-Branch Abstractions

Recent attention has been focused on the sensitivities of various tree reconstructing algorithms to sequence rate heterogeneity (long‐branch attraction). Phylogenetic conclusions from two recent empirical studies have been indicted as artifacts attributable to long‐branch attraction. Siddall et al.(1995) concluded that Myxozoa are cnidarians and sister group toPolypodium based on 18S rDNA and morphology. Haneltet al.(1996) argued that this result is due to long‐branch attraction. Whiting et al.(1997) concluded that the Strepsiptera are sister group to Diptera based on parsimony analysis of 18S rDNA, 28S rDNA, and morphology. Huelsenbeck (1997) argued that this result also is attributable to long‐branch attraction. We demonstrate that the analyses and arguments dismissing these results as the effects of long‐branch attraction are fundamentally flawed. The criteria employed by these authors were applied arbitrarily by them to the groups that they did not want, and yet using those same criteria, there is more reason to exclude other taxa besides Polypodium and there is more reason to disbelieve monophyly of Diptera than monophyly of Strepsiptera with Diptera. Moreover, it is asserted, long‐branch attraction cannot explain the presence of nematocysts in Myxozoa and halteres in Strepsiptera. For these reasons, and in light of the demonstration that long branches cannot attract each other in their mutual absence, we conclude that the monophyly of Myxozoa + Polypodium and Strepsiptera + Diptera is not due to long‐branch attraction. We suggest that maximum likelihood methods are extremely sensitive to taxon and character sampling and that these data sets are demonstrative of the long‐branch repulsion problem.

A mitochondrial genome phylogeny of Diptera: whole genome sequence data accurately resolve relationships over broad timescales with high precision

Abstract Mitochondrial genomes provide a promising new tool for understanding deep‐level insect phylogenetics, but have yet to be evaluated for their ability to resolve intraordinal relationships. We tested the utility of mitochondrial genome data for the resolution of relationships within Diptera, the insect order for which the most data are available. We sequenced an additional three genomes, from a syrphid, nemestrinid and tabanid, representing three additional dipteran clades, ‘aschiza’, non‐heteroneuran muscomorpha and ‘basal brachyceran’, respectively. We assessed the influence of optimality criteria, gene inclusion/exclusion, data recoding and partitioning strategies on topology and nodal support within Diptera. Our consensus phylogeny of Diptera was largely consistent with previous phylogenetic hypotheses of the order, except that we did not recover a monophyletic Muscomorpha (Nesmestrinidae grouped with Tabanidae) or Acalyptratae (Drosophilidae grouped with Calliphoridae). The results were very robust to optimality criteria, as parsimony, likelihood and Bayesian approaches yielded very similar topologies, although nodal support varied. The addition of ribosomal and transfer RNA genes to the protein coding genes traditionally used in mitochondrial genome phylogenies improved the resolution and support, contrary to previous suggestions that these genes would evolve too quickly or prove too difficult to align to provide phylogenetic signal at deep nodes. Strategies to recode data, aimed at reducing homoplasy, resulted in a decrease in tree resolution and branch support. Bayesian analyses were highly sensitive to partitioning strategy: biologically realistic partitions into codon groups produced the best results. The implications of this study for dipteran systematics and the effective approaches to using mitochondrial genome data are discussed. Mitochondrial genomes resolve intraordinal relationships within Diptera accurately over very wide time ranges (1–200 million years ago) and genetic distances, suggesting that this may be an excellent data source for deep‐level studies within other, less studied, insect orders.

A Comparative Analysis of Mitochondrial Genomes in Coleoptera (Arthropoda: Insecta) and Genome Descriptions of Six New Beetles

Coleoptera is the most diverse group of insects with over 360,000 described species divided into four suborders: Adephaga, Archostemata, Myxophaga, and Polyphaga. In this study, we present six new complete mitochondrial genome (mtgenome) descriptions, including a representative of each suborder, and analyze the evolution of mtgenomes from a comparative framework using all available coleopteran mtgenomes. We propose a modification of atypical cox1 start codons based on sequence alignment to better reflect the conservation observed across species as well as findings of TTG start codons in other genes. We also analyze tRNA-Ser(AGN) anticodons, usually GCU in arthropods, and report a conserved UCU anticodon as a possible synapomorphy across Polyphaga. We further analyze the secondary structure of tRNA-Ser(AGN) and present a consensus structure and an updated covariance model that allows tRNAscan-SE (via the COVE software package) to locate and fold these atypical tRNAs with much greater consistency. We also report secondary structure predictions for both rRNA genes based on conserved stems. All six species of beetle have the same gene order as the ancestral insect. We report noncoding DNA regions, including a small gap region of about 20 bp between tRNA-Ser(UCN) and nad1 that is present in all six genomes, and present results of a base composition analysis.

Mitochondrial genome data alone are not enough to unambiguously resolve the relationships of Entognatha, Insecta and Crustacea sensu lato (Arthropoda)

An analysis of the relationships of the major arthropod groups was undertaken using mitochondrial genome data to examine the hypotheses that Hexapoda is polyphyletic and that Collembola is more closely related to branchiopod crustaceans than insects. We sought to examine the sensitivity of this relationship to outgroup choice, data treatment, gene choice and optimality criteria used in the phylogenetic analysis of mitochondrial genome data. Additionally we sequenced the mitochondrial genome of an archaeognathan, Nesomachilis australica, to improve taxon selection in the apterygote insects, a group poorly represented in previous mitochondrial phylogenies. The sister group of the Collembola was rarely resolved in our analyses with a significant level of support. The use of different outgroups (myriapods, nematodes, or annelids + mollusks) resulted in many different placements of Collembola. The way in which the dataset was coded for analysis (DNA, DNA with the exclusion of third codon position and as amino acids) also had marked affects on tree topology. We found that nodal support was spread evenly throughout the 13 mitochondrial genes and the exclusion of genes resulted in significantly less resolution in the inferred trees. Optimality criteria had a much lesser effect on topology than the preceding factors; parsimony and Bayesian trees for a given data set and treatment were quite similar. We therefore conclude that the relationships of the extant arthropod groups as inferred by mitochondrial genomes are highly vulnerable to outgroup choice, data treatment and gene choice, and no consistent alternative hypothesis of Collembolas relationships is supported. Pending the resolution of these identified problems with the application of mitogenomic data to basal arthropod relationships, it is difficult to justify the rejection of hexapod monophyly, which is well supported on morphological grounds.

The problem with “the Paleoptera Problem:” sense and sensitivity

While the monophyly of winged insects (Pterygota) is well supported, phylogenetic relationships among the most basal extant pterygote lineages are problematic. Ephemeroptera (mayflies) and Odonata (dragonflies) represent the two most basal extant lineages of winged insects, and determining their relationship with regard to Neoptera (remaining winged insects) is a critical step toward understanding insect diversification. A recent molecular analysis concluded that Paleoptera (Odonata Ephemeroptera) is monophyletic. However, we demonstrate that this result is supported only under a narrow range of alignment parameters. We have further tested the monophyly of Paleoptera using additional sequence data from 18SrDNA, 28S rDNA, and Histone 3 for a broader selection of taxa and a wider range of analytical methodologies. Our results suggest that the current suite of molecular data ambiguously resolve the three basal winged insect lineages and do not provide independent confirmation of Odonata + Neoptera as supported via morphological data.

A preliminary mitochondrial genome phylogeny of Orthoptera (Insecta) and approaches to maximizing phylogenetic signal found within mitochondrial genome data.

The phylogenetic utility of mitochondrial genomes (mtgenomes) is examined using the framework of a preliminary phylogeny of Orthoptera. This study presents five newly sequenced genomes from four orthopteran families. While all ensiferan and polyneopteran taxa retain the ancestral gene order, all caeliferan lineages including the newly sequenced caeliferan species contain a tRNA rearrangement from the insect ground plan tRNA(Lys)(K)-tRNA(Asp)(D) swapping to tRNA(Asp) (D)-tRNA(Lys) (K) confirming that this rearrangement is a possible molecular synapomorphy for this suborder. The phylogenetic signal in mtgenomes is rigorously examined under the analytical regimens of parsimony, maximum likelihood and Bayesian inference, along with how gene inclusion/exclusion, data recoding, gap coding, and different partitioning schemes influence the phylogenetic reconstruction. When all available data are analyzed simultaneously, the monophyly of Orthoptera and its two suborders, Caelifera and Ensifera, are consistently recovered in the context of our taxon sampling, regardless of the optimality criteria. When protein-coding genes are analyzed as a single partition, nearly identical topology to the combined analyses is recovered, suggesting that much of the signals of the mtgenome come from the protein-coding genes. Transfer and ribosomal RNAs perform poorly when analyzed individually, but contribute signal when analyzed in combination with the protein-coding genes. Inclusion of third codon position of the protein-coding genes does not negatively affect the phylogenetic reconstruction when all genes are analyzed together, whereas recoding of the protein-coding genes into amino acid sequences introduces artificial resolution. Over-partitioning in a Bayesian framework appears to have a negative effect in achieving convergence. Our findings suggest that the best phylogenetic inferences are made when all available nucleotide data from the mtgenome are analyzed simultaneously, and that the mtgenome data can resolve over a wide time scale from the Permian (approximately 260 MYA) to the Tertiary (approximately 50 MYA).