Unicircular structure of the Brassica hirta mitochondrial genome

Highly active repeat-mediated recombination in the mitogenome of the aquatic grass Hygroryza aristata

BACKGROUND

Floating bamboo (Hygroryza aristata) is an endangered species with a narrow native distribution and is renowned for its unique aesthetic qualities, which holds significant ecological and ornamental value. However, the lack of genetic information research, with only one complete plastome available, significantly hampers conservation efforts and further research for this species.

RESULTS

In this research, we sequenced and assembled the organelle genomes of floating bamboo, including the mitogenome (587,847 bp) and plastome (135,675 bp). The mitogenome can recombine into various configurations, which are mediated by 25 repeat pairs (13 SRs, 6 MRs, 1 LR, and 5 CRs). LR1 and SR5 are particularly notable as they have the ability to combine with other contigs, forming complex repeat units that facilitate further homologous recombination. The rate of homologous recombination varies significantly among species, yet there is still a pronounced positive correlation observed between the length of these repeat pairs and the rate of recombination they mediate. The mitogenome integrates seven intact protein-coding genes from the chloroplast. The codon usage patterns in both organelles are similar, with a noticeable bias towards C and T on the third codon. The gene map of Poales shows the entire loss of rpl6, succinate dehydrogenase subunits (sdh3 and sdh4). Additionally, the BOP clade retained more variable genes compared to the PACMAD clade.

CONCLUSIONS

We provided a high-quality and well-annotated mitogenome for floating bamboo and demonstrated the presence of diverse configurations. Our study has revealed the correlation between repeat length and their corresponding recombination rate despite variations among species. Although the mitogenome can potentially exist in the form of a unicircular in vivo, this occurrence is rare and may not be stable.

The first two whole mitochondrial genomes for the genus Dactylis species: assembly and comparative genomics analysis

BACKGROUND

Orchardgrass (Dactylis glomerata L.), a perennial forage, has the advantages of rich leaves, high yield, and good quality and is one of the most significant forage for grassland animal husbandry and ecological management in southwest China. Mitochondrial (mt) genome is one of the major genetic systems in plants. Studying the mt genome of the genus Dactylis could provide more genetic information in addition to the nuclear genome project of the genus.

RESULTS

In this study, we sequenced and assembled two mitochondrial genomes of Dactylis species of D. glomerata (597, 281 bp) and D. aschersoniana (613, 769 bp), based on a combination of PacBio and Illumina. The gene content in the mitochondrial genome of D. aschersoniana is almost identical to the mitochondrial genome of D. glomerata, which contains 22-23 protein-coding genes (PCGs), 8 ribosomal RNAs (rRNAs) and 30 transfer RNAs (tRNAs), while D. glomerata lacks the gene encoding the Ribosomal protein (rps1) and D. aschersoniana contains one pseudo gene (atp8). Twenty-three introns were found among eight of the 30 protein-coding genes, and introns of three genes (nad 1, nad2, and nad5) were trans-spliced in Dactylis aschersoniana. Further, our mitochondrial genome characteristics investigation of the genus Dactylis included codon usage, sequences repeats, RNA editing and selective pressure. The results showed that a large number of short repetitive sequences existed in the mitochondrial genome of D. aschersoniana, the size variation of two mitochondrial genomes is due largely to the presence of a large number of short repetitive sequences. We also identified 52-53 large fragments that were transferred from the chloroplast genome to the mitochondrial genome, and found that the similarity was more than 70%. ML and BI methods used in phylogenetic analysis revealed that the evolutionary status of the genus Dactylis.

CONCLUSIONS

Thus, this study reveals the significant rearrangements in the mt genomes of Pooideae species. The sequenced Dactylis mt genome can provide more genetic information and improve our evolutionary understanding of the mt genomes of gramineous plants.

Comparative analysis of mitochondrial genomes provides insights into the mechanisms underlying an S-type cytoplasmic male sterility (CMS) system in wheat (Triticum aestivum L.)

Cytoplasmic male sterility (CMS) has been widely used in crop cross breeding. There has been much research on wheat CMS. However, the correlation between S-type CMS and mitochondrial genome remains elusive. Herein, we sequenced the mitochondrial genome of wheat CMS line and compared it with the maintainer line. The results showed that the mitochondrial genome of CMS line encoded 26 tRNAs, 8 rRNAs, and 35 protein-coding genes, and the cob encoding complex III in which the protein coding gene is mutated. This protein is known to affect reactive oxygen (ROS) production. The analysis of ROS metabolism in developing anthers showed that the deficiency of antioxidants and antioxidant enzymes in the sterile system aggravated membrane lipid oxidation, resulting in ROS accumulation, and influencing the anther development. Herein, cob is considered as a candidate causative gene sequence for CMS.

Complete Mitogenomes of Two Aragoa Species and Phylogeny of Plantagineae (Plantaginaceae, Lamiales) Using Mitochondrial Genes and the Nuclear Ribosomal RNA Repeat

Aragoa, comprising 19 high-altitude North Andean species, is one of three genera in the Plantagineae (Plantaginaceae, Lamiales), along with Littorella and Plantago. Based primarily on plastid data and nuclear ITS, Aragoa is sister to a clade of Littorella + Plantago, but Plantagineae relationships have yet to be assessed using multigene datasets from the nuclear and mitochondrial genomes. Here, complete mitogenomes were assembled for two species of Aragoa (A. abietina and A. cleefii). The mitogenomes of both species have a typical suite of genes for 34 proteins, 17 tRNAs, and three rRNAs. The A. abietina mitogenome assembled into a simple circular map, with no large repeats capable of producing alternative isoforms. The A. cleefii mitogenomic map was more complex, involving two circular maps bridged by a substoichiometric linear fragment. Phylogenetics of three mitochondrial genes or the nuclear rRNA repeat placed Aragoa as sister to Littorella + Plantago, consistent with previous studies. However, P. nubicola, the sole representative of subg. Bougueria, was nested within subg. Psyllium based on the mitochondrial and nuclear data, conflicting with plastid-based analyses. Phylogenetics of the nuclear rRNA repeat provided better resolution overall, whereas relationships from mitochondrial data were hindered by extensive substitution rate variation among lineages.

