Complete nucleotide sequence of the Coturnix chinensis (blue-breasted quail) mitochondrial genome and a phylogenetic analysis with related species

Coturnix chinensis (blue-breasted quail) has been classically grouped in Galliformes Phasianidae Coturnix, based on morphologic features and biochemical evidence. Since the blue-breasted quail has the smallest body size among the species of Galliformes, in addition to a short generation time and an excellent reproductive performance, it is a possible model fowl for breeding and physiological studies of the Coturnix japonica (Japanese quail) and Gallus gallus domesticus (chicken), which are classified in the same family as blue-breasted quail. However, since its phylogenetic position in the family Phasianidae has not been determined conclusively, the sequence of the entire blue-breasted quail mitochondria (mt) genome was obtained to provide genetic information for phylogenetic analysis in the present study. The blue-breasted quail mtDNA was found to be a circular DNA of 16,687 base pairs (bp) with the same genomic structure as the mtDNAs of Japanese quail and chicken, though it is smaller than Japanese quail and chicken mtDNAs by 10 bp and 88 bp, respectively. The sequence identity of all mitochondrial genes, including those for 12S and 16S ribosomal RNAs, between blue-breasted quail and Japanese quail ranged from 84.5% to 93.5%; between blue-breasted quail and chicken, sequence identity ranged from 78.0% to 89.6%. In order to obtain information on the phylogenetic position of blue-breasted quail in Galliformes Phasianidae, the 2,184 bp sequence comprising NADH dehydrogenase subunit 2 and cytochrome b genes available for eight species in Galliformes [Japanese quail, chicken, Gallus varius (green junglefowl), Bambusicola thoracica (Chinese bamboo partridge), Pavo cristatus (Indian peafowl), Perdix perdix (gray partridge), Phasianus colchicus (ring-neck pheasant), and Tympanchus phasianellus (sharp-tailed grouse)] together with that of Aythya americana (redhead) were examined using a maximum likelihood (ML) method. The ML analyses on the first/second codon positions, the third codon positions, and amino acid sequence consistently demonstrated that blue-breasted quail and Japanese quail are in the same phylogenetic cluster.     

Phylogenetic analyses of the harvest mouse,Micromys minutus (Rodentia: Muridae) based on the complete mitogenome sequences

Harvest Mouse (Micromys minutus) has a very wide range of distribution in Asia and Europe. However, the phylogenetic relationship of M. minutus is still uncertain. In this study, we determined the complete mitochondrial (mt) genome sequences of M. minutus, and used the complete mitochondrial genome sequences constructed the phylogenetic tree of Muroidea. The size of the genome is 16,232 bp in length and has a base composition of 33.6% A, 29.1% T, 24.8% C, and 12.5% G. The mitogenome structure was similar to that of typical vertebrate and other rodents' mitochondrial genomes, includes 13 protein-coding genes, 2 rRNA genes (12S rRNA and 16S rRNA), 22 tRNA genes, and 1 control region. We suggested a new initiation codon for ND5 (NADH dehydrogenase subunit), which has been never reported in the mitochondrial genome of vertebrate. The ML and BI phylogenetic trees, which based on the combination of the 12 protein-coding genes, supported strongly that the genus Micromys was represent an early offshoot within the Muridae with high support values (BI = 1.00, ML = 100).     

<Emphasis Type="Italic">Gallus gallus</Emphasis> aggrecan gene-based phylogenetic analysis of selected avian taxonomic groups

Mitochondrial DNA (mtDNA) sequences remain the most widely used for phylogenetic analysis in birds. A major limitation of mtDNA sequences, however, is that mitochondria genes are inherited as a single linkage group. Here we describe the use of a 540-bp DNA sequence corresponding to the G3 domain of Gallus gallus nuclear aggrecan gene (AGC1) for phylogenetic analysis of the main groups of Galliformes including Phasianidae, Numididae, and Odontophoridae. We also included species from Cracidae and Megapodiidae which are considered by some as Craciformes and others, including here as Galliformes. The uncorrected sequence divergence of the G3 fragments ranges from 1 among the grouses to 36% between some of the distant groups within Galliformes. These sequences contain 39–48% AT nucleotides and the ratios of transition versus transversion are above 1.5 in majority of the comparisons. Using G3 sequences from an Anseriform, Oxyura jamaicensis, as out-groups, phylogenetic trees were obtained using maximum parsimony and distance algorithms and bootstrap analyses. These trees were consistent with those described using Avian sarcoma and leucosis virus gag genes and those from amino acid sequences of hemoglobin and lysozyme c. Our data also support relationships among Galliformes which were defined using mtDNA sequences. In addition to the general support of the five main families of Galliformes, our data are also consistent with previous work that showed Francolinus africanus and Gallus gallus are in the same clade and that Tetraoninae is a well-supported monophyletic subfamily within Phasianidae. The results presented here suggest that the AGC1 sequences meet the criterion of novel nuclear DNA sequences that can be used to help resolve the relationships among Galliformes.     

