The first mitogenome of the subfamily Stenoponiinae (Siphonaptera: Ctenophthalmidae) and implications for its phylogenetic position

Fleas are the most important insect vectors that parasitize warm-blooded animals and are known vectors of zoonotic pathogens. A recent study showed that Stenoponia polyspina parasitizing Eospalax baileyi in Zoige County have carried Bartonella spp. and Spotted fever group Rickettsia (SFGR). Accurate identification and differentiation of fleas are essential for prevention and control of zoonotic pathogens. To understand phylogenetic relationship of the subfamily Stenoponiinae, we described morphological characteristics of S. polyspina and sequenced its mitogenome with 14,933 bp in size and high A + T content (~ 79%). The S. polyspina mitogenome retained the ancestral pattern of mitochondrial gene arrangement of arthropods without rearrangement. The start codons of 13 protein-coding genes (PCGs) are traditional ATN and the stop codons are TAA or TAG. Anticodon loops of all tRNA genes were 7 bp except for trnL2 and trnD had anticodon loops with 9 bp and the abnormal anticodon loops may be associated with frameshifting mutation. Genetic distance and Ka/Ks ratios indicated that all 13 PCGs of S. polyspina were subjected to purifying selection, with cox1 at the slowest rate and atp8 at the fastest rate. The mitogenomes of 24 species representing 7 families in the order Siphonaptera were selected to reconstruct phylogenetic tree based on concatenated nucleotide sequences of two datasets (PCGRNA matrix and PCG12RNA matrix) using Bayesian inference (BI) and Maximum likelihood (ML) methods. Phylogenetic tree supported that the superfamilies Ceratophylloidea, Vermipsylloidea, Pulicoidea were monophyletic and the superfamily Hystrichopsylloidea was paraphyletic. The family Ctenophthalmidae was monophyletic in PCGRNA-ML (codon partition) and paraphyletic in the remain trees. S. polyspina belongs to the subfamily Stenoponiinae was closely more related to the subfamily Rhadinopsyllinae. This paper explored phylogenetic position of diverse clades within the order Siphonaptera based on morphological and mitogenome data of S. polyspina. Our research enriched NCBI database of the order Siphonaptera.

The mitochondrial genome of Bottapotamon fukienense (Brachiura: Potamidae) is fragmented into two chromosomes

BACKGROUND

China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies.

RESULTS

Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae.

CONCLUSIONS

This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures.

Substantial rearrangements, single nucleotide frameshift deletion and low diversity in mitogenome of Wolbachia-infected strepsipteran endoparasitoid in comparison to its tephritid hosts

Insect mitogenome organisation is highly conserved, yet, some insects, especially with parasitic life cycles, have rearranged mitogenomes. Furthermore, intraspecific mitochondrial diversity can be reduced by fitness-affecting bacterial endosymbionts like Wolbachia due to their maternal coinheritance with mitochondria. We have sequenced mitogenomes of the Wolbachia-infected endoparasitoid Dipterophagus daci (Strepsiptera: Halictophagidae) and four of its 22 known tephritid fruit fly host species using total genomic extracts of parasitised flies collected across > 700 km in Australia. This halictophagid mitogenome revealed extensive rearrangements relative to the four fly mitogenomes which exhibited the ancestral insect mitogenome pattern. Compared to the only four available other strepsipteran mitogenomes, the D. daci mitogenome had additional transpositions of one rRNA and two tRNA genes, and a single nucleotide frameshift deletion in nad5 requiring translational frameshifting or, alternatively, resulting in a large protein truncation. Dipterophagus daci displays an almost completely endoparasitic life cycle when compared to Strepsiptera that have maintained the ancestral state of free-living adults. Our results support the hypothesis that the transition to extreme endoparasitism evolved together with increased levels of mitogenome changes. Furthermore, intraspecific mitogenome diversity was substantially smaller in D. daci than the parasitised flies suggesting Wolbachia reduced mitochondrial diversity because of a role in D. daci fitness.

Species Diversity and Phylogenetic Relationships of Olive Lace Bugs (Hemiptera: Tingidae) Found in South Africa

Olive lace bugs (Hemiptera: Tingidae) are small sap-sucking insects that feed on wild and cultivated Olea europaea. The diversity of olive lace bug species in South Africa, the most important olive producer on the continent, has been incompletely surveyed. Adult specimens were collected in the Western Cape province for morphological and DNA-based species identification, and sequencing of complete mitogenomes. Cysteochila lineata, Plerochila australis, Neoplerochila paliatseasi and Neoplerochila sp. were found at 12 sites. Intra- and interspecific genetic divergences and phylogenetic clustering in 30 species in 18 genera of Tingidae using new and publicly available DNA barcodes showed high levels of congruity between taxonomic and genetic data. The phylogenetic position of the four species found in South Africa was inferred using new and available mitogenomes of Tingidae. Notably, olive lace bugs formed a cluster of closely related species. However, Cysteochila was non-monophyletic as C. lineata was recovered as a sister species to P. australis whereas Cysteochila chiniana, the other representative of the genus, was grouped with Trachypeplus jacobsoni and Tingis cardui in a different cluster. This result suggests that feeding on O. europaea may have a common origin in Tingidae and warrants future research on potential evolutionary adaptations of olive lace bugs to this plant host.

