The Complete Mitochondrial Genome of Gossypium hirsutum and Evolutionary Analysis of Higher Plant Mitochondrial Genomes

Background

Mitochondria are the main manufacturers of cellular ATP in eukaryotes. The plant mitochondrial genome contains large number of foreign DNA and repeated sequences undergone frequently intramolecular recombination. Upland Cotton (Gossypium hirsutum L.) is one of the main natural fiber crops and also an important oil-producing plant in the world. Sequencing of the cotton mitochondrial (mt) genome could be helpful for the evolution research of plant mt genomes.

Methodology/Principal Findings

We utilized 454 technology for sequencing and combined with Fosmid library of the Gossypium hirsutum mt genome screening and positive clones sequencing and conducted a series of evolutionary analysis on Cycas taitungensis and 24 angiosperms mt genomes. After data assembling and contigs joining, the complete mitochondrial genome sequence of G. hirsutum was obtained. The completed G.hirsutum mt genome is 621,884 bp in length, and contained 68 genes, including 35 protein genes, four rRNA genes and 29 tRNA genes. Five gene clusters are found conserved in all plant mt genomes; one and four clusters are specifically conserved in monocots and dicots, respectively. Homologous sequences are distributed along the plant mt genomes and species closely related share the most homologous sequences. For species that have both mt and chloroplast genome sequences available, we checked the location of cp-like migration and found several fragments closely linked with mitochondrial genes.

Conclusion

The G. hirsutum mt genome possesses most of the common characters of higher plant mt genomes. The existence of syntenic gene clusters, as well as the conservation of some intergenic sequences and genic content among the plant mt genomes suggest that evolution of mt genomes is consistent with plant taxonomy but independent among different species.     

Decoding and analysis of organelle genomes of Indian tea (Camellia assamica) for phylogenetic confirmation

The NCBI database has >15 chloroplast (cp) genome sequences available for different Camellia species but none for C. assamica. There is no report of any mitochondrial (mt) genome in the Camellia genus or Theaceae family. With the strong believes that these organelle genomes can play a great tool for taxonomic and phylogenetic analysis, we successfully assembled and analyzed cp and mt genome of C. assamica. We assembled the complete mt genome of C. assamica in a single circular contig of 707,441 bp length comprising of a total of 66 annotated genes, including 35 protein-coding genes, 29 tRNAs and two rRNAs. The first ever cp genome of C. assamica resulted in a circular contig of 157,353 bp length with a typical quadripartite structure. Phylogenetic analysis based on these organelle genomes showed that C. assamica was closely related to C. sinensis and C. leptophylla. It also supports Caryophyllales as Superasterids.     

Comparative Mitochondrial Genomics within and among Yeast Species of the Lachancea Genus

Yeasts are leading model organisms for mitochondrial genome studies. The explosion of complete sequence of yeast mitochondrial (mt) genomes revealed a wide diversity of organization and structure between species. Recently, genome-wide polymorphism survey on the mt genome of isolates of a single species, Lachancea kluyveri, was also performed. To compare the mitochondrial genome evolution at two hierarchical levels: within and among closely related species, we focused on five species of the Lachancea genus, which are close relatives of L. kluyveri. Hence, we sequenced the complete mt genome of L. dasiensis, L. nothofagi, L. mirantina, L. fantastica and L. meyersii. The phylogeny of the Lachancea genus was explored using these data. Analysis of intra- and interspecific variability across the whole Lachancea genus led to the same conclusions regarding the mitochondrial genome evolution. These genomes exhibit a similar architecture and are completely syntenic. Nevertheless, genome sizes vary considerably because of the variations of the intergenic regions and the intron content, contributing to mitochondrial genome plasticity. The high variability of the intergenic regions stands in contrast to the high level of similarity of protein sequences. Quantification of the selective constraints clearly revealed that most of the mitochondrial genes are under purifying selection in the whole genus.     

Complete Mitochondrial Genome of Echinostoma hortense (Digenea: Echinostomatidae)

