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.     

小长蝽线粒体基因组全序列测定与长蝽总科系统发育分析

为了解小长蝽Nysius ericae(Schilling)线粒体基因组结构及长蝽总科的分子系统发育关系。本试验采用Illumina MiSeq测序平台对小长蝽线粒体基因进行测序,对基因组序列进行拼装、注释和特征分析,利用最大似然法和贝叶斯法构建基于12种长蝽总科昆虫线粒体全基因组核苷酸序列的系统发育树。小长蝽线粒体基因组全长为16 330 bp(GenBank登录号：MW465654),基因组包括13个蛋白编码基因(PCGs),22个tRNA基因,2个rRNA基因和1段非编码控制区。11个蛋白质编码基因的起始密码子为典型的ATN;cox1,nad4l的起始密码子为TTG。cob的终止密码子为TAG,其余蛋白编码基因的终止密码子为TAA。只有trnS1缺少DHU臂,其余tRNA基因均能形成典型的三叶草结构。12种长蝽总科昆虫线粒体全基因组序列构建的昆虫系统发育树结果显示,小长蝽与Nysius plebeius具有更近的亲缘关系,且与传统形态学分类基本一致。小长蝽线粒体基因组符合长蝽总科线粒体基因组的一般特征。结果表明小长蝽与N.plebeius的亲缘关系更近。     

The complete mitochondrial genome of the nudibranch Roboastra europaea (Mollusca: Gastropoda) supports the monophyly of opisthobranchs

The complete nucleotide sequence (14,472 bp) of the mitochondrial genome of the nudibranch Roboastra europaea (Gastropoda: Opisthobranchia) was determined. This highly compact mitochondrial genome is nearly identical in gene organization to that found in opisthobranchs and pulmonates (Euthyneura) but not to that in prosobranchs (a paraphyletic group including the most basal lineages of gastropods). The newly determined mitochondrial genome differs only in the relative position of the trnC gene when compared with the mitochondrial genome of Pupa strigosa, the only opisthobranch mitochondrial genome sequenced so far. Pupa and Roboastra represent the most basal and derived lineages of opisthobranchs, respectively, and their mitochondrial genomes are more similar in sequence when compared with those of pulmonates. All phylogenetic analyses (maximum parsimony, minimum evolution, maximum likelihood, and Bayesian) based on the deduced amino acid sequences of all mitochondrial protein-coding genes supported the monophyly of opisthobranchs. These results are in agreement with the classical view that recognizes Opisthobranchia as a natural group and contradict recent phylogenetic studies of the group based on shorter sequence data sets. The monophyly of opisthobranchs was further confirmed when a fragment of 2,500 nucleotides including the mitochondrial cox1, rrnL, nad6, and nad5 genes was analyzed in several species representing five different orders of opisthobranchs with all common methods of phylogenetic inference. Within opisthobranchs, the polyphyly of cephalaspideans and the monophyly of nudibranchs were recovered. The evolution of mitochondrial tRNA rearrangements was analyzed using the cox1+rrnL+nad6+nad5 gene phylogeny. The relative position of the trnP gene between the trnA and nad6 genes was found to be a synapomorphy of opisthobranchs that supports their monophyly.     

The complete nucleotide sequence and gene organization of the mitochondrial genome of the oriental mole cricket, Gryllotalpa orientalis (Orthoptera: Gryllotalpidae)

