Two rounds of whole genome duplication in the ancestral vertebrate

The hypothesis that the relatively large and complex vertebrate genome was created by two ancient, whole genome duplications has been hotly debated, but remains unresolved. We reconstructed the evolutionary relationships of all gene families from the complete gene sets of a tunicate, fish, mouse, and human, and then determined when each gene duplicated relative to the evolutionary tree of the organisms. We confirmed the results of earlier studies that there remains little signal of these events in numbers of duplicated genes, gene tree topology, or the number of genes per multigene family. However, when we plotted the genomic map positions of only the subset of paralogous genes that were duplicated prior to the fish–tetrapod split, their global physical organization provides unmistakable evidence of two distinct genome duplication events early in vertebrate evolution indicated by clear patterns of four-way paralogous regions covering a large part of the human genome. Our results highlight the potential for these large-scale genomic events to have driven the evolutionary success of the vertebrate lineage.     

Calculation of Evolutionary Correlation between Individual Genes and Full-Length Genome: A Method Useful for Choosing Phylogenetic Markers for Molecular Epidemiology

Individual genes or regions are still commonly used to estimate the phylogenetic relationships among viral isolates. The genomic regions that can faithfully provide assessments consistent with those predicted with full-length genome sequences would be preferable to serve as good candidates of the phylogenetic markers for molecular epidemiological studies of many viruses. Here we employed a statistical method to evaluate the evolutionary relationships between individual viral genes and full-length genomes without tree construction as a way to determine which gene can match the genome well in phylogenetic analyses. This method was performed by calculation of linear correlations between the genetic distance matrices of aligned individual gene sequences and aligned genome sequences. We applied this method to the phylogenetic analyses of porcine circovirus 2 (PCV2), measles virus (MV), hepatitis E virus (HEV) and Japanese encephalitis virus (JEV). Phylogenetic trees were constructed for comparisons and the possible factors affecting the method accuracy were also discussed in the calculations. The results revealed that this method could produce results consistent with those of previous studies about the proper consensus sequences that could be successfully used as phylogenetic markers. And our results also suggested that these evolutionary correlations could provide useful information for identifying genes that could be used effectively to infer the genetic relationships.     

From gene trees to organismal phylogeny in prokaryotes: the case of the gamma-Proteobacteria

The rapid increase in published genomic sequences for bacteria presents the first opportunity to reconstruct evolutionary events on the scale of entire genomes. However, extensive lateral gene transfer (LGT) may thwart this goal by preventing the establishment of organismal relationships based on individual gene phylogenies. The group for which cases of LGT are most frequently documented and for which the greatest density of complete genome sequences is available is the gamma-Proteobacteria, an ecologically diverse and ancient group including free-living species as well as pathogens and intracellular symbionts of plants and animals. We propose an approach to multigene phylogeny using complete genomes and apply it to the case of the gamma-Proteobacteria. We first applied stringent criteria to identify a set of likely gene orthologs and then tested the compatibilities of the resulting protein alignments with several phylogenetic hypotheses. Our results demonstrate phylogenetic concordance among virtually all (203 of 205) of the selected gene families, with each of the exceptions consistent with a single LGT event. The concatenated sequences of the concordant families yield a fully resolved phylogeny. This topology also received strong support in analyses aimed at excluding effects of heterogeneity in nucleotide base composition across lineages. Our analysis indicates that single-copy orthologous genes are resistant to horizontal transfer, even in ancient bacterial groups subject to high rates of LGT. This gene set can be identified and used to yield robust hypotheses for organismal phylogenies, thus establishing a foundation for reconstructing the evolutionary transitions, such as gene transfer, that underlie diversity in genome content and organization.     

Mitochondrial gene rearrangements: new paradigm in the evolutionary biology and systematics

Mitochondrial (mt) genomic study may reveal significant insight into the molecular evolution and several other aspects of genome evolution such as gene rearrangements evolution, gene regulation, and replication mechanisms. Other questions such as patterns of gene expression mechanism evolution, genomic variation and its correlation with physiology, and other molecular and biochemical mechanisms can be addressed by the mt genomics. Rare genomic changes have attracted evolutionary biology community for providing homoplasy free evidence of phylogenetic relationships. Gene rearrangements are considered to be rare evolutionary events and are being used to reconstruct the phylogeny of diverse group of organisms. Mt gene rearrangements have been established as a hotspot for the phylogenetic and evolutionary analysis of closely as well as distantly related organisms.     

