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.)

De novo assembly and characterization of the complete mitochondrial genome of Phellodendron amurense reveals three repeat-mediated recombination

Phellodendron amurense Rupr., a rare herb renowned for its medicinal and ecological significance, has remained genetically unexplored at the mitochondrial level until now. This study presents the first-ever systematic assembly and annotation of the complete mitochondrial genome of P. amurense, achieved through a hybrid strategy combining Illumina and Nanopore sequencing data. The mitochondrial genome spans 566,285 bp with a GC content of 45.51 %, structured into two circular molecules. Our comprehensive analysis identified 32 protein-coding genes (PCGs), 33 tRNA genes, and 3 rRNA genes, alongside 181 simple sequence repeats, 19 tandem repeats, and 310 dispersed repeats. Notably, multiple genome conformations were predicted due to repeat-mediated homologous recombination. Additionally, we assembled the chloroplast genome, identifying 21 mitochondrial plastid sequences that provide insights into organelle genome interactions. A total of 380 RNA-editing sites within the mitochondrial PCGs were predicted, enhancing our understanding of gene regulation and function. Phylogenetic analysis using mitochondrial PCGs from 30 species revealed evolutionary relationships, confirming the homology between P. amurense and Citrus species. This foundational study offers a valuable genetic resource for the Rutaceae family, facilitating further research into genetic evolution and molecular diversity in plant mitochondrial genomes.

Past and future of cytoplasmic male sterility and heterosis breeding in crop plants

Plant breeding needs to embrace genetic innovations to ensure stability in crop yields under fluctuating climatic conditions. Development of commercial hybrid varieties has proven to be a sustainable and economical alternative to deliver superior yield, quality and resistance with uniformity in a number of food crops. Cytoplasmic male sterility (CMS), a maternally inherited inability to produce functional pollen, facilitates a three-line system for efficient hybrid seed production strategies in crops. The CMS system has illustrated its potential as a robust pollination control mechanism to support the billion-dollar seed industry. In plants, CMS arises due to a genomic conflict between mitochondrial open reading frames (orfs) and nuclear-encoding restoration-of-fertility (Rf) genes, leading to floral abnormalities and pollen sterility. Research on pollen sterility and fertility restoration provides deeper insights into cytoplasmic-nuclear interplay in plants and elucidates key molecular targets for hybrid breeding in crops. More recently, programmable gene editing (e.g., TALEN, CRISPR-Cas) has emerged as a promising tool to functionally validate CMS and Rf genes and obviate the need for pollen donors or Rf-genes for hybrid breeding. Modern genomic prediction models have allowed establishment of high-performing heterotic groups and patterns for sustaining long-term gain in hybrid breeding. This article reviews latest discoveries elucidating the molecular mechanisms behind CMS and fertility restoration in plants. We then present our perspective on how evolving genetic technologies are contributing to advance fundamental knowledge of the CMS-Rf genetic system for producing crop hybrids with high heterosis.

Features and evolutionary adaptations of the mitochondrial genome of Garuga forrestii W. W. Sm

Introduction

Garuga forrestii W. W. Sm. is a tree species of the Burseraceae family, endemic to China, found in hot/warm-dry valleys. This species plays a crucial role in maintaining biodiversity in these ecosystems.

Methods

We performed de novo assembly of the Garuga forrestii mitochondrial genome using PMAT (v.1.5.4), resulting in a typical circular molecule of 606,853 bp. The genome consists of 31 tRNA genes, 3 rRNA genes, 35 protein-coding genes, and 1 pseudogene. The study also investigates RNA editing sites and evolutionary patterns.

Results

The mitochondrial genome exhibits a low proportion of repetitive sequences (3.30%), suggesting a highly conserved structure. A high copy number of the trnM-CAT gene (4 copies) is noted, which may contribute to genomic rearrangement and adaptive evolution. Among the 476 RNA editing sites, hydrophilic-hydrophobic and hydrophobic-hydrophobic editing events are most common, accounting for 77.10%. Negative selection predominates among most genes (Ka/Ks < 1), while a few genes (e.g., matR, nad3, rps1, rps12, and rps4) show signs of positive selection (Ka/Ks > 1), potentially conferring evolutionary advantages. Additionally, a significant A/T bias is observed at the third codon position. Phylogenomic analysis supports the APG IV classification, with no evidence of horizontal gene transfer.

Discussion

This mitochondrial genome offers valuable insights into the adaptive mechanisms and evolutionary processes of Garuga forrestii. It enhances our understanding of the species' biogeography in tropical Southeast Asia and Southwest China, providing key information on the evolutionary history of this genus.

A breeding method for Ogura CMS restorer line independent of restorer source in Brassica napus

The Ogura cytoplasmic male sterility (CMS) line of Brassica napus has gained significant attention for its use in harnessing heterosis. It remains unaffected by temperature and environment and is thorough and stable. The Ogura cytoplasmic restorer line of Brassica napus is derived from the distant hybridization of Raphanus sativus L. and B. napus, but it carried a large number of radish fragments into Brassica napus, because there is no homologous allele of the restorer gene in B. napus, transferring it becomes challenging. In this study, the double haploid induction line in B. napus was used as the male parent for hybridization with the Ogura CMS of B. napus. Surprisingly, fertile plants appeared in the offspring. Further analysis revealed that the cytoplasmic type, ploidy, and chromosome number of the fertile offspring were consistent with the sterile female parent. Moreover, the mitochondrial genome similarity between the fertile offspring and the sterile female parent was 97.7% indicates that the cytoplasm of the two is the same, while the nuclear gene difference between fertile offspring and sterile female parent was only 10.33%, indicates that new genes appeared in the offspring. To further investigate and locate the restorer gene, the BSA method was employed to construct extreme mixed pools. As a result, the restorer gene was mapped to three positions: A09 chromosome 10.99-17.20 Mb, C03 chromosome 5.07-5.34 Mb, and C09 chromosome 18.78-36.60 Mb. The experimental results have proved that induction does produce restorer genes. The induction of the Ogura CMS restorer gene through DH induction line provides a promising new approach for harnessing heterosis in B. napus.

Methods for assembling complex mitochondrial genomes in land plants

The large size and complex structural rearrangements inherent in the mitochondrial genomes of land plants pose challenges for their sequencing. Originally, the assembly of these genomes required the cloning of mitochondrial DNA fragments followed by Sanger sequencing. Subsequently, the advent of next-generation sequencing significantly expedited the process. This review highlights examples of plant mitochondrial genome assembly employing various technologies, including 454 sequencing, Illumina short sequencing reads, and Pacific Biosciences or Oxford Nanopore Technology long sequencing reads. The combination of short and long reads in hybrid assembly has proven to be the most efficient approach for achieving reliable assemblies of land plant mitochondrial genomes.

Assembly and comparative analysis of the complete mitochondrial genome of Fritillaria ussuriensis Maxim. (Liliales: Liliaceae), an endangered medicinal plant

BACKGROUND

Fritillaria ussuriensis is an endangered medicinal plant known for its notable therapeutic properties. Unfortunately, its population has drastically declined due to the destruction of forest habitats. Thus, effectively protecting F. ussuriensis from extinction poses a significant challenge. A profound understanding of its genetic foundation is crucial. To date, research on the complete mitochondrial genome of F. ussuriensis has not yet been reported.

RESULTS

The complete mitochondrial genome of F. ussuriensis was sequenced and assembled by integrating PacBio and Illumina sequencing technologies, revealing 13 circular chromosomes totaling 737,569 bp with an average GC content of 45.41%. A total of 55 genes were annotated in this mitogenome, including 2 rRNA genes, 12 tRNA genes, and 41 PCGs. The mitochondrial genome of F. ussuriensis contained 192 SSRs and 4,027 dispersed repeats. In the PCGs of F. ussuriensis mitogenome, 90.00% of the RSCU values exceeding 1 exhibited a preference for A-ended or U-ended codons. In addition, 505 RNA editing sites were predicted across these PCGs. Selective pressure analysis suggested negative selection on most PCGs to preserve mitochondrial functionality, as the notable exception of the gene nad3 showed positive selection. Comparison between the mitochondrial and chloroplast genomes of F. ussuriensis revealed 20 homologous fragments totaling 8,954 bp. Nucleotide diversity analysis revealed the variation among genes, and gene atp9 was the most notable. Despite the conservation of GC content, mitogenome sizes varied significantly among six closely related species, and colinear analysis confirmed the lack of conservation in their genomic structures. Phylogenetic analysis indicated a close relationship between F. ussuriensis and Lilium tsingtauense.