Complete sequence of mitochondrial genome of Cucumis melo L

Cucumis melo L. is one of the most important fruit-type vegetables in the world. This genome is divided into a main loop and two small loops. The length of the main loop is 2,709,526 bp, and the two small loops are 149,555 bp and 47,592 bp long, respectively. There are 88 coding genes in the melon mitochondrial genome, including 40 protein-coding genes (which accounted for about 1.23% of the whole genome), 8 rRNAs, and 40 tRNAs. The total length of rRNAs and tRNAs spans 0.31% of the total genome sequence. Among the 88 mitochondrial coding genes, only 5 tRNAs were located into the second largest circular DNA molecule. The complete mitogenome sequence provided herein would help understand C. melo evolution.

Complete Sequence, Multichromosomal Architecture and Transcriptome Analysis of the Solanum tuberosum Mitochondrial Genome

Mitochondrial genomes (mitogenomes) in higher plants can induce cytoplasmic male sterility and be somehow involved in nuclear-cytoplasmic interactions affecting plant growth and agronomic performance. They are larger and more complex than in other eukaryotes, due to their recombinogenic nature. For most plants, the mitochondrial DNA (mtDNA) can be represented as a single circular chromosome, the so-called master molecule, which includes repeated sequences that recombine frequently, generating sub-genomic molecules in various proportions. Based on the relevance of the potato crop worldwide, herewith we report the complete mtDNA sequence of two S. tuberosum cultivars, namely Cicero and Désirée, and a comprehensive study of its expression, based on high-coverage RNA sequencing data. We found that the potato mitogenome has a multi-partite architecture, divided in at least three independent molecules that according to our data should behave as autonomous chromosomes. Inter-cultivar variability was null, while comparative analyses with other species of the Solanaceae family allowed the investigation of the evolutionary history of their mitogenomes. The RNA-seq data revealed peculiarities in transcriptional and post-transcriptional processing of mRNAs. These included co-transcription of genes with open reading frames that are probably expressed, methylation of an rRNA at a position that should impact translation efficiency and extensive RNA editing, with a high proportion of partial editing implying frequent mis-targeting by the editing machinery.

Current understanding of male sterility systems in vegetable Brassicas and their exploitation in hybrid breeding

Overview of the current status of GMS and CMS systems available in Brassica vegetables, their molecular mechanism, wild sources of sterile cytoplasm and exploitation of male sterility in hybrid breeding. The predominantly herbaceous family Brassicaceae (crucifers or mustard family) encompasses over 3700 species, and many of them are scientifically and economically important. The genus Brassica is an economically important genus within the tribe Brassicaceae that comprises important vegetable, oilseed and fodder crops. Brassica vegetables display strong hybrid vigor, and heterosis breeding is the integral part in their improvement. Commercial production of F₁ hybrid seeds in Brassica vegetables requires an effective male sterility system. Among the available male sterility systems, cytoplasmic male sterility (CMS) is the most widely exploited in Brassica vegetables. This system is maternally inherited and studied intensively. A limited number of reports about the genic male sterility (GMS) are available in Brassica vegetables. The GMS system is reported to be dominant, recessive and trirecessive in nature in different species. In this review, we discuss the available male sterility systems in Brassica vegetables and their potential use in hybrid breeding. The molecular mechanism of mt-CMS and causal mitochondrial genes of CMS has been discussed in detail. Finally, the exploitation of male sterility system in heterosis breeding of Brassica vegetables, future prospects and need for further understanding of these systems are highlighted.

Distribution of plastids and mitochondria during male gametophyte formation in Tinantia erecta (Jacq.) Fenzl

During meiosis in microsporogenesis, autonomous cellular organelles, i.e., plastids and mitochondria, move and separate into daughter cells according to a specific pattern. This process called chondriokinesis is characteristic for a given plant species. The key criterion for classification of the chondriokinesis types was the arrangement of cell organelles during two meiosis phases: metaphase I and telophase I. The autonomous organelles participate in cytoplasmic inheritance; therefore, their precise distribution to daughter cells determines formation of identical viable microspores. In this study, the course of chondriokinesis during the development of the male gametophyte in Tinantia erecta was analyzed. The study was conducted using optical and transmission electron microscopes. During microsporogenesis in T. erecta, autonomous cell organelles moved in a manner defined as a neutral-equatorial type of chondriokinesis. Therefore, metaphase I plastids and mitochondria were evenly dispersed around the metaphase plate and formed an equatorial plate between the daughter nuclei in early telophase I. Changes in the ultrastructure of plastids and mitochondria during pollen microsporogenesis were also observed.

Genomic technologies for Hevea breeding

The commercial production of high quality natural rubber (NR) solely depends on Hevea brasiliensis Muell. Arg, (Para rubber tree) and accounts for >98% of total production worldwide. NR with its unique properties is an essential commodity for the automobile industry and its synthetic counterparts are in no way substitute to it. The rubber tree genome is very complex and plays an important role in delivering the unique properties of Hevea. But a lack of knowledge on the molecular mechanisms of rubber biosynthesis, disease resistance, etc., in elite clones of rubber still persists. Marker-assisted selection and transgenic techniques were proved to be advantageous in improving the breeding efficiency for latex yield, disease resistance, etc. The suppression subtractive hybridization (SSH), in the form of subtracted cDNA libraries and microarrays, can assist in searching the functions of expressed genes (candidate gene approach). Expressed sequence tags (ESTs) related to various metabolic aspects are well utilized to create EST banks that broadly represent the genes expressed in one tissue, such as latex cells, that assists in the study of gene function and regulation. Transcriptome analysis and gene mapping have been accomplished in Hevea at various stages. However, a selection criterion to delineate high yielding genotypes at the juvenile stage has not been accomplished so far. This is the main pit fall for rubber breeding apart from stock-scion interactions leading to yield differences among a clonally multiplied population. At least four draft genome sequences have been published on Hevea rubber, and all give different genome size and contig lengths-a comprehensive and acceptable genomic map remains unfulfilled. The progress made in molecular markers, latex biosynthesis genes, transcriptome analysis, chloroplast and mitochondrial DNA diversity, paternity identification through Breeding without Breeding (BwB), stimulated latex production and its molecular intricacies, molecular biology of tapping panel dryness, genomics for changed climates and genome mapping are discussed in this review. These information can be utilized to improvise the molecular breeding programs of Hevea in future.