Determination of whole prokaryotic phylogeny by the development of a random extraction method

The construction of accurate prokaryotic phylogeny is important not only in the field of evolutionary biology, but also in microbiology and pathology. However, in constructing a phylogenetic tree to trace prokaryotic evolution, the phylogenetic relationship is often changed by the choice of species. For the estimation of the accurate lineage of prokaryotes, a new method, named the "random extraction method", was developed. In this method, 16S rRNA sequence data were randomly extracted 1000 times from each closely-related taxa such as seven phyla of Eubacteria and one domain of Archaea and phylogenetic trees were constructed by the data to clarify the relationship of those groups. Next, the tree topology was counted and the most supported tree topology was found as the most plausible phylogenetic tree. To evaluate the reliability of each node, we developed the "Branching rate" (BR) and calculated for every tree. And also, computational simulation analysis was carried out to confirm these methods. On the assumption that the root of life is between Archaea and Eubacteria, the obtained phylogenetic relationships of phyla are the following. At first, Archaea (Euryarchaeota, Crenarchaeota and Korarchaeota) diverged, and Thermotogales, Cyanobacteria and Chlamydiales diverged in this order, then Firmicutes (Actinobacteria and Bacillus/Clostridium group cluster) and Proteobacteria (alpha and beta/gamma cluster) diverged. In addition, it was shown by the BR that the position of the node of Firmicutes Actinobacteria and Firmicutes Bacillus/Clostridium was changeable for each extraction. Therefore, it was suggested that the differences among the phylogenetic trees of prokaryotes were caused by the influence of these phyla.     

Grass Evolution Inferred from Chromosomal Rearrangements and Geometrical and Statistical Features in RNA Structure

The grasses (Poaceae) represent a monophyletic lineage that arose about 70 million years ago. The lineage contains about 10,000 species that differ widely in morphology and physiology. Species show striking differences in genome size, a feature important in the context of conservation of gene content and order (synteny and colinearity) and in the extension of genomic information directly from one grass species to another using comparative approaches. Grass diversification has been a contentious issue, as the exact branching order of the various subfamilies has been difficult to establish with standard methods. This motivated an evolutionary study of deep phylogenetic relationships based on the structure of coding and non-coding RNA molecules and on chromosomal rearrangements. Phylogenetic relationships in the grass family were inferred directly from the structure of RNA using cladistic principles and considerations in statistical mechanics. Coded attributes describing topological and thermodynamic information embedded in RNA molecules were treated as linearly ordered multi-state characters and were polarized by fixing the direction of character transformation toward molecular order. Intrinsically rooted phylogenies derived from the structure of signal recognition particle (SRP) RNA, the mRNA encoded by the early nodulation gene enod40, the small subunit of ribosomal RNA (rRNA), and the internal transcribed spacer ITS1 of rRNA established an order for the diversification of major grass lineages, suggesting a sister relationship of the Pooideae and the PACCAD clade. This same conclusion was reached when large-scale chromosomal rearrangements derived from the comparative genetic mapping of cereal genomes were studied. Chromosomal complements aligned in the most parsimonious manner allowed identification and coding of characters depicting chromosomal translocations, insertions, and linkage block arrangements and the reconstruction of phylogenetic trees based on large-scale chromosomal structure. Congruent reconstruction of deep branching relationships using geometrical and statistical features of RNA structure and orthology and large scale chromosomal recombination events support assumptions of polarization in character argumentation, and fail to falsify the claim that extant grass chromosomes can be considered combinations of linkage blocks of an ancestor of the rice genome. Congruence also suggests that the universal tendency toward order in RNA and the search for the most parsimonious organization of be genome architecture appear to be mutually supported drivers of molecular evolution. The study clarifies the relationship of major clades in the grasses, shows that phylogenetic history can be reconstructed effectively from the combinatorial exchange of chromosomal linkage blocks, and reveals considerable phylogenetic signal embedded in the structure of signal polypeptide-coding mRNA molecules, describing an instance where mRNA structure is the subject of strong evolutionary constraint.Reviewing Editor: Dr. David Pollock     

From Gene Trees to Organismal Phylogeny in Prokaryotes:The Case of the γ-Proteobacteria

