Characterization of four mitochondrial genomes of Crambidae (Lepidoptera,Pyraloidea) and phylogenetic implications

Loxostege turbidalis, Loxostege aeruginalis, Pyrausta despicata, and Crambus perlellus belong to Crambidae, Pyraloidea. Their mitochondrial genomes (mitogenomes) were successfully sequenced. The mitogenomes of L. turbidalis, L. aeruginalis, P. despicata, and C. perlellus are 15 240 bp, 15 339 bp, 15 389 bp, and 15 440 bp. The four mitogenomes all have a typical insect mitochondrial gene order, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and one A + T rich region (control region). The PCGs are initiated by the typical ATN codons, except CGA for the cox1 gene. Most PCGs terminate with common codon TAA or TAG, the incomplete codon T is found as the stop codon for cox2, nad4, and nad5. Most tRNA genes exhibit typical cloverleaf structure, except trnS1 (AGN) lacking the dihydrouridine (DHU) arm. The secondary structure of rRNA of four mitogenomes were predicted. Poly-T structure and micro-satellite regions are conserved in control regions. The phylogenetic analyses based on 13 PCGs showed the relationships of subfamilies in Pyraloidea. Pyralidae, and Crambidae are monophyletic, respectively. Pyralidae comprises four subfamilies, which form the following topology with high support values: (Galleriinae + ((Pyralinae + Epipaschiinae)+ Phycitinae)). Crambidae includes seven subfamilies and is divided into two lineages. Pyraustinae and Spilomelinae are sister groups of each other, and form the “PS clade.” Other five subfamilies (Crambinae, Acentropinae, Scopariinae, Schoenobiinae, and Glaphyriinae) form the “non-PS clade” in the Bayesian inference tree. However, Schoenobiinae is not grouped with the other four subfamilies and located at the base of Crambidae in two maximum likelihood trees.     

The mitochondrial genome of the butterfly Papilio xuthus (Lepidoptera: Papilionidae) and related phylogenetic analyses

The nearly complete mitochondrial genome of the butterfly Papilio xuthus (Lepidoptera: Papilionidae) was sequenced for its nucleotide sequence of 13,964 bp. The genome has a typical gene order identical to other lepidopteran species. All tRNAs showed same stable canonical clover-leaf structure as those of other insects, except for tRNA^Ser (AGN), in which the dihydrouracil arm (DHU arm) could not form stable stem–loop structure. Anomalous initiation codons have been observed for the cox1 gene, where the ATTACG hexa-nucleotide was believed to be involved in the initiation signaling. Twelve mitochondrial protein-coding gene sequence data were used to infer the phylogenetic relationships among the insect orders. Even though the number of insect orders represented by complete mitochondrial genomes is still limited, several well-established relationships are evident in the phylogenetic analysis of the complete sequences. Monophyly of the Homometabola was not supported in this paper. Phylogenetic analyses of the available species of Bombycoidea, Pyraloidea, Papilionoidea and Tortricidea bolstered the current morphology-based hypothesis that Bombycoidea, Pyraloidea and Papilionoidea are monophyletic (Obtectomera). Bombycoidea (Bombyx mandarina and Antheraea pernyi) and Papilionoidea (P. xuthus and Coreana raphaelis) formed a sister group.     

Description of new mitochondrial genomes (Spodoptera litura, Noctuoidea and Cnaphalocrocis medinalis, Pyraloidea) and phylogenetic reconstruction of Lepidoptera with the comment on optimization schemes

