The Mitochondrial Genome of Soybean Reveals Complex Genome Structures and Gene Evolution at Intercellular and Phylogenetic Levels

Determining mitochondrial genomes is important for elucidating vital activities of seed plants. Mitochondrial genomes are specific to each plant species because of their variable size, complex structures and patterns of gene losses and gains during evolution. This complexity has made research on the soybean mitochondrial genome difficult compared with its nuclear and chloroplast genomes. The present study helps to solve a 30-year mystery regarding the most complex mitochondrial genome structure, showing that pairwise rearrangements among the many large repeats may produce an enriched molecular pool of 760 circles in seed plants. The soybean mitochondrial genome harbors 58 genes of known function in addition to 52 predicted open reading frames of unknown function. The genome contains sequences of multiple identifiable origins, including 6.8 kb and 7.1 kb DNA fragments that have been transferred from the nuclear and chloroplast genomes, respectively, and some horizontal DNA transfers. The soybean mitochondrial genome has lost 16 genes, including nine protein-coding genes and seven tRNA genes; however, it has acquired five chloroplast-derived genes during evolution. Four tRNA genes, common among the three genomes, are derived from the chloroplast. Sizeable DNA transfers to the nucleus, with pericentromeric regions as hotspots, are observed, including DNA transfers of 125.0 kb and 151.6 kb identified unambiguously from the soybean mitochondrial and chloroplast genomes, respectively. The soybean nuclear genome has acquired five genes from its mitochondrial genome. These results provide biological insights into the mitochondrial genome of seed plants, and are especially helpful for deciphering vital activities in soybean.     

Mitochondrial genome structure of rice suspension culture from cytoplasmic male-sterile line (A-58CMS): reappraisal of the master circle

Summary The mitochondrial DNA (mtDNA) from the cultured cells of a cytoplasmic male-sterile line (A-58CMS) of rice (Oryza sativa) was cloned and its physical map was constructed. There was structural alteration on the mitochondrial genome during the cell culture. Detailed restriction analysis of cosmid clones having mtDNA fragments suggested either that the master genome has a 100-kb duplication (the genome size becomes 450 kb) or that a master circle is not present in the genome (the net structural complexity becomes 350 kb). The physical map of plant mitochondrial genomes thus far reported is illustrated in a single circle, namely a master circle. However, no circular DNA molecule corresponding to a master circle has yet been proved. In the present report, representation of plant mitochondrial genomes and a possibility for mitochondrial genome without a master circle are discussed.     

Rapid evolution of the compact and unusual mitochondrial genome in the ctenophore, Pleurobrachia bachei

Ctenophores are one of the most basally branching lineages of metazoans with the largest mitochondrial organelles in the animal kingdom. We sequenced the mitochondrial (mtDNA) genome from the Pacific cidipid ctenophore, Pleurobrachia bachei. The circular mitochondrial genome is 11,016 nts, with only 12 genes, and one of the smallest metazoan mtDNA genomes recorded. The protein coding genes are intronless cox1-3, cob, nad1, 3, 4, 4L and 5. The nad2 and 6 genes are represented as short fragments whereas the atp6 gene was found in the nuclear genome. Only the large ribosomal RNA subunit and two tRNAs were present with possibly the small subunit unidentifiable due to extensive fragmentation. The observed unique features of this mitochondrial genome suggest that nuclear and mitochondrial genomes have evolved at very different rates. This reduced mtDNA genome sharply contrasts with the very large sizes of mtDNA found in other basal metazoans including Porifera (sponges), and Placozoa (Trichoplax).     

Mitochondrial DNA of the blowflyPhormia regina: restriction analysis and gene localization

A study of an invertebrate mitochondrial genome, that of the blowflyPhormia regina, has been initiated to compare its structural and functional relatedness to other metazoan mitochondrial genomes. A restriction map of mitochondrial DNA (mtDNA) isolated from sucrose gradient-purified mitochondria has been established using a combination of single and double restriction endonuclease digestions and hybridizations with isolated mtDNA fragments, revealing a genome size of 17.5 kilobases (kb). A number of mitochondrial genes including those encoding the 12 S and 16 S ribosomal RNA, the cytochromec oxidase I subunit (COI) and an unidentified open reading frame (URF2) have been located on thePhormia mtDNA by Southern blot analysis using as probes both isolated mtDNA fragments and oligonucleotides derived from the sequences of previously characterized genes from rat andDrosophila yakuba mtDNAs. These data indicate that for those regions examined, the mitochondrial genome organization of blowfly mtDNA is the same as that ofDrosophila yakuba, the order being COI-URF2-12 S-16 S. These data also report the presence of an A + T-rich region, located as a 2.5-kb region between the URF2 and the 12 S rRNA genes, and its amplification by the polymerase chain reaction is described.     

