Complete mitochondrial genomes of five skippers (Lepidoptera: Hesperiidae) and phylogenetic reconstruction of Lepidoptera

We sequenced mitogenomes of five skippers (family Hesperiidae, Lepidoptera) to obtain further insight into the characteristics of butterfly mitogenomes and performed phylogenetic reconstruction using all available gene sequences (PCGs, rRNAs, and tRNAs) from 85 species (20 families in eight superfamilies). The general genomic features found in the butterflies also were found in the five skippers: a high A + T composition (79.3%–80.9%), dominant usage of TAA stop codon, similar skewness pattern in both strands, consistently length intergenic spacer sequence between tRNA^Gln and ND2 (64–87 bp), conserved ATACTAA motif between tRNA^{Ser (UCN)} and ND1, and characteristic features of the A + T-rich region (the ATAGA motif, varying length of poly-T stretch, and poly-A stretch). The start codon for COI was CGA in four skippers as typical, but Lobocla bifasciatus evidently possessed canonical ATG as start codon. All species had the ancestral arrangement tRNA^Asn/tRNA^{Ser (AGN)}, instead of the rearrangement tRNA^{Ser (AGN)}/tRNA^Asn, found in another skipper species (Erynnis). Phylogenetic analyses using all available genes (PCGs, rRNAS, and tRNAs) yielded the consensus superfamilial relationships ((((((Bombycoidea + Noctuoidea + Geometroidea) + Pyraloidea) + Papilionoidea) + Tortricoidea) + Yponomeutoidea) + Hepialoidea), confirming the validity of Macroheterocera (Bombycoidea, Noctuoidea, and Geometroidea in this study) and its sister relationship to Pyraloidea. Within Rhopalocera (butterflies and skippers) the familial relationships (Papilionidae + (Hesperiidae + (Pieridae + ((Lycaenidae + Riodinidae) + Nymphalidae)))) were strongly supported in all analyses (0.98–1 by BI and 96–100 by ML methods), rendering invalid the superfamily status for Hesperioidea. On the other hand, current mitogenome-based phylogeny did not find consistent superfamilial relationships among Noctuoidea, Geometroidea, and Bombycoidea and the familial relationships within Bombycoidea between analyses, requiring further taxon sampling in future studies.     

Complete mitochondrial genome of the stonefly Cryptoperla stilifera Sivec (Plecoptera: Peltoperlidae) and the phylogeny of Polyneopteran insects

We present the complete mitogenome of a stonefly, Cryptoperla stilifera Sivec (Plecoptera; Peltoperlidae). The mitogenome was a circular molecule consisting of 15,633 nucleotides, 37 genes and a A + T-rich region. C. stilifera mitogenome was similar to Pteronarcys princeps mitogenome (Plecoptera; Pteronarcyidae). All transfer RNA genes (tRNAs) had typical cloverleaf secondary structures except for trnSer (AGN), where the stem-loop structure of the dihydrouridine (DHU) arm was missing. The A + T-rich region of C. stilifera had two stem-loops and each had two interlink. Three conserved sequence blocks (CSBs) were present in the A + T-rich regions of C. stilifera, Peltoperla tarteri and Peltoperla arcuata. Moreover, many polynucleotide stretches (Poly N, N = A, T and C) in the A + T-rich region of C. stilifera Phylogenetic relationships of Polyneopteran species were constructed based on the nucleotide sequences of 13 protein coding genes (PCGs). Both maximum likelihood (ML) and Bayesian inference (BI) analyses supported Grylloblattodea as the sister group to Plecoptera + Dermaptera and Embiidina and Phasmatodea as sister groups.     

The complete mitochondrial genome of the sycamore lace bug Corythucha ciliata (Hemiptera: Tingidae)

The complete mitochondrial genome of the sycamore lace bug, Corythucha ciliata, was sequenced in this study. It represents the first sequenced mitogenome of family Tingidae in Heteroptera. The mitogenome of C. ciliata is 15,257 bp and contains 37 genes including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes and a large non-coding region. Gene arrangement, nucleotide content, codon usage, and amino acid composition and asymmetry indicate a high degree of conservation with six other species of Cimicomorpha. The 13 PCGs initiated with ATN as the start codon and terminated with TAA, TA or T as stop codon. The evolutionary rate of each PCG was different, among which ATP8 showed the highest rate while ATP6 indicated the lowest rate. The 22 tRNAs genes apparently fold into a typical cloverleaf structure; however, the anticodon (TTC) of trnSer (AGN) differs from other Heteropteran insects. Secondary structure modeling of rRNA genes revealed similarity to other insects, except for two incomplete helices (H1648 and H2735) in lrRNA. The predicted secondary structure of lrRNA indicates 45 helices in six domains, whereas srRNA has 27 helices in three domains. Three potential stem–loops and two tandem repeats (–TCTAAT–) were identified in the A+T-rich region. Phylogenetic analysis indicated that C. ciliata is a sister group to other Heteroptera species based on analysis of the 13 PCGs.     

