Phylogenetic analysis of the true water bugs (Insecta: Hemiptera: Heteroptera: Nepomorpha): evidence from mitochondrial genomes

Complete mitochondrial genome of Saldoida armata Horváth, 1911 (Heteroptera: Saldidae) and phylogenetic analysis

The complete mitochondrial genome of Saldoida armata (Heteroptera: Saldidae) is 16,049 bp in length, comprising 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and a control region. All the PCGs are initially encoded by ATN, TTG or GTG, and terminated coding with TAA or a single T. With the exception of trnS(AGN), all tRNAs exhibit a typical cloverleaf secondary structure. Phylogenetic analysis reveals the sister relationship of S. armata with other Saldidae members. The complete mitogenome of S. armata will provide useful genetic information for species identification, phylogenetic analysis and conservation of this species.

The complete mitochondrial genome of Sigara lateralis (Leach, 1817) (Nepomorpha: Corixidae)

The Sigara lateralis (Leach, 1817) is a small aquatic insect belonging to the family Corixidae. The study aims to reveal the genomic data of the mitochondrial genome of S. lateralis. The length of its circular mitochondrial genome is 15,725 bp long with an A + T bias (75.96%). The mitogenome comprises 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes, and two rRNA genes. The phylogenetic analyses showed that the S. lateralis is the closest to S. septemlineata. These findings will help the conservation of Corixidae from the perspective of genetic evolution.

Analysis of Thaumatotibia leucotreta (Lepidoptera: Tortricidae: Olethreutinae) mitochondrial genomes in the context of a recent host range expansion

BACKGROUND

The false codling moth (FCM), Thaumatotibia leucotreta (Meyrick, 1913), is a significant pest of various important economic crops and is a EU quarantine pest. In the last decade the pest has been reported on Rosa spp. In this study we determined whether this shift occurred within specific FCM populations across seven eastern sub-Saharan countries or whether the species opportunistically switches to this novel host as it presents itself. To achieve this, we assessed the genetic diversity of complete mitogenomes of T. leucotreta specimens intercepted at import and analysed potential linkages with the geographical origin and host species.

RESULTS

Genomic, geographical and host information were integrated into a T. leucotreta Nextstrain build which contains 95 complete mitogenomes generated from material intercepted at import between January 2013 and December 2018. Samples represented seven sub-Saharan countries and mitogenomic sequences grouped in six main clades.

DISCUSSION

If host strains of FCM would exist, specialization from a single haplotype towards the novel host is expected. Instead, we find specimens intercepted on Rosa spp. in all six clades. The absence of linkage between genotype and host suggests opportunistic expansion to the new host plant. This underlines risks of introducing new plant species to an area as the effect of pests already present on the new plant might be unpredictable with current knowledge.

Novel Gene Rearrangement Pattern in Pachycrepoideus vindemmiae Mitochondrial Genome: New Gene Order in Pteromalidae (Hymenoptera: Chalcidoidea)

The mitochondrial genomes of Muscidifurax similadanacus, M. sinesensilla, Nasonia vitripennis, and Pachycrepoideus vindemmiae were sequenced to better understand the structural evolution of Pteromalidae mitogenomes. These newly sequenced mitogenomes all contained 37 genes. Nucleotide composition was AT-biased and the majority of the protein-coding genes exhibited a negative AT skew. All 13 protein-coding genes (PCGs) initiated with the standard start codon of ATN, excepted for nad1 of N. vitripennis, which started with TTG, and terminated with a typical stop codon TAA/TAG or an incomplete stop codon T. All transfer RNA (tRNA) genes were predicted to fold into the typical clover-leaf secondary structures, except for trnS1, which lacks the DHU arm in all species. In P. vindemmiae, trnR and trnQ lack the DHU arm and TΨC arm, respectively. Although most genes evolved under a strong purifying selection, the Ka/Ks value of the atp8 gene of P. vindemmiae was greater than 1, indicating putative positive selection. A novel transposition of trnR in P. vindemmiae was revealed, which was the first of this kind to be reported in Pteromalidae. Two kinds of datasets (PCG12 and AA) and two inference methods (maximum likelihood and Bayesian inference) were used to reconstruct a phylogenetic hypothesis for the newly sequenced mitogenomes of Pteromalidae and those deposited in GenBank. The topologies obtained recovered the monophyly of the three subfamilies included. Pachyneurinae and Pteromalinae were recovered as sister families, and both appeared sister to Sycophaginae. The pairwise breakpoint distances of mitogenome rearrangements were estimated to infer phylogeny among pteromalid species. The topology obtained was not totally congruent with those reconstructed using the ML and BI methods.

The Leg Sensilla of Insects from Different Habitats-Comparison of Strictly Aquatic and Riparian Bugs (Corixidae, Ochteridae, Gelastocoridae: Nepomorpha: Insecta: Heteroptera)

The aim of this study was to analyze morphological types and arrangement of the leg sensilla of Corixidae, Ochteridae and Gelastocoridae, in relation to their various habitats. The leg sensilla of four species of Corixidae, six of Gelastocoridae and two of Ochteridae were studied. Eight main types of sensilla with six subtypes of sensilla trichodea and four subtypes of sensilla chaetica were found and described. The greatest variability was observed among mechanoreceptive sensilla. The study showed differences in the shape of the legs between strictly aquatic and terrestrial taxa. It is the first attempt to describe leg sensilla among nepomorphan taxa.

Phylogeny and historical biogeography of the water boatmen (Insecta: Hemiptera: Heteroptera: Nepomorpha: Corixoidea)

The water boatmen of Corixoidea, a group of aquatic bugs with more than 600 extant species, is one of the largest superfamilies of Nepomorpha. Contrary to the other nepomorphan lineages, the Corixoidea are most diverse in the Laurasian remnant Holarctic region. To explicitly test whether the present-day Holarctic distribution of diverse corixids is associated with the arising of the Laurasian landmass that was separated from Gondwana, we investigated the phylogeny, divergence times and historical biogeography of Corixoidea based on morphological and molecular characters sampled from 122 taxa representing all families, subfamilies, tribes and approximately 54 % of the genera. Our results were largely congruent with the phylogenetic relationships within the established nepomorphan phylogenetic context. The fossil calibrated chronogram, diversification analysis and ancestral ranges reconstruction indicated that Corixoidea began to diversify in Gondwana in the late Triassic approximately at 224 Ma and the arising of the most diverse subfamily Corixinae in Corixidae in the Holarctic region was largely congruent with the time of separation of Laurasia from Gondwana. The large-scale expansion of the temperate and cold zones on the northward-moving Laurasian landmass after the breakup of the Pangea provided new aquatic niches and ecological opportunities for promoting rapid diversification for the Holarctic corixid lineage.

