Sequence comparison of the mitochondrial genomes of five brackish water species of the family Neritidae: Phylogenetic implications and divergence time estimation

Neritids are ancient gastropod species which can live in marine, brackish water, and freshwater environments. In this study, we sequenced and annotated the mitochondrial genomes of five brackish water neritids (i.e., Clithon corona, Clithon lentiginosum, Clithon squarrosum, Neritina iris, and Septaria lineata). The mitogenomes ranged from 15,618 to 15,975 bp, and all contain 13 protein‐coding genes (PCGs), 22 tRNA genes, and two rRNA genes, with a closed ring structure. We calculated the Ka/Ks values of all 13 PCGs of Neritidae species, all ratios are less than 1, under purification selection. Phylogenetic analysis of the 13 PCGs showed that Neritimorpha is a sister group with Vetigastropoda and Caenogastopoda, genus Clithon is a sister group with Neritina and Septaria. Estimation of divergence time for all species of Neritidae showed that the main differentiation of Neritidae occurred in Cenozoic period (65 Mya), C. corona and C. lentiginosum were differentiated in the Cenozoic Neogene, the other three species diverged in the Cenozoic Paleogene. These results will help to better understand the evolutionary position of Neritidae and provide reference for further phylogenetic research on Neritidae species.     

Characteristics of the mitochondrial genome of Rana omeimontis and related species in Ranidae: Gene rearrangements and phylogenetic relationships

The Omei wood frog (Rana omeimontis), endemic to central China, belongs to the family Ranidae. In this study, we achieved detail knowledge about the mitogenome of the species. The length of the genome is 20,120 bp, including 13 protein‐coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a noncoding control region. Similar to other amphibians, we found that only nine genes (ND6 and eight tRNA genes) are encoded on the light strand (L) and other genes on the heavy strand (H). Totally, The base composition of the mitochondrial genome included 27.29% A, 28.85% T, 28.87% C, and 15.00% G, respectively. The control regions among the Rana species were found to exhibit rich genetic variability and A + T content. R. omeimontis was clustered together with R. chaochiaoensis in phylogenetic tree. Compared to R. amurensis and R. kunyuensi, it was more closely related to R. chaochiaoensis, and a new way of gene rearrangement (ND6‐trnE‐Cytb‐D‐loop‐trnL2 (CUN)‐ND5‐D‐loop) was also found in the mitogenome of R. amurensis and R. kunyuensi. Our results about the mitochondrial genome of R. omeimontis will contribute to the future studies on phylogenetic relationship and the taxonomic status of Rana and related Ranidae species.     

Characterization of the Complete Mitochondrial Genome of Cerura menciana and Comparison with Other Lepidopteran Insects

The complete mitochondrial genome (mitogenome) of Cerura menciana (Lepidoptera: Notodontidae) was sequenced and analyzed in this study. The mitogenome is a circular molecule of 15,369 bp, containing 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and a A+T-rich region. The positive AT skew (0.031) indicated that more As than Ts were present. All PCGs were initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which was initiated by CAG. Two of the 13 PCGs contained the incomplete termination codon T or TA, while the others were terminated with the stop codon TAA. The A+T-rich region was 372 bp in length and consisted of an ‘ATAGA’ motif followed by an 18 bp poly-T stretch, a microsatellite-like (AT)₈ and a poly-A element upstream of the trnM gene. Results examining codon usage indicated that Asn, Ile, Leu2, Lys, Tyr and Phe were the six most frequently occurring amino acids, while Cys was the rarest. Phylogenetic relationships, analyzed based on the nucleotide sequences of the 13 PCGs from other insect mitogenomes, confirmed that C. menciana belongs to the Notodontidae family.     

The first complete mitochondrial genome from Bostrychus genus (Bostrychus sinensis) and partitioned Bayesian analysis of Eleotridae fish phylogeny

