# Sequencing of the mitochondrial genome of the avocado lace bug Pseudacysta perseae (Heteroptera, Tingidae) using a genome skimming approach

Arthur Kocher, Eric Guilbert, Emeline Lhuillier, Jérôme Murienne

► To cite this version:

Arthur Kocher, Eric Guilbert, Emeline Lhuillier, Jérôme Murienne. Sequencing of the mitochondrial genome of the avocado lace bug Pseudacysta perseae (Heteroptera, Tingidae) using a genome skimming approach. Comptes Rendus. Biologies, 2015, 338 (3), pp.149-160. 10.1016/j.crvi.2014.12.004. hal-02636710

HAL Id: hal-02636710

https://hal.inrae.fr/hal-02636710v1

Submitted on 27 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Genetics/Génétique

# Sequencing of the mitochondrial genome of the avocado lace bug *Pseudacysta perseae* (Heteroptera, Tingidae) using a genome skimming approach

<figure><img src="image_2.png" /></figure>

CrossMark

## Séquençage du génome mitochondrial du tigre de l'avocatier Pseudacysta perseae (Heteroptera, Tingidae) par une approche d'écrémage de génome

Arthur Kocher<sup>a,*</sup>, Éric Guilbert<sup>b</sup>, Émeline Lhuillier<sup>c</sup>, Jerôme Murienne<sup>a</sup>

## ARTICLE INFO

**Article history:**
Received 1st September 2014
Accepted after revision 10 December 2014
Available online 27 January 2015

**Keywords:**
French Guiana
Illumina Hiseq
Mitogenome
Phylogenomics
Shotgun sequencing

**Mots clés :**
Guyane
Illumina Hiseq
Mitogenome
Phylogénomique
Séquençage aléatoire

## ABSTRACT

Lace bugs (Tingidae) are a family of phytophagous heteropterans, some of which are important agricultural and forestry pests. They currently comprise around 2500 species distributed worldwide, for which only one mitochondrial genome has been described so far. We sequenced the complete mitochondrial genome and the nuclear ribosomal gene segment of the avocado lace bug *Pseudacysta perseae* using a genome skimming approach on an Illumina Hiseq 2000 platform. Fifty-four additional heteropteran mitogenomes, including the one of the sycamore lace bug *Corythucha ciliata*, were retrieved to allow for comparisons and phylogenetic analyses. *P. perseae* mitochondrial genome was determined to be 15,850 bp long, and presented the typical organisation of insect mitogenomes. The phylogenetic analysis placed *P. perseae* as a sister to *C. ciliata* but did not confirm the monophyly of Miroidae including Tingidae. Our results contradicted widely accepted phylogenetic hypothesis, which highlights the limits of analyses based on mitochondrial data only. Shotgun sequencing approaches should provide substantial improvements in harmonizing mitochondrial and nuclear databases.

© 2014 Académie des sciences. Published by Elsevier Masson SAS. All rights reserved.

## RÉSUMÉ

Les tingidés (Tingidae) sont une famille d'hétéroptères phytophages, dont certains sont d'importants ravageurs agricoles et forestiers. Ils comprennent actuellement environ 2500 espèces distribuées à travers le monde, pour lesquelles un seul génome mitochondrial a été décrit. Nous avons séquencé le génome mitochondrial complet et les gènes ribosomiques nucléaires du tigre de l'avocatier *Pseudacysta perseae* avec une approche d'écrémage de génome sur une plateforme Illumina Hiseq 2000. Cinquante-quatre autres

<footnote><sup>a</sup>Laboratoire EDB (UMR 5174, CNRS/UPS/ENFA), université Paul-Sabatier, bâtiment 4R1, 118, route de Narbonne, 31062 Toulouse cedex 9, France</footnote>

<footnote><sup>b</sup>UMR 7205 CNRS, CP50, département de systématique et évolution, Muséum national d'histoire naturelle, 57, rue Cuvier, 75005 Paris, France</footnote>

<footnote><sup>c</sup>INRA, UAR 1209 département de génétique animale, GeT-PlaGe, Campus INRA, Auzeville, 31326 Castanet-Tolosan, France</footnote>

<footnote>* Corresponding author.</footnote>

<footnote>E-mail address: arthur.kocher@gmail.com (A. Kocher).</footnote>
mitogénomes d'hétéroptères, dont celui du tigre du platane *Corythucha ciliata*, ont été récupérés pour permettre des comparaisons et une analyse phylogénétiques. Le génome mitochondrial de *P. perseae* est long de 15 850 paires de bases et présente l'organisation typique des mitogénomes d'insectes. L'analyse phylogénétique place *P. perseae* comme groupe frère de *C. ciliata*, mais ne confirme pas la monophylie de Miroidae, incluant Tingidae. Notre analyse contredit des hypothèses phylogénétiques largement acceptées, ce qui souligne les limites des analyses basées uniquement sur de l'ADN mitochondrial. Les approches de séquençage aléatoire devraient apporter d'importantes améliorations dans l'harmonisation des bases de données mitochondriales et nucléaires.

© 2014 Académie des sciences. Publié par Elsevier Masson SAS. Tous droits réservés.

# 1. Introduction

Mitochondrial DNA has been increasingly used for comparative genomic and phylogenetic studies in the past decade. It has various interesting properties such as abundance in animal tissues, small size, simple genomic structure, fast rate of evolution (in animals) and exclusive maternal inheritance with low level of recombination [1,2]. Until recently, the sequencing of complete mitochondrial genomes has been relying on long range Polymerase Chain Reactions (PCR), which are difficult to perform and time consuming. The immense yield now provided by Next Generation Sequencing (NGS) allows alleviating this necessity. The high copy number of mitochondrial genomes per cell allows the use of a relatively low sequencing capacity compare to what is necessary to recover nuclear genomes. This "genome skimming" approach, originally developed for plant organelles assembly [3] has been successfully used to assemble a wide variety of animal mitochondrial genomes [4-7]. However, mitochondrial studies have been unbalanced among taxa, and the amount of available data for insects is still small compared to that for vertebrates [8,9].

In most metazoans, including insects, the mitogenome is a circular double-stranded molecule of 14-20 kb in size, which contains a set of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs) and two ribosomal RNA genes (rRNAs) [10,11]. However, variations in gene content exist [8,12-14]. It also exhibits a large non-coding region known as the control region (CR), which is implicated in the initiation of transcription and replication processes [15-18].

Cimicomorpha is the largest infraorder of Heteroptera with more than 20,000 species currently placed in 17 families [19]. Tingidae is a family of Cimicomorpha distributed worldwide. They are sister group of the Miridae into the super-family Miroidea, according to the last phylogenetical analysis [20,21]. They are commonly named "lace bugs" because of the highly reticulate lace-like patterning of the pronotum and forewings of adults. All Tingidae are phytophagous, and can be destructive to plants. Some of them are important economical pests in agriculture and forestry [22]. Tingidae is currently comprised of approximately 2500 species [23] for which only one mitochondrial genome has been described so far (the sycamore lace bug, *Corythucha ciliata*; [24]). The avocado lace bug, *Pseudacysta perseae* (Heidemann, 1908) was considered a minor pest of avocado for several years

[25] but since the mid-1990s, population outbreaks of *P. perseae* have been observed in Florida and the Caribbean. It is now considered as a serious threat for avocado in these regions [26] and several studies have focused on biological and chemical solutions for its management [27-31]. The known geographic distribution of *P. perseae* includes southern states of the USA (Florida, Georgia, Louisiana, Texas and, more recently, California) as well as various locations in the Caribbean, Mexico and Guatemala [30]. In South America, it has recently been reported in Venezuela and French Guiana [32,33].

In this paper, we use a genome skimming approach based on Illumina technology to assemble the complete mitochondrial genome of *P. perseae*. Its organization and features are compared to the mitogenome of *C. ciliata*. A set of 55 heteropteran mitochondrial genomes was used to perform a phylogenetic analysis.

# 2. Material and method

## 2.1. Specimen, DNA extraction and sequencing

Specimens of *P. perseae* were collected on an avocado in Remire-Montjoly, French Guiana on May 6th 2009. Total genomic DNA was extracted from leg muscle tissues using the DNeasy Blood and Tissue kit (Qiagen, Valencia, CA, USA) following a protocol adapted from the manufacturer's instructions. The quality and quantity of extracted genomic DNA was evaluated using the NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and the PicoGreen double stranded DNA quantitation assay kit (Life Technologies, Carlsbad, CA, USA).

The genomic DNA was sent for library construction and sequencing to the GeT-PlaGe core facilities of Genotoul (Toulouse, France). One hundred and sixty-two ng of DNA were used for library construction using the Illumina TruSeq Nano DNA Sample Prep Kit (Illumina Inc. San Diego, CA) following the instructions of the supplier. After shearing by ultrasonication with a Covaris M220 (Covaris Inc, Woburn, MA), purified fragments were A-tailed and ligated to sequencing indexed adapters. Fragments with an insert size around 450 bp were selected with Agencourt Ampure XP beads (Beckman Coulter, Inc.), and enriched with 8 cycles of PCR before library quantification and validation. The library was multiplexed with 23 other libraries (generated in other projects). The pool of libraries was then hybridized on one lane of Illumina Hiseq
2000 flow cell using the Illumina TruSeq PE Cluster Kit v.3, and paired-end reads of 100 nucleotides were collected on the HiSeq 2000 sequencer using the Illumina TruSeq SBS Kit v.3 (200 cycles). Quality filtering was performed by the Consensus Assessment of Sequence and Variation (CASAVA) pipeline (Illumina Inc., San Diego, CA, USA). Sequence data were stored on the NG6 platform [34] and all the computations were performed on the computer cluster of the Genotoul bioinformatic platform (Toulouse, France).