The complete mitochondrial genome of Sasakia funebris (Leech) (Lepidoptera: Nymphalidae) and comparison with other Apaturinae insects

Sasakia funebris, a member of the lepidopteran family, Nymphalidae (superfamily Papilionoidea) is a rare species and is found only in some areas of South China. In this study, the 15,233 bp long complete mitochondrial genome of S. funebris was determined, and harbors the gene arrangement identical to all other sequenced lepidopteran insects. The nucleotide composition of the genome is highly A + T biased, accounting for 81.2%. All protein-coding genes (PCGs) start with typical ATN codons, except for COI which begins with the CGA codon. All tRNAs have a typical clover-leaf secondary structure, except for tRNA^Ser(AGN), the dihydrouridine (DHU) arm of which forms a simple loop. The S. funebris A + T-rich region of 370 bp contains several features common to the Lepidoptera insects, including the motif ATAGA followed by a 19 bp poly-T stretch, and two tandem repeats consisting of 18 bp repeat units and 14 bp repeat units. The phylogenetic analyses of Apaturinae based on mitogenome sequences showed: (S. funebris + Sasakia charonda) + (Apatura metis + Apatura ilia). This result is consistent with the morphological classification.     

Complete mitochondrial genome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae) and the phylogenetic relationship of Saturniidae species

Mitochondrial genome (mitogenome) can provide information for genomic structure as well as for phylogenetic analysis and evolutionary biology. In this study, we present the complete mitogenome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae), a well-known silk-producing and ornamental insect with the largest wing surface area of all moths. The mitogenome of A. atlas is a circular molecule of 15,282 bp long, and its nucleotide composition shows heavily biased towards As and Ts, accounting for 79.30%. This genome comprises 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A + T-rich region. It is of note that this genome exhibits a slightly positive AT skew, which is different from the other known Saturniidae species. All PCGs are initiated by ATN codons, except for COI with CGA instead. Only six PCGs use a common stop codon of TAA or TAG, whereas the remaining seven use an incomplete termination codon T or TA. All tRNAs have the typical clover-leaf structure, with an exception for tRNA^Ser(AGN). The A. atlas A + T-rich region contains non-repetitive sequences, but harbors several features common to the Bombycoidea insects. The phylogenetic relationships based on Maximum Likelihood method provide a well-supported outline of Saturniidae, which is in accordance with the traditional morphological classification and recent molecular works.     

Characterization of the complete mitochondrial genome of Tryporyza incertulas, in comparison with seven other Pyraloidea moths

The complete 15,223-bp mitochondrial genome (mitogenome) of Tryporyza incertulas (Walker) (Lepidoptera: Pyraloidea: Crambidae) was determined, characterized and compared with seven other species of superfamily Pyraloidea. The order of 37 genes was typical of insect mitochondrial DNA sequences described to date. Compared with other moths of Pyraloidea, the A + T biased (77.0%) of T. incertulas was the lowest. Eleven protein-coding genes (PCGs) utilized the standard ATN, but cox1 used CGA and nad4 used AAT as the initiation codons. Ten protein-coding genes had the common stop codon TAA, except nad3 having TAG as the stop codon, and cox2, nad4 using T, TA as the incomplete stop codons, respectively. All of the tRNA genes had typical cloverleaf secondary structures except trnS1(AGN), in which the dihydrouridine (DHU) arm did not form a stable stem-loop structure. There was a spacer between trnQ and nad2, which was common in Lepidoptera moths. A 6-bp motif ‘ATACTA’ between trnS2(UCN) and nad1, a 7-bp motif “AGC(T)CTTA” between trnW and trnC and a 6-bp motif “ATGATA” of overlapping region between atp8 and atp6 were found in Pyraloidea moths. The A + T-rich region contained an ‘ATAGT(A)’-like motif followed by a poly-T stretch. In addition, two potential stem-loop structures, a duplicated 19-bp repeat element, and two microsatellites ‘(TA)₁₂’ and ‘(TA)₉’ were observed in the A + T-rich region of T. incertulas mitogenome. Finally, the phylogenetic relationships of Pyraloidea species were constructed based on amino acid sequences of 13 PCGs of mitogenomes using Bayesian inference (BI) and maximum likelihood (ML) methods. These molecular-based phylogenies supported the morphological classification on relationships within Pyraloidea species.     