To understand the phylogenetic position of Bostrychus sinensis in Eleotridae and the phylogenetic relationships of the family, we determined the nucleotide sequence of the mitochondrial (mt) genome of Bostrychus sinensis. It is the first complete mitochondrial genome sequence of Bostrychus genus. The entire mtDNA sequence was 16508 bp in length with a standard set of 13 protein-coding genes, 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a noncoding control region. The mitochondrial genome of B. sinensis had common features with those of other bony fishes with respect to gene arrangement, base composition, and tRNA structures. Phylogenetic hypotheses within Eleotridae fish have been controversial at the genus level. We used the mitochondrial cytochrome b (cytb) gene sequence to examine phylogenetic relationships of Eleotridae by using partitioned Bayesian method. When the specific models and parameter estimates were presumed for partitioning the total data, the harmonic mean –lnL was improved. The phylogenetic analysis supported the monophyly of Hypseleotris and Gobiomorphs. In addition, the Bostrychus were most closely related to Ophiocara, and the Philypnodon is also the sister to Microphlypnus, based on the current datasets. Further, extensive taxonomic sampling and more molecular information are needed to confirm the phylogenetic relationships in Eleotridae.     

The complete mitochondrial genome of the black mud crab, Scylla serrata (Crustacea: Brachyura: Portunidae) and its phylogenetic position among (pan)crustaceans

The black mud crab, Scylla serrata (Forsk?l 1775), is the most economically important edible crab in South-East Asia. In the present study, the complete mitochondrial genome of black mud crab, S.?serrata, was determined with the sequential polymerase chain reaction and primer walking sequencing. The complete mitochondrial genome was 15,721?bp in length with an A+T content of 69.2?% and contained 37 mitochondrial genes (13 protein coding genes (PCGs), 2 ribosomal RNA genes and 22 transfer RNA genes) and a control region (CR). The analysis of the CR sequence shows that it contains a multitude of repetitive fragments which can fold into hairpin-like or secondary structures and conserved elements as in other arthropods. The gene order of S. serrata mainly retains as the pancrustacean ground pattern, except for a single translocation of trnH. The gene arrangement of S. serrata appears to be a typical feature of portunid crabs. Phylogenetic analyses with concatenated amino acid sequences of 12 PCGs establishes that S. serrata in a well-supported monophyletic Portunidae and is consistent with previous morphological classification. Moreover, the phylogenomic results strongly support monophyletic Pancrustacea (Hexapoda plus “Crustaceans”). Within Pancrustacea, this study identifies Malacostraca?+?Entomostraca and Branchiopoda as the sister group to Hexapoda, which confirms that “Crustacea” is not monophyletic. Cirripedia?+?Remipedia appear to be a basal lineage of Pancrustacea. The present study also provides considerable data for the application of both population and phylogenetic studies of other crab species.     

Comparative plastome analysis of Blumea,with implications for genome evolution and phylogeny of Asteroideae

The genus Blumea (Asteroideae, Asteraceae) comprises about 100 species, including herbs, shrubs, and small trees. Previous studies have been unable to resolve taxonomic issues and the phylogeny of the genus Blumea due to the low polymorphism of molecular markers. Therefore, suitable polymorphic regions need to be identified. Here, we de novo assembled plastomes of the three Blumea species B. oxyodonta, B. tenella, and B. balsamifera and compared them with 26 other species of Asteroideae after correction of annotations. These species have quadripartite plastomes with similar gene content, genome organization, and inverted repeat contraction and expansion comprising 113 genes, including 80 protein‐coding, 29 transfer RNA, and 4 ribosomal RNA genes. The comparative analysis of codon usage, amino acid frequency, microsatellite repeats, oligonucleotide repeats, and transition and transversion substitutions has revealed high resemblance among the newly assembled species of Blumea. We identified 10 highly polymorphic regions with nucleotide diversity above 0.02, including rps16‐trnQ, ycf1, ndhF‐rpl32, petN‐psbM, and rpl32‐trnL, and they may be suitable for the development of robust, authentic, and cost‐effective markers for barcoding and inference of the phylogeny of the genus Blumea. Among these highly polymorphic regions, five regions also co‐occurred with oligonucleotide repeats and support use of repeats as a proxy for the identification of polymorphic loci. The phylogenetic analysis revealed a close relationship between Blumea and Pluchea within the tribe Inuleae. At tribe level, our phylogeny supports a sister relationship between Astereae and Anthemideae rooted as Gnaphalieae, Calenduleae, and Senecioneae. These results are contradictory to recent studies which reported a sister relationship between “Senecioneae and Anthemideae” and “Astereae and Gnaphalieae” or a sister relationship between Astereae and Gnaphalieae rooted as Calenduleae, Anthemideae, and then Senecioneae using nuclear genome sequences. The conflicting phylogenetic signals observed at the tribal level between plastidt and nuclear genome data require further investigation.     

