Mitochondrial genomes of stick insects (Phasmatodea) and phylogenetic considerations

Insight into the Phylogenetic Relationships of Phasmatodea and Selection Pressure Analysis of Phraortes liaoningensis Chen & He, 1991 (Phasmatodea: Lonchodidae) Using Mitogenomes

Stick and leaf insects are a group among the Insecta that are famous for their extraordinary mimicry ability. Since the establishment of the Phasmatodea, their internal classification has been constantly revised. Mitochondrial genes as molecular markers have been widely used for species classification, but the phylogenetic relationships within the Phasmatodea remain to be thoroughly discussed. In the present study, five mitogenomes of Phasmatodea ranging from 15,746 bp to 16,747 bp in length were sequenced. Bayesian inference (BI) and maximum likelihood (ML) analyses were carried out based on a 13 PCGs data matrix (nt123) and a combined matrix of 13 PCGs and two rRNA genes (nt123_rRNA). The present study supports the conclusion that Phylliidae was the basal group of Neophasmatodea and confirms the monophyly of Lonchodinae and Necrosciinae, but it shows that Lonchodidae was polyphyletic. A sister group of Bacillidae and Pseudophasmatidae was also recovered. The phylogenetic tree based on the nt_123 dataset showed higher node support values. The construction of a divergent time tree in this study supported the conclusion that extant Phasmatodea originated in the Jurassic (170 Mya) and most lineages diverged after the Cretaceous-Paleogene extinction event. To explore whether the mitochondrial genes of Phraortes liaoningensis collected from high latitudes where low temperatures occur for eight months of the year are under selection pressure, this study used the branch-site model and the branch model to analyze the selection pressure on the 13 mitochondria protein-coding genes (PCGs). We found that both ND2 and ND4L of Ph. liaoningensis exhibited positive selection sites using the branch-site model. This study shows that a low-temperature environment causes mitochondrial genes to be selected to meet the energy requirements for survival.

A comprehensive sampling of mitogenomes shows the utility to infer phylogeny of termites (Blattodea: Termitoidae)

The mitogenome sequence data have been widely used in inferring the phylogeny of insects. In this study, we determined the complete mitogenome for Macrotermes sp. (Termitidae, Macrotermitinae) using next-generation sequencing. Macrotermes sp. possesses a typical insect mitogenome, displaying an identical gene order and gene content to other existing termite mitogenomes. We present the first prediction of the secondary structure of ribosomal RNA genes in termites. The rRNA secondary structures of Macrotermes sp. exhibit similarities to closely related insects and also feature distinctive characteristics in their helical structures. Together with 321 published mitogenomes of termites as ingroups and 8 cockroach mitogenomes as outgroups, we compiled the most comprehensive mitogenome sequence matrix for Termitoidae to date. Phylogenetic analyses were conducted using datasets employing different data coding strategies and various inference methods. Robust relationships were recovered at the family or subfamily level, demonstrating the utility of comprehensive mitogenome sampling in resolving termite phylogenies. The results supported the monophyly of Termitoidae, and consistent relationships within this group were observed across different analyses. Mastotermitidae was consistently recovered as the sister group to all other termite families. The families Hodotermitidae, Stolotermitidae, and Archotermopsidae formed the second diverging clade, followed by the Kalotermitidae. The Neoisoptera was consistently supported with strong node support, with Stylotermitidae being sister to the remaining families. Rhinotermitidae was found to be non-monophyletic, and Serritermitidae nested within the basal clades of Rhinotermitidae and was sister to Psammotermitinae. Overall, our phylogenetic results are largely consistent with earlier mitogenome studies.

Nine Mitochondrial Genomes of Phasmatodea with Two Novel Mitochondrial Gene Rearrangements and Phylogeny