Complete mitochondrial genomes of Bittacus strigosus and Panorpa debilis and genomic comparisons of Mecoptera

Mitochondrial genomes play a significant role in reconstructing phylogenetic relationships and revealing molecular evolution of insects. However, only four mitochondrial genomes were reported in Mecoptera to date. Here, we obtained two new complete mitochondrial genomes of the hangingfly Bittacus strigosus Hagen, 1861 and the scorpionfly Panorpa debilis Westwood, 1846. The results show that the complete mitogenome sequences of B. strigosus and P. debilis are 15,825 and 17,018 bp, respectively, both containing 37 genes and one control region. The mecopteran mitogenomes are highly similar in A + T bias, AT-skew, and GC-skew. Tandem repeats of the control region were discovered in Mecoptera for the first time. The sliding window, genetic distance, and Ka/Ks ratio analyses indicate the purifying selection of 13 protein-coding genes, the lowest evolutionary rate of cox1, and the highest sequence variability of atp8. Considering the sufficiently large size, fast evolution, and high ratio of Ka/Ks, nad4L and nad6 are regarded as potential markers for future phylogenetic analyses, population genetics, and species delimitations in Mecoptera. The phylogenetic relationships were reconstructed for four families of Mecoptera based on all six available mitogenomes using Bayesian inference and maximum likelihood methods. The phylogeny is presented as Boreidae + (Nannochoristidae + (Bittacidae + Panorpidae)).

Mitochondrial DNA and their nuclear copies in the parasitic wasp Pteromalus puparum: A comparative analysis in Chalcidoidea

Chalcidoidea (chalcidoid wasps) are an abundant and megadiverse insect group with both ecological and economical importance. Here we report a complete mitochondrial genome in Chalcidoidea from Pteromalus puparum (Pteromalidae). Eight tandem repeats followed by 6 reversed repeats were detected in its 3308 bp control region. This long and complex control region may explain failures of amplifying and sequencing of complete mitochondrial genomes in some chalcidoids. In addition to 37 typical mitochondrial genes, an extra identical isoleucine tRNA (trnI) was detected at the opposite end of the control region. This recent mitochondrial gene duplication indicates that gene arrangements in chalcidoids are ongoing. A comparison among available chalcidoid mitochondrial genomes reveals rapid gene order rearrangements overall and high protein substitution rates in most chalcidoid taxa. In addition, we identified 24 nuclear sequences of mitochondrial origin (NUMTs) in P. puparum, summing up to 9989 bp, with 3617 bp of these NUMTs originating from mitochondrial coding regions. NUMTs abundance in P. puparum is only one-twelfth of that in its relative, Nasonia vitripennis. Based on phylogenetic analysis, we provide evidence that a faster nuclear degradation rate contributes to the reduced NUMT numbers in P. puparum. Overall, our study shows unusually high rates of mitochondrial evolution and considerable variation in NUMT accumulation in Chalcidoidea.

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.

Comparative Mt Genomics of the Tipuloidea (Diptera: Nematocera: Tipulomorpha) and Its Implications for the Phylogeny of the Tipulomorpha

A traditionally controversial taxon, the Tipulomorpha has been frequently discussed with respect to both its familial composition and relationships with other Nematocera. The interpretation of internal relationships within the Tipuloidea, which include the Tipulidae sensu stricto, Cylindrotomidae, Pediciidae and Limoniidae, is also problematic. We sequenced the first complete mitochondrial (mt) genome of Symplecta hybrida (Meigen, 1804), which belongs to the subfamily Chioneinae of family Limoniidae, and another five nearly complete mt genomes from the Tipuloidea. We did a comparative analysis of these mt genomics and used them, along with some other representatives of the Nematocera to construct phylogenetic trees. Trees inferred by Bayesian methods strongly support a sister-group relationship between Trichoceridae and Tipuloidea. Tipulomorpha are not supported as the earliest branch of the Diptera. Furthermore, phylogenetic trees indicate that the family Limoniidae is a paraphyletic group.

Mitochondrial genomes of two phlebotomine sand flies, Phlebotomus chinensis and Phlebotomus papatasi (Diptera: Nematocera), the first representatives from the family Psychodidae

BACKGROUND