The complete mitogenome of Apocheima cinerarius (Lepidoptera: Geometridae: Ennominae) and comparison with that of other lepidopteran insects

The complete mitochondrial genome (mitogenome) of a female flightless geometrid moth Apocheima cinerarius was found to be 15,722 bp in length, containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a control region. The A + T content of the complete mitogenome is 80.83%. The AT skew value ([A − T] / [A + T]) is 0.027. The 13 PCGs of the mitogenome start with typical ATN codons, except for cox1 with the start codon CGA. All the tRNA genes have typical cloverleaf secondary structures, except for trnSer(AGN). The secondary structures of rrnL and rrnS were predicted. Six structural domains including conserved regions (IV, V) and variable regions (I, II, III, VI) were identified in the secondary structure of rrnL. The secondary structure of rrnS consists of 3 structural domains. The control region of A. cinerarius begins with conserved motifs of “ATAGA” + 19-bp poly T. It also contains a microsatellite-like (TA)₂₆, a stem-and-loop structure, and a poly-A stretch. Phylogenetic analysis showed that Geometroidea is more closely related to Bombycoidea than to Noctuoidea. A. cinerarius is more closely related to Biston panterinaria than to Phthonandria atrilineata, which is in accordance with the conventional morphology-based classification.     

A unique tRNA gene family and a novel,highly expressed ORF in the mitochondrial genome of the silver-lip pearl oyster, Pinctada maxima (Bivalvia: Pteriidae)

Characteristics of mitochondrial (mt) DNA such as gene content and arrangement, as well as mt tRNA secondary structure, are frequently used in comparative genomic analyses because they provide valuable phylogenetic information. However, most analyses do not characterize the relationship of tRNA genes from the same mt genome and, in some cases, analyses overlook possible novel open reading frames (ORFs) when the 13 expected protein-coding genes are already annotated. In this study, we describe the sequence and characterization of the complete mt genome of the silver-lip pearl oyster, Pinctada maxima. The 16,994-bp mt genome contains the same 13 protein-coding genes (PCGs) and two ribosomal RNA genes typical of metazoans. The gene arrangement, however, is completely distinct from that of all other available bivalve mt genomes, and a unique tRNA gene family is observed in this genome. The unique tRNA gene family includes two trnS^− AGY and trnQ genes, a trnM isomerism, but it lacks trnS^− CUN. We also report the first clear evidence of alloacceptor tRNA gene recruitment (trnP → trnS^− AGY) in mollusks. In addition, a novel ORF (orfUR1) expressed at high levels is present in the mt genome of this pearl oyster. This gene contains a conserved domain, “Oxidored_q1_N”, which is a member of Complex I and thus may play an important role in key biological functions. Because orfUR1 has a very similar nucleotide composition and codon bias to that of other genes in this genome, we hypothesize that this gene may have been moved to the mt genome via gene transfer from the nuclear genome at an early stage of speciation of P. maxima, or it may have evolved as a result of gene duplication, followed by rapid sequence divergence. Lastly, a 319-bp region was identified as the possible control region (CR) even though it does not correspond to the longest non-coding region in the genome. Unlike other studies of mt genomes, this study compares the evolutionary patterns of all available bivalve mt tRNA and atp8 genes.     

Complete nucleotide sequence and gene organization of the mitochondrial genome of Paa spinosa (Anura: Ranoidae)

The mt genome of Paa spinosa (Anura: Ranoidae) is a circular molecule of 18,012 bp in length, containing 38 genes (including an extra copy of tRNA-Met gene). This mt genome is characterized by three distinctive features: a cluster of rearranged tRNA genes (LTPF tRNA gene cluster), a tandem duplication of tRNA-Met gene (Met1 and Met2), and distinct repeat regions at both 5′ and 3′-sides in the control region. Comparing the locations and the sequences of all tRNA-Met genes among Ranoidae, and constructing NJ tree of the nucleotide of those tRNA-Met genes, we suggested a tandem duplication of tRNA-Met gene can be regarded as a synapomorphy of Dicroglossinae. To further investigate the phylogenetic relationships of anurans, phylogenetic analyses (BI, ML and MP) based on the nucleotide dataset and the corresponding amino acid dataset of 11 protein-coding genes (except ND5 and ATP8) arrived at the similar topology.     

