Involvement of two differenturf-s related mitochondrial sequences in the molecular evolution of the CMS-specificS-Pcf locus in petunia

In petunia, a mitochondrial (mt) locus,S-Pcf, has been found to be strongly associated with cytoplasmic male sterility (CMS). TheS-Pcf locus consists of three open reading frames (ORF) that are co-transcribed. The first ORF,Pcf, contains parts of theatp9 andcoxII genes and an unidentified reading frame,urf-s. The second and third ORFs contain NADH dehydrogenase subunit 3 (nad3) and ribosomal protein S12 (rps12) sequences, respectively. Thenad3 andrps12 sequences included in theS-Pcf locus are identical to the corresponding sequences on the mt genome of fertile petunia. In both CMS and fertile petunia, only a single copy ofnad3 andrps12 has been detected on the physical map of the main mt genome. The origin of theurf-s sequence and the molecular events leading to the formation of the chimericS-Pcf locus are not known. This paper presents evidence indicating that two different mt sequences, related tourf-s and found in fertile petunia lines (orf-h and Rf-1), might have been involved in the molecular evolution of theS-Pcf locus. Southern analysis of mtDNA derived from both fertile and sterile petunia plants suggests that one of theseurf-s related sequences (showing 100% homology tourf-s and termedorf-h) is located on a sublimon. An additional, low-homologyurf-s related sequence (Rf-1) is shown to be located on the main mt genome 5′ to thenad3 gene. It is, thus, suggested that the sequence of events leading to the generation of theS-Pcf locus might have involved introduction of theorf-h sequence, via homologous recombination, into the main mt genome 5′ tonad3 at the region where the Rf-1 sequence is located.     

Complete mitogenome of Nycteribia allotopa Speiser, 1901 (Diptera,Hippoboscoidea, Nycteribiidae) and comparative analysis of mitochondrial genomes of Nycteribiidae

In recent years, the global pandemic of bat-associated pathogens has led to increasing attention on bat ectoparasites. Numerous studies have identified human-associated pathogens in Nycteribiidae, indicating their potential as vectors. In this study, the first complete sequencing of the mitochondrial genome of Nycteribia allotopa Speiser, 1901 was sequenced and analyzed. We also compared the mitochondrial sequences of N. allotopa with those available in the database for other Nycteribiidae species. The complete mitochondrial genome of N. allotopa was found to be 15,161 bp in size with an A + T content of 82.49%. Nucleotide polymorphism analysis of 13 protein-coding genes from five species of Nycteribiidae showed that nad6 exhibited the most significant variation, while cox1 was the most conserved. Furthermore, selection pressure analysis revealed cox1 to exhibit the strongest purifying selection, while atp8, nad2, nad4L, and nad5 showed slightly looser purifying selection. Pairwise genetic distances indicated that cox1 and cox2 were evolving comparatively slowly, whereas atp8, nad2, and nad6 were evolving comparatively quickly. Phylogenetic trees constructed using Bayesian inference and maximum likelihood methods demonstrated that all four families within the superfamily Hippoboscoidea clustered into one branch each, indicating their monophyly. N. allotopa was found to be most closely related to the same genus N. parvula. This study significantly enriches the molecular database for Nycteribiidae and provides invaluable reference data for future species identification, phylogenetic analysis, and exploration of their potential as vectors for human-associated pathogens.     

The Australian fresh water isopod (Phreatoicidea: Isopoda) allows insights into the early mitogenomic evolution of isopods

