The first complete mitochondrial genome sequence of Nanorana parkeri and Nanorana ventripunctata (Amphibia: Anura: Dicroglossidae), with related phylogenetic analyses

Members of the Nanorana genus (family Dicroglossidae) are often referred to as excellent model species with which to study amphibian adaptations to extreme environments and also as excellent keystone taxa for providing insights into the evolution of the Dicroglossidae. However, a complete mitochondrial genome is currently only available for Nanorana pleskei. Thus, we analyzed the complete mitochondrial genomes of Nanorana parkeri and Nanorana ventripunctata to investigate their evolutionary relationships within Nanorana and their phylogenetic position in the family Dicroglossidae. Our results showed that the genomes of N. parkeri (17,837 bp) and N. ventripunctata (18,373 bp) encode 13 protein‐coding genes (PCGs), two ribosomal RNA genes, 23 transfer RNA (tRNA) genes, and a noncoding control region. Overall sequences and genome structure of the two species showed high degree of similarity with N. pleskei, although the motif structures and repeat sequences of the putative control region showed clear differences among these three Nanorana species. In addition, a tandem repeat of the tRNA‐Met gene was found located between the tRNA‐Gln and ND2 genes. On both the 5′ and 3′‐sides, the control region possessed distinct repeat regions; however, the CSB‐2 motif was not found in N. pleskei. Based on the nucleotide sequences of 13 PCGs, our phylogenetic analyses, using Bayesian inference and maximum‐likelihood methods, illustrate the taxonomic status of Nanorana with robust support showing that N. ventripunctata and N. pleskei are more closely related than they are to N. parkeri. In conclusion, our analyses provide a more robust and reliable perspective on the evolutionary history of Dicroglossidae than earlier analyses, which used only a single species (N. pleskei).     

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 mitochondrial genome of the plant bug Apolygus lucorum （Hemiptera： Miridae）＂ Presently known as the smallest in Heteroptera

The complete mitochondrial （mt） genome of the plant bug, Apolygus lucorum, an important cotton pest, has been sequenced and annotated in this study. The entire circular genome is 14 768 bp in size and represents the smallest in presently known heteropteran mt genomes. The mt genome is encoding for two ribosomal RNA （rRNA） genes, 22 transfer RNA （tRNA） genes, 13 protein coding genes and a control region, and the order, content, codon usage and base organization show similarity to a great extent to the hypothetical ancestral model. All protein coding genes use standard initiation codons ATN. Conventional stop codons TAA and TAG have been assigned to the most protein coding genes; however, COIII, ND4 and ND5 genes show incomplete terminator signal （T）. All tRNA genes possess the typical clover leaf structure, but the dihydrouridine arm of tRNAser（A6N） only forms a simple loop. Secondary structure models of rRNA genes are generally in accordance with the former models, although some differences exist in certain parts. Three intergenic spacers have never been found in sequenced mt genomes of Heteroptera. The phylogenetic study based on protein coding genes is largely congruent with previous phylogenetic work. Both Bayesian inference and maximum likelihood analyses highly support the sister relationship ofA. lucorum and Lygus lineolaris, and Miridae presents a sister position to Anthocoridae.     

Divergence and species discrimination in freshwater bryozoans (Bryozoa: Phylactolaemata)

We explored the suitability of nuclear and mitochondrial ribosomal markers [small subunit nuclear ribosomal RNA gene, large subunit nuclear ribosomal RNA gene, and a region spanning partial small mitochondrial ribosomal RNA subunit, four transfer RNA genes, and partial large mitochondrial ribosomal RNA subunit (referred to as rrnS‐rrnL)] for resolving patterns of diversification of 27 freshwater bryozoan species (class: Phylactolaemata) and evaluated the utility of statoblast ultrastructural features and molecular phylogenies for species discrimination in the Fredericellidae and Plumatellidae. Molecular data identified Plumatella fruticosa as distinct from the rest of the plumatellids, rendering the latter polyphyletic. rrnS‐rrnL was the most suitable marker for species discrimination and identified two undescribed species of Plumatella and at least two undescribed species of Fredericella. Lack of wide dispersal by fredericellid statoblasts may underlie the observed propensity for cryptic speciation and phylogeographical structure in Fredericella. Conversely, the strong dispersal potential of plumatellid statoblasts may mediate efficient gene flow between distant populations and explain the relatively low intraspecific divergence and lack of evidence for cryptic speciation. We show that species identification based on external features of statoblasts can be problematic in both genera, including for a putatively highly invasive, biofouling species, Plumatella vaihiriae, thereby highlighting the utility of rrnS‐rrnL sequences for species barcoding. © 2013 The Linnean Society of London     

Species‐specific mitochondrial gene rearrangements in biting midges and vector species identification