Decoding and analysis of organelle genomes of Indian tea (Camellia assamica) for phylogenetic confirmation

The NCBI database has >15 chloroplast (cp) genome sequences available for different Camellia species but none for C. assamica. There is no report of any mitochondrial (mt) genome in the Camellia genus or Theaceae family. With the strong believes that these organelle genomes can play a great tool for taxonomic and phylogenetic analysis, we successfully assembled and analyzed cp and mt genome of C. assamica. We assembled the complete mt genome of C. assamica in a single circular contig of 707,441 bp length comprising of a total of 66 annotated genes, including 35 protein-coding genes, 29 tRNAs and two rRNAs. The first ever cp genome of C. assamica resulted in a circular contig of 157,353 bp length with a typical quadripartite structure. Phylogenetic analysis based on these organelle genomes showed that C. assamica was closely related to C. sinensis and C. leptophylla. It also supports Caryophyllales as Superasterids.

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.

PacBio-Based Mitochondrial Genome Assembly of Leucaena trichandra (Leguminosae) and an Intrageneric Assessment of Mitochondrial RNA Editing

Reconstructions of vascular plant mitochondrial genomes (mt-genomes) are notoriously complicated by rampant recombination that has resulted in comparatively few plant mt-genomes being available. The dearth of plant mitochondrial resources has limited our understanding of mt-genome structural diversity, complex patterns of RNA editing, and the origins of novel mt-genome elements. Here, we use an efficient long read (PacBio) iterative assembly pipeline to generate mt-genome assemblies for Leucaena trichandra (Leguminosae: Caesalpinioideae: mimosoid clade), providing the first assessment of non-papilionoid legume mt-genome content and structure to date. The efficiency of the assembly approach facilitated the exploration of alternative structures that are common place among plant mitochondrial genomes. A compact version (729 kbp) of the recovered assemblies was used to investigate sources of mt-genome size variation among legumes and mt-genome sequence similarity to the legume associated root holoparasite Lophophytum. The genome and an associated suite of transcriptome data from select species of Leucaena permitted an in-depth exploration of RNA editing in a diverse clade of closely related species that includes hybrid lineages. RNA editing in the allotetraploid, Leucaena leucocephala, is consistent with co-option of nearly equal maternal and paternal C-to-U edit components, generating novel combinations of RNA edited sites. A preliminary investigation of L. leucocephala C-to-U edit frequencies identified the potential for a hybrid to generate unique pools of alleles from parental variation through edit frequencies shared with one parental lineage, those intermediate between parents, and transgressive patterns.

Multichromosomal structure of the onion mitochondrial genome and a transcript analysis

The structures of plant mitochondrial genomes are more complex than those of animals. One of the reasons for this is that plant mitochondrial genomes typically have many long and short repeated sequences and intra- and intermolecular recombination may create various DNA molecules in this organelle. Recombination may sometimes create a novel gene that causes cytoplasmic male sterility (CMS). The onion has several cytoplasm types, with some causing CMS while others do not. The complete mitochondrial genome sequence of the onion was reported for an inbred line with CMS-S cytoplasm; however, the number of differences between onion strains remains unclear, and studies on purified mitochondrial DNA (mtDNA) have not yet been performed. Furthermore, analyses of transcripts in the mitochondrial genome have not been conducted. In the present study, we examined the mitochondrial genome of the onion variety "Momiji-3" (Allium cepa L.) possessing CMS-S-type cytoplasm using next-generation sequencing (NGS). The "Momiji-3" mitochondrial genome mainly exists as three circles as a result of recombination through repeated sequences and we herein succeeded for the first time in visualizing its structure using pulsed field gel electrophoresis (PFGE). The ability to clarify the structure of the mitochondrial genome is rare in plant mitochondria; therefore, "Momiji-3" represents a good example for elucidating complex plant mitochondrial genomes. We also mapped transcript data to the mitochondrial genome in order to identify the RNA-editing positions in all gene-coding regions and estimate the expression levels of genes. We identified 635 editing positions in gene-coding regions. Start and stop codons were created by RNA editing in six genes (nad1, nad4L, atp6, atp9, ccmFC, and orf725). The transcript amounts of novel open reading frames (ORFs) were all markedly lower than those of functional genes. These results suggest that a new functional gene was not present in the mitochondrial genome of "Momiji-3", and that the candidate gene for CMS is orf725, as previously reported.

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.

Analysis of the Complete Mitochondrial Genome Sequence of the Diploid Cotton Gossypium raimondii by Comparative Genomics Approaches

