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

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

Comparative Mitogenomic Analysis of Damsel Bugs Representing Three Tribes in the Family Nabidae (Insecta: Hemiptera)

Background

Nabidae, a family of predatory heteropterans, includes two subfamilies and five tribes. We previously reported the complete mitogenome of Alloeorhynchus bakeri, a representative of the tribe Prostemmatini in the subfamily Prostemmatinae. To gain a better understanding of architecture and evolution of mitogenome in Nabidae, mitogenomes of five species representing two tribes (Gorpini and Nabini) in the subfamily Nabinae were sequenced, and a comparative mitogenomic analysis of three nabid tribes in two subfamilies was carried out.

Methodology/Principal Findings

Nabid mitogenomes share a similar nucleotide composition and base bias, except for the control region, where differences are observed at the subfamily level. In addition, the pattern of codon usage is influenced by the GC content and consistent with the standard invertebrate mitochondrial genetic code and the preference for A+T-rich codons. The comparison among orthologous protein-coding genes shows that different genes have been subject to different rates of molecular evolution correlated with the GC content. The stems and anticodon loops of tRNAs are extremely conserved, and the nucleotide substitutions are largely restricted to TψC and DHU loops and extra arms, with insertion-deletion polymorphisms. Comparative analysis shows similar rates of substitution between the two rRNAs. Long non-coding regions are observed in most Gorpini and Nabini mtDNAs in-between trnI-trnQ and/or trnS2-nad1. The lone exception, Nabis apicalis, however, has lost three tRNAs. Overall, phylogenetic analysis using mitogenomic data is consistent with phylogenies constructed mainly form morphological traits.

Conclusions/Significance

This comparative mitogenomic analysis sheds light on the architecture and evolution of mitogenomes in the family Nabidae. Nucleotide diversity and mitogenomic traits are phylogenetically informative at subfamily level. Furthermore, inclusion of a broader range of samples representing various taxonomic levels is critical for the understanding of mitogenomic evolution in damsel bugs.     

Complex phylogeographic history of central African forest elephants and its implications for taxonomy

Background

The mitochondrial genomes of snakes are characterized by an overall evolutionary rate that appears to be one of the most accelerated among vertebrates. They also possess other unusual features, including short tRNAs and other genes, and a duplicated control region that has been stably maintained since it originated more than 70 million years ago. Here, we provide a detailed analysis of evolutionary dynamics in snake mitochondrial genomes to better understand the basis of these extreme characteristics, and to explore the relationship between mitochondrial genome molecular evolution, genome architecture, and molecular function. We sequenced complete mitochondrial genomes from Slowinski's corn snake (Pantherophis slowinskii) and two cottonmouths (Agkistrodon piscivorus) to complement previously existing mitochondrial genomes, and to provide an improved comparative view of how genome architecture affects molecular evolution at contrasting levels of divergence.

Results

We present a Bayesian genetic approach that suggests that the duplicated control region can function as an additional origin of heavy strand replication. The two control regions also appear to have different intra-specific versus inter-specific evolutionary dynamics that may be associated with complex modes of concerted evolution. We find that different genomic regions have experienced substantial accelerated evolution along early branches in snakes, with different genes having experienced dramatic accelerations along specific branches. Some of these accelerations appear to coincide with, or subsequent to, the shortening of various mitochondrial genes and the duplication of the control region and flanking tRNAs.

Conclusion

Fluctuations in the strength and pattern of selection during snake evolution have had widely varying gene-specific effects on substitution rates, and these rate accelerations may have been functionally related to unusual changes in genomic architecture. The among-lineage and among-gene variation in rate dynamics observed in snakes is the most extreme thus far observed in animal genomes, and provides an important study system for further evaluating the biochemical and physiological basis of evolutionary pressures in vertebrate mitochondria.     

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.     

Mitochondrial genome of Thais clavigera (Mollusca: Gastropoda): Affirmation of the conserved,ancestral gene pattern within the mollusks