Abstract Partial mitochondrial gene sequences of 16 Culicoides species were determined to elucidate phylogenetic relations among species and to develop a molecular identification method for important virus vector species. In addition, the analysis found mitochondrial gene rearrangement in several species. Sequences of the mitochondrial genome region, cox1–trnL2–cox2 (1940–3785 bp) of 16 Culicoides and additional sequences were determined in some species, including whole mitochondrial genome sequences of Culicoides arakawae. Nine species showed common organization in this region, with three genes cox1–trnL2–cox2 and a small or no intergenic region (0–30 bp) between them. The other seven species showed translocation of tRNA and protein‐coding genes and/or insertion of AT‐rich non‐coding sequences (65–1846 bp) between the genes. The varied gene rearrangements among species within a genus is very rare for mitochondrial genome organization. Phylogenetic analyses based on the sequences of cox1+cox2 suggest a few clades among Japanese Culicoides species. No relationships between phylogenetic closeness and mitochondrial gene rearrangements were observed. Sequence data were used to establish a polymerase chain reaction tool to distinguish three important vector species from other Culicoides species, for which classification during larval stages is not advanced and identification is difficult.     

Complete mitochondrial genome of a forensically important muscid, Hydrotaea chalcogaster (Diptera: Muscidae), with notes on its phylogenetic position

Muscidae are a dipteran family which is important for forensic investigations. However, it has received limited attention in forensic entomological experiments as a reason of identification issues. It is hard to identify specimens by morphological methods, especially in developmental stages. Therefore, complete mitochondrial genome sequences can be important tool in forensic entomology for identifying species. In this study we sequenced and analyzed the first complete mitochondrial genome from a forensically important Muscidae species Hydrotaea (=Ophyra) chalcogaster by next-generation sequencing. The mitochondrial genome of the sequenced species is circular molecules of 17,076?bp which have the typical mitochondrial genome complement of 13 protein-coding genes, 22 tRNAs, two ribosomal RNA genes and a control region. Rearrangements of gene positions are identical with the ancestral insect genome. Furthermore, phylogenetic relationships of the family Muscidae were evaluated in regard to mitochondrial protein coding genes. The inferred trees indicate that the Muscidae is a paraphyletic family. These data provide additional information for molecular identification of muscid species.     

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.     

Evolution of mushroom mitochondrial DNA:Suillus and related genera

Summary Mapping studies were performed with 18 cloned probes on mitochondrial DNA (mtDNA) from 15 species ofSuillus and four species from three related genera of fleshy pore mushrooms (Boletaceae). WithinSuillus, mtDNAs vary in size from 36 to 121 kb, differ in gene order by only one major rearrangement, and have diverged in nucleotide sequence within the large subunit ribosomal RNA gene region by up to 2.9%. Three additional gene orders exist in related genera. Two of the three can be transformed into the predominantSuillus order by either one or two rearrangements. The fourth requires two to three rearrangements to be converted to any of the others. The minimum estimates of nucleotide divergence within the large subunit ribosomal RNA gene region vary from 8.3% to 11% in comparisons betweenSuillus and these related species. Trees based on restriction-site and size differences within the mitochondrial ribosomal RNA genes were consistent with the hypothesized sequence of genome rearrangements and provide suggestive evidence for a major expansion of the mitochondrial genome withinSuillus. Structural and sequence changes in mtDNA provided information about phylogenetic relationships within the Boletaceae.     

Analysis of the mitochondrial genome of Schistocerca gregaria gregaria (Orthoptera: Acrididae)

In the present study, we present the full sequence of the mitochondrial genome of the African desert locust Schistocerca gregaria gregaria. The size of 15625 bp reported matches very well with mitochondrial genomes of other Orthopteriodea. The mitochondrial genome comprises 13 protein‐coding genes, two ribosomal RNAs and 22 t‐RNAs with two t‐RNA (trnD and trnK) rearrangements that are typical for the taxon Caelifera. We compared the sequence with 12 mitochondrial genes of Schistocerca gregaria flaviventris and Schistocerca americana and used some of these data to construct phylogenetic trees, which confirm the close relationship between the two subspecies S. g. flaviventris and S. g. gregaria. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 296–305.     

红足壮异蝽Urochela quadrinotata Reuter线粒体基因组的测定与分析

研究测定并分析了红足壮异蝽Urochela quadrinotata Reuter的线粒体基因组全序列。该线粒体基因组全长16585bp(GenBank登录号为JQ743678),A+T含量为75.4%,共编码35个基因,包括13个蛋白质基因、20个tRNA基因(两个tRNA基因,即tRNAIle和tRNAGln未被检测到)、2个rRNA基因及一段较长的非编码区(控制区,亦称A+T富含区)。基因排序与大部分昆虫的线粒体基因排列方式相同,没有发生基因重排。除tRNASer(AGN)的DHU臂无法形成典型的茎环结构,其余tRNA基因均能稳定形成典型的三叶草二级结构。预测了红足壮异蝽16S和12S rRNA的二级结果,分别包括6个结构域43个茎环和3个结构域27茎环。控制区含一个长1652bp的串联重复区域,由16个串联重复单元组成。     