A synopsis of the numbers of testicular follicles and ovarioles in true bugs (Heteroptera, Hemiptera) - sixty-five years of progress after J. Pendergrast's review

The structure of testes and ovaries can be described in its simplest form by the number of follicles and ovarioles they contain. Sixty-five years after the last review of the internal reproductive systems in true bugs (Heteroptera), the data accumulated today on the number of testicular follicles and ovarioles in their gonads are summarized. In addition, data on the number and type (mesadenia/ectadenia) of accessory glands are given. The hemipteran suborder Heteroptera constitutes one of the most diverse groups of non-homometabolous ('Hemimetabola') insects, comprising more than 40,000 described species worldwide and approximately 100 families, classified into seven infraorders. Data are available for all infraorders; however, more than 90% of studied species belong to the largest and most evolutionarily derived infraorders Cimicomorpha and Pentatomomorpha. In true bugs, in general, the number of follicles varies from one to nine (in a testis), and the number of ovarioles varies from two to 24 (in an ovary). Seven follicles per testis and seven ovarioles per ovary prevail being found in approximately 43.5% (307 species) and 24.4% (367 species) of studied species, respectively. Such a structure of testes and ovaries is considered an ancestral character state in the Heteroptera. In the evolution of this group, the number of follicles and ovarioles both increased and decreased, but the trend towards a decrease clearly prevailed.

The first complete mitochondrial genome of Dufouriellini (Hemiptera: Anthocoridae) and implications for its phylogenetic position

The mitochondrial genome (mitogenome) is extensively used to better understand the phylogenetic relationships within the family level, but there are still limited representations at the tribe level of Anthocoridae. Here we describe the first complete mitogenome of Dufouriellini. The mitogenome of Cardiastethus sp. is 15,209 bp in size, containing 13 typical protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and a control region. All genes are arranged in the same gene order as the most other known cimicomorphan mitogenomes. The phylogenetic relationships based on mitogenomes using Bayesian inference and maximum likelihood methods show that Dufouriellini is sister to Anthocorini, and then both of them together form sister group with Oriini. The monophyly of each superfamily of Cimicomorpha is generally well supported. Reduvioidea is basal within Cimicomorpha. The topology of the remaining superfamily is as follows: (Miroidea + (Cimicoidea + (Velocipedoidea + Nabioidea))). This study will help to enhance our understanding of mitochondrial genomic evolution and phylogenetic relationships in the tribe level of Anthocoridae and also superfamily level of Cimicomorpha.

Genome Organization and Comparative Evolutionary Mitochondriomics of Brown Planthopper, Nilaparvata lugens Biotype 4 Using Next Generation Sequencing (NGS)

Nilaparvata lugens is the main rice pest in India. Until now, the Indian N. lugens mitochondrial genome has not been sequenced, which is a very important basis for population genetics and phylogenetic evolution studies. An attempt was made to sequence two examples of the whole mitochondrial genome of N. lugens biotype 4 from the Indian population for the first time. The mitogenomes of N. lugens are 16,072 and 16,081 bp long with 77.50% and 77.45% A + T contents, respectively, for both of the samples. The mitochondrial genome of N. lugens contains 37 genes, including 13 protein-coding genes (PCGs) (cox1-3, atp6, atp8, nad1-6, nad4l, and cob), 22 transfer RNA genes, and two ribosomal RNA (rrnS and rrnL) subunits genes, which are typical of metazoan mitogenomes. However, both samples of N. lugens mitogenome in the present study retained one extra copy of the trnC gene. Additionally, we also found 93 bp lengths for the atp8 gene in both of the samples, which were 60–70 bp less than that of the other sequenced mitogenomes of hemipteran insects. The phylogenetic analysis of the 19 delphacids mitogenome dataset yielded two identical topologies when rooted with Ugyops sp. in one clade, and the remaining species formed another clade with P. maidis and M. muiri being sisters to the remaining species. Further, the genus Nilaparvata formed a separate subclade with the other genera (Sogatella, Laodelphax, Changeondelphax, and Unkanodes) of Delphacidae. Additionally, the relationship among the biotypes of N. lugens was recovered as the present study samples (biotype-4) were separated from the three biotypes reported earlier. The present study provides the reference mitogenome for N. lugens biotype 4 that may be utilized for biotype differentiation and molecular-aspect-based future studies of N. lugens.

Variation in the Mitochondrial Genome of the Chagas Disease Vector Triatoma infestans (Hemiptera: Reduviidae)

Chagas' disease is transmitted mainly by members of the subfamily Triatominae (Hemiptera: Reduviidae). Among them, Triatoma infestans (Klug) is the main vector of the disease in Southern Cone of Latin America. In order to contribute to knowledge of the genetic variation between triatomine vectors, in the present study, we analyzed the intraspecific and interspecific variations of the seven mitogenomes available from Triatominae. In addition, in order to examine their evolutionary relationships with others species of Reduviidae and to estimate the divergence time of the main lineages, we constructed phylogenetic trees including mitogenome sequences of 30 species from Reduviidae. Comparative analysis between mitochondrial DNA sequences from two specimens of T. infestans revealed a total of 54 variable sites. Triatoma infestans, Triatoma dimidiata (Latreille), Triatoma rubrofasciata (De Geer), Triatoma migrans (Breddin), Rhodnius pictipes (Stål), and Panstrongylus rufotuberculatus (Champion) present similar mitogenome organization and the length differences observed among these species are primarily caused by variations in control region (CR) and intergenic spacers (IGS). The relative synonymous codon usage values (RSCU) were similar in the six species of Triatominae, and in agreement with the observed in other insects, a biased use of A and C nucleotides in the majority strand was detected. The monophyly of five subfamilies was strongly supported (Phymatinae, Peiratinae, Triatominae, Stenopodainae, and Harpactorinae), while the sampled species of Reduviinae were grouped with one specie from the Salyavatinae subfamily. The oldest subfamily is Phymatinae at 100.3 Mya (99.6-102.2 Mya) and the youngest is Triatominae and Stenopodainae at 52.6 Mya (42.5-63.7 Mya). The estimated diversification time for the Triatominae subfamily agrees with the Andean uplift geological event. An analysis with more mitogenomes from more Triatominae species would be necessary to provide sufficient evidence to support this finding.

Morphology of the Antennal Sensilla of Notonectoidea and Comparison of Evolutionary Changes in Sensilla Types and Distribution in Infraorder Nepomorpha (Insecta: Heteroptera)

Simple Summary

