Complete mitochondrial genome sequence of the Himalayan Griffon,Gyps himalayensis (Accipitriformes: Accipitridae): Sequence,structure, and phylogenetic analyses

This is the first study to describe the mitochondrial genome of the Himalayan Griffon, Gyps himalayensis, which is an Old World vulture belonging to the family Accipitridae and occurring along the Himalayas and the adjoining Tibetan Plateau. Its mitogenome is a closed circular molecule 17,381 bp in size containing 13 protein‐coding genes, 22 tRNA coding genes, two rRNA‐coding genes, a control region (CR), and an extra pseudo‐control region (CCR) that are conserved in most Accipitridae mitogenomes. The overall base composition of the G. himalayensis mitogenome is 24.55% A, 29.49% T, 31.59% C, and 14.37% G, which is typical for bird mitochondrial genomes. The alignment of the Accipitridae species control regions showed high levels of genetic variation and abundant AT content. At the 5′ end of the domain I region, a long continuous poly‐C sequence was found. Two tandem repeats were found in the pseudo‐control regions. Phylogenetic analysis with Bayesian inference and maximum likelihood based on 13 protein‐coding genes indicated that the relationships at the family level were (Falconidae + (Cathartidae + (Sagittariidae + (Accipitridae + Pandionidae))). In the Accipitridae clade, G. himalayensis is more closely related to Aegypius monachus than to Spilornis cheela. The complete mitogenome of G. himalayensis provides a potentially useful resource for further exploration of the taxonomic status and phylogenetic history of Gyps species.     

The first complete mitogenome of Picumnus innominatus (Aves,Piciformes, Picidae) and phylogenetic inference within the Picidae

The avian family Picidae, which is nearly global in distribution, contains the piculets (Picumninae and Nesoctitinae), the woodpeckers (Picinae), and the wrynecks (Jynginae). However, the phylogenetic relationships within the Picidae remain obscure for most genera. In the present study, the complete mitochondrial genome of Picumnus innominatus was determined and described, which was the first complete mitogenome reported in the Picumnus. The circular mitogenome of P. innominatus was 17,180 bp in size and consisted of 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNA genes, a control region (CR), and a noncoding region. The gene order and arrangement of the P. innominatus mitogenome were identical to other mitogenomes of the Picidae. Moreover, strikingly large tandem repeats were found in the noncoding region of the P. innominatus mitogenome, which have not yet been covered in other picid species to date. At the family level (Picidae), the highest dN/dS ratio was detected for the ND1 gene (1.38726) among 13 PCGs, indicating that positive selection was powerful for this gene. Bayesian and Maximum Likelihood phylogenetic analyses based on the combination of 12S rRNA and CYTB gene supported strongly that the Picumninae is monophyletic.     

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 Japanese sparrowhawk (Accipiter gularis) and the phylogenetic relationships among some predatory birds

Phylogenetic relationships among raptors, especially various groups are rather complex and controversial. We determined the complete mtDNA of Japanese sparrowhawk, and estimated phylogenetic trees based on the complete mtDNA alignment of it and 36 other raptor species, to clarify raptor phylogenetics. Phylogenetic trees were also estimated using a multiple sequence alignment of 12S rRNA and 16S rRNA from 81 typical species in GenBank, to further clarify the phylogenetic relationships of several groups among the raptors. The new mtDNA is a circular molecule, 17 917 bp in length, containing the 37 typical genes, with a pseudo-control region. ATG is generally the start codon, TAA is the most frequent stop codon. All tRNAs can be folded into canonical cloverleaf secondary structures except for tRNA^{Ser (AGY)} and tRNA^{Leu (CUN)}, which are missing the “DHU” arm. Phylogenetic relationships demonstrate that raptors can be divided into four branches: Accipitriformes, Falconiformes, Strigiformes and Caprimulgiformes in this study. We suggest that Accipitriformes should to be an independent order, Accipitriformes. The results also indicate that Accipitriformes contains three clades: Accipitridae, Pandionidae and Sagittariidae. Strigiformes includes species from Tytonidae and Strigidae. Caprimulgiformes contains Aegothelidae and Caprimulgidae.     

Complete mtDNA of Meretrix lusoria (Bivalvia: Veneridae) reveals the presence of an atp8 gene, length variation and heteroplasmy in the control region

The complete nucleotide sequence of the mitochondrial genome of the clam Meretrix lusoria (Bivalvia: Veneridae) was determined. It comprises 20,268 base pairs (bp) and contains 13 protein-coding genes, including ATPase subunit 8 (atp8), two ribosomal RNAs, 22 transfer RNAs, and a non-coding control region. The atp8 encodes a protein of 39 amino acids. All genes are encoded on the same strand. A putative control region (CR or D-loop) was identified in the major non-coding region (NCR) between the tRNA^Gly and tRNA^Gln. A 1087 bp tandem repeat fragment was identified that comprises nearly 11 copies of a 101 bp motif and accounts for approximately 41% of the NCR. The 101 bp tandem repeat motif of the NCR can be folded into a stem–loop secondary structure. Samples of eight individuals from Hainan and Fujian provinces were collected and their NCR regions were successfully amplified and sequenced. The data revealed a highly polymorphic VNTR (variable number of tandem repeats) associated with high levels of heteroplasmy in the D-loop region. The size of the CR ranged from 1942 to 3354 bp depending upon the copy number of the repeat sequence.     

Two novel mitogenomes of Dipodidae species and phylogeny of Rodentia inferred from the complete mitogenomes

Two novel mitogenomes of Eozapus setchuanus (KJ648495) and Sicista concolor (KJ648496) were reported and their total lengths were 16,630 bp and 16,493 bp, respectively. Both mitogenomes which were analogous to other rodent mitogenomes, contained 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and a control region. Specifically, the ND2 gene of S. concolor has three amino acids lesser than that of two other Dipodidae species (E. setchuanus and Jaculus jaculus) due to a premature termination codon in the 3′ end. We detected a tandem repeat cluster of 221 bp and 274 bp in the control region of S. concolor and E. setchuanus, respectively. Along with phylogenetic relationship analysis, we speculated that the tandem repeats in control regions might be common in Dipodinae species. Our phylogenetic analysis using concatenated mitochondrial gene datasets suggested five suborder and 16 family monophyletic groups in 54 rodent taxa sampled and strongly supported a basal position of the squirrel-related clade (PP = 1; BP = 100). Dipodidae had a sister-group relationship with Muroidea, and Sicistinae was in the base of Dipodidae clade. The complete mitochondrial genomes showed high resolution in deep-level phylogenetic relationship reconstructions of Rodentia.     

Complete Chloroplast Genome Sequence of Poisonous and Medicinal Plant Datura stramonium: Organizations and Implications for Genetic Engineering

Datura stramonium is a widely used poisonous plant with great medicinal and economic value. Its chloroplast (cp) genome is 155,871 bp in length with a typical quadripartite structure of the large (LSC, 86,302 bp) and small (SSC, 18,367 bp) single-copy regions, separated by a pair of inverted repeats (IRs, 25,601 bp). The genome contains 113 unique genes, including 80 protein-coding genes, 29 tRNAs and four rRNAs. A total of 11 forward, 9 palindromic and 13 tandem repeats were detected in the D. stramonium cp genome. Most simple sequence repeats (SSR) are AT-rich and are less abundant in coding regions than in non-coding regions. Both SSRs and GC content were unevenly distributed in the entire cp genome. All preferred synonymous codons were found to use A/T ending codons. The difference in GC contents of entire genomes and of the three-codon positions suggests that the D. stramonium cp genome might possess different genomic organization, in part due to different mutational pressures. The five most divergent coding regions and four non-coding regions (trnH-psbA, rps4-trnS, ndhD-ccsA, and ndhI-ndhG) were identified using whole plastome alignment, which can be used to develop molecular markers for phylogenetics and barcoding studies within the Solanaceae. Phylogenetic analysis based on 68 protein-coding genes supported Datura as a sister to Solanum. This study provides valuable information for phylogenetic and cp genetic engineering studies of this poisonous and medicinal plant.     

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.     

Complete mitochondrial genome of Oxuderces dentatus (Perciformes, Gobioidei)

The complete mitochondrial genome (mitogenome) of Oxuderces dentatus was determined first. The genome was 17,116?bp in length and consisted of 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA genes, and 2 main non-coding regions [the control region (CR) and the origin of the light strand replication], the gene composition and order of which was similar to most other vertebrates. The overall base composition of the heavy strand was T 27.9%, C 26.8%, A 30.2%, and G 15.1%, with a slight A+T bias of 58.1%. In addition to the discrete and conserved sequence blocks, unusual long tandem repeat unit (three 150-bp tandem repeat units and an incomplete copy of 146?bp) was also detected within CR. This mitogenome sequence data would play an important role in population genetics and phylogenetic analysis of the Gobioidei.     

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 Chloroplast Genome Sequence of the Medicinal Plant Salvia miltiorrhiza

Salvia miltiorrhiza is an important medicinal plant with great economic and medicinal value. The complete chloroplast (cp) genome sequence of Salvia miltiorrhiza, the first sequenced member of the Lamiaceae family, is reported here. The genome is 151,328 bp in length and exhibits a typical quadripartite structure of the large (LSC, 82,695 bp) and small (SSC, 17,555 bp) single-copy regions, separated by a pair of inverted repeats (IRs, 25,539 bp). It contains 114 unique genes, including 80 protein-coding genes, 30 tRNAs and four rRNAs. The genome structure, gene order, GC content and codon usage are similar to the typical angiosperm cp genomes. Four forward, three inverted and seven tandem repeats were detected in the Salvia miltiorrhiza cp genome. Simple sequence repeat (SSR) analysis among the 30 asterid cp genomes revealed that most SSRs are AT-rich, which contribute to the overall AT richness of these cp genomes. Additionally, fewer SSRs are distributed in the protein-coding sequences compared to the non-coding regions, indicating an uneven distribution of SSRs within the cp genomes. Entire cp genome comparison of Salvia miltiorrhiza and three other Lamiales cp genomes showed a high degree of sequence similarity and a relatively high divergence of intergenic spacers. Sequence divergence analysis discovered the ten most divergent and ten most conserved genes as well as their length variation, which will be helpful for phylogenetic studies in asterids. Our analysis also supports that both regional and functional constraints affect gene sequence evolution. Further, phylogenetic analysis demonstrated a sister relationship between Salvia miltiorrhiza and Sesamum indicum. The complete cp genome sequence of Salvia miltiorrhiza reported in this paper will facilitate population, phylogenetic and cp genetic engineering studies of this medicinal plant.     

The chloroplast genome sequence from <Emphasis Type="Italic">Eugenia uniflora</Emphasis>, a Myrtaceae from Neotropics

Eugenia uniflora is a plant native to tropical America that holds great ecological and economic importance. The complete chloroplast (cp) genome sequence of Eugenia uniflora, a member of the Neotropical Myrtaceae family, is reported here. The genome is 158,445 bp in length and exhibits a typical quadripartite structure of the large (LSC, 87,459 bp) and small (SSC, 18,318 bp) single-copy regions, separated by a pair of inverted repeats (IRs, 26,334 bp). It contains 111 unique genes, including 77 protein-coding genes, 30 tRNAs and 4 rRNAs. The genome structure, gene order, GC content and codon usage are similar to the typical angiosperm cp genomes. Comparison of the entire cp genomes of E. uniflora L. and three other Myrtaceae revealed an expansion of 43 bp in the intergenic spacer located between the IRA/large single-copy (LSC) border and the first gene of LSC region. Simple sequence repeat (SSR) analysis revealed that most SSRs are AT rich, which contribute to the overall AT richness of the cp genome. Additionally, fewer SSRs are distributed in the protein-coding sequences compared to the noncoding regions. Phylogenetic analysis among 58 species based on 57 cp genes demonstrated a closer relationship between E. uniflora L. and Syzygium cumini (L). Skeels compared to the Eucalyptus clade in the Myrtaceae family. The complete cp genome sequence of E. uniflora reported here has importance for population genetics, as well as phylogenetic and evolutionary studies in this species and other Myrtaceae species from Neotropical regions.     

Complete mitochondrial genome of Coelomactra antiquata (Mollusca: Bivalvia): The first representative from the family Mactridae with novel gene order and unusual tandem repeats

The complete mitochondrial genome plays an important role in the accurate inference of phylogenetic relationships among metazoans. Mactridae, also known as trough shells or duck clams, is an important family of marine bivalve clams in the order Veneroida. Here we present the complete mitochondrial genome sequence of the Xishishe Coelomactra antiquata (Mollusca: Bivalvia), which is the first representative from the family Mactridae. The mitochondrial genome of C. antiquata is of 17,384 bp in length, and encodes 35 genes, including 12 protein-coding, 21 transfer RNA, and 2 ribosomal RNA genes. Compared with the typical gene content of animal mitochondrial genomes, atp8 and tRNAS₂ are missing. Gene order of the mitochondrial genome of C. antiquata is unique compared with others from Veneroida. In the mitochondrial genome of the C. antiquata, a total of 2189 bp of non-coding nucleotides are scattered among 26 non-coding regions. The largest non-coding region contains one section of tandem repeats (99 bp × 11), which is the second largest tandem repeats found in the mitochondrial genomes from Veneroida. The phylogenetic trees based on mitochondrial genomes support the monophyly of Veneridae and Lucinidae, and the relationship at the family level: ((Veneridae + Mactridae) + (Cardiidae + Solecurtidae)) + Lucinidae. The phylogenetic result is consistent with the morphological classification. Meanwhile, bootstrap values are very high (BP = 94–100), suggesting that the evolutionary relationship based on mitochondrial genomes is very reliable.     

Characterization of the chloroplast genome sequence of oil palm (Elaeis guineensis Jacq.)

Oil palm (Elaeis guineensis Jacq.) is an economically important crop, which is grown for oil production. To better understand the molecular basis of oil palm chloroplasts, we characterized the complete chloroplast (cp) genome sequence obtained from 454 pyrosequencing. The oil palm cp genome is 156,973 bp in length consisting of a large single-copy region of?85,192 bp flanked on each side by inverted repeats of 27,071 bp with a small single-copy region of 17,639 bp joining the?repeats. The genome contains 112 unique genes: 79 protein-coding genes, 4 ribosomal RNA genes and 29 tRNA genes. By aligning the cp?genome sequence with oil palm cDNA sequences, we observed 18 non-silent and 10 silent RNA editing events among 19 cp protein-coding genes. Creation of an initiation codon by RNA editing in rpl2 has been reported in several monocots and was also found in the oil palm cp genome. Fifty common chloroplast protein-coding genes from 33 plant taxa were used to construct ML and MP?phylogenetic trees. Their topologies are similar and strongly support for the position of E. guineensis as the sister of closely related species Phoenix dactylifera in Arecaceae (palm families) of monocot subtrees.     

Complete mitochondrial genome of Trypauchen vagina (Perciformes, Gobioidei)

The complete mitochondrial genome of Trypauchen vagina was determined first. The genome is 16,686 bp in length and consists of 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA genes, and 2 main non-coding regions [the control region (CR) and the origin of the light strand replication], the gene composition and order of which was similar to most other vertebrates. The overall base composition of T. vagina is T 27.6%, C 27.6%, A 29.5%, and G 15.3%, with a slight A+T bias of 57.1%. In addition to the discrete and conserved sequence blocks, an incomplete tandem repeat unit is detected within the CR. This mitogenome sequence data would play an important role in population genetics and phylogenetic analysis of the Gobioidei.     

Complete mitochondrial genome of Odontobutis haifengensis (Perciformes,Odontobutiae): A unique rearrangement of tRNAs and additional non-coding regions identified in the genus Odontobutis

Herein, the complete mitochondrial genome of Odontobutis haifengensis was sequenced for the first time. The O. haifengensis mitogenome was 17,016 bp in length and included 13 protein-coding genes, 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region (CR). The genome organization, base composition, codon usage, and gene rearrangement was similar to other Odontobutis species. Furthermore, a tRNA gene rearrangement within the SLH cluster was found to be identical to other Odontobutis species. Moreover, the gene order and the positions of additional intergenic non-coding regions suggests that the observed unique gene rearrangement resulted from a tandem duplication and random loss of large-scale gene regions. Additionally, phylogenetic analysis showed that Odontobutis species form a monophyletic clade due to the conserved mitochondrial gene rearrangement. This study provides useful information that aids in a better understanding of mitogenomic diversity and evolutionary patterns of Odontobutidae species.     

Complete sequence of the chloroplast genome from pear (Pyrus pyrifolia): genome structure and comparative analysis

The chloroplast genome of Pyrus was found to be 159,922?bp in length which included a pair of inverted repeats (IRs) of 26,392?bp, separated by a small single-copy region of 19,237?bp and a large single-copy region (LSC) of 87,901?bp. A total of 130 predicted genes (113 unique genes and 17 genes, which were duplicated in the IR) including 79 protein-coding genes, four ribosomal RNA genes and 30 tRNA genes were identified based on similarity to homologs from the chloroplast genome of Nicotiana tabacum. Genome organization was very similar to the inferred ancestral angiosperm chloroplast genome. Comparisons between Pyrus, Malus, and Prunus in Rosaceae revealed 220 indels (??10?bp). Excluding ycf1 and ycf2, which contained deletions in the coding region, all of these were detected in the spacer or intron regions. Three insertions and 13 deletions were detected in Pyrus compared to the same loci in Malus and Prunus. After comparing 89 noncoding chloroplast DNA regions in Pyrus and Malus, highly variable regions such as ndhC-trnV and trnR-atpA were identified. In Pyrus and Malus, the IR/LSC borders were 62?bp shorter than those of Prunus. In addition, there were length mutations at the IRa/LSC junction and in trnH. A total of 67 simple sequence repeats (more than 10 repeated motifs) were identified in the Pyrus chloroplast genome. The indels and simple sequence repeats will be useful evolutionary tools at both intra- and interspecific levels. Phylogenetic analysis demonstrated a close relationship between Pyrus and Prunus in the Rosaceae.     

Complete chloroplast genome sequences of Dipterygium glaucum and Cleome chrysantha and other Cleomaceae Species,comparative analysis and phylogenetic relationships

This current study presents, for the first time, the complete chloroplast genome of two Cleomaceae species: Dipterygium glaucum and Cleome chrysantha in order to evaluate the evolutionary relationship. The cp genome is 158,576 bp in length with 35.74% GC content in D. glaucum and 158,111 bp with 35.96% GC in C. chrysantha. Inverted repeats IR 26,209 bp, 26,251 bp each, LSC of 87,738 bp, 87,184 bp and SSC of 18,420 bp, 18,425 bp respectively. There are 136 genes in the genome, which includes 80 protein coding genes, 31 tRNA genes and four rRNA genes were observed in both chloroplast genomes. 117 genes are unique while the remaining 19 genes are duplicated in IR regions. The analysis of repeats shows that the cp genome includes all types of repeats with more frequent occurrences of palindromic; Also, this analysis indicates that the total number of simple sequence repeats (SSR) were 323 in D. glaucum, and 313 in C. chrysantha, of which the majority of the SSRs in these plastid genomes were mononucleotide repeats A/T which are located in the intergenic spacer. Moreover, the comparative analysis of the four cp sequences revealed four hotspot genes (atpF, rpoC2, rps19, and ycf1), these variable regions could be used as molecular makers for the species authentication as well as resources for inferring phylogenetic relationships of the species. All the relationships in the phylogenetic tree are with high support, this indicate that the complete chloroplast genome is a useful data for inferring phylogenetic relationship within the Cleomaceae and other families. The simple sequence repeats identified will be useful for identification, genetic diversity, and other evolutionary studies of the species. This study reported the first cp genome of the genus Dipterygium and Cleome. The finding of this study will be beneficial for biological disciplines such as evolutionary and genetic diversity studies of the species within the core Cleomaceae.