Leishmaniasis is a worldwide but neglected disease of humans and animal transmitted by sand flies, vectors that also transmit other important diseases. Mitochondrial genomes contain abundant information for population genetic and phylogenetic studies, important in disease management. However, the available mitochondrial sequences of these crucial vectors are limited, emphasizing the need for developing more mitochondrial genetic markers.

METHODS

The complete mitochondrial genome of Phlebotomus chinensis was amplified in eight fragments and sequenced using primer walking. The mitochondrial genome of Phlebotomus papatasi was reconstructed from whole-genome sequencing data available on Genbank. The phylogenetic relationship of 24 selected representatives of Diptera was deduced from codon positions 1 and 2 for 13 protein coding genes, using Bayesian inference (BI) and maximum likelihood (ML) methods.

RESULTS

We provide the first Phlebotomus (P. chinensis and P. papatasi) mitochondrial genomes. Both genomes contain 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and an A + T-rich region. The gene order of Phlebotomus mitochondrial genomes is identical with the ancestral gene order of insect. Phylogenetic analyses demonstrated that Psychodidae and Tanyderidae are sister taxa. Potential markers for population genetic study of Phlebotomus species were also revealed.

CONCLUSION

The generated mitochondrial genomes of P. chinensis and P. papatasi represent a useful resource for comparative genomic studies and provide valuable future markers for the population genetic study of these important Leishmania vectors. Our results also preliminary demonstrate the phylogenetic placement of Psychodidae based on their mitochondrial genomes.

A resourceful genome: updating the functional repertoire and evolutionary role of animal mitochondrial DNAs

Recent data from mitochondrial genomics and proteomics research demonstrate the existence of several atypical mitochondrial protein-coding genes (other than the standard set of 13) and the involvement of mtDNA-encoded proteins in functions other than energy production in several animal species including humans. These results are of considerable importance for evolutionary and cellular biology because they indicate that animal mtDNAs have a larger functional repertoire than previously believed. This review summarizes recent studies on animal species with a non-standard mitochondrial functional repertoire and discusses how these genetic novelties represent promising candidates for studying the role of the mitochondrial genome in speciation.

Complete mitochondrial genomes of Ceratobaeus sp. and Idris sp. (Hymenoptera: Scelionidae): shared gene rearrangements as potential phylogenetic markers at the tribal level

We sequenced the complete mitochondrial genomes of two sceliond taxa (Ceratobaeus sp. and Idris sp.). An atypical tRNA-Arg which lacks a D-stem was identified in both taxa, and represents a potentially derived character of sceliond wasps. A number of tRNA genes have rearranged in the two mitochondrial genomes compared with the ancestral organization. Some of these derived genome organizations are shared, and thus have much potential as phylogenetic markers at the tribal level in the subfamily Scelioninae. We test the influence of third codon inclusion/exclusion, alignment methods and partition schemes on the reconstruction of phylogenetic relationships. The results show that inclusion of third codon positions does not appear to be problematic when investigating the phylogeny of closely related taxa. Muscle and PartitionFinder schemes significantly improve the likelihood scores.

Insect mitochondrial genomics: implications for evolution and phylogeny

The mitochondrial (mt) genome is, to date, the most extensively studied genomic system in insects, outnumbering nuclear genomes tenfold and representing all orders versus very few. Phylogenomic analysis methods have been tested extensively, identifying compositional bias and rate variation, both within and between lineages, as the principal issues confronting accurate analyses. Major studies at both inter- and intraordinal levels have contributed to our understanding of phylogenetic relationships within many groups. Genome rearrangements are an additional data type for defining relationships, with rearrangement synapomorphies identified across multiple orders and at many different taxonomic levels. Hymenoptera and Psocodea have greatly elevated rates of rearrangement offering both opportunities and pitfalls for identifying rearrangement synapomorphies in each group. Finally, insects are model systems for studying aberrant mt genomes, including truncated tRNAs and multichromosomal genomes. Greater integration of nuclear and mt genomic studies is necessary to further our understanding of insect genomic evolution.

The complete mitochondrial genomes of two ghost moths, Thitarodes renzhiensis and Thitarodes yunnanensis: the ancestral gene arrangement in Lepidoptera

Background

Lepidoptera encompasses more than 160,000 described species that have been classified into 45–48 superfamilies. The previously determined Lepidoptera mitochondrial genomes (mitogenomes) are limited to six superfamilies of the lineage Ditrysia. Compared with the ancestral insect gene order, these mitogenomes all contain a tRNA rearrangement. To gain new insights into Lepidoptera mitogenome evolution, we sequenced the mitogenomes of two ghost moths that belong to the non-ditrysian lineage Hepialoidea and conducted a comparative mitogenomic analysis across Lepidoptera.

Results