The complete mitochondrial (mt) genome sequence of the Australian fresh water isopod Eophreatoicus sp.-14 has been determined. The new species is a member of the taxon Phreatoicidea, a clade of particular interest, as it is often regarded as the sister group to all other Isopoda. Although the overall genome organization of Eophreatoicus sp.-14 conforms to the typical state of Metazoa—it is a circular ring of DNA hosting the usual 37 genes and one major non-coding region—it bears a number of derived characters that fall within the scope of “genome morphology”. Earlier studies have indicated that the isopod mitochondrial gene order is not as conserved as that of other crustaceans. Indeed, the mt genome of Eophreatoicus sp.-14 shows an inversion of seven genes (including cox1), which is as far as we know unique. Even more interesting is the derived arrangement of nad1, trnL(CUN), rrnS, control region, cob, trnT, nad5 and trnF that is shared by nearly all available isopod mt genomes. A striking feature is the close proximity of the rearranged genes to the mt control region. Inferable gene translocation events are, however, more suitable to trace the evolution of mt genomes. Genes like nad1/trnL(CUN) and nad5/trnF, which retained their adjacent position after being rearranged, were most likely translocated together. A very good example for the need to understand the mechanisms of translocations is the remolding of trnL(UUR) to trnL(CUN). Both tRNA genes are adjacent and have a high sequence similarity, probably the result of a gene duplication and subsequent anticodon mutation. Modified secondary structures were found in three tRNAs of Eophreatoicus sp.-14, which are all characterized by the loss of the DHU-arm. This is common to crustaceans for tRNA Serine(AGY), while the arm-loss in tRNA Cysteine within Malacostraca is only shared by other isopods. Modification of the third tRNA, Isoleucine, is not known from any other related species. Nucleotide frequencies of genes have been found to be indirectly correlated to the orientation of the mitochondrial replication process. In Eophreatoicus sp.-14 and in other Isopoda the associated nucleotide bias is inversed to the state of other Malacostraca. This is a strong indication for an inversion of the control region that most likely evolved in the isopod ancestor.     

Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera)

The complete arrangement of genes in the mitochondrial (mt) genome is known for 12 species of insects, and part of the gene arrangement in the mt genome is known for over 300 other species of insects. The arrangement of genes in the mt genome is very conserved in insects studied, since all of the protein-coding and rRNA genes and most of the tRNA genes are arranged in the same way. We sequenced the entire mt genome of the wallaby louse, Heterodoxus macropus, which is 14,670 bp long and has the 37 genes typical of animals and some noncoding regions. The largest noncoding region is 73 bp long (93% A+T), and the second largest is 47 bp long (92% A+T). Both of these noncoding regions seem to be able to form stem-loop structures. The arrangement of genes in the mt genome of this louse is unlike that of any other animal studied. All tRNA genes have moved and/or inverted relative to the ancestral gene arrangement of insects, which is present in the fruit fly Drosophila yakuba. At least nine protein-coding genes (atp6, atp8, cox2, cob, nad1-nad3, nad5, and nad6) have moved; moreover, four of these genes (atp6, atp8, nad1, and nad3) have inverted. The large number of gene rearrangements in the mt genome of H. macropus is unprecedented for an arthropod.     

Making the most of mitochondrial genomes--markers for phylogeny, molecular ecology and barcodes in Schistosoma (Platyhelminthes: Digenea)

An increasing number of complete sequences of mitochondrial (mt) genomes provides the opportunity to optimise the choice of molecular markers for phylogenetic and ecological studies. This is particularly the case where mt genomes from closely related taxa have been sequenced; e.g., within Schistosoma. These blood flukes include species that are the causative agents of schistosomiasis, where there has been a need to optimise markers for species and strain recognition. For many phylogenetic and population genetic studies, the choice of nucleotide sequences depends primarily on suitable PCR primers. Complete mt genomes allow individual gene or other mt markers to be assessed relative to one another for potential information content, prior to broad-scale sampling. We assess the phylogenetic utility of individual genes and identify regions that contain the greatest interspecific variation for molecular ecological and diagnostic markers. We show that variable characters are not randomly distributed along the genome and there is a positive correlation between polymorphism and divergence. The mt genomes of African and Asian schistosomes were compared with the available intraspecific dataset of Schistosoma mansoni through sliding window analyses, in order to assess whether the observed polymorphism was at a level predicted from interspecific comparisons. We found a positive correlation except for the two genes (cox1 and nad1) adjoining the putative control region in S. mansoni. The genes nad1, nad4, nad5, cox1 and cox3 resolved phylogenies that were consistent with a benchmark phylogeny and in general, longer genes performed better in phylogenetic reconstruction. Considering the information content of entire mt genome sequences, partial cox1 would not be the ideal marker for either species identification (barcoding) or population studies with Schistosoma species. Instead, we suggest the use of cox3 and nad5 for both phylogenetic and population studies. Five primer pairs designed against Schistosoma mekongi and Schistosoma malayensis were tested successfully against Schistosoma japonicum. In combination, these fragments encompass 20-27% of the variation amongst the genomes (average total length approximately 14,000bp), thus providing an efficient means of encapsulating the greatest amount of variation within the shortest sequence. Comparative mitogenomics provides the basis of a rational approach to molecular marker selection and optimisation.     