## 2.2. Sequence assembly

The mitochondrial genome was assembled using a previously used strategy [35–37] in essence similar to that of Hahn et al. [38]. Reads aligning with the mitochondrial protein sequences of *C. ciliata* [24] were retrieved using the PLAST program [39]. Reads with a match of at least 90% were assembled into contigs using the Velvet assembler [40] with a k-mer length of 69 and all the remaining parameters left at their default values. The resulting contigs were used as seeds to initiate the genome walking approach using the extractreads2 program (included in the Obitools package; http://metabarcoding.org/obitools). Reads sharing at least 80 consecutive bp with the seeds were selected and subsequently used as seed to repeat the operation until no new read was identified. The newly selected reads were assembled with the Velvet assembler. The few resulting contigs were assembled using Geneious 6.0.6 Pro (Biomatters, Auckland, New Zealand). Additionally, reads aligning with the 28S and 18S rRNA genes of *Eurydema maracandica* [41] were identified and the sequence of the complete nuclear ribosomal gene segment was recovered with the same method.

## 2.3. Genome annotation and sequence analysis

Coverage statistics were computed with Geneious 6.0.6 Pro (Biomatters, Auckland, New Zealand) by mapping the Illumina reads to the assembled genome. Reads were set as paired and mapped to the reference sequence using a custom sensitivity. The following mapping parameter were used: a minimum overlap of 10 bp, a minimum overlap identity of 90%, a maximum mismatch per read of 10%, a maximum ambiguity of four bases and no gaps allowed.

The mitochondrial genome was first annotated using the MITOS web server [42] applying the invertebrate mitochondrial genetic code (NCBI code table 5). The annotations of tRNA genes were kept unchanged. The annotations of protein-coding genes (PCGs) were refined by checking manually for consistent start/stop codons and reading frames. The annotations of rRNA genes were extended until adjacent tRNAs when existing, following the punctuation model of mtDNA transcription [43,44]. Annotation of srRNA's 3' end was further adjusted by mapping the srRNA of *C. ciliata* for which the secondary structure was predicted and presented all expected domains [24]. We used Geneious 6.0.6 Pro with the following parameters: a word length of seven, maximum gap size of 15 and maximum mismatches of 40%. This approach conducted to extend IrRNA's annotation by 13 bp until tRNA<sub>Leu</sub> at its 3' end and by

194 bp until tRNA<sub>Val</sub> at its 5' end. srRNA's annotation was extended by 1 bp until tRNA<sub>Val</sub> at its 3' end and by 35 bp at its 5' end. We further verified the consistency of these new annotations by mapping the rRNAs of *C. ciliata* as described above. The large remaining non-annotated sequence was annotated as the control region in concordance with other insect mitogenomes. The tandem repeats copy number was inferred using a coverage analysis.

The 18S and 28S rRNAs were annotated in comparison with that of *E. maracandica* [41]. The 5.8S rRNA was annotated in comparison with that of *Triatoma dimidiata* (accession number: KF142517).

## 2.4. Sequence analyses and phylogenetics

Base composition and codon usage where computed with MEGA6 [45]. AT-skew ([A-T]/[A+T]) and GC-skew ([G-C]/[G+C]) where used to measure nucleotide compositional differences between genes [46]. Relative synonymous codon usage (RSCU) where used to described bias in synonymous codon composition. Tandem repeats were identified using Tandem Repeat Finder web server [47]. tRNA's secondary structure were inferred via the MITOS web server pipeline. We looked for putative secondary structures within the 100 bp of the control region flanking tRNA<sub>Ile</sub> where stem loops structures have already been reported in Tingidae and other heteropterans [24,36,48,49].

A set of 54 additional heteropteran mitogenomes was downloaded from Genbank (Table 1). Two mitogenomes of Auchenorrhyncha were used as outgroups. The 13 PCGs were used for the analysis to allow for comparison with previous studies [24,50]. They were first aligned separately based on amino-acids translation with translatorX [51]. Divergent regions were removed with Gblocks.0.91b before back-translation in order to conserve reading frames. All resulting alignments were then concatenated using FASconCAT [52]. The best partitioning scheme and substitution model were inferred with PartitionFinder.1.1.1 [53], using the greedy algorithm for scheme search and the Bayesian information criterion for scheme selection. A maximum-likelihood (ML) analysis was performed with RAxML 8.0 [54] and nodal support was estimated using a rapid bootstrap procedure [55] with the majority-rule tree based bootstopping criteria. A Bayesian analysis was conducted using Mr. Bayes 3.2 [56], starting from four random trees with ten Markov chains (nine heated chain and one cold chain), 2,000,000 generations and all other parameters set to default. Each set was sampled every 200 generations with a burn-in of 25% of the sampled trees. At the end of the analysis, the average standard deviation of split frequencies was below the recommended 0.01.

# 3. Results and discussion

## 3.1. Genome sequencing, assembly and annotation

After filtering 5.08% of the initial reads, raw sequence data represent a total of 6,965,443 paired-end reads (13,930,886 reads in total). Fifty thousands and eight
<table><caption>**Table 1**
Complete or near complete mitochondrial genomes used in this study.</caption><thead><tr><th>Suborder</th><th>Infra-order/Superfamily</th><th>Family</th><th>Species</th><th>Accession num.</th><th>Reference</th></tr></thead><tbody><tr><td rowspan="2">Auchenorryncha</td><td rowspan="2">Fulgoroidea</td><td>Fulgoridae</td><td>Lycorma delicatula</td><td>EU909203</td><td>[83]</td></tr><tr><td>Flatidae</td><td>Geisha distinctissima</td><td>NC_012617</td><td>[70]</td></tr><tr><td rowspan="10">Heteroptera</td><td>Cimicomorpha</td><td></td><td>Orius niger</td><td>EU427341</td><td>[74]</td></tr><tr><td>Cimicoidea</td><td>Anthocoridae</td><td>Adelphocoris fasciaticollis</td><td>NC_023796</td><td>[84]</td></tr><tr><td rowspan="4">Miroidea</td><td rowspan="4">Miridae</td><td>Apolygus lucorum</td><td>NC_023083</td><td>[85]</td></tr><tr><td>Lygus lineolaris</td><td>NC_021975</td><td>Unpublished</td></tr><tr><td>Nesidiocoris tenuis</td><td>NC_022677</td><td>[86]</td></tr><tr><td>Corythucha ciliata</td><td>NC_022922</td><td>[24]</td></tr><tr><td rowspan="6">Naboidea</td><td rowspan="6">Nabidae</td><td>Pseudacysta perseae</td><td></td><td>This study</td></tr><tr><td>Alloehrhynchus bakeri</td><td>HM235722</td><td>[87]</td></tr><tr><td>Gorpis annulatus</td><td>JF907591</td><td>[87]</td></tr><tr><td>Gorpis humeralis</td><td>JF927830</td><td>[87]</td></tr><tr><td>Himacerus apterus</td><td>JF927831</td><td>[87]</td></tr><tr><td>Himacerus nodipes</td><td>JF927832</td><td>[87]</td></tr><tr><td rowspan="5">Reduviidea</td><td rowspan="5">Reduviidae</td><td>Nabis apicalis</td><td>JF907590</td><td>[87]</td></tr><tr><td>Agriosphodrus dohrni</td><td>NC 015842</td><td>[50]</td></tr><tr><td>Oncocephalus breviscutum</td><td>NC 022816</td><td>[88]</td></tr><tr><td>Sirthenea flavipes</td><td>HQ645959</td><td>[49]</td></tr><tr><td>Triatoma dimidiata</td><td>NC 002609</td><td>[48]</td></tr><tr><td></td><td></td><td>Valentia hoffmanni</td><td>NC 012823</td><td>[69]</td></tr><tr><td>Enicocephalomorpha</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Enicocephaloidea</td><td>Enicocephalidae</td><td></td><td>Stenopirates sp.</td><td>NC_016017</td><td>[89]</td></tr><tr><td>Gerromorpha</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gerroidea</td><td>Gerridae</td><td></td><td>Aquarius paludum</td><td>NC_012841</td><td>[69]</td></tr><tr><td>Hydrometroidea</td><td>Hydrometridae</td><td></td><td>Hydrometra sp.</td><td>NC_012842</td><td>[69]</td></tr><tr><td>Leptopodomorpha</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Leptopodoidea</td><td>Leptopodidae</td><td></td><td>Leptopus sp.</td><td>FJ456946</td><td>[69]</td></tr><tr><td>Saldoidae</td><td>Saldidae</td><td></td><td>Saldula arsenjevi</td><td>EU427345</td><td>[74]</td></tr><tr><td>Nepomorpha</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Corixoidea</td><td>Corixidae</td><td></td><td>Sigara septemlineata</td><td>FJ456941</td><td>[69]</td></tr><tr><td rowspan="2">Naucoroidea</td><td>Aphelocheiridae</td><td></td><td>Aphelocheirus ellipsoides</td><td>FJ456939</td><td>[69]</td></tr><tr><td>Naucoridae</td><td></td><td>Ilyocoris cimicoides</td><td>NC_012845</td><td>[69]</td></tr><tr><td rowspan="2">Nepoidea</td><td>Belostomatidae</td><td></td><td>Diplonychus rusticus</td><td>FJ456940</td><td>[69]</td></tr><tr><td>Nepidae</td><td></td><td>Laccotrephes robustus</td><td>NC_012817</td><td>[69]</td></tr><tr><td>Notonectoidea</td><td>Notonectidae</td><td></td><td>Enithares tibialis</td><td>NC_012819</td><td>[69]</td></tr><tr><td rowspan="2">Ochteroidea</td><td>Gelastocoridae</td><td></td><td>Neothra sp.</td><td>NC_012838</td><td>[69]</td></tr><tr><td>Ochteridae</td><td></td><td>Ochterus marginatus</td><td>NC_012820</td><td>[69]</td></tr><tr><td>Pleoidea</td><td></td><td></td><td>Helotrephes sp.</td><td>FJ456951</td><td>[69]</td></tr><tr><td>Pentatomomorpha</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Aradoidea</td><td>Aradidae</td><td></td><td>Aradacanthia heissi</td><td>HQ441233</td><td>[90]</td></tr><tr><td rowspan="6">Coreoidea</td><td rowspan="2">Alydidae</td><td></td><td>Brachyrhynchus hsiaoii</td><td>NC_022670</td><td>[91]</td></tr><tr><td></td><td>Neuroctenus parus</td><td>EU427340</td><td>[74]</td></tr><tr><td rowspan="2">Coreidae</td><td></td><td>Riptortus pedestris</td><td>EU427344</td><td>[74]</td></tr><tr><td></td><td>Hydaropsis longirostris</td><td>EU427337</td><td>[74]</td></tr><tr><td rowspan="5">Rhopalidae</td><td></td><td>Aeschyntelus notatus</td><td>EU427333</td><td>[74]</td></tr><tr><td></td><td>Stictopleurus subviridis</td><td>NC_012888</td><td>[69]</td></tr><tr><td>Berytidae</td><td></td><td>Yemmalysus parallelus</td><td>EU427346</td><td>[74]</td></tr><tr><td>Colobathristidae</td><td></td><td>Phaenacantha marcida</td><td>EU427342</td><td>[74]</td></tr><tr><td>Geocoridae</td><td></td><td>Geocoris pallidipennis</td><td>EU427336</td><td>[74]</td></tr><tr><td rowspan="5">Lygaeoidea</td><td></td><td></td><td>Chauliops fallax</td><td>NC_020772</td><td>[74]</td></tr><tr><td></td><td></td><td>Malcus inconspicus</td><td>EU427339</td><td>[74]</td></tr><tr><td>Cydnidae</td><td></td><td>Macroscytus gibbulus</td><td>EU427338</td><td>[74]</td></tr><tr><td>Dinidoridae</td><td></td><td>Coridius chinensis</td><td>JQ739179</td><td>[92]</td></tr><tr><td>Pentatomidae</td><td></td><td>Dolycoris baccarum</td><td>NC_020373</td><td>[93]</td></tr><tr><td rowspan="2">Pentatomoidea</td><td></td><td></td><td>Halyomorpha halys</td><td>NC_013272</td><td>[94]</td></tr><tr><td></td><td></td><td>Nezara viridula</td><td>NC_011755</td><td>[74]</td></tr><tr><td rowspan="5">Plataspidae</td><td></td><td></td><td>Coptosoma bifaria</td><td>EU427334</td><td>[74]</td></tr><tr><td></td><td></td><td>Megacopta cribraria</td><td>NC_015842</td><td>[73]</td></tr><tr><td>Tessaratomidae</td><td></td><td>Eusthenes cupreus</td><td>NC_022449</td><td>[95]</td></tr><tr><td>Urostylididae</td><td></td><td>Urochela quadrinotata</td><td>NC_020144</td><td>[96]</td></tr><tr><td>Largidae</td><td></td><td>Physopelta gutta</td><td>EU427343</td><td>[74]</td></tr><tr><td rowspan="2">Pyrrhocoroidea</td><td></td><td></td><td>Dysdercus cingulatus</td><td>EU427335</td><td>[74]</td></tr><tr><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
hundred reads were assembled into a 15,850-bp circular sequence (Genbank accession number: KM278221), representing the complete mitochondrial genome with an average sequencing depth of 320.5. A circular map of the mitogenome and the assembly coverage are presented on Fig. 1. The typical set of 37 genes (13 PCGs, 22 tRNAs, two rRNAs) and one control region were identified. Twenty-three genes are encoded on the majority strand and the others located on the minority strand. Eight gene overlaps where observed, the longest being a 14-bp region between *ATP6* and *COX3* (Table 2). Apart from the control region, nine non-coding regions ranging from 1 bp to 19 bp where identified.

The complete nuclear ribosomal sequence was recovered. A total of 45,600 reads were assembled into a 7588-bp sequence comprising of 18S rRNA (1899 bp), ITS1 (232 bp), 5.8S rRNA (155 bp), ITS2 (1220 bp) and 28S rRNA (4082 bp). The sequence was deposited on Genbank with the accession number KM278220.

*P. perseae* mitogenome shares the same architecture and orientation as *C. ciliata*. This gene arrangement is shared with many insects [57,58]. It is also found in crustaceans and was the first to be determined differing by a single tRNA translocation from that of the Chelicerata *Limulus polyphemus*, considered as ancestral for Arthropoda [59,60]. Therefore, this mitogenome organisation is thought to be ancestral for the insect-crustacean clade [10,61]. *P. perseae* mitogenome is 593 longer than *C. ciliata* mitogenome which is almost entirely due to a variation in

the size of the control region (1426 bp in *P. perseae* vs 818 in *C. ciliata*). More generally, variations in insect mitogenomes size are mainly due to the control region, which ranges from 70 bp in *Ruspolia dubia* (Orthoptera) [62] to 4599 bp in *Drosophila melanogaster* (Diptera) [62,63].

## 3.2. Protein-coding genes

The total length of the 13 PCGs was 10,992 bp. Their nucleotide composition is strongly biased toward AT with an overall AT content of 76.4% (Table 3). All PCGs have an ATN start codon (Table 2). Six PCGs initiated with ATT (COX2, COX3, ND3, ND5, ND6 and CYTB), four initiated with ATA (ND2, ATP8, ND4L and ND1) and three initiated with ATG (COX1, ATP6 and ND4). Six genes initiated with the same ATG start codon between both lace bugs mitogenomes (*COX1*, *ATP8*, *ATP6*, *ND4*, *ND6* and *ND1*). Other unconventional start codons were found in insects such as TTG in heteropterans [24], CGA and TTAG in lepidopterans [64,65], or ATAA, GTAA and TTAA in dipterans [57,66] but that was not the case in *P. perseae*.

The majority of PCGs have a usual TAA stop codon, but four T and one TA stop codons were identified (ND2, COX2, COX3, CYTB and ND5 respectively). These incomplete stop codons are immediately adjacent to tRNA genes encoded on the same strand which is consistent with the punctuation model for primary transcripts processing followed by 3' polyadenylation of mature mRNA that will allow the

<figure><img src="image_3.png" /><figcaption>**Fig. 1.** (Color online.) Schematic representation of *Pseudacysta perseae* mitogenome. tRNAs are labelled according to the IUPAC-IUB single-letter amino acid codes. Arrows indicate directions of genes (clockwise = majority strand).</figcaption></figure>
<caption>**Table 2**</caption>

Summary of the mitochondrial genome of *Pseudacysta perseae*. Strand + –refers to major/minor strand. Negative values of IGS (intergenic spacer) refer to overlapping nucleotides.

<table><thead><tr><th>Locus</th><th>Strand</th><th>Location (bp)</th><th>Size (bp)</th><th>Anticodon</th><th>Start codon</th><th>Stop codon</th><th>IGS</th></tr></thead><tbody><tr><td>tRNA<sub>Ile</sub></td><td>+</td><td>1–63</td><td>63</td><td>GAT</td><td></td><td></td><td>0</td></tr><tr><td>tRNA<sub>Gln</sub></td><td>–</td><td>65–132</td><td>68</td><td>TTG</td><td></td><td></td><td>1</td></tr><tr><td>tRNA<sub>Met</sub></td><td>+</td><td>134–196</td><td>63</td><td>CAT</td><td></td><td></td><td>1</td></tr><tr><td>ND2</td><td>+</td><td>197–1178</td><td>982</td><td></td><td>ATA</td><td>T–</td><td>0</td></tr><tr><td>tRNA<sub>Trp</sub></td><td>+</td><td>1179–1245</td><td>67</td><td>TCA</td><td></td><td></td><td>0</td></tr><tr><td>tRNA<sub>Cys</sub></td><td>–</td><td>1238–1301</td><td>64</td><td>GCA</td><td></td><td></td><td>-8</td></tr><tr><td>tRNATyr</td><td>–</td><td>1304–1367</td><td>64</td><td>GTA</td><td></td><td></td><td>2</td></tr><tr><td>COX1</td><td>+</td><td>1368–2903</td><td>1536</td><td></td><td>ATG</td><td>TAA</td><td>0</td></tr><tr><td>tRNA<sub>Leu</sub> (UUR)</td><td>+</td><td>2906–2969</td><td>64</td><td>TAA</td><td></td><td></td><td>2</td></tr><tr><td>COX2</td><td>+</td><td>2970–3648</td><td>679</td><td></td><td>ATT</td><td>T–</td><td>0</td></tr><tr><td>tRNA<sub>lys</sub></td><td>+</td><td>3649–3720</td><td>72</td><td>CTT</td><td></td><td></td><td>0</td></tr><tr><td>tRNA<sub>Asp</sub></td><td>+</td><td>3721–3781</td><td>61</td><td>GTC</td><td></td><td></td><td>0</td></tr><tr><td>ATP8</td><td>+</td><td>3782–3940</td><td>159</td><td></td><td>ATA</td><td>TAA</td><td>0</td></tr><tr><td>ATP6</td><td>+</td><td>3934–4599</td><td>666</td><td></td><td>ATG</td><td>TAA</td><td>-7</td></tr><tr><td>COX3</td><td>+</td><td>4586–5372</td><td>787</td><td></td><td>ATT</td><td>T–</td><td>-14</td></tr><tr><td>tRNA<sub>Gly</sub></td><td>+</td><td>5373–5435</td><td>63</td><td>TCC</td><td></td><td></td><td>0</td></tr><tr><td>ND3</td><td>+</td><td>5436–5789</td><td>354</td><td></td><td>ATT</td><td>TAA</td><td>0</td></tr><tr><td>tRNA<sub>Ala</sub></td><td>+</td><td>5793–5853</td><td>61</td><td>TGC</td><td></td><td></td><td>3</td></tr><tr><td>tRNA<sub>Arg</sub></td><td>+</td><td>5854–5914</td><td>61</td><td>TCG</td><td></td><td></td><td>0</td></tr><tr><td>tRNA<sub>Asn</sub></td><td>+</td><td>5915–5979</td><td>65</td><td>GTT</td><td></td><td></td><td>0</td></tr><tr><td>tRNA<sub>Ser</sub> (AGN)</td><td>+</td><td>5979–6047</td><td>69</td><td>GCT</td><td></td><td></td><td>-1</td></tr><tr><td>tRNA<sub>Glu</sub></td><td>+</td><td>6048–6112</td><td>65</td><td>TTC</td><td></td><td></td><td>0</td></tr><tr><td>tRNA<sub>Phe</sub></td><td>R</td><td>6111–6174</td><td>64</td><td>GAA</td><td></td><td></td><td>-2</td></tr><tr><td>ND5</td><td>–</td><td>6175–7853</td><td>1679</td><td></td><td>ATT</td><td>TA–</td><td>0</td></tr><tr><td>tRNA<sub>His</sub></td><td>–</td><td>7854–7918</td><td>65</td><td>GTG</td><td></td><td></td><td>0</td></tr><tr><td>ND4</td><td>–</td><td>7921–9240</td><td>1320</td><td></td><td>ATG</td><td>TAA</td><td>2</td></tr><tr><td>ND4L</td><td>–</td><td>9234–9515</td><td>282</td><td></td><td>ATA</td><td>TAA</td><td>-7</td></tr><tr><td>tRNA<sub>Thr</sub></td><td>+</td><td>9517–9578</td><td>62</td><td>TGT</td><td></td><td></td><td>1</td></tr><tr><td>tRNA<sub>Pro</sub></td><td>–</td><td>9578–9643</td><td>66</td><td>TGG</td><td></td><td></td><td>-1</td></tr><tr><td>ND6</td><td>+</td><td>9645–10,124</td><td>480</td><td></td><td>ATT</td><td>TAA</td><td>1</td></tr><tr><td>CYTB</td><td>+</td><td>10,124–11,264</td><td>1141</td><td></td><td>ATT</td><td>T–</td><td>-1</td></tr><tr><td>tRNA<sub>Ser</sub> (TCN)</td><td>+</td><td>11,265–11,333</td><td>69</td><td>TGA</td><td></td><td></td><td>0</td></tr><tr><td>ND1</td><td>–</td><td>11,353–12,279</td><td>927</td><td></td><td>ATA</td><td>TAA</td><td>19</td></tr><tr><td>tRNA<sub>Leu</sub> (CUN)</td><td>–</td><td>12,281–12,344</td><td>64</td><td>TAG</td><td></td><td></td><td>0</td></tr><tr><td>IrRNA</td><td>–</td><td>12,345–13,564</td><td>1220</td><td></td><td></td><td></td><td>0</td></tr><tr><td>tRNA<sub>Val</sub></td><td>–</td><td>13,565–13,632</td><td>68</td><td>TAC</td><td></td><td></td><td>0</td></tr><tr><td>srRNA</td><td>–</td><td>13,633–14,424</td><td>792</td><td></td><td></td><td></td><td>0</td></tr><tr><td>Control region</td><td></td><td>14,425–15,850</td><td>1426</td><td></td><td></td><td></td><td>0</td></tr></tbody></table>

completion of termination codons [43,44,67]. Incomplete stop codons can also be found in *C. ciliata* and are shared with many arthropods [68].

### 3.3. Ribosomal and transfer RNA genes

tRNA genes locations and lengths are similar to those of other insect mitogenomes. IrRNA is located between tRNA<sub>Leu</sub>(CUN) and tRNA<sub>Val</sub> and is 1220 bp-long. srRNA is located between tRNA<sub>Val</sub> and the control region and is 792 bp-long. Their AT content is respectively 81.0% and 79.8%.

The classical set of 22 tRNAs found in arthropods is present in *P. perseae*. Their lengths vary between 61 bp (tRNA<sub>Ala</sub> and tRNA<sub>Arg</sub>) and 72 bp (tRNA<sub>lys</sub>). Secondary structures of tRNAs are schematized in Fig. 2. The classical clover leaf structure was observed for each of them, except for tRNA<sub>Ser</sub>(AGN), in which the D arm is reduced to a simple loop, as in many insects, and more generally, in most bilaterians [11,15]. The unusual TTC anticodon observed in *C. ciliata* for tRNA<sub>Ser</sub>(AGN) was not found in *P. perseae* which exhibits the classical GCT anticodon.

<caption>**Table 3**</caption>

Nucleotide composition of the mitogenomes of Tingidae: *Pseudacysta perseae* (P) and *Corythucha ciliata* (C).

<table><thead><tr><th rowspan="2"></th><th colspan="2">A + T%</th><th colspan="2">AT skew</th><th colspan="2">GC skew</th></tr><tr><th>P</th><th>C</th><th>P</th><th>C</th><th>P</th><th>C</th></tr></thead><tbody><tr><td>Whole genome</td><td>79.8</td><td>77.0</td><td>0.05</td><td>0.13</td><td>-0.23</td><td>-0.16</td></tr><tr><td>Protein coding genes</td><td>76.4</td><td>76.4</td><td>-0.12</td><td>0.16</td><td>-0.01</td><td>-0.17</td></tr><tr><td>1st codon position</td><td>72.7</td><td>74.2</td><td>0.18</td><td>0.36</td><td>-0.07</td><td>-0.02</td></tr><tr><td>2nd codon position</td><td>69.8</td><td>73.5</td><td>-0.03</td><td>0.46</td><td>-0.15</td><td>-0.22</td></tr><tr><td>3rd codon position</td><td>86.6</td><td>80.7</td><td>0.20</td><td>0.24</td><td>-0.54</td><td>-0.33</td></tr><tr><td>Protein coding genes-M<sup>a</sup></td><td>75.0</td><td>75.3</td><td>0.00</td><td>0.44</td><td>-0.16</td><td>-0.15</td></tr><tr><td>1st codon position</td><td>70.6</td><td>78.9</td><td>0.16</td><td>0.06</td><td>0.09</td><td>-0.30</td></tr><tr><td>2nd codon position</td><td>68.6</td><td>76.5</td><td>-0.35</td><td>0.16</td><td>-0.23</td><td>-0.07</td></tr><tr><td>3rd codon position</td><td>85.7</td><td>70.5</td><td>0.16</td><td>-0.10</td><td>-0.53</td><td>-0.11</td></tr><tr><td>Protein coding genes-m</td><td>78.6</td><td>78.0</td><td>0.31</td><td>0.35</td><td>-0.26</td><td>-0.21</td></tr><tr><td>1st codon position</td><td>76.1</td><td>77.4</td><td>0.21</td><td>0.33</td><td>-0.39</td><td>-0.19</td></tr><tr><td>2nd codon position</td><td>71.6</td><td>79.8</td><td>0.45</td><td>0.37</td><td>-0.03</td><td>-0.23</td></tr><tr><td>3rd codon position</td><td>88.0</td><td>77.2</td><td>0.27</td><td>0.33</td><td>-0.55</td><td>-0.23</td></tr><tr><td>tRNAs</td><td>80.0</td><td>78.3</td><td>0.04</td><td>0.08</td><td>-0.11</td><td>-0.11</td></tr><tr><td>tRNAs-M</td><td>79.9</td><td>78.2</td><td>0.05</td><td>0.66</td><td>0.03</td><td>0.01</td></tr><tr><td>tRNAs-m</td><td>80.1</td><td>78.5</td><td>0.03</td><td>0.11</td><td>-0.37</td><td>-0.31</td></tr><tr><td>rRNAs</td><td>80.5</td><td>79.2</td><td>0.05</td><td>0.11</td><td>-0.27</td><td>-0.25</td></tr><tr><td>Control region</td><td>82.2</td><td>74.0</td><td>-0.07</td><td>-0.07</td><td>-0.01</td><td>0.10</td></tr></tbody></table>

<footnote><sup>a</sup> M/m: genes encoded on the major/minor strand.</footnote>
<figure><img src="image_4.png" /><figcaption>Fig. 2. Secondary structure of *Pseudacysta perseae* mitochondrial tRNAs.</figcaption></figure>

### 3.4. Non-coding regions

Nine short intergenic spacers (IGS) and a long control region were identified, matching the usual organisation of insect mitogenomes. Most of IGS are very short, with less than four base pairs, and eight overlapping sequences are found (Table 2). The longest IGS are the 19 bp found between ND1 and tRNA<sub>Ser2</sub>. The control region is 1426 bp-long and is located between srRNA and tRNA<sub>Ile</sub> genes. In *P. persea*, it

exhibits one of the highest A + T content of the genome. The alternative name “A + T rich region” that has been used for insect control region could therefore be appropriate in this case [18]. However, this designation is hardly generalizable for heteropterans. Indeed, in *C. ciliata* and other bugs such as *T. dimidiata* [48] and *Valentia hoffmanni* [69], the control region exhibits a lower A + T content compared to that the whole mitogenome. In *P. persea* and *C. ciliata*, the control regions can be divided in a similar way (Fig. 3a):
<figure><img src="image_5.png" /><figcaption>Fig. 3. (Color online.) A. Structural organization of the mitochondrial control region of Tingidae. The boxes with roman numerals indicate tandem repeat units; B. Putative stem loop structures found in the mitochondrial control region of *Pseudacysta perseae*.</figcaption></figure>

* a region with a G + C content higher than the average;

* a region with a high A + T content;

* a region with a G + C content higher than the average, and;

* a region containing tandem repeats.

Tandem repeats in *P. perseae* control region consisted of six 38-bp repeated units and one 31-bp unit (corresponding to a partial copy of the previous ones) followed by two 180 bp repeated units and one 84-bp partial copy. In *C. ciliata*, it consisted of two 189-bp repeated units followed by a 48-bp partial copy. However, the sequencing technology used in this study does not allow inferring repeats copy numbers with certainty. We identified three putative stem loop structures in the control region of *P. perseae* (Fig. 3b). Such features may be involved in the mitochondrial DNA replication mechanism [70,71], but further studies are needed to assess this hypothesis.

In *Drosophila* and other insect species, a conserved string of about 20 Ts was identified in the control region and suspected to be involved in the initiation of replication [17]. Similarly, a conserved sequence block including a “G element” was observed in Reduviidae [49]. None of these were found in *P. perseae*. The alignment of the two lace bugs mitogenomes revealed various identical segments all along the control region, but more data would be necessary to know if some of them are well conserved among Tingidae.

### 3.5. Nucleotide content and codon usage

The nucleotide composition is strongly biased toward adenine and thymine in *P. perseae* mitogenome, with A + T representing 79.8% of the whole sequence and ranging from 76.4% in the protein-coding genes, 80.0% in tRNA genes, and 80.5% in rRNA genes to 82.2% in the control region. The nucleotide bias is reflected in the codon usage. AT-rich codons are predominant in protein-coding genes, with the most prevalent being in order TTA (Leu), ATT (Ile), ATA (Met) and TTT (Phe). The relative synonymous codon usage (RSCU) clearly indicates that AT rich codons are favoured among synonymous codons (Table 4). At the third codon position,

AT content is particularly high (86.6%), and G nucleotides are under-represented (GC skew = -0.54). AT content, as well as A-T and G-C skews patterns are similar between the two lace bugs mitochondrial genomes (Table 3).

<table><caption>Table 4
Codon usage of *Pseudacysta perseae* mitogenome protein-coding genes.</caption><thead><tr><th>Amino Acid<sup>a</sup></th><th>Codon</th><th>N<sup>b</sup></th><th>RSCU<sup>c</sup></th><th>Amino Acid</th><th>Codon</th><th>N</th><th>RSCU</th></tr></thead><tbody><tr><td rowspan="2">F</td><td>UUU</td><td>291</td><td>1.69</td><td rowspan="2">Y</td><td>UAU</td><td>150</td><td>1.72</td></tr><tr><td>UUC</td><td>54</td><td>0.31</td><td>UAC</td><td>24</td><td>0.28</td></tr><tr><td rowspan="6">L</td><td>UUA</td><td>363</td><td>4.02</td><td rowspan="2">H</td><td>CAU</td><td>51</td><td>1.59</td></tr><tr><td>UUG</td><td>61</td><td>0.68</td><td>CAC</td><td>13</td><td>0.41</td></tr><tr><td>CUU</td><td>52</td><td>0.58</td><td rowspan="2">Q</td><td>CAA</td><td>45</td><td>1.7</td></tr><tr><td>CUC</td><td>8</td><td>0.09</td><td>CAG</td><td>8</td><td>0.3</td></tr><tr><td>CUA</td><td>54</td><td>0.6</td><td rowspan="2">N</td><td>AAU</td><td>159</td><td>1.58</td></tr><tr><td>CUG</td><td>4</td><td>0.04</td><td>AAC</td><td>42</td><td>0.42</td></tr><tr><td rowspan="2">I</td><td>AUU</td><td>313</td><td>1.67</td><td rowspan="2">K</td><td>AAA</td><td>108</td><td>1.74</td></tr><tr><td>AUC</td><td>62</td><td>0.33</td><td>AAG</td><td>16</td><td>0.26</td></tr><tr><td rowspan="2">M</td><td>AUA</td><td>304</td><td>1.77</td><td rowspan="2">D</td><td>GAU</td><td>51</td><td>1.55</td></tr><tr><td>AUG</td><td>39</td><td>0.23</td><td>GAC</td><td>15</td><td>0.45</td></tr><tr><td rowspan="4">V</td><td>GUU</td><td>81</td><td>1.89</td><td rowspan="2">E</td><td>GAA</td><td>77</td><td>1.75</td></tr><tr><td>GUC</td><td>10</td><td>0.23</td><td>GAG</td><td>11</td><td>0.25</td></tr><tr><td>GUA</td><td>70</td><td>1.64</td><td rowspan="2">C</td><td>UGU</td><td>43</td><td>1.91</td></tr><tr><td>GUG</td><td>10</td><td>0.23</td><td>UGC</td><td>2</td><td>0.09</td></tr><tr><td rowspan="5">S</td><td>UCU</td><td>115</td><td>2.67</td><td rowspan="2">W</td><td>UGA</td><td>80</td><td>1.72</td></tr><tr><td>UCC</td><td>14</td><td>0.33</td><td>UGG</td><td>13</td><td>0.28</td></tr><tr><td>UCA</td><td>85</td><td>1.98</td><td rowspan="2">R</td><td>CGU</td><td>13</td><td>1.04</td></tr><tr><td>UCG</td><td>10</td><td>0.23</td><td>CGC</td><td>1</td><td>0.08</td></tr><tr><td rowspan="4">P</td><td>CCU</td><td>69</td><td>2.28</td><td rowspan="2">S</td><td>CGA</td><td>33</td><td>2.64</td></tr><tr><td>CCC</td><td>10</td><td>0.33</td><td>CGG</td><td>3</td><td>0.24</td></tr><tr><td>CCA</td><td>39</td><td>1.29</td><td rowspan="2">G</td><td>AGU</td><td>36</td><td>0.84</td></tr><tr><td>CCG</td><td>3</td><td>0.1</td><td>AGC</td><td>1</td><td>0.02</td></tr><tr><td rowspan="5">T</td><td>ACU</td><td>64</td><td>1.39</td><td rowspan="5">G</td><td>AGA</td><td>83</td><td>1.93</td></tr><tr><td>ACC</td><td>15</td><td>0.33</td><td>AGG</td><td>0</td><td>0</td></tr><tr><td>ACA</td><td>103</td><td>2.24</td><td>GGU</td><td>71</td><td>1.65</td></tr><tr><td>ACG</td><td>2</td><td>0.04</td><td>GGC</td><td>4</td><td>0.09</td></tr><tr><td>GCU</td><td>47</td><td>1.9</td><td>GGA</td><td>74</td><td>1.72</td></tr><tr><td rowspan="3">A</td><td>GCC</td><td>10</td><td>0.4</td><td rowspan="3"></td><td>GGG</td><td>23</td><td>0.53</td></tr><tr><td>GCA</td><td>39</td><td>1.58</td><td></td><td></td><td></td></tr><tr><td>GCG</td><td>3</td><td>0.12</td><td></td><td></td><td></td></tr></tbody></table>

<footnote><sup>a</sup> Amino acids are labeled according to the IUPAC-IUB single-letter codes.</footnote>

<footnote><sup>b</sup> N: total number in all PCGs.</footnote>

<footnote><sup>c</sup> RSCU: relative synonymous codon usage.</footnote>
3.6. Phylogenetic analysis

Bayesian inference and Maximum Likelihood analysis (ML) generated phylogenetic trees with very similar topologies (Figs. 4 and 5). *P. perseae* is placed as a sister group of *C. ciliata*. Miroidea is paraphyletic in our analysis, in contrast with previous results based on morphological and molecular data [20,72], but in accordance with the results of Tian et al. [73] inferred from mitochondrial and nuclear genes. The 18 other superfamilies represented in our dataset are monophyletic. At the infra-order level, Pentatomorpha is monophyletic with the following relationships between the five superfamilies: Aradoidea + (Pentatomoidea + [Lygaeoidea + [Pyrrhocoroidea + Coreoidea]]). These relationships are strongly supported in our analyses and do not confirm the results of previous studies based on a subset of the present mitogenomic data that placed Coreoidea and Lygaeoidea as sister groups [24,74]. On the tree generated by Bayesian inference, Nepomorpha is monophyletic with the following relationships between the six superfamilies: (Pleoidea + Corixoidea) + ([Notonectoidea + Naucoroidea] + [Nepoidea + Ochteroidea]). ML analysis differs by positioning Pleoidea as a sister group of the remaining Nepomorpha. These results are not consistent with those of a study based on molecular

and morphological data [75]. Interestingly, they also contradict a previous analysis of a subset of the present mitogenomic dataset by confirming the monophyly of Nepomorpha including Pleoidea for which an infraordinal status was proposed [69]. Gerromorpha and Leptopodomorpha are also monophyletic but only two species of each were included in the study. Enicocephalomorpha was only represented by *Stenopirates* sp.

Our analysis supports the paraphyly of Cimicomorpha that consisted of four different clades: (Cimicoidea + Na-boidea), Reduvioidea, Miridae and Tingidae. Tingidae is placed as a sister group to all remaining Heteroptera. However, the monophyly of Cimicomorpha has been widely accepted and is supported by various analyses. Schuh et al. [20] conducted a total-evidence phylogenetic analysis of Cimicomorpha based on 73 morphological characters, nuclear and mitochondrial DNA for 92 taxa including eight outgroups belonging to Pentatomomorpha, Leptopodomorpha and Nepomorpha. Their result con-firmed the monophyly of the group except for the inclusion of Thaumastocoridae, which is not represented in our dataset. Similarly, Tian et al. [73] inferred the monophyly of Cimicomorpha from the analysis of nuclear and mitochondrial data on 46 taxa including three outgroups belonging to Pentatomomorpha and Leptopodomorpha.

<figure><img src="image_6.png" /><figcaption>Fig. 4. (Color online.) Phylogenetic tree inferred by Bayesian analysis from 55 heteropteran mitogenomes. Bayesian posterior probabilities expressed in percentages are indicated on the nodes.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>Fig. 5. (Color online.) Phylogenetic tree inferred by ML analysis from 55 heteropteran mitogenomes. Bootstrap support values expressed in percentages are indicated on the nodes.</figcaption></figure>

More recently, Li et al. [76] performed a phylogenetic analysis of Heteroptera with a dataset of nuclear and mitochondrial DNA for 83 species including 22 Cimicomorpha, which also confirmed the monophyly of the group.

Infraordinal relationships are conserved in both ML and Bayesian analyses: Tingidae + (Miridae + (Pentatomomorpha + ((Cimicoidea + Naboidea) + (Reduviodea + (Leptopodomorpha + ((Enicocephalomorpha + Gerromorpha) + Nepomorpha)))))). These results are in contrast with the general consensus that considers Enicocephalomorpha as the early infraorder of Heteroptera [19]. However, the relationships between the infraorders of Heteroptera remain controversial. Only few phylogenetic studies have addressed the question and those have only included a small number of taxa [77–79]. Mitogenomic data provide a new insight in this regard, but more taxa should be added to the current database, especially in the poorly represented infra-orders Enicocephalomorpha, Gerromorpha, Leptopodomorpha and Dipsocoromorpha.

Our results are incongruent with generally accepted hypothesis [19]. On the other hand, they are in accordance with previous studies based on a subset of the present mitogenomic data [24,50]. Analyses performed on individual genes by Tian et al. [73] and Schuh et al. [20] indicate that the monophyly of Cimicomorpha is only supported by

nuclear rDNA partitions, among a dataset comprising 16 rDNA, 18S rDNA, 28S rDNA and COI. Incongruences of phylogenetic analyses among different genomic regions are a well-known issue that can have various biological causes such as incomplete lineage sorting or evolution rates variations among sites [80]. For instance, Kopp and True [81] reported that single-gene datasets of *D. melanogaster* species group produce several strongly supported conflicting clades. Lin and Danforth [82] studied the differences in the pattern of nucleotide substitution among nuclear and mitochondrial genes and concluded that insect phylogenetic studies should increasingly focus on nuclear data. This raises the limitations of phylogenetic inferences from complete mitochondrial genomes only. While reducing stochastic errors by providing large molecular datasets, this approach is susceptible to potential site-specific bias and would benefit from a conjoint analysis with other genes. Unfortunately, there is a lack of correspondence between publicly available mitogenomic and nuclear data. Currently, on Genbank, nuclear ribosomal sequences are available for only 14 species out of 55 represented in our mitochondrial dataset. Moreover, these are mainly partial sequences that do not necessarily overlap. The genome skimming-approach employed in the present study could provide an interesting improvement in this regard. In a single experimentation, it allows to recover the full mitogenome sequence as well as
nuclear genes that are classically used for phylogenetic
inference.

Acknowledgements

We would like to thank Maria Kamilari for help with
DNA extraction and Eric Coissac for help with genome
assembly. This work was supported by « Investissement
d'avenir » grants managed by the Agence Nationale de la
Recherche (CEBA, ref. ANR-10-LABX-25-01; TULIP, ANR-
10-LABX-41, ANR-11-IDEX-0002-02) as well as project
METABAR (ANR-11-BSV7-0020). We are grateful to the
Genotoul bioinformatics platform Toulouse Midi-Pyrénées
for providing computing and storage resources. The
material studied was collected in the frame of the
CAFOTROP-Muséum project.

References

[1] J.C. Avise, J. Arnold, R.M. Ball, E. Bermingham, T. Lamb, J.E. Neigel, et al., Intraspecific phylogeography: the mitochondrial DNA bridge between population genetics and systematics, Annu. Rev. Ecol. Syst. 18 (1987) 489–522.

[2] C. Moritz, T.E. Dowling, W.M. Brown, Evolution of animal mitochondrial DNA: relevance for population biology and systematics, Annu. Rev. Ecol. Syst. (1987) 269–292.

[3] S.C. Straub, M. Parks, K. Weitemier, M. Fishbein, R.C. Cronn, A. Liston, Navigating the tip of the genomic iceberg: next-generation sequencing for plant systematics, Am. J. Bot. 99 (2012) 349–364.

[4] M. Berman, C.M. Austin, A.D. Miller, Characterisation of the complete mitochondrial genome and 13 microsatellite loci through next-generation sequencing for the New Caledonian spider-ant *Leptomyrmex pallens*, Mol. Biol. Rep. 41 (2014) 1179–1187.

[5] S.R. Doyle, I.S. Griffith, N.P. Murphy, J.M. Strugnell, Low-coverage MiSeq next generation sequencing reveals the mitochondrial genome of the Eastern Rock Lobster, *Sagmariusus verreauxi*, Mitochondrial DNA (2014) 1–2.

[6] K.F. Thompson, S. Patel, L. Williams, P. Tsai, R. Constantine, C.S. Baker, et al., High coverage of the complete mitochondrial genome of the rare Gray's beaked whale (*Mesoplodon grayi*) using Illumina next generation sequencing, Mitochondrial DNA (2014) 1–2.

[7] A.J. Veale, L. Williams, P. Tsai, V. Thakur, S. Lavery, The complete mitochondrial genomes of two chiton species (*Sypharochiton pellisperpentis* and *Sypharochiton sinclairi*) obtained using Illumina next generation sequencing, Mitochondrial DNA (2014) 1–2.

[8] C. Gissi, F. Iannelli, G. Pesole, Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species, Heredity 101 (2008) 301–320.

[9] P. Salvato, M. Simonato, A. Battisti, E. Negrisolo, The complete mitochondrial genome of the bag-shelter moth *Ochrogaster lunifer* (Lepidoptera, Notodontidae), BMC Genomics 9 (2008) 331.

[10] J.L. Boore, Animal mitochondrial genomes, Nucleic Acids Res. 27 (1999) 1767–1780.

[11] D.R. Wolstenholme, Animal mitochondrial DNA: structure and evolution, Int. Rev. Cytol. 141 (1992) 173–216.

[12] A.C.M. Junqueira, A.C. Lessinger, T.T. Torres, F.R. da Silva, A.L. Vettore, P. Arruda, et al., The mitochondrial genome of the blowfly *Chrysomya chloropyga* (Diptera: Calliphoridae), Gene 339 (2004) 7–15.

[13] R. Shao, S.C. Barker, The highly rearranged mitochondrial genome of the plague thrips, *Thrips imaginis* (Insecta: Thysanoptera): convergence of two novel gene boundaries and an extraordinary arrangement of rRNA genes, Mol. Biol. Evol. 20 (2003) 362–370.

[14] M.L. Thao, L. Baumann, P. Baumann, Organization of the mitochondrial genomes of whiteflies, aphids, and psyllids (Hemiptera, Sternorrhyncha), BMC Evol. Biol. 4 (2004) 25.

[15] M. Bernt, A. Braband, B. Schierwater, P.F. Stadler, Genetic aspects of mitochondrial genome evolution, Mol. Phylogenet. Evol. 69 (2013) 328–338.

[16] D.A. Clayton, Transcription and replication of animal mitochondrial DNAs, Int. Rev. Cytol. 141 (1992) 217–232.

[17] S. Saito, Replication origin of mitochondrial DNA in insects, Genetics 171 (2005) 1695–1705.

[18] D.-X. Zhang, G.M. Hewitt, Insect mitochondrial control region: a review of its structure, evolution and usefulness in evolutionary studies, Biochem. Syst. Ecol. 25 (1997) 99–120.

[19] C. Weirauch, R.T. Schuh, Systematics and evolution of heteroptera: 25 years of progress, Annu. Rev. Entomol. 56 (2011) 487–510.

[20] R.T. Schuh, C. Weirauch, W.C. Wheeler, Phylogenetic relationships within the Cimicorompha (Hemiptera: Heteroptera): a total-evidence analysis, Syst. Entomol. 34 (2009) 15–48.

[21] R.T. Schuh, G. Cassis, E. Guilbert, Description of the first recent macroporous species of Vianaidinae (Heteroptera: Tingidae) with comments on the phylogenetic relationships of the family within the Cimicorompha, J. N. Y. Entomol. Soc. 114 (2006) 38–53.

[22] L.T. Miller, Lace Bugs (Hemiptera: Tingidae), in: J.L. Capinera (Ed.), Encyclopedia of Entomology, Springer, Netherlands, 2008, pp. 2099–2102.

[23] E. Guilbert, J. Damgaard, C.A. D'haese, Phylogeny of the lacebugs (Insecta: Heteroptera: Tingidae) using morphological and molecular data, Syst. Entomol. 39 (2014) 431–441.

[24] W. Yang, W. Yu, Y. Du, The complete mitochondrial genome of the sycamore lace bug *Corythucha ciliata* (Hemiptera: Tingidae), Gene 532 (2013) 27–40.

[25] F. Mead, J.E. Pena, Avocado lace bug, *Pseudacysta perseae* (Hemiptera: Tingidae), Florida Dept. Agr. Consum. Serv. Div. Plant. Ind. Entomol. Circ. 346 (1991).

[26] J.E. Peña, Pests of avocado in Florida, in: Proc. Vth World Avocado Congress, 2003, pp. 487–494.

[27] F.J. Byrne, E.C. Humeres, A.A. Urena, M.S. Hoddle, J.G. Morse, Field evaluation of systemic imidacloprid for the management of avocado thrips and avocado lace bug in California avocado groves, Pest Manage. Sci. 66 (2010) 1129–1136.

[28] R.J. Gagné, J.E. Peña, F.E. Acevedo, A new lestodipolsine (Diptera: Cecidomyiidae) preying on the avocado lace bug, *Pseudacysta perseae* (Heteroptera: Tingidae) in Southern Florida, Fla. Entomol. 91 (2008) 43–48.

[29] C.M. Holguín, J.E. Peña, T.J. Henry, F. Acevedo, Biology of *Stethoconus praefectus* (Distant) (Heteroptera: Miridae), a newly established predator of the avocado lace bug, *Pseudacysta perseae* (Heteroptera: Tingidae), in Florida, Fla. Entomol. 92 (2009) 54–57.

[30] E.C. Humeres, J.G. Morse, R. Stouthamer, W. Roltsch, M.S. Hoddle, Evaluation of natural enemies and insecticides for control of *Pseudacysta perseae* (Hemiptera: Tingidae) on avocados in Southern California, Fla. Entomol. 92 (2009) 35–42.

[31] J.E. Peña, S.V. Triapitsyn, D. Long, G.A. Evans, W. Roltsch, First record of Erythmelus klopomor (Hymenoptera: Mymaridae) as a parasitoid of the avocado lace bug, *Pseudacysta perseae* Heteroptera: Tingidae), Fla. Entomol. 92 (2009) 394–395.

