Evolution of plant mitochondrial genomes via substoichiometric intermediates

The mitochondrial orf117Sha gene desynchronizes pollen development and causes pollen abortion in Arabidopsis Sha cytoplasmic male sterility

Cytoplasmic male sterility (CMS) is of major agronomical relevance in hybrid breeding. In gametophytic CMS, abortion of pollen is determined by the grain genotype, while in sporophytic CMS, it is determined by the mother plant genotype. While several CMS mechanisms have been dissected at the molecular level, gametophytic CMS has not been straightforwardly accessible. We used the gametophytic Sha-CMS in Arabidopsis to characterize the cause and process of pollen abortion by implementing in vivo biosensing in single pollen and mitoTALEN mutagenesis. We obtained conclusive evidence that orf117Sha is the CMS-causing gene, despite distinct characteristics from other CMS genes. We measured the in vivo cytosolic ATP content in single pollen, followed pollen development, and analyzed pollen mitochondrial volume in two genotypes that differed only by the presence of the orf117Sha locus. Our results showed that the Sha-CMS is not triggered by ATP deficiency. Instead, we observed desynchronization of a pollen developmental program. Pollen death occurred independently in pollen grains at diverse stages and was preceded by mitochondrial swelling. We conclude that pollen death is grain-autonomous in Sha-CMS and propose that mitochondrial permeability transition, which was previously described as a hallmark of developmental and environmental-triggered cell death programs, precedes pollen death in Sha-CMS.

A Systematic Review and Developmental Perspective on Origin of CMS Genes in Crops

Cytoplasmic male sterility (CMS) arises from the incompatibility between the nucleus and cytoplasm as typical representatives of the chimeric structures in the mitochondrial genome (mitogenome), which has been extensively applied for hybrid seed production in various crops. The frequent occurrence of chimeric mitochondrial genes leading to CMS is consistent with the mitochondrial DNA (mtDNA) evolution. The sequence conservation resulting from faithfully maternal inheritance and the chimeric structure caused by frequent sequence recombination have been defined as two major features of the mitogenome. However, when and how these chimeric mitochondrial genes appear in the context of the highly conserved reproduction of mitochondria is an enigma. This review, therefore, presents the critical view of the research on CMS in plants to elucidate the mechanisms of this phenomenon. Generally, distant hybridization is the main mechanism to generate an original CMS source in natural populations and in breeding. Mitochondria and mitogenomes show pleomorphic and dynamic changes at key stages of the life cycle. The promitochondria in dry seeds develop into fully functioning mitochondria during seed imbibition, followed by massive mitochondria or mitogenome fusion and fission in the germination stage along with changes in the mtDNA structure and quantity. The mitogenome stability is controlled by nuclear loci, such as the nuclear gene Msh1. Its suppression leads to the rearrangement of mtDNA and the production of heritable CMS genes. An abundant recombination of mtDNA is also often found in distant hybrids and somatic/cybrid hybrids. Since mtDNA recombination is ubiquitous in distant hybridization, we put forward a hypothesis that the original CMS genes originated from mtDNA recombination during the germination of the hybrid seeds produced from distant hybridizations to solve the nucleo-cytoplasmic incompatibility resulting from the allogenic nuclear genome during seed germination.

Cytoplasmic male sterility and abortive seed traits generated through mitochondrial genome editing coupled with allotopic expression of atp1 in tobacco

Allotopic expression is the term given for the deliberate relocation of gene function from an organellar genome to the nuclear genome. We hypothesized that the allotopic expression of an essential mitochondrial gene using a promoter that expressed efficiently in all cell types except those responsible for male reproduction would yield a cytoplasmic male sterility (CMS) phenotype once the endogenous mitochondrial gene was inactivated via genome editing. To test this, we repurposed the mitochondrially encoded atp1 gene of tobacco to function in the nucleus under the transcriptional control of a CaMV 35S promoter (construct 35S:nATP1), a promoter that has been shown to be minimally expressed in early stages of anther development. The endogenous atp1 gene was eliminated (Δatp1) from 35S:nATP1 tobacco plants using custom-designed meganucleases directed to the mitochondria. Vegetative growth of most 35S:nATP1/Δatp1 plants appeared normal, but upon flowering produced malformed anthers that failed to shed pollen. When 35S:nATP1/Δatp1 plants were cross-pollinated, ovary/capsule development appeared normal, but the vast majority of the resultant seeds were small, largely hollow and failed to germinate, a phenotype akin to the seedless trait known as stenospermocarpy. Characterization of the mitochondrial genomes from three independent Δatp1 events suggested that spontaneous recombination over regions of microhomology and substoichiometric shifting were the mechanisms responsible for atp1 elimination and genome rearrangement in response to exposure to the atp1-targeting meganucleases. Should the results reported here in tobacco prove to be translatable to other crop species, then multiple applications of allotopic expression of an essential mitochondrial gene followed by its elimination through genome editing can be envisaged. Depending on the promoter(s) used to drive the allotopic gene, this technology may have potential application in the areas of: (1) CMS trait development for use in hybrid seed production; (2) seedless fruit production; and (3) transgene containment.

RADA-dependent branch migration has a predominant role in plant mitochondria and its defect leads to mtDNA instability and cell cycle arrest

Mitochondria of flowering plants have large genomes whose structure and segregation are modulated by recombination activities. The post-synaptic late steps of mitochondrial DNA (mtDNA) recombination are still poorly characterized. Here we show that RADA, a plant ortholog of bacterial RadA/Sms, is an organellar protein that drives the major branch-migration pathway of plant mitochondria. While RadA/Sms is dispensable in bacteria, RADA-deficient Arabidopsis plants are severely impacted in their development and fertility, correlating with increased mtDNA recombination across intermediate-size repeats and accumulation of recombination-generated mitochondrial subgenomes. The radA mutation is epistatic to recG1 that affects the additional branch migration activity. In contrast, the double mutation radA recA3 is lethal, underlining the importance of an alternative RECA3-dependent pathway. The physical interaction of RADA with RECA2 but not with RECA3 further indicated that RADA is required for the processing of recombination intermediates in the RECA2-depedent recombination pathway of plant mitochondria. Although RADA is dually targeted to mitochondria and chloroplasts we found little to no effects of the radA mutation on the stability of the plastidial genome. Finally, we found that the deficient maintenance of the mtDNA in radA apparently triggers a retrograde signal that activates nuclear genes repressing cell cycle progression.

In flowering plants, the mitochondrial genome is very large and dynamic, and its stability influences plant fitness and development. Rearrangements by recombination drive its very rapid evolution and can lead to valuable agronomic traits such as cytoplasmic sterility, used by breeders for the production of hybrid seeds. Here we describe RADA, a DNA helicase essential for the stability of the mitochondrial DNA in Arabidopsis. We demonstrate that RADA has branch migrating activity, accelerating the processing of recombination intermediates. radA mutants are severely affected in development and fertility. They display mitochondrial genome instability that results in uncoordinated replication of subgenomes created by recombination. Furthermore, we found that an important component of the growth defects of radA mutants is apparently a cellular response triggered by the sensing of damages to the mitochondrial genome, resulting in the activation of genes that suppress the progression of the cell cycle. Our results underline the importance of better understanding the plant mitochondrial recombination pathways and their cross-talk with nuclear gene expression.

TALEN-mediated depletion of the mitochondrial gene orf312 proves that it is a Tadukan-type cytoplasmic male sterility-causative gene in rice

Cytoplasmic male sterility (CMS) is a trait that causes pollen or anther dysfunctions, resulting in the lack of seed setting. CMS is considered to be caused by the expression of a unique mitochondrial open reading frame referred to as CMS-associated gene. orf312 has been reported as a CMS-associated gene of Tadukan-type CMS (TAA) in rice (Oryza sativa L.), which exhibits impaired anther dehiscence; however, evidence thereof has not yet been reported. Here, we took a loss-of-function approach, using a mitochondria-targeted transcription activator-like effector nuclease (mitoTALEN) designed to knock out orf312 in TAA, to prove that orf312 indeed is a CMS-causative gene. Out of 28 transgenic TAA plants harboring the mitoTALEN expression vector, deletion of orf312 was detected in 24 plants by PCR, Southern blot, and sequencing analyses. The 24 plants were grouped into three groups based on the deleted regions. All orf312-depleted TAA plants exhibited recovery of anther dehiscence and seed setting. The depletion of orf312 and fertility restoration was maintained in the next generation, even in mitoTALEN expression cassette null segregants. In contrast, orf312-retaining plants were sterile. These results provide robust evidence that orf312 is a Tadukan-type CMS-causative gene.