Linear mitochondrial DNAs of yeasts: frequency of occurrence and general features.

In most yeast species, the mitochondrial DNA (mtDNA) has been reported to be a circular molecule. However, two cases of linear mtDNA with specific termini have previously been described. We examined the frequency of occurrence of linear forms of mtDNA among yeasts by pulsed-field gel electrophoresis. Among the 58 species from the genera Pichia and Williopsis that we examined, linear mtDNA was found with unexpectedly high frequency. Thirteen species contained a linear mtDNA, as confirmed by restriction mapping, and labeling, and electron microscopy. The mtDNAs from Pichia pijperi, Williopsis mrakii, and P. jadinii were studied in detail. In each case, the left and right terminal fragments shared homologous sequences. Between the terminal repeats, the order of mitochondrial genes was the same in all of the linear mtDNAs examined, despite a large variation of the genome size. This constancy of gene order is in contrast with the great variation of gene arrangement in circular mitochondrial genomes of yeasts. The coding sequences determined on several genes were highly homologous to those of the circular mtDNAs, suggesting that these two forms of mtDNA are not of distant origins.     

Increase of mitochondrial DNA copies with low level of DNA repair in tissue cells of gamma-irradiated mice

The damage and the change in the number of mitochondrial DNA (mtDNA) copies in brain and spleen tissues of gamma-irradiated mice were studied. The changes in the number of mitochondrial DNA (mtDNA) copies were assayed by the comparative analysis of the density values of long-extension PCR products of the mtDNA fragments (16 kb) and the cluster nuclear gene of beta-globin (8.7 kb). PCRs of mtDNA fragments and the nuclear gene of beta-globin were carried out simultaneously in one test-tube within total DNA. Our results showed that in brain and in spleen cells of mice exposed to gamma-radiation an increase in copy number (polyploidization) of mtDNA with regard to the nuclear gene beta-globin took place. The induction of polyploidization of mtDNA observed in cells of gamma-irradiated animals is regarded as the development of a compensatory reaction because of the energy deficiency due to the increased ATP consumption and structural alteration of genes controlling OXPHOS. The data enabled the assumption that because of the low efficiency of repair systems in mitochondria the induction of synthesis of new mtDNA copies on intact or little affected mtDNA templates may be the major mechanism for the retention of the mitochondrial genome which is constantly damaged by the endogenous ROS and is affected by ionizing radiation and/or other exogenous factors.     

Next generation molecular diagnosis of mitochondrial disorders

Mitochondrial disorders are by far the most genetically heterogeneous group of diseases, involving two genomes, the 16.6 kb mitochondrial genome and ~ 1500 genes encoded in the nuclear genome. For maternally inherited mitochondrial DNA disorders, a complete molecular diagnosis requires several different methods for the detection and quantification of mtDNA point mutations and large deletions. For mitochondrial disorders caused by autosomal recessive, dominant, and X-linked nuclear genes, the diagnosis has relied on clinical, biochemical, and molecular studies to point to a group of candidate genes followed by stepwise Sanger sequencing of the candidate genes one-by-one. The development of Next Generation Sequencing (NGS) has revolutionized the diagnostic approach. Using massively parallel sequencing (MPS) analysis of the entire mitochondrial genome, mtDNA point mutations and deletions can be detected and quantified in one single step. The NGS approach also allows simultaneous analyses of a group of genes or the whole exome, thus, the mutations in causative gene(s) can be identified in one-step. New approaches make genetic analyses much faster and more efficient. Huge amounts of sequencing data produced by the new technologies brought new challenges to bioinformatics, analytical pipelines, and interpretation of numerous novel variants. This article reviews the clinical utility of next generation sequencing for the molecular diagnoses of complex dual genome mitochondrial disorders.     

中华卵索线虫线粒体基因组多态性分析