[32] M.F.S. Cabrera, M. Cermeli, Presencia de *Pseudacysta perseae* (Heidemann, 1908) (Insecta: Hemiptera: Tingidae) en Venezuela, Entomotropica 20 (2008) 271–273.

[33] J.C. Streito, Y. Morival, Première capture en Guyane Française de *Pseudacysta perseae* (Heidemann), 1908, un ravageur de l'avocatier (Heteroptera: Tingidae), Nouv. Rev. Entomol. 22 (2005) 191–192.

[34] J. Mariette, F. Escudié, N. Allias, G. Salin, C. Noirot, S. Thomas, et al., NG6: Integrated next generation sequencing storage and processing environment, BMC Genomics 13 (2012) 462.

[35] G. Besnard, F. Jühling, É. Chapuis, L. Zedane, É. Lhuillier, T. Mateille, et al., Fast assembly of the mitochondrial genome of a plant parasitic nematode (*Meloidogyne graminicola*) using next generation sequencing, C.R. Biologies 337 (2014) 295–301.

[36] A. Kocher, M. Kamilari, E. Lhuillier, E. Coissac, J. Péneau, J. Chave, et al., Shotgun assembly of the assassin bug *Brontostoma colossus* mitochondrial genome (Heteroptera, Reduviidae), Gene 552 (2014) 184–194.