Structure and Phylogeny of the Curly Birch Chloroplast Genome

Curly birch [Betula pendula var. carelica (Merckl.) Hämet-Ahti] is a relatively rare variety of silver birch (B. pendula Roth) that occurs mainly in Northern Europe and northwest part of Russia (Karelia). It is famous for the beautiful decorative texture of wood. Abnormal xylogenesis underlying this trait is heritable, but its genetic mechanism has not yet been fully understood. The high number of potentially informative genetic markers can be identified through sequencing nuclear and organelle genomes. Here, the de novo assembly, complete nucleotide sequence, and annotation of the chloroplast genome (plastome) of curly birch are presented for the first time. The complete plastome length is 160,523 bp. It contains 82 genes encoding structural and enzymatic proteins, 37 transfer RNAs (tRNAs), and eight ribosomal RNAs (rRNAs). The chloroplast DNA (cpDNA) is AT-rich containing 31.5% of A and 32.5% of T nucleotides. The GC-rich regions represent inverted repeats IR1 and IR2 containing genes of rRNAs (5S, 4.5S, 23S, and 16S) and tRNAs (trnV, trnI, and trnA). A high content of GC was found in rRNA (55.2%) and tRNA (53.2%) genes, but only 37.0% in protein-coding genes. In total, 384 microsatellite or simple sequence repeat (SSR) loci were found, mostly with mononucleotide motifs (92% of all loci) and predominantly A or T motifs (94% of all mononucleotide motifs). Comparative analysis of cpDNA in different plant species revealed high structural and functional conservatism in organization of the angiosperm plastomes, while the level of differences depends on the phylogenetic relationship. The structural and functional organization of plastome in curly birch was similar to cpDNA in other species of woody plants. Finally, the identified cpDNA sequence variation will allow to develop useful genetic markers.

Mitochondrial genes are involved in the fertility transformation of the thermosensitive male-sterile line YS3038 in wheat

Heterosis can improve the stress resistance, quality, and yield of crops, and the male sterility of wheat can be utilized to accelerate the breeding process of hybrid. To determine whether mitochondrial genes are involved in the fertility of K-type cytoplasmic male-sterile (CMS) line and the YS-type thermosensitive male-sterile (TMS) line in wheat, we sequenced and assembled the mitochondrial genomes of K519A, 519B, and YS3038 by next-generation sequencing (NGS). The non-synonymous mutations were analyzed, and the first-generation sequencing was conducted to verify the non-synonymous mutation sites. Furthermore, the expression patterns of genes with non-synonymous mutations were analyzed. Finally, the candidate genes were silenced by barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) to test the functions of the candidate genes. The results revealed that the mitochondrial genomes of K519A, 519B, and YS3038 were 420,543, 433,560, and 452,567 bp in length, respectively. Besides, 33, 31, and 37 protein-coding genes were identified in K519A, 519B, and YS3038, respectively. There were 14 protein-coding genes and 83 open reading frame (ORF) sequences that differed between K519A and 519B and 10 protein-coding genes and 122 ORF sequences that differed between K519A and YS3038. At the binucleate stage, seven genes (nad6, ORF256, ORF216, ORF138, atp6, nad3, and cox1) were downregulated in K519A compared with 519B, and 10 genes (nad6, atp6, cox3, atp8, nad3, cox1, rps3, ORF216, ORF138, and ORF224) were downregulated in YS3038 compared with K519A. Besides, six genes (nad6, ORF138, cox3, cox1, rps3, and ORF224) were downregulated under fertile conditions relative to sterile conditions in YS3038. Gene silencing analysis showed that the silencing of cox1 significantly reduced the seed setting rate of YS3038, indicating that the cox1 gene may be involved in the fertility transformation of YS3038.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01252-x.

Assessment of Mitochondrial DNA Copy Number, Stability, and Repair in Arabidopsis

Mitochondrial functions depend on the proper maintenance and expression of the mitochondrial genome (mtDNA). Therefore, understanding mtDNA replication and repair requires methods to assess its integrity. Mutations or chemical treatments that affect processes involved in the maintenance or stability of the mtDNA can affect its global copy number, but also the relative abundance of different genomic regions or the frequency of illegitimate recombination across repeated sequences. These can be conveniently tested by quantitative PCR (qPCR). Arabidopsis thaliana offers several advantages for studying these processes, because of the extensive collections of mutants, natural accessions and other genetic resources available from stock centers. Here we describe protocols we routinely use to explore changes in mtDNA copy number and relative stoichiometry in Arabidopsis mutants of genes involved in the replication, repair and recombination of the mtDNA.

In and out: Evolution of viral sequences in the mitochondrial genomes of legumes (Fabaceae)

Plant specific mitoviruses in the 'genus' Mitovirus (Narnaviridae) and their integrated sequences (non-retroviral endogenous RNA viral elements or NERVEs) have been recently identified in various plant lineages. However, the sparse phylogenetic coverage of complete plant mitochondrial genome (mitogenome) sequences and the non-conserved nature of mitochondrial intergenic regions have hindered comparative studies on mitovirus NERVEs in plants. In this study, 10 new mitogenomes were sequenced from legumes (Fabaceae). Based on comparative genomic analysis of 27 total mitogenomes, we identified mitovirus NERVEs and transposable elements across the family. All legume mitogenomes included NERVEs and total NERVE length varied from ca. 2 kb in the papilionoid Trifolium to 35 kb in the mimosoid Acacia. Most of the NERVE integration sites were in highly variable intergenic regions, however, some were positioned in six cis-spliced mitochondrial introns. In the Acacia mitogenome, there were L1-like transposon sequences including an almost full-length copy with target site duplications (TSDs). The integration sites of NERVEs in four introns showed evidence of L1-like retrotransposition events. Phylogenetic analysis revealed that there were multiple instances of precise deletion of NERVEs between TSDs. This study provides clear evidence that a L1-like retrotransposition mechanism has a long history of contributing to the integration of viral RNA into plant mitogenomes while microhomology-mediated deletion can restore the integration site.

Effects of mitoTALENs-Directed Double-Strand Breaks on Plant Mitochondrial Genomes

Mitochondrial genomes in flowering plants differ from those in animals and yeasts in several ways, including having large and variable sizes, circular, linear and branched structures, long repeat sequences that participate in homologous recombinations, and variable genes orders, even within a species. Understanding these differences has been hampered by a lack of genetic methods for transforming plant mitochondrial genomes. We recently succeeded in disrupting targeted genes in mitochondrial genomes by mitochondria-targeted transcription activator-like effector nucleases (mitoTALENs) in rice, rapeseed, and Arabidopsis. Double-strand breaks created by mitoTALENs were repaired not by non-homologous end-joining (NHEJ) but by homologous recombination (HR) between repeats near and far from the target sites, resulting in new genomic structures with large deletions and different configurations. On the other hand, in mammals, TALENs-induced DSBs cause small insertions or deletions in nuclear genomes and degradation of mitochondrial genomes. These results suggest that the mitochondrial and nuclear genomes of plants and mammals have distinct mechanisms for responding to naturally occurring DSBs. The different responses appear to be well suited to differences in size and copy numbers of each genome.

Comparative analysis of mitochondrial genomes of soybean cytoplasmic male-sterile lines and their maintainer lines

In soybean, only one mitochondrial genome of cultispecies has been completely obtained. To explore the effect of mitochondrial genome on soybean cytoplasmic male sterility (CMS), two CMS lines and three maintainer lines were used for sequencing. Comparative analysis showed that mitochondrial genome of the CMS line was more compact than that of its maintainer line, but genes were highly conserved. Conserved and unique sequence coexisted in the genomes. Mitochondrial genomes contained different sequence lengths and copy numbers of repeats between CMS line and maintainer line. Large and short repeats mediated intramolecular and intermolecular recombination in mitochondria. Unique sequences and genes were also involved in recombination process and constituted a complex network. orf178 and orf261 were identified as CMS-associated candidate genes. They had sequence characteristics of reported CMS genes in other crops and could be transcribed in CMS lines but not in maintainer lines. This report reveals mitochondrial genome of soybean CMS lines and compares complete mitochondrial sequence between CMS lines and their maintainer lines. The information will be helpful in further understanding the characteristics of soybean mitochondrial genome and the mechanism underlying CMS.