Antennal sensilla are important sensory organs for insects. According to their morphological structures, they respond to different chemical or mechanical stimuli. The antennae of the studied families of water insects (Notonectidae, Pleidae and Helotrephidae) are short and concealed under the head, leaving a small amount of space for the existence of sensory structures. Nevertheless, six main types of sensilla have been discovered on the surfaces of these antennae. The morphological types described in this study were further compared with other studies on the antennal sensilla of water bugs (Nepomorpha) in order to compare their evolutionary changes within the group.

Abstract

This article introduces the results of a study of three families of Nepomorpha and is the last part of a series of studies that sums up our work on the morphologies of the antennal sensory structures in this taxon. The morphologies and distribution of the sensilla in the families Notonectidae, Pleidae and Helotrephidae were studied under a scanning electron microscope. Six main types (sensilla trichodea, chaetica, campaniformia, basiconica, ampullacea and coeloconica) and ten subtypes (five subtypes of sensilla trichodea and five subtypes of sensilla basiconica) were described. The results were compared with other studies on the antennal sensilla of Nepomorpha in order to assess evolutionary changes within the infraorder. With the use of cladistics analysis, the monophyly of the families Nepidae, Micronectidae, Corixidae and Gelastocoridae was supported. On the other hand, the occurrence of some clades forming superfamilies was weakly supported by bootstrap analysis. These results, supported by presence of the numerous autapomorphies, suggest that antennal sensilla evolved within inner groups.

Complete mitochondrial genome sequencing of Oxycarenus laetus (Hemiptera: Lygaeidae) from two geographically distinct regions of India

Oxycarenus laetus is a seed-sap sucking pest affecting a variety of crops, including cotton plants. Rising incidence and pesticide resistance by O. laetus have been reported from India and neighbouring countries. In this study, O. laetus samples were collected from Bhatinda and Coimbatore (India). Pure mtDNA was isolated and sequenced using Illumina MiSeq. Both the samples were found to be identical species (99.9%), and the complete genome was circular (15,672 bp), consisting of 13 PCGs, 2 rRNA, 23 tRNA genes, and a 962 bp control region. The mitogenome is 74.1% AT-rich, 0.11 AT, and - 0.19 GC skewed. All the genes had ATN as the start codon except cox1 (TTG), and an additional trnT was predicted. Nearly all tRNAs folded into the clover-leaf structure, except trnS1 and trnV. The intergenic space between trnH and nad4, considered as a synapomorphy of Lygaeoidea, was displaced. Two 5 bp motifs AATGA and ACCTA, two tandem repeats, and a few microsatellite sequences, were also found. The phylogenetic tree was constructed using 36 mitogenomes from 7 super-families of Hemiptera by employing rigorous bootstrapping and ML. Ours is the first study to sequence the complete mitogenome of O. laetus or any Oxycarenus species. The findings from this study would further help in the evolutionary studies of Lygaeidae.

Comparative analysis of twelve mitogenomes of Caliscelidae (Hemiptera: Fulgoromorpha) and their phylogenetic implications

Here, the complete mitochondrial genomes (mitogenomes) of 12 Caliscelidae species, Augilina tetraina, Augilina triaina, Symplana brevistrata, Symplana lii, Neosymplana vittatum, Pseudosymplanella nigrifasciata, Symplanella brevicephala, Symplanella unipuncta, Augilodes binghami, Cylindratus longicephalus, Caliscelis shandongensis, and Peltonotellus sp., were determined and comparatively analyzed. The genomes varied from 15,424 to 16,746 bp in size, comprising 37 mitochondrial genes and an A+T-rich region. The typical gene content and arrangement were similar to those of most Fulgoroidea species. The nucleotide compositions of the mitogenomes were biased toward A/T. All protein-coding genes (PCGs) started with a canonical ATN or GTG codon and ended with TAN or an incomplete stop codon, single T. Among 13 PCGs in 16 reported Caliscelidae mitogenomes, cox1 and atp8 showed the lowest and highest nucleotide diversity, respectively. All PCGs evolved under purifying selection, with atp8 considered a comparatively fast-evolving gene. Phylogenetic relationships were reconstructed based on 13 PCGs in 16 Caliscelidae species and five outgroups using maximum likelihood and Bayesian inference analyses. All species of Caliscelidae formed a steadily monophyletic group with high support. Peltonotellini was present at the basal position of the phylogenetic tree. Augilini was the sister group to Caliscelini and Peltonotellini.

Complete Mitochondrial Genomes of Metcalfa pruinosa and Salurnis marginella (Hemiptera: Flatidae): Genomic Comparison and Phylogenetic Inference in Fulgoroidea

The complete mitochondrial genomes (mitogenomes) of two DNA barcode-defined haplotypes of Metcalfa pruinosa and one of Salurnis marginella (Hemiptera: Flatidae) were sequenced and compared to those of other Fulgoroidea species. Furthermore, the mitogenome sequences were used to reconstruct phylogenetic relationships among fulgoroid families. The three mitogenomes, including that of the available species of Flatidae, commonly possessed distinctive structures in the 1702-1836 bp A+T-rich region, such as two repeat regions at each end and a large centered nonrepeat region. All members of the superfamily Fulgoroidea, including the Flatidae, consistently possessed a motiflike sequence (TAGTA) at the ND1 and trnS₂ junction. The phylogenetic analyses consistently recovered the familial relationships of (((((Ricaniidae + Issidae) + Flatidae) + Fulgoridae) + Achilidae) + Derbidae) in the amino acid-based analysis, with the placement of Cixiidae and Delphacidae as the earliest-derived lineages of fulgoroid families, whereas the monophyly of Delphacidae was not congruent between tree-constructing algorithms.

The complete mitochondrial genome of Pentatoma rufipes (Hemiptera, Pentatomidae) and its phylogenetic implications

Pentatoma rufipes (Linnaeus, 1758) is an important agroforestry pest widely distributed in the Palaearctic region. In this study, we sequence and annotate the complete mitochondrial genome of P. rufipes and reconstruct the phylogenetic trees for Pentatomoidea using existing data for eight families published in the National Center for Biotechnology Information database. The mitogenome of P. rufipes is 15,887-bp-long, comprising 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a control region, with an A+T content of 77.7%. The genome structure, gene order, nucleotide composition, and codon usage of the mitogenome of P. rufipes were consistent with those of typical Hemiptera insects. Among the protein-coding genes of Pentatomoidea, the evolutionary rate of ATP8 was the fastest, and COX1 was found to be the most conservative gene in the superfamily. Substitution saturation assessment indicated that neither transition nor transversion substitutions were saturated in the analyzed datasets. Phylogenetic analysis using the Bayesian inference method showed that P. rufipes belonged to Pentatomidae. The node support values based on the dataset concatenated from protein-coding and RNA genes were the highest. Our results enrich the mitochondrial genome database of Pentatomoidea and provide a reference for further studies of phylogenetic systematics.