Evolution of multipartite mitochondrial genomes in the booklice of the genus Liposcelis (Psocoptera)

Background

The genus Liposcelis (Psocoptera: Troctomorpha) has more than 120 species with a worldwide distribution and they pose a risk for global food security. The organization of mitochondrial (mt) genomes varies between the two species of booklice investigated in the genus Liposcelis. Liposcelis decolor has its mt genes on a single chromosome, like most other insects; L. bostrychophila, however, has a multipartite mt genome with genes on two chromosomes.

Results

To understand how multipartite mt genome organization evolved in the genus Liposcelis, we sequenced the mt genomes of L. entomophila and L. paeta in this study. We found that these two species of booklice also have multipartite mt genomes, like L. bostrychophila, with the mt genes we identified on two chromosomes. Numerous pseudo mt genes and non-coding regions were found in the mt genomes of these two booklice, and account for 30% and 10% respectively of the entire length we sequenced. In L. bostrychophila, the mt genes are distributed approximately equally between the two chromosomes. In L. entomophila and L. paeta, however, one mt chromosome has most of the genes we identified whereas the other chromosome has largely pseudogenes and non-coding regions. L. entomophila and L. paeta differ substantially from each other and from L. bostrychophila in gene content and gene arrangement in their mt chromosomes.

Conclusions

Our results indicate unusually fast evolution in mt genome organization in the booklice of the genus Liposcelis, and reveal different patterns of mt genome fragmentation among L. bostrychophila, L. entomophila and L. paeta.

Electronic supplementary material

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

Complete sequence of the mitochondrial genome of Tetrahymena thermophila and comparative methods for identifying highly divergent genes

The complete sequence of the mitochondrial genome of Tetrahymena thermophila has been determined and compared with the mitochondrial genome of Tetrahymena pyriformis. The sequence similarity clearly indicates homology of the entire T.thermophila and T.pyriformis mitochondrial genomes. The T.thermophila genome is very compact, most of the intergenic regions are short (only three are longer than 63 bp) and comprise only 3.8% of the genome. The nad9 gene is tandemly duplicated in T.thermophila. Long terminal inverted repeats and the nad9 genes are undergoing concerted evolution. There are 55 putative genes: three ribosomal RNA genes, eight transfer RNA genes, 22 proteins with putatively assigned functions and 22 additional open reading frames of unknown function. In order to extend indications of homology beyond amino acid sequence similarity we have examined a number of physico-chemical properties of the mitochondrial proteins, including theoretical pI, molecular weight and particularly the predicted transmembrane spanning regions. This approach has allowed us to identify homologs to ymf58 (nad4L), ymf62 (nad6) and ymf60 (rpl6).     

Complete mitochondrial genome of Bactrocera ritsemai (Insecta: Tephritidae) and phylogenetic relationship with its congeners and related tephritid taxa

Bactrocera ritsemai is a dacine fruit fly found in Indonesia. We report here the complete mitogenome of this fruit fly from Lombok, Indonesia determined by Illumina MiSeq sequencing and its phylogenetic relationship with its congeners and related tephritid taxa. The whole mitogenome of B. ritsemai had a total length of 15,927 bp, comprising 37 genes – 13 protein-coding genes (PCGs), 2 ribosomal ribonucleic acid (rRNA) and 22 transfer ribonucleic acid (tRNA) genes – and a control region (D-loop). Of the PCGs, 6 (atp6, cob, cox2, cox3, nad4, nad4l) had ATG start codon, 4 (nad2, nad3, nad5, nad6) had ATT, and one each had ATA (nad1), GTG (atp8) and TCG (cox1). Seven PCGs (atp6, atp8, cox2, cox3, nad2, nad4l, nad6) had TAA stop codon, 3 (cob, nad3, nad4) had TAG, and 3 had incomplete stop codon (cox1 – TA; nad1, nad5 – T). The TΨC-loop of tRNA was absent in trnF while trnS1 lacked the DHU-loop. Phylogenetic analysis based on 15 mt-genes (13 PCGs + 2 rRNA genes) indicated B. ritsemai forming a sister group with B. umbrosa and the subgenus Bactrocera was monophyletic. The genera Bactrocera and Zeugodacus were monophyletic while the subfamilies Dacinae and Tephritinae were paraphyletic. A broader taxa sampling of the Tephritidae is needed to better elucidate the phylogenetics and systematics of the tribes and subfamilies of tephritid fruit flies.     