DNA Repair and the Stability of the Plant Mitochondrial Genome

The mitochondrion stands at the center of cell energy metabolism. It contains its own genome, the mtDNA, that is a relic of its prokaryotic symbiotic ancestor. In plants, the mitochondrial genetic information influences important agronomic traits including fertility, plant vigor, chloroplast function, and cross-compatibility. Plant mtDNA has remarkable characteristics: It is much larger than the mtDNA of other eukaryotes and evolves very rapidly in structure. This is because of recombination activities that generate alternative mtDNA configurations, an important reservoir of genetic diversity that promotes rapid mtDNA evolution. On the other hand, the high incidence of ectopic recombination leads to mtDNA instability and the expression of gene chimeras, with potential deleterious effects. In contrast to the structural plasticity of the genome, in most plant species the mtDNA coding sequences evolve very slowly, even if the organization of the genome is highly variable. Repair mechanisms are probably responsible for such low mutation rates, in particular repair by homologous recombination. Herein we review some of the characteristics of plant organellar genomes and of the repair pathways found in plant mitochondria. We further discuss how homologous recombination is involved in the evolution of the plant mtDNA.

Curing cytoplasmic male sterility via TALEN-mediated mitochondrial genome editing

Sequence-specific nucleases are commonly used to modify the nuclear genome of plants. However, targeted modification of the mitochondrial genome of land plants has not yet been achieved. In plants, a type of male sterility called cytoplasmic male sterility (CMS) has been attributed to certain mitochondrial genes, but none of these genes has been validated by direct mitochondrial gene-targeted modification. Here, we knocked out CMS-associated genes (orf79 and orf125) of CMS varieties of rice and rapeseed, respectively, using transcription activator-like effector nucleases (TALENs) with mitochondria localization signals (mitoTALENs). We demonstrate that knocking out these genes cures male sterility, strongly suggesting that these genes are causes of CMS. Sequencing revealed that double-strand breaks induced by mitoTALENs were repaired by homologous recombination, and that during this process, the target genes and surrounding sequences were deleted. Our results show that mitoTALENs can be used to stably and heritably modify the mitochondrial genome in plants.

Elucidating genomic patterns and recombination events in plant cybrid mitochondria

KEY MESSAGE

Cybrid plant mitochondria undergo homologous recombination, mainly BIR, keep a single allele for each gene, and maintain exclusive sequences of each parent and a single copy of the homologous regions. The maintenance of a dynamic equilibrium between the mitochondrial and nuclear genomes requires continuous communication and a high level of compatibility between them, so that alterations in one genetic compartment need adjustments in the other. The co-evolution of nuclear and mitochondrial genomes has been poorly studied, even though the consequences and effects of this interaction are highly relevant for human health, as well as for crop improvement programs and for genetic engineering. The mitochondria of plants represent an excellent system to understand the mechanisms of genomic rearrangements, chimeric gene formation, incompatibility between nucleus and cytoplasm, and horizontal gene transfer. We carried out detailed analyses of the mtDNA of a repeated cybrid between the solanaceae Nicotiana tabacum and Hyoscyamus niger. The mtDNA of the cybrid was intermediate between the size of the parental mtDNAs and the sum of them. Noticeably, most of the homologous sequences inherited from both parents were lost. In contrast, the majority of the sequences exclusive of a single parent were maintained. The mitochondrial gene content included a majority of N. tabacum derived genes, but also chimeric, two-parent derived, and H. niger-derived genes in a tobacco nuclear background. Any of these alterations in the gene content could be the cause of CMS in the cybrid. The parental mtDNAs interacted through 28 homologous recombination events and a single case of illegitimate recombination. Three main homologous recombination mechanisms were recognized in the cybrid mitochondria. Break induced replication (BIR) pathway was the most frequent. We propose that BIR could be one of the mechanisms responsible for the loss of the majority of the repeated regions derived from H. niger.

The complete organelle genomes of Physochlaina orientalis: Insights into short sequence repeats across seed plant mitochondrial genomes

Short repeats (SR) play an important role in shaping seed plant mitochondrial genomes (mtDNAs). However, their origin, distribution, and relationships across the different plant lineages remain unresolved. We focus on the angiosperm family Solanaceae that shows great variation in repeat content and extend the study to a wide diversity of seed plants. We determined the complete nucleotide sequences of the organellar genomes of the medicinal plant Physochlaina orientalis (Solanaceae), member of the tribe Hyoscyameae. To understand the evolution of the P. orientalis mtDNA we made comparisons with those of five other Solanaceae. P. orientalis mtDNA presents the largest mitogenome (∼685 kb in size) among the Solanaceae and has an unprecedented 8-copy repeat family of ∼8.2 kb in length and a great number of SR arranged in tandem-like structures. We found that the SR in the Solanaceae share a common origin, but these only expanded in members of the tribe Hyoscyameae. We discuss a mechanism that could explain SR formation and expansion in P. orientalis and Hyoscyamus niger. Finally, the great increase in plant mitochondrial data allowed us to systematically extend our repeat analysis to a total of 136 seed plants to characterize and analyze for the first time families of SR among seed plant mtDNAs.

Identification of a gene responsible for cytoplasmic male-sterility in onions (Allium cepa L.) using comparative analysis of mitochondrial genome sequences of two recently diverged cytoplasms

Almost identical mitochondrial genome sequences of two recently diverged male-fertile normal and male-sterile CMS-T-like cytoplasms were obtained in onions. A chimeric gene, orf725 , was found to be a CMS-inducing gene. In onions (Allium cepa L.), cytoplasmic male-sterility (CMS) has been widely used in hybrid seed production. Two types of CMS (CMS-S and CMS-T) have been reported in onions. A complete mitochondrial genome sequence of the CMS-S cytoplasm has been reported in our previous study. Draft mitochondrial genome sequences of male-fertile normal and CMS-T-like cytoplasms are reported in this study. Raw reads obtained from normal and CMS-T-like cytoplasms were assembled into eight and nine almost identical contigs, respectively. After connection and reorganization of contigs by PCR amplification and genome walking, four scaffold sequences with total length of 339 and 180 bp were produced for the normal cytoplasm. A mitochondrial genome sequence of the CMS-T-like cytoplasm was obtained by mapping trimmed reads of CMS-T onto scaffold sequences of the normal cytoplasm. Compared with the CMS-S mitochondrial genome, the normal mitochondrial genome was highly rearranged with 31 syntenic blocks. A total of 499 single nucleotide polymorphisms (SNPs) or insertions/deletions were identified in these syntenic regions. On the other hand, normal and CMS-T-like mitochondrial genome sequences were almost identical except for orf725, a chimeric gene consisting of cox1 with other sequences. Only three SNPs were identified between normal and CMS-T-like syntenic sequences. These results indicate that orf725 is likely to be the casual gene for CMS induction in onions and that CMS-T-like cytoplasm has recently diverged from the normal cytoplasm by introduction of orf725.

Candidate Gene Selection for Cytoplasmic Male Sterility in Pepper (Capsicum annuum L.) through Whole Mitochondrial Genome Sequencing

Cytoplasmic male sterility (CMS), which is controlled by mitochondrial genes, is an important trait for commercial hybrid seed production. So far, genes controlling this trait are still not clear in pepper. In this study, complete mitochondrial genomes were sequenced and assembled for the CMS line 138A and its maintainer line 138B. The genome size of 138A is 504,210 bp, which is 8618 bp shorter than that of 138B. Meanwhile, more than 214 and 215 open reading frames longer than 100 amino acids (aas) were identified in 138A and 138B, respectively. Mitochondrial genome structure of 138A was quite different from that of 138B, indicating the existence of recombination and rearrangement events. Based on the mitochondrial genome sequence and structure variations, mitochondrion of 138A and FS4401, a Korean origin CMS line, may have inherited from a common female ancestor, but their CMS traits did originate separately. Candidate gene selection was performed according to the published characteristics of the CMS genes, including the presence SNPs and InDels, located in unique regions, their chimeric structure, co-transcription, and transmembrane domain. A total of 35 ORFs were considered as potential candidate genes and 14 of these were selected, with orf300a and 0rf314a as strong candidates. A new marker, orf300a, was developed which did co-segregate with the CMS trait.