Complete mitochondrial genome of <Emphasis Type="Italic">Otis tarda</Emphasis> (Gruiformes: Otididae) and phylogeny of Gruiformes inferred from mitochondrial DNA sequences

The complete nucleotide sequence of mitochondrial genome of the Great bustard (Otis tarda) was determined by using polymerase chain reaction (PCR) method. The genome is 16,849 bp in size, containing 13 protein-coding, 2 ribosomal and 22 transfer RNA genes. Sequences of the tRNA genes can be folded into canonical cloverleaf secondary structure except for tRNA-Cys and tRNA-Ser (AGY), which lose “DHU” arm. Sequence analysis showed that the O. tarda mitochondrial control region (mtCR) contained many elements in common with other avian mtCRs. A microsatellite repeat was found in the 3′-peripheral domain of the O. tarda mtCR. Based on the mitochondrial DNA sequences of 12S rRNA, 16S rRNA and tRNA-Val, a phylogenetic study of Gruiformes was performed. The result showed that Otididae was a sister group to “core Gruiformes” and Charadriiformes with strong support (97% posterior probability values) in Bayesian analysis. The taxonomic status of Rhynochetidae, Mesitornithidae, Pedionomidae and Turnicidae that traditionally belonged to Gruiformes was also discussed in this paper.     

The complete mitochondrial genome of Chrysolophus pictus (Galliformes: Phasianidae) and a phylogenetic analysis with related species

The 16,678 bp mitochondrial genome of the Chrysolophus pictus has been sequenced in this paper. To determine the phylogentic position of C. pictus with related species within Phasianidae, the phylogenetic tree was reconstructed with the concatenated nucleotide dataset of the 12 heavy-strand-encoded protein genes. The phylogenetic analysis was carried out using maximum parsimony (MP) and Bayesian inference (BI) methods. MP and BI phylogenetic trees here showed similar topology and consistently suggested that C. pictus shared a close relationship with Phasianus versicolor. The results also showed that the Meleagris gallopavo possessed a basal phylogenetic position within Phasianidae, which may imply that it should be classified into the Phasianidae.     

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.     

Complete sequence of the Japanese quail (Coturnix japonica) mitochondrial genome and its genetic relationship with related species

The Japanese quail (Coturnix japonica; JQ) is one of the domesticated fowl species of Japan. To provide DNA sequence information for examination of its phylogenetic position in the order Galliformes, the complete sequence of the JQ mitochondria was determined. Sequence analysis revealed that the JQ mitochondrial genome is a circular DNA of 16 697 basepairs (bp), which is smaller than the chicken mitochondrial DNA of 16 775 bp, but the genomic structure of JQ mitochondria was the same as that of the chicken. The sequence homologies of all mitochondrial genes including those for 12S and 16S ribosomal RNA (rRNA), between Japanese quail and chicken ranged from 78.0 to 89.9%. Because the sequences of NADH dehydrogenase subunit 2 and cytochrome b genes had been reported in five species [Phasianus colchicus (ring-neck pheasant: RP), Gallus gallus domesticus (chicken: CH), Perdix perdix (grey partridge: GP), Bambusicola thoracia (Chinese bamboo partridge: CP), and Aythya americana (redhead: RH)], the concatenated nucleotide sequences (2184 bp) and amino acid sequences of these two genes were used in a phylogenetic analysis of JQ against these five species using a maximum likelihood (ML) method. Using the first and second bases of the codons, and the third base of the codons indicated a phylogenic tree of [RH, (RP, GP), (JQ, (CH, CP))]. A phylogenic tree of [RH, JQ, (RP, GP), (CH, CP)] was determined using amino acid sequences. Because the local bootstrap values for the JQ branch in these trees are not high, additional sequence is necessary for construction of a reliable tree.     