Cotton is one of the most important economic crops and the primary source of natural fiber and is an important protein source for animal feed. The complete nuclear and chloroplast (cp) genome sequences of G. raimondii are already available but not mitochondria. Here, we assembled the complete mitochondrial (mt) DNA sequence of G. raimondii into a circular genome of length of 676,078 bp and performed comparative analyses with other higher plants. The genome contains 39 protein-coding genes, 6 rRNA genes, and 25 tRNA genes. We also identified four larger repeats (63.9 kb, 10.6 kb, 9.1 kb, and 2.5 kb) in this mt genome, which may be active in intramolecular recombination in the evolution of cotton. Strikingly, nearly all of the G. raimondii mt genome has been transferred to nucleus on Chr1, and the transfer event must be very recent. Phylogenetic analysis reveals that G. raimondii, as a member of Malvaceae, is much closer to another cotton (G. barbadense) than other rosids, and the clade formed by two Gossypium species is sister to Brassicales. The G. raimondii mt genome may provide a crucial foundation for evolutionary analysis, molecular biology, and cytoplasmic male sterility in cotton and other higher plants.

Pentatricopeptide repeat 336 as the candidate gene for paternal sorting of mitochondria (Psm) in cucumber

Pentatricopeptide repeat (PPR) 336 was identified as the candidate gene for Paternal Sorting of Mitochondria ( Psm ), a nuclear locus that affects the predominant mitochondria transmitted to progenies. Cucumber (Cucumis sativus L.) is a useful plant to study organellar-nuclear interactions because its organelles show differential transmission, maternal for chloroplasts and paternal for mitochondria. The mitochondrial DNA (mtDNA) of cucumber is relatively large due in part to accumulation of repetitive DNAs and recombination among these repetitive regions produces structurally polymorphic mtDNAs associated with paternally transmitted mosaic (MSC) phenotypes. The mitochondrial mutant MSC16 possesses an under-representation of ribosomal protein S7 (rps7), a key component of the small ribosomal subunit in the mitochondrion. A nuclear locus, Paternal Sorting of Mitochondria (Psm), affects the predominant mitochondria transmitted to progenies generated from crosses with MSC16 as the male parent. Using single nucleotide polymorphisms, Psm was mapped to a 170 kb region on chromosome 3 of cucumber and pentatricopeptide repeat (PPR) 336 was identified as the likely candidate gene. PPR336 stabilizes mitochondrial ribosomes in Arabidopsis thaliana and because MSC16 shows reduced transcription of rps7, the cucumber homolog of PPR336 (CsPPR336) as the candidate for Psm is consistent with a nuclear effect on ribosome assembly or stability in the mitochondrion. We used polymorphisms in CsPPR336 to genotype progenies segregating at Psm and recovered only one Psm -/- plant with the MSC phenotype, indicating that the combination of the Psm- allele with mitochondria from MSC16 is almost always lethal. This research illustrates the usefulness of the MSC mutants of cucumber to reveal and study unique interactions between the mitochondrion and nucleus.

The two chromosomes of the mitochondrial genome of a sugarcane cultivar: assembly and recombination analysis using long PacBio reads

Sugarcane accounts for a large portion of the worlds sugar production. Modern commercial cultivars are complex hybrids of S. officinarum and several other Saccharum species. Historical records identify New Guinea as the origin of S. officinarum and that a small number of plants originating from there were used to generate all modern commercial cultivars. The mitochondrial genome can be a useful way to identify the maternal origin of commercial cultivars. We have used the PacBio RSII to sequence and assemble the mitochondrial genome of a South East Asian commercial cultivar, known as Khon Kaen 3. The long read length of this sequencing technology allowed for the mitochondrial genome to be assembled into two distinct circular chromosomes with all repeat sequences spanned by individual reads. Comparison of five commercial hybrids, two S. officinarum and one S. spontaneum to our assembly reveals no structural rearrangements between our assembly, the commercial hybrids and an S. officinarum from New Guinea. The S. spontaneum, from India, and one sample of S. officinarum (unknown origin) are substantially rearranged and have a large number of homozygous variants. This supports the record that S. officinarum plants from New Guinea are the maternal source of all modern commercial hybrids.

Mitochondrial Genome Analysis of Wild Rice (Oryza minuta) and Its Comparison with Other Related Species

Oryza minuta (Poaceae family) is a tetraploid wild relative of cultivated rice with a BBCC genome. O. minuta has the potential to resist against various pathogenic diseases such as bacterial blight (BB), white backed planthopper (WBPH) and brown plant hopper (BPH). Here, we sequenced and annotated the complete mitochondrial genome of O. minuta. The mtDNA genome is 515,022 bp, containing 60 protein coding genes, 31 tRNA genes and two rRNA genes. The mitochondrial genome organization and the gene content at the nucleotide level are highly similar (89%) to that of O. rufipogon. Comparison with other related species revealed that most of the genes with known function are conserved among the Poaceae members. Similarly, O. minuta mt genome shared 24 protein-coding genes, 15 tRNA genes and 1 ribosomal RNA gene with other rice species (indica and japonica). The evolutionary relationship and phylogenetic analysis revealed that O. minuta is more closely related to O. rufipogon than to any other related species. Such studies are essential to understand the evolutionary divergence among species and analyze common gene pools to combat risks in the current scenario of a changing environment.

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.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

The Mosaic Mutants of Cucumber: A Method to Produce Knock-Downs of Mitochondrial Transcripts

Cytoplasmic effects on plant performance are well-documented and result from the intimate interaction between organellar and nuclear gene products. In plants, deletions, mutations, or chimerism of mitochondrial genes are often associated with deleterious phenotypes, as well as economically important traits such as cytoplasmic male sterility used to produce hybrid seed. Presently, genetic analyses of mitochondrial function and nuclear interactions are limited because there is no method to efficiently produce mitochondrial mutants. Cucumber (Cucumis sativus L.) possesses unique attributes useful for organellar genetics, including differential transmission of the three plant genomes (maternal for plastid, paternal for mitochondrial, and bi-parental for nuclear), a relatively large mitochondrial DNA in which recombination among repetitive motifs produces rearrangements, and the existence of strongly mosaic (MSC) paternally transmitted phenotypes that appear after passage of wild-type plants through cell cultures and possess unique rearrangements in the mitochondrial DNA. We sequenced the mitochondrial DNA from three independently produced MSC lines and revealed under-represented regions and reduced transcription of mitochondrial genes carried in these regions relative to the wild-type parental line. Mass spectrometry and Western blots did not corroborate transcriptional differences in the mitochondrial proteome of the MSC mutant lines, indicating that post-transcriptional events, such as protein longevity, may compensate for reduced transcription in MSC mitochondria. Our results support cucumber as a model system to produce transcriptional "knock-downs" of mitochondrial genes useful to study mitochondrial responses and nuclear interactions important for plant performance.