Class Gastropoda includes a large number of described species, many with extensively rearranged mitochondrial genomes. We sequenced the mitogenome of the rock shell, Thais clavigera (Gastropoda: Muricidae), an intertidal snail, using long PCR with primers designed on the basis of expressed sequence tags. The mitogenome of T. clavigera consists of 2 rRNAs, 22 tRNAs, and 13 protein-coding genes, but no control region. Structural comparisons revealed that the order Sorbeoconcha, including T. clavigera, have nearly identical mitochondrial gene patterns. However, they have an inversion between a tRNA^Phe–tRNA^Glu cluster that comprises 21 genes, but most of the remaining structure is similar to the putative mollusk ground pattern. These findings will provide a better insight into mitochondrial gene rearrangement over the course of gastropod evolution.     

Evolution of selenium utilization traits

Background  

The essential trace element selenium is used in a wide variety of biological processes. Selenocysteine (Sec), the 21st amino acid, is co-translationally incorporated into a restricted set of proteins. It is encoded by an UGA codon with the help of tRNA^Sec (SelC), Sec-specific elongation factor (SelB) and a cis-acting mRNA structure (SECIS element). In addition, Sec synthase (SelA) and selenophosphate synthetase (SelD) are involved in the biosynthesis of Sec on the tRNA^Sec. Selenium is also found in the form of 2-selenouridine, a modified base present in the wobble position of certain tRNAs, whose synthesis is catalyzed by YbbB using selenophosphate as a precursor.     

The complete mitochondrial genome of Pardosa pusiola (Araneae,Lycosidae) and its phylogenetic implications

The mitochondrial genome (mitogenome) is useful for identification and phylogenetic analyses among arthropods, but there are no sufficient mitogenome data for wolf spiders. To enrich the mitogenome database of wolf spiders, the complete mitogenome of Pardosa pusiola was sequenced by high-throughput sequencing. It is 14,284 bp, comprising 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a control region (CR). It represents a high bias toward A and T nucleotides with an A + T content of 76.49%. The mitogenome exhibited a negative AT skew (−0.13) and a positive GC skew (0.32). Most PCGs started with ATN codons and ended with TAA, TAG, or an incomplete T. In addition, most tRNAs had aberrant secondary structures with the absence of DHU arm or TΨC arm. Analysis performed with CREx software demonstrated that large-scale rearrangements of tRNAs were observed in the mitogenome of P. pusiola as compared with the putative ancestral mitogenome. The Bayesian inference (BI) and maximum likelihood (ML) phylogenetic trees based on the 13 PCGs of 25 spiders had the same topology, which could be presented as (Araneidae + (Agelenidae + (Dictynidae + Desidae)) + (Salticidae + (Thomisidae + (Oxyopidae + (Pisauridae + Lycosidae))))). This study offers a useful genetic resource for the taxonomy and phylogeny of spiders.     

Characterization of Small RNAs Derived from tRNAs,rRNAs and snoRNAs and Their Response to Heat Stress in Wheat Seedlings

Small RNAs (sRNAs) derived from non-coding RNAs (ncRNAs), such as tRNAs, rRNAs and snoRNAs, have been identified in various organisms. Several observations have indicated that cleavage of tRNAs and rRNAs is induced by various stresses. To clarify whether sRNAs in wheat derived from tRNAs (stRNAs), rRNAs (srRNAs) and snoRNAs (sdRNAs) are produced specifically in association with heat stress responses, we carried out a bioinformatic analysis of sRNA libraries from wheat seedlings and performed comparisons between control and high-temperature-treated samples to measure the differential abundance of stRNAs, srRNAs and sdRNAs. We found that the production of sRNAs from tRNAs, 5.8S rRNAs, and 28S rRNAs was more specific than that from 5S rRNAs and 18S rRNAs, and more than 95% of the stRNAs were processed asymmetrically from the 3’ or 5’ ends of mature tRNAs. We identified 333 stRNAs and 8,822 srRNAs that were responsive to heat stress. Moreover, the expression of stRNAs derived from tRNA-Val-CAC, tRNA-Thr-UGU, tRNA-Tyr-GUA and tRNA-Ser-UGA was not only up-regulated under heat stress but also induced by osmotic stress, suggesting that the increased cleavage of tRNAs might be a mechanism that developed in wheat seedlings to help them cope with adverse environmental conditions.     