The classification of stick and leaf insects (Order Phasmatodea) is flawed at various taxonomic ranks due to a lack of robust phylogenetic relationships and convergent morphological characteristics. In this study, we sequenced nine new mitogenomes that ranged from 15,011 bp to 17,761 bp in length. In the mitogenome of Carausis sp., we found a translocation of trnR and trnA, which can be explained by the tandem duplication/random loss (TDRL) model. In the Stheneboea repudiosa Brunner von Wattenwyl, 1907, a novel mitochondrial structure of 12S rRNA-CR1-trnI-CR2-trnQ-trnM was found for the first time in Phasmatodea. Due to the low homology of CR1 and CR2, we hypothesized that trnI was inverted through recombination and then translocated into the middle of the control region. Control region repeats were frequently detected in the newly sequenced mitogenomes. To explore phylogenetic relationships in Phasmatodea, mtPCGs from 56 Phasmatodean species (composed of 9 stick insects from this study, 31 GenBank data, and 16 data derived from transcriptome splicing) were used for Bayesian inference (BI), and maximum likelihood (ML) analyses. Both analyses supported the monophyly of Lonchodinae and Necrosciinae, but Lonchodidae was polyphyletic. Phasmatidae was monophyletic, and Clitumninae was paraphyletic. Phyllidae was located at the base of Neophasmatodea and formed a sister group with the remaining Neophasmatodea. Bacillidae and Pseudophasmatidae were recovered as a sister group. Heteroptergidae was monophyletic, and the Heteropteryginae sister to the clade (Obriminae + Dataminae) was supported by BI analysis and ML analysis.

Cytochrome P450 Genes Expressed in Phasmatodea Midguts

Cytochrome P450s (CYPs) are xenobiotic detoxification genes found in most eukaryotes, and linked in insects to the tolerance of plant secondary chemicals and insecticide resistance. The number and diversity of CYP clans, families, and subfamilies that an organism produces could correlate with its dietary breadth or specialization. This study examined the CYP diversity expressed in the midguts of six species of folivorous stick insects (Phasmatodea), to identify their CYP complement and see if any CYPs correlate with diet toxicity or specialization, and see what factors influenced their evolution in this insect order. CYP genes were mined from six published Phasmatodea transcriptomes and analyzed phylogenetically. The Phasmatodea CYP complement resembles that of other insects, though with relatively low numbers, and with significant expansions in the CYP clades 6J1, 6A13/14, 4C1, and 15A1. The CYP6 group is known to be the dominant CYP family in insects, but most insects have no more than one CYP15 gene, so the function of the multiple CYP15A1 genes in Phasmatodea is unknown, with neofunctionalization following gene duplication hypothesized. No correlation was found between CYPs and diet specialization or toxicity, with some CYP clades expanding within the Phasmatodea and others likely inherited from a common ancestor.

Two new stick insect species of Sosibia Stål (Phasmatodea: Lonchodidae: Necrosciinae) from China and the first report on mitochondrial genomes of this genus

We describe and illustrate two new species of Sosibia from China: Sosibia gibba sp. nov. and Sosibia ovata sp. nov. This report includes a key to Sosibia species from China and a description of the distribution area in China. The two mitochondrial genomes of these new Sosibia species were sequenced and annotated for the first time. The compositional biases, codon usage, nucleotide composition, and construct tRNA secondary structures of the two mitogenomes were analyzed. The phylogenetic relationships based on the mitogenomes using Bayesian inference and maximum likelihood methods supported the monophyly of Necrosciinae and divided it into two distinct clades: A: (Sipyloidea + [Sosibia + Calvisia]); and B: (Neohirasea + Micadina).

Organization of the mitochondrial genome of Ramulus irregulatiter dentatus (Phasmatidae: Phasmatidae)

Recently, the species of the insect order Phasmatodea, have attracted the interest of more and more enthusiasts. Here, we obtained the complete mitochondrial genome of Ramulus irregulatiter dentatus (R. irregulatiter dentatus), which belongs to the subfamily of Phasmatidae, detected by Illumina next-generation sequencing. The entire mitochondrial genome is 16,060 bp in length and contains a standard set of 13 protein-coding genes, 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and a putative A + T-rich region. The base composition and codon usage were typical of Phasmatodea species. The mitochondrial gene organization (37 genes) was consistent with that of other Phasmatidae. A phylogenetic tree was built from the sequence information of the 13 protein-coding genes by Bayesian analyses. The newly sequenced R. irregulatiter dentatus was most closely related to the family Phasmatidae. The complete mitochondrial genome of R. irregulatiter dentatus also provides valuable molecular information for future studies on Phasmatidae insect taxonomy and a framework to unveil more of their cryptic and unknown diversity, so that it can be used to control forest pests and protect crops.