The Mitochondrial Genome of a Deep-Sea Bamboo Coral (Cnidaria,Anthozoa, Octocorallia,Isididae): Genome Structure and Putative Origins of Replication Are Not Conserved Among Octocorals

Octocoral mitochondrial (mt) DNA is subject to an exceptionally low rate of substitution, and it has been suggested that mt genome content and structure are conserved across the subclass, an observation that has been supported for most octocorallian families by phylogenetic analyses using PCR products spanning gene boundaries. However, failure to recover amplification products spanning the nad4L-msh1 gene junction in species from the family Isididae (bamboo corals) prompted us to sequence the complete mt genome of a deep-sea bamboo coral (undescribed species). Compared to the "typical" octocoral mt genome, which has 12 genes transcribed on one strand and 5 genes on the opposite (cox2, atp8, atp6, cox3, trnM), in the bamboo coral genome a contiguous string of 5 genes (msh1, rnl, nad2, nad5, nad4) has undergone an inversion, likely in a single event. Analyses of strand-specific compositional asymmetry suggest that (i) the light-strand origin of replication was also inverted and is adjacent to nad4, and (ii) the orientation of the heavy-strand origin of replication (OriH) has reversed relative to that of previously known octocoral mt genomes. Comparative analyses suggest that intramitochondrial recombination and errors in replication at OriH may be responsible for changes in gene order in octocorals and hexacorals, respectively. Using primers flanking the regions at either end of the inverted set of five genes, we examined closely related taxa and determined that the novel gene order is restricted to the deep-sea subfamily Keratoisidinae; however, we found no evidence for strand-specific mutational biases that may influence phylogenetic analyses that include this subfamily of bamboo corals.     

The complete mitochondrial genome of Evania appendigaster (Hymenoptera: Evaniidae) has low A+T content and a long intergenic spacer between atp8 and atp6

The apocritan Hymenoptera show extraordinary features in mitochondrial genomes, but no complete sequence has been reported for the basal lineage, Evanioidea. Here, we sequenced the complete mitochondrial genome of Evania appendigaster. This genome is 17,817 bp long; with low A+T content, 77.8%, compared with other hymenopteran species. Four tRNA genes were rearranged, among which remote inversion is the dominant gene rearrangement event. Gene shuffling is caused by tandem duplication-random loss while remote inversion is best explained by recombination. The start codon of nad1 was found as TTG, which might be common across Hymenoptera. trnS2 and trnK use abnormal anticodons TCT and TTT, respectively, and the D-stem pairings in trnS2 are absent. The secondary structure of two rRNA genes are predicted and compared with those in other insects. Five long intergenic spacers were present, including a long intergenic spacer between atp8 and atp6, where these two genes overlap in the previously reported animal genomes. A conserved motif was found between trnS1 and nad1, which is proposed to be associated with mtTERM. The A+T-rich region is 2,325 bp long, among the longest in insects, and contains a tandem repeat region.     

Complete mitochondrial DNA sequence of the ark shell Scapharca broughtonii: An ultra-large metazoan mitochondrial genome

The complete mitochondrial (mt) genome of the ark shell Scapharca broughtonii was determined using long PCR and a genome walking sequencing strategy with genus-specific primers. The S. broughtonii mt genome (GenBank accession number AB729113) contained 12 protein-coding genes (the atp8 gene is missing, as in most bivalves), 2 ribosomal RNA genes, and 42 transfer tRNA genes, in a length of 46,985 nucleotides for the size of mtDNA with only one copy of the heteroplasmic tandem repeat (HTR) unit. Moreover the S. broughtonii mt genome shows size variation; these genomes ranged in size from about 47 kb to about 50 kb because of variation in the number of repeat sequences in the non-coding region. The mt-genome of S. broughtonii is, to date, the longest reported metazoan mtDNA sequence. Sequence duplication in non-coding region and the formation of HTR arrays were two of the factors responsible for the ultra-large size of this mt genome. All the tRNA genes were found within the S. broughtonii mt genome, unlike the other bivalves usually lacking one or more tRNA genes. Twelve additional specimens were used to analyze the patterns of tandem repeat arrays by PCR amplification and agarose electrophoresis. Each of the 12 specimens displayed extensive heteroplasmy and had 8–10 length variants. The motifs of the HTR arrays are about 353–362 bp and the number of repeats ranges from 1 to 11.     