Pathology-related substitutions in human mitochondrial tRNA(Ile) reduce precursor 3' end processing efficiency in vitro

The human mitochondrial genome encodes 22 tRNAs interspersed among the two rRNAs and 11 mRNAs, often without spacers, suggesting that tRNAs must be efficiently excised. Numerous maternally transmitted diseases and syndromes arise from mutations in mitochondrial tRNAs, likely due to defect(s) in tRNA metabolism. We have systematically explored the effect of pathogenic mutations on tRNA^Ile precursor 3′ end maturation in vitro by 3′-tRNase. Strikingly, four pathogenic tRNA^Ile mutations reduce 3′-tRNase processing efficiency (V_max / K_M) to ~10-fold below that of wild-type, principally due to lower V_max. The structural impact of mutations was sought by secondary structure probing and wild-type tRNA^Ile precursor was found to fold into a canonical cloverleaf. Among the mutant tRNA^Ile precursors with the greatest 3′ end processing deficiencies, only G4309A displays a secondary structure substantially different from wild-type, with changes in the T domain proximal to the substitution. Reduced efficiency of tRNA^Ile precursor 3′ end processing, in one case associated with structural perturbations, could thus contribute to human mitochondrial diseases caused by mutant tRNAs.     

A manually curated database of tetrapod mitochondrially encoded tRNA sequences and secondary structures

Background  

Mitochondrial tRNAs have been the subject of study for structural biologists interested in their secondary structure characteristics, evolutionary biologists have researched patterns of compensatory and structural evolution and medical studies have been directed towards understanding the basis of human disease. However, an up to date, manually curated database of mitochondrially encoded tRNAs from higher animals is currently not available.     

Complete sequence and gene organization of the mitochondrial genome of scaly-sided merganser (<Emphasis Type="Italic">Mergus squamatus</Emphasis>) and phylogeny of some Anatidae species

The scaly-sided merganser (Mergus squamatus) is an endangered bird species on the IUCN Red List with the estimated global population of less than 2,500 individuals at present. In the present study, we studied the complete mitochondrial genome (mtDNA) and the phylogenetic of M. squamatus by PCR amplification and GenBank data. The genome was 16,595 bp in length and contained 37 genes (13 protein coding genes, two rRNAs, and 22 tRNAs) and a non-coding control region (D-loop). All protein-coding genes of M. squamatus mtDNA start with a typical ATG codon, except ND1, COI, and COII uses GTG as their initial codon. TAA, T- and TAG as the terminate codon occurred very commonly in the sequence. All tRNA genes can be folded into canonical cloverleaf secondary structure except for tRNA^Ser (AGY) and tRNA^Leu (CUN), which lose ‘‘DHU’’ arm. The genome sequences had been deposited in GenBank under accession number HQ833701. Based on the concatenated nucleotide sequences of mtDNA genes (Cyt b and D-loop), we reconstructed phylogenetic trees and discussed the phylogenetic relationships among ten Anatidae species. The results are different from the present classification, and we support Lophodytes cucullatus and Mergullus albellus to be members of the genus Mergus.     

Description of Aureibaculum luteum sp. nov. and Aureibaculum flavum sp. nov. isolated from Antarctic intertidal sediments

Two Gram-stain-negative, aerobic, non-motile, and rod-shaped bacterial strains, designated SM1352^T and A20^T, were isolated from intertidal sediments collected from King George Island, Antarctic. They shared 99.8% 16S rRNA gene sequence similarity with each other and had the highest sequence similarity of 98.1% to type strain of Aureibaculum marinum but?<?93.4% sequence similarity to those of other known bacterial species. The genomes of strains SM1352^T and A20^T consisted of 5,108,092 bp and 4,772,071 bp, respectively, with the G?+?C contents both being 32.0%. They respectively encoded 4360 (including 37 tRNAs and 6 rRNAs) and 4032 (including 36 tRNAs and 5 rRNAs) genes. In the phylogenetic trees based on 16S rRNA gene and single-copy orthologous clusters (OCs), both strains clustered with Aureibaculum marinum and together formed a separate branch within the family Flavobacteriaceae. The ANI and DDH values between the two strains and Aureibaculum marinum BH-SD17^T were all below the thresholds for species delineation. The major cellular fatty acids (>?10%) of the two strains included iso-C_15:0, iso-C_15:1 G, iso-C_17:0 3-OH. Their polar lipids predominantly included phosphatidylethanolamine, one unidentified aminophospholipid, one unidentified aminolipid, and two unidentified lipids. Genomic comparison revealed that both strains possessed much more glycoside hydrolases and sulfatase-rich polysaccharide utilization loci (PULs) than Aureibaculum marinum BH-SD17^T. Based on the above polyphasic evidences, strains SM1352^T and A20^T represent two novel species within the genus Aureibaculum, for which the names Aureibaculum luteum sp. nov. and Aureibaculum flavum sp. nov. are proposed. The type strains are SM1352^T (=?CCTCC AB 2014243 ^T?=?JCM 30335 ^T) and A20^T (=?CCTCC AB 2020370 ^T?=?KCTC 82503 ^T), respectively.

     