The first draft genome of Lophophorus: A step forward for Phasianidae genomic diversity and conservation

The monal genus (Lophophorus) is a branch of Phasianidae and its species inhabit the high-altitude mountains of the Qinghai-Tibet Plateau. The Chinese monal, L. lhuysii, is a threatened endemic bird of China that possesses high-altitude adaptability, diversity of plumage color and potentially low reproductive life history. This is the first study to describe the monal genome using next generation sequencing technology. The Chinese monal genome size is 1.01 Gb, with 16,940 protein-coding genes. Gene annotation yielded 100.93 Mb (9.97%) repeat elements, 785 ncRNA, 5,465,549 bp (0.54%) SSR and 15,550 (92%) genes in public databases. Compared to other birds and mammals, the genome evolution analysis showed numerous expanded gene families and positive selected genes involved in high-altitude adaptation, especially related to the adaptation of low temperature and hypoxia. Consequently, this gene data can be used to investigate the molecular evolution of high-altitude adaptation in future bird research. Our first published genome of the genus Lophophorus will be integral for the study of monal population genetic diversity and conservation, genomic evolution and Galliformes species differentiation in the Qinghai-Tibetan Plateau.     

Complete mitochondrial genome sequence of <Emphasis Type="Italic">Cucullaea labiata</Emphasis> (Arcoida: Cucullaeidae) and phylogenetic implications

The complete mitochondrial genome of Cucullaea labiata (Arcoida: Cucullaeidae) was firstly determined in this study in order to better understand the phylogenetic relationship between Cucullaeidae and Arcidae. The C. labiata mitochondrial genome was 25,845 bp in size and contained 12 protein-coding genes, 2 rRNA and 22 tRNA genes. The number and the location of the tRNA genes were different from three Arcidae species (Scapharca broughtonii, Scapharca kagoshimensis and Tegillarca granosa). Gene arrangement also differed dramatically. The length of the non-coding regions was 10,559 bp, in which the largest one (6057 bp) included eight point nine copies of a 659 bp repeat motif. The number of repeated sequences was different in different individuals, similar to the findings from the mitochondrial genome of S. broughtonii and Placopecten magellanicus. One intron was found in cox1 gene both in CL_98 and in CL_99 individuals of C. labiata. The reason why mitochondrial introns are retained so scarcely in bivalve taxa needs further research. Phylogenetic analyses based on 12 concatenated amino acid sequences of protein-coding genes supported Cucullaeidae was the sister group of Arcidae.     

Characterization of the complete mitochondrial genome of Hydrotaea spinigera (Diptera: Muscidae) with phylogenetic implications

Hydrotaea spinigera Stein is a muscid species with high forensic importance. To further supplement the genome-level features of related species, the complete mitochondrial genome of H. spinigera is amplified, sequenced, annotated, analyzed, and compared with 13 other species of the family Muscidae. This mitogenome is 15,517?bp in length, with a standard set of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a typical non-coding control region. The gene organization, base composition, and secondary structures of tRNA genes of this mitogenome are similar to other dipteran species. Phylogenetic analyses based on all the mitochondrial PCGs and rRNA genes for 14 muscid and one outgroup taxa are performed using Bayesian inference (BI) and maximum likelihood (ML) analysis. The inferred trees indicate that the Muscidae, the subfamily Muscinae, and the tribes Azeliini and Reinwardtiini are monophyletic, whereas the monophyly of the subfamily Azeliinae and the tribe Muscini is not supported.     

The Complete Sequence of the Mitochondrial Genome of Butomus umbellatus – A Member of an Early Branching Lineage of Monocotyledons

In order to study the evolution of mitochondrial genomes in the early branching lineages of the monocotyledons, i.e., the Acorales and Alismatales, we are sequencing complete genomes from a suite of key taxa. As a starting point the present paper describes the mitochondrial genome of Butomus umbellatus (Butomaceae) based on next-generation sequencing data. The genome was assembled into a circular molecule, 450,826 bp in length. Coding sequences cover only 8.2% of the genome and include 28 protein coding genes, four rRNA genes, and 12 tRNA genes. Some of the tRNA genes and a 16S rRNA gene are transferred from the plastid genome. However, the total amount of recognized plastid sequences in the mitochondrial genome is only 1.5% and the amount of DNA transferred from the nucleus is also low. RNA editing is abundant and a total of 557 edited sites are predicted in the protein coding genes. Compared to the 40 angiosperm mitochondrial genomes sequenced to date, the GC content of the Butomus genome is uniquely high (49.1%). The overall similarity between the mitochondrial genomes of Butomus and Spirodela (Araceae), the closest relative yet sequenced, is low (less than 20%), and the two genomes differ in size by a factor 2. Gene order is also largely unconserved. However, based on its phylogenetic position within the core alismatids Butomus will serve as a good reference point for subsequent studies in the early branching lineages of the monocotyledons.     