Molecular phylogeny and the underestimated species diversity of the endemic white‐bellied rat (Rodentia: Muridae: Niviventer) in Southeast Asia and China

The white‐bellied rat, Niviventer, is a genus endemic to Southeast Asia and China. However, the interspecific phylogenetic relationships and species diversity of this genus remain poorly understood. In the present study, single and multi‐locus analyses were performed. Phylogenetic reconstruction on Cytochrome b (512 individuals, including data from Genbank) revealed five major clades with approximately 35 operational taxonomic units (OTUs), a number twice the existing taxonomy. The first clade (N. langbianis species group) was the earliest diverged. The second clade (N. fulvescens species group) diverged in Southeast Asia, the south and lower altitude regions of the Hengduan Mountains, and Southeast China. The third clade (the N. eha species group) is endemic to high altitudes in Northwest Yunnan and the central region of Himalaya. The fourth clade (the N. andersoni species group), is mainly confined to alpine regions of the Hengduan Mountains. The fifth clade (N. confucianus species group) is mainly distributed in the north and higher altitude regions of eastern Himalaya, the Hengduan Mountains and Taiwan, with the complex also invading central and northern China. Results from the combined dataset of four genes (Cytochrome b, Cytochrome oxidase subunit I, the D‐loop sequence of the mitochondrial genome and the first exon of the nuclear interphotoreceptor retinoid binding protein) for 82 representative individuals from China generally coincide with the result of the single gene, with 12 OTUs identified. These results provide a preliminary framework for the existing classification of this highly diversified genus. The divergence time of Niviventer based on the four gene topology was dated to the late Miocene ~6.41 Ma. Significant differences were detected in the general body form changes among these units based on voucher specimens. Moreover, geometric morphometric analysis of the cranium shape of voucher specimens indicated significant differences among five major species groups. Shape divergence of the cranium among several OTUs within the N. confucinaus complex is also significant. Our results provide further evidence for rapid and highly underestimated diversification of Niviventer both in genetics and morphology.     

Characterization of the complete mitochondrial genome of Drawida gisti (Metagynophora,Moniligastridae) and comparison with other Metagynophora species

Here, the complete mitochondrial genome (mitogenome) of Drawida gisti was sequenced and compared with the mitogenomes of other Metagynophora species. The circular mitogenome was 14,648 bp in length and contained two ribosomal RNA genes (rRNAs), 13 protein-coding genes (PCGs), and 22 transfer RNA genes (tRNAs). The types of constitutive genes and the direction of the coding strand that appeared in Drawida mitogenome were identical to those observed in other Metagynophora species, except for a missing lengthy non-coding region. The conservative relationships between Drawida species were supported by the overall analyses of 13 PCGs, two rRNAs, and 22 tRNAs. A comparison of the Metagynophora mitogenomes revealed that the ATP8 gene possessed the highest polymorphism among the 13 PCGs and two rRNAs. Phylogenetic analysis suggested that the Moniligastridae contained Drawida, which is a primitive Metagynophora group. Our study provides a step forward toward elucidating the evolutionary linkages within Drawida and even Metagynophora.     