CONCLUSIONS

In this study, we sequenced and annotated the mitogenome of F. ussuriensis and compared it with the mitogenomes of other closely related species. In addition to genomic features and evolutionary position, this study also provides valuable genomic resources to further understand and utilize this medicinal plant.

Cytological and Molecular Characterization of a New Ogura Cytoplasmic Male Sterility Restorer of Brassica napus L

Ogura cytoplasmic male sterility (CMS) is considered the rapeseed (Brassica napus L.) with the most potential to be utilized as a heterosis system worldwide, but it lacks sufficient restorers. In this study, root tip cell (RTC) mitotic and pollen mother cell (PMC) meiosis observations were compared to ensure the number of chromosomes and the formation of a chromosomal bridge using restorer lines R2000, CLR650, and Zhehuhong (a new restorer) as the experimental material. Further, molecular markers of exogenous chromosomal fragments were detected and the sequence and expression differences of restorer genes in the three lines were determined to identify the distinctive characteristics of Zhehuhong. The results showed that the number of chromosomes in Zhehuhong was stable (2n = 38), indicating that the exogenous radish chromosome segment had been integrated into the chromosome of Zhehuhong. Molecular marker detection revealed that Zhehuhong was detected at most loci, with only the RMA05 locus being missed. The exogenous radish chromosome segment of Zhehuhong differed from R2000 and CLR650. The pollen mother cells of Zhehuhong showed chromosome lagging in the meiotic metaphase I, meiotic anaphase I, and meiotic anaphase II, which was consistent with R2000 and CLR650. The restorer gene PPRB in Zhehuhong had 85 SNPs compared with R2000 and 119 SNPs compared with CLR650, indicating the distinctive characteristic of PPRB in Zhehuhong. In terms of the spatial expression of PPRB, the highest level was detected in the anthers in the three restorer lines. In addition, in terms of temporal expression, the PPRB gene expression of Zhehuhong was highest at a bud length of 4 mm. Our results clearly indicated that Zhehuhong is a new restorer line for the Ogura CMS system, which can be used further in rapeseed heterosis utilization.

Assembly and comparative analysis of the complete mitochondrial genome of Brassica rapa var. Purpuraria

BACKGROUND

Purple flowering stalk (Brassica rapa var. purpuraria) is a widely cultivated plant with high nutritional and medicinal value and exhibiting strong adaptability during growing. Mitochondrial (mt) play important role in plant cells for energy production, developing with an independent genetic system. Therefore, it is meaningful to assemble and annotate the functions for the mt genome of plants independently. Though there have been several reports referring the mt genome of in Brassica species, the genome of mt in B. rapa var. purpuraria and its functional gene variations when compared to its closely related species has not yet been addressed.

RESULTS

The mt genome of B. rapa var. purpuraria was assembled through the Illumina and Nanopore sequencing platforms, which revealed a length of 219,775 bp with a typical circular structure. The base composition of the whole B. rapa var. purpuraria mt genome revealed A (27.45%), T (27.31%), C (22.91%), and G (22.32%). 59 functional genes, composing of 33 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes, were annotated. The sequence repeats, codon usage, RNA editing, nucleotide diversity and gene transfer between the cp genome and mt genome were examined in the B. rapa var. purpuraria mt genome. Phylogenetic analysis show that B. rapa var. Purpuraria was closely related to B. rapa subsp. Oleifera and B. juncea. Ka/Ks analysis reflected that most of the PCGs in the B. rapa var. Purpuraria were negatively selected, illustrating that those mt genes were conserved during evolution.

CONCLUSIONS

The results of our findings provide valuable information on the B.rapa var. Purpuraria genome, which might facilitate molecular breeding, genetic variation and evolutionary researches for Brassica species in the future.

Editing of ORF138 restores fertility of Ogura cytoplasmic male sterile broccoli via mitoTALENs

Cytoplasmic male sterility (CMS), encoded by the mitochondrial open reading frames (ORFs), has long been used to economically produce crop hybrids. However, the utilization of CMS also hinders the exploitation of sterility and fertility variation in the absence of a restorer line, which in turn narrows the genetic background and reduces biodiversity. Here, we used a mitochondrial targeted transcription activator-like effector nuclease (mitoTALENs) to knock out ORF138 from the Ogura CMS broccoli hybrid. The knockout was confirmed by the amplification and re-sequencing read mapping to the mitochondrial genome. As a result, knockout of ORF138 restored the fertility of the CMS hybrid, and simultaneously manifested a cold-sensitive male sterility. ORF138 depletion is stably inherited to the next generation, allowing for direct use in the breeding process. In addition, we proposed a highly reliable and cost-effective toolkit to accelerate the life cycle of fertile lines from CMS-derived broccoli hybrids. By applying the k-mean clustering and interaction network analysis, we identified the central gene networks involved in the fertility restoration and cold-sensitive male sterility. Our study enables mitochondrial genome editing via mitoTALENs in Brassicaceae vegetable crops and provides evidence that the sex production machinery and its temperature-responsive ability are regulated by the mitochondria.

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.

Male sterility in plants: an overview of advancements from natural CMS to genetically manipulated systems for hybrid seed production

KEY MESSAGE

The male sterility system in plants has traditionally been utilized for hybrid seed production. In last three decades, genetic manipulation for male sterility has revolutionized this area of research related to hybrid seed production technology. Here, we have surveyed some of the natural cytoplasmic male sterility (CMS) systems that existed/ were developed in different crop plants for developing male sterility-fertility restoration systems used in hybrid seed production and highlighted some of the recent biotechnological advancements in the development of genetically engineered systems that occurred in this area. We have indicated the possible future directions toward the development of engineered male sterility systems. Cytoplasmic male sterility (CMS) is an important trait that is naturally prevalent in many plant species, which has been used in the development of hybrid varieties. This is associated with the use of appropriate genes for fertility restoration provided by the restorer line that restores fertility on the corresponding CMS line. The development of hybrids based on a CMS system has been demonstrated in several different crops. However, there are examples of species, which do not have usable cytoplasmic male sterility and fertility restoration systems (Cytoplasmic Genetic Male Sterility Systems-CGMS) for hybrid variety development. In such plants, it is necessary to develop usable male sterile lines through genetic engineering with the use of heterologous expression of suitable genes that control the development of male gametophyte and fertile male gamete formation. They can also be developed through gene editing using the recently developed CRISPR-Cas technology to knock out suitable genes that are responsible for the development of male gametes. The present review aims at providing an insight into the development of various technologies for successful production of hybrid varieties and is intended to provide only essential information on male sterility systems starting from naturally occurring ones to the genetically engineered systems obtained through different means.

Genetic Mechanisms for Hybrid Breeding in Vegetable Crops

To address the complex challenges faced by our planet such as rapidly changing climate patterns, food and nutritional insecurities, and the escalating world population, the development of hybrid vegetable crops is imperative. Vegetable hybrids could effectively mitigate the above-mentioned fundamental challenges in numerous countries. Utilizing genetic mechanisms to create hybrids not only reduces costs but also holds significant practical implications, particularly in streamlining hybrid seed production. These mechanisms encompass self-incompatibility (SI), male sterility, and gynoecism. The present comprehensive review is primarily focused on the elucidation of fundamental processes associated with floral characteristics, the genetic regulation of floral traits, pollen biology, and development. Specific attention is given to the mechanisms for masculinizing and feminizing cucurbits to facilitate hybrid seed production as well as the hybridization approaches used in the biofortification of vegetable crops. Furthermore, this review provides valuable insights into recent biotechnological advancements and their future utilization for developing the genetic systems of major vegetable crops.

Comparative Transcriptome Analysis Reveals a Potential Regulatory Network for Ogura Cytoplasmic Male Sterility in Cabbage (Brassica oleracea L.)

Ogura cytoplasmic male sterility (CMS) lines are widely used breeding materials in cruciferous crops and play important roles in heterosis utilization; however, the sterility mechanism remains unclear. To investigate the microspore development process and gene expression changes after the introduction of orf138 and Rfo, cytological observation and transcriptome analysis were performed using a maintainer line, an Ogura CMS line, and a restorer line. Semithin sections of microspores at different developmental stages showed that the degradation of tapetal cells began at the tetrad stage in the Ogura CMS line, while it occurred at the bicellular microspore stage to the tricellular microspore stage in the maintainer and restorer lines. Therefore, early degradation of tapetal cells may be the cause of pollen abortion. Transcriptome analysis results showed that a total of 1287 DEGs had consistent expression trends in the maintainer line and restorer line, but were significantly up- or down-regulated in the Ogura CMS line, indicating that they may be closely related to pollen abortion. Functional annotation showed that the 1287 core DEGs included a large number of genes related to pollen development, oxidative phosphorylation, carbohydrate, lipid, and protein metabolism. In addition, further verification elucidated that down-regulated expression of genes related to energy metabolism led to decreased ATP content and excessive ROS accumulation in the anthers of Ogura CMS. Based on these results, we propose a transcriptome-mediated induction and regulatory network for cabbage Ogura CMS. Our research provides new insights into the mechanism of pollen abortion and fertility restoration in Ogura CMS.