Complete mitochondrial genome of Malaysian Mahseer (Tor tambroides)

This is the first documentation of the complete mitochondrial genome sequence of the Malaysian Mahseer, Tor tambroides. The 16,690 bp mitogenome with GenBank accession number JX444718 contains 13 protein genes, 22 tRNAs, two rRNAs, and a noncoding control region (D-loop) as is typical of most vertebrates. The phylogenomic reconstruction of this newly generated data with 21 Cypriniformes GenBank accession ID concurs with the recognized status of T. tambroides within the subfamily Cyprininae. This is in agreement with previous hypotheses based on morphological and partial mitochondrial analyses.     

The complete mitochondrial genome of the Ailanthus silkmoth, Samia cynthia cynthia (Lepidoptera: Saturniidae)

The complete mitochondrial genome (mitogenome) of the Ailanthus silkmoth, Samia cynthia cynthia (Lepidoptera: Saturniidae) was determined. The circular genome is 15,345 bp long, and presents a typical gene organization and order for sequenced mitogenomes of Bombycidea species. The nucleotide composition of the genome is highly A+T biased, accounting for 79.86%. The AT skew of the genome is slightly negative, indicating the occurrence of more Ts than As, as found in other Saturniidae species. All protein-coding genes (PCGs) are initiated by ATN codons, except for COI and COII, which are tentatively designated by CGA and GTG, respectively, as observed in other insects. Four of 13 PCGs, including COI, COII, ATP6, and ND3, harbor the incomplete termination codons, T or TA. With an exception for tRNA^Ser(AGN), all other tRNAs can form a typical clover-leaf structure of mitochondrial tRNA. The 359 bp A+T-rich region of S. c. cynthia contains non-repetitive sequences, but harbors several features common to the Bombycidea insects, including the motif ATAGA followed by a poly-T stretch of 19 bp, a microsatellite-like (AT)₇ element preceded by the ATTTA motif, and a poly-A element upstream tRNA^Met. The phylogenetic analyses support the morphology-based current hypothesis that Bombycidae and Saturniidae are monophyletic. Our result confirms that Saturniini and Attacini form a reciprocal monophyletic group within Saturniidae.     

The complete mitochondrial genome of the Chinese hook snout carp Opsariichthys bidens (Actinopterygii: Cypriniformes) and an alternative pattern of mitogenomic evolution in vertebrate

The complete mitochondrial genome sequence of the Chinese hook snout carp, Opsariichthys bidens, was newly determined using the long and accurate polymerase chain reaction method. The 16,611-nucleotide mitogenome contains 13 protein-coding genes, two rRNA genes (12S, 16S), 22 tRNA genes, and a noncoding control region. We use these data and homologous sequence data from multiple other ostariophysan fishes in a phylogenetic evaluation to test hypothesis pertaining to codon usage pattern of O. bidens mitochondrial protein genes as well as to re-examine the ostariophysan phylogeny. The mitochondrial genome of O. bidens reveals an alternative pattern of vertebrate mitochondrial evolution. For the mitochondrial protein genes of O. bidens, the most frequently used codon generally ends with either A or C, with C preferred over A for most fourfold degenerate codon families; the relative synonymous codon usage of G-ending codons is greatly elevated in all categories. The codon usage pattern of O. bidens mitochondrial protein genes is remarkably different from the general pattern found previously in the relatively closely related zebrafish and most other vertebrate mitochondria. Nucleotide bias at third codon positions is the main cause of codon bias in the mitochondrial protein genes of O. bidens, as it is biased particularly in favor of C over A. Bayesian analysis of 12 concatenated mitochondrial protein sequences for O. bidens and 46 other teleostean taxa supports the monophyly of Cypriniformes and Otophysi and results in a robust estimate of the otophysan phylogeny.     

The complete mitochondrial genome sequence and gene organization of the mud crab (Scylla paramamosain) with phylogenetic consideration