The mitochondrial genome of Priapulus caudatus Lamarck (Priapulida: Priapulidae)

We sequenced and annotated the complete mitochondrial (mt) genome of the priapulid Priapulus caudatus in order to provide a source of phylogenetic characters including an assessment of gene order arrangement. The genome was 14,919 bp in its entirety with few, short non-coding regions. A number of protein-coding and tRNA genes overlapped, making the genome relatively compact. The gene order was: cox1, cox2, trnK, trnD, atp8, atp6, cox3, trnG, nad3, trnA, trnR, trnN, rrnS, trnV, rrnL, trnL(yaa), trnL(nag), nad1, -trnS(nga), -cob, -nad6, trnP, -trnT, nad4L, nad4, trnH, nad5, trnF, -trnE, -trnS(nct), trnI, -trnQ, trnM, nad2, trnW, -trnC, -trnY; where '-' indicates genes transcribed on the opposite strand. The gene order, although unique amongst Metazoa, shared the greatest number of gene boundaries and the longest contiguous fragments with the chelicerate Limulus polyphemus. The mt genomes of these taxa differed only by a single inversion of 18 contiguous genes bounded by rrnS and trnS(nct). Other arthropods and nematodes shared fewer gene boundaries but considerably more than the most similar non-ecdysozoan.     

Comparative analyses of the complete mitochondrial genomes of the two ruminant hookworms Bunostomum trigonocephalum and Bunostomum phlebotomum

Bunostomum trigonocephalum and Bunostomum phlebotomum are blood-feeding hookworms of sheep and cattle, causing considerable economic losses to the live stock industries. Studying genetic variability within and among hookworm populations is critical to addressing epidemiological and ecological questions. Mitochondrial (mt) DNA is known to provide useful markers for investigations of population genetics of hookworms, but mt genome sequence data are scant. In the present study, the complete mitochondrial DNA (mtDNA) sequences of the sheep and goat hookworm B. trigonocephalum were determined for the first time, and the mt genome of B. phlebotomum from yak in China was also sequenced for comparative analyses of their gene contents and genome organizations. The lengths of mt DNA sequences of B. trigonocephalum sheep isolate, B.trigonocephalum goat isolate and B. phlebotomum China yak isolate were 13,764 bp, 13,771 bp and 13,803 bp in size, respectively. The identity of the mt genomes was 99.7% between B. trigonocephalum sheep isolate and B. trigonocephalum goat isolate. The identity of B. phlebotomum China yak isolate mt genomes was 85.3% with B. trigonocephalum sheep isolate, and 85.2% with B. trigonocephalum goat isolate. All the mt genes of the two hookworms were transcribed in the same direction and gene arrangements were consistent with those of the GA3 type, including 12 protein-coding genes, 2 rRNA genes and 22 tRNA genes, but lacking ATP synthetase subunit 8 gene. The mt genomes of B. trigonocephalum and B. phlebotomum were similar to prefer bases A and T, the contents of A + T are 76.5% (sheep isolate), 76.4% (goat isolate) and 76.9% (China yak isolate), respectively. Phylogenetic relationships reconstructed using concatenated amino acid sequences of 12 protein-coding genes with three methods (maximum likelihood, Bayesian inference and neighbor joining) revealed that the B. trigonocephalum and B. phlebotomum represent distinct but closely-related species. These data provide novel and useful genetic markers for studying the systematics, and population genetics of the two ruminant hookworms.     

The mitochondrial genome of Baylisascaris procyonis

Background

Baylisascaris procyonis (Nematoda: Ascaridida), an intestinal nematode of raccoons, is emerging as an important helminthic zoonosis due to serious or fatal larval migrans in animals and humans. Despite its significant veterinary and public health impact, the epidemiology, molecular ecology and population genetics of this parasite remain largely unexplored. Mitochondrial (mt) genomes can provide a foundation for investigations in these areas and assist in the diagnosis and control of B. procyonis. In this study, the first complete mt genome sequence of B. procyonis was determined using a polymerase chain reaction (PCR)-based primer-walking strategy.