The complete nucleotide sequences of the mitochondrial genome (mitogenome) of the oriental mole cricket, Gryllotalpa orientalis (Orthoptera: Gryllotalpidae), were determined. The 15,521-bp-long G. orientalis mitogenome contains typical gene complement, base composition, and codon usage found in metazoan mitogenomes. The G. orientalis mitogenome contains the third lowest A+T content (70.5%) among the complete insects mt genome sequences. The initiation codon for the G. orientalis COI gene appears to be ATG, instead of the tetranucleotides, which have been postulated to act as initiation codon for Locusta migratoria and some lepidopteran COI genes. The initiation codon for ND2 appears to be GTG, which is rare, but has been designated as an initiator of Tricholepidion gertschi ND2. All anticodons of G. orientalis tRNAs were identical to Drosophila yakuba and L. migratoria. The tRNA(Ser)(AGN) could not form a stable stem loop structure in the DHU arm as shown in many other insect tRNA(Ser)(AGN). Phylogenetic analysis of nucleotide sequence information from all mt genes supported a monophyletic Diptera, a monophyletic Lepidoptera, a monophyletic Coleoptera, a monophyletic Mecopterida (Diptera+Lepidoptera), and a monophyletic Endopterygota (Diptera+Lepidoptera+Coleoptera), suggesting that the complete insect mitogenome sequence has a resolving power to the diversification events within Endopterygota. However, the relationships of ancient insect orders were unstable, indicating the limited use of mitogenome information at deeper phylogenetic depth.     

The mitochondrial genome of the thermal dimorphic fungus Penicillium marneffei is more closely related to those of molds than yeasts

We report the complete sequence of the mitochondrial genome of Penicillium marneffei, the first complete mitochondrial DNA sequence of a thermal dimorphic fungus. This 35 kb mitochondrial genome contains the genes encoding ATP synthase subunits 6, 8, and 9 (atp6, atp8, and atp9), cytochrome oxidase subunits I, II, and III (cox1, cox2, and cox3), apocytochrome b (cob), reduced nicotinamide adenine dinucleotide ubiquinone oxireductase subunits (nad1, nad2, nad3, nad4, nad4L, nad5, and nad6), ribosomal protein of the small ribosomal subunit (rps), 28 tRNAs, and small and large ribosomal RNAs. Analysis of gene contents, gene orders, and gene sequences revealed that the mitochondrial genome of P. marneffei is more closely related to those of molds than yeasts.     

Species‐specific mitochondrial gene rearrangements in biting midges and vector species identification

Abstract Partial mitochondrial gene sequences of 16 Culicoides species were determined to elucidate phylogenetic relations among species and to develop a molecular identification method for important virus vector species. In addition, the analysis found mitochondrial gene rearrangement in several species. Sequences of the mitochondrial genome region, cox1–trnL2–cox2 (1940–3785 bp) of 16 Culicoides and additional sequences were determined in some species, including whole mitochondrial genome sequences of Culicoides arakawae. Nine species showed common organization in this region, with three genes cox1–trnL2–cox2 and a small or no intergenic region (0–30 bp) between them. The other seven species showed translocation of tRNA and protein‐coding genes and/or insertion of AT‐rich non‐coding sequences (65–1846 bp) between the genes. The varied gene rearrangements among species within a genus is very rare for mitochondrial genome organization. Phylogenetic analyses based on the sequences of cox1+cox2 suggest a few clades among Japanese Culicoides species. No relationships between phylogenetic closeness and mitochondrial gene rearrangements were observed. Sequence data were used to establish a polymerase chain reaction tool to distinguish three important vector species from other Culicoides species, for which classification during larval stages is not advanced and identification is difficult.     

入侵害虫甘薯凹胫跳甲的鉴定及线粒体基因组分析