Whole-genome sequencing of Brassica oleracea var. capitata reveals new diversity of the mitogenome

Plant mitochondrial genomes (mtDNAs) vary in sequence structure. We assembled the Brassica oleracea var. capitata mtDNA using a mean coverage depth of 25X whole genome sequencing (WGS) and confirmed the presence of eight contigs/fragments by BLASTZ using the previously reported KJ820683 and AP012988 mtDNA as reference. Assembly of the mtDNA sequence reads resulted in a circular structure of 219,975 bp. Our assembled mtDNA, NCBI acc. no. KU831325, contained 34 protein-coding genes, 3 rRNA genes, and 19 tRNA genes with similarity to the KJ820683 and AP012988 reference mtDNA. No large repeats were found in the KU831325 assembly. However, KU831325 showed differences in the arrangement of bases at different regions compared to the previously reported mtDNAs. In the reference mtDNAs KJ820683 and AP012988, contig/fragment number 4 is partitioned into two contigs/fragments, 4a and 4b. However, contig/fragment number 4 was a single contig/fragment with 29,661 bp in KU831325. PCR and qRT-PCR using flanking markers from separate parts of contig/fragment number 4 confirmed it to be a single contig/fragment. In addition, genome re-alignment of the plastid genome and mtDNAs supported the presence of heteroplasmy and reverse arrangement of the heteroplasmic blocks within the other mtDNAs compared to KU831325 that might be one of the causal factors for its diversity. Our results thus confirm the existence of different mtDNAs in diverse B. oleracea subspecies.

Plant Mitochondrial Genomes: Dynamics and Mechanisms of Mutation

The large mitochondrial genomes of angiosperms are unusually dynamic because of recombination activities involving repeated sequences. These activities generate subgenomic forms and extensive genomic variation even within the same species. Such changes in genome structure are responsible for the rapid evolution of plant mitochondrial DNA and for the variants associated with cytoplasmic male sterility and abnormal growth phenotypes. Nuclear genes modulate these processes, and over the past decade, several of these genes have been identified. They are involved mainly in pathways of DNA repair by homologous recombination and mismatch repair, which appear to be essential for the faithful replication of the mitogenome. Mutations leading to the loss of any of these activities release error-prone repair pathways, resulting in increased ectopic recombination, genome instability, and heteroplasmy. We review the present state of knowledge of the genes and pathways underlying mitochondrial genome stability.

Comparative mitochondrial genome analysis reveals the evolutionary rearrangement mechanism in Brassica

The genus Brassica has many species that are important for oil, vegetable and other food products. Three mitochondrial genome types (mitotype) originated from its common ancestor. In this paper, a B. nigra mitochondrial main circle genome with 232,407 bp was generated through de novo assembly. Synteny analysis showed that the mitochondrial genomes of B. rapa and B. oleracea had a better syntenic relationship than B. nigra. Principal components analysis and development of a phylogenetic tree indicated maternal ancestors of three allotetraploid species in Us triangle of Brassica. Diversified mitotypes were found in allotetraploid B. napus, in which napus-type B. napus was derived from B. oleracea, while polima-type B. napus was inherited from B. rapa. In addition, the mitochondrial genome of napus-type B. napus was closer to botrytis-type than capitata-type B. oleracea. The sub-stoichiometric shifting of several mitochondrial genes suggested that mitochondrial genome rearrangement underwent evolutionary selection during domestication and/or plant breeding. Our findings clarify the role of diploid species in the maternal origin of allotetraploid species in Brassica and suggest the possibility of breeding selection of the mitochondrial genome.

Completion of the mitochondrial genome sequence of onion (Allium cepa L.) containing the CMS-S male-sterile cytoplasm and identification of an independent event of the ccmF N gene split

Cytoplasmic male-sterility (CMS) conferred by the CMS-S cytoplasm has been most commonly used for onion (Allium cepa L.) F₁ hybrid seed production. We first report the complete mitochondrial genome sequence containing CMS-S cytoplasm in this study. Initially, seven contigs were de novo assembled from 150-bp paired-end raw reads produced from the total genomic DNA using the Illumina NextSeq500 platform. These contigs were connected into a single circular genome consisting of 316,363 bp (GenBank accession: KU318712) by PCR amplification. Although all 24 core protein-coding genes were present, no ribosomal protein-coding genes, except rps12, were identified in the onion mitochondrial genome. Unusual trans-splicing of the cox2 gene was verified, and the cox1 gene was identified as part of the chimeric orf725 gene, which is a candidate gene responsible for inducing CMS. In addition to orf725, two small chimeric genes were identified, but no transcripts were detected for these two open reading frames. Thirteen chloroplast-derived sequences, with sizes of 126-13,986 bp, were identified in the intergenic regions. Almost 10 % of the onion mitochondrial genome was composed of repeat sequences. The vast majority of repeats were short repeats of <100 base pairs. Interestingly, the gene encoding ccmF_N was split into two genes. The ccmF _N gene split is first identified outside the Brassicaceae family. The breakpoint in the onion ccmF _N gene was different from that of other Brassicaceae species. This split of the ccmF _N gene was also present in 30 other Allium species. The complete onion mitochondrial genome sequence reported in this study would be fundamental information for elucidation of onion CMS evolution.

MutS HOMOLOG1 silencing mediates ORF220 substoichiometric shifting and causes male sterility in Brassica juncea

Cytoplasmic male sterility (CMS) has consistently been associated with the expression of mitochondrial open reading frames (ORFs) that arise from genomic rearrangements. Spontaneous fertility reversion in CMS has been observed in several cases, but a clear understanding of fertility reversion controlled by nuclear genetic influences has been lacking. Here, we identified spontaneous fertile revertant lines for Brassica juncea CMS cytoplasm in which the mitochondrial genome has undergone substoichiometric shifting (SSS) to suppress ORF220 copy number. We placed ORF220, with or without a mitochondrial targeting presequence, under the control of the CaMV35S and AP3 promoters in Arabidopsis to confirm that ORF220 causes male sterility when mitochondrially localized. We found that copy number of the ORF220 gene was altered under conditions that suppress MSH1, a nuclear gene that controls illegitimate recombination in plant mitochondria. MSH1-RNAi lines with increased ORF220 copy number were male sterile compared with wild type. We found that a wide range of genes involved in anther development were up- and down-regulated in revertant and MSH1-RNAi lines, respectively. The system that we have developed offers valuable future insight into the interplay of MSH1 and SSS in CMS induction and fertility reversion as a mediator of nuclear-mitochondrial crosstalk.

Identification of candidate genes associated with fertility restoration of cytoplasmic male-sterility in onion (Allium cepa L.) using a combination of bulked segregant analysis and RNA-seq

A combination of BSA and RNA-seq was performed to identify candidates for the restorer-of-fertility gene in onion. The AcPMS1 involved in DNA mismatch repair was identified as the best candidate. To identify candidate genes of the restorer-of-fertility gene (Ms) responsible for fertility restoration of onion cytoplasmic male-sterility, a combined approach of bulked segregant analysis and RNA-seq was employed. From 32,674 de novo assembled contigs, 430 perfectly homozygous SNPs between male-fertile (MF) and male-sterile (MS) bulks were identified in 141 contigs. After verifying the homozygosity of the SNPs by PCR amplification and sequencing, the SNPs on 139 of the contigs were genotypes for the two recombinants which contained crossover events between the Ms locus and two tightly linked molecular markers. As a result, 30 contigs showing perfect linkage with the Ms locus in the large-sized segregating population were identified. Among them, 14 showed perfect linkage disequilibrium (LD) with the Ms locus, as determined by genotyping 251 domestic breeding lines. Furthermore, molecular markers tagging the 14 contigs also showed almost perfect LD with each other in 124 exotic accessions introduced from 21 countries, except for one accession which contained a crossover event by which the 14 markers were divided into two groups. After sequencing of the full-length cDNA of the 14 contigs showing perfect LD, the deduced amino acids sequences of the MF and MS alleles were compared. Four genes were shown to harbor putative critical amino acid changes in the known domains. Among them, the gene encoding PMS1, involved in the DNA mismatch repair pathway, was assumed to be the best candidate gene responsible for fertility restoration of male-sterility in onion.