Comparative analyses of three complete Primula mitogenomes with insights into mitogenome size variation in Ericales

BACKGROUND

Although knowledge of the sizes, contents, and forms of plant mitochondrial genomes (mitogenomes) is increasing, little is known about the mechanisms underlying their structural diversity. Evolutionary information on the mitogenomes of Primula, an important ornamental taxon, is more limited than the information on their nuclear and plastid counterparts, which has hindered the comprehensive understanding of Primula mitogenomic diversity and evolution. The present study reported and compared three Primula mitogenomes and discussed the size expansion of mitogenomes in Ericales.

RESULTS

Mitogenome master circles were sequenced and successfully assembled for three Primula taxa and were compared with publicly available Ericales mitogenomes. The three mitogenomes contained similar gene contents and varied primarily in their structures. The Primula mitogenomes possessed relatively high nucleotide diversity among all examined plant lineages. In addition, high nucleotide diversity was found among Primula species between the Mediterranean and Himalaya-Hengduan Mountains. Most predicted RNA editing sites appeared in the second amino acid codon, increasing the hydrophobic character of the protein. An early stop in atp6 caused by RNA editing was conserved across all examined Ericales species. The interfamilial relationships within Ericales and interspecific relationships within Primula could be well resolved based on mitochondrial data. Transfer of the two longest mitochondrial plastid sequences (MTPTs) occurred before the divergence of Primula and its close relatives, and multiple independent transfers could also occur in a single MTPT sequence. Foreign sequence [MTPTs and mitochondrial nuclear DNA sequences (NUMTs)] uptake and repeats were to some extent associated with changes in Ericales mitogenome size, although none of these relationships were significant overall.

CONCLUSIONS

The present study revealed relatively conserved gene contents, gene clusters, RNA editing, and MTPTs but considerable structural variation in Primula mitogenomes. Relatively high nucleotide diversity was found in the Primula mitogenomes. In addition, mitogenomic genes, collinear gene clusters, and locally collinear blocks (LCBs) all showed phylogenetic signals. The evolutionary history of MTPTs in Primula was complicated, even in a single MTPT sequence. Various reasons for the size variation observed in Ericales mitogenomes were found.

Complete mitochondrial genome sequencing and identification of candidate genes responsible for C5-type cytoplasmic male sterility in cabbage (B. oleracea var. capitata)

Cytoplasmic male sterility (CMS) is widely used in cruciferous vegetables hybrid breeding. The C5-type CMS cabbage line exhibits stable male sterility and offers great value for cabbage breeding. However, the underlying CMS mechanism remains unclear. Here, the complete mitochondrial genome was sequenced and assembled for this line. The genome size was 221,862 bp. Mitochondrial genome comparison showed that the mitochondrial genome was likely generated by recombination with a nap-type CMS B. napus strain. Sixty-seven unknown-function open reading frames (ORFs) were identified. Seven orfs, orf114a, orf123a, orf188a, orf222a, orf261a, orf286a, and orf322a, were specifically identified in this genome. The presence of these candidate CMS genes decreased ATPase activity and ATP content by affecting the transcript levels of energy metabolism-related genes and F₁F₀-ATP synthase assembly. Among them, orf188a, orf222a, orf261a, orf286a, and orf322a possessed a transmembrane structure, and orf188a was cotranscribed with rps7 and trnfM. orf222a was partially homologous to atp8 and coexpressed with nad5. orf261a and orf322a were cotranscribed with cox1 and atp9, respectively. Additionally, orf114a was cotranscribed with atp8. Yeast two-hybrid assays showed that the ORF222a protein interacts with a B. oleracea ATP17 homolog (Bo7g114140) during F₀-type ATP synthase assembly, reducing the quantity and activity of assembled F₁F₀-ATP synthase. Cytological sections showed that premature separation of the tapetum from the connective tissue and delayed tapetal programmed cell death (PCD) might be the immediate causes of CMS in C5-type CMS cabbage lines. Our results may help uncover the molecular mechanism of C5-type CMS in B. oleracea from the perspectives of the whole mitochondrial genome and cytology of anther development.

Mechanism and Utilization of Ogura Cytoplasmic Male Sterility in Cruciferae Crops

Hybrid production using lines with cytoplasmic male sterility (CMS) has become an important way to utilize heterosis in vegetables. Ogura CMS, with the advantages of complete pollen abortion, ease of transfer and a progeny sterility rate reaching 100%, is widely used in cruciferous crop breeding. The mapping, cloning, mechanism and application of Ogura CMS and fertility restorer genes in Brassica napus, Brassica rapa, Brassica oleracea and other cruciferous crops are reviewed herein, and the existing problems and future research directions in the application of Ogura CMS are discussed.

Organelle Comparative Genome Analysis Reveals Novel Alloplasmic Male Sterility with orf112 in Brassica oleracea L

B. oleracea Ogura CMS is an alloplasmic male-sterile line introduced from radish by interspecific hybridization and protoplast fusion. The introduction of alien cytoplasm resulted in many undesirable traits, which affected the yield of hybrids. Therefore, it is necessary to identify the composition and reduce the content of alien cytoplasm in B. oleracea Ogura CMS. In the present study, we sequenced, assembled, and compared the organelle genomes of Ogura CMS cabbage and its maintainer line. The chloroplast genome of Ogura-type cabbage was completely derived from normal-type cabbage, whereas the mitochondrial genome was recombined from normal-type cabbage and Ogura-type radish. Nine unique regions derived from radish were identified in the mitochondrial genome of Ogura-type cabbage, and the total length of these nine regions was 35,618 bp, accounting for 13.84% of the mitochondrial genome. Using 32 alloplasmic markers designed according to the sequences of these nine regions, one novel sterile source with less alien cytoplasm was discovered among 305 materials and named Bel CMS. The size of the alien cytoplasm in Bel CMS was 21,587 bp, accounting for 8.93% of its mtDNA, which was much less than that in Ogura CMS. Most importantly, the sterility gene orf138 was replaced by orf112, which had a 78-bp deletion, in Bel CMS. Interestingly, Bel CMS cabbage also maintained 100% sterility, although orf112 had 26 fewer amino acids than orf138. Field phenotypic observation showed that Bel CMS was an excellent sterile source with stable 100% sterility and no withered buds at the early flowering stage, which could replace Ogura CMS in cabbage heterosis utilization.

Comparative analysis of the complete mitochondrial genome sequences and anther development cytology between maintainer and Ogura-type cytoplasm male-sterile cabbage (B. oleracea Var. capitata)

Background

Cytoplasmic male sterility (CMS) has been widely used for commercial F1 hybrid seeds production. CMS is primarily caused by chimeric genes in mitochondrial genomes. However, which specific stages of anther development in cabbage are affected by the chimeric genes remain unclear.

Results

In the present study, the complete mitochondrial genomes were sequenced and assembled for the maintainer and Ogura CMS cabbage lines. The genome size of the maintainer and Ogura CMS cabbage are 219,962 bp and 236,648 bp, respectively. There are 67 and 69 unknown function ORFs identified in the maintainer and Ogura CMS cabbage mitochondrial genomes, respectively. Four orfs, orf102a, orf122b, orf138a and orf154a were specifically identified in the Ogura CMS mitochondrial genome, which were likely generated by recombination with Ogura type radish during breeding process. Among them, ORF138a and ORF154a possessed a transmembrane structure, and orf138a was co-transcribed with the atp8 and trnfM genes. orf154a is partially homologous to the ATP synthase subunit 1 (atpA) gene. Both these genes were likely responsible for the CMS phenotype. In addition, cytological sections showed that the abnormal proliferation of tapetal cells might be the immediate cause of cytoplasmic male-sterility in Ogura CMS cabbage lines. RNA-seq results showed that orf138a and orf154a in Ogura CMS might influence transcript levels of genes in energy metabolic pathways.

Conclusions

The presence of orf138a and orf154a lead to increased of ATPase activity and ATP content by affecting the transcript levels of genes in energy metabolic pathways, which could provide more energy for the abnormal proliferation of tapetal cells. Our data provides new insights into cytoplasmic male-sterility from whole mitochondrial genomes, cytology of anther development and transcriptome data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07963-x.

