Rapid evolution of enormous, multichromosomal genomes in flowering plant mitochondria with exceptionally high mutation rates

MitoSatPlant: mitochondrial microsatellites database of viridiplantae

Microsatellites also known as simple sequence repeats (SSRs) consist of 1-6 nucleotide long repeating units. The importance of mitochondrial SSRs (mtSSRs) in fields like population genetics, plant phylogenetics and genome mapping motivated us to develop MitoSatPlant, a repository of plant mtSSRs. It contains information for perfect, imperfect and compound SSRs mined from 92 mitochondrial genomes of green plants, available at NCBI (as of 1 Feb 2014). A total of 72,798 SSRs were found, of which PCR primers were designed for 72,495 SSRs. Among all sequences, tetranucleotide repeats (26,802) were found to be most abundant whereas hexanucleotide repeats (2751) were detected with least frequency. MitoSatPlant was developed using SQL server 2008 and can be accessed through a front end designed in ASP.Net. It is an easy to use, user-friendly database and will prove to be a useful resource for plant scientists. To the best of our knowledge MitoSatPlant is the only database available for plant mtSSRs and can be freely accessed at http://compubio.in/mitosatplant/.

Assembly and analysis of a male sterile rubber tree mitochondrial genome reveals DNA rearrangement events and a novel transcript

BACKGROUND

The rubber tree, Hevea brasiliensis, is an important plant species that is commercially grown to produce latex rubber in many countries. The rubber tree variety BPM 24 exhibits cytoplasmic male sterility, inherited from the variety GT 1.

RESULTS

We constructed the rubber tree mitochondrial genome of a cytoplasmic male sterile variety, BPM 24, using 454 sequencing, including 8 kb paired-end libraries, plus Illumina paired-end sequencing. We annotated this mitochondrial genome with the aid of Illumina RNA-seq data and performed comparative analysis. We then compared the sequence of BPM 24 to the contigs of the published rubber tree, variety RRIM 600, and identified a rearrangement that is unique to BPM 24 resulting in a novel transcript containing a portion of atp9.

CONCLUSIONS

The novel transcript is consistent with changes that cause cytoplasmic male sterility through a slight reduction to ATP production efficiency. The exhaustive nature of the search rules out alternative causes and supports previous findings of novel transcripts causing cytoplasmic male sterility.

Transcription could be the key to the selection advantage of mitochondrial deletion mutants in aging

The mitochondrial theory of aging is widely popular but confronted by several apparent inconsistencies. On the one hand, mitochondrial energy production is of central importance to the health and proper functioning of cells, and single-cell studies have shown that mtDNA deletion mutants accumulate in a clonal fashion in various mammalian species, displacing the wild-type mtDNAs. On the other hand, no explanation exists yet for the clonal expansion of mtDNA mutants that is compatible with experimental observations. We present here a new idea based on the distinctive connection between transcription and replication of metazoan mtDNA. Bioinformatic analysis of mtDNA deletion spectra strongly supports the predictions of this hypothesis and identifies specific candidates for proteins involved in transcriptional control of mtDNA replication. Computer simulations show the mechanism to be compatible with the available data from short- and long-lived mammalian species.

Idiosyncrasies in decoding mitochondrial genomes

Mitochondria originate from the α-proteobacterial domain of life. Since this unique event occurred, mitochondrial genomes of protozoans, fungi, plants and metazoans have highly derived and diverged away from the common ancestral DNA. These resulting genomes highly differ from one another, but all present-day mitochondrial DNAs have a very reduced coding capacity. Strikingly however, ATP production coupled to electron transport and translation of mitochondrial proteins are the two common functions retained in all mitochondrial DNAs. Paradoxically, most components essential for these two functions are now expressed from nuclear genes. Understanding how mitochondrial translation evolved in various eukaryotic models is essential to acquire new knowledge of mitochondrial genome expression. In this review, we provide a thorough analysis of the idiosyncrasies of mitochondrial translation as they occur between organisms. We address this by looking at mitochondrial codon usage and tRNA content. Then, we look at the aminoacyl-tRNA-forming enzymes in terms of peculiarities, dual origin, and alternate function(s). Finally we give examples of the atypical structural properties of mitochondrial tRNAs found in some organisms and the resulting adaptive tRNA-protein partnership.

Horizontal transfer of entire genomes via mitochondrial fusion in the angiosperm Amborella

We report the complete mitochondrial genome sequence of the flowering plant Amborella trichopoda. This enormous, 3.9-megabase genome contains six genome equivalents of foreign mitochondrial DNA, acquired from green algae, mosses, and other angiosperms. Many of these horizontal transfers were large, including acquisition of entire mitochondrial genomes from three green algae and one moss. We propose a fusion-compatibility model to explain these findings, with Amborella capturing whole mitochondria from diverse eukaryotes, followed by mitochondrial fusion (limited mechanistically to green plant mitochondria) and then genome recombination. Amborella's epiphyte load, propensity to produce suckers from wounds, and low rate of mitochondrial DNA loss probably all contribute to the high level of foreign DNA in its mitochondrial genome.

The plastid genomes of flowering plants

The plastid genome (plastome) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations are allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and tools of plastid genetic engineering.

Plant mitochondrial genome evolution can be explained by DNA repair mechanisms