Complete mitochondrial genome sequence of black mustard (Brassica nigra; BB) and comparison with Brassica oleracea (CC) and Brassica carinata (BBCC)

Crop species of Brassica (Brassicaceae) consist of three monogenomic species and three amphidiploid species resulting from interspecific hybridizations among them. Until now, mitochondrial genome sequences were available for only five of these species. We sequenced the mitochondrial genome of the sixth species, Brassica nigra (nuclear genome constitution BB), and compared it with those of Brassica oleracea (CC) and Brassica carinata (BBCC). The genome was assembled into a 232 145 bp circular sequence that is slightly larger than that of B. oleracea (219 952 bp). The genome of B. nigra contained 33 protein-coding genes, 3 rRNA genes, and 17 tRNA genes. The cox2-2 gene present in B. oleracea was absent in B. nigra. Although the nucleotide sequences of 52 genes were identical between B. nigra and B. carinata, the second exon of rps3 showed differences including an insertion/deletion (indel) and nucleotide substitutions. A PCR test to detect the indel revealed intraspecific variation in rps3, and in one line of B. nigra it amplified a DNA fragment of the size expected for B. carinata. In addition, the B. carinata lines tested here produced DNA fragments of the size expected for B. nigra. The results indicate that at least two mitotypes of B. nigra were present in the maternal parents of B. carinata.

History of plastid DNA insertions reveals weak deletion and at mutation biases in angiosperm mitochondrial genomes

Angiosperm mitochondrial genomes exhibit many unusual properties, including heterogeneous nucleotide composition and exceptionally large and variable genome sizes. Determining the role of nonadaptive mechanisms such as mutation bias in shaping the molecular evolution of these unique genomes has proven challenging because their dynamic structures generally prevent identification of homologous intergenic sequences for comparative analyses. Here, we report an analysis of angiosperm mitochondrial DNA sequences that are derived from inserted plastid DNA (mtpts). The availability of numerous completely sequenced plastid genomes allows us to infer the evolutionary history of these insertions, including the specific nucleotide substitutions and indels that have occurred because their incorporation into the mitochondrial genome. Our analysis confirmed that many mtpts have a complex history, including frequent gene conversion and multiple examples of horizontal transfer between divergent angiosperm lineages. Nevertheless, it is clear that the majority of extant mtpt sequence in angiosperms is the product of recent transfer (or gene conversion) and is subject to rapid loss/deterioration, suggesting that most mtpts are evolving relatively free from functional constraint. The evolution of mtpt sequences reveals a pattern of biased mutational input in angiosperm mitochondrial genomes, including an excess of small deletions over insertions and a skew toward nucleotide substitutions that increase AT content. However, these mutation biases are far weaker than have been observed in many other cellular genomes, providing insight into some of the notable features of angiosperm mitochondrial architecture, including the retention of large intergenic regions and the relatively neutral GC content found in these regions.

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.

Complete mitochondrial genome of Eruca sativa Mill. (Garden rocket)

Eruca sativa (Cruciferae family) is an ancient crop of great economic and agronomic importance. Here, the complete mitochondrial genome of Eruca sativa was sequenced and annotated. The circular molecule is 247 696 bp long, with a G+C content of 45.07%, containing 33 protein-coding genes, three rRNA genes, and 18 tRNA genes. The Eruca sativa mitochondrial genome may be divided into six master circles and four subgenomic molecules via three pairwise large repeats, resulting in a more dynamic structure of the Eruca sativa mtDNA compared with other cruciferous mitotypes. Comparison with the Brassica napus MtDNA revealed that most of the genes with known function are conserved between these two mitotypes except for the ccmFN2 and rrn18 genes, and 27 point mutations were scattered in the 14 protein-coding genes. Evolutionary relationships analysis suggested that Eruca sativa is more closely related to the Brassica species and to Raphanus sativus than to Arabidopsis thaliana.

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 11 Brassicales mitochondrial genomes and the mitochondrial transcriptome of Brassica oleracea

To elucidate the evolution of mitochondrial genomic diversity within a single order of angiosperms, we sequenced seven Brassicales genomes and the transcriptome of Brassica oleracea. In the common ancestor of Brassicaceae, several genes of known function were lost and the ccmFN gene was split into two independent genes, which also coincides with a trend of genome reduction towards the smallest sequenced angiosperm genomes of Brassica. For most ORFs of unknown function, the lack of conservation throughout Brassicales and the generally low expression and absence of RNA editing in B. oleracea argue against functionality. However, two chimeric ORFs were expressed and edited in B. oleracea, suggesting a potential role in cytoplasmic male sterility in certain nuclear backgrounds. These results demonstrate how frequent shifts in size, structure, and content of plant mitochondrial genomes can occur over short evolutionary time scales.

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.

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 complete mitochondrial genome of Gossypium hirsutum and evolutionary analysis of higher plant mitochondrial genomes

BACKGROUND

Mitochondria are the main manufacturers of cellular ATP in eukaryotes. The plant mitochondrial genome contains large number of foreign DNA and repeated sequences undergone frequently intramolecular recombination. Upland Cotton (Gossypium hirsutum L.) is one of the main natural fiber crops and also an important oil-producing plant in the world. Sequencing of the cotton mitochondrial (mt) genome could be helpful for the evolution research of plant mt genomes.