The complete mitochondrial genome is of great importance for better understanding the genome-level characteristics and phylogenetic relationships among related species. In the present study, we determined the complete mitochondrial genome DNA sequence of the mud crab (Scylla paramamosain) by 454 deep sequencing and Sanger sequencing approaches. The complete genome DNA was 15,824 bp in length and contained a typical set of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a putative control region (CR). Of 37 genes, twenty-three were encoded by the heavy strand (H-strand), while the other ones were encoded by light strand (L-strand). The gene order in the mitochondrial genome was largely identical to those obtained in most arthropods, although the relative position of gene tRNA^His differed from other arthropods. Among 13 protein-coding genes, three (ATPase subunit 6 (ATP6), NADH dehydrogenase subunits 1 (ND1) and ND3) started with a rare start codon ATT, whereas, one gene cytochrome c oxidase subunit I (COI) ended with the incomplete stop codon TA. All 22 tRNAs could fold into a typical clover-leaf secondary structure, with the gene sizes ranging from 63 to 73 bp. The phylogenetic analysis based on 12 concatenated protein-coding genes showed that the molecular genetic relationship of 19 species of 11 genera was identical to the traditional taxonomy.     

The complete mitochondrial genome of Macrobrachium nipponense

The complete mitochondrial (mt) genome sequence plays an important role in the accurate determination of phylogenetic relationships among metazoans. Herein, we determined the complete mt genome sequence, structure and organization of Macrobrachium nipponense (M. nipponense) (GenBank ID: NC_015073.1) and compared it to that of Macrobrachium lanchesteri (M. lanchesteri) and Macrobrachium rosenbergii (M. rosenbergii). The 15,806 base pair (bp) M. nipponense mt genome, which is comprised of 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs) and 2 ribosomal RNAs (rRNAs), is slightly larger than that of M. lanchesteri (15,694 bp, GenBank ID: NC_012217.1) and M. rosenbergii (15,772 bp, GenBank ID: NC_006880.1). The M. nipponense genome contains a high AT content (66.0%), which is a common feature among metazoan mt genomes. Compared with M. lanchesteri and M. rosenbergii, we found a peculiar non-coding region of 950 bp with a microsatellite-like (TA)₆ element and many hairpin structures. The 13 PCGs are comprised of a total of 3707 codons, excluding incomplete termination codons, and the most frequently used amino acid is Leu (16.0%). The predicted start codons in the M. nipponense mt genome include ATG, ATC and ATA. Seven PCGs use TAA as a stop codon, whereas two use TAG, three use T and only one uses TA. Twenty-three of the genes are encoded on the L strand, and ND1, ND4, ND5, ND4L, 12S rRNA, 16S rRNA, tRNA^His, tRNA^Pro, tRNA^Phe, tRNA^Val, tRNA^Gln, tRNA^Cys, tRNA^Tyr and a tRNA^Leu are encoded on the H strand. The two rRNAs of M. nipponense and M. rosenbergii are encoded on the H strand, whereas the M. lanchesteri rRNAs are encoded on the L stand.     

Characterization of the complete mitochondrial genome of the giant silkworm moth, Eriogyna pyretorum (Lepidoptera: Saturniidae)

The complete mitochondrial genome (mitogenome) of Eriogyna pyretorum (Lepidoptera: Saturniidae) was determined as being composed of 15,327 base pairs (bp), including 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The arrangement of the PCGs is the same as that found in the other sequenced lepidopteran. The AT skewness for the E. pyretorum mitogenome is slightly negative (-0.031), indicating the occurrence of more Ts than As. The nucleotide composition of the E. pyretorum mitogenome is also biased toward A + T nucleotides (80.82%). All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 and 2 (cox1 and cox2). Two of the 13 PCGs harbor the incomplete termination codon by T. All tRNA genes have a typical clover-leaf structure of mitochondrial tRNA, with the exception of trnS1(AGN) and trnS2(UCN). Phylogenetic analysis among the available lepidopteran species supports the current morphology-based hypothesis that Bombycoidea, Geometroidea, Notodontidea, Papilionoidea and Pyraloidea are monophyletic. As has been previously suggested, Bombycidae (Bombyx mori and Bombyx mandarina), Sphingoidae (Manduca sexta) and Saturniidae (Antheraea pernyi, Antheraea yamamai, E. pyretorum and Caligula boisduvalii) formed a group.     

The complete mitochondrial genomes of three grasshoppers, Asiotmethis zacharjini, Filchnerella helanshanensis and Pseudotmethis rubimarginis (Orthoptera: Pamphagidae)