Echinostoma hortense (Digenea: Echinostomatidae) is one of the intestinal flukes with medical importance in humans. However, the mitochondrial (mt) genome of this fluke has not been known yet. The present study has determined the complete mt genome sequences of E. hortense and assessed the phylogenetic relationships with other digenean species for which the complete mt genome sequences are available in GenBank using concatenated amino acid sequences inferred from 12 protein-coding genes. The mt genome of E. hortense contained 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 non-coding region. The length of the mt genome of E. hortense was 14,994 bp, which was somewhat smaller than those of other trematode species. Phylogenetic analyses based on concatenated nucleotide sequence datasets for all 12 protein-coding genes using maximum parsimony (MP) method showed that E. hortense and Hypoderaeum conoideum gathered together, and they were closer to each other than to Fasciolidae and other echinostomatid trematodes. The availability of the complete mt genome sequences of E. hortense provides important genetic markers for diagnostics, population genetics, and evolutionary studies of digeneans.     

The mitochondrial genomes of Amphiascoides atopus and Schizopera knabeni (Harpacticoida: Miraciidae) reveal similarities between the copepod orders Harpacticoida and Poecilostomatoida

Members of subclass Copepoda are abundant, diverse, and—as a result of their variety of ecological roles in marine and freshwater environments—important, but their phylogenetic interrelationships are unclear. Recent studies of arthropods have used gene arrangements in the mitochondrial (mt) genome to infer phylogenies, but for copepods, only seven complete mt genomes have been published. These data revealed several within-order and few among-order similarities. To increase the data available for comparisons, we sequenced the complete mt genome (13,831 base pairs) of Amphiascoides atopus and 10,649 base pairs of the mt genome of Schizopera knabeni (both in the family Miraciidae of the order Harpacticoida). Comparison of our data to those for Tigriopus japonicus (family Harpacticidae, order Harpacticoida) revealed similarities in gene arrangement among these three species that were consistent with those found within and among families of other copepod orders. Comparison of the mt genomes of our species with those known from other copepod orders revealed the arrangement of mt genes of our Harpacticoida species to be more similar to that of Sinergasilus polycolpus (order Poecilostomatoida) than to that of T. japonicus. The similarities between S. polycolpus and our species are the first to be noted across the boundaries of copepod orders and support the possibility that mt-gene arrangement might be used to infer copepod phylogenies. We also found that our two species had extremely truncated transfer RNAs and that gene overlaps occurred much more frequently than has been reported for other copepod mt genomes.     

The complete mitochondrial genomes of six heterodont bivalves (Tellinoidea and Solenoidea): variable gene arrangements and phylogenetic implications

Background

Taxonomy and phylogeny of subclass Heterodonta including Tellinoidea are long-debated issues and a complete agreement has not been reached yet. Mitochondrial (mt) genomes have been proved to be a powerful tool in resolving phylogenetic relationship. However, to date, only ten complete mitochondrial genomes of Heterodonta, which is by far the most diverse major group of Bivalvia, have been determined. In this paper, we newly sequenced the complete mt genomes of six species belonging to Heterodonta in order to resolve some problematical relationships among this subclass.

Principal Findings

The complete mt genomes of six species vary in size from 16,352 bp to 18,182. Hairpin-like secondary structures are found in the largest non-coding regions of six freshly sequenced mt genomes, five of which contain tandem repeats. It is noteworthy that two species belonging to the same genus show different gene arrangements with three translocations. The phylogenetic analysis of Heterodonta indicates that Sinonovacula constricta, distant from the Solecurtidae belonging to Tellinoidea, is as a sister group with Solen grandis of family Solenidae. Besides, all five species of Tellinoidea cluster together, while Sanguinolaria diphos has closer relationship with Solecurtus divaricatus, Moerella iridescens and Semele scaba rather than with Sanguinolaria olivacea.

Conclusions/Significance

