“…To illuminate the phylogenetic relationships of Phasmatodea, we first performed maximum likelihood (ML) and Bayesian inference (BI) analyses based on data from 85 species, including the three newly determined sequences, sixty-five previously sequenced mitochondrial genomes, and 13 PCGs of seventeen species of Phasmatodea assembled from transcriptome data [ 14 , 33 , 34 , 36 , 55 , 56 , 57 , 58 , 59 , 60 , 61 , 62 , 63 , 64 , 65 , 66 , 67 , 68 , 69 , 70 , 71 , 72 , 73 , 74 , 75 , 76 , 77 ] ( Table 1 , Table 2 and Table S2 ). Long-branch attraction between Embioptera and Zoraptera appeared in both ML and BI phylogenetic trees and affected the stability of the topology.…”