METHODOLOGY/PRINCIPAL FINDINGS

We utilized 454 technology for sequencing and combined with Fosmid library of the Gossypium hirsutum mt genome screening and positive clones sequencing and conducted a series of evolutionary analysis on Cycas taitungensis and 24 angiosperms mt genomes. After data assembling and contigs joining, the complete mitochondrial genome sequence of G. hirsutum was obtained. The completed G.hirsutum mt genome is 621,884 bp in length, and contained 68 genes, including 35 protein genes, four rRNA genes and 29 tRNA genes. Five gene clusters are found conserved in all plant mt genomes; one and four clusters are specifically conserved in monocots and dicots, respectively. Homologous sequences are distributed along the plant mt genomes and species closely related share the most homologous sequences. For species that have both mt and chloroplast genome sequences available, we checked the location of cp-like migration and found several fragments closely linked with mitochondrial genes.

CONCLUSION

The G. hirsutum mt genome possesses most of the common characters of higher plant mt genomes. The existence of syntenic gene clusters, as well as the conservation of some intergenic sequences and genic content among the plant mt genomes suggest that evolution of mt genomes is consistent with plant taxonomy but independent among different species.

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.

The genome and transcriptome of perennial ryegrass mitochondria

Background

Perennial ryegrass (Lolium perenne L.) is one of the most important forage and turf grass species of temperate regions worldwide. Its mitochondrial genome is inherited maternally and contains genes that can influence traits of agricultural importance. Moreover, the DNA sequence of mitochondrial genomes has been established and compared for a large number of species in order to characterize evolutionary relationships. Therefore, it is crucial to understand the organization of the mitochondrial genome and how it varies between and within species. Here, we report the first de novo assembly and annotation of the complete mitochondrial genome from perennial ryegrass.

Results

Intact mitochondria from perennial ryegrass leaves were isolated and used for mtDNA extraction. The mitochondrial genome was sequenced to a 167-fold coverage using the Roche 454 GS-FLX Titanium platform, and assembled into a circular master molecule of 678,580 bp. A total of 34 proteins, 14 tRNAs and 3 rRNAs are encoded by the mitochondrial genome, giving a total gene space of 48,723 bp (7.2%). Moreover, we identified 149 open reading frames larger than 300 bp and covering 67,410 bp (9.93%), 250 SSRs, 29 tandem repeats, 5 pairs of large repeats, and 96 pairs of short inverted repeats. The genes encoding subunits of the respiratory complexes – nad1 to nad9, cob, cox1 to cox3 and atp1 to atp9 – all showed high expression levels both in absolute numbers and after normalization.

Conclusions

The circular master molecule of the mitochondrial genome from perennial ryegrass presented here constitutes an important tool for future attempts to compare mitochondrial genomes within and between grass species. Our results also demonstrate that mitochondria of perennial ryegrass contain genes crucial for energy production that are well conserved in the mitochondrial genome of monocotyledonous species. The expression analysis gave us first insights into the transcriptome of these mitochondrial genes in perennial ryegrass.

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.

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.

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.

Comparative analysis of mitochondrial genomes between a wheat K-type cytoplasmic male sterility (CMS) line and its maintainer line

BACKGROUND

Plant mitochondria, semiautonomous organelles that function as manufacturers of cellular ATP, have their own genome that has a slow rate of evolution and rapid rearrangement. Cytoplasmic male sterility (CMS), a common phenotype in higher plants, is closely associated with rearrangements in mitochondrial DNA (mtDNA), and is widely used to produce F1 hybrid seeds in a variety of valuable crop species. Novel chimeric genes deduced from mtDNA rearrangements causing CMS have been identified in several plants, such as rice, sunflower, pepper, and rapeseed, but there are very few reports about mtDNA rearrangements in wheat. In the present work, we describe the mitochondrial genome of a wheat K-type CMS line and compare it with its maintainer line.

RESULTS

The complete mtDNA sequence of a wheat K-type (with cytoplasm of Aegilops kotschyi) CMS line, Ks3, was assembled into a master circle (MC) molecule of 647,559 bp and found to harbor 34 known protein-coding genes, three rRNAs (18 S, 26 S, and 5 S rRNAs), and 16 different tRNAs. Compared to our previously published sequence of a K-type maintainer line, Km3, we detected Ks3-specific mtDNA (> 100 bp, 11.38%) and repeats (> 100 bp, 29 units) as well as genes that are unique to each line: rpl5 was missing in Ks3 and trnH was absent from Km3. We also defined 32 single nucleotide polymorphisms (SNPs) in 13 protein-coding, albeit functionally irrelevant, genes, and predicted 22 unique ORFs in Ks3, representing potential candidates for K-type CMS. All these sequence variations are candidates for involvement in CMS. A comparative analysis of the mtDNA of several angiosperms, including those from Ks3, Km3, rice, maize, Arabidopsis thaliana, and rapeseed, showed that non-coding sequences of higher plants had mostly divergent multiple reorganizations during the mtDNA evolution of higher plants.

CONCLUSION

The complete mitochondrial genome of the wheat K-type CMS line Ks3 is very different from that of its maintainer line Km3, especially in non-coding sequences. Sequence rearrangement has produced novel chimeric ORFs, which may be candidate genes for CMS. Comparative analysis of several angiosperm mtDNAs indicated that non-coding sequences are the most frequently reorganized during mtDNA evolution in higher plants.

The mitochondrial genome of the legume Vigna radiata and the analysis of recombination across short mitochondrial repeats