Analysis of mitochondrial recombination in the male sterile Brassica juncea cybrid Og1 and identification of the molecular basis of fertility reversion

Recombinations between the parental genomes produced a novel mitochondrial genome in the cytoplasmic male sterile Brassica juncea cybrid Og1. A mitochondrial stoichiometric shift greatly reduced the molecule containing male-sterility-inducing orf138 gene leading to reversion to male fertility. An improved, chlorosis-corrected, cytoplasmic male sterile Brassica juncea cybrid Og1 derived from Ogura cytoplasm shows frequent reversion to male fertility. To determine the nature of mitochondrial recombination in the cybrid and to uncover the molecular mechanism of male fertility reversion, we sequenced the mitochondrial genomes of Og1, its isonuclear parental lines (OgRLM and Brassica juncea RLM198) and the revertant line (Og1-rt). Assembly of Og1 mitochondrial genome gave two circular molecules, Og1a (250.999 kbp) and Og1b (96.185 kbp) sharing two large direct repeat regions capable of recombining to form a single circular molecule. Og1a contains all essential mitochondrial genes, but the male-sterility-causing orf138 was uniquely present in Og1b along with 16 other complete or partial genes already represented in Og1a. Eleven and four recombinations between the parental mitochondrial genomes produced the Og1a and the Og1b molecules, respectively. Five genes were duplicated within Og1a, of which trnfM was inherited from both the parents while the other four genes, atp4, cox1 nad4L and trnM, were inherited from RLM198. RFLP analysis revealed that orf138-containing molecules were less abundant than Og1a in the male-sterile plants while og1b bearing molecules were undetectable in the revertant line. However, orf138 transcripts were amplified in RT-PCR and were also detected in northern blots revealing that Og1b molecules are not completely lost in the revertant plants. This is the first report where the mitochondrial genome of a cybrid is compared with its actual parents. The findings are discussed in the light of previous reports on mitochondrial genome recombination in cybrids.

Elucidating Mitochondrial DNA Markers of Ogura-Based CMS Lines in Indian Cauliflowers (Brassica oleracea var. botrytis L.) and Their Floral Abnormalities Due to Diversity in Cytonuclear Interactions

Mitochondrial markers can be used to differentiate diverse mitotypes as well as cytoplasms in angiosperms. In cauliflower, cultivation of hybrids is pivotal in remunerative agriculture and cytoplasmic male sterile lines constitute an important component of the hybrid breeding. In diversifying the source of male sterility, it is essential to appropriately differentiate among the available male sterile cytoplasms in cauliflower. PCR polymorphism at the key mitochondrial genes associated with male sterility will be instrumental in analyzing, molecular characterization, and development of mitotype-specific markers for differentiation of different cytoplasmic sources. Presence of auto- and alloplasmic cytonuclear combinations result in complex floral abnormalities. In this context, the present investigation highlighted the utility of organelle genome-based markers in distinguishing cytoplasm types in Indian cauliflowers and unveils the epistatic effects of the cytonuclear interactions influencing floral phenotypes. In PCR-based analysis using a set of primers targeted to orf-138, 76 Indian cauliflower lines depicted the presence of Ogura cytoplasm albeit the amplicons generated exhibited polymorphism within the ofr-138 sequence. The polymorphic fragments were found to be spanning over 200-280 bp and 410-470 bp genomic regions of BnTR4 and orf125, respectively. Sequence analysis revealed that such cytoplasmic genetic variations could be attributed to single nucleotide polymorphisms and insertion or deletions of 31/51 nucleotides. The cytoplasmic effects on varying nuclear-genetic backgrounds rendered an array of floral abnormalities like reduction in flower size, fused flowers, splitted style with the exposed ovule, absence of nonfunctional stamens, and petaloid stamens. These floral malformations caused dysplasia of flower structure affecting female fertility with inefficient nectar production. The finding provides an important reference to ameliorate understanding of mechanism of cytonuclear interactions in floral organ development in Brassicas. The study paves the way for unraveling developmental biology of CMS phenotypes in eukaryotic organisms and intergenomic conflict in plant speciation.

Advances in Research and Application of Male Sterility in Brassica oleracea

Brassica oleracea is an important vegetable species which belongs to the genus Brassica and the mustard family Brassicaceae Burnett. Strong heterosis in B. oleracea is displayed in yield, quality, disease resistance, and stress tolerance. Heterosis breeding is the main way to improve B. oleracea varieties. Male sterile mutants play an important role in the utilization of heterosis and the study of development and regulation in plant reproduction. In this paper, advances in the research and application of male sterility in B. oleracea were reviewed, including aspects of the genetics, cytological characteristics, discovery of genes related to male sterility, and application of male sterility in B. oleracea. Moreover, the main existing problems and prospect of male sterility application in B. oleracea were addressed and a new hybrids’ production strategy with recessive genic male sterility is introduced.

Physical mapping of introgressed chromosome fragment carrying the fertility restoring (Rfo) gene for Ogura CMS in Brassica juncea L. Czern & Coss

Rfo is located on a radish chromosome fragment (~ 108 Kb), which is seated in the middle of a pretty large C genome translocation at the distal region of chromosome A09 of B. juncea. Ogura cytoplasmic male sterility (CMS) is used to produce hybrids in Indian mustard (Brassica juncea L.). Fertility restorers for this CMS were developed by cross-hybridizing B. juncea (AABB; 2n = 36) with B. napus (AACC; 2n = 38) carrying radish Rfo gene. This hybrid production system is normally stable, but many commercial mustard hybrids show male sterile contaminants. We aimed to identify linkage drag associated with Rfo by comparing hybridity levels of 295 handmade CMS x Rfo crosses. Although Rfo was stably inherited, hybridity was < 85 percent in several combinations. Genome re-sequencing of five fertility restorers, mapping sequencing reads to B. juncea reference and synteny analysis with Raphanus sativus D81Rfo genomic region (AJ550021.2) helped to detect ~ 108 Kb of radish chromosome (R) fragment substitution in all fertility restorers. This radish segment substitution was itself located amidst a large C genome translocation on the terminal region of chromosome A09 of B. juncea. The size of alien segment substitution varied from 11.3 (NTCN-R9) to 22.0 Mb (NAJR-102B-R). We also developed an in silico SSR map for chromosome A09 and identified many homoeologous A to the C genome exchanges in the introgressed region. A to the R genome exchanges were rare. Annotation of the substituted fragment showed the gain of many novel genes from R and C genomes and the loss of B. juncea genes from the corresponding region. We have developed a KASPar marker for marker-aided transfer of Rfo and testing hybridity levels in seed production lots.

The dynamic evolution of mobile open reading frames in plastomes of Hymenophyllum Sm. and new insight on Hymenophyllum coreanum Nakai

In this study, four plastomes of Hymenophyllum, distributed in the Korean peninsula, were newly sequenced and phylogenomic analysis was conducted to reveal (1) the evolutionary history of plastomes of early-diverging fern species at the species level, (2) the importance of mobile open reading frames in the genus, and (3) plastome sequence divergence providing support for H. coreanum to be recognized as an independent species distinct from H. polyanthos. In addition, 1C-values of H. polyanthos and H. coreanum were measured to compare the genome size of both species and to confirm the diversification between them. The rrn16-trnV intergenic regions in the genus varied in length caused by Mobile Open Reading Frames in Fern Organelles (MORFFO). We investigated enlarged noncoding regions containing MORFFO throughout the fern plastomes and found that they were strongly associated with tRNA genes or palindromic elements. Sequence identity between plastomes of H. polyanthos and H. coreanum is quite low at 93.35% in the whole sequence and 98.13% even if the variation in trnV-rrn16 intergenic spacer was ignored. In addition, different genome sizes were found for these species based on the 1C-value. Consequently, there is no reason to consider them as a conspecies.

Genome-wide analysis of mRNA and lncRNA expression and mitochondrial genome sequencing provide insights into the mechanisms underlying a novel cytoplasmic male sterility system, BVRC-CMS96, in Brassicarapa

