The complete mitochondrial genome of the grand jackknife clam, Solen grandis (Bivalvia: Solenidae): a novel gene order and unusual non-coding region

Molluscs in general, and bivalves in particular, exhibit an extraordinary degree of mitochondrial gene order variation when compared with other metazoans. The complete mitochondrial genome of Solen grandis (Bivalvia: Solenidae) was determined using long-PCR and genome walking techniques. The entire mitochondrial genome sequence of S. grandis is 16,784 bp in length, and contains 36 genes including 12 protein-coding genes (atp8 is absent), 2 ribosomal RNAs, and 22 tRNAs. All genes are encoded on the same strand. Compared with other species, it bears a novel gene order. Besides these, we find a peculiar non-coding region of 435 bp with a microsatellite-like (TA)₁₂ element, poly-structures and many hairpin structures. In contrast to the available heterodont mitochondrial genomes from GenBank, the complete mtDNA of S. grandis has the shortest cox3 gene, and the longest atp6, nad4, nad5 genes.     

Nuclear and mitochondrial subunits from the white shrimp Litopenaeus vannamei F0F1 ATP-synthase complex: cDNA sequence, molecular modeling, and mRNA quantification of atp9 and atp6

We studied for the first time the ATP-synthase complex from shrimp as a model to understand the basis of crustacean bioenergetics since they are exposed to endogenous processes as molting that demand high amount of energy. We analyzed the cDNA sequence of two subunits of the Fo sector from mitochondrial ATP-synthase in the white shrimp Litopenaeus vannamei. The nucleus encoded atp9 subunit presents a 773 bp sequence, containing a signal peptide sequence only observed in crustaceans, and the mitochondrial encoded atp6 subunit presents a sequence of 675 bp, and exhibits high identity with homologous sequences from invertebrate species. ATP9 and ATP6 protein structural models interaction suggest specific functional characteristics from both proteins in the mitochondrial enzyme. Differences in the steady-state mRNA levels of atp9 and atp6 from five different tissues correlate with tissue function. Moreover, significant changes in the mRNA levels of both subunits at different molt stages were detected. We discussed some insights about the enzyme structure and the regulation mechanisms from both ATP-synthase subunits related to the energy requirements of shrimp.     

Reductions in complexity of mitochondrial genomes in lichen‐forming fungi shed light on genome architecture of obligate symbioses

Symbioses among co‐evolving taxa are often marked by genome reductions such as a loss of protein‐coding genes in at least one of the partners as a means of reducing redundancy or intergenomic conflict. To explore this phenomenon in an iconic yet under‐studied group of obligate symbiotic organisms, mitochondrial genomes of 22 newly sequenced and annotated species of lichenized fungi were compared to 167 mitochondrial genomes of nonlichenized fungi. Our results demonstrate the first broad‐scale loss of atp9 from mitochondria of lichenized fungi. Despite key functions in mitochondrial energy production, we show that atp9 has been independently lost in three different lineages spanning 10 of the 22 studied species. A search for predicted, functional copies of atp9 among genomes of other symbionts involved in each lichen revealed the full‐length, presumably functional copies of atp9 in either the photosynthetic algal partner or in other symbiotic fungi in all 10 instances. Together, these data yield evidence of an obligate symbiotic relationship in which core genomic processes have been streamlined, likely due to co‐evolution.     

Complete mitochondrial genome of Tubulipora flabellaris (Bryozoa: Stenolaemata): the first representative from the class Stenolaemata with unique gene order

Mitochondrial genomes play a significant role in the reconstruction of phylogenetic relationships within metazoans. There are still many controversies concerning the phylogenetic position of the phylum Bryozoa. In this research, we have finished the complete mitochondrial genome of one bryozoan (Tubulipora flabellaris), which is the first representative from the class Stenolaemata. The complete mitochondrial genome of T. flabellaris is 13,763 bp in length and contains 36 genes, which lacks the atp8 gene in contrast to the typical metazoan mitochondrial genomes. Gene arrangement comparisons indicate that the mitochondrial genome of T. flabellaris has unique gene order when compared with other metazoans. The four known bryozoans complete mitochondrial genomes also have very different gene arrangements, indicates that bryozoan mitochondrial genomes have experienced drastic rearrangements. To investigate the phylogenetic relationship of Bryozoa, phylogenetic analyses based on amino acid sequences of 11 protein coding genes (excluding atp6 and atp8) from 26 metazoan complete mitochondrial genomes were made utilizing Maximum Likelihood (ML) and Bayesian methods, respectively. The results indicate the monopoly of Lophotrochozoa and a close relationship between Chaetognatha and Bryozoa. However, more evidences are needed to clarify the relationship between two groups. Lophophorate appeared to be polyphyletic according to our analyses. Meanwhile, neither analysis supports close relationship between Branchiopod and Phoronida. Four bryozoans form a clade and the relationship among them is T. flabellaris + (F. hispida + (B. neritina + W. subtorquata)), which is in coincidence with traditional classification system.     