The complete mitochondrial genome of Parasesarma pictum (Brachyura: Grapsoidea: Sesarmidae) and comparison with other Brachyuran crabs

Mitochondrial DNA (mtDNA) is an extrachromosomal genome which can provide important information for evolution and phylogenetic analysis. In this study, we assembled a complete mitogenome of a crab Parasesarma pictum (Brachyura: Grapsoidea: Sesarmidae) from next generation sequencing reads at the first time. P. pictum is a mudflat crab, belonging to the Sesarmidae family (subfamily Sesarminae), which is perched on East Asia. The 15,716 bp mitogenome covers 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and one control region (CR). The control region spanns 420 bp. The genome composition was highly A+T biased 75.60% and showed negative AT-skew (−0.03) and negative GC-skew (−0.47). Compared with the ancestor of Brachyura, the gene order of Sesarmidae has several differences and the gene order of P. pictum is typical for mitogenomes of Sesarmidae. Phylogenetic tree based on nucleotide sequences of mitochondrial 13 PCGs using BI and ML determined that P. pictum has a sister group relationship with Parasesarma tripectinis and belongs to Sesarmidae.     

The characterization of complete mitochondrial genome and phylogenetic relationship within Rattus genus (Rodentia: Muridae)

The genus Rattus is one of the main pest genus of rodent. Most species of the genus carry all kinds of pathogenic bacteria to human being. They are traditionally considered to be a least understood group. The complete mitochondrial genome of the White-Footed Indochinese Rat, Rattus nitidus was determined in this study. The characterization of mitochondrial genomes of Rattus genus was also analyzed based on comprehensive comparison. The result of evolutionary patterns of protein-coding genes (PCGs) suggested purifying selection was the predominant evolutionary forces in the mitochondrial genomes of Rattus genus. The NADH dehydrogenase 4 gene (ND4) showed a highly elevated Ka/Ks ratio compared to the other protein-coding genes, which indicated ND4 was most likely under relaxed selection pressure. Phylogenetic analysis provided a well-supported outline of Rattus genus, and revealed two groups in the genus. R. nitidus had a sister relationship with R. norvegicus.     

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 the Geophilomorph Centipede Strigamia maritima

Strigamia maritima (Myriapoda; Chilopoda) is a species from the soil-living order of geophilomorph centipedes. The Geophilomorpha is the most speciose order of centipedes with over a 1000 species described. They are notable for their large number of appendage bearing segments and are being used as a laboratory model to study the embryological process of segmentation within the myriapods. Using a scaffold derived from the recently published genome of Strigamia maritima that contained multiple mitochondrial protein-coding genes, here we report the complete mitochondrial genome of Strigamia, the first from any geophilomorph centipede. The mitochondrial genome of S. maritima is a circular molecule of 14,938 base pairs, within which we could identify the typical mitochondrial genome complement of 13 protein-coding genes and 2 ribosomal RNA genes. Sequences resembling 16 of the 22 transfer RNA genes typical of metazoan mitochondrial genomes could be identified, many of which have clear deviations from the standard ‘cloverleaf’ secondary structures of tRNA. Phylogenetic trees derived from the concatenated alignment of protein-coding genes of S. maritima and >50 other metazoans were unable to resolve the Myriapoda as monophyletic, but did support a monophyletic group of chilopods: Strigamia was resolved as the sister group of the scolopendromorph Scolopocryptos sp. and these two (Geophilomorpha and Scolopendromorpha), along with the Lithobiomorpha, formed a monophyletic group the Pleurostigmomorpha. Gene order within the S. maritima mitochondrial genome is unique compared to any other arthropod or metazoan mitochondrial genome to which it has been compared. The highly unusual organisation of the mitochondrial genome of Strigamia maritima is in striking contrast with the conservatively evolving nuclear genome: sampling of more members of this order of centipedes will be required to see whether this unusual organization is typical of the Geophilomorpha or results from a more recent reorganisation in the lineage leading to Strigamia.     