Miniaturized mitogenome of the parasitic plant Viscum scurruloideum is extremely divergent and dynamic and has lost all nad genes

Despite the enormous diversity among parasitic angiosperms in form and structure, life-history strategies, and plastid genomes, little is known about the diversity of their mitogenomes. We report the sequence of the wonderfully bizarre mitogenome of the hemiparasitic aerial mistletoe Viscum scurruloideum. This genome is only 66 kb in size, making it the smallest known angiosperm mitogenome by a factor of more than three and the smallest land plant mitogenome. Accompanying this size reduction is exceptional reduction of gene content. Much of this reduction arises from the unexpected loss of respiratory complex I (NADH dehydrogenase), universally present in all 300+ other angiosperms examined, where it is encoded by nine mitochondrial and many nuclear nad genes. Loss of complex I in a multicellular organism is unprecedented. We explore the potential relationship between this loss in Viscum and its parasitic lifestyle. Despite its small size, the Viscum mitogenome is unusually rich in recombinationally active repeats, possessing unparalleled levels of predicted sublimons resulting from recombination across short repeats. Many mitochondrial gene products exhibit extraordinary levels of divergence in Viscum, indicative of highly relaxed if not positive selection. In addition, all Viscum mitochondrial protein genes have experienced a dramatic acceleration in synonymous substitution rates, consistent with the hypothesis of genomic streamlining in response to a high mutation rate but completely opposite to the pattern seen for the high-rate but enormous mitogenomes of Silene. In sum, the Viscum mitogenome possesses a unique constellation of extremely unusual features, a subset of which may be related to its parasitic lifestyle.

A high-resolution linkage map of the Rfd1, a restorer-of-fertility locus for cytoplasmic male sterility in radish (Raphanus sativus L.) produced by a combination of bulked segregant analysis and RNA-Seq

We utilized a combination of BSA and RNA-Seq to identify SNPs linked to the Rfd1 locus, a restorer-of-fertility gene in radish. A high-density linkage map was constructed using this approach. Male fertility of cytoplasmic male sterility conditioned by the Dongbu cytoplasmic and genic male-sterility cytoplasm can be restored by a restorer-of-fertility locus, Rfd1, in radish. To construct a high-density linkage map and to identify a candidate gene for the Rfd1 locus, bulked segregant analysis and RNA-seq approaches were combined. A total of 26 and 28 million reads produced from male-fertile and male-sterile bulked RNA were mapped to the radish reference unigenes. After stringent screening of SNPs, 327 reliable SNPs of 109 unigenes were selected. Arabidopsis homologs for 101 of the 109 genes were clustered around the 4,000 kb region of Arabidopsis chromosome 3, which was syntenic to the Rfd1 flanking region. Since the reference unigene set was incomplete, the contigs were de novo assembled to identify 134 contigs harboring SNPs. Most of SNP-containing contigs were also clustered on the same syntenic region in Arabidopsis chromosome. A total of 21 molecular markers positioned within a 2.1 cM interval including the Rfd1 locus were developed, based on the selected unigenes and contigs. A segregating population consisting of 10,459 individuals was analyzed to identify recombinants containing crossovers within this interval. A total of 284 identified recombinants were then used to construct a high-density map, which delimited the Rfd1 locus into an 83-kb syntenic interval of Arabidopsis chromosome 3. Since no candidate gene, such as a pentatricopeptide repeat (PPR)-coding gene, was found in this interval, 231 unigenes and 491 contigs containing putative PPR motifs were analyzed further, but no PPR gene in linkage disequilibrium with the Rfd1 locus could be found.

Extensive structural variations between mitochondrial genomes of CMS and normal peppers (Capsicum annuum L.) revealed by complete nucleotide sequencing

BACKGROUND

Cytoplasmic male sterility (CMS) is an inability to produce functional pollen that is caused by mutation of the mitochondrial genome. Comparative analyses of mitochondrial genomes of lines with and without CMS in several species have revealed structural differences between genomes, including extensive rearrangements caused by recombination. However, the mitochondrial genome structure and the DNA rearrangements that may be related to CMS have not been characterized in Capsicum spp.

RESULTS

We obtained the complete mitochondrial genome sequences of the pepper CMS line FS4401 (507,452 bp) and the fertile line Jeju (511,530 bp). Comparative analysis between mitochondrial genomes of peppers and tobacco that are included in Solanaceae revealed extensive DNA rearrangements and poor conservation in non-coding DNA. In comparison between pepper lines, FS4401 and Jeju mitochondrial DNAs contained the same complement of protein coding genes except for one additional copy of an atp6 gene (ψatp6-2) in FS4401. In terms of genome structure, we found eighteen syntenic blocks in the two mitochondrial genomes, which have been rearranged in each genome. By contrast, sequences between syntenic blocks, which were specific to each line, accounted for 30,380 and 17,847 bp in FS4401 and Jeju, respectively. The previously-reported CMS candidate genes, orf507 and ψatp6-2, were located on the edges of the largest sequence segments that were specific to FS4401. In this region, large number of small sequence segments which were absent or found on different locations in Jeju mitochondrial genome were combined together. The incorporation of repeats and overlapping of connected sequence segments by a few nucleotides implied that extensive rearrangements by homologous recombination might be involved in evolution of this region. Further analysis using mtDNA pairs from other plant species revealed common features of DNA regions around CMS-associated genes.

CONCLUSIONS

Although large portion of sequence context was shared by mitochondrial genomes of CMS and male-fertile pepper lines, extensive genome rearrangements were detected. CMS candidate genes located on the edges of highly-rearranged CMS-specific DNA regions and near to repeat sequences. These characteristics were detected among CMS-associated genes in other species, implying a common mechanism might be involved in the evolution of CMS-associated genes.

The complete mitochondrial genome sequence of Brassica oleracea and analysis of coexisting mitotypes

The complete mitochondrial genome sequences of Brassica species have provided insight into inter- and intraspecific variation of plant mitochondrial genomes. However, the size of mitochondrial genome sequenced for Brassica oleracea hitherto does not match to its physical mapping data. This fact led us to investigate B. oleracea mitochondrial genome in detail. Here we report novel B. oleracea mitochondrial genome, derived from var. capitata, a cabbage cultivar ''Fujiwase''. The genome was assembled into a 219,952-bp circular sequence that is comparable to the mitochondrial genomes of other Brassica species (ca. 220-232 kb). This genome contained 34 protein-coding genes, 3 rRNA genes and 17 tRNA genes. Due to absence of a large repeat (140 kb), the mitochondrial genome of ''Fujiwase'' is clearly smaller than the previously reported mitochondrial genome of B. oleracea accession ''08C717'' (360 kb). In both mitotypes, all genes were identical, except cox2-2, which was present only in the Fujiwase type. At least two rearrangement events via large and small repeat sequences have contributed to the structural differences between the two mitotypes. PCR-based marker analysis revealed that the Fujiwase type is predominant, whereas the 08C717 type coexists at low frequency in all B. oleracea cultivars examined. Intraspecific variations in the mitochondrial genome in B. oleracea may occur because of heteroplasmy, coexistence of different mitotypes within an individual, and substoichiometric shifting. Our data indicate that the Fujiwase-type genome should be used as the representative genome of B. oleracea.

Comparative analysis of mitochondrial genomes between the hau cytoplasmic male sterility (CMS) line and its iso-nuclear maintainer line in Brassica juncea to reveal the origin of the CMS-associated gene orf288

Background