We newly sequenced mitochondrial genomes of Spodoptera litura and Cnaphalocrocis medinalis belonging to Lepidoptera to obtain further insight into mitochondrial genome evolution in this group and investigated the influence of optimal strategies on phylogenetic reconstruction of Lepidoptera. Estimation of p-distances of each mitochondrial gene for available taxonomic levels has shown the highest value in ND6, whereas the lowest values in COI and COII at the nucleotide level, suggesting different utility of each gene for different hierarchical group when individual genes are utilized for phylogenetic analysis. Phylogenetic analyses mainly yielded the relationships (((((Bombycoidea + Geometroidea) + Noctuoidea) + Pyraloidea) + Papilionoidea) + Tortricoidea), evidencing the polyphyly of Macrolepidoptera. The Noctuoidea concordantly recovered the familial relationships (((Arctiidae + Lymantriidae) + Noctuidae) + Notodontidae). The tests of optimality strategies, such as exclusion of third codon positions, inclusion of rRNA and tRNA genes, data partitioning, RY recoding approach, and recoding nucleotides into amino acids suggested that the majority of the strategies did not substantially alter phylogenetic topologies or nodal supports, except for the sister relationship between Lycaenidae and Pieridae only in the amino acid dataset, which was in contrast to the sister relationship between Lycaenidae and Nymphalidae in Papilionoidea in the remaining datasets.     

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.     

Complete mitochondrial genome of the Eurasian flying squirrel Pteromys volans (Sciuromorpha, Sciuridae) and revision of rodent phylogeny

In this study, the complete mitochondrial genome of the Eurasian flying squirrel Pteromys volans (Rodentia, Sciuromorpha, Sciuridae) was sequenced and characterized in detail. The entire mitochondrial genome of P. volans consisted of 16,513 bp and contained 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and two non-coding regions. Its gene arrangement pattern was consistent with the mammalian ground pattern. The overall base composition and AT contents were similar to those of other rodent mitochondrial genomes. The light-strand origin generally identified between tRNA ^Asn and tRNA ^Cys consisted of a secondary structure with an 11-bp stem and an 11-bp loop. The large control region was constructed of three characteristic domains, ETAS, CD, and CSB without any repeat sequences. Each domain contained ETAS1, subsequences A, B, and C, and CSB1, respectively. In order to examine phylogenetic contentious issues of the monophyly of rodents and phylogenetic relationships among five rodent suborders, here, phylogenetic analyses based on nucleotide sequence data of the 35 rodent and 3 lagomorph mitochondrial genomes were performed using the Bayesian inference and maximum likelihood method. The result strongly supported the rodent monophyly with high node confidence values (BP 100 % in ML and BPP 1.00 in BI) and also monophylies of four rodent suborders (BP 85–100 % in ML and BPP 1.00 in BI), except for Anomalumorpha in which only one species was examined here. Also, phylogenetic relationships among the five rodent suborders were suggested and discussed in detail.     

The complete mitochondrial genome of the leafminer <Emphasis Type="Italic">Liriomyza trifolii</Emphasis> (Diptera: Agromyzidae)

Liriomyza trifolii (Diptera: Agromyzidae) is one of the most economically significant pests in the world. In this paper we present sequence data for the complete mitochondrial genome of L. trifolii. The circular genome is 16,141 bp long and contains one encoding region including 37 genes and one non-coding A+T-rich region. Gene numbers and organization is similar to that of the typical insect mitochondrial genomes except that two additional tRNA genes are found in the A+T-rich region (tRNA^Thr and tRNA^Leu(UUR)). All of the protein initiation codons are ATN, except ND1 which begins with GTG and COI which is initiated by the quadruplet ATCA. The 22 tRNA anticodons of L. trifolii match those observed in Drosophila yakuba, and all of tRNAs form the typical cloverleaf structure except for tRNA^Ser(AGN), which has lost the DHU-arm. The A+T-rich region of L. trifolii also contains two previously noted Diperan features—a highly conserved polyT stretch and a (TA)_n stretch.     

The phylogeny of leaf beetles (Chrysomelidae) inferred from mitochondrial genomes