Methodology/Principal Findings

The circular mt genome (14781 bp) of B. procyonis contained 12 protein-coding, 22 transfer RNA and 2 ribosomal RNA genes congruent with other chromadorean nematodes. Interestingly, the B. procyonis mtDNA featured an extremely long AT-rich region (1375 bp) and a high number of intergenic spacers (17), making it unique compared with other secernentean nematodes characterized to date. Additionally, the entire genome displayed notable levels of AT skew and GC skew. Based on pairwise comparisons and sliding window analysis of mt genes among the available 11 Ascaridida mtDNAs, new primer pairs were designed to amplify specific short fragments of the genes cytb (548 bp fragment) and rrnL (200 bp fragment) in the B. procyonis mtDNA, and tested as possible alternatives to existing mt molecular beacons for Ascaridida. Finally, phylogenetic analysis of mtDNAs provided novel estimates of the interrelationships of Baylisasaris and Ascaridida.

Conclusions/Significance

The complete mt genome sequence of B. procyonis sequenced here should contribute to molecular diagnostic methods, epidemiological investigations and ecological studies of B. procyonis and other related ascaridoids. The information will be important in refining the phylogenetic relationships within the order Ascaridida and enriching the resource of markers for systematic, population genetic and evolutionary biological studies of parasitic nematodes of socio-economic importance.     

Fragmented mitochondrial genomes are present in both major clades of the blood-sucking lice (suborder Anoplura): evidence from two Hoplopleura rodent lice (family Hoplopleuridae)

Background

The suborder Anoplura contains 540 species of blood-sucking lice that parasitize over 840 species of eutherian mammals. Fragmented mitochondrial (mt) genomes have been found in the lice of humans, pigs, horses and rats from four families: Pediculidae, Pthiridae, Haematopinidae and Polyplacidae. These lice, eight species in total, are from the same major clade of the Anoplura. The mt genomes of these lice consist of 9–20 minichromosomes; each minichromosome is 1.5–4 kb in size and has 1–8 genes. To understand mt genome fragmentation in the other major clade of the Anoplura, we sequenced the mt genomes of two species of rodent lice in the genus Hoplopleura (family Hoplopleuridae).

Results

We identified 28 mt genes on 10 minichromosomes in the mouse louse, Ho. akanezumi; each minichromosome is 1.7–2.7 kb long and has 1–6 genes. We identified 34 mt genes on 11 minichromosomes in the rat louse, Ho. kitti; each minichromosome is 1.8–2.8 kb long and has 1–5 genes. Ho. akanezumi also has a chimeric minichromosome with parts of two rRNA genes and a full-length tRNA gene for tyrosine. These two rodent lice share the same pattern for the distribution of all of the protein-coding and rRNA genes but differ in tRNA gene content and gene arrangement in four minichromosomes. Like the four genera of blood-sucking lice that have been investigated in previous studies, the Hoplopleura species have four minichromosomes that are only found in this genus.

Conclusions

Our results indicate that fragmented mt genomes were present in the most recent common ancestor of the two major clades of the blood-sucking lice, which lived ~75 million years ago. Intra-genus variation in the pattern of mt genome fragmentation is common in the blood-sucking lice (suborder Anoplura) and genus-specific minichromosomes are potential synapomorphies. Future studies should expand into more species, genera and families of blood-sucking lice to explore further the phylogenetic utility of the novel features associated with fragmented mt genomes.

Electronic supplementary material

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

Involvement of two differenturf-s related mitochondrial sequences in the molecular evolution of the CMS-specificS-Pcf locus in petunia