Characterization of a novel and valuable CMS system in Brassicarapa. Cytoplasmic male sterility (CMS) is extensively used to produce F₁ hybrid seeds in a variety of crops. However, it has not been successfully used in Chinese cabbage (Brassicarapa L. ssp. pekinensis) because of degeneration or temperature sensitivity. Here, we characterize a novel CMS system, BVRC-CMS96, which originated in B.napus cybrid obtained from INRAE, France and transferred by us to B.rapa. Floral morphology and agronomic characteristics indicate that BVRC-CMS96 plants are 100% male sterile and show no degeneration in the BC₇ generation, confirming its suitability for commercial use. We also sequenced the BVRC-CMS96 and maintainer line 18BCM mitochondrial genomes. Genomic analyses showed the presence of syntenic blocks and distinct structures between BVRC-CMS96 and 18BCM and the other known CMS systems. We found that BVRC-CMS96 has one orf222 from 'Nap'-type CMS and two copies of orf138 from 'Ogu'-type CMS. We analyzed expression of orf222, orf138, orf261b, and the mitochondrial energy genes (atp6, atp9, and cox1) in flower bud developmental stages S1-S5 and in four floral organs. orf138 and orf222 were both highly expressed in S4, S5-stage buds, calyx, and the stamen. RNA-seq identified differentially expressed mRNAs and lncRNAs (long non-coding RNAs) that were significantly enriched in pollen wall assembly, pollen development, and pollen coat. Our findings suggest that an energy supply disorder caused by orf222/orf138/orf261b may inhibit a series of nuclear pollen development-related genes. Our study shows that BVRC-CMS96 is a valuable CMS system, and our detailed molecular analysis will facilitate its application in Chinese cabbage breeding.

Comparative Analysis of the Mitochondrial Genomes of Nicotiana tabacum: Hints Toward the Key Factors Closely Related to the Cytoplasmic Male Sterility Mechanism

Background

Cytoplasmic male sterility (CMS) is a complex phenomenon of plant sterility that can produce non-functional pollen. It is caused by mutation, rearrangement or recombination in the mitochondrial genome. So far, the systematic structural characteristics of the changes in the mitochondrial genome from the maintainer lines to the CMS lines have not been reported in tobacco.

Results

The mitochondrial genomes of the flower buds from both CMS lines and maintainer lines of two Nicotiana tabacum cultivars (YY85, sYY85, ZY90, and sZY90) were sequenced using the PacBio and Illumina Hiseq technology, and several findings were produced by comparative analysis based on the de novo sequencing. (1) The genomes of the CMS lines were larger, and the different areas were mostly non-coding regions. (2) A large number of rearrangement regions were detected in the CMS lines, with many translocation regions. (3) Thirteen gene clusters were shared by the four mitochondrial genomes, among which two of the gene clusters, nad2-sdh3 and nad6-rps4, were far from each other in the CMS lines. (4) Thirty-three protein-coding genes were conserved in four mitochondrial genomes. However, nad3 was detected one additional copy in the maintainer lines, and sequence differences were revealed in the four candidate genes (atp6, cox2, nad2, and sdh3). Importantly, the evolutionary tree based on the four genes could be used to distinguish the CMS lines and the maintainer lines well for the sequenced mitochondrial genomes of the tobacco. (5) Sixteen CMS-specific open reading frames (ORFs) were found, three of which (orf91, orf115b, and orf100) were previously reported. (6) The differences in intensity of the protein–protein (PPI) interaction in ATP6 were further verified using the yeast two-hybrid analysis.

Conclusion

Although the majority of the sequences, genes and gene clusters were shared by the mitochondrial genomes of the maintainer and the CMS lines in tobacco, extensive structural variations identified with comprehensive analysis based on the mitochondrial genomes, including rearrangement, gene order, the mitochondrial genome expansion and shrinkage events, might be related to CMS. Additionally, the candidate protein-coding genes and CMS-specific ORFs were closely associated with the CMS mechanism. Verification experiments of one of the candidate genes were performed, and the validity of our research results was supported.

Mitochondrial Genome Sequencing Reveals orf463a May Induce Male Sterility in NWB Cytoplasm of Radish

Radish (Raphanus sativus L.) is an important root vegetable worldwide. The development of F1 hybrids, which are extensively used for commercial radish production, relies on cytoplasmic male sterility (CMS). To identify candidate genes responsible for CMS in NWB cytoplasm, we sequenced the normal and NWB CMS radish mitochondrial genomes via next-generation sequencing. A comparative analysis revealed 18 syntenic blocks and 11 unique regions in the NWB CMS mitogenome. A detailed examination indicated that orf463a was the most likely causal factor for male sterility in NWB cytoplasm. Interestingly, orf463a was identical to orf463, which is responsible for CMS in Dongbu cytoplasmic and genic male sterility (DCGMS) radish. Moreover, only structural variations were detected between the NWB CMS and DCGMS mitochondrial genomes, with no nucleotide polymorphisms (SNPs) or meaningful indels. Further analyses revealed these two mitochondrial genomes are coexisting isomeric forms belonging to the same mitotype. orf463a was more highly expressed in flower buds than in vegetative organs and its expression was differentially regulated in the presence of restorer of fertility (Rf) genes. orf463a was confirmed to originate from Raphanus raphanistrum. In this study, we identified a candidate gene responsible for the CMS in NWB cytoplasm and clarified the relationship between NWB CMS and DCGMS.

Organelle genome composition and candidate gene identification for Nsa cytoplasmic male sterility in Brassica napus

Background

Nsa cytoplasmic male sterility (CMS) is a novel alloplasmic male sterility system derived from somatic hybridization between Brassica napus and Sinapis arvensis. Identification of the CMS-associated gene is a prerequisite for a better understanding of the origin and molecular mechanism of this CMS. With the development of genome sequencing technology, organelle genomes of Nsa CMS line and its maintainer line were sequenced by pyro-sequencing technology, and comparative analysis of the organelle genomes was carried out to characterize the organelle genome composition of Nsa CMS as well as to identify the candidate Nsa CMS-associated genes.

Results

Nsa CMS mitochondrial genome showed a higher collinearity with that of S. arvensis than B. napus, indicating that Nsa CMS mitochondrial genome was mainly derived from S. arvensis. However, mitochondrial genome recombination of parental lines was clearly detected. In contrast, the chloroplast genome of Nsa CMS was highly collinear with its B. napus parent, without any evidence of recombination of the two parental chloroplast genomes or integration from S. arvensis. There were 16 open reading frames (ORFs) specifically existed in Nsa CMS mitochondrial genome, which could not be identified in the maintainer line. Among them, three ORFs (orf224, orf309, orf346) possessing chimeric and transmembrane structure are most likely to be the candidate CMS genes. Sequences of all three candidate CMS genes in Nsa CMS line were found to be 100% identical with those from S. arvensis mitochondrial genome. Phylogenetic and homologous analysis showed that all the mitochondrial genes were highly conserved during evolution.

Conclusions

Nsa CMS contains a recombined mitochondrial genome of its two parental species with the majority form S. arvensis. Three candidate Nsa CMS genes were identified and proven to be derived from S. arvensis other than recombination of its two parental species. Further functional study of the candidate genes will help to identify the gene responsible for the CMS and the underlying molecular mechanism.

Mitotypes Based on Structural Variation of Mitochondrial Genomes Imply Relationships With Morphological Phenotypes and Cytoplasmic Male Sterility in Peppers

Plant mitochondrial genomes characteristically contain extensive structural variation that can be used to define and classify cytoplasm types. We developed markers based on structural variation in the mitochondrial genomes of fertile and cytoplasmic male sterility (CMS) pepper lines and applied them to a panel of Capsicum accessions. We designed a total of 20 sequence characterized amplified region (SCAR) markers based on DNA rearrangement junctions or cytoplasm-specific segments that did not show high similarity to any nuclear mitochondrial DNA segments. We used those markers to classify the mitotypes of 96 C. annuum accessions into 15 groups. Precise genotyping of other Capsicum species (C. frutescens, C. chinense, and C. baccatum) was hampered because of various stoichiometric levels of marker amplicons. We developed a multiplex PCR system based on four of the markers that efficiently classified the C. annuum accessions into five mitotype groups. Close relationships between specific mitotypes and morphological phenotypes implied that diversification or domestication of C. annuum germplasm might have been accompanied by structural rearrangements of mitochondrial DNA or the selection of germplasms with specific mitotypes. Meanwhile, CMS lines shared the same amplification profile of markers with another mitotype. Further analysis using mitochondrial DNA (mtDNA) markers based on single-nucleotide polymorphisms (SNPs) or insertions and deletions (InDels) and CMS-specific open reading frames (orfs) provided new information about the origin of the CMS-specific mitotype and evaluation of candidates for CMS-associated genes, respectively.

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.

High and Variable Rates of Repeat-Mediated Mitochondrial Genome Rearrangement in a Genus of Plants