The high-level classification of Chrysomelidae (leaf beetles) currently recognizes 12 or 13 well-established subfamilies, but the phylogenetic relationships among them remain ambiguous. Full mitochondrial genomes were newly generated for 27 taxa and combined with existing GenBank data to provide a dataset of 108 mitochondrial genomes covering all subfamilies. Phylogenetic analysis under maximum likelihood and Bayesian inference recovered the monophyly of all subfamilies, except that Timarcha was split from Chrysomelinae in some analyses. Three previously recognized major clades of Chrysomelidae were broadly supported: the ‘chrysomeline’ clade consisting of (Chrysomelinae (Galerucinae + Alticinae)); the ‘sagrine’ clade with internal relationships of ((Bruchinae + Sagrinae) + (Criocerinae + Donaciinae)), and the ‘eumolpine’ clade comprising (Spilopyrinae (Cassidinae (Eumolpinae (Cryptocephalinae + Lamprosomatinae)))). Relationships among these clades differed between data treatments and phylogenetic algorithms, and were complicated by two additional deep lineages, Timarcha and Synetinae. Various topological tests favoured the PhyloBayes software as the preferred inference method, resulting in the arrangement of (chrysomelines (eumolpines + sagrines)), with Timarcha placed as sister to the chrysomeline clade and Synetinae as a deep lineage splitting near the base. Whereas mitogenomes provide a solid framework for the phylogeny of Chrysomelidae, the basal relationships do not agree with the topology of existing molecular studies and remain one of the most difficult problems of Chrysomelidae phylogenetics.     

The complete mitochondrial genome of the yellow coaster, <Emphasis Type="Italic">Acraea</Emphasis> <Emphasis Type="Italic">issoria</Emphasis> (Lepidoptera: Nymphalidae: Heliconiinae: Acraeini): sequence,gene organization and a unique tRNA translocation event

In this paper, the complete mitochondrial genome of Acraea issoria (Lepidoptera: Nymphalidae: Heliconiinae: Acraeini) is reported; a circular molecule of 15,245 bp in size. For A. issoria, genes are arranged in the same order and orientation as the complete sequenced mitochondrial genomes of the other lepidopteran species, except for the presence of an extra copy of tRNA^Ile(AUR)b in the control region. All protein-coding genes of A. issoria mitogenome start with a typical ATN codon and terminate in the common stop codon TAA, except that COI gene uses TTG as its initial codon and terminates in a single T residue. All tRNA genes possess the typical clover leaf secondary structure except for tRNA^Ser(AGN), which has a simple loop with the absence of the DHU stem. The sequence, organization and other features including nucleotide composition and codon usage of this mitochondrial genome were also reported and compared with those of other sequenced lepidopterans mitochondrial genomes. There are some short microsatellite-like repeat regions (e.g., (TA)₉, polyA and polyT) scattered in the control region, however, the conspicuous macro-repeats units commonly found in other insect species are absent.     

Mitotype-specific sequences related to cytoplasmic male sterility in Oryza species

Plant mitochondrial genomes contain a large number of mitotype-specific sequences (MSS) which establish a mitochondrial genome structure distinct from other mitotypes. In rice, nine mitochondrial genomes have been sequenced, which provides us with the possibility of characterizing the MSS of rice and probing their relationship to cytoplasmic male sterility (CMS) in rice. We therefore analyzed the mitochondrial genomes of CW-CMS, LD-CMS, WA-CMS, N and Nipponbare lines, and found 57 MSS with sizes ranging from 102 to 5,745 bp, and with an aggregate length of 92.4 kb. The MSS account for more than 14.5 % of the rice mitochondrial genome and are a significant contributing factor in the variation of mitochondrial genome sizes. Of the MSS tested, 34 MSS exhibited polymorphism among rice lines, and 14 MSS were further confirmed as being specific to CMS. This includes nine MSS specific to sporophytic CMS, three specific to gametophytic CMS, and two shared by all types of CMS. Interestingly, except for CMS genes orf(H)79 and orf352 which are partly or fully overlapping with some MSS fragments, there are ten more open reading frames of unknown function that were detected in CMS-specific MSS, hinting at their possible roles in plant CMS. These novel findings provide us with potential new molecular tools to direct the breeding of CMS lines in hybrid rice breeding programs.     