The mitochondrial genomes of seed plants are exceptionally fluid in size, structure, and sequence content, with the accumulation and activity of repetitive sequences underlying much of this variation. We report the first fully sequenced mitochondrial genome of a legume, Vigna radiata (mung bean), and show that despite its unexceptional size (401,262 nt), the genome is unusually depauperate in repetitive DNA and "promiscuous" sequences from the chloroplast and nuclear genomes. Although Vigna lacks the large, recombinationally active repeats typical of most other seed plants, a PCR survey of its modest repertoire of short (38–297 nt) repeats nevertheless revealed evidence for recombination across all of them. A set of novel control assays showed, however, that these results could instead reflect, in part or entirely, artifacts of PCR-mediated recombination. Consequently, we recommend that other methods, especially high-depth genome sequencing, be used instead of PCR to infer patterns of plant mitochondrial recombination. The average-sized but repeat- and feature-poor mitochondrial genome of Vigna makes it ever more difficult to generalize about the factors shaping the size and sequence content of plant mitochondrial genomes.

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.

Extensive loss of translational genes in the structurally dynamic mitochondrial genome of the angiosperm Silene latifolia

Background

Mitochondrial gene loss and functional transfer to the nucleus is an ongoing process in many lineages of plants, resulting in substantial variation across species in mitochondrial gene content. The Caryophyllaceae represents one lineage that has experienced a particularly high rate of mitochondrial gene loss relative to other angiosperms.

Results

In this study, we report the first complete mitochondrial genome sequence from a member of this family, Silene latifolia. The genome can be mapped as a 253,413 bp circle, but its structure is complicated by a large repeated region that is present in 6 copies. Active recombination among these copies produces a suite of alternative genome configurations that appear to be at or near "recombinational equilibrium". The genome contains the fewest genes of any angiosperm mitochondrial genome sequenced to date, with intact copies of only 25 of the 41 protein genes inferred to be present in the common ancestor of angiosperms. As observed more broadly in angiosperms, ribosomal proteins have been especially prone to gene loss in the S. latifolia lineage. The genome has also experienced a major reduction in tRNA gene content, including loss of functional tRNAs of both native and chloroplast origin. Even assuming expanded wobble-pairing rules, the mitochondrial genome can support translation of only 17 of the 61 sense codons, which code for only 9 of the 20 amino acids. In addition, genes encoding 18S and, especially, 5S rRNA exhibit exceptional sequence divergence relative to other plants. Divergence in one region of 18S rRNA appears to be the result of a gene conversion event, in which recombination with a homologous gene of chloroplast origin led to the complete replacement of a helix in this ribosomal RNA.

Conclusions

These findings suggest a markedly expanded role for nuclear gene products in the translation of mitochondrial genes in S. latifolia and raise the possibility of altered selective constraints operating on the mitochondrial translational apparatus in this lineage.

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.

Heteroplasmy and stoichiometric complexity of plant mitochondrial genomes--though this be madness, yet there's method in't

Mitochondrial heteroplasmy is defined as the coexistence of divergent mitochondrial genotypes in a cell. The ratio of the alternative genomes may be variable, but in plants, the usually prevalent main genome is accompanied by sublimons--substoichiometric mitochondrial DNA (mtDNA) molecules. Plant mitochondrial heteroplasmy was originally viewed as being associated with pathological mutations or was found in non-natural plant populations. Currently, it is considered to be a common situation in plants. Recent years have changed the previous view on the role of homologous recombination, small-scale mutations, and paternal leakage of mtDNA in the generation of heteroplasmy. Newly developed sensitive techniques have allowed the precise estimation of mtDNA stoichiometry. Mechanisms of maintenance and transmission of heteroplasmic genomes, including DNA recombination and replication, as well as mitochondrial fusion and fission, have been studied. This review describes the high level of plant mitochondrial genome complication--the 'madness' resulting from the heteroplasmic state and explains the method hidden in this madness. Heteroplasmy is described as the evolutionary strategy of uniparentally inherited plant mitochondrial genomes which do not undergo sexual recombination. In order to compensate for this deficiency, alternative types of mtDNA are substoichiometrically accumulated as a reservoir of genetic variability and may undergo accelerated evolution. Occasionally, sublimons are selected and amplified in the process called substoichiometric shifting, to take over the role of the main genome. Alternative mitochondrial genomes may recombine, yielding new mtDNA variants, or segregate during plant growth resulting in plants with mosaic phenotypes. Two opposite roles of mitochondrial heteroplasmy with respect to acceleration or counteracting of mutation accumulation are also discussed. Finally, nuclear control of heteroplasmy and substoichiometric shifting is described.

Comparison of mitochondrial and chloroplast genome segments from three onion (Allium cepa L.) cytoplasm types and identification of a trans-splicing intron of cox2

To study genetic relatedness of two male sterility-inducing cytotypes, the phylogenetic relationship among three cytotypes of onions (Allium cepa L.) was assessed by analyzing polymorphisms of the mitochondrial DNA organization and chloroplast sequences. The atp6 gene and a small open reading frame, orf22, did not differ between the normal and CMS-T cytotypes, but two SNPs and one 4-bp insertion were identified in CMS-S cytotype. Partial sequences of the chloroplast ycf2 gene were integrated in the upstream sequence of the cob gene via short repeat sequence-mediated recombination. However, this chloroplast DNA-integrated organization was detected only in CMS-S. Interestingly, disruption of a group II intron of cox2 was identified for the first time in this study. Like other trans-splicing group II introns in mitochondrial genomes, fragmentation of the intron occurred in domain IV. Two variants of each exon1 and exon2 flanking sequences were identified. The predominant types of four variants were identical in both the normal and the CMS-T cytotypes. These predominant types existed as sublimons in CMS-S cytotypes. Altogether, no differences were identified between normal and CMS-T, but significant differences in gene organization and nucleotide sequences were identified in CMS-S, suggesting recent origin of CMS-T male-sterility from the normal cytotype.

Identification of highly variable chloroplast sequences and development of cpDNA-based molecular markers that distinguish four cytoplasm types in radish (Raphanus sativus L.)