The mitogenomes of Thitarodes renzhiensis and T. yunnanensis are 16,173 bp and 15,816 bp long with an A + T content of 81.28 % and 82.34 %, respectively. Both mitogenomes include 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and the A + T-rich region. Different tandem repeats in the A + T-rich region mainly account for the size difference between the two mitogenomes. All the protein-coding genes start with typical mitochondrial initiation codons, except for cox1 (CGA) and nad1 (TTG) in both mitogenomes. The anticodon of trnS(AGN) in T. renzhiensis and T. yunnanensis is UCU instead of the mostly used GCU in other sequenced Lepidoptera mitogenomes. The 1,584-bp sequence from rrnS to nad2 was also determined for an unspecified ghost moth (Thitarodes sp.), which has no repetitive sequence in the A + T-rich region. All three Thitarodes species possess the ancestral gene order with trnI-trnQ-trnM located between the A + T-rich region and nad2, which is different from the gene order trnM-trnI-trnQ in all previously sequenced Lepidoptera species. The formerly identified conserved elements of Lepidoptera mitogenomes (i.e. the motif ‘ATAGA’ and poly-T stretch in the A + T-rich region and the long intergenic spacer upstream of nad2) are absent in the Thitarodes mitogenomes.

Conclusion

The mitogenomes of T. renzhiensis and T. yunnanensis exhibit unusual features compared with the previously determined Lepidoptera mitogenomes. Their ancestral gene order indicates that the tRNA rearrangement event(s) likely occurred after Hepialoidea diverged from other lepidopteran lineages. Characterization of the two ghost moth mitogenomes has enriched our knowledge of Lepidoptera mitogenomes and contributed to our understanding of the mechanisms underlying mitogenome evolution, especially gene rearrangements.

Novel gene rearrangements in the mitochondrial genome of a webspinner, Aposthonia japonica (Insecta: Embioptera)

Webspinners (order Embioptera) are polyneopteran insects characterized by enlarged foretarsi with silk glands, whose silk is used to produce galleries in which the insects live gregariously. The phylogenetic position of webspinners has been debated. In the present study, an almost complete mitochondrial DNA (mtDNA) sequence of Embioptera is reported for the first time. The mtDNA of a webspinner, Aposthonia japonica , has the 13 protein-coding genes (PCGs) generally found in metazoan mtDNA sequences. There is a translocation of a large region including atp6, atp8, cox3, nad3, and nad5 as well as a duplication of the 12S rRNA gene. The rearrangement does not seem to affect nucleotide composition, although amino acid composition in some parts of the mtDNA is biased compared with other Polyneoptera species. Based on phylogenetic analyses using nucleotide sequences of all PCGs concatenated with two rRNA genes and the amino acid sequences of all PCGs, A. japonica is sister to Verophasmatodea, a suborder of typical stick and leaf insects.

The first mitochondrial genome for the wasp superfamily Platygastroidea: the egg parasitoid Trissolcus basalis

The nearly complete mitochondrial (mt) genome of an egg parasitoid, Trissolcus basalis (Wollaston), was sequenced using both 454 and Illumina next-generation sequencing technologies. A portion of the noncoding region remained unsequenced, possibly owing to the presence of repeats. The sequenced portion of the genome is 15,768 bp and has a high A+T content (84.2%), as is typical for hymenopteran mt genomes. A total of 36 of the 37 genes normally present in animal mt genomes were located. The one exception was trnR; a truncated version of this gene is present between trnS(1) and nd5, but it is unclear whether this gene fragment could code for the entire trnR gene. The mt gene arrangement of T. basalis is different from other Proctotrupomorpha mt genomes, with a number of trn genes in different positions. However, no shared derived gene rearrangements were identified in the present study. Bayesian analyses of mt genomes from 29 hymenopteran taxa and seven other orders of holometabolous insects support some uncontroversial evolutionary relationships, but indicate that much higher levels of taxonomic sampling are necessary for the resolution of family and superfamily relationships.

Mitochondrial genome sequences of Nematocera (lower Diptera): evidence of rearrangement following a complete genome duplication in a winter crane fly

The complete mitochondrial DNA sequences of eight representatives of lower Diptera, suborder Nematocera, along with nearly complete sequences from two other species, are presented. These taxa represent eight families not previously represented by complete mitochondrial DNA sequences. Most of the sequences retain the ancestral dipteran mitochondrial gene arrangement, while one sequence, that of the midge Arachnocampa flava (family Keroplatidae), has an inversion of the trnE gene. The most unusual result is the extensive rearrangement of the mitochondrial genome of a winter crane fly, Paracladura trichoptera (family Trichocera). The pattern of rearrangement indicates that the mechanism of rearrangement involved a tandem duplication of the entire mitochondrial genome, followed by random and nonrandom loss of one copy of each gene. Another winter crane fly retains the ancestral diperan gene arrangement. A preliminary mitochondrial phylogeny of the Diptera is also presented.