For 30 years, it has been clear that angiosperm mitochondrial genomes evolve rapidly in sequence arrangement (i.e., synteny), yet absolute rates of rearrangement have not been measured in any plant group, nor is it known how much these rates vary. To investigate these issues, we sequenced and reconstructed the rearrangement history of seven mitochondrial genomes in Monsonia (Geraniaceae). We show that rearrangements (occurring mostly as inversions) not only take place at generally high rates in these genomes but also uncover significant variation in rearrangement rates. For example, the hyperactive mitochondrial genome of Monsonia ciliata has accumulated at least 30 rearrangements over the last million years, whereas the branch leading to M. ciliata and its sister species has sustained rearrangement at a rate that is at least ten times lower. Furthermore, our analysis of published data shows that rates of mitochondrial genome rearrangement in seed plants vary by at least 600-fold. We find that sites of rearrangement are highly preferentially located in very close proximity to repeated sequences in Monsonia. This provides strong support for the hypothesis that rearrangement in angiosperm mitochondrial genomes occurs largely through repeat-mediated recombination. Because there is little variation in the amount of repeat sequence among Monsonia genomes, the variable rates of rearrangement in Monsonia probably reflect variable rates of mitochondrial recombination itself. Finally, we show that mitochondrial synonymous substitutions occur in a clock-like manner in Monsonia; rates of mitochondrial substitutions and rearrangements are therefore highly uncoupled in this group.

Mitochondrial genome and transcriptome analysis of five alloplasmic male-sterile lines in Brassica juncea

Background

Alloplasmic lines, in which the nuclear genome is combined with wild cytoplasm, are often characterized by cytoplasmic male sterility (CMS), regardless of whether it was derived from sexual or somatic hybridization with wild relatives. In this study, we sequenced and analyzed the mitochondrial genomes of five such alloplasmic lines in Brassica juncea.

Results

The assembled and annotated mitochondrial genomes of the five alloplasmic lines were found to have virtually identical gene contents. They preserved most of the ancestral mitochondrial segments, and the same candidate male sterility gene (orf108) was found harbored in mitotype-specific sequences. We also detected promiscuous sequences of chloroplast origin that were conserved among plants of the Brassicaceae, and found the RNA editing profiles to vary across the five mitochondrial genomes.

Conclusions

On the basis of our characterization of the genetic nature of five alloplasmic mitochondrial genomes, we speculated that the putative candidate male sterility gene orf108 may not be responsible for the CMS observed in Brassica oxyrrhina and Diplotaxis catholica. Furthermore, we propose the potential coincidence of CMS in alloplasmic lines. Our findings lay the foundation for further elucidation of male sterility gene.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5721-2) contains supplementary material, which is available to authorized users.

Characterization of the mitochondrial genome of the MAX1 type of cytoplasmic male-sterile sunflower

BACKGROUND

More than 70 cytoplasmic male sterility (CMS) types have been identified in Helianthus, but only for less than half of them, research of mitochondrial organization has been conducted. Moreover, complete mitochondrion sequences have only been published for two CMS sources - PET1 and PET2. It has been demonstrated that other sunflower CMS sources like MAX1, significantly differ from the PET1 and PET2 types. However, possible molecular causes for the CMS induction by MAX1 have not yet been proposed. In the present study, we have investigated structural changes in the mitochondrial genome of HA89 (MAX1) CMS sunflower line in comparison to the fertile mitochondrial genome.

RESULTS

Eight significant major reorganization events have been determined in HA89 (MAX1) mtDNA: one 110 kb inverted region, four deletions of 439 bp, 978 bp, 3183 bp and 14,296 bp, respectively, and three insertions of 1999 bp, 5272 bp and 6583 bp. The rearrangements have led to functional changes in the mitochondrial genome of HA89 (MAX1) resulting in the complete elimination of orf777 and the appearance of new ORFs - orf306, orf480, orf645 and orf1287. Aligning the mtDNA of the CMS sources PET1 and PET2 with MAX1 we found some common reorganization features in their mitochondrial genome sequences.

CONCLUSION

The new open reading frame orf1287, representing a chimeric atp6 gene, may play a key role in MAX1 CMS phenotype formation in sunflower, while the contribution of other mitochondrial reorganizations seems to appear negligible for the CMS development.

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.

Evolution of six novel ORFs in the plastome of Mankyua chejuense and phylogeny of eusporangiate ferns

In this paper, three plastomes of Mankyua chejuense, Helminthostachys zeylanica, and Botrychium ternatum in Ophioglossaceae were completely sequenced in order to investigate the plastome evolution and phylogeny of eusporangiate ferns. They were similar to each other in terms of length and the gene orders; however, six unknown open reading frames (ORFs) were found between rps4 and trnL-UAA genes in M. chejuense. Similar sequence regions of six ORFs of M. chejuense were found at the plastomes of Ophioglossum californicum and H. zeylanica, as well as the mitochondrial genome (mitogenome) of H. zeylanica, but not in B. ternatum. Interestingly, the translated amino acid sequences of three ORFs were more similar to the proteins of distantly related taxa such as algae and bacteria than they were to proteins in land plants. It is likely that the six ORFs region arose from endosymbiotic gene transfer (EGT) or horizontal gene transfer (HGT), but further study is needed to verify this. Phylogenetic analyses suggested that Mankyua was resolved as the earliest diverging lineage and that Ophioglossum was subsequently diverged in Ophioglossaceae. This result supports why the plastome of M. chejuense have contained the most ancestral six ORFs in the family.

The comparison of four mitochondrial genomes reveals cytoplasmic male sterility candidate genes in cotton

BACKGROUND

The mitochondrial genomes of higher plants vary remarkably in size, structure and sequence content, as demonstrated by the accumulation and activity of repetitive DNA sequences. Incompatibility between mitochondrial genome and nuclear genome leads to non-functional male reproductive organs and results in cytoplasmic male sterility (CMS). CMS has been used to produce F₁ hybrid seeds in a variety of plant species.

RESULTS

Here we compared the mitochondrial genomes (mitogenomes) of Gossypium hirsutum sterile male lines CMS-2074A and CMS-2074S, as well as their restorer and maintainer lines. First, we noticed the mitogenome organization and sequences were conserved in these lines. Second, we discovered the mitogenomes of 2074A and 2074S underwent large-scale substitutions and rearrangements. Actually, there were five and six unique chimeric open reading frames (ORFs) in 2074A and 2074S, respectively, which were derived from the recombination between unique repetitive sequences and nearby functional genes. Third, we found out four chimeric ORFs that were differentially transcribed in sterile line (2074A) and fertile-restored line.

CONCLUSIONS

These four novel and recombinant ORFs are potential candidates that confer CMS character in 2074A. In addition, our observations suggest that CMS in cotton is associated with the accelerated rates of rearrangement, and that novel expression products are derived from recombinant ORFs.

Organellar genome analysis reveals endosymbiotic gene transfers in tomato

We assembled three complete mitochondrial genomes (mitogenomes), two of Solanum lycopersicum and one of Solanum pennellii, and analyzed their intra- and interspecific variations. The mitogenomes were 423,596-446,257 bp in length. Despite numerous rearrangements between the S. lycopersicum and S. pennellii mitogenomes, over 97% of the mitogenomes were similar to each other. These mitogenomes were compared with plastid and nuclear genomes to investigate genetic material transfers among DNA-containing organelles in tomato. In all mitogenomes, 9,598 bp of plastome sequences were found. Numerous nuclear copies of mitochondrial DNA (NUMTs) and plastid DNA (NUPTs) were observed in the S. lycopersicum and S. pennellii nuclear genomes. Several long organellar DNA fragments were tightly clustered in the nuclear genome; however, the NUMT and NUPT locations differed between the two species. Our results demonstrate the recent occurrence of frequent endosymbiotic gene transfers in tomato genomes.

Integrated analysis of transcriptome and proteome changes related to the Ogura cytoplasmic male sterility in cabbage