Characterization of Three Complete Mitogenomes of Flatidae (Hemiptera: Fulgoroidea) and Compositional Heterogeneity Analysis in the Planthoppers' Mitochondrial Phylogenomics

Although sequences of mitogenomes have been widely used for investigating phylogenetic relationship, population genetics, and biogeography in many members of Fulgoroidea, only one complete mitogenome of a member of Flatidae has been sequenced. Here, the complete mitogenomes of Cerynia lineola, Cromna sinensis, and Zecheuna tonkinensis are sequenced. The gene arrangements of the three new mitogenomes are consistent with ancestral insect mitogenomes. The strategy of using mitogenomes in phylogenetics remains in dispute due to the heterogeneity in base composition and the possible variation in evolutionary rates. In this study, we found compositional heterogeneity and variable evolutionary rates among planthopper mitogenomes. Phylogenetic analysis based on site-homogeneous models showed that the families (Delphacidae and Derbidae) with high values of Ka/Ks and A + T content tended to fall together at a basal position on the trees. Using a site-heterogeneous mixture CAT + GTR model implemented in PhyloBayes yielded almost the same topology. Our results recovered the monophyly of Fulgoroidea. In this study, we apply the heterogeneous mixture model to the planthoppers’ phylogenetic analysis for the first time. Our study is based on a large sample and provides a methodological reference for future phylogenetic studies of Fulgoroidea.

The complete mitochondrial genome sequence of Asclepios apicalis (Gerromorpha: Gerridae)

In this study, the complete mitochondrial genome of Asclepios apicalis was sequenced and assembled, which was first reported in Asclepios. The mitogenome of Asclepios apicalis was 15,391 bp in length, and it contained 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, 13 protein-coding genes (PCGs), and a control region (D-loop), the overall base nucleotide compositions encoded was 42.9% A, 14.3% C, 10.0% G, and 32.8% T.

The mitochondrial genome of a giant water bug Lethocerus deyrollei (Hemiptera: Belostomatidae) from South Korea

A giant water bug Lethocerus deyrollei (Hemiptera: Belostomatidae) is a large, predatory, and nocturnal hemipteran insect, which has been considered threatened and thus enrolled as an endangered species in South Korea and Japan. Here, we characterized the complete mitochondrial genome of L. deyrollei, which has a circular form with 19,295 bp in length, which is the longest when compared to those of the 111 hemipteran species reported so far. Its longest genome size is due to the extremely extended CR (4686 bp), which is much longer than those of China and Japan. It consisted of a total of 37 genes (13 PCGs, 22 tRNA genes, and two rRNA genes) and one control region (CR). The genome composition and gene order were identical to those previously reported from the same species of China and Japan with over 99.7% sequence similarities except for CR and trnI. The nucleotide composition was highly A + T biased, accounting for 71% of the whole mitochondrial genome, as in other species of Nepoidea. Based on the aa sequences of 13 PCGs, we reconstructed a maximum likelihood tree, which indicated that the three mitochondrial genomes of L. deyrollei from South Korea, China, and Japan are grouped, and also Lethocerus, Belostomatidae, Nepoidea, Nepomorpha, Heteroptra are strong monophyletic groups, respectively.

The complete mitochondrial genome sequence of Aquarius elongatus (Hemiptera: Gerridae)

In this study, we report the complete mitochondrial genome of Aquarius elongatus. The mitogenome was 15,370 bp in length, comprising 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a control region. Maximum-likelihood (ML) phylogenetic tree indicated that Aquarius elongatus has a close relationship with Aquarius paludum. In general, this study provides meaningful genetic information for A. elongatus.

Comparative Study of Antennal Sensilla of Corixidae and Micronectidae (Hemiptera: Heteroptera: Nepomorpha: Corixoidea)

The goal of this study was to analyze the types and distributional patterns of sensilla in Corixoidea, which is part of the approach to the phylogeny study of Nepomorpha, based on the morphological characters of sensilla. This paper presents the results of the study, with the use of a scanning electron microscope (SEM), on the antennae of species from the families Corixidae and Micronectidae. The antennal sensilla of eleven species from Corixidae and two species from Micronectidae were studied. Five main types of sensilla with several subtypes of sensilla trichodea were found and described. The study has shown that the family Corixidae has a strong uniformity when it comes to antennal sensilla (similar patterns of sensilla trichodea and basiconica), and a similarity to the types and distributions of sensilla in two species of the family Micronectidae. However, significant differences between the families were also discovered (differences in sensilla presence on the first and second antennomeres, lack of sensilla coeloconica on the third antennomere in Micronectidae), which leads to a supportive conclusion of the systematic position of Micronectidae as a family.

Novel tRNA gene rearrangements in the mitochondrial genome of Camarochiloides weiweii (Hemiptera: Pachynomidae)

Pachynomidae is a poorly known heteropteran group currently considered as a distinct family closely related to Reduviidae. Taxonomic position of Pachynomidae has been debated for a long time. Knowledge about mitochondrial genome (mitogenome) of this family also remain limited. In the present study, we describe the complete mitogenome of Camarochiloides weiweii, the first mitogenome sequenced for this family. This mitogenome consists of 13 protein-coding genes, two rRNA genes and 21 tRNA genes. Two tRNA gene rearrangements in this genome, including the translocation of trnA and trnR as well as the transposition of trnH, were observed and could be explained by the tandem duplication-random loss (TDRL) model. Bayesian and maximum likelihood analyses based on sequences of three different mitogenomic datasets recovered the monophyly of Reduvioidea, and the sister relationship between Reduvioidea and the remaining Cimicomorpha lineages. This study enlightens the future application of mitogenomic data in the phylogenetic study of Cimicomorpha based on the large-scale sampling.

Phylogeny and diversification of the true water bugs (Insecta: Hemiptera: Heteroptera: Nepomorpha)

Climate fluctuations and tectonic reconfigurations associated with environmental changes play large roles in determining patterns of adaptation and diversification, but studies documenting how such drivers have shaped the evolutionary history and diversification dynamics of limnic organisms during the Mesozoic are scarce. Members of the heteropteran infraorder Nepomorpha, or aquatic bugs, are ideal for testing the effects of these determinants on their diversification pulses because most species are confined to aquatic environments during their entire life. The group has a relatively mature taxonomy and is well represented in the fossil record. We investigated the evolution of Nepomorpha based on phylogenetic analyses of morphological and molecular characters sampled from 115 taxa representing all 13 families and approximately 40% of recognized genera. Our results were largely congruent with the phylogenetic relationships inferred from morphology. A divergence dating analysis indicated that Nepomorpha began to diversify in the late Permian (approximately 263 Ma), and diversification analyses suggested that palaeoecological opportunities probably promoted lineage diversification in this group.