By comparative study of gene order rearrangements and phylogenetic relationships of the five species belonging to Tellinoidea, our results support that comparisons of mt gene order rearrangements, to some extent, are a useful tool for phylogenetic studies. Based on phylogenetic analyses of multiple protein-coding genes, we prefer classifying the genus Sinonovacula within the superfamily Solenoidea and not the superfamily Tellinoidea. Besides, both gene order and sequence data agree that Sanguinolaria (Psammobiidae) is not monophyletic. Nevertheless, more studies based on more mt genomes via combination of gene order and phylogenetic analysis are needed to further understand the phylogenetic relationships in subclass Heterodonta.     

The Complete Maternally and Paternally Inherited Mitochondrial Genomes of the Endangered Freshwater Mussel Solenaia carinatus (Bivalvia: Unionidae) and Implications for Unionidae Taxonomy

Doubly uniparental inheritance (DUI) is an exception to the typical maternal inheritance of mitochondrial (mt) DNA in Metazoa, and found only in some bivalves. In species with DUI, there are two highly divergent gender-associated mt genomes: maternal (F) and paternal (M), which transmit independently and show different tissue localization. Solenaia carinatus is an endangered freshwater mussel species exclusive to Poyang Lake basin, China. Anthropogenic events in the watershed greatly threaten the survival of this species. Nevertheless, the taxonomy of S. carinatus based on shell morphology is confusing, and the subfamilial placement of the genus Solenaia remains unclear. In order to clarify the taxonomic status and discuss the phylogenetic implications of family Unionidae, the entire F and M mt genomes of S. carinatus were sequenced and compared with the mt genomes of diverse freshwater mussel species. The complete F and M mt genomes of S. carinatus are 16716 bp and 17102 bp in size, respectively. The F and M mt genomes of S. carinatus diverge by about 40% in nucleotide sequence and 48% in amino acid sequence. Compared to F counterparts, the M genome shows a more compact structure. Different gene arrangements are found in these two gender-associated mt genomes. Among these, the F genome cox2-rrnS gene order is considered to be a genome-level synapomorphy for female lineage of the subfamily Gonideinae. From maternal and paternal mtDNA perspectives, the phylogenetic analyses of Unionoida indicate that S. carinatus belongs to Gonideinae. The F and M clades in freshwater mussels are reciprocal monophyly. The phylogenetic trees advocate the classification of sampled Unionidae species into four subfamilies: Gonideinae, Ambleminae, Anodontinae, and Unioninae, which is supported by the morphological characteristics of glochidia.     

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.     

Rapid sequencing of the bamboo mitochondrial genome using Illumina technology and parallel episodic evolution of organelle genomes in grasses

Background

Compared to their counterparts in animals, the mitochondrial (mt) genomes of angiosperms exhibit a number of unique features. However, unravelling their evolution is hindered by the few completed genomes, of which are essentially Sanger sequenced. While next-generation sequencing technologies have revolutionized chloroplast genome sequencing, they are just beginning to be applied to angiosperm mt genomes. Chloroplast genomes of grasses (Poaceae) have undergone episodic evolution and the evolutionary rate was suggested to be correlated between chloroplast and mt genomes in Poaceae. It is interesting to investigate whether correlated rate change also occurred in grass mt genomes as expected under lineage effects. A time-calibrated phylogenetic tree is needed to examine rate change.

Methodology/Principal Findings

We determined a largely completed mt genome from a bamboo, Ferrocalamus rimosivaginus (Poaceae), through Illumina sequencing of total DNA. With combination of de novo and reference-guided assembly, 39.5-fold coverage Illumina reads were finally assembled into scaffolds totalling 432,839 bp. The assembled genome contains nearly the same genes as the completed mt genomes in Poaceae. For examining evolutionary rate in grass mt genomes, we reconstructed a phylogenetic tree including 22 taxa based on 31 mt genes. The topology of the well-resolved tree was almost identical to that inferred from chloroplast genome with only minor difference. The inconsistency possibly derived from long branch attraction in mtDNA tree. By calculating absolute substitution rates, we found significant rate change (∼4-fold) in mt genome before and after the diversification of Poaceae both in synonymous and nonsynonymous terms. Furthermore, the rate change was correlated with that of chloroplast genomes in grasses.

