The complete mitochondrial genome of the Sichuan taimen (Hucho bleekeri): repetitive sequences in the control region and phylogenetic implications for Salmonidae

The complete mitochondrial DNA genome of the Sichuan taimen (Hucho bleekeri) was determined by the long and accurate polymerase chain reaction (LA-PCR) and primer walking sequence method. The entire mitochondrial genome of this species is 16,997 bp in length, making it the longest among the completely sequenced Salmonidae mitochondrial genomes. It consists of two ribosomal RNA (rRNA) genes, 13 protein-coding genes, 22 transfer RNA (tRNA) genes, and one control region (CR). The gene arrangement, nucleotide composition, and codon usage pattern of the mitochondrial genome are similar to those of other teleosts. A T-type mononucleotide microsatellite and an 82 bp tandem repeat were identified in the control region, which were almost identical among the three H. bleekeri individuals examined. Both phylogenetic analyses based on 12 concatenated protein-coding genes of the heavy strand and on just the control region show that H. bleekeri is a basal species in Salmoninae. In addition, Salmo, Salvelinus and Oncorhynchus all represent monophyletic groups, respectively. All freshwater species occupied basal phylogenetic positions, and also possessed various tandem repeats in their mitochondrial control regions. These results support established phylogenetic relationships among genera in Salmonidae based on morphological and molecular analyses, and are consistent with the hypothesis that Salmonidae evolved from freshwater species.     

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.     

Comparison of complete mitochondrial genomes of pine wilt nematode Bursaphelenchus xylophilus and Bursaphelenchus mucronatus (Nematoda: Aphelenchoidea) and development of a molecular tool for species identification

We determined the complete mitochondrial genome sequences for Bursaphelenchus mucronatus, one species of pinewood nematode. The genome is a circular-DNA molecule of 14,583 bp (195 bp smaller than its congener Bursaphelenchus xylophilus) and contains 12 protein-coding genes (lacking atp8), 22 tRNA genes, and 2 rRNA genes encoded in the same direction, consistent with most other nematodes. Based on sequence comparison of mtDNA genomes, we developed a PCR-based molecular assay to differentiate B. xylophilus (highly pathogenic) and B. mucronatus (relatively less virulent) using species-specific primers. The molecular identification system employs multiplex-PCR and is very effective and reliable for discriminating these Bursaphelenchus species, which are economically important, but difficult to distinguish based on morphology. The comparison of the mitochondrial genomes and molecular identification system of the two species of Bursaphelenchus spp. should provide a rich source of genetic information to support the effective control and management (quarantine) of the pine wilt disease caused by pinewood nematodes.     

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.     

The complete mitochondrial genome of Arctic Calanus hyperboreus (Copepoda,Calanoida) reveals characteristic patterns in calanoid mitochondrial genome

Copepoda is the most diverse and abundant group of crustaceans, but its phylogenetic relationships are ambiguous. Mitochondrial (mt) genomes are useful for studying evolutionary history, but only six complete Copepoda mt genomes have been made available and these have extremely rearranged genome structures. This study determined the mt genome of Calanus hyperboreus, making it the first reported Arctic copepod mt genome and the first complete mt genome of a calanoid copepod. The mt genome of C. hyperboreus is 17,910 bp in length and it contains the entire set of 37 mt genes, including 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. It has a very unusual gene structure, including the longest control region reported for a crustacean, a large tRNA gene cluster, and reversed GC skews in 11 out of 13 protein-coding genes (84.6%). Despite the unusual features, comparing this genome to published copepod genomes revealed retained pan-crustacean features, as well as a conserved calanoid-specific pattern. Our data provide a foundation for exploring the calanoid pattern and the mechanisms of mt gene rearrangement in the evolutionary history of the copepod mt genome.     

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.     

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.     

The complete mitochondrial genome of Biston panterinaria (Lepidoptera: Geometridae), with phylogenetic utility of mitochondrial genome in the Lepidoptera