Complete mitochondrial genome of the spotted lanternfly, Lycorma delicatula White, 1845 (Hemiptera: Fulgoridae)

The spotted lanternfly, Lycorma delicatula White, 1845 (Hemiptera: Fulgoridae), is an invasive pest that attacks forest as well as agricultural trees. We sequenced the 15,798-bp long complete mitochondrial genome (mitogenome) of this species; it consists of a typical set of genes (13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes) and one major non-coding A + T-rich region. The orientation and gene order of the L. delicatula mitogenome are identical to that of the ancestral type found in majority of the insects. Bayesian inference (BI) and maximum-likelihood (ML) phylogeny placed the L. delicatula examined in our study, together with other geographical samples of the species in a group with the highest nodal support, forming the subfamily Aphaeninae to which L. delicatula belongs.

Aquatic Hemiptera in Southwest Cameroon: Biodiversity of Potential Reservoirs of Mycobacterium ulcerans and Multiple Wolbachia Sequence Types Revealed by Metagenomics

Buruli ulcer (BU), caused by Mycobacterium ulcerans, is a neglected tropical disease associated with freshwater habitats. A variety of limnic organisms harbor this pathogen, including aquatic bugs (Hemiptera: Heteroptera), which have been hypothesized to be epidemiologically important reservoirs. Aquatic Hemiptera exhibit high levels of diversity in the tropics, but species identification remains challenging. In this study, we collected aquatic bugs from emerging foci of BU in the Southwest Region of Cameroon, which were identified using morphological and molecular methods. The bugs were screened for mycobacterial DNA and a selection of 20 mycobacteria-positive specimens from the families Gerridae and Veliidae were subjected to next-generation sequencing. Only one individual revealed putative M. ulcerans DNA, but all specimens contained sequences from the widespread alpha-proteobacterial symbiont, Wolbachia. Phylogenetic analysis placed the Wolbachia sequences into supergroups A, B, and F. Circularized mitogenomes were obtained for seven gerrids and two veliids, the first from these families for the African continent. This study suggests that aquatic Hemiptera may have a minor role (if any) in the spread of BU in Southwest Cameroon. Our metagenomic analysis provides new insights into the incursion of Wolbachia into aquatic environments and generated valuable resources to aid molecular taxonomic studies of aquatic Hemiptera.

Tracking outbreak populations of the pepper weevil Anthonomus eugenii (Coleoptera; Curculionidae) using complete mitochondrial genomes

The pepper weevil, Anthonomus eugenii, is a major pest on Capsicum species. Apart from natural spread, there is a risk of spread via international pepper trade. In the Netherlands, a pepper weevil outbreak occurred in 2012 and affected six greenhouses producing different sweet pepper varieties. The following year, a pepper weevil outbreak occurred in Italy. To trace the origin of the Dutch outbreak and to establish if the Dutch and Italian outbreaks were linked, we determined the mitogenomes of A. eugenii specimens collected at outbreak locations, and compared these with specimens from the native area, and other areas where the pest was introduced either by natural dispersal or via trade. The circular 17,257 bp A. eugenii mitogenome comprises thirteen mitochondrial genes typically found in insect species. Intra-species variation of these mitochondrial genes revealed four main mitochondrial lineages encompassing 41 haplotypes. The highest diversity was observed for specimens from its presumed native area (i.e. Mexico). The Dutch outbreak specimens represented three highly similar haplotypes, suggesting a single introduction of the pest. The major Dutch haplotype was also found in two specimens from New Jersey. As the Netherlands does not have pepper trade with New Jersey, it is likely that the specimens sampled in New Jersey and those sampled in the Netherlands originate from a shared source that was not included in this study. In addition, our analysis shows that the Italian and Dutch outbreaks were not linked. The mitochondrial genome is a useful tool to trace outbreak populations and the methodology presented in this paper could prove valuable for other invasive pest species, such as the African fruit moth Thaumatotibia leucotreta and emerald ash borer Agrilus planipennis.

Contribution to the mitogenome diversity in Delphacinae: Phylogenetic and ecological implications

We document the complete (or nearly complete) mitogenomes of 20 Delphacidae taxa, and together with 17 other delphacid mitogenomes currently in GenBank, to reconstruct the phylogeny of the Delphacinae and to investigate mitogenome differences among members of Delphacini, Tropidocephalini and Saccharosydnini. The mitogenomes of the 20 species encode the complete set of 37 genes usually found in animal mitogenomes. The length of complete mitogenomes in Delphacinae ranges from 15,531 to 16,231 bp. The gene order of all newly sequenced mitogenomes are identical, and the mitogenome gene order of Stenocranus matsumurai Metcalf in Stenocraninae has a transposition of tRNA^Thr. The two-clade system in Tropidocephalini was supported with high value (PP = 1, BS = 100), and the monophyly of Bambusiphaga was recovered in this study. Finally, we found that the host shift from plants with a C₃ to a C₄ photosynthetic pathway appears to have occurred independently in several clades.

The complete mitochondrial genome of Tetraphleps aterrimus (Hemiptera: Anthocoridae): Genomic comparisons and phylogenetic analysis of Cimicomorpha

The mitochondrial genome (mitogenome) provides important information for better understanding the phylogenetic relationships within heteropteran infraorder Cimicomorpha (Hemiptera: Heteroptera), but there are still limited representations at the family level of Anthocoridae. Here we sequenced the complete mitogenome of Tetraphleps aterrimus. It is 15,803 bp in size, and contains the expected 37 genes (13 PCGs, 22 tRNAs and 2 rRNAs) and control region. Gene order is identical to that of typical cimicomorphans. In comparison with other cimicomorphans, the ratios of Ka/Ks are increasing from 0.17 for COI to 0.85 for ATP8, which demonstrates COI shows the lowest evolutionary rate, while ATP8 appears to be the highest. The ratios of conserved sites of COI is the highest, while ATP8 is the lowest, suggesting that the evolutionary rate of ATP8 is higher than COI. The phylogenetic relationships based on mitogenomes using Bayesian inference (BI) and Maximum likelihood (ML) methods show that Tetraphleps aterrimus is sister to (Orius niger + Orius sauteri), suggesting that Tetraphleps aterrimus belongs to Anthocoridae. The monophyly of each superfamily is generally well supported and Reduvioidea is placed as basal branch in Cimicomorpha. The results support the remaining superfamily groupings (Miroidea + (Cimicoidea + (Velocipedoidea + Nabioidea))).