The First Mitochondrial Genome for Caddisfly (Insecta: Trichoptera) with Phylogenetic Implications

The Trichoptera (caddisflies) is a holometabolous insect order with 14,300 described species forming the second most species-rich monophyletic group of animals in freshwater. Hitherto, there is no mitochondrial genome reported of this order. Herein, we describe the complete mitochondrial (mt) genome of a caddisfly species, Eubasilissa regina (McLachlan, 1871). A phylogenomic analysis was carried out based on the mt genomic sequences of 13 mt protein coding genes (PCGs) and two rRNA genes of 24 species belonging to eight holometabolous orders. Both maximum likelihood and Bayesian inference analyses highly support the sister relationship between Trichoptera and Lepidoptera.     

The First Mitochondrial Genome for the Superfamily Hagloidea and Implications for Its Systematic Status in Ensifera

Hagloidea Handlirsch, 1906 was an ancient group of Ensifera, that was much more diverse in the past extending at least into the Triassic, apparently diminishing in diversity through the Cretaceous, and now only represented by a few extant species. In this paper, we report the complete mitochondrial genome (mitogenome) of Tarragoilus diuturnus Gorochov, 2001, representing the first mitogenome of the superfamily Hagloidea. The size of the entire mitogenome of T. diuturnus is 16144 bp, containing 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one control region. The order and orientation of the gene arrangement pattern is identical to that of D. yakuba and most ensiferans species. A phylogenomic analysis was carried out based on the concatenated dataset of 13 PCGs and 2 rRNA genes from mitogenome sequences of 15 ensiferan species, comprising four superfamilies Grylloidea, Tettigonioidae, Rhaphidophoroidea and Hagloidea. Both maximum likelihood and Bayesian inference analyses strongly support Hagloidea T. diuturnus and Rhaphidophoroidea Troglophilus neglectus as forming a monophyletic group, sister to the Tettigonioidea. The relationships among four superfamilies of Ensifera were (Grylloidea, (Tettigonioidea, (Hagloidea, Rhaphidophoroidea))).     

Complete mitochondrial genome of Dacus conopsoides (Insecta: Tephritidae) with tRNA gene duplication and molecular phylogeny of Dacini tribe

To date there is only a single report on the complete mitochondrial genome of the Dacus fruit flies. We report here the whole mitogenome of Dacus conopsoides with first report of tRNA gene duplication in tephritid fruit flies determined using next-generation sequencing and discuss the molecular phylogeny of Dacini tribe. It had a total length of 15,852 bp, comprising 13 protein coding genes, 2 rRNA genes, 23 tRNA genes, and a non-coding region (A + T-rich control region). The 65-bp trnF gene was duplicated, and the 68-bp trnE gene was partially duplicated resulting in a 31-bp pseudogene. The cloverleaf structure for trnN, trnH, and trnF lacked the TΨC-loop, while trnS lacked the D-stem. The start codons for the protein coding genes included 6 ATG, 3 ATC, 2 ATA, and 1 each of ATT and TCG. Seven PCGs had TAA stop codon, two had TAG and four had incomplete T stop codon. Molecular phylogeny based on 15 mt-genes (13 PCGs +2 rRNA genes) and 30 taxa of Tephritidae indicated D. conopsoides forming a monophyletic sister group with D. longicornis supported by high bootstrap value. The lineage containing also the monophyletic genus Zeugodacus. The Dacini and Ceratitidini tribes of the subfamily Dacinae were monophyletic but the subfamilies Dacinae and Trypetinae were paraphyletic. A broader taxa sampling of the Tephritidae is needed to better elucidate the phylogenetics and systematics of the tribes and subfamilies of tephritid fruit flies.     

Mitochondrial genomes of five Hyphessobrycon tetras and their phylogenetic implications