In petunia, a mitochondrial (mt) locus,S-Pcf, has been found to be strongly associated with cytoplasmic male sterility (CMS). TheS-Pcf locus consists of three open reading frames (ORF) that are co-transcribed. The first ORF,Pcf, contains parts of theatp9 andcoxII genes and an unidentified reading frame,urf-s. The second and third ORFs contain NADH dehydrogenase subunit 3 (nad3) and ribosomal protein S12 (rps12) sequences, respectively. Thenad3 andrps12 sequences included in theS-Pcf locus are identical to the corresponding sequences on the mt genome of fertile petunia. In both CMS and fertile petunia, only a single copy ofnad3 andrps12 has been detected on the physical map of the main mt genome. The origin of theurf-s sequence and the molecular events leading to the formation of the chimericS-Pcf locus are not known. This paper presents evidence indicating that two different mt sequences, related tourf-s and found in fertile petunia lines (orf-h and Rf-1), might have been involved in the molecular evolution of theS-Pcf locus. Southern analysis of mtDNA derived from both fertile and sterile petunia plants suggests that one of theseurf-s related sequences (showing 100% homology tourf-s and termedorf-h) is located on a sublimon. An additional, low-homologyurf-s related sequence (Rf-1) is shown to be located on the main mt genome 5′ to thenad3 gene. It is, thus, suggested that the sequence of events leading to the generation of theS-Pcf locus might have involved introduction of theorf-h sequence, via homologous recombination, into the main mt genome 5′ tonad3 at the region where the Rf-1 sequence is located. Contribution [No. 1581-E (1995 series)] from the Agricultural Research Organization, The Volcani Center, Bet Dagan, Israel 50 250     

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 genome of Mahanta tanyae compared with other zygaenoid moths (Lepidoptera: Zygaenoidea)

The complete mitochondrial genome (mitogenome) of Mahanta tanyae was sequenced and extensively compared with all seven additionally reported zygaenoid mitogenomes. The M. tanyae mitogenome is circular, double-stranded, and 15,323 bp long. Gene content, gene order, and orientation are all typical of Lepidoptera, despite the existence of gene rearrangements for some other zygaenoid mitogenomes. Comparative analyses further showed that the incomplete termination codon T is consistently recognized in the mitochondrial cox1, cox2 and nad4 genes of all zygaenoid species, as well as in the nad5 gene in two limacodid species. Among 13 protein-coding genes, nad6 exhibits the highest evolutionary rate. The structure for each tRNA is highly conserved, including loss of the dihydorouidine (DHU) arm in trnS1 (AGN), but remarkable nucleotide variation exists, primarily in the pseudouridine (TψC) loops. Interestingly, in four species of Zygaenidae, the anticodons for trnS1 (AGN) are consistently UCU, instead of the routinely used codon GCU, in all three species of Limacodidae. In the intergenic region between trnS2 and nad1, a short sequence before the motif “ATACTAA” is present in the M. tanyae mitogenome that is unique among reported zygaenoid mitogenomes. In the A + T-rich region between the motif “ATTTA” and the microsatellite (AT)_n element, some nucleotides were present for most zygaenoid mitogenomes, which is, to our knowledge, rare even in reported lepidopteran mitogenomes. Phylogenetic analyses based on the combined 37 mitochondrial genes confirmed the position of M. tanyae in Limacodidae of the Zygaenoidea.     

A novel mitochondrial genome architecture in thrips (Insecta: Thysanoptera): extreme size asymmetry among chromosomes and possible recent control region duplication

Background

Multipartite mitochondrial genomes are very rare in animals but have been found previously in two insect orders with highly rearranged genomes, the Phthiraptera (parasitic lice), and the Psocoptera (booklice/barklice).

Results

We provide the first report of a multipartite mitochondrial genome architecture in a third order with highly rearranged genomes: Thysanoptera (thrips). We sequenced the complete mitochondrial genomes of two divergent members of the Scirtothrips dorsalis cryptic species complex. The East Asia 1 species has the single circular chromosome common to animals while the South Asia 1 species has a genome consisting of two circular chromosomes. The fragmented South Asia 1 genome exhibits extreme chromosome size asymmetry with the majority of genes on the large, 14.28 kb, chromosome and only nad6 and trnC on the 0.92 kb mini-circle chromosome. This genome also features paralogous control regions with high similarity suggesting a very recent origin of the nad6 mini-circle chromosome in the South Asia 1 cryptic species.

Conclusions

Thysanoptera, along with the other minor paraenopteran insect orders should be considered models for rapid mitochondrial genome evolution, including fragmentation. Continued use of these models will facilitate a greater understanding of recombination and other mitochondrial genome evolutionary processes across eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1672-4) contains supplementary material, which is available to authorized users.