A mitochondrial genome of Micronectidae and implications for its phylogenetic position

The mitochondrial genome (mitogenome) has been extensively used to better understand the phylogenetic relationships within the heteropteran infraorder Nepomorpha (Hemiptera), but no mitogenome in Micronectidae has been sequenced to date. Here we describe the first complete mitogenome of Micronecta sahlbergii (Jakovlev, 1881). The mitogenome is 15,005 bp in size, containing 13 typical PCGs, 22 tRNAs, two rRNAs and a control region (CR). All genes are arranged in the same gene order as the most other known heteropteran mitogenome. The phylogenetic relationships based on mitogenomes using Bayesian inference and Maximum likelihood methods showed that Micronecta sahlbergii was sister to Sigara septemlineata, suggesting that Micronecta sahlbergii belongs to Corixoidea. Corixoidea was basal within Nepomorpha. The PCG12 and PCG12RT matrices of BI and ML analyses yielded the consistent topology, respectively. Whereas there was no consistent conclusions in PCG123 and PCG123RT matrices. Saturation tests showed that PCG12 and PCG12RT were free of saturation in evaluation of transition and transversion substitution, while PCG123 and PCG123RT exhibited a plateau revealing saturation of transition suggesting that the third codon positions of PCGs were not suitable for addressing relationships at the superfamily level in Nepomorpha. So our results supported the phylogenetic analysis of PCG12 and PCG12RT in Nepomorpha.

The Complete Mitochondrial Genome of Ugyops sp. (Hemiptera: Delphacidae)

The complete mitochondrial genome (mitogenome) of Ugyops sp. (Hemiptera: Delphacidae) was sequenced, making it the first determined mitogenome from the subfamily Asiracinae, the basal clade of the family Delphacidae. The mitogenome was 15,259 bp in length with A + T content of 77.65% and contained 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a control region. The gene order was identical with that of the ancestral insect. The nucleotide composition analysis indicated that the whole mitogenome was strongly A-skewed (0.288) and highly C-skewed (-0.270). For PCGs on the J-strand, the AT skew was positive, and the GC skew was negative. All PCGs started with canonical ATN codons, except for cox1 and nad5, which used CTG and GTG as start codon, respectively. All tRNAs could fold into typical cloverleaf secondary structures, with the exception of trnS1 (AGN), in which the dihydrouridine arm was reduced to a simple loop. The control region included a poly-T stretch downstream of the small rRNA gene (rrnS), a subregion of higher A + T content and tandemly repeated sequence near trnI. The mitogenome of Ugyops sp. could be very helpful in exploring the diversity and evolution of mitogenomes in Delphacidae.

The first complete mitochondrial genome from the subfamily Phyllocephalinae (Heteroptera: Pentatomidae) and its phylogenetic analysis

The true bug, Gonopsis affinis, is an important pest in China. Here, we determined the complete mitogenome of G. affinis, which is the first for the subfamily Phyllocephalinae. This 16,011-basepair (bp) mitogenome comprises of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes, and an A + T-rich region. The gene order and the orientation are similar to those of other sequenced Hemipteran species. All PCGs start with ATN codons, except COI, and end with TAA. The 22 tRNAs have a typical cloverleaf secondary structure except trnS-Ser (AGN), which lacks a dihydrouridine (DHU) arm. Phylogenetic analyses highly supported the monophyly of each family, and confirmed the reasonable placement of G. affinis.

The first mitochondrial genome for the subfamily Asopinae (Heteroptera: Pentatomidae) and its phylogenetic implications

Here, we determined the complete mitogenome of Picromerus griseus, the first for the subfamily Asopinae. This mitochondrial genome contains 16,338 bp, with an A + T content of 71.69%, and contains 37 typical mitochondrial genes (13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and a control region). The genome size, gene arrangement, A + T content, codon usage and secondary structures of 22 transfer RNA genes of the P. griseus mitogenome were similar to those of other sequenced pentatomoids. Bayesian analyses performed using the mitogenomes of P. griseus and its relatives, including 14 taxa, confirmed the reasonable placement of P. griseus.

Complete mitochondrial genomes of Anopheles stephensi and An. dirus and comparative evolutionary mitochondriomics of 50 mosquitoes

To better understand the phylogeny and evolution of mosquitoes, the complete mitochondrial genome (mitogenome) of Anopheles stephensi and An. dirus were sequenced and annotated, and a total of 50 mosquito mitogenomes were comparatively analyzed. The complete mitogenome of An. stephensi and An. dirus is 1,5371 bp and 1,5406 bp long, respectively. The main features of the 50 mosquito mitogenomes are conservative: 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, positive AT-skew and negative GC-skew. The gene order trnA-trnR in ancestral insects is rearranged. All tRNA genes have the typical clover leaf secondary structure but tRNA Ser . The control regions are highly variable in size. PCGs show signals of purifying selection, but evidence for positive selection in ND2, ND4 and ND6 is found. Bayesian and Maximum Likelihood phylogenetic analyses based on all PCG nucleotides produce an identical tree topology and strongly support the monophyly of subgenera Cellia, Anopheles, Keterszia and Nyssorhynchus, the sister relationship of the subgenera Nyssorhynchus and Keterszia, and Cellia and Anopheles. The most recent ancestor of the genus Anopheles and Culicini + Aedini exited ~145 Mya ago. This is the first comprehensive study of mosquito mitogenomes, which are effective for mosquito phylogeny at various taxonomic levels.

Complete mitochondrial genome of Triatoma infestans (Hemiptera, Reduviidae, Triatominae), main vector of Chagas disease

The complete mitogenome of Triatoma infestans, the main vector of Chagas disease in South America, was obtained by combining Illumina and Sanger sequencing sequence data. The 17,301bp long genome contains 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and a control region. The number, order and orientation of mitochondrial genes are the same as in T. dimidiata, the only Triatominae mt genome published so far. The main differences between both mitogenomes are found in the control region and in the intergenic spacer between the nd1 gene and the tRNA-Ser. Comparative analysis with other Reduviidae species shows high conservation in the mt genome organization. Molecular phylogeny using all available complete mt genomes from Reduviidae species confirms the close relationship between Triatominae and Stenopodainae.

A Mitochondrial Genome of Rhyparochromidae (Hemiptera: Heteroptera) and a Comparative Analysis of Related Mitochondrial Genomes