The complete mitogenomes of Asiotmethis zacharjini, Filchnerella helanshanensis and Pseudotmethis rubimarginis are 15,660 bp, 15,657 bp and 15,661 bp in size, respectively. All three mitogenomes contain a standard set of 13 protein - coding genes, 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and an A + T-rich region in the same order as those of the other analysed caeliferan species, including the rearrangement of trnAsp and trnLys. The putative initiation codon for the cox1 gene in the three species is CCG. The long polythymine stretch (T-stretch) in the A + T-rich region of the three species is not adjacent to the trnIle but inside the stem–loop sequence in the majority strand. The mitogenomes of F. helanshanensis and P. rubimarginis have higher overall similarities. The characterization of the three mitogenomes will enrich our knowledge on the Pamphagidae mitogenome. The phylogenetic analyses indicated that within the Caelifera, Pyrgomorphoidea is a sister group to Acridoidea. The species from the Pamphagidae form a monophyletic group, as is the case for Acrididae. Furthermore, the two families cluster as sister groups, supporting the monophyly of Acridoidea. The relationships among eight acridid subfamilies were (Cyrtacanthacridinae + (Calliptaminae + (Catantopinae + (Oxyinae + (Melanopline + (Acridinae + (Oedipodinae + Gomphocerinae))))))).     

The Complete Mitochondrial Genome of the Beet Webworm,Spoladea recurvalis (Lepidoptera: Crambidae) and Its Phylogenetic Implications