The phylogenic position of aschiphasmatidae in euphasmatodea based on mitochondrial genomic evidence

The mitochondrial genome (mitogenome) has been regarded as significant source of data to better understand the phylogenetic relationships within the Euphasmatodea, but no mitogenome in Aschiphasmatoidea has been sequenced to date. In this study, two mitogenomes of Orthomeria smaragdinum and Nanhuaphasma hamicercum of Aschiphasmatidae were sequenced and annotated for the first time. The same mitochondrial gene rearrangement structure was present in the two mitogenomes sequenced, showing as the translocation of tRNA-Arg and tRNA-Asn, which conformed to the tandem duplication-random loss and could be used as a possible synapomorphy for Aschiphasmatidae. The phylogenetic results based on the maximum likelihood (ML) and bayesian inference (BI) methods both showed that Aschiphasmatidae and Neophasmatodea in Euphasmatodea are sister taxa. Although the monophyly of Oriophasmata, Occidophasmata, Diapheromeridae, Phasmatidae, Lonchodidae and Bacilloidea has not been solved, the monophyly of Neophasmatodea and Phyllioidea was well supported.

A mitochondrial genome phylogeny of voles and lemmings (Rodentia: Arvicolinae): Evolutionary and taxonomic implications

Arvicolinae is one of the most impressive placental radiations with over 150 extant and numerous extinct species that emerged since the Miocene in the Northern Hemisphere. The phylogeny of Arvicolinae has been studied intensively for several decades using morphological and genetic methods. Here, we sequenced 30 new mitochondrial genomes to better understand the evolutionary relationships among the major tribes and genera within the subfamily. The phylogenetic and molecular dating analyses based on 11,391 bp concatenated alignment of protein-coding mitochondrial genes confirmed the monophyly of the subfamily. While Bayesian analysis provided a high resolution across the entire tree, Maximum Likelihood tree reconstruction showed weak support for the ordering of divergence and interrelationships of tribal level taxa within the most ancient radiation. Both the interrelationships among tribes Lagurini, Ellobiusini and Arvicolini, comprising the largest radiation and the position of the genus Dinaromys within it also remained unresolved. For the first time complex relationships between genus level taxa within the species-rich tribe Arvicolini received full resolution. Particularly Lemmiscus was robustly placed as sister to the snow voles Chionomys in the tribe Arvicolini in contrast with a long-held belief of its affinity with Lagurini. Molecular dating of the origin of Arvicolinae and early divergences obtained from the mitogenome data were consistent with fossil records. The mtDNA estimates for putative ancestors of the most genera within Arvicolini appeared to be much older than it was previously proposed in paleontological studies.

The mitochondrial genome of Eurycantha calcarata Lucas, 1869 (Phasmatodea: Lonchodinae) and its phylogeny

The Lonchodinae (Phasmatodea: Phasmatidae) is rich in insect species with more than 330 species of 40 genera. The phylogenetic relationships within Lonchodinae have been under debate. We successfully sequenced the complete mitogenome of Eurycantha calcarata Lucas, 1869 (Phasmatodea: Lonchodinae) with a length of 16,280 bp, which had the same genes and gene arrangements as those of various published papers on stick insects. The whole mitogenome and control region of E. calcarata had a high AT content of 78.2 and 85.9%, respectively. All PCGs used ATN as the start codon, and most PCGs used TAA/TAG as the stop codons excluding COX2 (T), COX3 (TA), and ND5 (TA). To discuss the phylogeny of Lonchodinae, we reconstructed the phylogenetic relationships of 27 species of Phasmatodea including E. calcarata and two species of Embioptera used as outgroups. In BI and ML trees, the monophyly of Lonchodinae and Necrosciinae was well supported, whereas the monophyly of Clitumninae was not recovered. These results indicated that Lonchodinae was a sister clade to Phylliinae and E. calcarata was a sister clade to Phraortes genus.

Three Complete Mitochondrial Genomes of Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis (Insecta: Phasmatodea) and Their Phylogeny