Complete mitochondrial genome sequence of the spinyhead croaker <Emphasis Type="Italic">Collichthys lucidus</Emphasis> (Perciformes,Sciaenidae) with phylogenetic considerations

The complete mitochondrial genome of the spinyhead croaker Collichthys lucidus was determined using long-PCR and primer walking methods. It is a circular molecule of 16,451 bp in length with a standard set of 22 tRNAs, 2 rRNAs, 13 protein-coding genes as well as a non-coding control region in the same order as those of the other bony fishes. C. lucidus mitogenome exhibited a clear strand-specific bias in nucleotide composition, as evidence by a GC- skew of the H-strand of −0.319. The strand-specific bias was also reflected in the codon usage of genes oriented in opposite directions. All tRNA genes except for tRNA ^Ser(AGY) harbored the typical cloverleaf secondary structures and possessed anticodons that matched the vertebrate mitochondrial genetic code. Contrary to the typical structure of control region consistig of TAS, central, and CSB domains, there were no central conserved blocks available in C. lucidus mitogenome. Despite extensive studies based on both morphology and molecules, phylogenetic position of C. lucidus with Sciaenidae is still controversial. Our phylogenetic results provided more evidence to support previous morphological studies and consistently placed C. lucidus as a sister taxon to Collichthys niveatus, with both of these taxa forming the monophyletic Collichthys.     

Three LIF-dependent signatures and gene clusters with atypical expression profiles, identified by transcriptome studies in mouse ES cells and early derivatives

Background

The apparent scarcity of available sequence data has greatly impeded evolutionary studies in Acari (mites and ticks). This subclass encompasses over 48,000 species and forms the largest group within the Arachnida. Although mitochondrial genomes are widely utilised for phylogenetic and population genetic studies, only 20 mitochondrial genomes of Acari have been determined, of which only one belongs to the diverse order of the Sarcoptiformes. In this study, we describe the mitochondrial genome of the European house dust mite Dermatophagoides pteronyssinus, the most important member of this largely neglected group.

Results

The mitochondrial genome of D. pteronyssinus is a circular DNA molecule of 14,203 bp. It contains the complete set of 37 genes (13 protein coding genes, 2 rRNA genes and 22 tRNA genes), usually present in metazoan mitochondrial genomes. The mitochondrial gene order differs considerably from that of other Acari mitochondrial genomes. Compared to the mitochondrial genome of Limulus polyphemus, considered as the ancestral arthropod pattern, only 11 of the 38 gene boundaries are conserved. The majority strand has a 72.6% AT-content but a GC-skew of 0.194. This skew is the reverse of that normally observed for typical animal mitochondrial genomes. A microsatellite was detected in a large non-coding region (286 bp), which probably functions as the control region. Almost all tRNA genes lack a T-arm, provoking the formation of canonical cloverleaf tRNA-structures, and both rRNA genes are considerably reduced in size. Finally, the genomic sequence was used to perform a phylogenetic study. Both maximum likelihood and Bayesian inference analysis clustered D. pteronyssinus with Steganacarus magnus, forming a sistergroup of the Trombidiformes.

Conclusion

