Skip to main content

Comparative analysis of mitochondrial genomes of maize CMS-S subtypes provides new insights into male sterility stability



Cytoplasmic male sterility (CMS) is a trait of economic importance in the production of hybrid seeds. In CMS-S maize, exerted anthers appear frequently in florets of field-grown female populations where only complete male-sterile plants were expected. It has been reported that these reversions are associated with the loss of sterility-conferring regions or other rearrangements in the mitochondrial genome. However, the relationship between mitochondrial function and sterility stability is largely unknown.


In this study, we determined the ratio of plants carrying exerted anthers in the population of two CMS-S subtypes. The subtype with a high ratio of exerted anthers was designated as CMS-Sa, and the other with low ratio was designated as CMS-Sb. Through next-generation sequencing, we assembled and compared mitochondrial genomes of two CMS-S subtypes. Phylogenetic analyses revealed strong similarities between the two mitochondrial genomes. The sterility-associated regions, S plasmids, and terminal inverted repeats (TIRs) were intact in both genomes. The two subtypes maintained high transcript levels of the sterility gene orf355 in anther tissue. Most of the functional genes/proteins were identical at the nucleotide sequence and amino acid sequence levels in the two subtypes, except for NADH dehydrogenase subunit 1 (nad1). In the mitochondrial genome of CMS-Sb, a 3.3-kilobase sequence containing nad1-exon1 was absent from the second copy of the 17-kb repeat region. Consequently, we detected two copies of nad1-exon1 in CMS-Sa, but only one copy in CMS-Sb. During pollen development, nad1 transcription and mitochondrial biogenesis were induced in anthers of CMS-Sa, but not in those of CMS-Sb. We suggest that the impaired mitochondrial function in the anthers of CMS-Sb is associated with its more stable sterility.


Comprehensive analyses revealed diversity in terms of the copy number of the mitochondrial gene nad1-exon1 between two subtypes of CMS-S maize. This difference in copy number affected the transcript levels of nad1 and mitochondrial biogenesis in anther tissue, and affected the reversion rate of CMS-S maize. The results of this study suggest the involvement of mitochondrial robustness in modulation of sterility stability in CMS-S maize.

Peer Review reports


Cytoplasmic male sterility (CMS) is the maternally inherited inability of plants to produce viable pollen, while plant growth and female fertility are normal [1, 2]. In maize (Zea mays L.), there are three major groups of CMS: CMS-C, CMS-T, and CMS-S, which are defined according to the sterilizing factors in mitochondrion and the corresponding main restorers in the nucleus. Maize CMS is an economically valuable trait used for the production of hybrid seeds. The CMS-T system was first used with three-line system technology to produce hybrid maize, which accounted for about 85% of hybrid seeds in the U.S.A. until the 1970 southern corn leaf blight epidemic [3, 4]. Because only one single major restorer (Rf3) is required for pollen fertility in the S-system, some breeders transitioned to use CMS-S to produce hybrid maize.

CMS-S maize has a clearly elucidated mitochondria-nucleus interaction system, wherein the sterilizing factor orf355 initiates microspore degeneration while the single main restorer Rf3 cleaves orf355 transcripts to restore fertility. CMS-S mitochondria have two linear plasmids known as S1 and S2, which have exactly the same terminal inverted repeats (TIRs) as that in the CMS-associated region. These plasmids actively recombine with the circular genome to linearize the mitochondrial DNA (mtDNA) [5,6,7,8]. The gene encoding the sterilizing factor orf355 is located near the end of the linearized mtDNA and is expressed at the bi-cellular stage of microspore development, leading to the gametophytic sterility of maize [9, 10]. In the presence of Rf3, the 1.6-kb transcripts containing orf355 are cleaved via posttranscriptional modification, leading to reduced orf355 transcript levels and the restoration of fertility [11].

However, CMS-S maize is a relatively unstable system where revertants frequently arise as a result of genetic mutation [12]. In the absence of Rf3, sterility reversions are caused by less effective genetic restorer-of-fertility genes, or by the mtDNA changes that usually arise in de novo [8, 9, 13, 14]. Although the reversions usually take place in a single floret or a small sector of the tassel, they can potentially reduce the purity of hybrid seeds. This can hinder the use of CMS-S for hybrid seed production. Unstable sterility usually occurs in the absence of the main restorer Rf3. For example, nuclear-encoded restorer Rf9 does not cause the cleavage of the 1.6-kb transcripts but decreases the abundance of the linearized transcripts of the mtDNA template [8]. In addition, the loss of the free S1 and S2 plasmids from the mitochondria blocks the rearrangement with the sterility-associated region, leading to the restoration of fertility in CMS-S maize [9]. In some cases, mtDNA rearrangement disrupts the sterility-associated region so that the sterilizing factor is either lost or not expressed. For example, some revertants have a 7.3-kb inversion in the mtDNA that separates the TIR sequences from the CMS-associated region [13]. As a result, orf355 is no longer transcribed from mtDNA linear ends but is co-transcribed with cox2. Although orf355 transcripts can still be detected in these plants, they are not highly expressed and no longer confer sterility [13]. The restorer-of-fertility lethal 1 (rfl1) mutant disrupts mitochondrial gene expression and the accumulation of the α-subunit of ATP synthase, suggesting that the functional plasticity of mitochondria is linked with the male sterility stability of CMS-S [14].

The maize mitochondrial genome is about 500 kb in size, and contains a suite of relatively conserved protein-encoding genes within different cytotypes [15]. The relative placement of the genes and the intergenic spacer regions within the mitochondrial genome vary extensively among different maize subgroups [16]. A major reason for the highly variable structural organization of the maize mtDNA is the abundance of recombination-active repeated sequences [17]. The shuffling of mtDNA sequences by recombination plays a role in evolution, and changes gene organization and creates gene chimeras [16]. All three types of CMS in maize, as well as CMS systems in other species, are caused by chimeric open reading frames (ORFs) arising from mtDNA rearrangements. Mitochondrial rearrangements associated with the loss of portions of essential genes cause poor growth or lethality in maize [18, 19]. In contrast, rearrangements affecting only non-coding regions of mtDNA are usually neutral, but sometimes cause different phenotypes, such as fertility reversion in CMS-S maize [20,21,22]. Maize mitochondrial genomes have a special type of DNA sequence acquired from the plastid, which accounts for 4.4% of CMS-S mitochondrial genome [23]. The plastid-derived DNA sequences vary widely in their content and location among the five maize mitochondrial genomes [15]. However, the relationship between the plastid-derived DNA sequences and recombination is unclear.

