Complete mitochondrial genome sequence of Bonasa sewerzowi(Galliformes: Phasianidae) and phylogenetic analysis

Bonasa sewerzowi, the smallest and most southerly distributed grouse species in the world, is a bird endemic to China. The population of B. sewerzowi had shown a declining trend, which made it to be the endangered species in the China Red Data Book and Category I of nationally protected animals. So far, however, most studies about this species were mainly focused on the morphological and ecological aspects. In order to further study the feature of B. sewerzowi, the complete mitochondrial genome(mitogenome) of B. sewerzowi was sequenced by Illumina Hiseq 2000 high-throughput sequencing. Then, we focused on comparative genomics of two Bonasa species to find their characteristics. Finally, phylogenetic position of Bonasa was made based on the mitogenome dataset. Our results revealed that:(1) the mitogenome of B. sewerzowi, consisting of 16 658 bp, displayed typical genome organization and gene order found in other previously determined Galliformes mitogenomes;(2) the structure and composition of mitogenomes were similar between B. sewerzowi and B. bonasia;(3) the monophyly of Bonasa was well supported, which had a closer phylogenetic relationship with Meleagris gallopavo.     

Complete mitochondrial genome sequence of the Indian pipistrelle Pipistrellus coromandra (Vespertilioninae)

We characterized the complete mitogenome of Pipistrellus coromandra (Indian pipistrelle) for comparative analysis of mitogenomes and for resolving the phylogenetic relationship of four tribes in the subfamily Vespertilioninae. The mitogenome size of P. coromandra was 17,153?bp, with a control region and a typical set of 37 mitochondrial genes. The nucleotide composition of the P. coromandra mitogenome showed an AT bias with a nucleotide composition of 33.5% A, 30.7% T, 13.3% G, and 22.5% C. The mitochondrial protein-coding genes in P. coromandra use the standard start codon (ATN), two stop codons (TAA and AGA), and two incomplete stop codons (TA- and T--). The intertribal relationship of four tribes was highly resolved from the phylogenetic analysis of mitogenome sequences.     

Characterization of the complete mitochondrial genomes of Cnaphalocrocis medinalis and Chilo suppressalis (Lepidoptera: Pyralidae)

The complete mitochondrial genomes (mitogenomes) of Cnaphalocrocis medinalis and Chilo suppressalis (Lepidoptera: Pyralidae) were determined and analyzed. The circular genomes were 15,388 bp long for C. medinalis and 15,395 bp long for C. suppressalis. Both mitogenomes contained 37 genes, with gene order similar to that of other lepidopterans. Notably, 12 protein-coding genes (PCGs) utilized the standard ATN, but the cox1 gene used CGA as the initiation codon; the cox1, cox2, and nad4 genes in the two mitogenomes had the truncated termination codons T, T, and TA, respectively, but the nad5 gene was found to use T as the termination codon only in the C. medinalis mitogenome. Additionally, the codon distribution and Relative Synonymous Codon Usage of the 13 PCGs in the C. medinalis mitogenome were very different from those in other pyralid moth mitogenomes. Most of the tRNA genes had typical cloverleaf secondary structures. However, the dihydrouridine (DHU) arm of the trnS1(AGN) gene did not form a stable stem-loop structure. Forty-nine helices in six domains, and 33 helices in three domains were present in the secondary structures of the rrnL and rrnS genes of the two mitogenomes, respectively. There were four major intergenic spacers, except for the A+T-rich region, spanning at least 12 bp in the two mitogenomes. The A+T-rich region contained an 'ATAGT(A)'-like motif followed by a poly-T stretch in the two mitogenomes. In addition, there were a potential stem-loop structure, a duplicated 25-bp repeat element, and a microsatellite '(TA)(13)' observed in the A+T-rich region of the C. medinalis mitogenome. A poly-T motif, a duplicated 31-bp repeat element, and a 19-bp triplication were found in the C. suppressalis mitogenome. However, there are many differences in the A+T-rich regions between the C. suppressalis mitogenome sequence in the present study and previous reports. Finally, the phylogenetic relationships of these insects were reconstructed based on amino acid sequences of mitochondrial 13 PCGs using Bayesian inference and maximum likelihood methods. These molecular-based phylogenies support the traditional morphologically based view of relationships within the Pyralidae.     

Complete mitochondrial genomes of three skippers in the tribe Aeromachini (Lepidoptera: Hesperiidae: Hesperiinae) and their phylogenetic implications