Evaluation of a novel 16S rRNA/tRNA mitochondrial marker for the identification and phylogenetic analysis of shrimp species belonging to the superfamily Penaeoidea

In this study, we infer the phylogenetic relationships within commercial shrimp using sequence data from a novel mitochondrial marker consisting of an approximately 530-bp region of the 16S ribosomal RNA (rRNA)/transfer RNA (tRNA)^Val genes compared with two other mitochondrial genes: 16S rRNA and cytochrome c oxidase I (COI). All three mitochondrial markers were considerably AT rich, exhibiting values up to 78.2% for the species Penaeus monodon in the 16S rRNA/tRNA^Val genes, notably higher than the average among other Malacostracan mitochondrial genomes. Unlike the 16S rRNA and COI genes, the 16S rRNA/tRNA^Val marker evidenced that Parapenaeus is more closely related to Metapenaeus than to Solenocera, a result that seems to be more in agreement with the taxonomic status of these genera. To our knowledge, our study using the 16S rRNA/tRNA^Val gene as a marker for phylogenetic analysis offers the first genetic evidence to confirm that Pleoticus muelleri and Solenocera agassizi constitute a separate group and that they are more related to each other than to genera belonging to the family Penaeidae. The 16S rRNA/tRNA^Val region was also found to contain more variable sites (56%) than the other two regions studied (33.4% for the 16S rRNA region and 42.7% for the COI region). The presence of more variable sites in the 16S rRNA/tRNA^Val marker allowed the interspecific differentiation of all 19 species examined. This is especially useful at the commercial level for the identification of a large number of shrimp species, particularly when the lack of morphological characteristics prevents their differentiation.     

Characterization of the complete mitochondrial genome of Diploastrea heliopora and phylogeny of the scleractinia species which have group I introns in their COI genes

Mitochondrial genome DNA is a powerful marker for resolving phylogenetic relationships among scleractinian corals. Here, we decode the complete mitochondrial genome of Diploastrea heliopora (Lamarck, 1816) for the first time. The general features are 18 363 bp in length, and conventionally, with 13 protein coding genes, two ribosomal RNAs, and two transfer RNAs. Gene arrangement and distribution are similar to other scleractinian corals. Moreover, the COI gene of D. heliopora is broken up into two parts by a complex group I intron. This intron is 1076 bases in length and contains helical structures (P1-P10, except P2) and four conserved regions (P, Q, R, and S). The mitochondrial genome of D. heliopora has asymmetric base composition (13.03% C, 20.29% G, 25.91% A, and 40.77% for T). Based on concatenated protein coding genes, ML and BI trees show similar phylogenetic relationship: D. heliopora clustered closely with Sclerophyllia maxima and Echinophyllia aspera into the robust branch. The data and conclusion in this study are reference for further phylogenetic studies of corals.     

Karyotypic evolution in the Galliformes: an examination of the process of karyotypic evolution by comparison of the molecular cytogenetic findings with the molecular phylogeny

To define the process of karyotypic evolution in the Galliformes on a molecular basis, we conducted genome-wide comparative chromosome painting for eight species, i.e. silver pheasant (Lophura nycthemera), Lady Amherst's pheasant (Chrysolophus amherstiae), ring-necked pheasant (Phasianus colchicus), turkey (Meleagris gallopavo), Western capercaillie (Tetrao urogallus), Chinese bamboo-partridge (Bambusicola thoracica) and common peafowl (Pavo cristatus) of the Phasianidae, and plain chachalaca (Ortalis vetula) of the Cracidae, with chicken DNA probes of chromosomes 1-9 and Z. Including our previous data from five other species, chicken (Gallus gallus), Japanese quail (Coturnix japonica) and blue-breasted quail (Coturnix chinensis) of the Phasianidae, guinea fowl (Numida meleagris) of the Numididae and California quail (Callipepla californica) of the Odontophoridae, we represented the evolutionary changes of karyotypes in the 13 species of the Galliformes. In addition, we compared the cytogenetic data with the molecular phylogeny of the 13 species constructed with the nucleotide sequences of the mitochondrial cytochrome b gene, and discussed the process of karyotypic evolution in the Galliformes. Comparative chromosome painting confirmed the previous data on chromosome rearrangements obtained by G-banding analysis, and identified several novel chromosome rearrangements. The process of the evolutionary changes of macrochromosomes in the 13 species was in good accordance with the molecular phylogeny, and the ancestral karyotype of the Galliformes is represented.     