Cytoplasmic male sterility (CMS) is not only important for exploiting heterosis in crop plants, but also as a model for investigating nuclear-cytoplasmic interaction. CMS may be caused by mutations, rearrangement or recombination in the mitochondrial genome. Understanding the mitochondrial genome is often the first and key step in unraveling the molecular and genetic basis of CMS in plants. Comparative analysis of the mitochondrial genome of the hau CMS line and its maintainer line in B. juneca (Brassica juncea) may help show the origin of the CMS-associated gene orf288.

Results

Through next-generation sequencing, the B. juncea hau CMS mitochondrial genome was assembled into a single, circular-mapping molecule that is 247,903 bp in size and 45.08% in GC content. In addition to the CMS associated gene orf288, the genome contains 35 protein-encoding genes, 3 rRNAs, 25 tRNA genes and 29 ORFs of unknown function. The mitochondrial genome sizes of the maintainer line and another normal type line “J163-4” are both 219,863 bp and with GC content at 45.23%. The maintainer line has 36 genes with protein products, 3 rRNAs, 22 tRNA genes and 31 unidentified ORFs. Comparative analysis the mitochondrial genomes of the hau CMS line and its maintainer line allowed us to develop specific markers to separate the two lines at the seedling stage. We also confirmed that different mitotypes coexist substoichiometrically in hau CMS lines and its maintainer lines in B. juncea. The number of repeats larger than 100 bp in the hau CMS line (16 repeats) are nearly twice of those found in the maintainer line (9 repeats). Phylogenetic analysis of the CMS-associated gene orf288 and four other homologous sequences in Brassicaceae show that orf288 was clearly different from orf263 in Brassica tournefortii despite of strong similarity.

Conclusion

The hau CMS mitochondrial genome was highly rearranged when compared with its iso-nuclear maintainer line mitochondrial genome. This study may be useful for studying the mechanism of natural CMS in B. juncea, performing comparative analysis on sequenced mitochondrial genomes in Brassicas, and uncovering the origin of the hau CMS mitotype and structural and evolutionary differences between different mitotypes.

Electronic supplementary material

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

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.

Biparental inheritance of organelles in Pelargonium: evidence for intergenomic recombination of mitochondrial DNA

While uniparental transmission of mtDNA is widespread and dominating in eukaryotes leaving mutation as the major source of genotypic diversity, recently, biparental inheritance of mitochondrial genes has been demonstrated in reciprocal crosses of Pelargonium zonale and P. inquinans. The thereby arising heteroplasmy carries the potential for recombination between mtDNAs of different descent, i.e. between the parental mitochondrial genomes. We have analyzed these Pelargonium hybrids for mitochondrial intergenomic recombination events by examining differences in DNA blot hybridization patterns of the mitochondrial genes atp1 and cob. Further investigation of these genes and their flanking regions using nucleotide sequence polymorphisms and PCR revealed DNA segments in the progeny, which contained both P. zonale and P. inquinans sequences suggesting an intergenomic recombination in hybrids of Pelargonium. This turns Pelargonium into an interesting subject for studies of recombination and evolutionary dynamics of mitochondrial genomes.

Cytoplasmic male sterility in plants: molecular evidence and the nucleocytoplasmic conflict

A much-debated issue in plant evolutionary biology concerns the maintenance of a high frequency of male sterility in natural populations. For the past decade, a theoretical framework has been provided by the concept of nucleocytoplasmic conflict. Recent molecular studies on cytoplasmic male sterility indicate that novel chimeric genes, resulting from duplications and rearrangements of mitochondrial DNA sequences, are involved In its control. Thus, male sterility, which is phenotypically the loss of the male function, is encoded by a new mitochondrial function at the molecular level. Molecular data are in agreement with theoretical models that consider cytoplasmic male sterility as a stage in the coevolution between nucleus and mitochondria, and not simply as a deleterious mitochondrial mutation.

Construction of a high-resolution linkage map of Rfd1, a restorer-of-fertility locus for cytoplasmic male sterility conferred by DCGMS cytoplasm in radish (Raphanus sativus L.) using synteny between radish and Arabidopsis genomes

Cytoplasmic male sterility caused by Dongbu cytoplasmic and genic male-sterility (DCGMS) cytoplasm and its nuclear restorer-of-fertility locus (Rfd1) with a linked molecular marker (A137) have been reported in radish (Raphanus sativus L.). To construct a linkage map of the Rfd1 locus, linked amplified fragment length polymorphism (AFLP) markers were screened using bulked segregant analysis. A 220-bp linked AFLP fragment sequence from radish showed homology with an Arabidopsis coding sequence. Using this Arabidopsis gene sequence, a simple PCR marker (A220) was developed. The A137 and A220 markers flanked the Rfd1 locus. Two homologous Arabidopsis genes with both marker sequences were positioned on Arabidopsis chromosome-3 with an interval of 2.4 Mb. To integrate the Rfd1 locus into a previously reported expressed sequence tag (EST)-simple sequence repeat (SSR) linkage map, the radish EST sequences located in three syntenic blocks within the 2.4-Mb interval were used to develop single nucleotide polymorphism (SNP) markers for tagging each block. The SNP marker in linkage group-2 co-segregated with male fertility in an F(2) population. Using radish ESTs positioned in linkage group-2, five intron length polymorphism (ILP) markers and one cleaved amplified polymorphic sequence (CAPS) marker were developed and used to construct a linkage map of the Rfd1 locus. Two closely linked markers delimited the Rfd1 locus within a 985-kb interval of Arabidopsis chromosome-3. Synteny between the radish and Arabidopsis genomes in the 985-kb interval were used to develop three ILP and three CAPS markers. Two ILP markers further delimited the Rfd1 locus to a 220-kb interval of Arabidopsis chromosome-3.

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.

RecA-dependent DNA repair results in increased heteroplasmy of the Arabidopsis mitochondrial genome

Plant mitochondria have very active DNA recombination activities that are responsible for its plastic structures and that should be involved in the repair of double-strand breaks in the mitochondrial genome. Little is still known on plant mitochondrial DNA repair, but repair by recombination is believed to be a major determinant in the rapid evolution of plant mitochondrial genomes. In flowering plants, mitochondria possess at least two eubacteria-type RecA proteins that should be core components of the mitochondrial repair mechanisms. We have performed functional analyses of the two Arabidopsis (Arabidopsis thaliana) mitochondrial RecAs (RECA2 and RECA3) to assess their potential roles in recombination-dependent repair. Heterologous expression in Escherichia coli revealed that RECA2 and RECA3 have overlapping as well as specific activities that allow them to partially complement bacterial repair pathways. RECA2 and RECA3 have similar patterns of expression, and mutants of either display the same molecular phenotypes of increased recombination between intermediate-size repeats, thus suggesting that they act in the same recombination pathways. However, RECA2 is essential past the seedling stage and should have additional important functions. Treatment of plants with several DNA-damaging drugs further showed that RECA3 is required for different recombination-dependent repair pathways that significantly contribute to plant fitness under stress. Replication repair of double-strand breaks results in the accumulation of crossovers that increase the heteroplasmic state of the mitochondrial DNA. It was shown that these are transmitted to the plant progeny, enhancing the potential for mitochondrial genome evolution.

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.

A reevaluation of rice mitochondrial evolution based on the complete sequence of male-fertile and male-sterile mitochondrial genomes

Plant mitochondrial genomes have features that distinguish them radically from their animal counterparts: a high rate of rearrangement, of uptake and loss of DNA sequences, and an extremely low point mutation rate. Perhaps the most unique structural feature of plant mitochondrial DNAs is the presence of large repeated sequences involved in intramolecular and intermolecular recombination. In addition, rare recombination events can occur across shorter repeats, creating rearrangements that result in aberrant phenotypes, including pollen abortion, which is known as cytoplasmic male sterility (CMS). Using next-generation sequencing, we pyrosequenced two rice (Oryza sativa) mitochondrial genomes that belong to the indica subspecies. One genome is normal, while the other carries the wild abortive-CMS. We find that numerous rearrangements in the rice mitochondrial genome occur even between close cytotypes during rice evolution. Unlike maize (Zea mays), a closely related species also belonging to the grass family, integration of plastid sequences did not play a role in the sequence divergence between rice cytotypes. This study also uncovered an excellent candidate for the wild abortive-CMS-encoding gene; like most of the CMS-associated open reading frames that are known in other species, this candidate was created via a rearrangement, is chimeric in structure, possesses predicted transmembrane domains, and coopted the promoter of a genuine mitochondrial gene. Our data give new insights into rice mitochondrial evolution, correcting previous reports.