The complete mitochondrial genome (mitogenome) of the Chinese pistacia looper Biston panterinaria was sequenced and annotated (15,517 bp). It contains the typical 37 genes of animal mitogenomes and a high A + T content (79.5%). All protein coding genes (PCGs) use standard ATN initiation codons except for cytochrome c oxidase 1 (COX1) with CGA. Eleven PCGs use a common stop codon of TAA or TAG, whereas COX2 and NADH dehydrogenase 4 (ND4) use a single T. All transfer RNA (tRNA) genes have the typical clover-leaf structure with the exception of tRNA^Ser(AGN). We reconstructed a preliminary mitochondrial phylogeny of six ditrysian superfamilies and performed comparative analyses of inference methods (Bayesian Inference (BI), Maximum Likelihood (ML), and Maximum Parsimony (MP)), dataset compositions (including and excluding 3rd codon positions), and alignment methods (Muscle, Clustal W, and MAFFT). Our analyses indicated that inference methods and dataset compositions more significantly affected the phylogenetic results than alignment methods. BI analysis consistently revealed uncontroversial relationships with all dataset compositions. By contrast, ML analysis failed to reconstruct stable phylogeny at two nodes, whereas MP analysis had more difficulties in the tree resolution and nodal support. Distinct from most previous studies, our analyses revealed that Geometroidea had a closer lineage relationship with Bombycoidea than Noctuoidea. Similar to previous molecular studies, our analyses revealed that Hesperiidae were nested in the Papilionoidea clade, providing further evidence to the previous concept that Papilionoidea was paraphyletic, and none of the butterflies were associated with the Macroheterocera.     

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.     

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.     

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 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.     

Complete Mitochondrial Genome of Haplorchis taichui and Comparative Analysis with Other Trematodes

Mitochondrial genomes have been extensively studied for phylogenetic purposes and to investigate intra- and interspecific genetic variations. In recent years, numerous groups have undertaken sequencing of platyhelminth mitochondrial genomes. Haplorchis taichui (family Heterophyidae) is a trematode that infects humans and animals mainly in Asia, including the Mekong River basin. We sequenced and determined the organization of the complete mitochondrial genome of H. taichui. The mitochondrial genome is 15,130 bp long, containing 12 protein-coding genes, 2 ribosomal RNAs (rRNAs, a small and a large subunit), and 22 transfer RNAs (tRNAs). Like other trematodes, it does not encode the atp8 gene. All genes are transcribed from the same strand. The ATG initiation codon is used for 9 protein-coding genes, and GTG for the remaining 3 (nad1, nad4, and nad5). The mitochondrial genome of H. taichui has a single long non-coding region between trnE and trnG. H. taichui has evolved as being more closely related to Opisthorchiidae than other trematode groups with maximal support in the phylogenetic analysis. Our results could provide a resource for the comparative mitochondrial genome analysis of trematodes, and may yield genetic markers for molecular epidemiological investigations into intestinal flukes.     

The complete mitochondrial genome and phylogenetic analysis of Nyctereutes procyonoides

The complete mitochondrial genome sequence of the raccoon dog (Nyctereutes procyonoides) was determined by using the long and accurate polymerase chain reaction. The entire mitochondrial genome sequence is 16,713 bp in length contains two ribosomal RNA genes, 13 protein-coding genes, 22 transfer RNA genes and 1 control region. Most mitochondrial genes are encoded on the H strand, except for the ND6 gene and 8 tRNA genes. The base compositions of mitochondrial genomes present clearly A–T skew. All the transfer RNA genes can be folded into the typical cloverleaf-shaped structure except tRNA-Ser (AGY), which lacks the dihydrouridine arm. Protein-coding genes mainly initiate with ATG and terminate with TAA. Some reading frame intervals and overlaps are found in the mitochondrial genome. The control region can be divided into three domains: the extended termination associated sequences (ETASs) domain, the central conserved domain and the conserved sequence blocks (CSBs) domain. Three conserved sequence blocks (CSBs) and one extended termination associated sequences (ETAS-1) is found in the control region. The phylogenetic analysis based on the concatenated data set of 14 genes in the mitochondrial genome of Canidae shows that the raccoon dog has close phylogenetic position with the red fox (Vulpes vulpes) and they constitute a clade which has an equil evolutionary position with the clade formed by the genera Canis and Cuon.     

The mitochondrial genome of the Russian wheat aphid Diuraphis noxia: Large repetitive sequences between trnE and trnF in aphids