[37] S. Cally, E. Lhuillier, A. Iribar, I. Garzón-Orduña, E. Coissac, J. Murienne, Shotgun assembly of the complete mitochondrial genome of the neotropical cracker butterfly *Hamadryas* epine, Mitochondrial DNA (2014) 1–3.

[38] C. Hahn, L. Bachmann, B. Chevreux, Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads—a baiting and iterative mapping approach, Nucleic Acids Res. 41 (2013) e129.

[39] V.H. Nguyen, D. Lavenier, PLAST: parallel local alignment search tool for database comparison, BMC Bioinform. 10 (2009) 329.

[40] D.R. Zerbino, E. Birney, Velvet: algorithms for de novo short read assembly using de Bruijn graphs, Genome Res. 18 (2008) 821–829.

[41] S. Yu, Y. Wang, D. Rédei, Q. Xie, W. Bu, Secondary structure models of 18S and 28S rRNAs of the true bugs based on complete rDNA sequences of Eurydema maracandica Oshanin, 1871 (Heteroptera, Pentatomidae), ZooKeys 319 (2013) 363–377.

[42] M. Bernt, A. Donath, F. Jühling, F. Externbrink, C. Florentz, G. Fritzsch, et al., MITOS: improved de novo metazoan mitochondrial genome annotation, Mol. Phylogenet. Evol. 69 (2013) 313–319.
[43] D. Ojala, J. Montoya, G. Attardi, tRNA punctuation model of RNA processing in human mitochondria, Nature 290 (1981) 470–474.