Four types of cytoplasms (Ogura, DCGMS, DBRMF1, and DBRMF2) were identified in the previous studies using molecular markers based on mitochondrial genome variations in radish (Raphanus sativus L.). However, mtDNA markers have limitations in obtaining clear results due to complexity of radish mitochondrial genomes. To improve fidelity, molecular markers based on variation of chloroplast genome sequences were developed in this study. We searched for the sequence variations of chloroplast genome among the four cytoplasm types in 11 noncoding intergenic regions of ~8.7 kb. Highly variable intergenic regions between trnK and rps16 were identified, and a couple of 4-34 bp indels were used to develop a simple PCR-based marker that distinguished the four cytoplasm types based on the PCR product length polymorphism. Two additional cpDNA markers were developed by using a single nucleotide polymorphism and 17-bp insertion. Analysis of 90 accessions using both mtDNA and cpDNA markers showed the perfect match of results of both the markers, suggesting strict co-transmission of mitochondria and chloroplast in radish. Phylogenetic trees showed that two male-sterility inducing cytoplasms, Ogura and DCGMS, were closely related to DBRMF1 and DBRMF2, respectively. Analysis of 120 radish germplasms introduced from diverse countries showed that the frequency of male-sterility inducing mitotypes of Ogura and DCGMS was very low, and DCGMS was predominately detected in eastern European countries. Majority of accessions from Europe and Asia were shown to contain DBRMF2 and DBRMF1 mitotypes, respectively.

Identification of a novel chimeric gene, orf725, and its use in development of a molecular marker for distinguishing among three cytoplasm types in onion (Allium cepa L.)

A novel chimeric gene with a 5' end containing the nearly complete sequence of the coxI gene and a 3' end showing homology with chive orfA501 was isolated by genome walking from two cytoplasm types: CMS-S and CMS-T, both of which induce male-sterility in onion (Allium cepa L.). In addition, the normal active and variant inactive coxI genes were also isolated from onions containing the normal and CMS-S cytoplasms, respectively. The chimeric gene, designated as orf725, was nearly undetectable in normal cytoplasm, and the copy number of the normal coxI gene was significantly reduced in CMS-S cytoplasm. RT-PCR results showed that orf725 was not transcribed in normal cytoplasm. Meanwhile, the normal coxI gene, which is essential for normal mitochondrial function, was not expressed in CMS-S cytoplasm. However, both orf725 and coxI were transcribed in CMS-T cytoplasm. The expression of orf725, a putative male-sterility-inducing gene, was not affected by the presence of nuclear restorer-of-fertility gene(s) in male-fertility segregating populations originating from the cross between a male-sterile plant containing either CMS-T or CMS-S and a male-fertile plant whose genotypes of nuclear restorer gene(s) might be heterozygous. The specific stoichiometry of orf725 and coxI in the mtDNA of the three cytoplasm types was consistent among diverse germplasm. Therefore, a molecular marker based on the relative copy numbers of orf725 and coxI was designed for distinguishing among the three cytoplasm types by one simple PCR. The reliability and applicability of the molecular marker was shown by testing diverse onion germplasm.

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

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

Discovery of a novel cytoplasmic male-sterility and its restorer lines in radish (Raphanus sativus L.)

A male-sterile (MS) radish (Raphanus sativus L.) was found in an accession collected from Uzbekistan. Unlike Ogura MS radishes in which no pollen grain is typically visible during anthesis, a small number of pollen grains stuck together in the dehiscing anthers was observed in the newly identified MS radish. Fluorescein diacetate tests and scanning electron micrographs showed that pollen grains in the new MS radish were severely deformed and non-viable. Cytological examination of pollen development stages showed a clear difference in the defective stage from that seen in Ogura male-sterility. Reciprocal cross-pollination with diverse male-fertile lines indicated that pollen grains of the new MS radish were completely sterile, and the female organs were fully fertile. When the new MS radish and Ogura MS lines were cross-pollinated with a set of eight breeding lines, all F1 progeny originating from crosses with the new MS radish were male-sterile. In contrast, most of the F1 progeny resulting from crosses with Ogura MS lines were male-fertile. These results demonstrated that factors associated with induction of the newly identified male-sterility are different from those of Ogura male-sterility. The lack of restorer lines for the newly identified male-sterility led us to predict that it might be a complete cytoplasmic male-sterility without restorer-of-fertility genes in nuclear genomes. However, cross-pollination with more diverse radish germplasm identified one accession introduced from Russia that could completely restore fertility, proving the existence of restorer-of-fertility gene(s) for the new male-sterility. Meanwhile, the PCR amplification profile of molecular markers for the classification of radish mitochondrial genome types revealed that the new MS radish contained a novel mitotype.

Gene duplication and transfer events in plant mitochondria genome

Gene or genome duplication events increase the amount of genetic material available to increase the genomic, and thereby phenotypic, complexity of organisms during evolution. Gene duplication and transfer events have been important to molecular evolution in all three domains of life, and may be the first step in the emergence of new gene functions. Gene transfer events have been proposed as another accelerator of evolution. The duplicated gene or genome, mainly nuclear, has been the subject of several recent reviews. In addition to the nuclear genome, organisms have organelle genomes, including mitochondrial genome. In this review, we briefly summarize gene duplication and transfer events in the plant mitochondrial genome.

Angiosperm mitochondrial genomes and mutations

Flowering plants harbor the largest mitochondrial genomes reported so far. At present, the nucleotide sequences of 15 mitochondrial genomes from seven angiosperm species are available, making detailed comparative analysis feasible. The gene content is variable among the species, but the most striking feature is the fluidity of intergenic regions, where species-specific sequences predominate. Additionally, angiosperm mitochondrial genomes, even within a species, show a remarkable amount of rearrangement. We also review mitochondrial mutants in angiosperms from a genomic viewpoint, and discuss how they have arisen. The involvement of nuclear genes in mitochondrial genome stability and organization is currently being revealed through the analysis of mutants.