【目的】基于形态学鉴定和分子生物学技术确认甘薯凹胫跳甲Chaetocnema confinis是否入侵中国大陆,测定甘薯凹胫跳甲线粒体基因组序列,分析基因组结构及其系统发育关系。【方法】应用显微镜观察从广东不同地点采集的甘薯凹胫跳甲成虫的形态特征,并扩增cox1基因DNA序列进行分子鉴定;利用Illumina MiSeq测序平台对甘薯凹胫跳甲线粒体基因组进行测序、拼装、注释和特征分析;基于亲缘关系相近种属的线粒体基因组序列进行共线性分析和构建系统发育树,分析基因重排和系统发育关系。【结果】形态和分子鉴定结果表明大陆甘薯上发现的跳甲为甘薯凹胫跳甲。甘薯凹胫跳甲线粒体基因组序列大小为15 685 bp,包括有13个蛋白质编码基因、2个rRNA基因、22个tRNA基因和1个非编码控制区;这37个基因之间排列紧凑,间隔总长度101 bp,排列顺序与模式昆虫Drosophila yakuba线粒体基因排列顺序相同。甘薯凹胫跳甲线粒体基因组A+T含量为77.3%,具有明显的AT偏向性。13个蛋白质编码基因的起始密码子均为ATN。在22个tRNA基因中除trnS1的DHU臂缺失,trnD, trnG, trnN和trnT的二级结构中缺少TψC环外,其余17个都能形成典型的三叶草式二级结构,另trnK的反密码子突变为UUU,trnS1的反密码子突变为UCU。甘薯凹胫跳甲的控制区片段长度仅有60 bp,是目前已报道的昆虫线粒体基因组中最短的控制区。基于线粒体基因组的系统发育分析表明,甘薯凹胫跳甲与跳甲亚科(Alticinae)黄曲条跳甲Phyllotreta striolata亲缘关系最近。【结论】甘薯凹胫跳甲已经入侵到中国大陆。本研究获得了甘薯凹胫跳甲的线粒体基因组序列,为防控甘薯凹胫跳甲和分析叶甲科(Chrysomelidae)各种属间的系统发育关系奠定了基础。     

Evolutionary rates, divergence dates, and the performance of mitochondrial genes in Bayesian phylogenetic analysis

The mitochondrial genome is one of the most frequently used loci in phylogenetic and phylogeographic analyses, and it is becoming increasingly possible to sequence and analyze this genome in its entirety from diverse taxa. However, sequencing the entire genome is not always desirable or feasible. Which genes should be selected to best infer the evolutionary history of the mitochondria within a group of organisms, and what properties of a gene determine its phylogenetic performance? The current study addresses these questions in a Bayesian phylogenetic framework with reference to a phylogeny of plethodontid and related salamanders derived from 27 complete mitochondrial genomes; this topology is corroborated by nuclear DNA and morphological data. Evolutionary rates for each mitochondrial gene and divergence dates for all nodes in the plethodontid mitochondrial genome phylogeny were estimated in both Bayesian and maximum likelihood frameworks using multiple fossil calibrations, multiple data partitions, and a clock-independent approach. Bayesian analyses of individual genes were performed, and the resulting trees compared against the reference topology. Ordinal logistic regression analysis of molecular evolution rate, gene length, and the G-shape parameter a demonstrated that slower rate of evolution and longer gene length both increased the probability that a gene would perform well phylogenetically. Estimated rates of molecular evolution vary 84-fold among different mitochondrial genes and different salamander lineages, and mean rates among genes vary 15-fold. Despite having conserved amino acid sequences, cox1, cox2, cox3, and cob have the fastest mean rates of nucleotide substitution, and the greatest variation in rates, whereas rrnS and rrnL have the slowest rates. Reasons underlying this rate variation are discussed, as is the extensive rate variation in cox1 in light of its proposed role in DNA barcoding.     

绿眼赛茧蜂线粒体基因组全序列测定和分析