A Comparative Analysis of the <Emphasis Type="Italic">atp8</Emphasis> Gene Between a Cytoplasmic Male Sterile Line and Its Maintainer and Further Development of a Molecular Marker Specific to Male Sterile Cytoplasm in Kenaf (<Emphasis Type="Italic">Hibiscus cannabinus</Emphasis> L.)

In this study, the atp8 gene was cloned from the cytoplasmic male sterile (CMS) line UG93A and its maintainer line UG93B in kenaf. Its DNA sequence analysis showed that atp8 containing 480-bp, encoding 159 amino acid residues, and a 9-bp insertion was found at the 3′flanking sequence in UG93A compared with UG93B. The cDNA sequence of atp8 analyzed by RT-PCR indicated that there were five loci edited, but six loci edited in UG93B. The editing frequencies were higher in sterile cytoplasm than in fertile cytoplasm. The relative expression of atp8 analyzed by real-time PCR showed that the expressed level of atp8 in UG93A was lower than that of its maitainer UG93B and its F₁ hybrid UG93A/992 (a restore line). Furthermore, based on the difference of the 9-bp differences at the 3′flanking sequence of atp8 between UG93A and UG93B, a molecular marker specific to male sterile cytoplasm was developed, which can be used for indentifying whether any germplasm of kenaf is male sterile cytoplasm or male fertile cytoplasm.     

Phylogeny of Rhigonematomorpha based on the complete mitochondrial genome of Rhigonema thysanophora (Nematoda: Chromadorea)

We determined the complete mitochondrial genome sequence of Rhigonema thysanophora, the first representative of Rhigonematomorpha, and used this sequence along with 57 other nematode species for phylogenetic analyses. The R. thysanophora mtDNA is 15 015 bp and identical to all other chromadorean nematode mtDNAs published to date in that it contains 36 genes (lacking atp8) encoded in the same direction. Phylogenetic analyses of nucleotide and amino acid sequence data for the 12 protein‐coding genes recovered Rhigonematomorpha as the sister group to the heterakoid species, Ascaridia columbae (Ascaridomorpha). The organization of R. thysanophora mtDNA resembles the most common pattern for the Rhabditomorpha+Ascaridomorpha+Diplogasteromorpha clade in gene order, but with some substantial gene rearrangements. This similarity in gene order is in agreement with the sequence‐based analyses that indicate a close relationship between Rhigonematomorpha and Rhabditomorpha+Ascaridomorpha+Diplogasteromorpha. These results are consistent with certain analyses of nuclear SSU rDNA for R. thysanophora and some earlier classification systems that asserted phylogenetic affinity between Rhigonematomorpha and Ascaridomorpha, but inconsistent with morphology‐based phylogenetic hypotheses that suggested a close (taxonomic) relationship between rhigonematomorphs and oxyuridomorphs (pinworms). These observations must be tempered by noting that few rhigonematomorph species have been sequenced and included in phylogenetic analyses, and preliminary studies based on SSU rDNA suggest the group is not monophyletic. Additional mitochondrial genome sequences of rhigonematids are needed to characterize their phylogenetic relationships within Chromadorea, and to increase understanding of mitochondrial genome evolution.     

Characterization of a mitochondrial ORF from the gender-associated mtDNAs of Mytilus spp. (Bivalvia: Mytilidae): identification of the "missing" ATPase 8 gene