The complete chloroplast and mitochondrial genome sequences of Boea hygrometrica: insights into the evolution of plant organellar genomes

The complete nucleotide sequences of the chloroplast (cp) and mitochondrial (mt) genomes of resurrection plant Boea hygrometrica (Bh, Gesneriaceae) have been determined with the lengths of 153,493 bp and 510,519 bp, respectively. The smaller chloroplast genome contains more genes (147) with a 72% coding sequence, and the larger mitochondrial genome have less genes (65) with a coding faction of 12%. Similar to other seed plants, the Bh cp genome has a typical quadripartite organization with a conserved gene in each region. The Bh mt genome has three recombinant sequence repeats of 222 bp, 843 bp, and 1474 bp in length, which divide the genome into a single master circle (MC) and four isomeric molecules. Compared to other angiosperms, one remarkable feature of the Bh mt genome is the frequent transfer of genetic material from the cp genome during recent Bh evolution. We also analyzed organellar genome evolution in general regarding genome features as well as compositional dynamics of sequence and gene structure/organization, providing clues for the understanding of the evolution of organellar genomes in plants. The cp-derived sequences including tRNAs found in angiosperm mt genomes support the conclusion that frequent gene transfer events may have begun early in the land plant lineage.

The mitochondrial genomes of the early land plants Treubia lacunosa and Anomodon rugelii: dynamic and conservative evolution

Early land plant mitochondrial genomes captured important changes of mitochondrial genome evolution when plants colonized land. The chondromes of seed plants show several derived characteristics, e.g., large genome size variation, rapid intra-genomic rearrangement, abundant introns, and highly variable levels of RNA editing. On the other hand, the chondromes of charophytic algae are still largely ancestral in these aspects, resembling those of early eukaryotes. When the transition happened has been a long-standing question in studies of mitochondrial genome evolution. Here we report complete mitochondrial genome sequences from an early-diverging liverwort, Treubia lacunosa, and a late-evolving moss, Anomodon rugelii. The two genomes, 151,983 and 104,239 base pairs in size respectively, contain standard sets of protein coding genes for respiration and protein synthesis, as well as nearly full sets of rRNA and tRNA genes found in the chondromes of the liverworts Marchantia polymorpha and Pleurozia purpurea and the moss Physcomitrella patens. The gene orders of these two chondromes are identical to those of the other liverworts and moss. Their intron contents, with all cis-spliced group I or group II introns, are also similar to those in the previously sequenced liverwort and moss chondromes. These five chondromes plus the two from the hornworts Phaeoceros laevis and Megaceros aenigmaticus for the first time allowed comprehensive comparative analyses of structure and organization of mitochondrial genomes both within and across the three major lineages of bryophytes. These analyses led to the conclusion that the mitochondrial genome experienced dynamic evolution in genome size, gene content, intron acquisition, gene order, and RNA editing during the origins of land plants and their major clades. However, evolution of this organellar genome has remained rather conservative since the origin and initial radiation of early land plants, except within vascular plants.

Unique changes in mitochondrial genomes associated with reversions of S-type cytoplasmic male sterility in maizemar

Cytoplasmic male sterility (CMS) in plants is usually associated with the expression of specific chimeric regions within rearranged mitochondrial genomes. Maize CMS-S plants express high amounts of a 1.6-kb mitochondrial RNA during microspore maturation, which is associated with the observed pollen abortion. This transcript carries two chimeric open reading frames, orf355 and orf77, both unique to CMS-S. CMS-S mitochondria also contain free linear DNA plasmids bearing terminal inverted repeats (TIRs). These TIRs recombine with TIR-homologous sequences that precede orf355/orf77 within the main mitochondrial genome to produce linear ends. Transcription of the 1.6-kb RNA is initiated from a promoter within the TIRs only when they are at linear ends. Reversions of CMS-S to fertility occur in certain nuclear backgrounds and are usually associated with loss of the S plasmids and/or the sterility-associated region. We describe an unusual set of independently recovered revertants from a single maternal lineage that retain both the S plasmids and an intact orf355/orf77 region but which do not produce the 1.6-kb RNA. A 7.3-kb inversion resulting from illegitmate recombination between 14-bp microrepeats has separated the genomic TIR sequences from the CMS-associated region. Although RNAs containing orf355/orf77 can still be detected in the revertants, they are not highly expressed during pollen development and they are no longer initiated from the TIR promoter at a protein-stabilized linear end. They appear instead to be co-transcribed with cytochrome oxidase subunit 2. The 7.3-kb inversion was not detected in CMS-S or in other fertile revertants. Therefore, this inversion appears to be a de novo mutation that has continued to sort out within a single maternal lineage, giving rise to fertile progeny in successive generations.

Origin of the CMS gene locus in rapeseed cybrid mitochondria: active and inactive recombination produces the complex CMS gene region in the mitochondrial genomes of Brassicaceae

CMS (cytoplasmic male sterile) rapeseed is produced by asymmetrical somatic cell fusion between the Brassica napus cv. Westar and the Raphanus sativus Kosena CMS line (Kosena radish). The CMS rapeseed contains a CMS gene, orf125, which is derived from Kosena radish. Our sequence analyses revealed that the orf125 region in CMS rapeseed originated from recombination between the orf125/orfB region and the nad1C/ccmFN1 region by way of a 63 bp repeat. A precise sequence comparison among the related sequences in CMS rapeseed, Kosena radish and normal rapeseed showed that the orf125 region in CMS rapeseed consisted of the Kosena orf125/orfB region and the rapeseed nad1C/ccmFN1 region, even though Kosena radish had both the orf125/orfB region and the nad1C/ccmFN1 region in its mitochondrial genome. We also identified three tandem repeat sequences in the regions surrounding orf125, including a 63 bp repeat, which were involved in several recombination events. Interestingly, differences in the recombination activity for each repeat sequence were observed, even though these sequences were located adjacent to each other in the mitochondrial genome. We report results indicating that recombination events within the mitochondrial genomes are regulated at the level of specific repeat sequences depending on the cellular environment.

Substoichiometrically different mitotypes coexist in mitochondrial genomes of Brassica napus L

Cytoplasmic male sterility (CMS) has been identified in numerous plant species. Brassica napus CMS plants, such as Polima (pol), MI, and Shaan 2A, have been identified independently by different researchers with different materials in conventional breeding processes. How this kind of CMS emerges is unclear. Here, we report the mitochondrial genome sequence of the prevalent mitotype in the most widely used pol-CMS line, which has a length of 223,412 bp and encodes 34 proteins, 3 ribosomal RNAs, and 18 tRNAs, including two near identical copies of trnH. Of these 55 genes, 48 were found to be identical to their equivalents in the “nap” cytoplasm. The nap mitotype carries only one copy of trnH, and the sequences of five of the six remaining genes are highly similar to their equivalents in the pol mitotype. Forty-four open reading frames (ORFs) with unknown function were detected, including two unique to the pol mitotype (orf122 and orf132). At least five rearrangement events are required to account for the structural differences between the pol and nap sequences. The CMS-related orf224 neighboring region (∼5 kb) rearranged twice. PCR profiling based on mitotype-specific primer pairs showed that both mitotypes are present in B. napus cultivars. Quantitative PCR showed that the pol cytoplasm consists mainly of the pol mitotype, and the nap mitotype is the main genome of nap cytoplasm. Large variation in the copy number ratio of mitotypes was found, even among cultivars sharing the same cytoplasm. The coexistence of mitochondrial mitotypes and substoichiometric shifting can explain the emergence of CMS in B. napus.

Ribosomal RNA genes in Euglena gracilis mitochondrial DNA: fragmented genes in a seemingly fragmented genome