Plant mitochondrial genomes are notorious for their large and variable size, nonconserved open reading frames of unknown function, and high rates of rearrangement. Paradoxically, the mutation rates are very low. However, mutation rates can only be measured in sequences that can be aligned—a very small part of plant mitochondrial genomes. Comparison of the complete mitochondrial genome sequences of two ecotypes of Arabidopsis thaliana allows the alignment of noncoding as well as coding DNA and estimation of the mutation rates in both. A recent chimeric duplication is also analyzed. A hypothesis is proposed that the mechanisms of plant mitochondrial DNA repair account for these features and includes different mechanisms in transcribed and nontranscribed regions. Within genes, a bias toward gene conversion would keep measured mutation rates low, whereas in noncoding regions, break-induced replication (BIR) explains the expansion and rearrangements. Both processes are types of double-strand break repair, but enhanced second-strand capture in transcribed regions versus BIR in nontranscribed regions can explain the two seemingly contradictory features of plant mitochondrial genome evolution—the low mutation rates in genes and the striking expansions of noncoding sequences.

A recurring syndrome of accelerated plastid genome evolution in the angiosperm tribe Sileneae (Caryophyllaceae)

In flowering plants, plastid genomes are generally conserved, exhibiting slower rates of sequence evolution than the nucleus and little or no change in structural organization. However, accelerated plastid genome evolution has occurred in scattered angiosperm lineages. For example, some species within the genus Silene have experienced a suite of recent changes to their plastid genomes, including inversions, shifts in inverted repeat boundaries, large indels, intron losses, and rapid rates of amino acid sequence evolution in a subset of protein genes, with the most extreme divergence occurring in the protease gene clpP. To investigate the relationship between the rates of sequence and structural evolution, we sequenced complete plastid genomes from three species (Silene conoidea, S. paradoxa, and Lychnis chalcedonica), representing independent lineages within the tribe Sileneae that were previously shown to have accelerated rates of clpP evolution. We found a high degree of parallel evolution. Elevated rates of amino acid substitution have occurred repeatedly in the same subset of plastid genes and have been accompanied by a recurring pattern of structural change, including cases of identical inversions and intron loss. This "syndrome" of changes was not observed in the closely related outgroup Agrostemma githago or in the more slowly evolving Silene species that were sequenced previously. Although no single mechanism has yet been identified to explain the correlated suite of changes in plastid genome sequence and structure that has occurred repeatedly in angiosperm evolution, we discuss a possible mixture of adaptive and non-adaptive forces that may be responsible.

Cytonuclear interactions and relaxed selection accelerate sequence evolution in organelle ribosomes

Many mitochondrial and plastid protein complexes contain subunits that are encoded in different genomes. In animals, nuclear-encoded mitochondrial proteins often exhibit rapid sequence evolution, which has been hypothesized to result from selection for mutations that compensate for changes in interacting subunits encoded in mutation-prone animal mitochondrial DNA. To test this hypothesis, we analyzed nuclear genes encoding cytosolic and organelle ribosomal proteins in flowering plants. The model angiosperm genus Arabidopsis exhibits low organelle mutation rates, typical of most plants. Nevertheless, we found that (nuclear-encoded) subunits of organelle ribosomes in Arabidopsis have higher amino acid sequence polymorphism and divergence than their counterparts in cytosolic ribosomes, suggesting that organelle ribosomes experience relaxed functional constraint. However, the observed difference between organelle and cytosolic ribosomes was smaller than in animals and could be partially attributed to rapid evolution in N-terminal organelle-targeting peptides that are not involved in ribosome function. To test the role of organelle mutation more directly, we used transcriptomic data from an angiosperm genus (Silene) with highly variable rates of organelle genome evolution. We found that Silene species with unusually fast-evolving mitochondrial and plastid DNA exhibited increased amino acid sequence divergence in ribosomal proteins targeted to the organelles but not in those that function in cytosolic ribosomes. Overall, these findings support the hypothesis that rapid organelle genome evolution has selected for compensatory mutations in nuclear-encoded proteins. We conclude that coevolution between interacting subunits encoded in different genomic compartments within the eukaryotic cell is an important determinant of variation in rates of protein sequence evolution.

Paternal leakage, heteroplasmy, and the evolution of plant mitochondrial genomes

Plant mitochondrial genomes are usually transmitted to the progeny from the maternal parent. However, cases of paternal transmission are known and are perhaps more common than once thought. This review will consider recent evidence, both direct and indirect, of paternal transmission (leakage) of the mitochondrial genome of seed plants, especially in natural populations, and how this can result in offspring that carry a mixture of maternally and paternally derived copies of the genome; a type of heteroplasmy. It will further consider how this heteroplasmy facilitates recombination between genetically distinct partners; a process that can enhance mitochondrial genotypic diversity. This will then form the basis for a discussion of five evolutionary questions that arise from these observations. Questions include how plant mitochondrial genome evolution can be placed on a sexual to asexual continuum, whether cytoplasmic male sterility (CMS) facilitates the evolution of paternal leakage, whether paternal leakage is more likely in populations undergoing admixture, how leakage influences patterns of gene flow, and whether heteroplasmy occurs in natural populations at a frequency greater than predicted by crossing experiments. It is proposed that each of these questions offers fertile ground for future research on a diversity of plant species.

One ring to rule them all? Genome sequencing provides new insights into the 'master circle' model of plant mitochondrial DNA structure