【目的】测定绿眼赛茧蜂Zele chlorophthalmus线粒体基因组全序列,分析其基因组结构及茧蜂科(Braconidae)部分类群的系统发育关系。【方法】利用Illumina MiSeq二代测序技术对绿眼赛茧蜂的线粒体基因组进行测序,对基因组序列进行拼装、注释,分析其结构特点和碱基组成;基于22种茧蜂科昆虫的COX1蛋白编码基因序列,应用最大似然法(ML)和邻接法(NJ)构建系统发育树,分析绿眼赛茧蜂与茧蜂科其他昆虫的系统发育关系。【结果】绿眼赛茧蜂线粒体基因组全长16 661 bp(GenBank登录号: MG822749),包含13个蛋白质编码基因、22个tRNA基因和2个rRNA基因,共37个基因,以及1个控制区。线粒体基因组有明显的核苷酸组成的偏倚,AT偏正,GC偏负,其A+T含量为82.83%。基因排列顺序与推测的昆虫祖先的序列不完全一致,tRNA基因7处发生重排。13个蛋白质编码基因均以ATN为起始密码子,以TAA为终止密码子。在22个tRNA基因二级结构中,除tRNA^His(H)缺失TΨC环和tRNA^Cys(C)仅剩二氢尿嘧啶(DHU)臂和反密码子臂外,其余tRNA基因均能形成典型的三叶草结构。基于COX1蛋白编码序列的系统发育分析结果显示,与绿眼赛茧蜂亲缘关系最近的是同属于赛茧蜂属的雪跗赛茧蜂Z. niveitarsis。【结论】本研究首次获得绿眼赛茧蜂线粒体基因组全序列。结果表明绿眼赛茧蜂隶属于优茧蜂亚科(Euphorinae)赛茧蜂属,并支持赛茧蜂属的单系性。     

The Complete Mitochondrial Genomes of Three Bristletails (Insecta: Archaeognatha): The Paraphyly of Machilidae and Insights into Archaeognathan Phylogeny

The order Archaeognatha was an ancient group of Hexapoda and was considered as the most primitive of living insects. Two extant families (Meinertellidae and Machilidae) consisted of approximately 500 species. This study determined 3 complete mitochondrial genomes and 2 nearly complete mitochondrial genome sequences of the bristletail. The size of the 5 mitochondrial genome sequences of bristletail were relatively modest, containing 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one control region. The gene orders were identical to that of Drosophila yakuba and most bristletail species suggesting a conserved genome evolution within the Archaeognatha. In order to estimate archaeognathan evolutionary relationships, phylogenetic analyses were conducted using concatenated nucleotide sequences of 13 protein-coding genes, with four different computational algorithms (NJ, MP, ML and BI). Based on the results, the monophyly of the family Machilidae was challenged by both datasets (W12 and G12 datasets). The relationships among archaeognathan subfamilies seemed to be tangled and the subfamily Machilinae was also believed to be a paraphyletic group in our study.     

The complete mitochondrial genome of the black coral Chrysopathes formosa (Cnidaria:Anthozoa:Antipatharia) supports classification of antipatharians within the subclass Hexacorallia

Black corals comprise a globally distributed shallow- and deep-water taxon whose phylogenetic position within the Anthozoa has been debated. We sequenced the complete mitochondrial genome of the antipatharian Chrysopathes formosa to further evaluate its phylogenetic relationships. The circular mitochondrial genome (18,398 bp) consists of 13 energy pathway protein-coding genes and two ribosomal RNAs, but only two transfer RNA genes (trnM and trnW), as well as a group I intron within the nad5 gene that contains the only copies of nad1 and nad3. No novel genes were found in the antipatharian mitochondrial genome. Gene order and genome content are most similar to those of the sea anemone Metridium senile (subclass Hexacorallia), with differences being the relative location of three contiguous genes (cox2-nad4-nad6) and absence (from the antipatharian) of a group I intron within the cox1 gene. Phylogenetic analyses of multiple protein-coding genes support classifying the Antipatharia within the subclass Hexacorallia and not the subclass Ceriantipatharia; however, the sister-taxon relationships of black corals within Hexacorallia remain inconclusive.     

Complete Mitochondrial Genome of the Free-Living Earwig, Challia fletcheri (Dermaptera: Pygidicranidae) and Phylogeny of Polyneoptera