The mitochondrion is a semi-autonomous organelle, which is capable of responding to fluctuating energy demands and environmental stimuli [24]. Mitochondrion is capable of maintaining a constant ATP to ADP ratio over a wide range of conditions. This key property is referred to as mitochondrial robustness, and can explain why some even severe mutations can be tolerated within the oxidative phosphorylation system and be compatible with life. Studies on animals suggest that the tissue-specific control coefficients of different respiratory chain complexes counteract the influence of mitochondria, and probably vice versa [25]. Plant CMS might result from the inability of mitochondria to meet the increased energy demands of microspore development [2]. However, the relationship between mitochondrial robustness and sterility determination is largely unknown. Here we describe features of the mitochondrial genome of a CMS-S subtype that has retained an intact orf355/orf77 region but has lost the second copy of nad1-exon1. This change was found to reduce mitochondrial biogenesis and gene transcription in anther tissue, leading to stable male sterility in CMS-S maize.

Materials and methods

Plant materials and growth conditions

The subtype CMS-Sa was originated from male sterile line S-Mo17rf3rf3 [26] and were donated by professor Yonglian Zheng in Huazhong Agricultural University (Wuhan, China). The subtype CMS-Sb were originated from male fertile line S-SD13Rf3Rf3 [10] and were donated by professor Huabang Chen in Institute of Genetics and Developmental Biology, Chinese Academy of Sciences (Beijing, China). The sterility of both lines can be counteracted by the Rf3 restorer. Jing72464 is the female parent of an elite hybrid line in China, which was developed by the Maize Research Institute, Beijing Academy of Agricultural and Forestry Science (Beijing, China). Both subtypes were introduced into Jing72464 to develop the sterile inbred line. By crossing two cytoplasm donors with the fertile Jing72464 (NB), the sterile F1 was obtained, which was consecutively backcrossed with Jing72464 for more than 10 generations. Molecular markers were used to select the Jing72464 nuclear background in each backcross generation, so that the sterile families in the two CMS-S subtypes can be considered as near isogenic lines. The plant materials were cultivated at the Hainan maize propagation base in Yacheng, Hainan (HN-YC, 18.3°N, 109.5°E) and a farm near Beijing (40.1° N, 116.4° E).

Pollen viability test

A fertility survey was conducted at the Beijing farm (40.1° N, 116.4° E) in the summer of 2021. During flowering, every plant was observed every day to detect exerted anthers. Plants with even one exerted anther in any sector of the tassels were counted. Pollen viability was checked using the iodine-potassium iodide (I2-KI) staining method. The pre-exerted anthers were cut with tweezers in 50 μL I2-KI solution on a glass slide, and then large debris were carefully removed. The stained pollen grains without a cover slip were directly observed under a light microscope at 20× magnification as described previously [27]. At least 100 pollen grains per visual field were counted randomly under a light microscope. Successful seed set indicated that pollen was fertile. Pollen grains were carefully collected from the exerted anthers and crossed onto ears of CMS-Sb Jing72464 inbred plants to confirm pollen fertility. In some cases, anthers were cut with tweezers and the pollen was applied directly to filaments. Progeny families from the above tests were cultivated and their fertility was checked to determine whether the restoration of pollen fertility could be transmitted via the pollen.

Library construction and sequencing

Mitochondrial DNA was extracted as described previously [28]. Ten-day-old etiolated seedlings were homogenized in homogenization buffer (0.4 M mannitol, 10 mM TES pH 7.5, 5 mM EGTA, 0.05% cysteine, 0.2% BSA) in a blender, then the mixture was filtered through four layers of cheesecloth and one layer of Miracloth. Crude mitochondria were obtained by differential centrifugation (1000 g, 5 min; 2000 g, 10 min; 10,000 g, 20 min). The crude mitochondrial pellet was suspended in 2 mL homogenization buffer, then purified using sucrose step gradient centrifugation. The sucrose step gradient was established using sterile sucrose solutions (60, 47, 35, and 20% w/v sucrose in 10 mM tricine, 1 mM EGTA, pH 7.5) in ultraclear tubes. The resuspended mitochondria were carefully layered onto the sucrose gradient and centrifuged at 30,000 rpm (about 111,000 g) in a SW41 rotor for 60 min at 4 °C. The mitochondria were condensed at the interface between the 35 and 47% sucrose layers. The collected mitochondria were diluted with three volumes of homogenization buffer. Purified mitochondria were obtained by centrifuging the mixture at 10,000 g for 30 min at 4 °C. Then, DNA was extracted from the purified mitochondria using a DNeasy plant mini kit (QIAGEN, Hilden, Germany).

The libraries were constructed with the TruSeq Nano DNA LT Sample Preparation Kit (Illumina, San Diego, CA, USA). Briefly, the mitochondrial DNA was sheared into fragments with lengths of approximately 350 bp using a S220 Focused ultrasonicator (Covaris, San Diego, CA, USA). Adapters were ligated onto the 3′ end of the sheared fragments. After PCR amplification and purification, the final libraries were sequenced on the Illumina sequencing HiSeq X Ten platform (Illumina) and 150-bp paired-end reads were generated. The sequence data were de novo assembled using ABySS [29] with different k-mer values (k = 31 to 127, increasing step-wise by 2 each time) to choose the best k-mer value for assembly. To verify the quality and accuracy of our assemblies, the original reads were mapped to the corresponding mitogenomes using BWA [30]. Pilon was used to polish the final assembly. Mitochondrial protein-coding genes were predicted using GeSeq [31], and then the start/stop codons and the exon-intron boundaries of genes were revised manually. The GC content was determined with a shell script and the circular physical map of all mitogenomes was visualized using OGDRAW [32]. The Synteny and single nucleotide polymorphism (SNP) analyses of the different mitogenomes were conducted using MUMmer [33].

Phylogenetic analysis

To understand the phylogenetic position of two CMS-S subtypes, the mitochondrial genomes of CMS-Sa, CMS-Sb, and DQ490951 were used to build an evolutionary tree, mitochondrial genome of Sorghum (DQ984518) was used as the outgroup. Phylogenetic analysis was based on nucleotide sequences of the 4 mitochondrial genomes. These nucleotides were aligned using MAFFT [34] and gaps were removed from the alignments using the trimAL [35]. A maximum likelihood (ML) tree was constructed using IQ-TREE [36] with a Best-fit model from ModelFinder [37] and bootstrap consensus was inferred from 1000 replications. The tree was finally drawn using R script.

Analysis of mitochondrial plastid DNAs and the corresponding plastid homolog

A total of 11 pairs mitochondrial plastid DNA and plastid homologs were selected and compared based on previous study [15]. Fragments covering the target sequences were amplified from genomic DNA isolated from ten-days-old seedlings using suites of specific primers. Then the fragments were sequenced by Sanger sequencing. The primers for the amplification of these fragment are listed in Additional file 1: Table S1.