[44] J.B. Stewart, A.T. Beckenbach, Characterization of mature mitochondrial transcripts in Drosophila, and the implications for the tRNA punctuation model in arthropods, Gene 445 (2009) 49–57.

[45] K. Tamura, G. Stecher, D. Peterson, A. Filipski, S. Kumar, MEGA6: molecular evolutionary genetics analysis version 6.0, Mol. Biol. Evol. 30 (2013) 2725–2729.

[46] N.T. Perna, T.D. Kocher, Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes, J. Mol. Evol. 41 (1995) 353–358.

[47] G. Benson, Tandem repeats finder: a program to analyze DNA sequences, Nucleic Acids Res. 27 (1999) 573.

[48] E.M. Dotson, C.B. Beard, Sequence and organization of the mitochondrial genome of the Chagas disease vector, *Triatoma dimidiata*, Insect Mol. Biol. 10 (2001) 205–215.

[49] J. Gao, H. Li, X.L. Truong, X. Dai, J. Chang, W. Cai, Complete nucleotide sequence and organization of the mitochondrial genome of Sirtthenea flavipes (Hemiptera: Reduviidae: Peiratinae) and comparison with other assassin bugs, Zootaxa 3669 (2013) 1–16.

[50] H. Li, J. Gao, H. Liu, H. Liu, A. Liang, X. Zhou, et al., The architecture and complete sequence of mitochondrial genome of an assassin bug *Agriosphodrus dohri* (Hemiptera: Reduviidae), Int. J. Biol. Sci. 7 (2011) 792.