The in vivo molecular structure of plant mitochondrial DNA (mtDNA) has been a long-standing source of intrigue and controversy. Recent deep sequencing analyses of mitochondrial genomes from numerous plant species have provided the opportunity to revisit this decades-old question from a novel perspective. Whole-genome sequencing approaches have yielded new lines of evidence that the 'master circle' is not the predominant form of plant mtDNA and have revealed striking structural variation both within and among species. Here, I review these recent studies,including the discovery that at least two independent angiosperm lineages have evolved multichromosomal mitochondrial genome structures. These findings raise fascinating questions regarding the mechanisms of plant mtDNA replication and inheritance.

The plant mitochondrial genome: dynamics and maintenance

Plant mitochondria have a complex and peculiar genetic system. They have the largest genomes, as compared to organelles from other eukaryotic organisms. These can expand tremendously in some species, reaching the megabase range. Nevertheless, whichever the size, the gene content remains modest and restricted to a few polypeptides required for the biogenesis of the oxidative phosphorylation chain complexes, ribosomal proteins, transfer RNAs and ribosomal RNAs. The presence of autonomous plasmids of essentially unknown function further enhances the level of complexity. The physical organization of the plant mitochondrial DNA includes a set of sub-genomic forms resulting from homologous recombination between repeats, with a mixture of linear, circular and branched structures. This material is compacted into membrane-bound nucleoids, which are the inheritance units but also the centers of genome maintenance and expression. Recombination appears to be an essential characteristic of plant mitochondrial genetic processes, both in shaping and maintaining the genome. Under nuclear surveillance, recombination is also the basis for the generation of new mitotypes and is involved in the evolution of the mitochondrial DNA. In line with, or as a consequence of its complex physical organization, replication of the plant mitochondrial DNA is likely to occur through multiple mechanisms, potentially involving recombination processes. We give here a synthetic view of these aspects.

Complete mitochondrial genome sequence and identification of a candidate gene responsible for cytoplasmic male sterility in radish (Raphanus sativus L.) containing DCGMS cytoplasm

A novel cytoplasmic male sterility (CMS) conferred by Dongbu cytoplasmic and genic male-sterility (DCGMS) cytoplasm and its restorer-of-fertility gene (Rfd1) was previously reported in radish (Raphanus sativus L.). Its inheritance of fertility restoration and profiles of mitochondrial DNA (mtDNA)-based molecular markers were reported to be different from those of Ogura CMS, the first reported CMS in radish. The complete mitochondrial genome sequence (239,186 bp; GenBank accession No. KC193578) of DCGMS mitotype is reported in this study. Thirty-four protein-coding genes and three ribosomal RNA genes were identified. Comparative analysis of a mitochondrial genome sequence of DCGMS and previously reported complete sequences of normal and Ogura CMS mitotypes revealed various recombined structures of seventeen syntenic sequence blocks. Short-repeat sequences were identified in almost all junctions between syntenic sequence blocks. Phylogenetic analysis of three radish mitotypes showed that DCGMS was more closely related to the normal mitotype than to the Ogura mitotype. A single 1,551-bp unique region was identified in DCGMS mtDNA sequences and a novel chimeric gene, designated orf463, consisting of 128-bp partial sequences of cox1 gene and 1,261-bp unidentified sequences were found in the unique region. No other genes with a chimeric structure, a major feature of most characterized CMS-associated genes in other plant species, were found in rearranged junctions of syntenic sequence blocks. Like other known CMS-associated mitochondrial genes, the predicted gene product of orf463 contained 12 transmembrane domains. Thus, this gene product might be integrated into the mitochondrial membrane. In total, the results indicate that orf463 is likely to be a casual factor for CMS induction in radish containing the DCGMS cytoplasm.

Cytoplasmic male sterility-associated chimeric open reading frames identified by mitochondrial genome sequencing of four Cajanus genotypes

The hybrid pigeonpea (Cajanus cajan) breeding technology based on cytoplasmic male sterility (CMS) is currently unique among legumes and displays major potential for yield increase. CMS is defined as a condition in which a plant is unable to produce functional pollen grains. The novel chimeric open reading frames (ORFs) produced as a results of mitochondrial genome rearrangements are considered to be the main cause of CMS. To identify these CMS-related ORFs in pigeonpea, we sequenced the mitochondrial genomes of three C. cajan lines (the male-sterile line ICPA 2039, the maintainer line ICPB 2039, and the hybrid line ICPH 2433) and of the wild relative (Cajanus cajanifolius ICPW 29). A single, circular-mapping molecule of length 545.7 kb was assembled and annotated for the ICPA 2039 line. Sequence annotation predicted 51 genes, including 34 protein-coding and 17 RNA genes. Comparison of the mitochondrial genomes from different Cajanus genotypes identified 31 ORFs, which differ between lines within which CMS is present or absent. Among these chimeric ORFs, 13 were identified by comparison of the related male-sterile and maintainer lines. These ORFs display features that are known to trigger CMS in other plant species and to represent the most promising candidates for CMS-related mitochondrial rearrangements in pigeonpea.

The complete sequence of the mitochondrial genome of Butomus umbellatus--a member of an early branching lineage of monocotyledons