The mitochondrial genome of a sea anemone Bolocera sp. exhibits novel genetic structures potentially involved in adaptation to the deep‐sea environment

The deep sea is one of the most extensive ecosystems on earth. Organisms living there survive in an extremely harsh environment, and their mitochondrial energy metabolism might be a result of evolution. As one of the most important organelles, mitochondria generate energy through energy metabolism and play an important role in almost all biological activities. In this study, the mitogenome of a deep‐sea sea anemone (Bolocera sp.) was sequenced and characterized. Like other metazoans, it contained 13 energy pathway protein‐coding genes and two ribosomal RNAs. However, it also exhibited some unique features: just two transfer RNA genes, two group I introns, two transposon‐like noncanonical open reading frames (ORFs), and a control region‐like (CR‐like) element. All of the mitochondrial genes were coded by the same strand (the H‐strand). The genetic order and orientation were identical to those of most sequenced actiniarians. Phylogenetic analyses showed that this species was closely related to Bolocera tuediae. Positive selection analysis showed that three residues (31 L and 42 N in ATP6, 570 S in ND5) of Bolocera sp. were positively selected sites. By comparing these features with those of shallow sea anemone species, we deduced that these novel gene features may influence the activity of mitochondrial genes. This study may provide some clues regarding the adaptation of Bolocera sp. to the deep‐sea environment.     

Plant mitochondrial rps2 genes code for proteins with a C-terminal extension that is processed

A gene (rps2) coding for ribosomal protein S2 (RPS2) is present in the mitochondrial (mt) genome of several monocot plants, but absent from the mtDNA of dicots. Confirming that in dicot plants the corresponding gene has been transferred to the nucleus, a corresponding Arabidopsis thaliana nuclear gene was identified that codes for mitochondrial RPS2. As several yeast and mammalian genes coding for mt ribosomal proteins, the Arabidopsis RPS2 apparently has no N-terminal targeting sequence. In the maize mt genome, two rps2 genes were identified and both are transcribed, although at different levels. As in wheat and rice, the maize genes code for proteins with long C-terminal extensions, as compared to their bacterial counterparts. These extensions are not conserved in sequence. Using specific antibodies against one of the maize proteins we found that a large protein precursor is indeed synthesized, but it is apparently processed to give the mature RPS2 protein which is associated with the mitochondrial ribosome.     

Complete DNA sequence of the mitochondrial genome ofCepaea nemoralis (Gastropoda: Pulmonata)

The nucleotide sequence of a mitochondrial genome of the pulmonate gastropod molluscCepaea nemoralis has been determined. Contained within the 14,100 basepairs (bp) are the two ribosomal RNA genes and 13 protein coding genes typical of metazoan mitochondrial genomes. TheCepaea mtDNA does contain a gene for ATPase subunit 8, like the clausiliid pulmonate,Albinaria, and the chiton,Katharina, but unlike the bivalve mollusc,Mytilus. The mitochondrial genetic code ofCepaea is proposed to be the same as that ofMytilus, Katharina, andDrosophila. Only 14 putative tRNA genes are presented, although there is sufficient unassigned sequence to encode the remainder of the expected total of 22 tRNA genes. These 14 tRNA genes are a mixture of standard cloverleaf structures and nonstandard structures containing TV replacement loops as seen in nematode and mosquito mitochondrial genomes. If the eight unidentified tRNA genes are indeed present, very little unassigned sequence would remain to serve as a control region. Genes are transcribed from both strands of the molecule. Base composition is the least biased for any reported animal mitochondrial genome and is also very little skewed between strands using measures independent of base composition. TheCepaea mitochondrial gene order is quite unlike that of any other reported metazoan mtDNA, with the exception of the recently reported partial sequences ofAlbinaria. No gene bound-aries are shared among all the reported molluscan taxa, demonstrating a complete lack of conservation of mitochondrial gene order across the phylum Mollusca.     

Organization and closing of mitochondrial deoxyribonucleic acid from Paramecium tetraaurelia and Paramecium primaurelia.