The Rhyparochromidae, the largest family of Lygaeoidea, encompasses more than 1,850 described species, but no mitochondrial genome has been sequenced to date. Here we describe the first mitochondrial genome for Rhyparochromidae: a complete mitochondrial genome of Panaorus albomaculatus (Scott, 1874). This mitochondrial genome is comprised of 16,345 bp, and contains the expected 37 genes and control region. The majority of the control region is made up of a large tandem-repeat region, which has a novel pattern not previously observed in other insects. The tandem-repeats region of P. albomaculatus consists of 53 tandem duplications (including one partial repeat), which is the largest number of tandem repeats among all the known insect mitochondrial genomes. Slipped-strand mispairing during replication is likely to have generated this novel pattern of tandem repeats. Comparative analysis of tRNA gene families in sequenced Pentatomomorpha and Lygaeoidea species shows that the pattern of nucleotide conservation is markedly higher on the J-strand. Phylogenetic reconstruction based on mitochondrial genomes suggests that Rhyparochromidae is not the sister group to all the remaining Lygaeoidea, and supports the monophyly of Lygaeoidea.

The first complete mitochondrial genome of a Belostomatidae species, Lethocerus indicus, the giant water bug: An important edible insect

Lethocerus indicus of the family Belostomatidae is one of the most preferred and delicious edible insects in different parts of South-East Asia including North-East, India. The mitogenome of L. indicus represents the first complete mitogenome sequence of a Belostomatidae species in Heteroptera order. The mitogenome of L. indicus is 16,251bp and contains 37 genes including 13 protein coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a large non-coding region. The genome has a typical gene order which is identical to other Heteroptera species. All tRNAs exhibit the classic cloverleaf secondary structure except tRNASer (AGN). All the PCGs employ a complete translation termination codon either TAA or TAG except COII. The nucleotide composition showed heavy biased toward AT accounting to 70.9% of total mitogenome. The overall A+T content of L. indicus mitogenome was comparatively lower than some other Heteropteran bugs mitogenomes. The control region is divided into seven different parts which includes the putative stem loop, repeats, tandem repeats, GC and AT rich regions. The phylogenetic relationship based on maximum-likelihood method using all protein coding genes was congruent with the traditional morphological classification that Belostomatidae is closely related to Nepidae. The complete mitogenome sequence of L. indicus provides fundamental data useful in conservation genetics and aquaculture diversification.

Duplication and Remolding of tRNA Genes in the Mitochondrial Genome of Reduvius tenebrosus (Hemiptera: Reduviidae)

Most assassin bugs are predators that act as important natural enemies of insect pests. Mitochondrial (mt) genomes of these insects are double-strand circular DNAs that encode 37 genes. In the present study, we explore the duplication and rearrangement of tRNA genes in the mt genome of Reduvius tenebrosus, the first mt genome from the subfamily Reduviinae. The gene order rearranges from CR (control region)-trnI-trnQ-trnM-ND2 to CR-trnQ-trnI2-trnI1-trnM-ND2. We identified 23 tRNA genes, including 22 tRNAs commonly found in insects and an additional trnI (trnI2), which has high sequence similarity to trnM. We found several pseudo genes, such as pseudo-trnI, pseudo-CR, and pseudo-ND2, in the hotspot region of gene rearrangement (between the control region and ND2). These features provided evidence that this novel gene order could be explained by the tandem duplication/random loss (TDRL) model. The tRNA duplication/anticodon mutation mechanism further explains the presence of trnI2, which is remolded from a duplicated trnM in the TDRL process (through an anticodon mutation of CAT to GAT). Our study also raises new questions as to whether the two events proceed simultaneously and if the remolded tRNA gene is fully functional. Significantly, the duplicated tRNA gene in the mitochondrial genome has evolved independently at least two times within assassin bugs.

Molecular Species Delimitation and Morphology of Aquatic and Sub-Aquatic Bugs (Heteroptera) in Cameroon

Aquatic and semi-aquatic bugs (Heteroptera) represent a remarkable diversity and a resurging interest has been given to documenting at the species level these insects inhabiting Cameroon in Central Africa due to their potential implication in the transmission of the bacterium Mycobacterium ulcerans, the causal agent of Buruli ulcer, an emerging human disease. A survey was carried out over two years in Cameroon. Morphological analyses were done in two steps. A first step consisted in separating the specimens based on broadly shared characters into morphotypes. The specimens were then separated into two independent batches containing each the same representation of each morphotype. One batch (309 specimens) was used by taxonomy experts on aquatic bugs for species level identification and/or to reconcile nymph with their corresponding adult species. The second batch (188 specimens) was used to define species based on the COI DNA sequences (standard sequence used for "DNA barcoding") and using the Automatic Barcode Gap Discovery (ABGD) method. The first morphological analysis step separated the specimens into 63 different morphotypes (49 adults and 14 nymphs), which were then found to belong to 54 morphological species in the infra-orders Gerromorpha and Nepomorpha based on the species-level morphological identification, and 41-45 putative molecular species according to the gap value retained in the ABGD. Integrating morphology and "DNA barcoding" reconciled all the specimens into 62 aquatic bug species in Cameroon. Generally, we obtained a good congruence between species a priori identified based on morphology from adult morphotypes and molecular putative species. Moreover, molecular identification has allowed the association of 86% of nymphs with adults. This work illustrates the importance of integrative taxonomy.

Phylogenetic relationships of Hemiptera inferred from mitochondrial and nuclear genes

Here, we reconstructed the Hemiptera phylogeny based on the expanded mitochondrial protein-coding genes and the nuclear 18S rRNA gene, separately. The differential rates of change across lineages may associate with long-branch attraction (LBA) effect and result in conflicting estimates of phylogeny from different types of data. To reduce the potential effects of systematic biases on inferences of topology, various data coding schemes, site removal method, and different algorithms were utilized in phylogenetic reconstruction. We show that the outgroups Phthiraptera, Thysanoptera, and the ingroup Sternorrhyncha share similar base composition, and exhibit "long branches" relative to other hemipterans. Thus, the long-branch attraction between these groups is suspected to cause the failure of recovering Hemiptera under the homogeneous model. In contrast, a monophyletic Hemiptera is supported when heterogeneous model is utilized in the analysis. Although higher level phylogenetic relationships within Hemiptera remain to be answered, consensus between analyses is beginning to converge on a stable phylogeny.

The complete mitochondrial genome of Nepa hoffmanni (Hemiptera: Heteroptera: Nepidae)