The rapid increase in published genomic sequences for bacteria presents the first opportunity to reconstruct evolutionary events on the scale of entire genomes. However, extensive lateral gene transfer (LGT) may thwart this goal by preventing the establishment of organismal relationships based on individual gene phylogenies. The group for which cases of LGT are most frequently documented and for which the greatest density of complete genome sequences is available is the γ-Proteobacteria, an ecologically diverse and ancient group including free-living species as well as pathogens and intracellular symbionts of plants and animals. We propose an approach to multigene phylogeny using complete genomes and apply it to the case of the γ-Proteobacteria. We first applied stringent criteria to identify a set of likely gene orthologs and then tested the compatibilities of the resulting protein alignments with several phylogenetic hypotheses. Our results demonstrate phylogenetic concordance among virtually all (203 of 205) of the selected gene families, with each of the exceptions consistent with a single LGT event. The concatenated sequences of the concordant families yield a fully resolved phylogeny. This topology also received strong support in analyses aimed at excluding effects of heterogeneity in nucleotide base composition across lineages. Our analysis indicates that single-copy orthologous genes are resistant to horizontal transfer, even in ancient bacterial groups subject to high rates of LGT. This gene set can be identified and used to yield robust hypotheses for organismal phylogenies, thus establishing a foundation for reconstructing the evolutionary transitions, such as gene transfer, that underlie diversity in genome content and organization.     

The distribution of genes in the Drosophila genome

In the present work we show that in the Drosophila genome (which covers a 37-51% GC range at a DNA size of approx.50kb) a linear correlation holds between GC (or GC(3)50kb) genomic sequences embedding them. This correlation allows us to position the two compositional distributions of (a) coding sequences, and (b) of long DNA segments relative to each other and to calculate gene concentration across the compositional range of the Drosophila genome. Using this approach, we show that gene concentration increases with increasing GC of the regions embedding the genes, reaching a 7-fold higher level in the GC-richest regions compared with the GC-poorest regions. The gene distribution of the Drosophila genome is, therefore, similar to (although less striking than) that of the human genome, whereas it is very different from those of the Arabidopsis genome, which has about the same size as the Drosophila genome.     

From Gene Trees to Organismal Phylogeny in Prokaryotes:The Case of the γ-Proteobacteria

The rapid increase in published genomic sequences for bacteria presents the first opportunity to reconstruct evolutionary events on the scale of entire genomes. However, extensive lateral gene transfer (LGT) may thwart this goal by preventing the establishment of organismal relationships based on individual gene phylogenies. The group for which cases of LGT are most frequently documented and for which the greatest density of complete genome sequences is available is the γ-Proteobacteria, an ecologically diverse and ancient group including free-living species as well as pathogens and intracellular symbionts of plants and animals. We propose an approach to multigene phylogeny using complete genomes and apply it to the case of the γ-Proteobacteria. We first applied stringent criteria to identify a set of likely gene orthologs and then tested the compatibilities of the resulting protein alignments with several phylogenetic hypotheses. Our results demonstrate phylogenetic concordance among virtually all (203 of 205) of the selected gene families, with each of the exceptions consistent with a single LGT event. The concatenated sequences of the concordant families yield a fully resolved phylogeny. This topology also received strong support in analyses aimed at excluding effects of heterogeneity in nucleotide base composition across lineages. Our analysis indicates that single-copy orthologous genes are resistant to horizontal transfer, even in ancient bacterial groups subject to high rates of LGT. This gene set can be identified and used to yield robust hypotheses for organismal phylogenies, thus establishing a foundation for reconstructing the evolutionary transitions, such as gene transfer, that underlie diversity in genome content and organization.     

梅花草属叶绿体基因组进化分析

叶绿体基因组序列变异和基因组成等特征可有效反映植物类群间的系统发育和进化关系。本研究利用Illumina高通量测序平台对梅花草属（Parnassia）及其近缘属5种植物的叶绿体基因组进行测序和组装,同时基于已发表的近缘种叶绿体基因组信息,对梅花草属叶绿体基因组结构特征、序列遗传变异和蛋白编码基因密码子偏好性比对分析。结果显示：梅花草属叶绿体基因组整体结构较为保守,均为四分体结构;梅花草多个基因出现假基因化,而本属其他物种叶绿体基因组成一致,均编码115个基因;与近缘属物种相比,本属所有物种均丢失rpl16基因的内含子;蛋白质编码基因的非同义/同义替代率比值较低,叶绿体基因可能经历纯化选择作用;密码子偏好性聚类结果与蛋白编码序列重建的系统发育关系结果一致。本研究表明选择压力可能在梅花草属叶绿体基因组蛋白编码基因进化过程中发挥作用,有助于进一步理解梅花草属植物的进化和适应机制。     