Previously we showed that the mitochondrial deoxyribonucleic acid (DNA) from Paramecium aurelia consists of a linear genome and that replication of this genome is initiated at one terminus and proceeds unidirectionally to the other terminus. Analyses of mitochondria from four closely related species (1, 4, 5, and 7) indicated that the species 1, 5, and 7 DNAs are essentially completely homologous but that the species 4 mitochondrial DNA is only 40 to 50% homologous with that from species 1. The major regions of homology are those containing the genes for ribosomal ribonucleic acid (RNA). To understand the replication and organization of the linear mitochondrial genome better, we compared species 1 (Paramecium primaurelia) and 4 (Paramecium tetraaurelia) DNAs with regard to restriction fragment mapping and homology between initiation regions; we also identified the sites of the genes for ribosomal RNA. In general, the structures of the species 1 and 4 mitochondrial genomes were quite similar. Each ribosomal RNA gene was present in one copy per genome, with the large ribosomal RNA gene located near the terminal region of replication and the small ribosomal RNA gene located more centrally. These two genes were separated by about 10 kilobases in the species 1 genome and by about 12 kilobases in the species 4 genome. In contrast to our previous findings, by using nonstringent hybridization conditions we detected homology between the species 1 and 4 DNA fragments containing the initiation regions. We constructed recombinant DNA clones for many fragments, especially those containing the initiation region and the ribosomal RNA genes. We also constructed restriction enzyme maps for six enzymes for both P. primaurelia and P. tetraaurelia.     

Genome streamlining via complete loss of introns has occurred multiple times in lichenized fungal mitochondria

Reductions in genome size and complexity are a hallmark of obligate symbioses. The mitochondrial genome displays clear examples of these reductions, with the ancestral alpha‐proteobacterial genome size and gene number having been reduced by orders of magnitude in most descendent modern mitochondrial genomes. Here, we examine patterns of mitochondrial evolution specifically looking at intron size, number, and position across 58 species from 21 genera of lichenized Ascomycete fungi, representing a broad range of fungal diversity and niches. Our results show that the cox1gene always contained the highest number of introns out of all the mitochondrial protein‐coding genes, that high intron sequence similarity (>90%) can be maintained between different genera, and that lichens have undergone at least two instances of complete, genome‐wide intron loss consistent with evidence for genome streamlining via loss of parasitic, noncoding DNA, in Phlyctis boliviensisand Graphis lineola. Notably, however, lichenized fungi have not only undergone intron loss but in some instances have expanded considerably in size due to intron proliferation (e.g., Alectoria fallacina and Parmotrema neotropicum), even between closely related sister species (e.g., Cladonia). These results shed light on the highly dynamic mitochondrial evolution that is occurring in lichens and suggest that these obligate symbiotic organisms are in some cases undergoing recent, broad‐scale genome streamlining via loss of protein‐coding genes as well as noncoding, parasitic DNA elements.     

The complete chloroplast genome of Gracilariopsis lemaneiformis (Rhodophyta) gives new insight into the evolution of family Gracilariaceae

The complete chloroplast genome of Gracilariopsis lemaneiformis was recovered from a Next Generation Sequencing data set. Without quadripartite structure, this chloroplast genome (183,013 bp, 27.40% GC content) contains 202 protein‐coding genes, 34 tRNA genes, 3 rRNA genes, and 1 tmRNA gene. Synteny analysis showed plasmid incorporation regions in chloroplast genomes of three species of family Gracilariaceae and in Grateloupia taiwanensis of family Halymeniaceae. Combined with reported red algal plasmid sequences in nuclear and mitochondrial genomes, we postulated that red algal plasmids may have played an important role in ancient horizontal gene transfer among nuclear, chloroplast, and mitochondrial genomes. Substitution rate analysis showed that purifying selective forces maintaining stability of protein‐coding genes of nine red algal chloroplast genomes over long periods must be strong and that the forces acting on gene groups and single genes of nine red algal chloroplast genomes were similar and consistent. The divergence of Gp. lemaneiformis occurred ~447.98 million years ago (Mya), close to the divergence time of genus Pyropia and Porphyra (443.62 Mya).     

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.     

The mitochondrial genome of the water spider Argyroneta aquatica (Araneae: Cybaeidae)

We sequenced nearly the entire mitochondrial genome of Argyroneta aquatica, a wholly underwater‐living spider, thereby enhancing the available genomic information for Arachnida. The confirmed sequences contained the complete set of known genes present in other metazoan mitochondrial genomes. However, the mitochondrial gene order of A. aquatica was distinctly different from that of the most distant Chelicerata Limulus polyphemus (Xiphosura), probably because of a series of gene translocations and/or inversions. Comparison of arachnid mitochondrial gene orders for the purpose of phylogenetic inference is only minimally useful, but provides a strong signal in closely related lineages. To test the basal relationships and the evolutionary pattern of tRNA gene rearrangements among Arachnida, phylogenetic analyses using amino acid sequences of the 13 protein‐coding genes were performed. An interesting feature, the five 135‐bp tandem repeats and two 363‐bp tandem repeats, was identified in the putative control region. Although control region tandem repeats have been reported in many other arachnid and metazoan species, this is the first time it has been described in spiders.