The complete mitochondrial genome of Eterusia aedea (Lepidoptera,Zygaenidae) and comparison with other zygaenid moths

Zygaenidae comprises >1036 species, including many folivorous pests in agriculture. In the present study, the complete mitochondrial genome (mitogenome) of a major pest of tea trees, Eterusia aedea was determined. The 15,196-bp circular genome contained the common set of 37 mitochondrial genes (including 13 protein-coding genes, two rRNA genes, and 22 tRNA genes) and exhibited the similar genomic features to reported Zygaenidae mitogenome. Comparative analyses of Zygaenidae mitogenomes showed a typical evolutionary trend of lepidopteran mitogenomes. In addition, we also investigated the gene order of lepidopteran mitogenomes and proposed that the novel gene order trnA-trnR-trnN-trnE-trnS-trnF from Zygaenidae and Gelechiidae and most other gene rearrangements of this tRNA cluster evolved independently. Finally, the mitogenomic phylogeny of Lepidoptera was reconstructed based on multiple mitochondrial datasets. And all the phylogenetic results revealed the sister relationships of Cossoidea and Zygaenoidea with both BI and ML methods, which is the first stable mitogenomic evidence for this clade.     

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 nucleotide sequence of the mitochondrial genome of <Emphasis Type="Italic">Phthonandria atrilineata</Emphasis> (Lepidoptera: Geometridae)

Using long-polymerase chain reaction (Long-PCR) method, we determined the complete nucleotide sequence of the mitochondrial genome (mitogenome) of Phthonandria atrilineata. The complete mtDNA from P. atrilineata was 15,499 base pairs in length and contained 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The P. atrilineata genes were in the same order and orientation as the completely sequenced mitogenomes of other lepidopteran species. The nucleotide composition of P. atrilineata mitogenome was biased toward A + T nucleotides (81.02%), and the 13 PCGs show different A + T contents that range from 73.25% (cox1) to 92.12% (atp8). Phthonandria had the canonical set of 22 tRNA genes, that fold in the typical cloverleaf structure described for metazoan mt tRNAs, with the unique exception of trnS(AGN). The phylogenetic relationships were reconstructed with the concatenated sequences of the 13 PCGs of the mitochondrial genome, which confirmed that P. atrilineata is most closely related to the superfamily Bombycoidea.     

The mitochondrial genome sequence analysis of Ophidascaris baylisi from the Burmese python (Python molurus bivittatus)

Ophidascaris species are parasitic roundworms that inhabit the python gut, resulting in severe granulomatous lesions or even death. However, the classification and nomenclature of these roundworms are still controversial. Our study aims to identify a snake roundworm from the Burmese python (Python molurus bivittatus) and analyze the mitochondrial genome. We identified this roundworm as Ophidascaris baylisi based on the morphology and cytochrome c oxidase subunit I (cox1) sequence. Ophidascaris baylisi complete mitochondrial genome was 14,784 bp in length, consisting of two non-coding regions and 36 mitochondrial genes (12 protein-coding genes, 22 tRNA genes, and two rRNA genes). The protein-coding genes used TTG, ATG, ATT, or TTA as start codons and TAG, TAA, or T as stop codons. All tRNA genes showed a TV-loop structure, except trnS1^AGN and trnS2^UCN revealed a D-loop structure. The mitochondrial large ribosomal subunit 16S (rrnL) and small ribosomal subunit 12S (rrnS) were 956 bp and 700 bp long, respectively. Phylogenetic analysis based on O. baylisi mitochondrial protein-coding genes demonstrated that O. baylisi clustered with the family Ascarididae members and was most closely related to Ophidascaris wangi. These results may enhance the nematode mitochondrial genome database and provide valuable molecular markers for further research on the taxonomy, phylogeny, and genetic relationships of Ophidascaris nematodes.     

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.     

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

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