[51] F. Abascal, R. Zardoya, M.J. Telford, TranslatorX: multiple alignment of nucleotide sequences guided by amino acid translations, Nucleic Acids Res. 38 (2010) W7–W13.

[52] P. Kück, K. Meusemann, FASconCAT: convenient handling of data matrices, Mol. Phylogenet. Evol. 56 (2010) 1115–1118.

[53] R. Lanfear, B. Calcott, S.Y.W. Ho, S. Guindon, Partitionfinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses, Mol. Biol. Evol. 29 (2012) 1695–1701.

[54] A. Stamatakis, RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies, Bioinformatics (2014) [btu033].

[55] A. Stamatakis, P. Hoover, J. Rougemont, A rapid bootstrap algorithm for the RAxML web servers, Syst. Biol. 57 (2008) 758–771.

[56] F. Ronquist, J.P. Huelsenbeck, MrBayes 3: bayesian phylogenetic inference under mixed datasets, Bioinformatics 19 (2003) 1572–1574.

[57] D.O. Clary, D.R. Wolstenholme, The mitochondrial DNA molecule of *Drosophila yakuba*: nucleotide sequence, gene organization, and genetic code, J. Mol. Evol. 22 (1995) 252–271.

[58] M.H. de Bruijn, *Drosophila melanogaster* mitochondrial DNA, a novel organization and genetic code, Nature 304 (1983) 234–241.