To date, the taxonomic status and phylogenetic affinities within Hyphessobrycon, even among other genera in Characidae, remain unclear. Here, we determined five new mitochondrial genomes (mitogenomes) of Hyphessobrycon species (H. elachys, H. flammeus, H. pulchripinnis, H. roseus, and H. sweglesi). The mitogenomes were all classical circular structures, with lengths ranging from 16,008 to 17,224 bp. The type of constitutive genes and direction of the coding strand that appeared in the mitogenomes were identical to those of other species in Characidae. The highest value of the Ka/Ks ratio within 13 protein‐coding genes (PCGs) was found in ND2 with 0.83, suggesting that they were subject to purifying selection in the Hyphessobrycon genus. Comparison of the control region sequences among seven Hyphessobrycon fish revealed that repeat units differ in length and copy number across different species, which led to sharp differences in mitogenome sizes. Phylogenetic trees based on the 13 PCGs did not support taxonomic relationships, as the Hyphessobrycon fish mixed with those from other genera. These data were combined to explore higher level relationships within Characidae and could aid in the understanding of the evolution of this group.     

Characterization of the complete mitochondrial genome of Bombyx mori strain H9 (Lepidoptera: Bombycidae)

The complete mitochondrial genome (mitogenome) of Bombyx mori strain H9 (Lepidoptera: Bombycidae) is 15,670 base pairs (bp) in length, encoding 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. The nucleotide composition of the genome is highly A + T biased, accounting for 81.31%, with a slightly positive AT skewness (0.059). The arrangement of 13 PCGs is similar to that of other sequenced lepidopterans. All the PCGs are initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which is proposed by the TTAG sequence as observed in other lepidopterans. Unlike the other PCGs, the cox1 and cytochrome c oxidase subunit 2 (cox2) genes have incomplete stop codons consisting of just a T. All tRNAs have typical structures of insect mitochondrial tRNAs, which is different from other sequenced lepidopterans. The structure of A + T-rich region is similar to that of other sequenced lepidopterans, including non-repetitive sequences, the ATAGA binding domain, a 18 bp poly-T stretch and a poly-A element upstream of transfer RNA M (trnM) gene. Phylogenetic analysis shows that the domesticated silkmoth B. mori originated from the Chinese Bombyx mandarina.     

Complete mitochondrial genome of Pycnonotus xanthorrhous (Passeriformes,Pycnonotidae) and phylogenetic consideration

The complete mitochondrial genome (mitogenome) of Pycnonotus xanthorrhous was sequenced via next generation sequencing. The full length of the circular genome is 16,952 bp. It consists of 37 typical animal mitochondrial genes including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) and 2 ribosomal RNA (rRNA) genes. P. xanthorrhous also contains one control region (CR) and one pseudo-control region, and shares the identical gene arrangements with sequenced Pycnonotus spp. which differs from the typical vertebrates gene order. Phylogenetic analyses indicates that Passerida sensu stricto contains three major clades and the core Sylvioidea form a monophyletic group. Furthermore, we investigated the evolution of control region within this lineage and revealed the multiple independent origins of duplicate control region.     

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.     

Mitochondrial genome of the yellow catfish Pelteobagrus fulvidraco and insights into Bagridae phylogenetics

The mitochondrial genome (mitogenome) can provide important information for understanding phylogenetic analysis and molecular evolution. Herein, we amplified the complete mitogenome sequence of Pelteobagrus fulvidraco. The mitogenome was 16,526 bp in length and included 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes and a non-coding control region (D-loop). Both the organization and location of genes in the mitogenome were consistent with those from Siluriformes fishes previously published in GenBank. The phylogenetic relationships based on Bayesian inference (BI) and Maximum likelihood (ML) methods showed that P. fulvidraco has close relationships with Pelteobagrus eupogon and Tachysurus intermedius, suggesting that P. fulvidraco belongs to Tachysurus. This study provides evidence that Tachysurus, Pseudobagrus and Leiocassis do not form monophyly, but that these three genera form a monophyletic group. Our results provide reference for further phylogenetic research of the Bagridae species.