Trends in morphological evolution in homobasidiomycetes inferred using maximum likelihood: a comparison of binary and multistate approaches

The homobasidiomycetes is a diverse group of macrofungi that includes mushrooms, puffballs, coral fungi, and other forms. This study used maximum likelihood methods to determine if there are general trends (evolutionary tendencies) in the evolution of fruiting body forms in homobasidiomycetes, and to estimate the ancestral forms of the homobasidiomycetes and euagarics clade. Character evolution was modeled using a published 481-species phylogeny under two character-coding regimes: additive binary coding, using DISCRETE, and multistate (five-state) coding, using MULTISTATE. Inferences regarding trends in character evolution made under binary coding were often in conflict with those made under multistate coding, suggesting that the additive binary coding approach cannot serve as a surrogate for multistate methods. MULTISTATE was used to develop a"minimal"model of fruiting body evolution, in which the 20 parameters that specify rates of transformations among character states were grouped into the fewest possible rate categories. The minimal model required only four rate categories, one of which is approaching zero, and suggests the following conclusions regarding trends in evolution of homobasidiomycete fruiting bodies: (1) there is an active trend favoring the evolution of pileate-stipitate forms (those with a cap and stalk); (2) the hypothesis that the evolution of gasteroid forms (those with internal spore production, such as puffballs) is irreversible cannot be rejected; and (3) crustlike resupinate forms are not a particularly labile morphology. The latter finding contradicts the conclusions of a previous study that used binary character coding. Ancestral state reconstructions under binary coding suggest that the ancestor of the homobasidiomycetes was resupinate and the ancestor of the euagarics clade was pileate-stipitate, but ancestral state reconstructions under multistate coding did not resolve the ancestral form of either node. The results of this study illustrate the potential sensitivity of comparative analyses to character state definitions.     

Phylogenetic interpretation of ontogenetic change: sorting out the actual and artefactual in an empirical case study of centrarchid fishes

Hypothesized relationships between ontogenetic and phylogenetic change in morphological characters were empirically tested in centrarchid fishes by comparing observed patterns of character development with patterns of character evolution as inferred from a representative phylogenetic hypothesis. This phylogeny was based on 56–61 morphological characters that were polarized by outgroup comparison. Through these comparisons, evolutionary changes in character ontogeny were categorized in one of eight classes (terminal addition, terminal deletion, terminal substitution, non-terminal addition, non-terminal deletion, non-terminal substitution, ontogenetic reversal and substitution). The relative frequencies of each of these classes provided an empirical basis from which assumptions underlying hypothesized relationships between ontogeny and phylogeny were tested. In order to test hypothesized relationships between ontogeny and phylogeny that involve assumptions about the relative frequencies of terminal change (e.g. the use of ontogeny as a homology criterion), two additional phylogenies were generated in which terminal addition and terminal deletion were maximized and minimized for all characters. Character state change interpreted from these phylogenies thus represents the maxima and minima of the frequency range of terminal addition and terminal deletion for the 8.7 × 10³⁶ trees possible for centrarchids. It was found for these data that terminal change accounts for c. 75% of the character state change. This suggests either that early ontogeny is conserved in evolution or that interpretation and classification of evolutionary changes in ontogeny is biased in part by the way that characters are recognized, delimited and coded. It was found that ontogenetic interpretation is influenced by two levels of homology decision: an initial decision involving delimitation of the character (the ontogenetic sequence), and the subsequent recognition of homologous components of developmental sequences. Recognition of phylogenetic homology among individual components of developmental sequences is necessary for interpretation of evolutionary changes in ontogeny as either terminal or non-terminal. If development is the primary criterion applied in recognizing individual homologies among parts of ontogenetic sequences, the only possible interpretation of phylogenetic differences is that of terminal change. If homologies of the components cannot be ascertained, recognition of the homology of the developmental sequence as a whole will result in the interpretation of evolutionary differences as substitutions. Particularly when the objective of a study is to discover how ontogeny has evolved, criteria in addition to ontogeny must be used to recognize homology. Interpretation is also dependent upon delimitation within an ontogenetic sequence. This is in part a function of the way that an investigator ‘sees’ and codes characters. Binary and multistate characters influence interpretation differently and predictably. The use of ontogeny for determining phylogenetic polarity as previously proposed rests on the assumptions that ancestral ontogenies are conserved and that character evolution occurs predominantly through terminal addition. It was found for these data that terminal addition may comprise a maximum of 51.9% of the total character state change. It is concluded that the ontogenetic criterion is not a reliable indicator of phylogenetic polarity. Process and pattern data are collected simultaneously by those engaged in comparative morphological studies of development. The set of alternative explanatory processes is limited in the process of observing development. These form necessary starting points for the research of developmental biologists. Separating ‘empirical’ results from interpretational influences requires awareness of potential biases in the course of character selection, coding and interpretation. Consideration of the interpretational problems involved in identifying and classifying phylogenetic changes in ontogeny leads to a re-evaluation of the purpose, usefulness and information conveyed by the current classification system. It is recommended that alternative classification schemes be pursued.     

