Content uploaded by Min Li
Author content
All content in this area was uploaded by Min Li on Sep 26, 2017
Content may be subject to copyright.
Entomotaxonomia (2016) 38(2): 81–91, DOI: 10.11680/entomotax.2016021, ISSN 2095–8609 81
Reanalysis of the phylogenetic relationships of the
Pentatomomorpha (Hemiptera: Heteroptera) based on
ribosomal, Hox and mitochondrial genes
Min LI1,2,3, Yanhui WANG1, Qiang XIE1, Xiaoxuan TIAN1,4, Teng LI1, Hufang
ZHANG3, Wenjun BU1①
1. College of Life Sciences, Nankai University, Tianjin 300071, China
2. Department of Biology, Taiyuan Normal University, Taiyuan, Shanxi 030031, China
3. College of Agriculture, Shanxi Agricultural University, Taigu, Shanxi 030801, China
4. Tianjin State Key Laboratory of Modern Chinese Medicine, Tianjin University of Traditional Chinese
Medicine, Tianjin 300193, China
Abstract: Pentatomomorpha is one of the most biodiverse infraorders among the true bugs (Hemiptera:
Heteroptera). Phylogenetic relationships among the superfamilies within this infraorder have been uncertain,
especially for the Eutrichophora. The previous studies were based on morphological characters, or just
mitochondrial or nuclear genes, or only partial 18S rDNA and COI. In this study, we used maximum likelihood
(ML) and Bayesian inference (BI) based on massive molecular datasets (18S rDNA, 28S rDNA, Hox and
mitochondrial genes totaling 21 loci and 12,538 characters) to infer a robust phylogeny for this terrestrial group.
Results strongly support the monophyly of all superfamilies; the superfamily status of Aradoidea and the
following relationships: (Aradoidea + (Pentatomoidea + (Coreoidea + (Lygaeoidea + Pyrrhocoroidea)))) in
Pentatomomorpha, and (Coreoidea + (Lygaeoidea + Pyrrhocoroidea)) in Eutrichophora. Our results suggest that
sampling greater numbers of genes is an effective tool for resolving phylogenetic problems.
Key words: taxonomy; Eutrichophora; molecular phylogenetics; molecular datasets
基于核糖体、Hox 和线粒体基因对蝽次目系统发育关系的重新分析(半翅目:异翅亚目)
李敏 1,2, 3,王艳会 1,谢强 1,田晓轩 1,4 ,李腾 1,张虎芳 3,卜文俊 1①
1. 南开大学生命科学学院,天津 300071;2. 太原师范学院生物系,山西 太原 030031;3. 山西农业
大学农学院,山西 太谷 030801;4. 天津中医药大学天津市现代中药重点实验室-省部共建国家重点实
验室(培育),天津 300193
摘要:蝽次目是蝽类昆虫(半翅目:异翅亚目)中物种多样性最为丰富的次目之一,其总科之间的系统
发育关系尚未确定,尤其对于真毛点类。以往对于蝽次目系统发育关系的研究基于形态数据,或仅基于
线粒体或核基因,或为 18S rDNA 和COI 基因片段。本研究选取了大量分子数据(18S rDNA、28S
rDNA、Hox 和线粒体基因,共 21 个基因位点,12538 bp),利用最大似然法和贝叶斯方法推断蝽次目总
科间的系统发育关系。研究结果支持各总科的单系性,扁蝽总科的总科级地位,蝽次目各总科系统发育
关系(扁蝽总科 +(蝽总科 +(缘蝽总科 +(长蝽总科 + 红蝽总科))))和真毛点类系统发育关系(缘
蝽总科 +(长蝽总科 + 红蝽总科))。研究表明选取多个基因作为分子标记为解决系统发育问题的有
Received 23 March 2016. Published 25 June 2016. Published online 21 June 2016.
① Corresponding author, E-mail: wenjunbu@nankai.edu.cn
82 LI et al. Reanalysis of the phylogenetic relationships of Pentatomorpha
效途径。
关键词:分类;真毛点类;分子系统发育;分子数据集
Introduction
Pentatomomorpha (Hemiptera: Heteroptera) includes some of the largest terrestrial bugs.
This group causes economic damage to many species of crops because its members are chiefly
phytophagous (Schuh & Slater 1995). Pentatomomorphan bugs comprise more than 14,000
species (Weirauch & Schuh 2011). This large number of species has been variously grouped
into four (Schaefer 1964), five (Štys 1961; Schaefer 1993; Schuh & Slater 1995), six (Henry
1997), or seven (Schuh 1986; Henry & Froeschner 1988) superfamilies (Table 1) and 39
families (Schuh & Slater 1995; Henry 1997; Rider 2015). The superfamily Aradoidea was
raised to the infraorder Aradomorpha by Sweet (1996, 2006) based on morphological evidence,
together with ecological, fossil and biogeographical evidence, but this viewpoint was treated
as invalid in Cassis and Schuh (2010).
Table 1. Previous classification systems of Pentatomomorpha
Schaefer 1964 Štys 1961 Schaefer 1993
Schuh & Slater
1995
Schuh1986
Henry &
Froeschner 1988
Henry 1997
Aradoidea Aradoidea Aradoidea Aradoidea Aradoidea
Pentatomoidea Pentatomoidea Pentatomoidea Pentatomoidea Pentatomoidea
Coreoidea Coreoidea Coreoidea Coreoidea Coreoidea
Pyrrhocoroidea Pyrrhocoroidea Pyrrhocoroidea
Lygaeoidea Lygaeoidea Lygaeoidea
Piesmatoidea Piesmatoidea Piesmatoidea
Idiostoloidea Idiostoloidea Idiostoloidea
Phylogenetic relationships proposed previously at the superfamily level within
Pentatomomorpha are summarized in Fig. 1. The works of Henry (1997) (Fig. 1A) and Yao et
al. (2012) (Fig. 1E) were based on morphological characters. Li et al. (2005) (Fig. 1B)
analyzed the phylogeny of Pentatomomorpha using molecular sequence data (partial 18S
rDNA and COI mtDNA) for the first time, and supported the monophyly of Pentatomomorpha,
a sister-group relationship of Aradoidea and Trichophora, and a clade comprising Lygaeoidea,
Coreoidea, and Pyrrhocoroidea. Later, several other works based on molecular sequences were
published (Hua et al. 2008; Tian et al. 2011; Yuan et al. 2015). Of these, Hua et al. (2008 )
(Fig. 1C) and Yuan et al. (2015) (Fig. 1D) were both mitochondrial genome sequences-based
analyses, but reached different relationships for the Eutrichophora, a concept proposed in Xie
et al. (2005). In Tian et al. (2011) (Fig. 1D), 26 pentatomomorphan and four cimicomorphan
putative families and six Hox gene fragments were analyzed, and the relationship (Aradoidea
+ (Pentatomoidea + (Lygaeoidea + (Coreoidea + Pyrrhocoroidea)))) was supported.
Until now, the monophyly of Pentatomomorpha, and the close relationship between
Pentatomoidea and Eutrichophora have been generally accepted (Henry 1997; Xie et al. 2005;
Li et al. 2005; Hua et al. 2008; Tian et al. 2011; Yao et al. 2012; Yuan et al. 2015). The
phylogenetic relationships within the Pentatomomorpha, however, have not reached full
Entomotaxonomia (2016) 38(2): 81–93 83
agreement. The unsolved problems are: 1) the monophyly of the superfamilies (Xie et al. 2005;
Li et al. 2005); 2) whether Aradoidea should be treated as superfamily Aradoidea or infraorder
Aradomorpha (Sweet 1996, 2006); and 3) the phylogenetic relationships among the
superfamilies, especially in the Eutrichophora (Fig. 1).
Previous phylogenetic studies were mostly based on morphological characters, fossils or
several single source genes. Lack of a sufficient number of features, errors in accessing
homology among features of various organisms, and convergent evolutionary changes can
affect the accuracy of phylogenetic inference (Boore & Brown 1999). Multiple genes with
increased sequence lengths are favorable for accurate phylogenetic analyses (Wortley et al.
2005; Kawahara et al. 2011; Corl & Ellegren 2013). Sequences used to construct phylogenetic
trees of Pentatomomorpha were less than 1.5 kb (Li et al. 2005), or just from nuclear (Tian et
al. 2011) or mitochondria (Hua et al. 2008; Yuan et al. 2015). In this study, we took advantage
of nucleotide and genomic sequencing projects to assemble a data set of 21 genes and 37
species. Based on this data set, a preliminary phylogenetic framework of Pentatomomorpha is
proposed and the status of Aradoidea is analyzed.
Figure 1. Proposed phylogenetic hypotheses within the Pentatomomorpha. A. Henry (1997); B. Li et al.
(2005); C. Hua et al. (2008); D. Tian et al. (2011), Yuan et al. (2015); E. Yao et al. (2012).
Material and methods
Taxa sampling. Taxa were chosen to represent the major lineages according to the
classification of Schuh and Slater (1995), Henry (1997) and Rider (2015). The terminal taxa
were decided according to Campbell and Lapointe (2009). In total, 37 species were sampled,
of which 28 species were ingroups. Two gerromorphans, two nepomorphans, two
leptopodomorphans and three cimicomorphans were selected as successive outgroups (Table 2).
Gene sampling. A total of 15 fragments of Hox genes of Ptilomera tigrina Uhlen, 1860,
Microvelia douglasi Scott, 1874, Tingis ampliata (Herrich-Schaeffer, 1838) and Anthocoris
thibetanus Poppius, 1909 were initially sequenced in this study. All of the other Hox genes, all
84 LI et al. Reanalysis of the phylogenetic relationships of Pentatomorpha
18S rDNA and 28S rDNA sequences, and most mitochondrial genome genes were generated
by authors of the previously published studies (Hua et al. 2008; Tian et al. 2011; Li et al. 2012;
Wang et al. 2015) (Table 2).
Laboratory procedure and primers. The laboratory procedures were the same as those
described in Li et al. (2012). Our newly amplified five Hox gene exons are: the 5’ end of
abd-A, 3’end of Dfd, 5’ end of Dfd, 5’ end of Ubx and 3’ end of pb, with sequence lengths
ranging from 450 to 1,000 bp. Universal degenerate primers were used in this study. Some are
from Tian et al. (2011) and Li et al. (2012); some were newly designed as shown in Table 3.
Table 3. Universal degenerate primers used in this study
Upper primer (5’-3) Lower primer (5’-3)
5’ end of abd-A Tian et al. 2011
3’ end of Dfd
Tian et al. 2011
5’ end of Dfd
Tian et al. 2011
5’ end of Ubx
Li et al. 2012
3’ end of pb
Tian et al. 2011
3’ end of pb for
Gerromorpha
GTRAA RGTYT GGTTY CAAAA YCG TAGTA YTCSG GCGYR AAGTC GTT
5’ end of abd-A
for Gerromorpha
TACCANCCNCAAATGCTGCC
TGTCWGGCATWGGTTGKGAAG
Sequences of 18S rDNA, 28S rDNA, Hox and 13 protein coding genes of mitochondrial
genomes were aligned by Muscle, embedded within MEGA 6.06 (Tamura et al. 2013) under
the default conditions of the program. The alignments of 18S rDNA and 28S rDNA were
adjusted according to the secondary structure models of heteropteran 18S rRNA and 28S
rRNA (Xie et al. 2009; Wang et al. 2013), and length-variable regions of both genes were
eliminated before phylogenetic reconstructing (Wang et al. 2015). The six Hox genes were
also manually corrected according Tian et al. (2011). Finally, 18S rDNA (1937 bp),28S rDNA
(4163 bp), six Hox genes (2708 aa) and 13 mitochondrial protein coding genes (PCGs) (3730
aa) were selected as molecular markers. All of the four data sets (18S rDNA, 28S rDNA, Hox
and mt-genomes) were concatenated in DAMBE V5.3.48 (Xia 2013).
Phylogenetic analyses were conducted utilizing maximum likelihood (ML) and Bayesian
inference (BI). Best-fitting substitution models were selected by Modeltest 3.7 (Posada &
Crandall 1998) and Treefinder (Jobb 2011) for the nrDNAs and amino acids, respectively. ML
analyses were conducted using RAxML 8.0.12 in PThreads version (Stamatakis 2014). The
most appropriate substitution model was GTR+G+I for 18S rDNA and 28S rDNA, mtArt + G
+ I for ATP6, COI, COII, COIII and ND4L, mtZOA + G + I for ATP8, CYTB, ND1, ND2,
ND3, ND4, ND5, ND6 and 5’ end of Scr, Dayhoff + G + I for 5’ end of abd-A and 3’ end of
Dfd, JTT + G + I for 5’ end of Dfd and 5’ end of Ubx, MIX [Dayhoff, BLOSUM, JTT, LG,
PMB, VT, WAG, mtREV, mtMam, mtArt, mtZOA, cpREV, rtREV, betHIV, witHIV, FLU] + G
+ I for 3’ end of pb. The best ML tree was calculated by RAxML with 100 runs, followed by
2000 bootstrap replicates. The consensus tree was built based on the best ML tree and the
bootstrap tree.
Bayesian analyses were conducted using MrBayes 3.2.5 (Ronquist et al. 2012) on CPUs
with the gaps in the matrices treated as missing. The GTR substitution model for 18S rDNA
Entomotaxonomia (2016) 38(2): 81–93 85
and 28S rDNA, and the mixed substitution model for the amino acids were used under a
gamma distribution (+G) with a proportion of invariable sites (+I) to account for among-site
rate variation. Four chains were used and the datasets were run for 5,000,000 generations until
the average standard deviation of split frequencies between the two-parallel Markov chain
Monte Carlo runs dropped below 0.01. Trees were sampled every 1000th cycle, with a burn-in
of 25% of sampled trees.
Phylogenetic results
Two trees generated through two methods (ML and BI) recovered an identical topology
for the higher-level nodes (Figs. 2, 3). The monophyly of all superfamilies except Aradoidea of
Pentatomomorpha was supported. Aradoidea was the sister group of the remaining
pentatomomorphan taxa, and then Pentatomoidea came out. The relationships among
Eutrichophora were as follows: (Coreoidea + (Lygaeoidea + Pyrrhocoroidea)). The suggested
relationships among the nine outgroups and Pentatomomorpha in this paper were (Gerridae +
Veliidae + (Nepomorpha + (Leptopodomorpha + (Cimicomorpha + Pentatomomorpha)))).
Figure 2. ML phylogram inferred from the combined data. Bootstrap support values are indicated at each node
(> 50%).
Two trees generated through two methods (ML and BI) recovered an identical topology
for the higher-level nodes (Figs. 2, 3). The monophyly of all superfamilies except Aradoidea of
Pentatomomorpha was supported. Aradoidea was the sister group of the remaining
86 LI et al. Reanalysis of the phylogenetic relationships of Pentatomorpha
pentatomomorphan taxa, and then Pentatomoidea came out. The relationships among
Eutrichophora were as follows: (Coreoidea + (Lygaeoidea + Pyrrhocoroidea)). The suggested
relationships among the nine outgroups and Pentatomomorpha in this paper were (Gerridae +
Veliidae + (Nepomorpha + (Leptopodomorpha + (Cimicomorpha + Pentatomomorpha)))).
ML tree (Fig. 2). Branch support for the monophyly of all the superfamilies in
Trichophora was robust (BP = 100%), as well as for the basal branch, and clade (Coreoidea +
(Lygaeoidea + Pyrrhocoroidea)). The support indices of clades (Pentatomoidea + (Coreoidea +
(Lygaeoidea + Pyrrhocoroidea))) and (Lygaeoidea + Pyrrhocoroidea) were 91% and 55%,
respectively.
Figure 3. BI phylogram inferred from the combined data. Bayesian posterior probabilities values are indicated
at each node (> 50%).
BI tree (Fig. 3). Bayesian topologies were identical to those recovered in the ML tree (Fig.
2). The posterior probability values supporting the monophyly of all superfamilies and the
relationships among the superfamilies were all 100%. In addition, inter-family relationships
were also well-resolved.
Discussion
Pentatomomorpha is recognized as a monophyletic group by most researchers. Our
results were in accordance with this opinion with high support values in ML and Bayesian
trees.
The monophyly of Pentatomoidea and the sister relationship of Pentatomoidea and the
other Trichophora were strongly supported in both analyses (BP > 90%; PP = 1.00), views that
corroborate almost all of the recent studies (Henry 1997; Li et al. 2005; Hua et al. 2008; Yao et
al. 2012; Yuan et al., 2015). Although only one representative was chosen for Aradoidea, no
Entomotaxonomia (2016) 38(2): 81–93 87
one has questioned its monophyly up to now. So the monophyly of Aradoidea is also accepted
in this research.
The superfamily Aradoidea was raised to the infraorder Aradimorpha by Sweet (1996,
2006) based on morphological, ecological, fossil and biogeographical evidences. In the
mitochondrial genome analyses of Hua et al. (2008), they found a unique gene rearrangement
in Aradoidea as its autapomorphy. However, Aradoidea was consistently placed as an
independent branch of Pentatomomopha in our research. The sister relationship of Aradoidea
and Trichophora (Leston 1958; Henry 1997; Li et al. 2005; Hua et al. 2008; Yao et al. 2012;
Yu a n et al. 2015) was further confirmed in this research.
As for the monophyly and relationships of the other three superfamilies i.e. Coreoidea,
Lygaeoidea and Pyrrhocoroidea, there were different views: Henry (1997) gave the topology
((Coreoidea + Pyrrhocoroidea) + Lygaeoidea) based on morphological characters. Li et al.
(2005) recognized them as a paraphyletic group based on partial 18S rDNA and COI
sequences; Xie et al. (2005) performed a Bayesian analysis using 18S rDNA and reported the
hypothesis of (Pyrrhocoroidea + (Coreoidea + Lygaeoidea)); Hua et al. (2008) used mitochondrial
genomes as molecular markers and obtained the same results as Xie et al. (2005). Tian et al.
(2011) analyzed six Hox genes and found the same relationship as Henry (1997). Yao et al.
(2012) clarified that the remainder of the Trichophora (except for Pentatomoidea) should be
assigned to 5 superfamilies (Coreoidea, Pyrrhocoroidea, Idiostoloidea, Lygaeoidea, and
Piesmatoidea) based on a combination of fossil and extant morphological characters, but
didn’t give the relationships among them. Phylogenetic analyses based on mitogenomic data
of Yuan et al. (2015) supported a phylogenetic relationship of (Lygaeoidea + (Pyrrhocoroidea
+ Coreoidea)), the same result found by Henry (1997) and Tian et al. (2011).
In this study, the relationship ((Coreoidea + (Lygaeoidea + Pyrrhocoroidea)) was
demonstrated. The three superfamilies are united by four synapomorphies with high
consistency indices: presence of an m-chromosome, typically (ground plan) three trichobothria
on at least some of the abdominal segments, 3-lobed (or more) salivary glands, and loss of
dorsal abdominal scent gland in nymphs between tergites 3/4. Further, several studies have
suggested that the incomplete abdominal sutures found in some taxa (e.g., Physopeltinae
[Largidae] and some Pyrrhocoridae) link these bugs with the Rhyparochrominae in the
Lygaeoidea (Henry 1997).
For the relationship among outgroups, (Nepomorpha + (Leptopodomorpha +
(Cimicomorpha + Pentatomomorpha))) was supported and this is in congruence with results of
previous analyses of morphological data (Schuh 1979) and molecular data (Wang et al. 2015).
In this paper, we focused primarily on the higher taxa of Pentatomomorpha. Inclusion of
wider sampling to establish the relationships among lower taxa would be the expected next
step.
Acknowledgements
We thank anonymous referees for their useful comments and suggestions. This research
was supported by the National Natural Science Foundation of China (31372240, 31501840,
31440078), the Postdoctoral Science Foundation of China (134845) and the Taiyuan Normal
University Training Program of Innovation and Entrepreneurship for Undergraduates (CXCY1610).
88 LI et al. Reanalysis of the phylogenetic relationships of Pentatomorpha
References
Boore JL & Brown WM. 1999. Big trees from little genomes: mitochondrial gene order as a phylogenetic tool.
Current Opinion in Genetics & Development, 8(6): 668–674.
Campbell V & Lapointe FJ. 2009. The use and validity of composite taxa in phylogenetic analysis. Systematic
Biology, 58(6): 560–572.
Cassis G & Schuh RT. 2010. Systematic methods, fossils, and relationships within Heteroptera (Insecta).
Cladistics, 26(3): 262–280.
Corl A & Ellegren H. 2013. Sampling strategies for species of trees: the effects on phylogenetic inference of
the number of genes, number of individuals, and whether loci are mitochondrial, sex-linked, or autosomal.
Molecular Phylogenetics and Evolution, 67(2): 358–366.
Henry TJ & Froeschner RC. 1988. Catalog of the Heteroptera, or true bugs, of Canada and the Continental
United States. E. J. Brill, Leiden, the Netherlands, 958 pp.
Henry TJ. 1997. Phylogenetic analysis of family groups within the infraorder Pentatomomorpha (Hemiptera:
Heteroptera), with emphasis on the Lygaeoidea. Annals of the Entomological Society of America, 90(3):
275–301.
Hua J, Li M, Dong P, Cui Y, Qiang X & Bu WJ. 2008. Comparative and phylogenomic studies on the
mitochondrial genomes of Pentatomomorpha (Insecta: Hemiptera: Heteroptera). BMC Genomics, 9: 610.
doi:10.1186/1471-2164-9-610.
Jobb G. 2011. TREEFINDER version of March 2011. Munich, Germany. http://www.treefinder.de/ (Accessed
on 18 March 2011).
Kawahara AY, Ohshima I, Kawakita A, Regier JC, Mitter C, Cummings MP, Davis DR, Wagner DL, Prins JD
& Lopez-Vaamonde C. 2011. Increased gene sampling strengthens support for higher-level groups within
leaf-mining moths and relatives (Lepidoptera: Gracillariidae). BMC Evolutionary Biology, 11: 182.
doi:10.1186/1471-2148-11-182.
Leston D. 1958. Chromosome number and the systematics of Pentatomomorpha (Hemiptera). Proceedings of
the 10th International Congress Entomology, 2: 911–918.
Li HM, Deng RQ, Wang JW, Chen ZY, Jia FL & Wang XZ. 2005. A preliminary phylogeny of the
Pentatomomorpha (Hemiptera: Heteroptera) based on nuclear 18S rDNA and mitochondrial DNA
sequences. Molecular Phylogenetics and Evolution, 37(2): 313–326.
Li M, Wang J, Tian XX, Xie Qiang, Liu HX & Bu WJ. 2012. Phylogeny of the true water bugs
(Hemiptera-Heteroptera: Nepomorpha) based on four Hox genes. Entomotaxonomia, 34(1): 35–44.
Posada D & Crandall KA. 1998. Modeltest: testing the model of DNA substitution. Bioinformatics, 14(9):
817–818.
Rider DA. 2015. Pentatomoidea home page. Available from: https://www.ndsu.edu/ndsu/rider/Pentatomoidea/.
(Accessed on 15 February 2016)
Ronquist F, Teslenko M, vander Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA &
Huelsenbeck JP. 2012. MrBayes3.2: efficient Bayesian phylogenetic inference and model choice across
a large model space. Systematic Biology, 61(3): 539–542.
Schaefer CW. 1964. The morphology and higher classification of the Coreoidea (Hemiptera-Heteroptera):
Parts I and II. Annals of the Entomological Society of America, 57: 670–684.
Schaefer CW. 1993. The Pentatomomorpha (Hemiptera: Heteroptera): an annotated outline of its systematic
history. European Journal of Entomology, 90: 105–122.
Schuh RT. 1979. Evolutionary trends in Heteroptera. Part II. Mouthpart-structures and feeding strategies, by R.
H. Cobben. Systematic Zoology, 28: 653–656.
Schuh RT. 1986. The influence of cladistics on heteropteran classification. Annual Review of Entomology,
Entomotaxonomia (2016) 38(2): 81–93 89
31(31): 67–93.
Schuh RT & Slater JA. 1995. True bugs of the world (Hemiptera: Heteroptera): classification and natural
history. Cornell University Press, Ithaca, United States, 336 pp.
Stamatakis A. 2014. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies.
Bioinformatics, 30: 1312–1313.
Štys P. 1961. Morphology of the abdomen and female ectodermal genitalia of the trichophorus Heteroptera
and bearing on their classification. Proceedings of the 11th International Congress of Entomology, 1:
37–43.
Sweet MH. 1996. Comparative external morphology of the pregenital abdomen of the Hemiptera. In: Schaefer CW
(Ed.), Studies on Hemiptera phylogeny. Lanham, Maryland, pp. 119–158.
Sweet MH. 2006. Justification for the Aradimorpha as an infraorder of the suborder Heteroptera (Hemiptera,
Prosorrhyncha) with special reference to the pregenital abdominal structure. Denisia, 19: 225–248.
Swofford DL. 2003. PAUP* 4 0b10: Phylogenetic Analysis Using Parsimony (*and other methods). Sinauer
Associates, Sunderland, Massachusetts.
Tamura K, Stecher G, Peterson D, Filipski A & Kumar S. 2013. MEGA6: molecular evolutionary genetics
analysis Version 6.0. Molecular Biology and Evolution, 30: 2725–2729.
Thompson JD, Higgins DG & Gibson TJ. 1994. CLUSTAL W: improving the sensitivity of progressive
multiple sequence alignment through sequence weighting, position-specific gap penalties and weight
matrix choice. Nucleic Acids Research, 22(22): 4673–4680.
Tian XX, Xie Q, Li M, Gao CQ, Cui Y, Xi L & Bu WJ. 2011. Phylogeny of pentatomomorphan bugs
(Hemiptera-Heteroptera: Pentatomomorpha) based on six Hox gene fragments. Zootaxa, 2888: 57–68.
Wang YH, Engel MS, Rafael JA, Dang K, Wu H, Wang Y, Xie Q & Bu WJ. 2013. A unique box in 28S rRNA
is shared by the enigmatic insect orders Zoraptera and Dictyoptera. PLoS ONE, 8(1): e53679.
doi:10.1371/journal.pone.0053679.
Wang YH, Cui Y, Rédei DA, Baňař P, Xie Q, Štys P, Damgaard J, Chen PP, Yi WB, Wang Y, Dang K, Li CR &
Bu WJ. 2015. Phylogenetic divergences of the true bugs (Insecta: Hemiptera: Heteroptera), with
emphasis on the aquatic lineages: the last piece of the aquatic insect jigsaw originated in the Late
Permian/Early Triassic. Cladistics, 4663(1): 404–407.
Weirauch C & Schuh RT. 2011. Systematics and evolution of Heteroptera: 25 years of progress. Annual
Review of Entomology, 56: 487–510.
Wortley AH, Rudall PJ, Harris DJ & Scotland RW. 2005. How much data are needed to resolve a difficult
phylogeny? Case study in Lamiales. Systematic Biology, 54(5): 697–709.
Xia X. 2013. DAMBE5: A comprehensive software package for data analysis in molecular biology and
evolution. Molecular Biology and Evolution, 30: 1720–1728.
Xie Q, Bu WJ & Zheng LY. 2005. The Bayesian phylogenetic analysis of the 18S rRNA sequences from the
main lineages of Trichophora (Insecta: Heteroptera: Pentatomomorpha). Molecular Phylogenetics and
Evolution, 34(2): 448–451.
Xie Q, Tian XX, Qin Y & Bu WJ. 2009. 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. Molecular Phylogenetics and Evolution, 50: 310–316.
Yao YZ, Ren D, Rider DA & Cai WZ. 2012. Phylogeny of the infraorder Pentatomomorpha based on fossil
and extant morphology, with description of a new fossil family from China. PLoS ONE, 7(5): e37289.
doi:10.1371/journal.pone.0037289.
Yuan ML, Zhang QL, Guo ZL, Wang J & Shen YY. 2015. The complete mitochondrial genome of Corizus
tetraspilus (Hemiptera: Rhopalidae) and phylogenetic analysis of Pentatomomorpha. PLoS ONE, 10(6):
e0129003. doi: 10.1371/journal.pone.0129003.
90 LI et al. Reanalysis of the phylogenetic relationships of Pentatomorpha
Table 2 Taxa and Genes Used in this Analysis
Infraorder-level Family-level Species Accession
number of
18SrDNA
Accession
number of
28SrDNA
Species
Gerromorpha Gerroidea:Gerridae
P
tilomera tigrina KJ461233 JF719825 Ptilomera tigrina
Gerromorpha Gerroidea:Veliidae Microvelia douglasi KJ461281 KJ461178 Entomovelia sp.
Nepomorpha Nepoidea:Nepidae Ranatra chinensis KJ461186 JF719824
L
accotrephes robustus
Nepomorpha Ochteroidea:Ochteridae Ochterus marginatus KJ461251 KJ461315 Ochterus marginatus
Leptopodomorpha Saldoidea:Saldidae Saldula saltatoria KJ461170 KJ461213 Saldula arsenjevi
Leptopodomorpha Leptopodoidea:Leptopodidae Valleriola buenoi KJ461310 KJ461290
L
eptopus sp.
Cimicomorpha Naboidea:Nabidae
A
lloeorhynchus sp. KJ461245 KJ461293
A
lloeorhynchus bakeri
Cimicomorpha Miroidea:Tingidae Corythucha ciliata KJ461201 KJ461175 Corythucha ciliata
Cimicomorpha Cimicoidea:Anthocoridae
A
nthocoris pilosus KJ461226 KJ461268 Orius niger
Pentatomomorpha Aradoidea:Aradidae Mezira sp. KJ461285 KJ461189
N
euroctenus parus
Pentatomomorpha Pyrrhocoroidea:Largidae
P
hysopelta gutta gutta KJ461164 KJ461255 Physopelta gutta
Pentatomomorpha Pyrrhocoroidea:Pyrrhocoridae
D
ysdercus cingulatus KJ461263 KJ461235
D
ysdercus cingulatus
Pentatomomorpha Pyrrhocoroidea:Pyrrhocoridae Macrocheraia grandis KJ461199 KJ461172
/
Pentatomomorpha Lygaeoidea:Lygaeidae
A
rocatus melanocephalus KJ461273 KJ461308 Kleidocerys resedae
Pentatomomorpha Lygaeoidea:Lygaeidae Sinorsillus piliferus KJ461269 KJ461194
Pentatomomorpha Lygaeoidea:Artheneidae
A
rtheneis intricata KJ461299 KJ461203
/
Pentatomomorpha Lygaeoidea:Blissidae Cavelerius excavatus KJ461294 KJ461237
/
Pentatomomorpha Lygaeoidea:Cymidae Cymus koreanus KJ461232 KJ461312
/
Pentatomomorpha Lygaeoidea:Geocoridae Geocoris ater KJ461218 KJ461309 Geocoris pallidipennis
Pentatomomorpha Lygaeoidea:Piesmatidae
P
iesma maculatum KJ461171 KJ461169
/
Pentatomomorpha Lygaeoidea:Rhyparochromidae
N
eolethaeus assamensis KJ461280 KJ461225
Pentatomomorpha Coreoidea:Rhopalidae Stictopleurus punctatonervosus KJ461217 KJ461286
/
Pentatomomorpha Coreoidea:Alydidae Megalotomus ornaticeps KJ461221 KJ461236 Stictopleurus subviridis
Pentatomomorpha Coreoidea:Alydidae Riptortus pedestris AB725684 AB725684
/
Pentatomomorpha Coreoidea:Coreidae Cletus punctiger KJ461173 KJ461219 Riptortus pedestris
Pentatomomorpha Coreoidea:Stenocephalidae
D
icranocephalus alticolus KJ461228 KJ461267 Hydaropsis longirostris
Pentatomomorpha Pentatomoidea:Tessaratomidae Eurostus validus KJ461181 JF719828
/
Pentatomomorpha Pentatomoidea:Pentatomidae Rhaphigaster nebulosa X89495 EU426880 Eusthenes cupreus
Pentatomomorpha Pentatomoidea:Pentatomidae Eurydema maracandica JX997807 JX997806 Halyomorpha halys
Pentatomomorpha Pentatomoidea:Acanthosomatidae Elasmucha laeviventris KJ461262 KJ461208 Eurydema gebleri
Pentatomomorpha Pentatomoidea:Scutelleridae Cantao ocellatus KJ461182 KJ461230
/
Pentatomomorpha Pentatomoidea:Cydnidae Fromundus pygmaeus KJ461296 KJ461209
/
Pentatomomorpha Pentatomoidea:Plataspidae Coptosoma bifarium KJ461259 KJ461239 Macroscytus gibbulus
Pentatomomorpha Pentatomoidea:Plataspidae
P
aracopta maculata KJ535895 KJ535890 Coptosoma bifaria
Pentatomomorpha Pentatomoidea:Urostylididae Urochela luteovaria KJ461205 KJ461306 Megacopta cribraria
Pentatomomorpha
Pentatomomorpha
Pentatomoidea:Urostylididae
Pentatomoidea:Urostylididae
Tessaromerus licenti
Urostylis chinai
KJ461291
KJ461275
KJ461300
KJ461287
Urochela quadrinotata
/
Entomotaxonomia (2016) 38(2): 81–93 91
Continued Table 2 on the right
Accession
number of Mt
genomes
Species Accession
number of
5'end abd-A
3'end Dfd 5'end Dfd 3'end pb 5'end Scr 5'end Ubx
KP400583 Ptilomera tigrina 0 0 0 0
/
/
KP400582 Microvelia douglasi 0
/
0 0
/
/
NC_012817 Ranatra chinensis JQ409392 JQ409405 JQ409432 JQ409445
/
JQ409419
NC_012820 Ochterus marginatus JQ409394 JQ409407 JQ409434 JQ409447
/
JQ409421
NC_012463 Saldula pallipes JQ409389 JQ409402 JQ409429 JQ409442
/
JQ409416
FJ456946 Valleriola sp. JQ409390 JQ409403 JQ409430 JQ409443
/
JQ409417
NC_016432 Gorpis annulatus FJ851727 FJ851754 FJ851784 FJ851810 FJ851840 FJ851870
KC756280 Tingis ampliata 0 0 0 0
/
0
EU427341
A
nthocoris thibetanus 0 0 0
/
FJ851839 FJ851868
NC_012459 Carventus hainanensis FJ851698 FJ851729 FJ851756 FJ851785 FJ851811 FJ851842
NC_012432 Physopelta quadriguttata FJ851712 FJ851739 FJ851769 FJ851798 FJ851825 FJ851856
NC_012421 Pyrrhopeplus posthumus FJ851711 FJ851738 FJ851768 FJ851797 FJ851824 FJ851855
/
/
/
/
/
/
/
/
KJ584365
L
ygaeus equestris FJ851700 FJ851731 FJ851758 FJ851787 FJ851813 FJ851844
Oncopeltus fasciatus FJ851728 FJ851755 / / FJ851841 FJ851871
/
A
rtheneis intricata FJ851705 FJ851734 FJ851762 FJ851792 FJ851818 FJ851849
/
D
imorphopterus spinolae FJ851703 FJ851733 FJ851760 FJ851790 FJ851816 FJ851847
/ Cymus koreanus FJ851701 / / FJ851788 FJ851814 FJ851845
NC_012424 Geocoris pallidipennis FJ851704 / FJ851761 FJ851791 FJ851817 FJ851848
/ Piesma josifovi FJ851710 FJ851737 FJ851767 / FJ851823 FJ851854
/
D
ieuches femoralis FJ851707 / FJ851764 FJ851794 FJ851820 FJ851851
NC_012888
A
eschyntelus chinensis FJ851715 FJ851742 FJ851772 FJ851801 FJ851828 FJ851859
/
/
/
/
/
/
/
/
NC_012462 Riptortus pedestris FJ851713 FJ851740 FJ851770 FJ851799 FJ851826 FJ851857
NC_012456
A
canthocoris scabe
r
FJ851714 FJ851741 FJ851771 FJ851800 FJ851827 FJ851858
/
/
/
/
/
/
/
/
NC_022449 Eurostus validus FJ851722 FJ851749 FJ851779 FJ851808 FJ851835 FJ851866
NC_013272
N
ezara viridula FJ851723 FJ851750 FJ851780 FJ851809 FJ851836 FJ851867
KP207595
/
/
/
/
/
/
/
/
A
canthosoma crassicauda FJ851717 FJ851744 FJ851774 FJ851803 FJ851830 FJ851861
/ Poecilocoris lewisi FJ851720 FJ851747 FJ851777 FJ851806 FJ851833 FJ851864
NC_012457
A
ethus indicus FJ851718 FJ851745 FJ851775 FJ851804 FJ851831 FJ851862
NC_012449
/
/
/
/
/
/
/
NC_015342 Megacopta cribraria FJ851719 FJ851746 FJ851776 FJ851805 FJ851832 FJ851863
JQ743678 Urochela distincta FJ851716 FJ851743 FJ851773 FJ851802 FJ851829 FJ851860
/
/
/
/
/
/
/
/
/
/
/
/
/
/
/
/
Note: Sequences marked with 0 were newly amplified in this research and the accession numbers will
be supplied after assigned by the GenBank staff.