Conclusions/Significance

Our result demonstrates that it is a rapid and efficient approach to obtain angiosperm mt genome sequences using Illumina sequencing technology. The parallel episodic evolution of mt and chloroplast genomes in grasses is consistent with lineage effects.     

Analysis of organellar genomes in brown algae reveals an independent introduction of similar foreign sequences into the mitochondrial genome

Kelp species (Laminariales, Phaeophyceae) are globally widespread along temperate to Polar rocky coastal lines. Here we analyse the mitochondrial and chloroplast genomes of Laminaria rodriguezii, in comparison to the organellar genomes of other kelp species. We also provide the complete mitochondrial genome sequence of another endemic kelp species from a Polar habitat, the Arctic Laminaria solidungula. We compare phylogenetic trees derived from twenty complete mitochondrial and seven complete chloroplast kelp genomes. Interestingly, we found a stretch of more than 700 bp in the mitochondrial genome of L.rodriguezii, which is not present in any other yet sequenced member of the Phaeophyceae. This stretch matches a protein coding region in the mitochondrial genome from Desmarestia viridis, another brown seaweed. Their high similarity suggests that these sequences originated through independent introduction into the two species. Their origin could have been by infection by yet unknown similar mitoviruses, currently only known from fungi and plants.     

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

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

Mitochondrial genomes of the genus Ephydatia Lamouroux, 1816: can palindromic elements be used in species-level studies?

Poriferan mitochondrial DNA (mtDNA), especially large intergenic regions, is a target for the insertion of repetitive hairpin-forming elements. These elements are responsible for the large mt genome size differences observed even among closely related sponge taxa. In this study, we present the new, nearly complete, mt genome sequence of Ephydatia fluviatilis and compare it with previously published mt genomes of freshwater sponges. Special emphasis was placed on comparison with the closely related species Ephydatia muelleri, thereby comparing the only two species of the genus Ephydatia on the western Balkan Peninsula. In particular, we analyzed repetitive palindromic elements within the mitochondrial intergenic regions. The genomic distribution of these repetitive elements was analyzed and their potential role in the evolution of mt genomes discussed. We show here that palindromic elements are widespread through the whole mt genome, including the protein coding genes, thus introducing genetic variability into mt genomes.     

Mitochondrial genome evidence suggests Cooperia sp. from China may represent a distinct species from Cooperia oncophora from Australia

Cooperia spp. are parasitic nematodes parasitizing in small intestine of ruminants with a worldwide distribution. Infection of ruminants with Cooperia species can cause severe enteritis, causing significant socio-economic losses to the livestock industry. However, it is yet to know whether there is genetic diversity in mitochondrial (mt) DNA sequences of Cooperia nematodes from different geographic regions. The objective of the present study was to examine sequence difference in mt genomes between Cooperia sp. from China and other Cooperia species. We determined the sequences of the internal transcribed spacer (ITS-1 and ITS-2) of nuclear ribosomal DNA (rDNA) of 11 Cooperia specimens collected from the small intestine of a Tianzhu White yak in Gansu Province, northwestern China, which had 99% similarity with that of C. oncophora from Brazil (GenBank accession Number: AJ544290) in ITS-1, and 99% similarity with those from Denmark (AB245040), Scotland and Australia (AJ000032) in ITS-2, indicating that specimens used in the present study should at least represent parasites in Cooperia. We then determined the complete mt genome sequences of one representative specimen of Cooperia sp. from China (CspC), compared the mt DNA sequences with that of C. oncophora from Australia (COA, GQ888713), and conducted phylogenetic analysis with selected nematodes using both maximum likelihood (ML) and Bayesian inference (BI) methods based on both concatenated 12 PCGs, rrnL and rrnS sequences and partial cox2 sequences. The complete mt genome sequence of CspC (KY769271) is 13, 583 bp in length, which is 91 bp shorter than that from COA. The sequence difference over the entire mt genome between CspC and COA was 12.2% in nucleotide and 6.3% in inferred amino acids, with nad4L and nad1 being the most variable and the most conserved PCGs, respectively. Phylogenetic analysis indicated that CspC and COA were closely-related but distinct taxa. The determination of mt genome sequences for Cooperia sp. from China also provides novel resources for further studies of taxonomy, systematics and population genetics of Cooperia from different geographical locations.     