为完善昆虫病原索科线虫线粒体基因组全序列数据库, 更系统地研究其基因组特征和系统演化规律, 进而为发挥该线虫生防潜力打下基础, 我们开展了中华卵索线虫Ovomermis sinensis线粒体全基因组的研究。该研究通过线粒体基因组滚环复制及酶切图谱, 揭示了中华卵索线虫线粒体基因组具有种内遗传多态性, 即群体中单体线虫具有独特的酶切条带, 且条带累加之和变化范围较大, 为16.5～24.5 kb。为进一步了解线粒体基因组多态性特征及产生的分子机制, 采用两步长PCR方法对2条代表性成虫线粒体基因组进行了测序及拼接, 得其全长分别为18 864和16 777 bp。对这2条序列的比对表明, 线粒体基因组中位于ND2和ND4之间的可变区域, 不仅基因排列顺序不同, 且存在ND3基因重复现象, 这是导致中华卵索线虫线粒体基因组呈现多态性的主要原因。通过对以上研究结果的分析及与GenBank中已有的6种索科线虫线粒体基因组序列进行比对, 概括出其线粒体基因组基本特点: ①线粒体基因排列顺序各不相同;②部分线虫线粒体基因存在重复现象, 且重复次数不同;③线粒体基因组大小存在很大差异。     

Extensive variation in nuclear mitochondrial DNA content between the genomes of Phytophthora sojae and Phytophthora ramorum

Fragments of mitochondrial DNA (mtDNA) transferred to the nuclear genome are called nuclear mitochondrial DNAs (NUMTs). We report here a comparison of NUMT content between genomes from two species of the same genus. Analysis of the genomes of Phytophthora sojae and P. ramorum revealed large differences in the NUMT content of the two genomes: 16.27 x 10(-3) and 2.28 x 10(-3)% of each genome, respectively. Substantial differences also exist between the two species in the sizes of the NUMTs found in each genome, with ranges of 20 to 405 bp for P. sojae and 19 to 137 bp for P. ramorum. Furthermore, in P. sojae, fragments from the mitochondrial genes rns, rnl, coxl, and nad (various subunits) are found most frequently, whereas P. ramorum NUMTs most often originate from the cox3, rpsl4, nad4, and nad5 genes. The large differences in the presumptive mtDNA insertions suggest that the insertions occurred subsequent to the divergence of the two species, and this is supported by sequence comparisons among the NUMTs and the mtDNA sequences of the two species. P. sojae mtDNA sequences inserted in the nuclear genome appear to have been altered as a result of insertions, deletions, inversions, and translocations and provide insights into active mechanisms of sequence divergence in this plant pathogen. No clear examples were found of NUMTs forming functional nuclear genes or of NUMTs inserted into exons or introns of any nuclear gene.     

Small, repetitive DNAs contribute significantly to the expanded mitochondrial genome of cucumber

Closely related cucurbit species possess eightfold differences in the sizes of their mitochondrial genomes. We cloned mitochondrial DNA (mtDNA) fragments showing strong hybridization signals to cucumber mtDNA and little or no signal to watermelon mtDNA. The cucumber mtDNA clones carried short (30-53 bp), repetitive DNA motifs that were often degenerate, overlapping, and showed no homology to any sequences currently in the databases. On the basis of dot-blot hybridizations, seven repetitive DNA motifs accounted for >13% (194 kb) of the cucumber mitochondrial genome, equaling >50% of the size of the Arabidopsis mitochondrial genome. Sequence analysis of 136 kb of cucumber mtDNA revealed only 11.2% with significant homology to previously characterized mitochondrial sequences, 2.4% to chloroplast DNA, and 15% to the seven repetitive DNA motifs. The remaining 71.4% of the sequence was unique to the cucumber mitochondrial genome. There was <4% sequence colinearity surrounding the watermelon and cucumber atp9 coding regions, and the much smaller watermelon mitochondrial genome possessed no significant amounts of cucumber repetitive DNAs. Our results demonstrate that the expanded cucumber mitochondrial genome is in part due to extensive duplication of short repetitive sequences, possibly by recombination and/or replication slippage.     

Ionizing radiation can activate the insertion of mitochondrial DNA fragments in the nuclear genome

In analytical review is considered the possibility of the insertion of mitochondrial DNA (mtDNA) fragments into the nuclear genome of cells, exposed ionizing radiation (IR). Many studies show that integration fragment mtDNA in nuclear genome, as well as its fastening as NUMT-pseudogenes, proceed at ancient periods of the evolutions not only, but also at more late periods. The number of the investigations shows that under influence endogenous reactive oxygen species, chemical agent, UV-light and IR mtDNA is damaged with greater frequency, than nucleus DNA. Furthermore, the repair systems in mitochondria are low efficiency. In irradiated by IR cells mtDNA fragments can transition from the mitochondria to the cytoplasm. The binding of mtDNA fragment to a complex with proteins provides them the protection from nuclease destroying. Possibly, at such safe condition they and are carried to nucleus. At inductions of DNA double-strand breaks (under the action of IR and activated their reparation) mtDNA fragments may be inserted to nuclear genome. Such integration of mtDNA to nuclear genome, with shaping NUMT-pseudogenes de novo, may be proceed in irradiated cells in the course of the reparations DNA double-strand breaks by the nonhomologous end-joining pathway. These insertions of mtDNA can cardinally change the structure of nuclear genomes in area of their introduction and render the essential influence upon the realization of genetic information. Available information in literature also allows to suppose that integration mtDNA in nuclear genome can proceed and at raised genomic instability observed in cells at post radiation period. It in equal extent pertains and to malignant cells with raised by instability mitochondrial and nuclear genomes. As the most efficient agent, initiating insertion fragment mtDNA in nuclear genome, is considered ionizing radiation.     