Quantification of mitochondrial DNA

Mitochondrial DNA in maize seedlings was quantified as described previously [38] with slight modifications. Total DNA was extracted from a pool of five ten-day-old maize seedlings leaves, or from the anthers collected from pre-emerged tassels. A total of 40 ng DNA was analyzed by real-time PCR using a Light Cycler 480 (Roche, Basel, Switzerland). The reaction mixture (final volume, 20 μL) contained SYBR green master mix (Roche) and 1.25 mM forward and reverse primers. The proportion of mitochondrial DNA out of total DNA was determined using the primer pairs C33–C34 and C35–C36, then normalized to the internal standard nucleus-encoded single-copy fragment C39–C40 as described previously [38]. Each data point represents the average of at least three technical replicates. Primer sequences are shown in Additional file 1: Table S1.

RNA extraction and gene transcript level analysis

Whole seedlings with the first expanded true leaf were sampled from each subtype. Anther tissue was carefully collected from the florets of pre-emerged tassels. Total RNA was isolated using a Plant RNA Purification Kit (K0801, Thermo Fisher Scientific, Waltham, MA, USA). Full-length cDNA was generated from 2 μg total RNA per sample using random primers (cDNA synthesis kit, K1622, Thermo Fisher Scientific). The resulting cDNA was diluted 1:10 for quantitative RT-PCR (qRT-PCR) analyses, which were conducted using the TB Green® Premix Ex Taq™ II Real-Time PCR System (RR820A, Takara, Otsu, Japan). Each 20 μL PCR mixture consisted of 10 μL premix, 3 μL template cDNA, and 0.4 μM each PCR primer. The thermal cycling conditions were as follows: 95 °C for 2 min, followed by 40 cycles at 95 °C for 20 s, 56 °C for 1 min, and 72 °C for 30 s. Samples were run on the LightCycler 480 Real-Time PCR System (Roche). Melting curve data were used in every experiment to verify amplification of the appropriate PCR product. The mitochondrion-encoded gene rrn18 was used as the internal control to monitor sample uniformity of initial RNA input and reverse transcription efficiency. The 2−ΔΔCt method was used to calculate the relative transcript levels of mitochondrial genes. The primer sequences are listed in Additional file 1: Table S1.

Circular RT-PCR analyses