Complete mitochondrial genome sequences of three bats species and whole genome mitochondrial analyses reveal patterns of codon bias and lend support to a basal split in Chiroptera

Order Chiroptera is a unique group of mammals whose members have attained self-powered flight as their main mode of locomotion. Much speculation persists regarding bat evolution; however, lack of sufficient molecular data hampers evolutionary and conservation studies. Of ~ 1200 species, complete mitochondrial genome sequences are available for only eleven. Additional sequences should be generated if we are to resolve many questions concerning these fascinating mammals. Herein, we describe the complete mitochondrial genomes of three bats: Corynorhinus rafinesquii, Lasiurus borealis and Artibeus lituratus. We also compare the currently available mitochondrial genomes and analyze codon usage in Chiroptera. C. rafinesquii, L. borealis and A. lituratus mitochondrial genomes are 16438 bp, 17048 bp and 16709 bp, respectively. Genome organization and gene arrangements are similar to other bats. Phylogenetic analyses using complete mitochondrial genome sequences support previously established phylogenetic relationships and suggest utility in future studies focusing on the evolutionary aspects of these species. Comprehensive analyses of available bat mitochondrial genomes reveal distinct nucleotide patterns and synonymous codon preferences corresponding to different chiropteran families. These patterns suggest that mutational and selection forces are acting to different extents within Chiroptera and shape their mitochondrial genomes.     

The Mitochondrial Genome of an Aquatic Plant,Spirodela polyrhiza

Background

Spirodela polyrhiza is a species of the order Alismatales, which represent the basal lineage of monocots with more ancestral features than the Poales. Its complete sequence of the mitochondrial (mt) genome could provide clues for the understanding of the evolution of mt genomes in plant.

Methods

Spirodela polyrhiza mt genome was sequenced from total genomic DNA without physical separation of chloroplast and nuclear DNA using the SOLiD platform. Using a genome copy number sensitive assembly algorithm, the mt genome was successfully assembled. Gap closure and accuracy was determined with PCR products sequenced with the dideoxy method.

Conclusions

This is the most compact monocot mitochondrial genome with 228,493 bp. A total of 57 genes encode 35 known proteins, 3 ribosomal RNAs, and 19 tRNAs that recognize 15 amino acids. There are about 600 RNA editing sites predicted and three lineage specific protein-coding-gene losses. The mitochondrial genes, pseudogenes, and other hypothetical genes (ORFs) cover 71,783 bp (31.0%) of the genome. Imported plastid DNA accounts for an additional 9,295 bp (4.1%) of the mitochondrial DNA. Absence of transposable element sequences suggests that very few nuclear sequences have migrated into Spirodela mtDNA. Phylogenetic analysis of conserved protein-coding genes suggests that Spirodela shares the common ancestor with other monocots, but there is no obvious synteny between Spirodela and rice mtDNAs. After eliminating genes, introns, ORFs, and plastid-derived DNA, nearly four-fifths of the Spirodela mitochondrial genome is of unknown origin and function. Although it contains a similar chloroplast DNA content and range of RNA editing as other monocots, it is void of nuclear insertions, active gene loss, and comprises large regions of sequences of unknown origin in non-coding regions. Moreover, the lack of synteny with known mitochondrial genomic sequences shed new light on the early evolution of monocot mitochondrial genomes.     

The Complete Mitochondrial Genomes of Six Species of Tetranychus Provide Insights into the Phylogeny and Evolution of Spider Mites