Phylogenetic information from three mitochondrial genomes of Terebelliformia (Annelida) worms and duplication of the methionine tRNA

Mitochondrial genomes have been useful for inferring animal phylogeny across a wide range of clades, however they are still poorly sampled in some animal taxa, limiting our knowledge of mtDNA evolution. For example, despite being one of the most diverse animal phyla, only 5 complete annelid mitochrondial genomes have been published. To address this paucity of information, we obtained complete mitochondrial genomic sequences from Pista cristata (Terebellidae) and Terebellides stroemi (Trichobranchidae) as well as one nearly complete mitochondrial genome from Eclysippe vanelli (Ampharetidae). These taxa are within Terebelliformia (Annelida), which include spaghetti worms, icecream cone worms and their relatives. In contrast to the 37 genes found in most bilaterian metazoans, we recover 38 genes in the mitochondrial genomes of T. stroemi and P. cristata due to the presence of a second methionine tRNA (trnM). Interestingly, the two trnMs are located next to each other and are possibly a synapomorphy of these two taxa. The E. vanelli partial mitochondrial genome lacks this additional trnM at the same position, but it may be present in the region not sampled. Compared to other annelids, gene orders of these three mitochondrial genomes are generally conserved except for the atp6-mSSU region. Phylogenetic analyses reveal that mtDNA data strongly supports a Trichobranchidae/Terebellidae clade.     

Organization and expression of the mitochondrial genome in the Nicotiana sylvestris CMSII mutant.

Previous analyses suggested that the Nicotiana sylvestris CMSII mutant carried a large deletion in its mitochondrial genome. Here, we show by cosmid mapping that the deletion is 60 kb in length and contains several mitochondrial genes or ORFs, including the complex I nad7 gene. However, due to the presence of large duplications in the progenitor mitochondrial genome, the only unique gene that appears to be deleted is nad7. RNA gel blot data confirm the absence of nad7 expression, strongly suggesting that the molecular basis for the CMSII abnormal phenotype, poor growth and male sterility, is the altered complex I structure. The CMSII mitochondrial genome appears to consist essentially of one of two subgenomes resulting from recombination between direct short repeats. In the progenitor mitochondrial genome both recombination products are detected by PCR and, reciprocally, the parental fragments are detected at the substoichiometric level in the mutant. The CMSII mtDNA organization has been maintained through six sexual generations.     

Mitochondrial DNA-like sequences in the nuclear genome of the opossum genus Didelphis (Marsupialia: Didelphidae).

This study reports the occurrence of mtDNA-like sequences in the nuclear genome of the opossum genus Didelphis (Didelphidae, Marsupialia). A specific primer pair designed to amplify a region encompassing a 3' terminal 118 bp region of the cytochrome b gene, the Thr and Pro tRNA genes, and a 489 bp region of the D-loop of the D. virginiana mtDNA, was used in highly stringent PCR reactions. These PCR reactions resulted in several fragments per individual varying in size from 259 bp to 1 kb. The sequencing of some of these fragments showed the occurrence of paralogous mtDNA-like sequences among the PCR amplified fragments. Analyses of qualitative aspects of these sequences, their transition/transversion ratios, and phylogenetic relationships were conclusive in showing the occurrence of mtDNA-like sequences in the nuclear genome of the genus Didelphis. Comparisons and phylogenetic analysis of orthologous mtDNA from the four Didelphis species and paralogous nuclear sequences suggested that mtDNA migration to the nuclear genome occurred more than once in Didelphis evolution.     

Variation of the mitochondrial genome in the evolution of Drosophila