[59] J.L. Boore, T.M. Collins, D. Stanton, L.L. Daehler, W.M. Brown, Deducing the pattern of arthropod phylogeny from mitochondrial-DNA rearrangements, Nature 376 (1995) 163–165.

[60] D.V. Lavrov, J.L. Boore, W.M. Brown, The complete mitochondrial DNA sequence of the horseshoe crab *Limulus polyphemus*, Mol. Biol. Evol. 17 (2000) 813–824.

[61] J.L. Boore, D.V. Lavrov, W.M. Brown, Gene translocation links insects and crustaceans, Nature 392 (1998) 667–668. http://dx.doi.org/10.1038/33577.

[62] Z. Zhou, Y. Huang, F. Shi, The mitochondrial genome of *Ruspolia dubia* (Orthoptera: Conocephalidae) contains a short A + T-rich region of 70 bp in length, Genome 50 (2007) 855–866.

[63] R. Garesse, *Drosophila melanogaster* mitochondrial DNA: gene organization and evolutionary considerations, Genetics 118 (1988) 649–663.

[64] E.-S. Lee, K.S. Shin, M.-S. Kim, H. Park, S. Cho, C.-B. Kim, The mitochondrial genome of the smaller tea tortrix *Adoxophyes honmai* (Lepidoptera: Tortricidae), Gene 373 (2006) 52–57.