Many spider mites belonging to the genus Tetranychus are of agronomical importance. With limited morphological characters, Tetranychus mites are usually identified by a combination of morphological characteristics and molecular diagnostics. To clarify their molecular evolution and phylogeny, the mitochondrial genomes of the green and red forms of Tetranychus urticae as well as T. kanzawai, T. ludeni, T. malaysiensis, T. phaselus, T. pueraricola were sequenced and compared. The seven mitochondrial genomes are typical circular molecules of about 13,000 bp encoding and they are composed of the complete set of 37 genes that are usually found in metazoans. The order of the mitochondrial (mt) genes is the same as that in the mt genomes of Panonychus citri and P. ulmi, but very different from that in other Acari. The J-strands of the mitochondrial genomes have high (∼84%) A+T contents, negative GC-skews and positive AT-skews. The nucleotide sequence of the cox1 gene, which is commonly used as a taxon barcode and molecular marker, is more highly conserved than the nucleotide sequences of other mitochondrial genes in these seven species. Most tRNA genes in the seven genomes lose the D-arm and/or the T-arm. The functions of these tRNAs need to be evaluated. The mitochondrial genome of T. malaysiensis differs from the other six genomes in having a slightly smaller genome size, a slight difference in codon usage, and a variable loop in place of the T-arm of some tRNAs by a variable loop. A phylogenic analysis shows that T. malaysiensis first split from other Tetranychus species and that the clade of the family Tetranychoidea occupies a basal position in the Trombidiformes. The mt genomes of the green and red forms of T. urticae have limited divergence and short evolutionary distance.     

The First Mitochondrial Genomes of Antlion (Neuroptera: Myrmeleontidae) and Split-footed Lacewing (Neuroptera: Nymphidae), with Phylogenetic Implications of Myrmeleontiformia

In the holometabolous insect order Neuroptera (lacewings), the cosmopolitan Myrmeleontidae (antlions) are the most species-rich family, while the closely related Nymphidae (split-footed lacewings) are a small endemic family from the Australian-Malesian region. Both families belong to the suborder Myrmeleontiformia, within which controversial hypotheses on the interfamilial phylogenetic relationships exist. Herein, we describe the complete mitochondrial (mt) genomes of an antlion (Myrmeleon immanis Walker, 1853) and a split-footed lacewing (Nymphes myrmeleonoides Leach, 1814), representing the first mt genomes for both families. These mt genomes are relatively small (respectively composed of 15,799 and 15,713 bp) compared to other lacewing mt genomes, and comprise 37 genes (13 protein coding genes, 22 tRNA genes and two rRNA genes). The arrangement of these two mt genomes is the same as in most derived Neuroptera mt genomes previously sequenced, specifically with a translocation of trnC. The start codons of all PCGs are started by ATN, with an exception of cox1, which is ACG in the M. immanis mt genome and TCG in N. myrmeleonoides. All tRNA genes have a typical clover-leaf structure of mitochondrial tRNA, with the exception of trnS1(AGN). The secondary structures of rrnL and rrnS are similar with those proposed insects and the domain I contains nine helices rather than eight helices, which is common within Neuroptera. A phylogenetic analysis based on the mt genomic data for all Neuropterida sequenced thus far, supports the monophyly of Myrmeleontiformia and the sister relationship between Ascalaphidae and Myrmeleontidae.     

The First Mitochondrial Genome of the Sepsid Fly Nemopoda mamaevi Ozerov, 1997 (Diptera: Sciomyzoidea: Sepsidae), with Mitochondrial Genome Phylogeny of Cyclorrhapha

Sepsid flies (Diptera: Sepsidae) are important model insects for sexual selection research. In order to develop mitochondrial (mt) genome data for this significant group, we sequenced the first complete mt genome of the sepsid fly Nemopoda mamaevi Ozerov, 1997. The circular 15,878 bp mt genome is typical of Diptera, containing all 37 genes usually present in bilaterian animals. We discovered inaccurate annotations of fly mt genomes previously deposited on GenBank and thus re-annotated all published mt genomes of Cyclorrhapha. These re-annotations were based on comparative analysis of homologous genes, and provide a statistical analysis of start and stop codon positions. We further detected two 18 bp of conserved intergenic sequences from tRNA^Glu-tRNA^Phe and ND1-tRNA^Ser(UCN) across Cyclorrhapha, which are the mtTERM binding site motifs. Additionally, we compared automated annotation software MITOS with hand annotation method. Phylogenetic trees based on the mt genome data from Cyclorrhapha were inferred by Maximum-likelihood and Bayesian methods, strongly supported a close relationship between Sepsidae and the Tephritoidea.     