The evidence on mitochondrial genome variation and its role in evolution of the genus Drosophila are reviewed. The mitochondrial genome is represented by a circular double-stranded DNA molecule 16 to 19 kb in length. The genome contains no introns involved in recombination. The entire mitochondrial genome can be arbitrarily divided into three parts: (1) protein-coding genes; (2) genes encoding rRNA and tRNA; and (3) the noncoding regulatory region (A + T region). The selective importance of mutations within different mtDNA regions is therefore unequal. In Drosophila, the content of the A + T pairs in mtDNA is extremely low and a pattern of nucleotide substitution is characterized by a low transition/transversion ratio (and a low threshold of mutation saturation). The deletions and duplications are of common occurrence in the mitochondrial genome. However, this genome lacks such characteristic for the nuclear genome aberrations as the inversions and transpositions. The phenomena of introgression and heteroplasmy provide an opportunity to study the adaptive role of the mitochondrial genome and its role in speciation. Analysis of evidence concerning mtDNA variation in different species of the genus Drosophila made it possible to ascertain data on phylogenetic relationships among species obtained by studying nuclear genome variation. In some species, mtDNA variation may serve as a reliable marker for population differentiation within a species, although evidence on the population dynamics of the mtDNA variation is very scarce.     

Variation of the Mitochondral Genome in the Evolution of Drosophila

The evidence on mitochondrial genome variation and its role in evolution of the genus Drosophila are reviewed. The mitochondrial genome is represented by a circular double-stranded DNA molecule 16 to 19 kb in length. Mitochondrial genes lack introns and recombination. The entire mitochondrial genome can be arbitrarily divided into three parts: (1) protein-coding genes; (2) genes encoding rRNA and tRNA; and (3) the noncoding regulatory region (A + T region). The selective importance of mutations within different mtDNA regions is therefore unequal. In Drosophila, the content of the A + T pairs in mtDNA is extremely high and a pattern of nucleotide substitution is characterized by a low transition/transversion ratio (and a low threshold of mutation saturation). The deletions and duplications are of common occurrence in the mitochondrial genome. However, this genome lacks such characteristic for the nuclear genome aberrations as inversions and transpositions. The phenomena of introgression and heteroplasmy provide an opportunity to study the adaptive role of the mitochondrial genome and its role in speciation. Analysis of evidence concerning mtDNA variation in different species of the genus Drosophilamade it possible to ascertain data on phylogenetic relationships among species obtained by studying nuclear genome variation. In some species, mtDNA variation may serve as a reliable marker for population differentiation within a species, although evidence on the population dynamics of the mtDNA variation is very scarce.     

Physical mapping of the mitochondrial genome of Arabidopsis thaliana by cosmid and YAC clones

As part of the worldwide efforts at molecular analysis of Arabidopsis thaliana as a model plant the complete structure of the mitochondrial genome has been determined. The mitochondrial DNA molecules were mapped by restriction fragment analysis of more than 300 cosmid clones and purified mitochondrial DNA. The entire genome of 372 kb is contained in three different configurations of circular molecules and is split into two additional subgenomic molecules of 234 kb and 138 kb, respectively. These arrangements result from recombinations of the two sets of repeats present in combinations of inverted and/or direct orientation. Alignment of YAC clones confirms the in vivo presence of continuous DNA molecules of more than 300 kb in A. thaliana mitochondria. The presence of this comparatively large mitochondrial genome in a plant with one of the smallest nuclear genomes shows that different size constraints act upon the different genomes in plant cells.     

Capturing single-copy nuclear genes,organellar genomes,and nuclear ribosomal DNA from deep genome skimming data for plant phylogenetics: A case study in Vitaceae

With the decreasing cost and availability of many newly developed bioinformatics pipelines, next-generation sequencing (NGS) has revolutionized plant systematics in recent years. Genome skimming has been widely used to obtain high-copy fractions of the genomes, including plastomes, mitochondrial DNA (mtDNA), and nuclear ribosomal DNA (nrDNA). In this study, through simulations, we evaluated the optimal (minimum) sequencing depth and performance for recovering single-copy nuclear genes (SCNs) from genome skimming data, by subsampling genome resequencing data and generating 10 data sets with different sequencing coverage in silico. We tested the performance of four data sets (plastome, nrDNA, mtDNA, and SCNs) obtained from genome skimming based on phylogenetic analyses of the Vitis clade at the genus level and Vitaceae at the family level, respectively. Our results showed that optimal minimum sequencing depth for high-quality SCNs assembly via genome skimming was about 10× coverage. Without the steps of synthesizing baits and enrichment experiments, coupled with incredibly low sequencing costs, we showcase that deep genome skimming (DGS) is as effective for capturing large data sets of SCNs as the widely used Hyb-Seq approach, in addition to capturing plastomes, mtDNA, and entire nrDNA repeats. DGS may serve as an efficient and economical alternative and may be superior to the popular target enrichment/Hyb-Seq approach.