The insect order Dermaptera, belonging to Polyneoptera, includes ～2,000 extant species, but no dermapteran mitochondrial genome has been sequenced. We sequenced the complete mitochondrial genome of the free-living earwig, Challia fletcheri, compared its genomic features to other available mitochondrial sequences from polyneopterous insects. In addition, the Dermaptera, together with the other known polyneopteran mitochondrial genome sequences (protein coding, ribosomal RNA, and transfer RNA genes), were employed to understand the phylogeny of Polyneoptera, one of the least resolved insect phylogenies, with emphasis on the placement of Dermaptera. The complete mitochondrial genome of C. fletcheri presents the following several unusual features: the longest size in insects is 20,456 bp; it harbors the largest tandem repeat units (TRU) among insects; it displays T- and G-skewness on the major strand and A- and C-skewness on the minor strand, which is a reversal of the general pattern found in most insect mitochondrial genomes, and it possesses a unique gene arrangement characterized by a series of gene translocations and/or inversions. The reversal pattern of skewness is explained in terms of inversion of replication origin. All phylogenetic analyses consistently placed Dermaptera as the sister to Plecoptera, leaving them as the most basal lineage of Polyneoptera or sister to Ephemeroptera, and placed Odonata consistently as the most basal lineage of the Pterygota.     

The complete mitochondrial genome of an ectoparasitic monopisthocotylean fluke Benedenia hoshinai (Monogenea: Platyhelminthes)

An exponential growth of mitochondrial genome information has brought significant progress in understanding the organismal phylogeny and mitochondrial genome evolution for many metazoans including platyhelminth groups. In this study, we determined the complete mitochondrial genome sequence for Benedenia hoshinai, an ectoparasitic monogenean species, and compared it with its congener Benedenia seriolae. The complete mitochondrial genome is 13,554 bp in length and contains 12 protein-coding genes (lacking the atp8 gene), 2 rRNA genes, and 22 tRNA genes, all encoded in the same direction as found in all other platyhelminth species sequenced to date. The gene arrangement of B. hoshinai mtDNA is almost identical to B. seriolae, differing only by the translocation of trnT between cox1 and rrnL. It is unclear whether the shared position of trnT between B. hoshinai and Gyrodactylus represents evidence for their phylogenetic affinity; testing this hypothesis requires further mitogenomic evidence.     

The mitochondrial genome of the fission yeast Schizosaccharomyces pombe. 2. Localization of genes by interspecific hybridization in strain ade7-50h- and cloning of the genome in small fragments

A series of 18 small overlapping restriction fragments has been cloned, covering the complete mitochondrial genome of Schizosaccharomyces pombe. By hybridizing mitochondrial gene probes from Saccharomyces cerevisiae and Neurospora crassa with restriction fragments of Schizosaccharomyces pombe mitochondrial DNA, the following homologous genes were localized on the mitochondrial genome of S. pombe: cob, cox1, cox2 and cox3, ATPase subunit 6 and 9 genes, the large rRNA gene and both types of open reading frames occurring in mitochondrial introns of various ascomycetes. The region of the genome, hybridizing with cob exon probes is separated by an intervening sequence of about 2500 bp, which is homologous with the first two introns of the cox1 gene in Saccharomyces cerevisiae (class II introns according to Michel et al. 1982). Similarly, in the cox1 homologous region, which covers about 4000 bp, two regions were detected hybridizing with class I intron probes, suggesting the existence of two cox1 introns in Schizosaccharomyces pombe. Hybridization with several specific exon probes with a determined order has revealed that cob, cox1, cox3 and the large rRNA gene are all transcribed from the same DNA strand. The low intensities of hybridization signals suggest a large evolutionary distance between Schizosaccharomyces pombe and Saccharomyces cerevisiae or Neurospora crassa mitochondrial genes. Considering the length of the mitochondrial DNA of Schizosaccharomyces pombe (about 19.4 kbp) and the expected length of the localized genes and intron sequences there is enough space left for encoding the expected set of tRNAs and the small rRNA gene. The existence of leader-, trailer-, ori- and spacer sequences or further unassigned reading frames is then restricted to a total length of about 3000 bp only.     