[65] K. Yukuhiro, H. Sezutsu, M. Itoh, K. Shimizu, Y. Banno, Significant levels of sequence divergence and gene rearrangements have occurred between the mitochondrial genomes of the wild mulberry silkmoth, *Bombyx mandarina*, and its Close Relative, the Domesticated Silkmoth, *Bombyx mori*, Mol. Biol. Evol. 19 (2002) 1385–1389.

[66] J.W.O. Ballard, Comparative genomics of mitochondrial DNA in members of the *Drosophila melanogaster* subgroup, J. Mol. Evol. 51 (2000) 48–63.

[67] T. Nagaike, Human mitochondrial mrnas are stabilized with polyadenylation regulated by mitochondria-specific poly(a) polymerase and polynucleotide phosphorylase, J. Biol. Chem. 280 (2005) 19721–19727.

[68] J.L. Boore, The duplication/random loss model for gene rearrangement exemplified by mitochondrial genomes of deuterostome animals, in: D. Sankoff, J.H. Nadeau (Eds.), Comparative Genomics, Springer, Netherlands, 2000, pp. 133–147.

[69] J. Hua, M. Li, P. Dong, Y. Cui, Q. Xie, W. Bu, Phylogenetic analysis of the true water bugs (Insecta: Hemiptera: Heteroptera: Nepomorpha): evidence from mitochondrial genomes, BMC Evol. Biol. 9 (2009) 134.

[70] N. Song, A.-P. Liang, Complete mitochondrial genome of the small brown planthopper, *Laodelphax striatellus* (Delphacidae: Hemiptera), with a novel gene order, Zoolog. Sci. 26 (2009) 851–860.

[71] D.-X. Zhang, J.M. Szymura, G.M. Hewitt, Evolution and structural conservation of the control region of insect mitochondrial DNA, J. Mol. Evol. 40 (1995) 382–391.

[72] R.T. Schuh, P. Stys, Phylogenetic analysis of cimicoraphan family relationships (Heteroptera), J. N. Y. Entomol. Soc. (1991) 298–350.

[73] Y. Tian, W. Zhu, M. Li, Q. Xie, W. Bu, Influence of data conflict and molecular phylogeny of major clades in Cimicoraph true bugs (Insecta: Hemiptera: Heteroptera), Mol. Phylogenet. Evol. 47 (2008) 581–597.

[74] J. Hua, M. Li, P. Dong, Y. Cui, Q. Xie, W. Bu, Comparative and phylogenomic studies on the mitochondrial genomes of Pentatomomorpha (Insecta: Hemiptera: Heteroptera), BMC Genomics 9 (2008) 610.

[75] M.B. Hebsgaard, N.M. Andersen, J. Damgaard, Phylogeny of the true water bugs (Nepomorpha: Hemiptera-Heteroptera) based on 16S and 28S rDNA and morphology, Syst. Entomol. 29 (2004) 488–508.

[76] M. Li, Y. Tian, Y. Zhao, W. Bu, Higher level phylogeny and the first divergence time estimation of Heteroptera (Insecta: Hemiptera) based on multiple genes, PLoS ONE 7 (2012) e32152.

[77] M. Mahner, Systema Cryptoceratumus Phylogeneticum (Insecta, Heteroptera), Schweizerbart, Stuttgart, Germany, 1993.

[78] W.C. Wheeler, R.T. Schuh, R. Bang, Cladistic relationships among higher groups of Heteroptera: congruence between morphological and molecular data sets, Insect Syst. Evol. 24 (1993) 121–137.

[79] Q. Xie, Y. Tian, L. Zheng, W. Bu, 18S rRNA hyper-elongation and the phylogeny of Euhemiptera (Insecta: Hemiptera), Mol. Phylogenet. Evol. 47 (2008) 463–471.

[80] A. Som, Causes, consequences and solutions of phylogenetic incongruence, Brief. Bioinform. (2014).

[81] A. Kopp, J.R. True, Phylogeny of the Oriental *Drosophila melanogaster* species group: a multilocus reconstruction, Syst. Biol. 51 (2002) 786–805.

[82] C.-P. Lin, B.N. Danforth, How do insect nuclear and mitochondrial gene substitution patterns differ? Insights from Bayesian analyses of combined datasets, Mol. Phylogenet. Evol. 30 (2004) 686–702.

[83] N. Song, A.-P. Liang, C.-P. Bu, A molecular phylogeny of Hemiptera inferred from mitochondrial genome sequences, PloS One 7 (2012) e48778.

[84] Y. Wang, H. Li, H. Xun, W. Cai, Complete mitochondrial genome sequence of the plant bug *Adelphocoris fasciaticollis* (Hemiptera: Heteroptera: Miridae), Mitochondrial DNA (2014) 1–2.

[85] P. Wang, H. Li, Y. Wang, J.-H. Zhang, X. Dai, J. Chang, et al., The mitochondrial genome of the plant bug *Apolygus lucorum* (Hemiptera: Miridae): presently known as the smallest in Heteroptera, Insect Sci. 21 (2013) 159–173.

[86] X. Dai, H. Xun, J. Chang, J. Zhang, B. Hu, H. Li, et al., The complete mitochondrial genome of the plant bug *Nesidiocoris tenuis* (Reuter) (Hemiptera: Miridae: Bryocorinae: Dicyphini), Zootaxa (2012) 30–44.

[87] H. Li, H. Liu, L. Cao, A. Shi, H. Yang, W. Cai, The complete mitochondrial genome of the damsel bug *Alloehrhynchus bakeri* (Hemiptera: Nabidae), Int. J. Biol. Sci. 8 (2012) 93.

[88] H. Li, J. Gao, W. Cai, Complete mitochondrial genome of the assassin bug *Oncocaphalus breviscutum* (Hemiptera: Reduviidae), Mitochondrial DNA (2013) 1–2.

[89] H. Li, H. Liu, A. Shi, P. Štyš, X. Zhou, W. Cai, The complete mitochondrial genome and novel gene arrangement of the unique-headed bug *Steponipares* sp. (Hemiptera: Enicocephalidae), PloS One 7 (2012) e29419.

[90] A.M. Shi, H. Li, X.S. Bai, X. Dai, J. Chang, E. Guilbert, et al., The complete mitochondrial genome of the flat bug *Aradacanthina heissi* (Hemiptera: Aradidae), Zootaxa 3238 (2012) 23–38.

[91] H. Li, A. Shi, F. Song, W. Cai, Complete mitochondrial genome of the flat bug *Brachyrhynchus hsiaoi* (Hemiptera: Aradidae), Mitochondrial DNA (2014) 1–2.

[92] L. Liu, H. Li, F. Song, W. Song, X. Dai, J. Chang, et al., The mitochondrial genome of *Coridius chinensis* (Hemiptera: Dinidoridae), Zootaxa 3537 (2012) 29–40.

[93] Q.-L. Zhang, M.-L. Yuan, Y.-Y. Shen, The complete mitochondrial genome of *Dolycoris baccarum* (Insecta: Hemiptera: Pentatomidae), Mitochondrial DNA 24 (2013) 469–471.

[94] W. Lee, J. Kang, C. Jung, K. Hoelmer, S.H. Lee, S. Lee, Complete mitochondrial genome of brown marmorated stink bug *Halyomorpha halys* (Hemiptera: Pentatomidae), and phylogenetic relationships of hemipteran suborders, Mol. Cells 28 (2009) 155–165.

[95] W. Song, H. Li, F. Song, L. Liu, P. Wang, H. Xun, et al., The complete mitochondrial genome of a tessaratomid bug, *Eusthenes cupreus* (Hemiptera: Heteroptera: Pentatomomorpha: Tessaratomidae), Zootaxa 3620 (2013) 260–272.

[96] D. Yuting, P. Li, F. Jiang, Z. Song, X. Ye, Yuang, et al., Sequence and organization of the mitochondrial genome of an urostylid bug, *Urochela quadrinotata* Reuter (Hemiptera: Urostylididae), Entomotaxonomia 34 (2012) 613–623.