Although the mitochondrial genome of D. pteronyssinus shares different features with previously characterised Acari mitochondrial genomes, it is unique in many ways. Gene order is extremely rearranged and represents a new pattern within the Acari. Both tRNAs and rRNAs are truncated, corroborating the theory of the functional co-evolution of these molecules. Furthermore, the strong and reversed GC- and AT-skews suggest the inversion of the control region as an evolutionary event. Finally, phylogenetic analysis using concatenated mt gene sequences succeeded in recovering Acari relationships concordant with traditional views of phylogeny of Acari.     

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.     

Evolutionary relatedness of mackerels of the genus Scomber based on complete mitochondrial genomes: Strong support to the recognition of Atlantic Scomber colias and Pacific Scomber japonicus as distinct species

Mackerels of the genus Scomber are commercially important species, but their taxonomic status is still controversial. Although previous phylogenetic data support the recognition of Atlantic Scomber colias and Pacific Scomber japonicus as separate species, it is only based on the analysis of partial mitochondrial and nuclear DNA sequences. In an attempt to shed light on this relevant issue, we have determined the complete mitochondrial DNA sequence of S. colias, S. japonicus, and Scomber australasicus. The total length of the mitogenomes was 16,568 bp for S. colias and 16,570 bp for both S. japonicus and S. australasicus. All mitogenomes had a gene content (13 protein-coding, 2 rRNAs, and 22 tRNAs) and organization similar to that observed in Scomber scombrus and most other vertebrates. The major noncoding region (control region) ranged between 865 and 866 bp in length and showed the typical conserved blocks. Phylogenetic analyses revealed a monophyletic origin of Scomber species with regard to other scombrid fish. The major finding of this study is that S. colias and S. japonicus were significantly grouped in distinct lineages within Scomber cluster, which phylogenetically constitutes evidence that they may be considered as separate species. Additionally, molecular data here presented provide a useful tool for evolutionary as well as population genetic studies.     

Next-generation sequencing,phylogenetic signal and comparative mitogenomic analyses in Metacrangonyctidae (Amphipoda: Crustacea)

Background

Comparative mitochondrial genomic analyses are rare among crustaceans below the family or genus level. The obliged subterranean crustacean amphipods of the family Metacrangonyctidae, found from the Hispaniola (Antilles) to the Middle East, including the Canary Islands and the peri-Mediterranean region, have an evolutionary history and peculiar biogeography that can respond to Tethyan vicariance. Indeed, recent phylogenetic analysis using all protein-coding mitochondrial sequences and one nuclear ribosomal gene have lent support to this hypothesis (Bauzà-Ribot et al. 2012).

Results

We present the analyses of mitochondrial genome sequences of 21 metacrangonyctids in the genera Metacrangonyx and Longipodacrangonyx, covering the entire geographical range of the family. Most mitogenomes were attained by next-generation sequencing techniques using long-PCR fragments sequenced by Roche FLX/454 or GS Junior pyro-sequencing, obtaining a coverage depth per nucleotide of up to 281×. All mitogenomes were AT-rich and included the usual 37 genes of the metazoan mitochondrial genome, but showed a unique derived gene order not matched in any other amphipod mitogenome. We compare and discuss features such as strand bias, phylogenetic informativeness, non-synonymous/synonymous substitution rates and other mitogenomic characteristics, including ribosomal and transfer RNAs annotation and structure.

Conclusions

Next-generation sequencing of pooled long-PCR amplicons can help to rapidly generate mitogenomic information of a high number of related species to be used in phylogenetic and genomic evolutionary studies. The mitogenomes of the Metacrangonyctidae have the usual characteristics of the metazoan mitogenomes (circular molecules of 15,000-16,000 bp, coding for 13 protein genes, 22 tRNAs and two ribosomal genes) and show a conserved gene order with several rearrangements with respect to the presumed Pancrustacean ground pattern. Strand nucleotide bias appears to be reversed with respect to the condition displayed in the majority of crustacean mitogenomes since metacrangonyctids show a GC-skew at the (+) and (-) strands; this feature has been reported also in the few mitogenomes of Isopoda (Peracarida) known thus far. The features of the rRNAs, tRNAs and sequence motifs of the control region of the Metacrangonyctidae are similar to those of the few crustaceans studied at present.

Electronic supplementary material

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