In order to study the evolution of mitochondrial genomes in the early branching lineages of the monocotyledons, i.e., the Acorales and Alismatales, we are sequencing complete genomes from a suite of key taxa. As a starting point the present paper describes the mitochondrial genome of Butomus umbellatus (Butomaceae) based on next-generation sequencing data. The genome was assembled into a circular molecule, 450,826 bp in length. Coding sequences cover only 8.2% of the genome and include 28 protein coding genes, four rRNA genes, and 12 tRNA genes. Some of the tRNA genes and a 16S rRNA gene are transferred from the plastid genome. However, the total amount of recognized plastid sequences in the mitochondrial genome is only 1.5% and the amount of DNA transferred from the nucleus is also low. RNA editing is abundant and a total of 557 edited sites are predicted in the protein coding genes. Compared to the 40 angiosperm mitochondrial genomes sequenced to date, the GC content of the Butomus genome is uniquely high (49.1%). The overall similarity between the mitochondrial genomes of Butomus and Spirodela (Araceae), the closest relative yet sequenced, is low (less than 20%), and the two genomes differ in size by a factor 2. Gene order is also largely unconserved. However, based on its phylogenetic position within the core alismatids Butomus will serve as a good reference point for subsequent studies in the early branching lineages of the monocotyledons.

OrganellarGenomeDRAW--a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets

Mitochondria and plastids (chloroplasts) are cell organelles of endosymbiotic origin that possess their own genetic information. Most organellar DNAs map as circular double-stranded genomes. Across the eukaryotic kingdom, organellar genomes display great size variation, ranging from ∼15 to 20 kb (the size of the mitochondrial genome in most animals) to >10 Mb (the size of the mitochondrial genome in some lineages of flowering plants). We have developed OrganellarGenomeDraw (OGDRAW), a suite of software tools that enable users to create high-quality visual representations of both circular and linear annotated genome sequences provided as GenBank files or accession numbers. Although all types of DNA sequences are accepted as input, the software has been specifically optimized to properly depict features of organellar genomes. A recent extension facilitates the plotting of quantitative gene expression data, such as transcript or protein abundance data, directly onto the genome map. OGDRAW has already become widely used and is available as a free web tool (http://ogdraw.mpimp-golm.mpg.de/). The core processing components can be downloaded as a Perl module, thus also allowing for convenient integration into custom processing pipelines.

Mitochondrial genomes as living 'fossils'

The huge variation between mitochondrial genomes makes untangling their evolutionary histories difficult. Richardson et al. report on the remarkably unaltered 'fossil' genome of the tulip tree, giving us many clues as to how the mitochondrial genomes of flowering plants have evolved over the last 150 million years, and raising questions about how such extraordinary sequence conservation can be maintained.

The "fossilized" mitochondrial genome of Liriodendron tulipifera: ancestral gene content and order, ancestral editing sites, and extraordinarily low mutation rate

BACKGROUND

The mitochondrial genomes of flowering plants vary greatly in size, gene content, gene order, mutation rate and level of RNA editing. However, the narrow phylogenetic breadth of available genomic data has limited our ability to reconstruct these traits in the ancestral flowering plant and, therefore, to infer subsequent patterns of evolution across angiosperms.

RESULTS

We sequenced the mitochondrial genome of Liriodendron tulipifera, the first from outside the monocots or eudicots. This 553,721 bp mitochondrial genome has evolved remarkably slowly in virtually all respects, with an extraordinarily low genome-wide silent substitution rate, retention of genes frequently lost in other angiosperm lineages, and conservation of ancestral gene clusters. The mitochondrial protein genes in Liriodendron are the most heavily edited of any angiosperm characterized to date. Most of these sites are also edited in various other lineages, which allowed us to polarize losses of editing sites in other parts of the angiosperm phylogeny. Finally, we added comprehensive gene sequence data for two other magnoliids, Magnolia stellata and the more distantly related Calycanthus floridus, to measure rates of sequence evolution in Liriodendron with greater accuracy. The Magnolia genome has evolved at an even lower rate, revealing a roughly 5,000-fold range of synonymous-site divergence among angiosperms whose mitochondrial gene space has been comprehensively sequenced.

CONCLUSIONS

Using Liriodendron as a guide, we estimate that the ancestral flowering plant mitochondrial genome contained 41 protein genes, 14 tRNA genes of mitochondrial origin, as many as 7 tRNA genes of chloroplast origin, >700 sites of RNA editing, and some 14 colinear gene clusters. Many of these gene clusters, genes and RNA editing sites have been variously lost in different lineages over the course of the ensuing ∽200 million years of angiosperm evolution.

The "fossilized" mitochondrial genome of Liriodendron tulipifera: ancestral gene content and order, ancestral editing sites, and extraordinarily low mutation rate

BACKGROUND

The mitochondrial genomes of flowering plants vary greatly in size, gene content, gene order, mutation rate and level of RNA editing. However, the narrow phylogenetic breadth of available genomic data has limited our ability to reconstruct these traits in the ancestral flowering plant and, therefore, to infer subsequent patterns of evolution across angiosperms.

RESULTS

We sequenced the mitochondrial genome of Liriodendron tulipifera, the first from outside the monocots or eudicots. This 553,721 bp mitochondrial genome has evolved remarkably slowly in virtually all respects, with an extraordinarily low genome-wide silent substitution rate, retention of genes frequently lost in other angiosperm lineages, and conservation of ancestral gene clusters. The mitochondrial protein genes in Liriodendron are the most heavily edited of any angiosperm characterized to date. Most of these sites are also edited in various other lineages, which allowed us to polarize losses of editing sites in other parts of the angiosperm phylogeny. Finally, we added comprehensive gene sequence data for two other magnoliids, Magnolia stellata and the more distantly related Calycanthus floridus, to measure rates of sequence evolution in Liriodendron with greater accuracy. The Magnolia genome has evolved at an even lower rate, revealing a roughly 5,000-fold range of synonymous-site divergence among angiosperms whose mitochondrial gene space has been comprehensively sequenced.

CONCLUSIONS

Using Liriodendron as a guide, we estimate that the ancestral flowering plant mitochondrial genome contained 41 protein genes, 14 tRNA genes of mitochondrial origin, as many as 7 tRNA genes of chloroplast origin, >700 sites of RNA editing, and some 14 colinear gene clusters. Many of these gene clusters, genes and RNA editing sites have been variously lost in different lineages over the course of the ensuing ∽200 million years of angiosperm evolution.

Transcription of atp1 is influenced by both genomic configuration and nuclear background in the highly rearranged mitochondrial genomes of Silene vulgaris

An extraordinary variation in mitochondrial DNA sequence exists in angiosperm Silene vulgaris. The atp1 gene is flanked by very variable regions, as deduced from four completely sequenced mitochondrial genomes of this species. This diversity contributed to a highly variable transcript profile of this gene observed across S. vulgaris populations. We examined the atp1 transcript in the KOV mitochondrial genome and found three 5' ends, created most likely by the combination of transcription initiation and RNA processing. Most atp1 transcripts terminated about 70 bp upstream of the translation stop codon, which was present in only 10 % of them. Controlled crosses between a KOV mother and a geographically distant pollen donor (Krasnoyarsk, Russia) showed that nuclear background also affected atp1 transcription. The distant pollen donor introduced the factor(s) preventing the formation of a long 2,100 nt-transcript, because this long atp1 transcript reappeared in the progeny from self-crosses. The highly rearranged mitochondrial genomes with a variation in gene flanking regions make S. vulgaris an excellent model for the study of mitochondrial gene expression in plants.

Mitochondrial and plastid genomes of the colonial green alga Gonium pectorale give insights into the origins of organelle DNA architecture within the volvocales

Volvocalean green algae have among the most diverse mitochondrial and plastid DNAs (mtDNAs and ptDNAs) from the eukaryotic domain. However, nearly all of the organelle genome data from this group are restricted to unicellular species, like Chlamydomonas reinhardtii, and presently only one multicellular species, the ∼4,000-celled Volvox carteri, has had its organelle DNAs sequenced. The V. carteri organelle genomes are repeat rich, and the ptDNA is the largest plastome ever sequenced. Here, we present the complete mtDNA and ptDNA of the colonial volvocalean Gonium pectorale, which is comprised of ∼16 cells and occupies a phylogenetic position closer to that of V. carteri than C. reinhardtii within the volvocine line. The mtDNA and ptDNA of G. pectorale are circular-mapping AT-rich molecules with respective lengths and coding densities of 16 and 222.6 kilobases and 73 and 44%. They share some features with the organelle DNAs of V. carteri, including palindromic repeats within the plastid compartment, but show more similarities with those of C. reinhardtii, such as a compact mtDNA architecture and relatively low organelle DNA intron contents. Overall, the G. pectorale organelle genomes raise several interesting questions about the origin of linear mitochondrial chromosomes within the Volvocales and the relationship between multicellularity and organelle genome expansion.

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.

The mitochondrial genome of Raphanus sativus and gene evolution of cruciferous mitochondrial types

To explore the mitochondrial genes of the Cruciferae family, the mitochondrial genome of Raphanus sativus (sat) was sequenced and annotated. The circular mitochondrial genome of sat is 239,723 bp and includes 33 protein-coding genes, three rRNA genes and 17 tRNA genes. The mitochondrial genome also contains a pair of large repeat sequences 5.9 kb in length, which may mediate genome reorganization into two sub-genomic circles, with predicted sizes of 124.8 kb and 115.0 kb, respectively. Furthermore, gene evolution of mitochondrial genomes within the Cruciferae family was analyzed using sat mitochondrial type (mitotype), together with six other reported mitotypes. The cruciferous mitochondrial genomes have maintained almost the same set of functional genes. Compared with Cycas taitungensis (a representative gymnosperm), the mitochondrial genomes of the Cruciferae have lost nine protein-coding genes and seven mitochondrial-like tRNA genes, but acquired six chloroplast-like tRNAs. Among the Cruciferae, to maintain the same set of genes that are necessary for mitochondrial function, the exons of the genes have changed at the lowest rates, as indicated by the numbers of single nucleotide polymorphisms. The open reading frames (ORFs) of unknown function in the cruciferous genomes are not conserved. Evolutionary events, such as mutations, genome reorganizations and sequence insertions or deletions (indels), have resulted in the non-conserved ORFs in the cruciferous mitochondrial genomes, which is becoming significantly different among mitotypes. This work represents the first phylogenic explanation of the evolution of genes of known function in the Cruciferae family. It revealed significant variation in ORFs and the causes of such variation.

Organelle genome complexity scales positively with organism size in volvocine green algae

It has been argued that for certain lineages, noncoding DNA expansion is a consequence of the increased random genetic drift associated with long-term escalations in organism size. But a lack of data has prevented the investigation of this hypothesis in most plastid-bearing protists. Here, using newly sequenced mitochondrial and plastid genomes, we explore the relationship between organelle DNA noncoding content and organism size within volvocine green algae. By looking at unicellular, colonial, and differentiated multicellular algae, we show that organelle DNA complexity scales positively with species size and cell number across the volvocine lineage. Moreover, silent-site genetic diversity data suggest that the volvocine species with the largest cell numbers and most bloated organelle genomes have the smallest effective population sizes. Together, these findings support the view that nonadaptive processes, like random genetic drift, promote the expansion of noncoding regions in organelle genomes.

Plant mitochondrial genome evolution can be explained by DNA repair mechanisms

Plant mitochondrial genomes are notorious for their large and variable size, nonconserved open reading frames of unknown function, and high rates of rearrangement. Paradoxically, the mutation rates are very low. However, mutation rates can only be measured in sequences that can be aligned--a very small part of plant mitochondrial genomes. Comparison of the complete mitochondrial genome sequences of two ecotypes of Arabidopsis thaliana allows the alignment of noncoding as well as coding DNA and estimation of the mutation rates in both. A recent chimeric duplication is also analyzed. A hypothesis is proposed that the mechanisms of plant mitochondrial DNA repair account for these features and includes different mechanisms in transcribed and nontranscribed regions. Within genes, a bias toward gene conversion would keep measured mutation rates low, whereas in noncoding regions, break-induced replication (BIR) explains the expansion and rearrangements. Both processes are types of double-strand break repair, but enhanced second-strand capture in transcribed regions versus BIR in nontranscribed regions can explain the two seemingly contradictory features of plant mitochondrial genome evolution--the low mutation rates in genes and the striking expansions of noncoding sequences.

Intraspecific variation in mitochondrial genome sequence, structure, and gene content in Silene vulgaris, an angiosperm with pervasive cytoplasmic male sterility

In angiosperms, mitochondrial-encoded genes can cause cytoplasmic male sterility (CMS), resulting in the coexistence of female and hermaphroditic individuals (gynodioecy). We compared four complete mitochondrial genomes from the gynodioecious species Silene vulgaris and found unprecedented amounts of intraspecific diversity for plant mitochondrial DNA (mtDNA). Remarkably, only about half of overall sequence content is shared between any pair of genomes. The four mtDNAs range in size from 361 to 429 kb and differ in gene complement, with rpl5 and rps13 being intact in some genomes but absent or pseudogenized in others. The genomes exhibit essentially no conservation of synteny and are highly repetitive, with evidence of reciprocal recombination occurring even across short repeats (< 250 bp). Some mitochondrial genes exhibit atypically high degrees of nucleotide polymorphism, while others are invariant. The genomes also contain a variable number of small autonomously mapping chromosomes, which have only recently been identified in angiosperm mtDNA. Southern blot analysis of one of these chromosomes indicated a complex in vivo structure consisting of both monomeric circles and multimeric forms. We conclude that S. vulgaris harbors an unusually large degree of variation in mtDNA sequence and structure and discuss the extent to which this variation might be related to CMS.

Mitochondrial evolution

Viewed through the lens of the genome it contains, the mitochondrion is of unquestioned bacterial ancestry, originating from within the bacterial phylum α-Proteobacteria (Alphaproteobacteria). Accordingly, the endosymbiont hypothesis--the idea that the mitochondrion evolved from a bacterial progenitor via symbiosis within an essentially eukaryotic host cell--has assumed the status of a theory. Yet mitochondrial genome evolution has taken radically different pathways in diverse eukaryotic lineages, and the organelle itself is increasingly viewed as a genetic and functional mosaic, with the bulk of the mitochondrial proteome having an evolutionary origin outside Alphaproteobacteria. New data continue to reshape our views regarding mitochondrial evolution, particularly raising the question of whether the mitochondrion originated after the eukaryotic cell arose, as assumed in the classical endosymbiont hypothesis, or whether this organelle had its beginning at the same time as the cell containing it.

The mitochondrial genome of Malus domestica and the import-driven hypothesis of mitochondrial genome expansion in seed plants

Mitochondrial genomes of spermatophytes are the largest of all organellar genomes. Their large size has been attributed to various factors; however, the relative contribution of these factors to mitochondrial DNA (mtDNA) expansion remains undetermined. We estimated their relative contribution in Malus domestica (apple). The mitochondrial genome of apple has a size of 396 947 bp and a one to nine ratio of coding to non-coding DNA, close to the corresponding average values for angiosperms. We determined that 71.5% of the apple mtDNA sequence was highly similar to sequences of its nuclear DNA. Using nuclear gene exons, nuclear transposable elements and chloroplast DNA as markers of promiscuous DNA content in mtDNA, we estimated that approximately 20% of the apple mtDNA consisted of DNA sequences imported from other cell compartments, mostly from the nucleus. Similar marker-based estimates of promiscuous DNA content in the mitochondrial genomes of other species ranged between 21.2 and 25.3% of the total mtDNA length for grape, between 23.1 and 38.6% for rice, and between 47.1 and 78.4% for maize. All these estimates are conservative, because they underestimate the import of non-functional DNA. We propose that the import of promiscuous DNA is a core mechanism for mtDNA size expansion in seed plants. In apple, maize and grape this mechanism contributed far more to genome expansion than did homologous recombination. In rice the estimated contribution of both mechanisms was found to be similar.

A complete mitochondrial genome sequence of Ogura-type male-sterile cytoplasm and its comparative analysis with that of normal cytoplasm in radish (Raphanus sativus L.)

BACKGROUND

Plant mitochondrial genome has unique features such as large size, frequent recombination and incorporation of foreign DNA. Cytoplasmic male sterility (CMS) is caused by rearrangement of the mitochondrial genome, and a novel chimeric open reading frame (ORF) created by shuffling of endogenous sequences is often responsible for CMS. The Ogura-type male-sterile cytoplasm is one of the most extensively studied cytoplasms in Brassicaceae. Although the gene orf138 has been isolated as a determinant of Ogura-type CMS, no homologous sequence to orf138 has been found in public databases. Therefore, how orf138 sequence was created is a mystery. In this study, we determined the complete nucleotide sequence of two radish mitochondrial genomes, namely, Ogura- and normal-type genomes, and analyzed them to reveal the origin of the gene orf138.

RESULTS

Ogura- and normal-type mitochondrial genomes were assembled to 258,426-bp and 244,036-bp circular sequences, respectively. Normal-type mitochondrial genome contained 33 protein-coding and three rRNA genes, which are well conserved with the reported mitochondrial genome of rapeseed. Ogura-type genomes contained same genes and additional atp9. As for tRNA, normal-type contained 17 tRNAs, while Ogura-type contained 17 tRNAs and one additional trnfM. The gene orf138 was specific to Ogura-type mitochondrial genome, and no sequence homologous to it was found in normal-type genome. Comparative analysis of the two genomes revealed that radish mitochondrial genome consists of 11 syntenic regions (length >3 kb, similarity >99.9%). It was shown that short repeats and overlapped repeats present in the edge of syntenic regions were involved in recombination events during evolution to interconvert two types of mitochondrial genome. Ogura-type mitochondrial genome has four unique regions (2,803 bp, 1,601 bp, 451 bp and 15,255 bp in size) that are non-syntenic to normal-type genome, and the gene orf138 was found to be located at the edge of the largest unique region. Blast analysis performed to assign the unique regions showed that about 80% of the region was covered by short homologous sequences to the mitochondrial sequences of normal-type radish or other reported Brassicaceae species, although no homology was found for the remaining 20% of sequences.

CONCLUSIONS

Ogura-type mitochondrial genome was highly rearranged compared with the normal-type genome by recombination through one large repeat and multiple short repeats. The rearrangement has produced four unique regions in Ogura-type mitochondrial genome, and most of the unique regions are composed of known Brassicaceae mitochondrial sequences. This suggests that the regions unique to the Ogura-type genome were generated by integration and shuffling of pre-existing mitochondrial sequences during the evolution of Brassicaceae, and novel genes such as orf138 could have been created by the shuffling process of mitochondrial genome.

Genome reduction and co-evolution between the primary and secondary bacterial symbionts of psyllids

Genome reduction in obligately intracellular bacteria is one of the most well-established patterns in the field of molecular evolution. In the extreme, many sap-feeding insects harbor nutritional symbionts with genomes that are so reduced that it is not clear how they perform basic cellular functions. For example, the primary symbiont of psyllids (Carsonella) maintains one of the smallest and most AT-rich bacterial genomes ever identified and has surprisingly lost many genes that are thought to be essential for its role in provisioning its host with amino acids. However, our understanding of this extreme case of genome reduction is limited, as genomic data for Carsonella are available from only a single host species, and little is known about the functional role of "secondary" bacterial symbionts in psyllids. To address these limitations, we analyzed complete Carsonella genomes from pairs of congeneric hosts in three divergent genera within the Psyllidae (Ctenarytaina, Heteropsylla, and Pachypsylla) as well as complete secondary symbiont genomes from two of these host species (Ctenarytaina eucalypti and Heteropsylla cubana). Although the Carsonella genomes are generally conserved in size, structure, and GC content and exhibit genome-wide signatures of purifying selection, we found that gene loss has remained active since the divergence of the host species and had a particularly large impact on the amino acid biosynthesis pathways that define the symbiotic role of Carsonella. In some cases, the presence of additional bacterial symbionts may compensate for gene loss in Carsonella, as functional gene content indicates a high degree of metabolic complementarity between co-occurring symbionts. The genomes of the secondary symbionts also show signatures of long-term evolution as vertically transmitted, intracellular bacteria, including more extensive genome reduction than typically observed in facultative symbionts. Therefore, a history of co-evolution with secondary bacterial symbionts can partially explain the ongoing genome reduction in Carsonella. However, the absence of these secondary symbionts in other host lineages indicates that the relationships are dynamic and that other mechanisms, such as changes in host diet or functional coordination with the host genome, must also be at play.

Gene conversion occurs within the mating-type locus of Cryptococcus neoformans during sexual reproduction

Meiotic recombination of sex chromosomes is thought to be repressed in organisms with heterogametic sex determination (e.g. mammalian X/Y chromosomes), due to extensive divergence and chromosomal rearrangements between the two chromosomes. However, proper segregation of sex chromosomes during meiosis requires crossing-over occurring within the pseudoautosomal regions (PAR). Recent studies reveal that recombination, in the form of gene conversion, is widely distributed within and may have played important roles in the evolution of some chromosomal regions within which recombination was thought to be repressed, such as the centromere cores of maize. Cryptococcus neoformans, a major human pathogenic fungus, has an unusually large mating-type locus (MAT, >100 kb), and the MAT alleles from the two opposite mating-types show extensive nucleotide sequence divergence and chromosomal rearrangements, mirroring characteristics of sex chromosomes. Meiotic recombination was assumed to be repressed within the C. neoformans MAT locus. A previous study identified recombination hot spots flanking the C. neoformans MAT, and these hot spots are associated with high GC content. Here, we investigated a GC-rich intergenic region located within the MAT locus of C. neoformans to establish if this region also exhibits unique recombination behavior during meiosis. Population genetics analysis of natural C. neoformans isolates revealed signals of homogenization spanning this GC-rich intergenic region within different C. neoformans lineages, consistent with a model in which gene conversion of this region during meiosis prevents it from diversifying within each lineage. By analyzing meiotic progeny from laboratory crosses, we found that meiotic recombination (gene conversion) occurs around the GC-rich intergenic region at a frequency equal to or greater than the meiotic recombination frequency observed in other genomic regions. We discuss the implications of these findings with regards to the possible functional and evolutionary importance of gene conversion within the C. neoformans MAT locus and, more generally, in fungi.

Recombination has been thought to be repressed within sex chromosomes, as well as within the mating-type (MAT) loci in many fungi, due to the highly diverged and rearranged nature between alleles defining opposite sexes or mating-types. However, it has long been appreciated that recombination can occur within these presumptive recombinational “cold spots,” and recent studies reveal that recombination, including gene conversion, can occur at a frequency higher than previously appreciated and could play important roles in shaping evolution of these chromosomal regions. Here, we provide evidence that, during sexual reproduction of the human pathogenic fungus Cryptococcus neoformans, recombination (gene conversion) occurs across a GC-rich intergenic region within the MAT locus. The frequency of this gene conversion is comparable to those of typical meiotic recombination events observed in other chromosomal regions. This is in accord with population genetics analyses, which indicate homogenization between alleles of opposite mating-types within the intergenic region. Gene conversion within these highly rearranged chromosomal regions may serve to ensure proper meiosis and/or rejuvenate genes/chromosomal regions within MAT that are otherwise facing irreversible evolutionary decay. In conclusion, our study provides further experimental evidence that at least some recombinational “cold spots” are not that cold, after all.

Relative rates of evolution among the three genetic compartments of the red alga Porphyra differ from those of green plants and do not correlate with genome architecture

In photosynthetic eukaryotes, relative silent-site nucleotide substitution rates (which can be used to approximate relative mutation rates) among mitochondrial, plastid, and nuclear genomes (mtDNAs, ptDNAs, and nucDNAs) are estimated to be 1:3:10 respectively for seed plants and roughly equal for green algae. These estimates correlate with certain genome characteristics, such as size and coding density, and have therefore been taken to support a relationship between mutation rate and genome architecture. Plants and green algae, however, represent a small fraction of the major eukaryotic plastid-bearing lineages. Here, we investigate relative rates of mutation within the model red algal genus Porphyra. In contrast to plants, we find that the levels of silent-site divergence between the Porphyra purpurea and Porphyra umbilicalis mtDNAs are three times that of their ptDNAs and five times that of their nucDNAs. Moreover, relative mutation rates do not correlate with genome architecture: despite an estimated three-fold difference in their mutation rate, the mitochondrial and plastid genome coding densities are equivalent - an observation that extends to organisms with secondary red algal plastids. These findings are supported by within-species silent-site polymorphism data from P. purpurea.

The mitochondrial genome of the lycophyte Huperzia squarrosa: the most archaic form in vascular plants

Mitochondrial genomes have maintained some bacterial features despite their residence within eukaryotic cells for approximately two billion years. One of these features is the frequent presence of polycistronic operons. In land plants, however, it has been shown that all sequenced vascular plant chondromes lack large polycistronic operons while bryophyte chondromes have many of them. In this study, we provide the completely sequenced mitochondrial genome of a lycophyte, from Huperzia squarrosa, which is a member of the sister group to all other vascular plants. The genome, at a size of 413,530 base pairs, contains 66 genes and 32 group II introns. In addition, it has 69 pseudogene fragments for 24 of the 40 protein- and rRNA-coding genes. It represents the most archaic form of mitochondrial genomes of all vascular plants. In particular, it has one large conserved gene cluster containing up to 10 ribosomal protein genes, which likely represents a polycistronic operon but has been disrupted and greatly reduced in the chondromes of other vascular plants. It also has the least rearranged gene order in comparison to the chondromes of other vascular plants. The genome is ancestral in vascular plants in several other aspects: the gene content resembling those of charophytes and most bryophytes, all introns being cis-spliced, a low level of RNA editing, and lack of foreign DNA of chloroplast or nuclear origin.

Recent acceleration of plastid sequence and structural evolution coincides with extreme mitochondrial divergence in the angiosperm genus Silene

The angiosperm genus Silene exhibits some of the most extreme and rapid divergence ever identified in mitochondrial genome architecture and nucleotide substitution rates. These patterns have been considered mitochondrial specific based on the absence of correlated changes in the small number of available nuclear and plastid gene sequences. To better assess the relationship between mitochondrial and plastid evolution, we sequenced the plastid genomes from four Silene species with fully sequenced mitochondrial genomes. We found that two species with fast-evolving mitochondrial genomes, S. noctiflora and S. conica, also exhibit accelerated rates of sequence and structural evolution in their plastid genomes. The nature of these changes, however, is markedly different from those in the mitochondrial genome. For example, in contrast to the mitochondrial pattern, which appears to be genome wide and mutationally driven, the plastid substitution rate accelerations are restricted to a subset of genes and preferentially affect nonsynonymous sites, indicating that altered selection pressures are acting on specific plastid-encoded functions in these species. Indeed, some plastid genes in S. noctiflora and S. conica show strong evidence of positive selection. In contrast, two species with more slowly evolving mitochondrial genomes, S. latifolia and S. vulgaris, have correspondingly low rates of nucleotide substitution in plastid genes as well as a plastid genome structure that has remained essentially unchanged since the origin of angiosperms. These results raise the possibility that common evolutionary forces could be shaping the extreme but distinct patterns of divergence in both organelle genomes within this genus.