The mitochondrial genome is now widely used in the study of phylogenetics and molecular evolution due to its maternal inheritance, fast evolutionary rate, and highly conserved gene content. To explore the phylogenetic relationships of the tribe Aeromachini within the subfamily Hesperiinae at the mitochondrial genomic level, we sequenced and annotated the complete mitogenomes of 3 skippers: Ampittia virgata, Halpe nephele, and Onryza maga (new mitogenomes for 2 genera) with a total length of 15,333 bp, 15,291 bp, and 15,381 bp, respectively. The mitogenomes all contain 13 protein‐coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a noncoding A + T‐rich region and are consistent with other lepidopterans in gene order and type. In addition, we reconstructed the phylogenetic trees of Hesperiinae using maximum likelihood (ML) and Bayesian inference (BI) methods based on mitogenomic data. Results show that the tribe Aeromachini in this study robustly constitute a monophyletic group in the subfamily Hesperiinae, with the relationships Coeliadinae + (Euschemoninae + (Pyrginae + ((Eudaminae + Tagiadinae) + (Heteropterinae + ((Trapezitinae + Barcinae) + Hesperiinae))))). Moreover, our study supports the view that Apostictopterus fuliginosus and Barca bicolor should be placed out of the subfamily Hesperiinae.     

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.     

The first complete mitochondrial genome of the Indian Tent Turtle,Pangshura tentoria (Testudines: Geoemydidae): Characterization and comparative analysis

The characterization of a complete mitogenome is widely used in genomics studies for systematics and evolutionary research. However, the sequences and structural motifs contained within the mitogenome of Testudines taxa have rarely been examined. The present study decodes the first complete mitochondrial genome of the Indian Tent Turtle, Pangshura tentoria (16,657 bp) by using next‐generation sequencing. This denovo assembly encodes 37 genes: 13 protein‐coding genes (PCGs), 22 transfer RNA (tRNAs), two ribosomal RNA, and one control region (CR). Most of the genes were encoded on majority strand, except for one PCG (NADH dehydrogenase subunit 6) and eight tRNAs. Most of the PCGs were started with an ATG initiation codon, except for Cytochrome oxidase subunit 1 with “GTG” and NADH dehydrogenase subunit 5 with “ATA.” The termination codons, “TAA” and “AGA” were observed in two subunits of NADH dehydrogenase gene. The relative synonymous codon usage analysis revealed the maximum abundance of alanine, isoleucine, leucine, and threonine. The nonsynonymous/synonymous ratios were <1 in all PCGs, which indicates strong negative selection among all Geoemydid species. The study also found the typical cloverleaf secondary structure in most of the tRNA genes, except for serine with the lack of the conventional DHU arm. The comparative study of Geoemydid mitogenomes revealed the occurrence of tandem repeats was frequent in the 3′ end of CR. Further, two copies of a unique tandem repeat “TTCTCTTT” were identified in P. tentoria. The Bayesian and maximum‐likelihood phylogenetic trees using concatenation of 13 PCGs revealed the close relationships of P. tentoria with Batagur trivittata in the studied dataset. All the Geoemydid species showed distinct clustering with high bootstrap support congruent with previous evolutionary hypotheses. We suggest that the generations of more mitogenomes of Geoemydid species are required, to improve our understanding of their in‐depth phylogenetic and evolutionary relationships.     

Complete mitochondrial genome of the Iberian Mole <Emphasis Type="Italic">Talpa occidentalis</Emphasis> (Talpidae,Insectivora) and comparison with <Emphasis Type="Italic">Talpa europaea</Emphasis>

The complete mitogenome of Talpa occidentalis, the Iberian mole, was sequenced using a combination of the Illumina and Sanger methods. The 16,962 bp genome obtained contains 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and a control region. Thirty-seven identical repetitions of a 10-nucleotide (CACACGTACG) repeat element were identified in the non-coding control region (D-loop). The number, order, and orientation of the mitochondrial genes are the same as in T. europaea, the only mitogenome published so far for this genus. These two mitogenomes differ only at the repeat element included in the control region. The phylogeny obtained for the Talpidae species using the protein-coding genes of these mitogenomes agrees with the current classification of this family.     

The mitogenome of freshwater loach Homatula laxiclathra (Teleostei: Nemacheilidae) with phylogenetic analysis of Nemacheilidae