Bivalve species are characterized by extraordinary variability in terms of mitochondrial (mt) genome size, gene arrangement and tRNA gene number. Many species are thought to lack the mitochondrial protein-coding gene atp8. Of these species, the Mytilidae appears to be the only known taxon with doubly uniparental inheritance of mtDNA that does not possess the atp8 gene. This raises the question as to whether mytilids have completely lost the ATP8 protein, whether the gene has been transferred to the nucleus or whether they possess a highly modified version of the gene/protein that has led to its lack of annotation. In the present study, we re-investigated all complete (or nearly complete) F and M mytilid mt genomes previously sequenced for the presence of conserved open reading frames (ORFs) that might code for ATP8 and/or have other functional importance in these bivalves. We also revised the annotations of all available complete mitochondrial genomes of bivalves and nematodes that are thought to lack atp8 in an attempt to detect it. Our results indicate that a novel mytilid ORF of significant length (i.e., the ORF is >85 amino acids in length), with complete start and stop codons, is a candidate for the atp8 gene: (1) it possesses a pattern of evolution expected for a protein-coding gene evolving under purifying selection (i.e., the 3rd>1st>2nd codon pattern of evolution), (2) it is actively transcribed in Mytilus species, (3) it has one predicted transmembrane helix (as do other metazoan ATP8 proteins), (4) it has conserved functional motifs and (5), comparisons of its amino acid sequence with ATP8 sequences of other molluscan or bivalve species reveal similar hydropathy profiles. Furthermore, our revised annotations also confirmed the mt presence of atp8 in almost all bivalve species and in one nematode species. Our results thus support recognizing the presence of ATPase 8 in most bivalves mt genomes (if not all) rather than the continued characterization of these genomes as lacking this gene.     

Mitochondrial DNA data unveil highly divergent populations within the genus <Emphasis Type="Italic">Hynobius</Emphasis> (Caudata: Hynobiidae) in South Korea

Korean salamanders of the genus Hynobius are currently classified into 3 species, H. leechii, H. quelpaertensis, and H. yangi. To investigate the phylogenetic relationship of these species, we analyzed the partial sequence of mitochondrial cytochrome b gene (907 bp) of 197 specimens from 43 regions in South Korea. Of these specimens, 93 were additionally examined with 12S rRNA (799 bp). Based on the partial sequence of the mitochondrial cytochrome b gene and 12S rRNA, 89 and 36 haplotypes were defined, respectively, consisting of six subclades (H. leechii, H. quelpaertensis, H. yangi, HC1, HC2, and HC3). Among these subclades, the three subclades (HC1, HC2, and HC3) were clearly separated from the 3 previously reported species in the genus Hynobius. Pairwise sequence divergence between the six subclades ranged from 6.3 to 11.2% in cytochrome b gene and 2.0 to 4.3% in 12S rRNA. These results indicate there may be more divergent populations than the three currently described. Moreover, the estimation of divergence time revealed that the Hynobius species in South Korea diverged during the Miocene epoch, approximately 9 — 5 MYA. In addition, we confirmed the distribution of the three known species (H. leechii, H. quelpaertensis, and H. yangi) and determined the distributions of new, distinct groups (or subclades; HC1, HC1, and HC3). To more accurately establish the taxonomic status and population structure, further genetic, morphological, and ecological studies will be needed.     

Mitogenomic analysis of the genus <Emphasis Type="Italic">Panthera</Emphasis>

The complete sequences of the mitochondrial DNA genomes of Panthera tigris, Panthera pardus, and Panthera uncia were determined using the polymerase chain reaction method. The lengths of the complete mitochondrial DNA sequences of the three species were 16990, 16964, and 16773 bp, respectively. Each of the three mitochondrial DNA genomes included 13 protein-coding genes, 22 tRNA, two rRNA, one O_LR, and one control region. The structures of the genomes were highly similar to those of Felis catus, Acinonyx jubatus, and Neofelis nebulosa. The phylogenies of the genus Panthera were inferred from two combined mitochondrial sequence data sets and the complete mitochondrial genome sequences, by MP (maximum parsimony), ML (maximum likelihood), and Bayesian analysis. The results showed that Panthera was composed of Panthera leo, P. uncia, P. pardus, Panthera onca, P. tigris, and N. nebulosa, which was included as the most basal member. The phylogeny within Panthera genus was N. nebulosa (P. tigris (P. onca (P. pardus, (P. leo, P. uncia)))). The divergence times for Panthera genus were estimated based on the ML branch lengths and four well-established calibration points. The results showed that at about 11.3 MYA, the Panthera genus separated from other felid species and then evolved into the several species of the genus. In detail, N. nebulosa was estimated to be founded about 8.66 MYA, P. tigris about 6.55 MYA, P. uncia about 4.63 MYA, and P. pardus about 4.35 MYA. All these estimated times were older than those estimated from the fossil records. The divergence event, evolutionary process, speciation, and distribution pattern of P. uncia, a species endemic to the central Asia with core habitats on the Qinghai-Tibetan Plateau and surrounding highlands, mostly correlated with the geological tectonic events and intensive climate shifts that happened at 8, 3.6, 2.5, and 1.7 MYA on the plateau during the late Cenozoic period.     

A fragment of chloroplast DNA was transferred horizontally,probably from non-eudicots,to mitochondrial genome of <Emphasis Type="Italic">Phaseolus</Emphasis>