Exploring phylogenetic hypotheses via Gibbs sampling on evolutionary networks

Background

Phylogenetic networks are leaf-labeled graphs used to model and display complex evolutionary relationships that do not fit a single tree. There are two classes of phylogenetic networks: Data-display networks and evolutionary networks. While data-display networks are very commonly used to explore data, they are not amenable to incorporating probabilistic models of gene and genome evolution. Evolutionary networks, on the other hand, can accommodate such probabilistic models, but they are not commonly used for exploration.

Results

In this work, we show how to turn evolutionary networks into a tool for statistical exploration of phylogenetic hypotheses via a novel application of Gibbs sampling. We demonstrate the utility of our work on two recently available genomic data sets, one from a group of mosquitos and the other from a group of modern birds. We demonstrate that our method allows the use of evolutionary networks not only for explicit modeling of reticulate evolutionary histories, but also for exploring conflicting treelike hypotheses. We further demonstrate the performance of the method on simulated data sets, where the true evolutionary histories are known.

Conclusion

We introduce an approach to explore phylogenetic hypotheses over evolutionary phylogenetic networks using Gibbs sampling. The hypotheses could involve reticulate and non-reticulate evolutionary processes simultaneously as we illustrate on mosquito and modern bird genomic data sets.

     

2R or not 2R: Testing hypotheses of genome duplication in early vertebrates

The widely popular hypothesis that there were two rounds of genome duplication by polyploidization early in vertebrate history (the 2R hypothesis) has been difficult to test until recently. Among the lines of evidence adduced in support of this hypothesis are relative genome size, relative gene number, and the existence of genomic regions putatively duplicated during polyploidization. The availability of sequence for a substantial portion of the human genome makes possible the first rigorous tests of this hypothesis. Comparison of gene family size in the human genome and in invertebrate genomes shows no evidence of a 4:1 ratio between vertebrates and invertebrates. Furthermore, explicit phylogenetic tests for the topology expected from two rounds of polyploidization have revealed alternative topologies in a substantial majority of human gene families. Likewise, phylogenetic analyses have shown that putatively duplicated genomic regions often include genes duplicated at widely different times over the evolution of life. The 2R hypothesis thus can be decisively rejected. Rather, current evidence favors a model of genome evolution in which tandem duplication, whether of genomic segments or of individual genes, predominates.     

Amborella not a "basal angiosperm"? Not so fast

The sequence of the plastid genome of Amborella trichopoda, the putative sister to all other extant angiosperms, was recently reported (Molecular Biology and Evolution 20: 1499-1505). Goremykin et al. used sequence data for 61 plastid genes from Amborella and 12 other embryophytes in phylogenetic analyses and concluded that Amborella is not the sister to the remaining flowering plants; the monocots instead occupy this position. The authors attributed their results, which differ substantially from all recent phylogenetic analyses of angiosperms, to the increased character sampling (30?017 nucleotides in their aligned matrix) in their analysis relative to published studies that included fewer genes but more taxa. We hypothesized that the difference in topology is not due to limited character sampling in previous studies but to limited taxon sampling in the analysis by Goremykin et al. To test this, we conducted a series of phylogenetic analyses using a three-gene, 12 (or more)-taxon data set to evaluate the topological effects of (i) including three vs. 61 genes for (nearly) the same set of taxa, (ii) analyzing different codon positions, (iii) substituting representatives of other basal lineages for Amborella, (iv) replacing the grasses used to represent the monocots with other monocots, selected either for their phylogenetic position or randomly, and (v) adding other basal taxa-Nymphaea, Austrobaileya, magnoliids, and monocots-to the 12-taxon data set. Our results demonstrate that the "monocots basal" topology obtained by Goremykin et al. is not due to increased character sampling of the plastid genome; their topology was obtained using only two plastid genes or two plastid genes and one nuclear gene. This topology was also retained when either Nymphaea or Austrobaileya was substituted for Amborella, demonstrating that any of the three basal lineages will attach to Calycanthus for lack of any other close branch. Furthermore, the "monocots basal" topology is not robust to changes in sampling of monocots. Simply adding Oncidium, for example, places Amborella sister to the other angiosperms. Thus, limited taxon sampling, focusing on organisms with complete genome sequences, can lead to artifactual results.     