The complete mitogenome can provide valuable genetic information to reconstruct relationships between species. In this study, we sequenced a stone loach, Homatula laxiclathra (Teleostei: Nemacheilidae), which is found in the northern region of the Qinling Mountains in China. The size of the H. laxiclathra mitogenome is 16,570 bp, which contains 37 typical mitochondrial genes including 13 protein‐coding genes, 22 transfer RNAs, two ribosomal RNAs, and a control region (D‐loop) with a total AT content of 55.8%. This is similar to other Nemacheilidae sequences published in GenBank. Furthermore, a mito‐phylogenomic analysis of 46 Nemacheilidae species places H. laxiclathra in a robust monophyletic Homatula cluster with other Homatula species. Our results contribute toward a better understanding of a true phylogeny of these species based on large‐scale taxonomic samplings as well as to help grasp the evolution of fish mitogenomes.     

The complete mitogenome of the Chinese bush cricket, Gampsocleis gratiosa （Orthoptera： Tettigonioidea）

The complete mitochondrial genome (mitogenome) of Gampsocleis gratiosa was determined. The 15, 929 bp in the size of G. gratiosa mitogenome contains a typical gene content, base composition, and codon usage found in metazoan. All 13 protein coding genes (PCGs) of the G. gratiosa mitogenome start with a typical ATN codon. The usual termination codons (TAA and TAG) were found from 10 PCGs. However, the atp6, nad4, and nad5 had incomplete termination codon (T). The anticodons of all tRNAs are identical to those observed in Drosophila yakuba and Locusta migratoria, and can be folded in the form of a typical clover leaf structure except for trnS (AGN). The secondary structure of trnS (AGN) was drawn according with the Steinberg-Cedergren tertiary structure. The A T content (67.4%) of the A T-rich region is relatively lower among the mitogenome regions, in contrast, it usually contains the highest A T content for most insects. Two isolated sequence repeat regions (202 bp) were found in the A T-rich region with mapping and secondary structure information.     

The compact mitochondrial genome of Zorotypus medoensis provides insights into phylogenetic position of Zoraptera

Background

Zoraptera, generally regarded as a member of Polyneoptera, represents one of the most enigmatic insect orders. Although phylogenetic analyses based on a wide array of morphological and/or nuclear data have been performed, the position of Zoraptera is still under debate. Mitochondrial genome (mitogenome) information is commonly considered to be preferable to reconstruct phylogenetic relationships, but no efforts have been made to incorporate it in Zorapteran phylogeny. To characterize Zoraptera mitogenome features and provide insights into its phylogenetic placement, here we sequenced, for the first time, one complete mitogenome of Zoraptera and reconstructed the phylogeny of Polyneoptera.

Results

The mitogenome of Zorotypus medoensis with an A + T content of 72.50% is composed of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a noncoding A + T-rich region. The gene content and arrangement are identical to those considered ancestral for insects. This mitogenome shows a number of very unusual features. First, it is very compact, comprising 14,572 bp, and is the smallest among all known polyneopteran mitogenomes. Second, both noncoding sequences and coding genes exhibit a significant decrease in size compared with those of other polyneopterans. Third, Z. medoensis mitogenome has experienced an accelerated substitution rate. Fourth, truncated secondary structures of tRNA genes occur with loss of dihydrouridine (DHU) arm in trnC, trnR, and trnS(AGN) and loss of TΨC arm in trnH and trnT. The phylogenetic analyses based on the mitogenome sequence information indicate that Zoraptera, represented by Z. medoensis, is recovered as sister to Embioptera. However, both Zoraptera and Embioptera exhibit very long branches in phylogenetic trees.

Conclusions

Characterization of Z. medoensis mitogenome contributes to our understanding of the enigmatic Zoraptera. Mitogenome data demonstrate an overall strong resolution of deep-level phylogenies of Polyneoptera but not Insecta. It is preferable to expand taxon sampling of Zoraptera and other poorly represented orders in future to break up long branches.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1156) contains supplementary material, which is available to authorized users.     

The Complete Mitochondrial Genome of two Tetragnatha Spiders (Araneae: Tetragnathidae): Severe Truncation of tRNAs and Novel Gene Rearrangements in Araneae