The Mitochondrial Genome of Arctica islandica; Phylogeny and Variation

Arctica islandica is known as the longest-lived non-colonial metazoan species on earth and is therefore increasingly being investigated as a new model in aging research. As the mitochondrial genome is associated with the process of aging in many species and bivalves are known to possess a peculiar mechanism of mitochondrial genome inheritance including doubly uniparental inheritance (DUI), we aimed to assess the genomic variability of the A. islandica mitochondrial DNA (mtDNA). We sequenced the complete mitochondrial genomes of A. islandica specimens from three different sites in the Western Palaearctic (Iceland, North Sea, Baltic Sea). We found the A. islandica mtDNA to fall within the normal size range (18 kb) and exhibit similar coding capacity as other animal mtDNAs. The concatenated protein sequences of all currently known Veneroidea mtDNAs were used to robustly place A. islandica in a phylogenetic framework. Analysis of the observed single nucleotide polymorphism (SNP) patterns on further specimen revealed two prevailing haplotypes. Populations in the Baltic and the North Sea are very homogenous, whereas the Icelandic population, from which exceptionally old individuals have been collected, is the most diverse one. Homogeneity in Baltic and North Sea populations point to either stronger environmental constraints or more recent colonization of the habitat. Our analysis lays the foundation for further studies on A. islandica population structures, age research with this organism, and for phylogenetic studies. Accessions for the mitochondrial genome sequences: KC197241 Iceland; KF363951 Baltic Sea; KF363952 North Sea; KF465708 to KF465758 individual amplified regions from different speciemen     

Beyond barcoding: A mitochondrial genomics approach to molecular phylogenetics and diagnostics of blowflies (Diptera: Calliphoridae)

Members of the Calliphoridae (blowflies) are significant for medical and veterinary management, due to the ability of some species to consume living flesh as larvae, and for forensic investigations due to the ability of others to develop in corpses. Due to the difficulty of accurately identifying larval blowflies to species there is a need for DNA-based diagnostics for this family, however the widely used DNA-barcoding marker, cox1, has been shown to fail for several groups within this family. Additionally, many phylogenetic relationships within the Calliphoridae are still unresolved, particularly deeper level relationships. Sequencing whole mt genomes has been demonstrated both as an effective method for identifying the most informative diagnostic markers and for resolving phylogenetic relationships. Twenty-seven complete, or nearly so, mt genomes were sequenced representing 13 species, seven genera and four calliphorid subfamilies and a member of the related family Tachinidae. PCR and sequencing primers developed for sequencing one calliphorid species could be reused to sequence related species within the same superfamily with success rates ranging from 61% to 100%, demonstrating the speed and efficiency with which an mt genome dataset can be assembled. Comparison of molecular divergences for each of the 13 protein-coding genes and 2 ribosomal RNA genes, at a range of taxonomic scales identified novel targets for developing as diagnostic markers which were 117–200% more variable than the markers which have been used previously in calliphorids. Phylogenetic analysis of whole mt genome sequences resulted in much stronger support for family and subfamily-level relationships. The Calliphoridae are polyphyletic, with the Polleninae more closely related to the Tachinidae, and the Sarcophagidae are the sister group of the remaining calliphorids. Within the Calliphoridae, there was strong support for the monophyly of the Chrysomyinae and Luciliinae and for the sister-grouping of Luciliinae with Calliphorinae. Relationships within Chrysomya were not well resolved. Whole mt genome data, supported the previously demonstrated paraphyly of Lucilia cuprina with respect to L. sericata and allowed us to conclude that it is due to hybrid introgression prior to the last common ancestor of modern sericata populations, rather than due to recent hybridisation, nuclear pseudogenes or incomplete lineage sorting.