The complete mitochondrial genome of the butterfly <Emphasis Type="Italic">Apatura metis</Emphasis> (Lepidoptera: Nymphalidae)

As an important pest in the Slender Leaved Willow (Salix alba), Apatura metis is called Freyer’s purple emperor, and its mitochondrial genome is 15,236 bp long. The encoded genes for 22 tRNA genes, two ribosomal RNA (rrnL and rrnS) genes, and 13 protein-coding genes (PCGs), and a control region in the A. metis mitochondria are highly homologous to other lepidopteran species. The mitochondrial genome of A. metis is biased toward a high A + T content (A + T = 80.5%). All protein-coding genes, except for COI begins with the CGA codon as observed in other lepidopterans, start with a typical ATN initiation codon. All tRNAs show the classic clover-leaf structure, except that the dihydrouridine (DHU) arm of tRNA ^Ser(AGN) forms a simple loop. The A. metis A + T-rich region contains some conserved structures including a structure combining the motif ‘ATAGA’ and 19 bp poly (T) stretch, which is similar to those found in other lepidopteran mitogenomes. The phylogenetic analyses of lepidopterans based on mitogenomes sequences demonstrate that each of the six superfamilies is monophyletic, and the relationship among them is (((Noctuoidea + (Geometroidea + Bombycoidea)) + Pyraloidea) + Papilionoidea) + Tortricoidea. In Papilionoidea group, our conclusion argues that ((Lycaenidae + Pieridae) + Nymphalidae) + Papilionidae.     

Complete mitochondrial genomes of three Oxya grasshoppers (Orthoptera) and their implications for phylogenetic reconstruction

Oxya is a genus of grasshoppers (Orthoptera: Acridoidea) attacking rice and other gramineous plants in Africa and Asia. In the present study, we characterized complete mitochondrial genomes (mitogenomes) of three species, Oxya japonica japonica (15,427 bp), Oxya hainanensis (15,443 bp) and Oxya agavisa robusta (15,552 bp) collected from China. The three mitogenomes contained a typical gene set of metazoan mitogenomes and shared the same gene order with other Acridid grasshoppers, including the rearrangement of tRNA^Asp and tRNA^Lys. Analyses of pairwise genetic distances showed that ATP8 was the least conserved gene, while COI the most conserved. To determine the position of Oxya grasshoppers in the phylogeny of Acrididae, we reconstructed phylogenetic trees among 64 species from across 11 subfamilies using nucleotide sequences of mitogenomes. While the tree confirms traditional classifications of Acrididae at major higher-levels, it suggests a few modifications for classifications at lower-levels.     

The complete mitochondrial genomes of two brown algae (Laminariales, Phaeophyceae) and phylogenetic analysis within Laminaria

Two complete mitochondrial genome sequences for Laminaria longissima (=Saccharina longissima) and Laminaria hyperborea are reported in this study. They had circular mapping organization with slight difference in size (37,628 and 37,976 bp, respectively) and contained almost the same set of mitochondrial genes, including the genes for three rRNAs (23S, 16S, and 5S), 25 tRNAs, 35 known mitochondrial proteins, and three to four Open Reading Frame genes (ORFs). Both mitochondrial genomes exhibited typical gene content and organization of Laminaria mtDNAs except for the existence of ORF157 genes being located between rRNA large subunit gene 5 (rpl5) and ORF129-139 in L. hyperborea as found in that of Laminaria digitata. The phylogenetic analysis based on mitochondrial genomes supported the hypothesis of the split of the genus Laminaria, and the result of this study provided important information on the molecular evolution of Laminaria.     

Higher‐level phylogeny of longhorn beetles (Coleoptera: Chrysomeloidea) inferred from mitochondrial genomes