Mitogenomes can provide information for phylogenetic analysis and evolutionary biology. The Araneae is one of the largest orders of Arachnida with great economic importance. In order to develop mitogenome data for this significant group, we determined the complete mitogenomes of two long jawed spiders Tetragnatha maxillosa and T. nitens and performed the comparative analysis with previously published spider mitogenomes. The circular mitogenomes are 14578 bp long with A+T content of 74.5% in T. maxillosa and 14639 bp long with A+T content of 74.3% in T. nitens, respectively. Both the mitogenomes contain a standard set of 37 genes and an A+T-rich region with the same gene orientation as the other spider mitogenomes, with the exception of the different gene order by the rearrangement of two tRNAs (trnW and trnG). Most of the tRNAs lose TΨC arm stems and have unpaired amino acid acceptor arms. As interesting features, both trnS^AGN and trnS^UCN lack the dihydrouracil (DHU) arm and long tandem repeat units are presented in the A+T-rich region of both the spider mitogenomes. The phylogenetic relationships of 23 spider mitogenomes based on the concatenated nucleotides sequences of 13 protein-coding genes indicated that the mitogenome sequences could be useful in resolving higher-level relationship of Araneae. The molecular information acquired from the results of this study should be very useful for future researches on mitogenomic evolution and genetic diversities in spiders.     

Characterization of the complete mitochondrial genome of Diaphania pyloalis (Lepidoptera: Pyralididae)

The complete mitochondrial genome (mitogenome) of Diaphania pyloalis (Lepidoptera: Pyralididae) was determined to be 15,298 bp and has the typical gene organization of mitogenomes from lepidopteran insects. It consists of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. The A + T content of this mitogenome is 80.83% and the AT skew is slightly positive. All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which is initiated by CGA. Only the cox2 gene has an incomplete stop codon consisting of just a T. All the tRNA genes display a typical clover-leaf structure of mitochondrial tRNA. The A + T-rich region of the mitogenome is 332 bp in length, including several common features found in lepidopteran mitogenomes. Phylogenetic analysis showed that the D. pyloalis is close to Pyralididae.     

The time is ripe to investigate human centromeres by long-read sequencing

The complete sequencing of human centromeres, which are filled with highly repetitive elements, has long been challenging. In human centromeres, α-satellite monomers of about 171 bp in length are the basic repeating units, but α-satellite monomers constitute the higher-order repeat (HOR) units, and thousands of copies of highly homologous HOR units form large arrays, which have hampered sequence assembly of human centromeres. Because most HOR unit occurrences are covered by long reads of about 10 kb, the recent availability of much longer reads is expected to enable observation of individual HOR occurrences in terms of their single-nucleotide or structural variants. The time has come to examine the complete sequence of human centromeres.     

Mitogenomes provide new insights into the evolutionary history of Prodiamesinae (Diptera: Chironomidae)

Evolutionary analysis of Prodiamesinae has long been impeded by lack of information, and its phylogenetic relationship with Orthocladiinae remains questionable. Here, ten complete mitochondrial genomes (mitogenomes) of Orthocladiinae sensu lato were newly sequenced, including three Prodiamesinae species and seven Orthocladiinae species. Coupled with published mitogenomes, a total of 12 mitogenomes of Orthocladiinae sensu lato were selected for a comparative mitogenomic analysis and phylogenetic reconstruction. Mitogenomes of Orthocladiinae sensu lato are conserved in structure, and all genes arrange the same gene order as the ancestral insect mitogenome. Nucleotide composition is highly biased, and the control region displayed the highest A + T content. All protein-coding genes are under purifying selection, and the ATP8 evolves at the fastest rate. In addition, the mitogenomes of Orthocladiinae sensu lato are highly conserved, and they are practically useful for phylogenetic inference, suggesting a re-classification of Orthocladiinae by sinking Prodiamesinae as a subgroup of Orthocladiinae.     

肺形侧耳线粒体基因组多态性鉴定

线粒体是能量转化与ATP形成的重要场所,已有多种食用真菌的线粒体基因组被相继组装注释,但肺形侧耳的线粒体基因组鲜有报道。本研究对野生肺形侧耳X2菌株的线粒体基因组进行组装并注释,由野生菌株X2与主栽菌株JX线粒体大片段差异序列,构建分子标记来鉴别野生和主栽肺形侧耳菌株。结果显示X2的线粒体基因组为大小75 709bp的闭合环状结构,含有rRNA的大小亚基基因,25个负责携带氨基酸的tRNA基因以及14个常见蛋白编码基因,在cox1基因上含有9个内含子。内含子主要为IB型,包含LAGLIDADG_1 superfamily、GIY-YIG_Cterm superfamily保守结构域。位于菌株X2与JX内含子和基因间区上的大片段差异序列,是造成种内线粒体多态性的主要因素。根据两菌株的差异片段构建系统进化树的结果显示,内含子、rnl-trnX、trnD-atp6序列能够将两株野生菌株和其他主栽菌株区分开,可用于野生肺形侧耳种内鉴定,而且trnD-atp6效果最佳。     

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.