The complete mitochondrial genome (mt-genome) of Nepa hoffmanni has been reported in this study. This mitochondrial genome is 15 774 bp long, with an A + T content of 72.04%, containing the typical 37 genes (13 protein-coding genes (PCGs), 22 transfer RNA genes, and two ribosomal RNA genes) and a control region. All genes are arranged in the same gene order as most other known heteropteran mt-genome. This is the second completely sequenced mt-genome from the family Nepidae of Nepomorpha. Bayesian analyses were performed using the mt-genome of Nepa hoffmanni and its relatives, including 17 taxa, showing a reasonable placement of Nepa hoffmanni.

Efficient capture of natural history data reveals prey conservatism of cryptic termite predators

Stenophagy, specialization of a clade on a narrow range of taxa, has not been well studied in speciose clades of predators, principally due to the difficulty of obtaining adequate natural history data. The pantropical Salyavatinae (Hemiptera: Reduviidae; 17 genera, 107 species) contains members with enigmatic morphology and specialized behavior for feeding on termites. All Salyavatinae are suspected specialist termite predators; however, existing observations are limited to seven species. Prior analyses indicate that Salyavatinae may be paraphyletic with respect to another subfamily, Sphaeridopinae, also hypothesized to feed on termites. A molecular phylogeny of these putative termite assassins is here constructed using seven loci from 28 species in nine genera and is used in a dating analysis to shed light on the timing of Neotropical colonization by this primarily Old World clade. DNA extracted from gut contents of 50 individuals was assayed using PCR with prey-specific primers.Molecular assays, along with recent photographs and observations, provide substantial evidence that this clade feeds specifically upon termites, documenting 28 new individual associations. Our phylogeny supports a sister group relationship of the Neotropical genus Salyavata with Sphaeridopinae. Termite association data combined with our phylogeny provide evidence of previously unknown prey conservatism among clades of one of the most diverse groups of specialist termite predators.

Hemipteran mitochondrial genomes: features, structures and implications for phylogeny

The study of Hemipteran mitochondrial genomes (mitogenomes) began with the Chagas disease vector, Triatoma dimidiata, in 2001. At present, 90 complete Hemipteran mitogenomes have been sequenced and annotated. This review examines the history of Hemipteran mitogenomes research and summarizes the main features of them including genome organization, nucleotide composition, protein-coding genes, tRNAs and rRNAs, and non-coding regions. Special attention is given to the comparative analysis of repeat regions. Gene rearrangements are an additional data type for a few families, and most mitogenomes are arranged in the same order to the proposed ancestral insect. We also discuss and provide insights on the phylogenetic analyses of a variety of taxonomic levels. This review is expected to further expand our understanding of research in this field and serve as a valuable reference resource.

Oxygen-limited thermal tolerance is seen in a plastron-breathing insect and can be induced in a bimodal gas exchanger

Thermal tolerance has been hypothesized to result from a mismatch between oxygen supply and demand. However, the generality of this hypothesis has been challenged by studies on various animal groups, including air-breathing adult insects. Recently, comparisons across taxa have suggested that differences in gas exchange mechanisms could reconcile the discrepancies found in previous studies. Here, we test this suggestion by comparing the behaviour of related insect taxa with different gas exchange mechanisms, with and without access to air. We demonstrate oxygen-limited thermal tolerance in air-breathing adults of the plastron-exchanging water bug Aphelocheirus aestivalis. Ilyocoris cimicoides, a related, bimodal gas exchanger, did not exhibit such oxygen-limited thermal tolerance and relied increasingly on aerial gas exchange with warming. Intriguingly, however, when denied access to air, oxygen-limited thermal tolerance could also be induced in this species. Patterns in oxygen-limited thermal tolerance were found to be consistent across life-history stages in these insects, with nymphs employing the same gas exchange mechanisms as adults. These results advance our understanding of oxygen limitation at high temperatures; differences in the degree of respiratory control appear to modulate the importance of oxygen in setting tolerance limits.

Summary: Oxygen does not appear to universally limit thermal tolerance, but instead, oxygen-limited thermal tolerance is context dependent, being related to a species’ capacity to regulate oxygen consumption.

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

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.

Two nearly complete mitogenomes of wheat stem borers, Cephus pygmeus (L.) and Cephus sareptanus Dovnar-Zapolskij (Hymenoptera: Cephidae): an unusual elongation of rrnS gene

Two nearly complete mitochondrial genomes (mitogenomes) of wheat stem borers, Cephus pygmeus and Cephus sareptanus (Hymenoptera: Cephidae), were sequenced, characterised and compared with the previously known mitogenome of Cephus cinctus. The gene orders are mostly conserved, except for translocation of trnM and swapped position of trnI and trnQ. An A+T bias was found, but a deviation from strand asymmetry was also detected on the J strand. All protein coding genes (PCGs) are initiated by ATN codons, except for nad1, nad2 and atp8, and all are terminated with TAA, TA- or T- as a stop codon. The predicted secondary structures of rrnS and rrnL genes are mostly consistent with reported hymenopteran species. However, an unusual elongation in rrnS, not know elsewhere in the order, was discovered in Cephus species. Three autonomous sequences detected in domains I and II are mainly responsible for the length expansions.

Shotgun assembly of the assassin bug Brontostoma colossus mitochondrial genome (Heteroptera, Reduviidae)

The complete mitochondrial genome of the assassin bug Brontostoma colossus (Distant, 1902) (Heteroptera: Reduviidae) has been sequenced using a genome-skimming approach on an Illumina Hiseq 2000 platform. Fifty-four additional heteropteran mitogenomes, including five assassin bug species, were retrieved to allow for comparisons and phylogenetic analyses. The mitochondrial genome of B. colossus was determined to be 16,625 bp long, and consists of 13 protein-coding genes (PCGs), 23 transfer-RNA genes (tRNAs), two ribosomal-RNA genes (rRNAs), and one control region. The nucleotide composition is biased toward adenine and thymine (A+T=73.4%). Overall, architecture, nucleotide composition and genome asymmetry are similar among all available assassin bug mitogenomes. All PCGs have usual start-codons (Met and Ile). Three T and two TA incomplete termination codons were identified adjacent to tRNAs, which was consistent with the punctuation model for primary transcripts processing followed by 3' polyadenylation of mature mRNA. All tRNAs exhibit the classic clover-leaf secondary structure except for tRNASer(AGN) in which the DHU arm forms a simple loop. Two notable features are present in the B. colossus mitogenome: (i) a 131 bp duplicated unit including the complete tRNAArg gene, resulting in 23 potentially functional tRNAs in total, and (ii) a 857 bp duplicated region comprising 277 bp of the srRNA gene and 580 bp of the control region. A phylogenetic analysis based on 55 true bug mitogenomes confirmed that B. colossus belongs to Reduviidae, but contradicted a widely accepted hypothesis. This highlights the limits of phylogenetic analyses based on mitochondrial data only.