Cerambycidae (longhorn beetles) and related families in the superfamily Chrysomeloidea are important components of forest ecosystems and play a key role in nutrient cycling and pollination. Using full mitochondrial genomes and dense taxon sampling, the phylogeny of Chrysomeloidea with a focus on Cerambycidae and allied families was explored. We used 151 mitochondrial genomes (75 newly sequenced) covering all families and 29 subfamilies of Chrysomeloidea. Our results reveal that (i) Chrysomelidae (leaf beetles) are sister to all other chrysomeloid families; (ii) Cerambycidae sensu stricto (s. s.) is polyphyletic due to the inclusion of other families that split Cerambycidae into a ‘lamiine’ clade comprising Lepturinae sensu lato (s. l.) + (Lamiinae + Spondylidinae) and a ‘cerambycine’ clade comprising Dorcasominae + (Cerambycinae + Prioninae s. l.); (iii) the subfamilies within the two clades of Cerambycidae s. s. were monophyletic, except for the placement of Necydalinae nested in Lepturinae, and the placement of Parandrinae within Prioninae (now considered as tribes Necydalini and Parandrini, respectively); (iv) smaller families were grouped into two major clades: one composed of Disteniidae+Vesperidae and the other composed of Orsodacnidae + (Megalopodidae + Oxypeltidae); (v) relationships among the four major clades were poorly supported but were resolved as ((cerambycines + (Disteniidae + Vesperidae) + Orsodacnidae + (Megalopodidae + Oxypeltidae)) + lamiines. Divergence time analyses estimated that Chrysomeloidea originated ca. 154.1 Mya during the late Jurassic, and most subfamilies of Cerambycidae originated much earlier than subfamilies of Chrysomelidae. The diversification of families within Chrysomeloidea was largely coincident with the radiation of angiosperms during the Early Cretaceous.     

Complete sequence of the mitochondrial genome of Odontamblyopus rubicundus (Perciformes: Gobiidae): genome characterization and phylogenetic analysis

Odontamblyopus rubicundus is a species of gobiid fishes, inhabits muddy-bottomed coastal waters. In this paper, the first complete mitochondrial genome sequence of O. rubicundus is reported. The complete mitochondrial genome sequence is 17119 bp in length and contains 13 protein-coding genes, two rRNA genes, 22 tRNA genes, a control region and an L-strand origin as in other teleosts. Most mitochondrial genes are encoded on H-strand except for ND6 and seven tRNA genes. Some overlaps occur in protein-coding genes and tRNAs ranging from 1 to 7 bp. The possibly nonfunctional L-strand origin folded into a typical stem-loop secondary structure and a conserved motif (5^′-GCCGG-3^′) was found at the base of the stem within the tRNA ^Cys gene. The TAS, CSB-2 and CSB-3 could be detected in the control region. However, in contrast to most of other fishes, the central conserved sequence block domain and the CSB-1 could not be recognized in O. rubicundus, which is consistent with Acanthogobius hasta (Gobiidae). In addition, phylogenetic analyses based on different sequences of species of Gobiidae and different methods showed that the classification of O. rubicundus into Odontamblyopus due to morphology is debatable.     

The chloroplast and mitochondrial genomes of two Gomphonema parvulum (Bacillariophyta) environmental isolates from South Carolina (United States) and Virginia (United States)

Gomphonema parvulum is a cosmopolitan freshwater diatom that is used as an indicator in water quality biomonitoring. In this study, we report the culturing of two geographically separated isolates from southeastern North America, their morphology, and the sequencing and assembly of their mitochondrial and chloroplast genomes. Morphologically, both strains fit G. parvulum sensu lato, but the frustules from a protected habitat in South Carolina were smaller than those cited in the historic data of this species from the same location as well as a second culture from Virginia. Phylogenetic analyses using the rbcL gene placed both within a clade with G. parvulum. Genetic markers, including full chloroplast and mitochondrial genomes and the nuclear small subunit rRNA gene region were assembled from each isolate. The organellar genomes of the two strains varied slightly in size due to small differences in intergenic regions with chloroplast genomes of 121,035 bp and 121,482 bp and mitochondrial genomes of 34,639 bp and 34,654 bp. The intraspecific pairwise identities of the chloroplast and mitochondrial genomes of these two isolates were 97.9% and 95.4%, respectively. Multigene phylogenetic analysis demonstrated a close relationship between G. parvulum, Gomphoneis minuta, and Didymosphenia geminata.     

The complete mitogenome of Parapoynx crisonalis (Walker, 1859) (Lepidoptera: Crambidae), with phylogenetic relationships amongst three acentropine larval forms

The larvae of Parapoynx crisonalis (Walker, 1859) is an important pest in China, causing significant damage on the aquatic crops. To improve the management of P. crisonalis and demonstrate the phylogenetic relationships amongst Pyraloidea, the complete mitochondrial genome (mitogenome) sequences of four aquatic pyralid moths, which belong to three acentropine larval forms – Parapoynx-type, Nymphula-type and Potamomusa-type, and twenty-three terricolous pyralids were compared. The complete mitogenome P. crisonalis, a typical parapoynx larval species, was 15,374 base pairs in length and included 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and an A + T-rich region. The phylogenetic analysis revealed that all the aquatic acentropine moths formed a monophyletic clade with strong phylogenetic support via maximum likelihood and Bayesian methods with a bootstrap value of 100% and a posterior probability of 1.00. The Potamomusa-type was more closely related to the Nymphula-type than to the Parapoynx-type.     

The complete mitochondrial genome of the oriental fruit moth <Emphasis Type="Italic">Grapholita molesta</Emphasis> (Busck) (Lepidoptera: Tortricidae)

The oriental fruit moth, Grapholita molesta (Busck) (Lepidoptera: Tortricidae) currently is one of the economically most destructive pest species of stone and pome fruits worldwide. Here we sequenced the complete mitochondrial genome of this pest. This genome is 15,776 bp long, with an A + T content of 81.24%, containing 37 typical animal mitochondrial genes and an A + T-rich region. All gene are arranged as hypothesized ancestral gene order of insects except for trnM, which was shuffled from 3′ downstream of trnQ to 5′ upstream of trnI. cox1 gene uses unusual CGA start codon, as that in all other sequenced lepidopteran mitochondrial genome. The secondary structures for the two rRNA genes were predicted. All helices typically present in insect mitochondrial rRNA genes are generated. A microsatellite sequence was inserted into the region of H2347 in rrnL in G. molesta and two other sequenced tortricid mitochondrial genomes, indicating that the insertion event in this helix might occurred anciently in family Tortricidae. All of the 22 typical animal tRNA genes have a typical cloverleaf structure except for trnS2, in which the D-stem pairings in the DHU arm are absent. An intergenic sequence is present between trnQ and nad2 as well as in other sequenced lepidopteran mitochondrial genomes, which was presumed to be a remnant of trnM gene and its boundary sequences after the duplication of trnM to the upstream of trnI in Lepidoptera. The A + T-rich region is 836 bp, containing six repeat sequences of “TTATTATTATTATTAAATA(G)TTT.”     

Complete sequences of maternally inherited mitochondrial genomes in musselsUnio pictorum (Bivalvia, Unionidae)

Mitochondrial genomes are frequently used to infer phylogenetic relationships. Some taxa are, however, poorly represented. To facilitate better understanding of the potential of mitochondrial genome data in freshwater mussels, we present here, for the first time, the mitochondrial sequences of 4 complete F-type mitochondrial genomes from the European freshwater bivalveUnio pictorum (Unionidae). These genomes are very compact (15 761 bp) but have a typical gene complement for bilaterian mitochondrial genomes and a very similar organization to other unionid genomes available in databases. Very low nucleotide diversity within the species suggests a small effective population size of PolishU. pictorum, a phenomenon of potential importance for environmental management policies.     

A molecular phylogeny for the pyraloid moths (Lepidoptera: Pyraloidea) and its implications for higher‐level classification

Pyraloidea, one of the largest superfamilies of Lepidoptera, comprise more than 15 684 described species worldwide, including important pests, biological control agents and experimental models. Understanding of pyraloid phylogeny, the basis for a predictive classification, is currently provisional. We present the most detailed molecular estimate of relationships to date across the subfamilies of Pyraloidea, and assess its concordance with previous morphology‐based hypotheses. We sequenced up to five nuclear genes, totalling 6633 bp, in each of 42 pyraloids spanning both families and 18 of the 21 subfamilies, plus up to 14 additional genes, for a total of 14 826 bp, in 21 of those pyraloids plus all 24 outgroups. Maximum likelihood analyses yield trees that, within Pyraloidea, differ little among datasets and character treatments and are strongly supported at all levels of divergence (83% of nodes with bootstrap ≥80%). Subfamily relationships within Pyralidae, all very strongly supported (>90% bootstrap), differ only slightly from a previous morphological analysis, and can be summarized as Galleriinae + Chrysauginae (Phycitinae (Pyralinae + Epipaschiinae)). The main remaining uncertainty involves Chrysauginae, of which the poorly studied Australian genera may constitute the basal elements of Galleriinae + Chrysauginae or even of Pyralidae. In Crambidae the molecular phylogeny is also strongly supported, but conflicts with most previous hypotheses. Among the newly proposed groupings are a ‘wet‐habitat clade’ comprising Acentropinae + Schoenobiinae + Midilinae, and a provisional ‘mustard oil clade’ containing Glaphyriinae, Evergestinae and Noordinae, in which the majority of described larvae feed on Brassicales. Within this clade a previous synonymy of Dichogaminae with the Glaphyriinae is supported. Evergestinae syn. n. and Noordinae syn. n. are here newly synonymized with Glaphyriinae, which appear to be paraphyletic with respect to both. Pyraustinae and Spilomelinae as sampled here are each monophyletic but form a sister group pair. Wurthiinae n. syn. , comprising the single genus Niphopyralis Hampson, which lives in ant nests, are closely related to, apparently subordinate within, and here newly synonymized with, Spilomelinae syn. n.     

Characterization of mitochondrial genome of Chinese wild mulberry silkworm, <Emphasis Type="Italic">Bomyx mandarina</Emphasis> (Lepidoptera: Bombycidae)

The complete mitochondrial genome of Chinese Bombyx mandarina (ChBm) was determined. The circular genome is 15682 bp long, and contains a typical gene complement, order, and arrangement identical to that of Bombyx mori (B. mori) and Japanese Bombyx mandarina (JaBm) except for two additional tRNA-like structures: tRNA ^Ser(TGA)-like and tRN ^AIle(TAT)-like. All protein-coding sequences are initiated with a typical ATN codon except for the COI gene, which has a 4-bp TTAG putative initiator codon. Eleven of 13 protein-coding genes (PCGs) have a complete termination codon (all TAA), but the remaining two genes terminate with incomplete codons. All tRNAs have the typical clover-leaf structures of mitochondrial tRNAs, with the exception of tRNA ^Ser(TGA)-like, with a four stem-and-loop structure. The length of the A+T-rich region of ChBm is 484 bp, shorter than those of JaBm (747 bp) and B. mori (494–499 bp). Phylogenetic analysis among B. mori, ChBm, JaBm, and Antheraea pernyi (Anpe) showed that B. mori is more closely related to ChBm than JaBm. The earliest divergence time estimate for B. mori-ChBm and B. mori-JaBm is about 1.08±0.18–1.41±0.24 and 1.53±0.20–2.01±0.26 Mya, respectively. ChBm and JaBm diverged around 1.11±0.16–1.45±0.21 Mya.