Circular RT-PCR (cRT-PCR) analyses were conducted to detect the primary and processed transcripts of nad1 and nad2 from the 17-kb repeat region. The extracted RNAs were treated with RNA 5′ polyphosphatase (#M0356, NEB, Boston, MA, USA) before self-ligation. The circular cDNA was reverse-transcribed using gene-specific primers. Two-round PCR was conducted to amplify the sequences flanking the 5′ and 3′ extremities. The products of second-round PCR were separated by agarose gel electrophoresis, and the bands of expected size were sequenced by cloning into a vector. These procedures were conducted as described elsewhere [39] and according to the manufacturer’s instructions. The primers used for reverse transcription and cRT-PCR amplifications are listed in Additional file 1: Table S1.


Different sterility stability in two subtypes of CMS-S maize

In the hybrid breeding program, two subtypes of CMS-S (designated as CMS-Sa and CMS-Sb) were used as the cytoplasm donors to generate sterile female parents of elite maize hybrids (see methods). Of the 319 plants derived from subtype CMS-Sa, 188 had visible anthers exerted from the florets. These anthers were smaller and usually emerged from the basal region of the tassel sectors after natural pollination (Fig. 1 a-d). The male sterile stability was significantly higher in CMS-Sb than in CMS-Sa (only 4 of the 270 progenies produced visible anthers at the same stage). The anthers in subtype Sb were severely shrunken, implying that there was no viable pollen in the locules (Fig. 1a-d). Pollen development was checked via the I2-KI staining method. Among the 178 analyzed anthers from CMS-Sa, most (170/178) produced a small amount of stainable pollen (Fig. 1e). In contrast, none of the four analyzed anthers from CMS-Sb produced stainable pollen. To test the pollen viability of CMS-Sa, revertants were crossed as pollen parents onto ears of pollen-sterile testers. Twenty-three families were produced with successful seed set, indicating that the pollen was indeed viable (Fig. 1f). However, the seed setting rate was very low, which less than 10% in the representative ears. Notably, natural pollination of these plants had already finished by the time the anthers were exerted from the florets, implying that fertility reversion was only partial and the leaked pollen grains were able to meet the unfertilized silks to complete fertilization only rarely. We speculated that, because the experimental materials were near isogenic lines, variation in fertility reversion was likely due to non-nuclear differences. To test this idea, 15 of the above 23 families with good seed set were grown and their pollen fertility was evaluated. As expected, none of the plants were fertile or partially fertile, indicating that no nuclear restorer allele had been transmitted from the pollen to the next generation.

Fig. 1
figure 1

Different sterility stability in two subtypes of CMS-S maize. a Pictures of tassel sectors on cytoplasm N, CMS-Sa, and CMS-Sb. Scale bar, 1 cm. The lower left panel show full image of each tassel. Scale bar, 10 cm. The lower right panel show representative pollen grains that stained with I-Ki solution in each subtypes. b Spontaneous reversion rate of CMS-S subtype Sa and Sb, plant with exerted anther in floret was considered as spontaneous reverted individuals. c Image of anthers in CMS-Sa and CMS-Sb. Scale bar, 1 cm. d Anther size calculated using Image J. Data are mean ± s.e.m. (n = 6). e Fertile pollen ratios, calculated as the ratio of darkly stained pollens versus the lightly stained pollens in each plantlets. f Tests of pollen viability, as illustrated by seed setting rate in representative ears pollinated using pollen from each subtypes

Mitochondrial genomes of two CMS-S subtypes

Because mitochondrial genome structure, gene organization, and linearization all affect sterility, the mitochondrial genomes of the two CMS-S subtypes were sequenced from highly purified mitochondrial DNA. Previous studies have demonstrated that the majority of the CMS-S genome is present in the linear form [6]. However, in this paper, the assemblies of the subtypes’ genomes are depicted in circular form for better illustration. The two newly sequenced CMS-S subtype genomes showed high similarity with each other and with the reference genome [15] (Additional file 2: Fig. S1), but exhibited variable sizes with 557,050-bp in CMS-Sa and 553,762-bp in CMS-Sb (Fig. 2). Variations in large repeats (> 0.5-kb) accounted for the majority of the size differences between the mitochondrial genomes of subtype CMS-Sa and CMS-Sb. There are 22 predicted large repeats in the mitochondrial genomes of CMS-S maize [15]. The two newly sequenced mitochondrial genomes shared the same size and copy number of all repeats except for a 17-kb repeat region designated as R17. The first copy of R17 was identical among CMS-Sa, CMS-Sb, and the reference genome, whereas the second copy of R17 had a deletion of a 3.3-kb sequence from its internal region in the mitochondrial genome of CMS-Sb (Additional file 3: Fig. S2 and Additional file 4: Fig. S3). Most of the functional genes/proteins were identical at both the nucleotide and amino acid sequence levels in the two subtypes. However, apocytochrome b (cob) and NADH dehydrogenase subunit 1 (nad1) differed between the genomes of the two subtypes. The cob gene of CMS-Sa had a T to G substitution at position + 312, resulting in a predicted amino acid change from phenylalanine to leucine. There were two copies of nad1-exon1 in CMS-Sa, but only one copy of nad1-exon1 in CMS-Sb because of the deletion of the 3.3-kb sequence.

Fig. 2
figure 2

Graphical genome maps display major difference between CMS-Sa and CMS-Sb. Both genomes have identical sequence of 4.2-kb repeat regions that contains the sterilizing gene orf355. The second copy of 17-kb repeat of CMS-Sb lost 3.3-kb sequence, leading to the elimination of one copy of first exon of nad1. Yellow rectangles represent coding region of nad1. Repeat regions sequences were labeled in brown color

Nucleotide substitutions mainly occur in the plastid-derived sequence in newly sequenced CMS-S subtypes

Alignments with the reference genome revealed 117 and 104 nucleotide substitutions in CMS-Sa and CMS-Sb, respectively. The nucleotide substitution frequency in CMS-Sa (2.1 substitutions/10,000 bp) was higher than that in CMS-Sb (1.88 substitutions/10,000 bp) (Additional file 5: Table S2). It is thought that there is a lower frequency of indels than substitutions in plant mitochondrial genomes [15]. Consistent with this, the CMS-Sa and CMS-Sb mitochondrial genomes had a low ratio of indels to substitutions: 0.17 and 0.15 indels/substitution relative to the reference genome, respectively. More than half of the indels were five nucleotides in length. Length variations in simple sequence repeats (SSRs) accounted for 41% of the total indels (Additional file 6: Table S3 and Additional file 7: Table S4).

Compared with the reference genome, CMS-Sa and CMS-Sb shared 92 identical nucleotide substitutions. These substitution sites were not evenly distributed, but clustered in several hot spot regions (Fig. 3a). Notably, most of these hot spots were in plastid-derived sequences (mitochondrial plastid DNAs: mtpt). For example, 22 of the substitutions were within the 3726-bp mtpt sequence that was present in the mitochondrial genome of both subtypes. However, the corresponding 3681-bp homolog sequences in the plastid were almost identical among CMS-Sa, CMS-Sb, and the reference plastid genome (NC_001666.2) (Fig. 3a). To eliminate the possible interference of plastid genomes during sequencing and assembly, we amplified these mtpt and their corresponding plastid sequences separately from the mitochondrial genome and the plastid genome using suites of specific primers. The results indicated that the analyzed mtpt sequences in CMS-Sa and CMS-Sb were different from their ancestral homolog sequences in the plastid. Thus, nucleotide substitutions have occurred in the mitochondrial genomes of CMS-Sa and CMS-Sb, especially in plastid-derived sequences. The substitution rates in mtpt sequence were about 42.3 and 42.8 substitutions/10,000-bp in CMS-Sa and CMS-Sb, respectively, significantly higher than the rates in non-mtpt sequences (Fig. 3b). Most of the clustered substitution sites surrounded the functional plastid-originated tRNA coordinates (Additional file 5: Table S2).

Fig. 3
figure 3

Nucleotide substitutions mainly occurred in the plastid-derived sequence in new sequenced CMS-S subtypes. a Mutation bias in plastid-derived DNA sequences (mtpt) of CMS-S mitochondrial genome. Orange dots in upper panel represent nucleotides mutation in each subtypes that compared with referee. Mtpt sequences are labeled in green color. Lower panel shows individual nucleotide mutations in the representative 3658-bp mtpt sequence and chloroplast homolog among NB, CMS-Sa, and CMS-Sb. Mt, mitochondrial; Pt, plastid. b Mutation rate in mtpt and non-mtpt sequence of CMS-Sa and CMS-Sb, calculated as the number of substitution or indel per 10,000 nucleotide

Transcript levels of orf355 in the two CMS-S subtypes

Next, we focused on differences in the sterility-associated region between the two CMS-S subtypes. In the CMS-S mitochondrial genome, there are two copies of a 4.2-kb repeat (R4.2) containing orf355-orf77 sequences. The second copy of R4.2 located downstream of cox2 recombines with S plasmid to yield linearized CMS-S genomes, from which the sterilizing 1.6-kb RNA is transcribed [8]. Previous studies demonstrated that expression of orf355-orf77 is associated with the male sterility of maize CMS-S. Sequencing of the coding region revealed that the DNA sequence of R4.2 was identical in two subtypes (Fig. 2). Most spontaneous reversions in CMS-S maize are associated with the integrity of the free S plasmid and TIR sequence of main mitochondrial genome. However, the newly sequenced subtypes CMS-Sa and CMS-Sb shared identical intact S1, S2 and TIR sequences with the reference sequences [40, 41]. Next, we checked the transcript levels of orf355 and orf77 in the leaves and anthers of the two accessions, normalized against the mitochondrial rRNA gene, rrn18. Quantitative RT-PCR analyses demonstrated that the orf355 transcript levels were significantly higher in the anthers than in the leaves in both subtypes. The transcript abundance of orf355 was even slightly lower in CMS-Sb than in CMS-Sa (Fig. 4). Thus, other factors beside the sterilizing factor may also influence the stability of sterility in CMS-S maize.

Fig. 4
figure 4

Transcript level of sterilizing gene in CMS-Sa and CMS-Sb. Transcript level of orf355 and orf77 in leaf and anther in NB, CMS-Sa, and CMS-Sb. Data are mean ± s.e.m. (n = 3). Different letters denote significant differences determined by Tukey’s tests, P < 0.05

Similar post-transcriptional processes in the R17 region in CMS-Sa and CMS-Sb

The 3.3-kb deletion in the second copy of R17 accounted for the major differences between the two CMS-S subtypes. R17 contains several functional genes, including five tRNA genes, the first exon of nad1, and the fourth and fifth exons of nad2. Notably, a 411-bp mtpt sequence was located within the intergenic region between nad2-exon5 and nad1-exon1. The site of the deletion in CMS-Sb was between the mtpt sequence and trnP, leading to the elimination of nad1-exon1 (Fig. 5). In the maize mitochondrial genome, nad1 and nad2 were interrupted by group II introns (Fig. 6a). Mature functional RNA is generated through posttranscriptional modification of the corresponding precursors [23]. A previous study using NB mitochondrial maize illustrated that nad1-exon1 is co-transcribed with upstream tRNAs to generate polycistronic precursor RNAs [39]. Following endonucleolytic cleavages at specific sites adjacent to tRNAs, mature RNA of nad1-exon1 is released. R17 is present in all cytotypes of maize except for CMS-T [15], which implies that CMS-S probably has identical transcripts to those of NB maize. We performed cRT-PCR to explore the post-transcriptional processing of genes in R17. The results revealed that nad1-exon1 was processed downstream of trnP to generate the 5′ terminus and upstream of the mtpt sequence to generate the 3′ terminus (Fig. 6b). Because the first copy of R17 was intact in both subtypes, there was no difference in the fragment length and processing site of cRT-PCR products between them (Fig. 6b and Additional file 8: Fig. S4a). Similarly, cRT-PCR analysis only amplified transcripts flanking nad2-exon3, exon4, and exon5 adjacent to the first copy of R17. No nad2 transcripts specific to the second copy of R17 were identified (Fig. 6c-d and Additional file 8: Fig. S4b). Thus, according to our results, there was no obvious differences in post-transcriptional processing between CMS-Sa and CMS-Sb.

Fig. 5
figure 5

Schematic map of the second copy of 17-kb repeat in CMS-Sa and CMS-Sb. Coding genes are indicated by boxes of purple color. The mtpt sequence and its homolog in chloroplast genome are labeled in green color. The eliminated 3.3-kb sequence in CMS-Sb is labeled in gray color. The putative nad1 transcripts were cleaved at the sites adjacent to trnP to generate the 5′ end of precursor. 3′ end of nad1 precursor were tested at the sites adjacent to mtpt sequence

Fig. 6
figure 6

Similar post-transcriptional process in R17 region was observed in CMS-Sa and CMS-Sb. a Diagrammatic representation of nad1 and nad2 RNA transcripts that generated from post-transcriptionally processing of the precursor. b and c The cropped gels of cRT-PCR products of nad1-transcripts 1 (b) and nad2-transcripts 2 (c). “+” indicated RNA samples treated with RNA 5′-polyphosphatase. d Schematic map of two copies of R17 in CMS-Sb. Circular RT-PCR (cRT-PCR) in CMS-S revealed that nad1-exon1 were processed downstream of trnP to generate the 5′ terminus and upstream of mtpt sequence to generate the 3′ terminus. No obvious difference was observed in the fragment length and processing site of cRT-PCR products between subtype CMS-Sa and CMS-Sb, which due to the presence of intact first copy of R17 is both genomes. cRT-PCR analysis revealed similar precursor of nad2-exon3, exon4, and exon5 that flanking border of the first copy of R17. Consistent with previous work of Zhang et al, no specific precursor that flanking nad2-exon4 and exon5 in the second copy of R17 were identified, which probably due to low abundance of this precursor

Sterility-stable subtype CMS-Sb displayed reduced mitochondrial function in anthers as compared with CMS-Sa

Mitochondrial function increases in anthers to meet the large energy demands of developing pollen. We determined the ratio of the mitochondrial genome to the nuclear genome to estimate the average mitochondrial copy number per cell in different tissues. The results of qRT-PCR analyses suggested that, on average, there were more mitochondrial genome copies per cell in the anther tissues than in seedlings (Fig. 7a). The mitochondrial copy number in anthers was lower in CMS-Sb than in CMS-Sa, implying that mitochondrial biogenesis was impaired in CMS-Sb (Fig. 7a). Transcript levels of mitochondrial genes were also analyzed. The transcript levels of most mitochondrial genes were significantly increased in anthers of CMS-Sa, but not in anthers of CMS-Sb (Fig. 7d). For example, the transcript levels of nad1 and nad2 were more than four times higher in anthers than in seedlings of CMS-Sa, but not significantly different between anthers and seedlings in CMS-Sb (Fig. 7b, c). These results suggest that lower mitochondrial function in CMS-Sb maize during pollen development enhances the stability of its sterility.

Fig. 7
figure 7

The sterility-stable subtype CMS-Sb displayed reduced mitochondrial function in anther as compared CMS-Sa. a Mitochondria biogenesis in leaf and anther of fertile, CMS-Sa, and CMS-Sb, as illustrated by relative abundance of mitochondrial DNA versus nuclear DNA. Mitochondrial DNA levels were measured by two independent primer pairs of mitochondria, then normalized to the internal standard nucleus-encoded single-copy fragment. Data are mean ± s.e.m. (n = 6). b and c Transcripts abundances of representative mitochondria encoding genes nad1 (B) and nad2 (C) of Complex I. Abundances shown are relative to the level of CMS-sb (set to 1). Data are mean ± s.e.m. (n = 3). d Relative transcripts levels of mitochondria encoding genes in leaf and anther tissues. Transcripts levels were normalized to the minimum means among four groups (set to 1). The relative positions of genes as indicated by horizontal lines in the mitochondrial schematic map of CMS-S


The mitochondrion is a dynamic organelle that modulates their function and biogenesis in response to fluctuating energy demands triggered by developmental signals and environmental stimuli [24]. One key property of mitochondrion is its robustness, that is the ability to maintain the respiration rate and ATP synthesis over a wide regime, even in face of some severe mutations in the mitochondrial genome. Robustness of mitochondrion has important implications for understanding why some mutations in mitochondrial genomes are lethal and some only manifest in specific tissues [24]. Working models to explain such phenomena have been proposed on the basis of results from studies on animals [25]. One explanation is that the consequences of defects in oxidative phosphorylation (OXPHOS) complexes depend on severity of the defect and the biological threshold of the tissue. Thus, in a given tissue with high biochemical threshold, even a small quantity of normal mitochondria are sufficient to maintain a normal level of oxidative phosphorylation. Conversely, in another tissue with a low biochemical threshold, a very small decrease in OXPHOS function can induce mitochondrial collapse [25]. This hypothesis can explain the unstable sterility in CMS-S maize.

The number of mitochondrial copy number per diploid cell is large and highly variable during plant development, whereas the copy number in the nucleus remains essentially constant [42]. According to early proposal, high copy number of mtDNA reflects an increased demand for organellar ribosomes that can only be satisfied by increased rRNA gene number that results from genome amplification [43]. This probably reflect the capacity of mitochondrion to regulate gene expression in response to changed physiological condition. Studies on animals suggest that the increased mtDNA copy number is associated with ROS-generated oxidative stress [44]. Previous study detailed the ROS content change during pollen development [10]. Plant microsporogenesis is an energy-demanding process, but ATP required for pollen development is supplied mainly by mitochondria, which are non-photosynthetic [45]. Thus, large amount ROS were generated as inevitable byproduct of oxidative phosphorylation. Previous studies have reported that induction of mitochondrial biogenesis can result in a 20- to 40-fold increase in some anther cells in maize, and that this is accompanied by increased expression levels of several mitochondrial genes in microspores [45, 46]. We detected increased mitochondrial copy number and genes expression in CMS-Sa, although lower than those in the same tissues of NB. Therefore, anther tissue of CMS-Sa maintains considerable mitochondrial robustness. However, we cannot exclude the possibility that the increased mtDNA copy number is caused by orf355-mediated mitochondrial degradation, which also associated with the increase of mtDNA content. Further studies are needed to dissect this mechanism.

Although the precise molecular mechanism of CMS-S is as yet unknown, CMS-S maize is the only system where restorers are reported to arise in real time. This implies that the defect conferred by orf355 may be not as severe as other reported CMS genes in stable CMS system [47,48,49]. Orf355 expression increased during pollen development, leading to the gradual impairment of mitochondrial function [10]. The final degeneration of most microspores may result from the inability of mitochondria to counteract the deleterious effect of orf355. However, mitochondrial function varies among tissues and even among different cells in the same tissue. Thus, anthers in CMS-Sa with mitochondrial function that high enough to reach the threshold for normal pollen development can revert to being fertile. In this way, the unstable sterility in the CMS-Sa subtype probably reflects mitochondrial robustness that successfully counteract the deleterious effect of orf355. All reported maize mitochondrial genomes contain a basic suite of functional genes encoding components of macromolecular complexes, but the copy number of individual genes or exons varies among cytotypes [15, 23]. Variations in the copy number of these essential OXPHOS genes can affect transcription. We observed higher transcript levels of mitochondrial genes in the subtype with two copies of nad1-exon1, but lower levels in the subtype with one copy. We speculated that changes in OXPHOS gene transcription may affect translation or assembly efficiency of mitochondrial complex, ultimately affecting the functional robustness of the respiratory chain [50, 51]. The intact one copy of nad1-exon1 in CMS-Sb may be sufficient to meet the energy requirements for normal vegetative growth, but not to meet the increase energy requirements of pollen development. Thus, this defect is manifested only in specific situations such as sterility reversion where robust mitochondrial function is required (Fig. 8).

Fig. 8
figure 8

A hypothesis about the regulation of stability of sterility in CMS-S. Sterility in CMS-S is the consequence of inability of mitochondrial function to counteract the deleterious effect of orf355. Mitochondrial activity increases in anther tissue to meet the large amount energy required for pollen development, which is a representative example of mitochondrial robustness (indicated by blue arrows). However, mitochondrial activity varies among tissues and even among different cells of the same tissue. In anthers of CMS-Sa, microspores with moderate mitochondrial activity that higher enough to reach the threshold (dotted line) for normal pollen development can revert to being fertile. In CMS-Sb, the loss of one copy of nad1-exon1 reduces the mitochondrial robustness. The mitochondrial activity of all microspores could not reach the threshold required for normal pollen development. This defect may be sufficient to sustain the normal vegetative growth, but is manifested in specific situations such as sterility reversion

One point should be noted that the unstable sterility discussed here is different from the spontaneous reversion reported by previous studies. Usually, a fertile plant with CMS-S cytoplasm was caused by restorer-of-fertility gene such as Rf3, Rf9 and rfl1, or by S1- and S2-plasmid integration into mitochondrial genome. Here, a slight modification of mitochondrial genome that is not associated with CMS gene exerts additive effect on male sterility expression when combined with the CMS gene. In other words, the deletion of nad1-exon1 in the second copy of R17 plays as an enhancer in terms of genetics. Practical observations indicated that most of these plants with exerted anthers did not fully revert to being fertile, although a small portion of the pollen in the exerted anthers was stainable with I2-KI solution. In addition, the exerted anthers in CMS-Sa usually appeared several days after natural pollination, so they would not meet the functional silks to complete fertilization. Nevertheless, the presence of starch granules in these pollen grains implied that they had the potential to germinate, and represented the partial recovery of mitochondrial function in the supposedly sterile pollen.

Another notable feature of the two newly sequenced mitochondrial genomes is the nucleotide substitution bias in plastid-derived sequences. One important feature of angiosperm mitochondrial genomes is that they have acquired foreign DNA from diverse sources, especially by intracellular gene transfer from the plastids [52,53,54]. Although most of the foreign plastid DNA results in non-functional sequences, some of them play important role during plant development, e.g., providing tRNAs in mitochondrial protein synthesis, the creation of promoter regions and codons, and the participation in post-transcriptional RNA processing [55,56,57]. In this study, we detected higher nucleotide mutation rates in mtpt sequences than in non-mtpt regions. In addition, we detected a 3.3-kb sequence deletion adjacent to a 411-bp mtpt sequence. One explanation for this result is that the mtpt region is recombinationally active and the transfer of plastid DNA into the mitochondrion is a relatively common event [23]. In this scenario, mtpt sequences would recombine with homologous DNA from the plastid. In rare cases, replication of mtpt sequences may fail to utilize the recombinational repair machinery, resulting in the occasional introduction of mutations into the mitochondrial genome.

In this study, we found that mtpt sequences may serve as recombinationally active sites in mediating mitochondrial genome rearrangement, leading to the copy number variation of mitochondrial gene. In addition, our comparative analyses provide evidence for a relationship between the copy number of mitochondrial gene and functional plasticity, which gives new insights into the stability of male sterility in maize CMS-S.

Availability of data and materials

The obtained sequence data have been deposited into NCBI database with Accession No. PRJNA864113,


  1. Chase CD. Cytoplasmic male sterility: a window to the world of plant mitochondrial-nuclear interactions. Trends Genetics. 2007;23(2):81–90.

    Article  CAS  Google Scholar 

  2. Chen L, Liu YG. Male sterility and fertility restoration in crops. Annu Rev Plant Biol. 2014;65:579–606.

    Article  CAS  PubMed  Google Scholar 

  3. Ullstrup AJ. The impacts of the southern corn leaf blight epidemics of 1970-1971. Annu Rev Phytopathol. 1972;10(1):37–50.

    Article  Google Scholar 

  4. Levings CS. The Texas cytoplasm of maize: cytoplasmic male sterility and disease susceptibility. Science (New York, NY). 1990;250(4983):942–7.

    Article  CAS  Google Scholar 

  5. Levings CS 3rd, Kim BD, Pring DR, Conde MF, Mans RJ, Laughnan JR, et al. Cytoplasmic reversion of cms-S in maize: association with a Transpositional event. Science (New York, NY). 1980;209(4460):1021–3.

    Article  CAS  Google Scholar 

  6. Schardl CL, Lonsdale DM, Pring DR, Rose KR. Linearization of maize mitochondrial chromosomes by recombination with linear episomes. Nature. 1984;310(5975):292–6.

    Article  CAS  Google Scholar 

  7. Pring DR, Levings CS, Hu WW, Timothy DH. Unique DNA associated with mitochondria in the "S"-type cytoplasm of male-sterile maize. Proc Natl Acad Sci U S A. 1977;74(7):2904–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Gabay-Laughnan S, Kuzmin EV, Monroe J, Roark L, Newton KJ. Characterization of a novel thermosensitive restorer of fertility for cytoplasmic male sterility in maize. Genetics. 2009;182(1):91–103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Zabala G, Gabay-Laughnan S, Laughnan JR. The nuclear gene Rf3 affects the expression of the mitochondrial chimeric sequence R implicated in S-type male sterility in maize. Genetics. 1997;147(2):847–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Xiao S, Zang J, Pei Y, Liu J, Liu J, Song W, et al. Activation of mitochondrial orf355 gene expression by a nuclear-encoded DREB transcription factor causes cytoplasmic male sterility in maize. Mol Plant. 2020;13(9):1270–83.

    Article  CAS  PubMed  Google Scholar 

  11. Qin X, Tian S, Zhang W, Zheng Q, Wang H, Feng Y, et al. The main restorer Rf3 of maize S type cytoplasmic male sterility encodes a PPR protein that functions in reduction of the transcripts of orf355. Mol Plant. 2021;14(12):1961–4.

    Article  CAS  PubMed  Google Scholar 

  12. Gabay-Laughnan S, Settles AM, Hannah LC, Porch TG, Becraft PW, McCarty DR, et al. Chamusco KC et al: Restorer-of-Fertility Mutations Recovered in Transposon-Active Lines of S Male-Sterile Maize. G3 (Bethesda, Md). 2018;8(1):291–302.

    Article  CAS  Google Scholar 

  13. Matera JT, Monroe J, Smelser W, Gabay-Laughnan S, Newton KJ. Unique changes in mitochondrial genomes associated with reversions of S-type cytoplasmic male sterility in maizemar. PLoS One. 2011;6(8):e23405.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Wen L, Ruesch KL, Ortega VM, Kamps TL, Gabay-Laughnan S, Chase CD. A nuclear restorer-of-fertility mutation disrupts accumulation of mitochondrial ATP synthase subunit alpha in developing pollen of S male-sterile maize. Genetics. 2003;165(2):771–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Allen JO, Fauron CM, Minx P, Roark L, Oddiraju S, Lin GN, et al. Comparisons among two fertile and three male-sterile mitochondrial genomes of maize. Genetics. 2007;177(2):1173–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Gualberto JM, Newton KJ. Plant mitochondrial genomes: dynamics and mechanisms of mutation. Annu Rev Plant Biol. 2017;68:225–52.

    Article  CAS  PubMed  Google Scholar 

  17. Stern DB, Palmer JD. Recombination sequences in plant mitochondrial genomes: diversity and homologies to known mitochondrial genes. Nucleic Acids Res. 1984;12(15):6141–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Newton K. Stern D. Gabay-Laughnan S: Mitochondria and Chloroplasts; 2009. p. 481–503.

    Google Scholar 

  19. Marienfeld JR, Newton KJ. The maize NCS2 abnormal growth mutant has a chimeric nad4-nad7 mitochondrial gene and is associated with reduced complex I function. Genetics. 1994;138(3):855–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Bellaoui M, Martin-Canadell A, Pelletier G, Budar F. Low-copy-number molecules are produced by recombination, actively maintained and can be amplified in the mitochondrial genome of Brassicaceae: relationship to reversion of the male sterile phenotype in some cybrids. Mol Gen Genet. 1998;257(2):177–85.

    Article  CAS  PubMed  Google Scholar 

  21. Janska H, Sarria R, Woloszynska M, Arrieta-Montiel M, Mackenzie SA. Stoichiometric shifts in the common bean mitochondrial genome leading to male sterility and spontaneous reversion to fertility. Plant Cell. 1998;10(7):1163–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Chen J, Guan R, Chang S, Du T, Zhang H, Xing H. Substoichiometrically different mitotypes coexist in mitochondrial genomes of Brassica napus L. PLoS One. 2011;6(3):e17662.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Clifton SW, Minx P, Fauron CM, Gibson M, Allen JO, Sun H, et al. Sequence and comparative analysis of the maize NB mitochondrial genome. Plant Physiol. 2004;136(3):3486–503.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Vafai SB, Mootha VK. Mitochondrial disorders as windows into an ancient organelle. Nature. 2012;491(7424):374–83.

    Article  CAS  PubMed  Google Scholar 

  25. Rossignol R, Malgat M, Mazat JP, Letellier T. Threshold effect and tissue specificity. Implication for mitochondrial cytopathies. J Biol Chem. 1999;274(47):33426–32.

    Article  CAS  PubMed  Google Scholar 

  26. Xiao H, Zhang F, Zheng Y. The 5′ stem-loop and its role in mRNA stability in maize S cytoplasmic male sterility. Plant J. 2006;47(6):864–72.

    Article  CAS  PubMed  Google Scholar 

  27. Zhang D, Wu S, An X, Xie K, Dong Z, Zhou Y, et al. Construction of a multicontrol sterility system for a maize male-sterile line and hybrid seed production based on the ZmMs7 gene encoding a PHD-finger transcription factor. Plant Biotechnol J. 2018;16(2):459–71.

    Article  CAS  PubMed  Google Scholar 

  28. Newton KJ. Procedures for isolating mitochondria and mitochondrial DNA and RNA: The Maize Handbook; 1994.

  29. Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJ, Birol I. ABySS: a parallel assembler for short read sequence data. Genome Res. 2009;19(6):1117–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv e-prints 2013.

  31. Tillich M, Lehwark P, Pellizzer T, Ulbricht-Jones ES, Fischer A, Bock R, et al. GeSeq - versatile and accurate annotation of organelle genomes. Nucleic Acids Res. 2017;45(W1):W6–w11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Greiner S, Lehwark P, Bock R. OrganellarGenomeDRAW (OGDRAW) version 1.3.1: expanded toolkit for the graphical visualization of organellar genomes. Nucleic Acids Res. 2019;47(W1):W59–w64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Delcher AL, Salzberg SL, Phillippy AM. Using MUMmer to identify similar regions in large sequence sets. Current protocols in bioinformatics 2003, Chapter 10:Unit 10.13.

  34. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics (Oxford, England). 2009;25(15):1972–3.

    Article  Google Scholar 

  36. Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32(1):268–74.

    Article  CAS  PubMed  Google Scholar 

  37. Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14(6):587–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Ma J, Li XQ. Organellar genome copy number variation and integrity during moderate maturation of roots and leaves of maize seedlings. Curr Genet. 2015;61(4):591–600.

    Article  CAS  PubMed  Google Scholar 

  39. Zhang Y, Huang X, Zou J, Liao X, Liu Y, Lian T, et al. Major contribution of transcription initiation to 5′-end formation of mitochondrial steady-state transcripts in maize. RNA Biol. 2019;16(1):104–17.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Paillard M, Sederoff RR, Levings CS. Nucleotide sequence of the S-1 mitochondrial DNA from the S cytoplasm of maize. EMBO J. 1985;4(5):1125–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Levings CS, Sederoff RR. Nucleotide sequence of the S-2 mitochondrial DNA from the S cytoplasm of maize. Proc Natl Acad Sci U S A. 1983;80(13):4055–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Oldenburg DJ, Bendich AJ. DNA maintenance in plastids and mitochondria of plants. Front Plant Sci. 2015;6:883.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Bendich AJ. Why do chloroplasts and mitochondria contain so many copies of their genome? BioEssays. 1987;6(6):279–82.

    Article  CAS  PubMed  Google Scholar 

  44. Lee HC, Wei YH. Mitochondrial biogenesis and mitochondrial DNA maintenance of mammalian cells under oxidative stress. Int J Biochem Cell Biol. 2005;37(4):822–34.

    Article  CAS  PubMed  Google Scholar 

  45. Huang J, Struck F, Matzinger DF, Levings CS 3rd. Flower-enhanced expression of a nuclear-encoded mitochondrial respiratory protein is associated with changes in mitochondrion number. Plant Cell. 1994;6(3):439–48.

    CAS  PubMed  PubMed Central  Google Scholar 

  46. Warmke HE, Lee SL. Pollen abortion in T cytoplasmic male-sterile corn (Zea mays): a suggested mechanism. Science (New York, NY). 1978;200(4341):561–3.

    Article  CAS  Google Scholar 

  47. Rhoads DM, Levings CS 3rd, Siedow JN. URF13, a ligand-gated, pore-forming receptor for T-toxin in the inner membrane of cms-T mitochondria. J Bioenerg Biomembr. 1995;27(4):437–45.

    Article  CAS  PubMed  Google Scholar 

  48. Luo D, Xu H, Liu Z, Guo J, Li H, Chen L, et al. A detrimental mitochondrial-nuclear interaction causes cytoplasmic male sterility in rice. Nat Genet. 2013;45(5):573–7.

    Article  CAS  PubMed  Google Scholar 

  49. Duroc Y, Gaillard C, Hiard S, Defrance MC, Pelletier G, Budar F. Biochemical and functional characterization of ORF138, a mitochondrial protein responsible for Ogura cytoplasmic male sterility in Brassiceae. Biochimie. 2005;87(12):1089–100.

    Article  CAS  PubMed  Google Scholar 

  50. Vercellino I, Sazanov LA. The assembly, regulation and function of the mitochondrial respiratory chain. Nat Rev Mol Cell Biol. 2022;23(2):141–61.

    Article  CAS  PubMed  Google Scholar 

  51. Woodson JD, Chory J. Coordination of gene expression between organellar and nuclear genomes. Nat Rev Genet. 2008;9(5):383–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Palmer JD, Adams KL, Cho Y, Parkinson CL, Qiu YL, Song K. Dynamic evolution of plant mitochondrial genomes: mobile genes and introns and highly variable mutation rates. Proc Natl Acad Sci U S A. 2000;97(13):6960–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Stern DB, Lonsdale DM. Mitochondrial and chloroplast genomes of maize have a 12-kilobase DNA sequence in common. Nature. 1982;299(5885):698–702.

    Article  CAS  PubMed  Google Scholar 

  54. Ellis J. Promiscuous DNA--chloroplast genes inside plant mitochondria. Nature. 1982;299(5885):678–9.

    Article  CAS  PubMed  Google Scholar 

  55. Wang D, Rousseau-Gueutin M, Timmis JN. Plastid sequences contribute to some plant mitochondrial genes. Mol Biol Evol. 2012;29(7):1707–11.

    Article  CAS  PubMed  Google Scholar 

  56. Nakazono M, Nishiwaki S, Tsutsumi N, Hirai A. A chloroplast-derived sequence is utilized as a source of promoter sequences for the gene for subunit 9 of NADH dehydrogenase (nad9) in rice mitochondria. Mol Gen Genet. 1996;252(4):371–8.

    CAS  PubMed  Google Scholar 

  57. Miyata S, Nakazono M, Hirai A. Transcription of plastid-derived tRNA genes in rice mitochondria. Curr Genet. 1998;34(3):216–20.

    Article  CAS  PubMed  Google Scholar 

Download references


We thank Jennifer Smith, PhD, from Liwen Bianji (Edanz) ( for editing the English text of a draft of this manuscript. We thank Dr. Tianjun Xu and Fengchao Jiang for their valuable comments and discussion on previous versions of the manuscript. We are also grateful to anonymous reviewers for their helpful suggestions and comments.


This work was financially supported by Youth Research Fund (QNJJ202128) and The Beijing Scholars Program (BSP041).

Author information

Authors and Affiliations



SX, WS and JZ conceived the study. SX and TN performed the experiment. AS and YZ assembled, annotated, and analyzed the genome. SX drafted the manuscript. JX and RZ provided input to the manuscript. SX and WS revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Wei Song or Jiuran Zhao.

Ethics declarations

Ethics approval and consent to participate

Experimental research and field studies comply with relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declared that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information


Additional file 1.


Additional file 2.


Additional file 3.


Additional file 4.


Additional file 5.


Additional file 6.


Additional file 7.


Additional file 8.


Additional file 9.


Additional file 10.


Additional file 11.


Additional file 12.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xiao, S., Xing, J., Nie, T. et al. Comparative analysis of mitochondrial genomes of maize CMS-S subtypes provides new insights into male sterility stability. BMC Plant Biol 22, 469 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Maize
  • CMS-S
  • Sterility stability
  • Mitochondria
  • Robustness