The complete mitochondrial genome and novel gene arrangement of the unique-headed bug Stenopirates sp. (Hemiptera: Enicocephalidae)

Many of true bugs are important insect pests to cultivated crops and some are important vectors of human diseases, but few cladistic analyses have addressed relationships among the seven infraorders of Heteroptera. The Enicocephalomorpha and Nepomorpha are consider the basal groups of Heteroptera, but the basal-most lineage remains unresolved. Here we report the mitochondrial genome of the unique-headed bug Stenopirates sp., the first mitochondrial genome sequenced from Enicocephalomorpha. The Stenopirates sp. mitochondrial genome is a typical circular DNA molecule of 15, 384 bp in length, and contains 37 genes and a large non-coding fragment. The gene order differs substantially from other known insect mitochondrial genomes, with rearrangements of both tRNA genes and protein-coding genes. The overall AT content (82.5%) of Stenopirates sp. is the highest among all the known heteropteran mitochondrial genomes. The strand bias is consistent with other true bugs with negative GC-skew and positive AT-skew for the J-strand. The heteropteran mitochondrial atp8 exhibits the highest evolutionary rate, whereas cox1 appears to have the lowest rate. Furthermore, a negative correlation was observed between the variation of nucleotide substitutions and the GC content of each protein-coding gene. A microsatellite was identified in the putative control region. Finally, phylogenetic reconstruction suggests that Enicocephalomorpha is the sister group to all the remaining Heteroptera.     

The mitochondrial genome of the firefly, Pyrocoelia rufa: complete DNA sequence, genome organization, and phylogenetic analysis with other insects

The complete nucleotide sequences of the mt genome from the firefly, Pyrococelia rufa (Coeleoptera: Lampyridae) was determined. The circular genome is 17,739-bp long, and contains a typical gene complement, order, and arrangement identical to Drosophila yacuba. The presence of 1,724-bp long intergenic spacer in the P. rufa mt genome is unique. The putative initiation codon for ND1 gene appears to be TTG, instead of frequently found ATN. All tRNAs showed stable canonical clover-leaf structure of other mt tRNAs, except for tRNA(Ser) (AGN), DHU arm of which could not form stable stem-loop structure. Phylogenetic analysis among insect orders confirmed a monophyletic Endopterygota, a monophyletic Mecopterida, a monophyletic Diptera, a monophyletic Lepidoptera, and a monophyletic Coleoptera, suggesting that the complete insect mt genome sequence has a resolving power in the diversification events within Endopterygota. However, internal relationships among three coleopteran species are not clear, and the inclusion of some insect orders (i.e., apterygotan T. gertschi) in the analysis provided inconsistent results compared to other molecular studies.     

Complete mitochondrial genome of the darkling beetle Gonocephalum outreyi (Coleoptera: Tenebrionidae) with phylogenetic implications

The complete mitochondrial genome (mitogenome) of Gonocephalum outreyi was determined by using next-generation sequencing approach. The full length of this mitogenome is 15,836?bp, which consists of 37 typical metazoan mitochondrial genes with an identical genome organization to ancestral insects. The majority of the protein-coding genes begin with the codon ATN, except for cox1 and cox2 with AAT and AAA, respectively. To elucidate the phylogenetic position of G. outreyi, we used various sequence coding schemes for protein-coding genes and the combined nucleotide sequences of all mitochondrial genes for tree building under the Bayesian and Maximum Likelihood inferences. The phylogenetic results consistently supported G. outreyi as a member of the family Tenebrionidae. The monophyly of both Tenebrionoidea and Tenebrionidae were strongly supported. The Scraptiidae and Melandryidae were recovered to be non-monophyletic in regards to the Osphya. Within Tenebrionidae, the subfamilies Diaperinae and Tenebrioninae were found to be non-monophyletic.