Cabbage (Brassica oleracea L. var. capitata), an important vegetable crop in the Brassicaceae family, is economically important worldwide. In the process of hybrid seed production, Ogura cytoplasmic male sterility (OguCMS), controlled by the mitochondrial gene orf138, has been extensively used for cabbage hybrid production with complete and stable male sterility. To identify the critical genes and pathways involved in the sterility and to better understand the underlying molecular mechanisms, the anther of OguCMS line R2P2CMS and the fertile line R2P2 were used for RNA-seq and iTRAQ (Isobaric Tags for Relative and Absolute Quantitation) proteome analysis. RNA-seq analysis generated 13,037,109 to 13,066,594 SE50-clean reads, from the sterile and fertile lines, which were assembled into 36,890 unigenes. Among them, 1,323 differentially expressed genes (DEGs) were identified, consisting of 307 up- and 1016 down-regulated genes. For ITRAQ analysis, a total of 7,147 unique proteins were identified, and 833 were differentially expressed including 538 up- and 295 down-regulated proteins. These were mainly annotated to the ribosome, spliceosome and mRNA surveillance pathways. Combined transcriptomic and proteomic analyses identified 22 and 70 genes with the same and opposite expression profiles, respectively. Using KEGG analysis of DEGs, gibberellin mediated signaling pathways regulating tapetum programmed cell death and four different pathways involved in sporopollenin synthesis were identified. Secretion and translocation of the sporopollenin precursors were identified, and the key genes participating in these pathways were all significantly down-regulated in R2P2CMS. Light and transmission electron (TE) microscopy revealed fat abnormal tapetum rather than vacuolization and degradation at the tetrad and microspore stages of the OguCMS line. This resulted in the failed deposition of sporopollenin on the pollen resulting in sterility. This study provides a comprehensive understanding of the mechanism underlying OguCMS in cabbage.

Recombination Events Involving the atp9 Gene Are Associated with Male Sterility of CMS PET2 in Sunflower

Cytoplasmic male sterility (CMS) systems represent ideal mutants to study the role of mitochondria in pollen development. In sunflower, CMS PET2 also has the potential to become an alternative CMS source for commercial sunflower hybrid breeding. CMS PET2 originates from an interspecific cross of H. petiolaris and H. annuus as CMS PET1, but results in a different CMS mechanism. Southern analyses revealed differences for atp6, atp9 and cob between CMS PET2, CMS PET1 and the male-fertile line HA89. A second identical copy of atp6 was present on an additional CMS PET2-specific fragment. In addition, the atp9 gene was duplicated. However, this duplication was followed by an insertion of 271 bp of unknown origin in the 5' coding region of the atp9 gene in CMS PET2, which led to the creation of two unique open reading frames orf288 and orf231. The first 53 bp of orf288 are identical to the 5' end of atp9. Orf231 consists apart from the first 3 bp, being part of the 271-bp-insertion, of the last 228 bp of atp9. These CMS PET2-specific orfs are co-transcribed. All 11 editing sites of the atp9 gene present in orf231 are fully edited. The anther-specific reduction of the co-transcript in fertility-restored hybrids supports the involvement in male-sterility based on CMS PET2.

Organellar genome assembly methods and comparative analysis of horticultural plants

Although organellar genomes (including chloroplast and mitochondrial genomes) are smaller than nuclear genomes in size and gene number, organellar genomes are very important for the investigation of plant evolution and molecular ecology mechanisms. Few studies have focused on the organellar genomes of horticultural plants. Approximately 1193 chloroplast genomes and 199 mitochondrial genomes of land plants are available in the National Center for Biotechnology Information (NCBI), of which only 39 are from horticultural plants. In this paper, we report an innovative and efficient method for high-quality horticultural organellar genome assembly from next-generation sequencing (NGS) data. Sequencing reads were first assembled by Newbler, Amos, and Minimus software with default parameters. The remaining gaps were then filled through BLASTN search and PCR. The complete DNA sequence was corrected based on Illumina sequencing data using BWA (Burrows-Wheeler Alignment tool) software. The advantage of this approach is that there is no need to isolate organellar DNA from total DNA during sample preparation. Using this procedure, the complete mitochondrial and chloroplast genomes of an ornamental plant, Salix suchowensis, and a fruit tree, Ziziphus jujuba, were identified. This study shows that horticultural plants have similar mitochondrial and chloroplast sequence organization to other seed plants. Most horticultural plants demonstrate a slight bias toward A+T rich features in the mitochondrial genome. In addition, a phylogenetic analysis of 39 horticultural plants based on 15 protein-coding genes showed that some mitochondrial genes are horizontally transferred from chloroplast DNA. Our study will provide an important reference for organellar genome assembly in other horticultural plants. Furthermore, phylogenetic analysis of the organellar genomes of horticultural plants could accurately clarify the unanticipated relationships among these plants.

Rapid evolutionary divergence of diploid and allotetraploid Gossypium mitochondrial genomes

Background

Cotton (Gossypium spp.) is commonly grouped into eight diploid genomic groups and an allotetraploid genomic group, AD. The mitochondrial genomes supply new information to understand both the evolution process and the mechanism of cytoplasmic male sterility. Based on previously released mitochondrial genomes of G. hirsutum (AD₁), G. barbadense (AD₂), G. raimondii (D₅) and G. arboreum (A₂), together with data of six other mitochondrial genomes, to elucidate the evolution and diversity of mitochondrial genomes within Gossypium.

Results

Six Gossypium mitochondrial genomes, including three diploid species from D and three allotetraploid species from AD genome groups (G. thurberi D₁, G. davidsonii D_3-d and G. trilobum D₈; G. tomentosum AD₃, G. mustelinum AD₄ and G. darwinii AD₅), were assembled as the single circular molecules of lengths about 644 kb in diploid species and 677 kb in allotetraploid species, respectively. The genomic structures of mitochondrial in D group species were identical but differed from the mitogenome of G. arboreum (A₂), as well as from the mitogenomes of five species of the AD group. There mainly existed four or six large repeats in the mitogenomes of the A + AD or D group species, respectively. These variations in repeat sequences caused the major inversions and translocations within the mitochondrial genome. The mitochondrial genome complexity in Gossypium presented eight unique segments in D group species, three specific fragments in A + AD group species and a large segment (more than 11 kb) in diploid species. These insertions or deletions were most probably generated from crossovers between repetitive or homologous regions. Unlike the highly variable genome structure, evolutionary distance of mitochondrial genes was 1/6th the frequency of that in chloroplast genes of Gossypium. RNA editing events were conserved in cotton mitochondrial genes. We confirmed two near full length of the integration of the mitochondrial genome into chromosome 1 of G. raimondii and chromosome A03 of G. hirsutum, respectively, with insertion time less than 1.03 MYA.

Conclusion

Ten Gossypium mitochondrial sequences highlight the insights to the evolution of cotton mitogenomes.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4282-5) contains supplementary material, which is available to authorized users.

Complete mitochondrial genome sequences of Brassica rapa (Chinese cabbage and mizuna), and intraspecific differentiation of cytoplasm in B. rapa and Brassica juncea

The complete sequence of the mitochondrial genome was determined for two cultivars of Brassica rapa. After determining the sequence of a Chinese cabbage variety, 'Oushou hakusai', the sequence of a mizuna variety, 'Chusei shiroguki sensuji kyomizuna', was mapped against the sequence of Chinese cabbage. The precise sequences where the two varieties demonstrated variation were ascertained by direct sequencing. It was found that the mitochondrial genomes of the two varieties are identical over 219,775 bp, with a single nucleotide polymorphism (SNP) between the genomes. Because B. rapa is the maternal species of an amphidiploid crop species, Brassica juncea, the distribution of the SNP was observed both in B. rapa and B. juncea. While the mizuna type SNP was restricted mainly to cultivars of mizuna (japonica group) in B. rapa, the mizuna type was widely distributed in B. juncea. The finding that the two Brassica species have these SNP types in common suggests that the nucleotide substitution occurred in wild B. rapa before both mitotypes were domesticated. It was further inferred that the interspecific hybridization between B. rapa and B. nigra took place twice and resulted in the two mitotypes of cultivated B. juncea.

Identification of different cytoplasms based on newly developed mitotype-specific markers for marker-assisted selection breeding in Brassica napus L

Different mitotype-specific markers were developed to distinguish different cytoplasms in Brassica napus L. Mitotype-specific markers have been developed to distinguish different mitotypes in plant. And use of molecular markers to identify different mitotypes in Brassica napus would enhance breeding efficiency. Here, we comparatively analyzed six sequenced mitochondrial genomes in Brassica napus and identified collinear block sequences and mitotype-specific sequences (MSSs) of these mitochondrial genomes. The collinear block sequences between mitochondrial genomes of nap, cam, and pol cytoplasmic male sterility (CMS) lines were higher than those of other lines. After comparative analysis of the six sequenced mitochondrial genomes (cam, nap, ole, pol CMS, ogu CMS, and hau CMS), 90 MSSs with sizes ranging from 101 to 9981 bp and a total length of 103,756 bp (accounting for 6.77% of the mitochondrial genome sequences) were identified. Additionally, 12 mitotype-specific markers were developed based on the mitochondrial genome-specific sequences in order to distinguish among these different mitotypes. Cytoplasms of 570 different inbred lines collected across scientific research institutes in China were identified using the MSS markers developed in our study. In addition to confirming the accuracy of the cytoplasmic identification, we also identified mitotypes that have not been reported in Brassica napus. Our study may provide guidance for the classification of different mitotypes in B. napus breeding.