Insects of the order Phasmatodea are mainly distributed in the tropics and subtropics and are best known for their remarkable camouflage as plants. In this study, we sequenced three complete mitochondrial genomes from three different families: Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis. The lengths of the three mitochondrial genomes were 15,896 bp, 16,869 bp, and 17,005 bp, respectively, and the gene composition and structure of the three stick insects were identical to those of the most recent common ancestor of insects. The phylogenetic relationships among stick insects have been chaotic for a long time. In order to discuss the intra- and inter-ordinal relationship of Phasmatodea, we used the 13 protein-coding genes (PCGs) of 85 species for maximum likelihood (ML) and Bayesian inference (BI) analyses. Results showed that the internal topological structure of Phasmatodea had a few differences in both ML and BI trees and long-branch attraction (LBA) appeared between Embioptera and Zoraptera, which led to a non-monophyletic Phasmatodea. Consequently, after removal of the Embioptera and Zoraptera species, we re-performed ML and BI analyses with the remaining 81 species, which showed identical topology except for the position of Tectarchus ovobessus (Phasmatodea). We recovered the monophyly of Phasmatodea and the sister-group relationship between Phasmatodea and Mantophasmatodea. Our analyses also recovered the monophyly of Heteropterygidae and the paraphyly of Diapheromeridae, Phasmatidae, Lonchodidae, Lonchodinae, and Clitumninae. In this study, Peruphasma schultei (Pseudophasmatidae), Phraortes sp. YW-2014 (Lonchodidae), and species of Diapheromeridae clustered into the clade of Phasmatidae. Within Heteropterygidae, O. guangxiensis was the sister clade to O. mouhotii belonging to Dataminae, and the relationship of (Heteropteryginae + (Dataminae + Obriminae)) was recovered.

A tree of leaves: Phylogeny and historical biogeography of the leaf insects (Phasmatodea: Phylliidae)

The insect order Phasmatodea is known for large slender insects masquerading as twigs or bark. In contrast to these so-called stick insects, the subordinated clade of leaf insects (Phylliidae) are dorso-ventrally flattened and therefore resemble leaves in a unique way. Here we show that the origin of extant leaf insects lies in the Australasian/Pacific region with subsequent dispersal westwards to mainland Asia and colonisation of most Southeast Asian landmasses. We further hypothesise that the clade originated in the Early Eocene after the emergence of angiosperm-dominated rainforests. The genus Phyllium to which most of the ~100 described species pertain is recovered as paraphyletic and its three non-nominate subgenera are recovered as distinct, monophyletic groups and are consequently elevated to genus rank. This first phylogeny covering all major phylliid groups provides the basis for future studies on their taxonomy and a framework to unveil more of their cryptic and underestimated diversity.

Cladistic analysis of Paraphasma (Phasmatodea: Pseudophasmatidae) highlights the importance of the phallic organ for phasmid systematics

The internal male genitalia have been poorly investigated in Phasmatodea, remaining virtually unexplored in phylogenetic studies. Here we describe and illustrate the main phallic elements in several Neotropical stick insects, with emphasis on Paraphasma (Pseudophasmatidae), and present a phylogenetic analysis of this genus. The analysis included ten terminals in the ingroup and 18 in the outgroup, and was based on 32 characters of the phallic organ and 48 of external morphology. In order to compare these datasets in terms of phylogenetic signal and level of homoplasy, the consistency and retention indices of the cladogram were calculated separately for each of them, and partial analyses were also conducted using each dataset alone. The phylogenetic reconstruction revealed Paraphasma as polyphyletic and led us to propose a new, monotypic genus, Ecuadoriphasma gen. nov., three new combinations (Ecuadoriphasma cognatum, Paraphasma trianguliferum and Tithonophasma cancellatum) and place Oestrophora as a synonym of Paraphasma. Additionally, Olcyphides hopii and Paraphasma dentatum are synonymized with Paraphasma laterale. Both external and phallic characters were determinant for the topology obtained, and the latter were less homoplastic in the phylogenetic tree. Our results highlight the usefulness of phallic morphology for inferring phylogenetic relationships in Phasmatodea, especially among closely related genera and species.