The mitochondrial genomes of some Phaseolus species contain a fragment of chloroplast trnA gene intron, named pvs-trnA for its location within the Phaseolus vulgaris sterility sequence (pvs). The purpose of this study was to determine the type of transfer (intracellular or horizontal) that gave rise to pvs-trnA. Using a PCR approach we could not find the respective portion of the trnA gene as a part of pvs outside the Phaseolus genus. However, a BLAST search revealed longer fragments of trnA present in the mitochondrial genomes of some Citrus species, Helianthus annuus and Zea mays. Basing on the identity or near-identity between these mitochondrial sequences and their chloroplast counterparts we concluded that they had relocated from chloroplasts to mitochondria via recent, independent, intracellular DNA transfers. In contrast, pvs-trnA displayed a relatively higher sequence divergence when compared with its chloroplast counterpart from Phaseolus vulgaris. Alignment of pvs-trnA with corresponding trnAfragments from 35 plant species as well as phylogenetic analysis revealed that pvs-trnA grouped with non-eudicot sequences and was well separated from all Fabalessequences. In conclusion, we propose that pvs-trnA arose via horizontal transfer of a trnA intron fragment from chloroplast of a non-eudicot plant to Phaseolus mitochondria. This is the first example of horizontal transfer of a chloroplast sequence to the mitochondrial genome in higher plants.     

Mitochondrial DNA sequence-based phylogenetic relationship among flesh flies of the genus <Emphasis Type="Italic">Sarcophaga</Emphasis> (Sarcophagidae: Diptera)

The phylogenetic relationships among flesh flies of the family Sarcophagidae has been based mainly on the morphology of male genitalia. However, the male genitalic character-based relationships are far from satisfactory. Therefore, in the present study mitochondrial DNA has been used as marker to unravel genetic relatedness and to construct phylogeny among five sympatric species of the genus Sarcophaga. Two mitochondrial genes viz., cytochrome oxidase subunit 1 (COI) and NAD dehydrogenase subunit 5 (ND5) were sequenced and genetic distance values were calculated on the basis of sequence differences in both the mitochondrial genes. The data revealed very few genetic difference among the five species for the COI and ND5 gene sequences.     

The complete nucleotide sequence of the mitochondrial genome of <Emphasis Type="Italic">Phthonandria atrilineata</Emphasis> (Lepidoptera: Geometridae)

Using long-polymerase chain reaction (Long-PCR) method, we determined the complete nucleotide sequence of the mitochondrial genome (mitogenome) of Phthonandria atrilineata. The complete mtDNA from P. atrilineata was 15,499 base pairs in length and contained 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The P. atrilineata genes were in the same order and orientation as the completely sequenced mitogenomes of other lepidopteran species. The nucleotide composition of P. atrilineata mitogenome was biased toward A + T nucleotides (81.02%), and the 13 PCGs show different A + T contents that range from 73.25% (cox1) to 92.12% (atp8). Phthonandria had the canonical set of 22 tRNA genes, that fold in the typical cloverleaf structure described for metazoan mt tRNAs, with the unique exception of trnS(AGN). The phylogenetic relationships were reconstructed with the concatenated sequences of the 13 PCGs of the mitochondrial genome, which confirmed that P. atrilineata is most closely related to the superfamily Bombycoidea.     

REEXAMINATION OF PHYLOGENETIC RELATIONSHIPS WITHIN THE COLONIAL VOLVOCALES (CHLOROPHYTA): AN ANALYSIS OF atpB AND rbcL GENE SEQUENCES

The chloroplast-encoded atp B gene was sequenced from 33 strains representing 28 species of the colonial Volvocales (the Volvocaceae and its relatives) to reexamine phylogenetic relationships as previously deduced by morphological data and rbc L gene sequence data.1128 base pairs in the coding regions of the atp B gene were analyzed by MP, NJ, and ML analyses. Although supported with relatively low bootstrap values (75% and 65% in the NJ and ML analyses, respectively), three anisogamous/oogamous volvocacean genera— Eudorina, Pleodorina, and Volvox, excluding the section Volvox (= Euvolvox, illegitimate name), constituted a large monophyletic group (Eudorina group). Outside the Eudorina group, a robust lineage composed of three species of Volvox sect. Volvox was resolved as in the rbc L gene trees, rejecting the hypothesis of the previous cladistic analysis based on morphological data that the genus Volvox is monophyletic. In addition, the NJ and ML trees suggested that Eudorina is a nonmonophyletic genus as inferred from the morphological data and rbc L gene sequences. Although phylogenetic status of the genus Gonium is ambiguous in the rbc L gene trees and the paraphyly of this genus is resolved in the cladistic analysis based on morphological data, the atp B gene sequence data suggest monophyly of Gonium with relatively low bootstrap values (56–61%) in the NJ and ML trees. On the basis of the combined sequence data (2256 base pairs) from atp B and rbc L genes, Gonium was resolved as a robust monophyletic genus in the NJ and ML trees (with 68–86% bootstrap values), and Eudorina elegans Ehrenberg represented a paraphyletic species positioned most basally within the Eudorina group. However, phylogenetic status and relationships of the families of the colonial Volvocales were still almost ambiguous even in the combined analysis.     