Entire nucleotide sequences of Gossypium raimondii and G. arboreum mitochondrial genomes revealed A-genome species as cytoplasmic donor of the allotetraploid species

Cotton (Gossypium spp.) is commonly grouped into eight diploid genomic groups, designated A-G and K, and an allotetraploid genomic group, AD. Gossypium raimondii (D₅ ) and G. arboreum (A₂ ) are the putative contributors to the progenitor of G. hirsutum (AD₁ ), the economically important fibre-producing cotton species. Mitochondrial DNA from week-old etiolated seedlings was extracted from isolated organelles using discontinuous sucrose density gradient method. Mitochondrial genomes were sequenced, assembled, annotated and analysed in orderly. Gossypium raimondii (D₅ ) and G. arboreum (A₂ ) mitochondrial genomes were provided in this study. The mitochondrial genomes of two diploid species harboured circular genome of 643,914 bp (D₅ ) and 687,482 bp (A₂ ), respectively. They differ in size and number of repeat sequences, both contain illuminating triplicate sequences with 7317 and 10,246 bp, respectively, demonstrating dynamic difference and rearranged genome organisations. Comparing the D₅ and A₂ mitogenomes with mitogenomes of tetraploid Gossypium species (AD₁ , G. hirsutum; AD₂ , G. barbadense), a shared 11 kbp fragment loss was detected in allotetraploid species, three regions shared by G. arboreum (A₂ ), G. hirsutum (AD₁ ) and G. barbadense (AD₂ ), while eight regions were specific to G. raimondii (D₅ ). The presence/absence variations and gene-based phylogeny supported that A-genome is a cytoplasmic donor to the progenitor of allotetraploid species G. hirsutum and G. barbadense. The results present structure variations and phylogeny of Gossypium mitochondrial genome evolution.

A Novel Cytoplasmic Male Sterility in Brassica napus (inap CMS) with Carpelloid Stamens via Protoplast Fusion with Chinese Woad

A novel cytoplasmic male sterility (CMS) in Brassica napus (inap CMS) was selected from the somatic hybrid with Isatis indigotica (Chinese woad) by recurrent backcrossing. The male sterility was caused by the conversion of tetradynamous stamens into carpelloid structures with stigmatoid tissues at their tips and ovule-like tissues in the margins, and the two shorter stamens into filaments without anthers. The feminized development of the stamens resulted in the complete lack of pollen grains, which was stable in different years and environments. The pistils of inap CMS displayed normal morphology and good seed-set after pollinated by B. napus. Histological sections showed that the developmental alteration of the stamens initiated at the stage of stamen primordium differentiation. AFLP analysis of the nuclear genomic composition with 23 pairs of selective primers detected no woad DNA bands in inap CMS. Twenty out of 25 mitochondrial genes originated from I. indigotica, except for cox2-2 which was the recombinant between cox2 from woad and cox2-2 from rapeseed. The novel cox2-2 was transcribed in flower buds of inap CMS weakly and comparatively with the fertile B. napus addition line Me harboring one particular woad chromosome. The restorers of other autoplasmic and alloplasmic CMS systems in rapeseed failed to restore the fertility of inap CMS and the screening of B. napus wide resources found no fertility restoration variety, showing its distinct origin and the related mechanism of sterility. The reasons for the mitochondrial rearrangements and the breeding of the restorer for the novel CMS system were discussed.

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.

Development of a novel allele-specific Rfo marker and creation of Ogura CMS fertility-restored interspecific hybrids in Brassica oleracea

A novel allele-specific Rfo marker was developed and proved to be effective for MAS of Rfo gene in B. oleracea background and six Ogu-CMS fertility-restored interspecific hybrids were created for the first time. Ogura cytoplasmic male sterility (Ogu-CMS) has been extensively used for Brassica oleracea hybrid production. However, because of maternal inheritance, all the hybrids produced by CMS lines are male sterile and cannot be self-pollinated, which prohibits germplasm maintenance and innovation. This problem can be overcome by using the Ogu-CMS restorer line, but restorer material is absent in B. oleracea crops. Here, Rfo, a fertility-restored gene of Ogu-CMS, was transferred from rapeseed restorer lines into a Chinese kale Ogu-CMS line using interspecific hybridization combined with embryo rescue. Nine interspecific, triploid plant progenies were identified at morphological and ploidy level, with phenotypes intermediate between those of rapeseed and Chinese kale. Because the Rfo marker (Hu et al., Mol Breeding 22:663-674, 2008) cannot distinguish the Rfo and its homologies under a B. oleracea background, a novel allele-specific Rfo marker was developed based on the BLAST analysis of highly homologous Rfo sequences in B. oleracea. Screening using the novel Rfo marker found that six interspecific hybrids carrying Rfo were also fertile, although fertility varied during different flowering periods. Furthermore, BC1 offsprings with the Rfo gene were selected with the allele-specific Rfo marker and showed restored fertility. These results indicated that the novel allele-specific marker could be used for the MAS of Rfo gene in B. oleracea, and this study lays the foundation for the development of Ogu-CMS restorer material in cabbage and its related other subspecies.

The Whole Genome Assembly and Comparative Genomic Research of Thellungiella parvula (Extremophile Crucifer) Mitochondrion

The complete nucleotide sequences of the mitochondrial (mt) genome of an extremophile species Thellungiella parvula (T. parvula) have been determined with the lengths of 255,773 bp. T. parvula mt genome is a circular sequence and contains 32 protein-coding genes, 19 tRNA genes, and three ribosomal RNA genes with a 11.5% coding sequence. The base composition of 27.5% A, 27.5% T, 22.7% C, and 22.3% G in descending order shows a slight bias of 55% AT. Fifty-three repeats were identified in the mitochondrial genome of T. parvula, including 24 direct repeats, 28 tandem repeats (TRs), and one palindromic repeat. Furthermore, a total of 199 perfect microsatellites have been mined with a high A/T content (83.1%) through simple sequence repeat (SSR) analysis and they were distributed unevenly within this mitochondrial genome. We also analyzed other plant mitochondrial genomes' evolution in general, providing clues for the understanding of the evolution of organelles genomes in plants. Comparing with other Brassicaceae species, T. parvula is related to Arabidopsis thaliana whose characters of low temperature resistance have been well documented. This study will provide important genetic tools for other Brassicaceae species research and improve yields of economically important plants.

Comparative Transcriptome Profile of the Cytoplasmic Male Sterile and Fertile Floral Buds of Radish (Raphanus sativus L.)

Radish cytoplasmic male sterility (CMS) has been widely used for breeding in Raphanus and Brassica genera. However, the detailed regulation network of the male sterility remains to be determined. Our previous work has shown that the abnormalities in a CMS radish appeared shortly after the tetrad stage when microspores were malformed and the tapetal cells grew abnormally large. In this work, histological analysis shows that anthers are at the tetrad stage when the radish buds are about 1.5 mm in length. Furthermore, a high throughput RNA sequencing technology was employed to characterize the transcriptome of radish buds with length about 1.5 mm from two CMS lines possessing the CMS-inducing orf138 gene and corresponding near-isogenic maintainer lines. A total of 67,140 unigenes were functionally annotated. Functional terms for these genes are significantly enriched in 55 Gene Ontology (GO) groups and 323 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The transcriptome detected transcripts for 72 out of a total of 79 protein genes encoded in the chloroplast genome from radish. In contrast, the radish mitochondrial genome contains 34 protein genes, but only 16 protein transcripts were detected from the transcriptome. The transcriptome comparison between CMS and near-isogenic maintainer lines revealed 539 differentially expressed genes (DEGs), indicating that the false positive rate for comparative transcriptome profiling was clearly decreased using two groups of CMS/maintainer lines with different nuclear background. The level of 127 transcripts was increased and 412 transcripts were decreased in the CMS lines. No change in levels of transcripts except CMS-inducing orf138 was identified from the mitochondrial and chloroplast genomes. Some DEGs which would be associated with the CMS, encoding MYB and bHLH transcription factors, pentatricopeptide repeat (PPR) proteins, heat shock transcription factors (HSFs) and heat shock proteins (HSPs), are discussed. The transcriptome dataset and comparative analysis will provide an important resource for further understanding anther development, the CMS mechanism and to improve molecular breeding in radish.