To characterize aphid mitochondrial genome (mitogenome) features, we sequenced the complete mitogenome of the Russian wheat aphid, Diuraphis noxia. The 15,784-bp mitogenome with a high A + T content (84.76%) and strong C skew (− 0.26) was arranged in the same gene order as that of the ancestral insect. Unlike typical insect mitogenomes, D. noxia possessed a large tandem repeat region (644 bp) located between trnE and trnF. Sequencing partial mitogenome of the cotton aphid (Aphis gossypii) further confirmed the presence of the large repeat region in aphids, but with different repeat length and copy number. Another motif (58 bp) tandemly repeated 2.3 times in the control region of D. noxia. All repeat units in D. noxia could be folded into stem-loop secondary structures, which could further promote an increase in copy numbers. Characterization of the D. noxia mitogenome revealed distinct mitogenome architectures, thus advancing our understanding of insect mitogenomic diversities and evolution.     

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.     

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))))))).     

Variation in Coding (NADH Dehydrogenase Subunits 2, 3, and 6) and Noncoding Intergenic Spacer Regions of the Mitochondrial Genome in Octocorallia (Cnidaria: Anthozoa)

Low rates of evolution in cnidarian mitochondrial genes such as COI and 16S rDNA have hindered molecular systematic studies in this important invertebrate group. We sequenced fragments of 3 mitochondrial protein-coding genes (NADH dehydrogenase subunits ND2, ND3 and ND6) as well as the COI-COII intergenic spacer, the longest noncoding region found in the octocoral mitochondrial genome, to determine if any of these regions contain levels of variation sufficient for reconstruction of phylogenetic relationships among genera of the anthozoan subclass Octocorallia. Within and between the soft coral families Alcyoniidae and Xeniidae, sequence divergence in the genes ND2 (539 bp), ND3 (102 bp), and ND6 (444 bp) ranged from 0.5% to 12%, with the greatest pairwise distances between the 2 families. The COI-COII intergenic spacer varied in length from 106 to 122 bp, and pairwise sequence divergence values ranged from 0% to 20.4%. Phylogenetic trees constructed using each region separately were poorly resolved. Better phylogenetic resolution was obtained in a combined analysis using all 3 protein-coding regions (1085 bp total). Although relationships among some pairs of species and genera were well supported in the combined analysis, the base of the alcyoniid family tree remained an unresolved polytomy. We conclude that variation in the NADH subunit coding regions is adequate to resolve phylogenetic relationships among families and some genera of Octocorallia, but insufficient for most species - or population-level studies. Although the COI-COII intergenic spacer exhibits greater variability than the protein-coding regions and may contain useful species-specific markers, its short length limits its phylogenetic utility.     

The complete mitochondrial genome of Microtus fortis calamorum (Arvicolinae, Rodentia) and its phylogenetic analysis

Microtus fortis is a special resource of rodent in China. It is a promising experimental animal model for the study on the mechanism of Schistosome japonicum resistance. The first complete mitochondrial genome sequence for Microtus fortis calamorum, a subspecies of M. fortis (Arvicolinae, Rodentia), was reported in this study. The mitochondrial genome sequence of M. f. calamorum (Genbank: JF261175) showed a typical vertebrate pattern with 13 protein coding genes, 2 ribosomal RNAs, 22 transfer RNAs and one major noncoding region (CR region).The extended termination associated sequences (ETAS-1 and ETAS-2) and conserved sequence block 1 (CSB-1) were found in the CR region. The putative origin of replication for the light strand (O(L)) of M. f. calamorum was 35bp long and showed high conservation in stem and adjacent sequences, but the difference existed in the loop region among three species of genus Microtus. In order to investigate the phylogenetic position of M. f. calamorum, the phylogenetic trees (Maximum likelihood and Bayesian methods) were constructed based on 12 protein-coding genes (except for ND6 gene) on H strand from 16 rodent species. M. f. calamorum was classified into genus Microtus, Arvcicolinae for the highly phylogenetic relationship with Microtus kikuchii (Taiwan vole). Further phylogenetic analysis results based on the cytochrome b gene ranged from M. f. calamorum to one of the subspecies of M. fortis, which formed a sister group of Microtus middendorfii in the genus Microtus.