Because relatively little information is available about mtDNA in the euglenid protozoa, distant relatives of the kinetoplastid protozoa, we investigated mitochondrial genome structure and expression in Euglena gracilis. We found that isolated E. gracilis mtDNA comprises a heterodisperse collection of short molecules (modal size approximately 4 kbp) and that the mitochondrial large subunit (LSU) and small subunit (SSU) rRNAs are each split into two pieces. For the two halves of the SSU rRNA, we identified separate, non-contiguous coding modules that are flanked by a complex array of (primarily direct) A + T-rich repeats. The potential secondary structure of the bipartite SSU rRNA displays the expected conserved elements implicated in ribosome function. Label from [α-(32)P]GTP was incorporated in the presence of guanylyltransferase into each of the separate SSU and LSU rRNA fragments, confirming that these RNAs are primary transcripts, separately expressed from non-contiguous rRNA modules. In addition to authentic genes for SSU rRNA, we discovered numerous short fragments of protein-coding and rRNA genes dispersed throughout the E. gracilis mitochondrial genome. We propose that antisense transcripts of gene fragments of this type could have been the evolutionary precursors of the guide RNAs that mediate U insertion/deletion editing in the kinetoplastid relatives of the euglenids. To the extent that E. gracilis mtDNA is a representative euglenid mitochondrial genome, it differs radically in structure and organization from that of its kinetoplastid relatives, instead more closely resembling the mitochondrial genome of dinoflagellates in many of its features, an apparent evolutionary convergence.

Fewer genes than organelles: extremely low and variable gene copy numbers in mitochondria of somatic plant cells

Plant mitochondrial genomes are split into sub-genomes, i.e. genes are distributed across various sub-genomic molecules. To investigate whether copy numbers vary between individual mitochondrial genes, we used quantitative real-time PCR in combination with flow cytometric determination of nuclear DNA quantities to determine absolute per-cell-copy numbers of four mitochondrial genes in various Arabidopsis organs and the leaves of tobacco (Nicotiana tabacum) and barley (Hordeum vulgare). The copy numbers of the investigated mitochondrial genes (atp1, rps4, nad6 and cox1) not only differed from each other, but also varied between organs and changed during the development of cotyledons and leaves in Arabidopsis. We found no correlation between altered gene copy numbers, transcript levels and O(2) consumption. However, per cell, both the number of mitochondria and the number of gene copies increased with growing cell size. Gene copy numbers varied from approximately 40 (cox1 in young leaves) to approximately 280 (atp1 in mature leaves), and the mean number of mitochondria was approximately 300 in young leaves and 450 in mature leaves. Thus, cells are polyploid with respect to their mitochondrial genomes, but individual mitochondria may contain only part of the genome or even no DNA at all. Our data supports structural models of the mitochondrial genome in non-dividing cells of angiosperms that predict localization of the genes on sub-genomic molecules rather than master chromosomes. The data indicate control of the number of individual genes according to the genotype and developmental program(s) via amplification and/or degradation of sub-genomic molecules.

Cost of Having the Largest Mitochondrial Genome: Evolutionary Mechanism of Plant Mitochondrial Genome

The angiosperm mitochondrial genome is the largest and least gene-dense among the eukaryotes, because its intergenic regions are expanded. There seems to be no functional constraint on the size of the intergenic regions; angiosperms maintain the large mitochondrial genome size by a currently unknown mechanism. After a brief description of the angiosperm mitochondrial genome, this review focuses on our current knowledge of the mechanisms that control the maintenance and alteration of the genome. In both processes, the control of homologous recombination is crucial in terms of site and frequency. The copy numbers of various types of mitochondrial DNA molecules may also be controlled, especially during transmission of the mitochondrial genome from one generation to the next. An important characteristic of angiosperm mitochondria is that they contain polypeptides that are translated from open reading frames created as byproducts of genome alteration and that are generally nonfunctional. Such polypeptides have potential to evolve into functional ones responsible for mitochondrially encoded traits such as cytoplasmic male sterility or may be remnants of the former functional polypeptides.

A scenario of mitochondrial genome evolution in maize based on rearrangement events

BACKGROUND

Despite their monophyletic origin, animal and plant mitochondrial genomes have been described as exhibiting different modes of evolution. Indeed, plant mitochondrial genomes feature a larger size, a lower mutation rate and more rearrangements than their animal counterparts. Gene order variation in animal mitochondrial genomes is often described as being due to translocation and inversion events, but tandem duplication followed by loss has also been proposed as an alternative process. In plant mitochondrial genomes, at the species level, gene shuffling and duplicate occurrence are such that no clear phylogeny has ever been identified, when considering genome structure variation.

RESULTS

In this study we analyzed the whole sequences of eight mitochondrial genomes from maize and teosintes in order to comprehend the events that led to their structural features, i.e. the order of genes, tRNAs, rRNAs, ORFs, pseudogenes and non-coding sequences shared by all mitogenomes and duplicate occurrences. We suggest a tandem duplication model similar to the one described in animals, except that some duplicates can remain. This model enabled us to develop a manual method to deal with duplicates, a recurrent problem in rearrangement analyses. The phylogenetic tree exclusively based on rearrangement and duplication events is congruent with the tree based on sequence polymorphism, validating our evolution model.

CONCLUSIONS

This study suggests more similarity than usually reported between plant and animal mitochondrial genomes in their mode of evolution. Further work will consist of developing new tools in order to automatically look for signatures of tandem duplication events in other plant mitogenomes and evaluate the occurrence of this process on a larger scale.

Recombination and the maintenance of plant organelle genome stability

Like their nuclear counterpart, the plastid and mitochondrial genomes of plants have to be faithfully replicated and repaired to ensure the normal functioning of the plant. Inability to maintain organelle genome stability results in plastid and/or mitochondrial defects, which can lead to potentially detrimental phenotypes. Fortunately, plant organelles have developed multiple strategies to maintain the integrity of their genetic material. Of particular importance among these processes is the extensive use of DNA recombination. In fact, recombination has been implicated in both the replication and the repair of organelle genomes. Revealingly, deregulation of recombination in organelles results in genomic instability, often accompanied by adverse consequences for plant fitness. The recent identification of four families of proteins that prevent aberrant recombination of organelle DNA sheds much needed mechanistic light on this important process. What comes out of these investigations is a partial portrait of the recombination surveillance machinery in which plants have co-opted some proteins of prokaryotic origin but have also evolved whole new factors to keep their organelle genomes intact. These new features presumably optimized the protection of plastid and mitochondrial genomes against the particular genotoxic stresses they face.

Discovery of global genomic re-organization based on comparison of two newly sequenced rice mitochondrial genomes with cytoplasmic male sterility-related genes

Background

Plant mitochondrial genomes are known for their complexity, and there is abundant evidence demonstrating that this organelle is important for plant sexual reproduction. Cytoplasmic male sterility (CMS) is a phenomenon caused by incompatibility between the nucleus and mitochondria that has been discovered in various plant species. As the exact sequence of steps leading to CMS has not yet been revealed, efforts should be made to elucidate the factors underlying the mechanism of this important trait for crop breeding.

Results

Two CMS mitochondrial genomes, LD-CMS, derived from Oryza sativa L. ssp. indica (434,735 bp), and CW-CMS, derived from Oryza rufipogon Griff. (559,045 bp), were newly sequenced in this study. Compared to the previously sequenced Nipponbare (Oryza sativa L. ssp. japonica) mitochondrial genome, the presence of 54 out of 56 protein-encoding genes (including pseudo-genes), 22 tRNA genes (including pseudo-tRNAs), and three rRNA genes was conserved. Two other genes were not present in the CW-CMS mitochondrial genome, and one of them was present as part of the newly identified chimeric ORF, CW-orf307. At least 12 genomic recombination events were predicted between the LD-CMS mitochondrial genome and Nipponbare, and 15 between the CW-CMS genome and Nipponbare, and novel genetic structures were formed by these genomic rearrangements in the two CMS lines. At least one of the genomic rearrangements was completely unique to each CMS line and not present in 69 rice cultivars or 9 accessions of O. rufipogon.

Conclusion

Our results demonstrate novel mitochondrial genomic rearrangements that are unique in CMS cytoplasm, and one of the genes that is unique in the CW mitochondrial genome, CW-orf307, appeared to be the candidate most likely responsible for the CW-CMS event. Genomic rearrangements were dynamic in the CMS lines in comparison with those of rice cultivars, suggesting that 'death' and possible 'birth' processes of the CMS genes occurred during the breeding history of rice.