The complete mitochondrial genome (mitogenome) of the beet webworm, Spoladea recurvalis has been sequenced. The circular genome is 15,273 bp in size, encoding 13 protein-coding genes (PCGs), two rRNA genes, and 22 tRNA genes and containing a control region with gene order and orientation identical to that of other ditrysian lepidopteran mitogenomes. The nucleotide composition of the mitogenome shows a high A+T content of 80.9%, and the AT skewness is slightly negative (-0.023). All PCGs start with the typical ATN codons, except for COX1, which may start with the CGA codon. Nine of 13 PCGs have the common stop codon TAA; however, COX1, COX2 and ND5 utilize the T nucleotide and ND4 utilizes TA nucleotides as incomplete termination codons. All tRNAs genes are folded into the typical cloverleaf structure of mitochondrial tRNAs, except for the tRNA^Ser(AGY) gene, in which the DHU arm fails to form a stable stem-loop structure. A total of 157 bp intergenic spacers are scattered in 17 regions. The overlapping sequences are 42 bp in total and found in eight different locations. The 329 bp AT-rich region is comprised of non-repetitive sequences, including the motif ATAG, which is followed by a 14 bp poly-T stretch, a (AT₁₁ microsatellite-like repeat, which is adjacent to the motif ATTTA, and a 9 bp poly-A, which is immediately upstream from the tRNA^Met gene. Phylogenetic analysis, based on 13 PCGs and 13 PCGs+2 rRNAs using Bayesian inference and Maximum likelihood methods, show that the classification position of Pyraloidea is inconsistent with the traditional classification. Hesperioidea is placed within the Papilionoidea rather than as a sister group to it. The Pyraloidea is placed within the Macrolepidoptera with other superfamilies instead of the Papilionoidea.     

The complete mitochondrial genome of the leafminer Liriomyza sativae (Diptera: Agromyzidae): great difference in the A+T-rich region compared to Liriomyza trifolii

The complete mitochondrial genome sequence of Liriomyza sativae Blanchard (15,551 bp) was determined and analyzed in this study. The circular genome contained 37 genes including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and an A + T-rich region. The initiation codons of COI and ND1 were ‘ATCA’ and ‘GTG’, respectively. ND2 gene used the truncated termination codon ‘T’. All the tRNA genes had the typical cloverleaf secondary structures except for tRNA^Ser(AGN) gene, which was found with the absence of a DHU arm. In addition, a tRNA-like secondary structure (tRNA^Met) was found in the A + T-rich region. The great difference was that the length of L. sativae A + T-rich region was 597 bp shorter than that of Liriomyza trifolii (Burgess). Meanwhile, some minor differences such as ‘TATA’ block were also observed in L. sativae in contrast to ‘TACA’ block in L. trifolii. There were also some essential structure elements such as ‘TATA’ block, ‘G(A)_nT’ block, poly-T stretch and stem-and-loop structure in the A + T-rich region of L. sativae mitochondrial genome.     

Comparative analyses of the complete mitochondrial genomes of the two ruminant hookworms Bunostomum trigonocephalum and Bunostomum phlebotomum

Bunostomum trigonocephalum and Bunostomum phlebotomum are blood-feeding hookworms of sheep and cattle, causing considerable economic losses to the live stock industries. Studying genetic variability within and among hookworm populations is critical to addressing epidemiological and ecological questions. Mitochondrial (mt) DNA is known to provide useful markers for investigations of population genetics of hookworms, but mt genome sequence data are scant. In the present study, the complete mitochondrial DNA (mtDNA) sequences of the sheep and goat hookworm B. trigonocephalum were determined for the first time, and the mt genome of B. phlebotomum from yak in China was also sequenced for comparative analyses of their gene contents and genome organizations. The lengths of mt DNA sequences of B. trigonocephalum sheep isolate, B.trigonocephalum goat isolate and B. phlebotomum China yak isolate were 13,764 bp, 13,771 bp and 13,803 bp in size, respectively. The identity of the mt genomes was 99.7% between B. trigonocephalum sheep isolate and B. trigonocephalum goat isolate. The identity of B. phlebotomum China yak isolate mt genomes was 85.3% with B. trigonocephalum sheep isolate, and 85.2% with B. trigonocephalum goat isolate. All the mt genes of the two hookworms were transcribed in the same direction and gene arrangements were consistent with those of the GA3 type, including 12 protein-coding genes, 2 rRNA genes and 22 tRNA genes, but lacking ATP synthetase subunit 8 gene. The mt genomes of B. trigonocephalum and B. phlebotomum were similar to prefer bases A and T, the contents of A + T are 76.5% (sheep isolate), 76.4% (goat isolate) and 76.9% (China yak isolate), respectively. Phylogenetic relationships reconstructed using concatenated amino acid sequences of 12 protein-coding genes with three methods (maximum likelihood, Bayesian inference and neighbor joining) revealed that the B. trigonocephalum and B. phlebotomum represent distinct but closely-related species. These data provide novel and useful genetic markers for studying the systematics, and population genetics of the two ruminant hookworms.     

林氏按蚊线粒体全基因组序列的测定及基于线粒体基因组的按蚊属系统发育分析(英文)

【目的】对林氏按蚊Anopheles lindesayi完整的线粒体基因组进行测序及分析,依据已知的线粒体基因组构建并讨论按蚊属蚊虫的分子系统发育关系。【方法】对林氏按蚊线粒体基因组进行测序、注释,并对其基本特征和基本组成进行分析。基于串联的13个蛋白质编码基因的核苷酸序列和氨基酸序列,用ML法和贝叶斯法构建林氏按蚊和按蚊属其他32种蚊虫的系统发育树,据此探讨按蚊属蚊虫的系统发育关系和系统分类。【结果】林氏按蚊线粒体基因组全长为15 366 bp,包含13个蛋白质编码基因,22个tRNA基因,2个rRNA基因和一段控制区。林氏按蚊线粒体基因组呈现明显的AT偏斜和GC偏斜,AT偏斜为正,GC偏斜为负。除了COX1使用TCG和ND5使用GTG作为起始密码子以外,其他蛋白质编码基因的起始密码子均遵循ATN原则;终止密码子为TAA或者T。除了tRNA^Ser(AGN)以外,其他的tRNA基因均呈现典型的三叶草二级结构。控制区AT含量最高,为94.54%。滑窗分析显示蛋白质编码基因是用于构建亚属或属水平系统发育关系的最佳分子标记。系统发育树强烈支持塞蚊亚属Cellia、按蚊亚属Anopheles、徕蚊亚属Nyssorhynchus和柯特蚊亚属Kerteszia均为单系群。小五斑按蚊An. atroparvus和四斑按蚊An. quadrimaculatus A这两个种聚到一起,从传统的形态分类上讲,它们和林氏按蚊均属于按蚊亚属按蚊系蚊虫。但本研究构建的4个系统发育树均显示,(小五斑按蚊An. atroparvus+四斑按蚊An. quadrimaculatus A)和林氏按蚊被属于迈蚊系的中华按蚊分开,这为两个系的分类提供了新的论点。【结论】本研究获得了林氏按蚊的完整的线粒体基因组,探析了按蚊属的线粒体基因组特征和系统发育关系,为进一步研究蚊科线粒体基因组和系统发育关系提供了依据。