THE EFFECT OF ORDERED CHARACTERS ON PHYLOGENETIC RECONSTRUCTION

Abstract Morphological structures are likely to undergo more than a single change during the course of evolution. As a result, multistate characters are common in systematic studies and must be dealt with. Particularly interesting is the question of whether or not multistate characters should be treated as ordered (additive) or unordered (non-additive). In accepting a particular hypothesis of order, numerous others are necessarily rejected. We review some of the criteria often used to order character states and the underlying assumptions inherent in these criteria.
The effects that ordered multistate characters can have on phylogenetic reconstruction are examined using 27 data sets. It has been suggested that hypotheses of character state order are more informative then hypotheses of unorder and may restrict the number of equally parsimonious trees as well as increase tree resolution. Our results indicate that ordered characters can produce more, equal or less equally parsimonious trees and can increase, decrease or have no effect on tree resolution. The effect on tree resolution can be a simple gain in resolution or a dramatic change in sister-taxa relationships. In cases where several outgroups are included in the data matrix, hypotheses of order can change character polarities by altering outgroup topology. Ordered characters result in a different topology from unordered characters only when the hierarchy of the cladogram disagrees with the investigator's a priori hypothesis of order. If the best criterion for assessing character evolution is congruence with other characters, the practice of ordering multistate characters is inappropriate.     

Supervised Lowess normalization of comparative genome hybridization data – application to lactococcal strain comparisons

Background

Visualising the evolutionary history of a set of sequences is a challenge for molecular phylogenetics. One approach is to use undirected graphs, such as median networks, to visualise phylogenies where reticulate relationships such as recombination or homoplasy are displayed as cycles. Median networks contain binary representations of sequences as nodes, with edges connecting those sequences differing at one character; hypothetical ancestral nodes are invoked to generate a connected network which contains all most parsimonious trees. Quasi-median networks are a generalisation of median networks which are not restricted to binary data, although phylogenetic information contained within the multistate positions can be lost during the preprocessing of data. Where the history of a set of samples contain frequent homoplasies or recombination events quasi-median networks will have a complex topology. Graph reduction or pruning methods have been used to reduce network complexity but some of these methods are inapplicable to datasets in which recombination has occurred and others are procedurally complex and/or result in disconnected networks.

Results

We address the problems inherent in construction and reduction of quasi-median networks. We describe a novel method of generating quasi-median networks that uses all characters, both binary and multistate, without imposing an arbitrary ordering of the multistate partitions. We also describe a pruning mechanism which maintains at least one shortest path between observed sequences, displaying the underlying relations between all pairs of sequences while maintaining a connected graph.

Conclusion

Application of this approach to 5S rDNA sequence data from sea beet produced a pruned network within which genetic isolation between populations by distance was evident, demonstrating the value of this approach for exploration of evolutionary relationships.     

Insufficient resolving power of mitogenome data in deciphering deep phylogeny of Holometabola

The mitochondrial genome (mitogenome) is one of the most widely used markers for phylogenetic analysis. Compared with whole-genome data, mitogenome data are less expensive to obtain and easier to manipulate. However, compositional bias and accelerated evolutionary rate reduce the effectiveness of the mitogenome in determining insect phylogeny. This study shows that mitogenome data are not suitable to reconstruct deep holometabolan evolution, even with a most comprehensive data coding scheme and the more realistic CAT model. For the deep levels of divergence within Holometabola, protein-coding genes only retain weak phylogenetic signals, leading to peculiar interordinal relationships. Consensus relationships in the Holometabola phylogeny, such as the monophyly of Holometabola, the most basal position of Hymenoptera, and the sister group relationship between the Strepsiptera and Coleoptera were rarely resolved in our analyses. The relationships of the holometabolan groups as inferred by mitogenomes are highly vulnerable to gene types, data coding regimes, model choice, and optimality criteria, and no consistent alternative hypothesis of Holometabola's relationships is supported. Thus, we suggest that the slowly evolving nuclear genes or genome-scale approaches may be better options for resolving deep-level phylogeny of Holometabola.