Identification of mitochondrial genome rearrangements unique to novel cytoplasmic male sterility in radish (Raphanus sativus L.)

A novel cytoplasmic male-sterility (CMS) radish (Raphanus sativus L.) and its associated mitotype (DCGMS) were previously identified; however, no mtDNA fragments flanking the atp6 gene were found in the DCGMS mitotype. Unlike three other mitotypes in this study, a unique mtDNA organization, atp6–nad3–rps12, was found to be the major mtDNA structure associated with this mitotype. This organization may have arisen from short repeat sequence-mediated recombination events. The short repeat clusters involved in the mtDNA rearrangement around the atp6 gene also exist as repetitive sequences in the complete mitochondrial genomes of other members of the Brassicaceae family, including rapeseed and Arabidopsis. These sequences do not exist as repetitive elements in the mtDNA of tobacco, sugar beet, or rice. While studying the regions flanking atp6, we identified a truncated atp6 mtDNA fragment which consists of the 5′ part of the atp6 gene linked to an unidentified sequence. This mtDNA structure was present in all mitotypes; however, a single nucleotide insertion mutation leading to a frame-shift was identified only in the DCGMS mitotype. Although this truncated atp6 organization was transcribed, there was no significantly different expression between male-sterile and fertile segregating individuals from the BC₁F₁ population originating from a cross between male-sterile and restorer parents. Comprehensive survey of the single base-pair insertion showed that it was maternally inherited and unique to the DCGMS mitotype. Therefore, this single nucleotide polymorphism (SNP) in the coding sequence of the mtDNA will be a useful molecular marker for the detection of the DCGMS mitotype. Young-Pyo Lee and Sunggil Kim equally contributed to this work.     

Complete mitochondrial genomes of the Japanese pink coral (Corallium elatius) and the Mediterranean red coral (Corallium rubrum): a reevaluation of the phylogeny of the family Coralliidae based on molecular data

Precious corals are soft corals belonging to the family Coralliidae (Anthozoa: Octocorallia: Alcyonacea) and class Anthozoa, whose skeletal axes are used for jewelry.The family Coralliidae includes ca. 40 species and was originally thought to comprise of the single genus Corallium. In 2003, Corallium was split into two genera, Corallium and Paracorallium, and seven species were moved to this newly identified genus on the bases of morphological features. Previously, we determined the complete mitochondrial genome sequence of two precious corals Paracorallium japonicum and Corallium konojoi, in order to clarify their systematic positions. The two genomes showed high nucleotide sequence identity, but their gene order arrangements were not identical. Here, we determined three complete mitochondrial genome sequences from the one specimen of Mediterranean Corallium rubrum and two specimens of Corallium elatius coming from Kagoshima (South Japan). The circular mitochondrial genomes of C. rubrum and C. elatius are 18,915 bp and 18,969–18,970 bp in length, respectively, and encode 14 typical octocorallian protein-coding genes (nad1–6, nad4L, cox1–3, cob, atp6, atp8, and mtMutS, which is an octocoral-specific mismatch repair gene homologue), two ribosomal RNA genes (rns and rnl), and one transfer RNA (trnM). The overall nucleotide differences between C. konojoi and each C. elatius haplotype (T2007 and I2011) are only 10 and 11 nucleotides, respectively; this degree of similarity indicates that C. elatius and C. konojoi are very closely related species. Notably, the C. rubrum mitochondrial genome shows more nucleotide sequence identity to P. japonicum (99.5%) than to its congeneric species C. konojoi (95.3%) and C. elatius (95.3%). Moreover, the gene order arrangement of C. rubrum was the same as that of P. japonicum, while that of C. elatius was the same as C. konojoi. Phylogenetic analysis based on three mitochondrial genes from 24 scleraxonian species shows that the family Coralliidae is separated into two distinct groups, recovering Corallium as a paraphyletic genus. Our results indicate that the currently accepted generic classification of Coralliidae should be reconsidered.