- Research article
- Open Access
Wide crossing diversify mitogenomes of rice
BMC Plant Biology volume 20, Article number: 159 (2020)
In most angiosperms, the inheritance of the mitochondria takes place in a typical maternal manner. However, very less information is available about if the existence of structural variations or not in mitochondrial genomes (mitogenomes) between maternal parents and their progenies.
In order to find the answer, a stable rice backcross inbred line (BIL) population was derived from the crosses of Oryza glaberrima/Oryza sativa//Oryza sativa. The current study presents a comparative analysis of the mitogenomes between maternal parents and five BILs. There were recorded universal structural variations such as reversal, translocation, fusion, and fission among the BILs. The repeat-mediated recombination and non-homologous end-joining contributed virtually equal to the rearrangement of mitogenomes. Similarly, the relative order, copy-number, expression level, and RNA-editing rate of mitochondrial genes were also extensively varied among BILs.
These novel findings unraveled an unusual mystery of the maternal inheritance and possible cause for heterogeneity of mitogenomes in rice population. The current piece of work will greatly develop our understanding of the plant nucleo-cytoplasmic interaction and their potential role in plant growth and developmental processes.
A coordinated nucleo-cytoplasmic interaction plays a crucial role in the plant development [1, 2]. However, inter- or intraspecies hybridization in plants may lead to the reshuffling of the nuclear genome and eventually might disturb the equilibrium of cytonuclear interactions. Therefore, maintenance of a coordinated nucleo-cytoplasmic interaction would determine the growth, development, and even survival of the plants [3, 4]. Unlike the nuclear genome, mitochondria and chloroplast generally follow a uniparental inheritance in plants [5, 6]. It is considered that except some special cases, maternal inheritance appears to be the predominant mode of mitogenome transmission in plants especially angiosperms, and cycads and gnetophytes (gymnosperms) [7,8,9]. Therefore, it becomes very essential to understand the mechanism of acclimatization of the maternally inherited mitogenomes with the reshuffled nuclear genome to retain a determined nucleo-cytoplasmic interaction during plant evolution and crop domestication. Here, we speculate that the mitogenome may have evolved a special genetic mechanism in order to cope with the frequent rearrangement of nuclear genomes in nature. Interestingly, previous reports have also documented the apparent differentiation of mitogenome in plant populations including Silene vulgaris, Cucurbitaceae, and Actinidia families [10,11,12,13]. It is proposed that mitogenome possibly transmits in a dynamic way in the filial progenies to reinstate a balance of mitonuclear interaction in plants. Otherwise, it may be hard to explain the rich diversity of mitogenomes within a plant species population only by means of traditional maternal inheritance.
Rice, as a model cereal crop, shows a typical example of maternal inheritance. It is generally considered that indica rice was primarily originated from wild rice Oryza rufipogon in the middle of the Yangtz river, which was later domesticated into japonica rice in Northern China [14, 15]. However, recent genomic research suggests that both indica, as well as japonica rice, were originated from O. rufipogon simultaneously [16, 17]. It indicates that the mitogenomes of indica and japonica varieties are well-differentiated and show a narrow diversity in each subspecies if their mitogenomes strictly persist in stable maternal inheritance. However, rice population cannot only be differentiated into a group of subpopulations at nuclear level  but also significantly at mitogenome level. The mitogenomes of 224 worldwide landraces of rice were analyzed using 32 mitotype-specific molecular markers , all the landraces were classified into 20 subgroups, and the indica and japonica lines were found evenly distributed in each clade (Supplementary Fig. S1), implying a strong differentiation of rice mitogenomes [15, 17]. Such characteristics make rice as an apt model plant for the investigation of mitogenome inheritance and evolution patterns.
In order to prove this hypothesis, we developed a stable rice BC2F12 backcross inbred line (BIL) population which inherited their mitochondrial genomes from the same maternal parent (O. glaberrima × O. sativa 93–11). Subsequently, high-quality draft assemblies of mitogenomes were constructed and assayed for the genomic structure, gene content, DNA rearrangement, repeat sequences, as well as gene expression and RNA-editing among the maternal line and BILs with diverse genotypes. It was found that structure and organization of mitogenome, gene copy-number, expression, and RNA editing all varied to certain degrees among the BILs. The current novel findings revamp the traditional concept of maternal inheritance and greatly extend our understanding of the plant nucleo-cytoplasmic interactions.
Backcross inbred lines derived from crosses of O. glaberrima and O. sativa showed different mitotypes from the maternal line
In the quest for the existence of dynamic inheritance of rice mitogenome, two sets of BC2F12 backcross inbred lines (BIL) were constructed by crossing O. sativa 93–11 with O. glaberrima line 675 as maternal line, and then introgressed using 93–11 and O. glaberrima as the paternal line, respectively, in an attempt to simulate the natural evolutionary process of rice (Fig. 1a). Finally, 99 BILs from O. glaberrima/93–11//93–11 (P) and 82 BILs from O. glaberrima/93–11//O. glaberrima (M) were developed. These BILs were further analyzed for DNA polymorphism in the mitogenomes using mitochondria-specific molecular markers . Although most of the paternal, as well as maternal BILs, were found to be grouped together with O. glaberrima, about 35.4% of the paternal (Fig. 1b) and 12.2% of the maternal BILs (Fig. 1c) were apparently distinct from their maternal parent O. glaberrima. It indicates that the mitogenomes are differentiated among different BILs, despite their inheritance from the same maternal line.
Assembly of the mitogenomes of O. glaberrima and BILs
In order to decipher the structural variation of the mitogenomes, the BILs from the same cross (O. glaberrima/93–11//93–11) were selected for subsequent analysis owing to their distinct plant stature (Supplementary Fig. S2), differentiated mitotypes (Fig. 1b), and genotypes (Supplementary Fig. S3) verified with over 20-fold genomic sequencing (Supplementary Table S1).
In order to ascertain the accuracy of assembly, initially, a draft assembly of mitogenome of O. glaberrima was constructed using PacBio along with Illumina sequencing (Supplementary Table S2, S3). It contained two scaffolds, with minimum estimates of genome size of 402,174 bp, and 44.0% GC content (Table 1). The mitogenome of O. glaberrima was compared with that of Nipponbare (a high-quality reference mitogenome), and the boundaries of each bin were found (Supplementary Fig. S4a). PCR amplification and appearance of target boundaries (Supplementary Fig. S4b, Supplementary information 1) further advocate high reliability of O. glaberrima mitogenome.
Likewise, the draft assembly of mitogenomes of five BILs was constructed (Supplementary Table S2, S3), and each mitogenome contained 2–3 scaffolds, with nearly stable GC content (~ 43%) as the maternal line, and genome size varied from 409,553 to 441,199 bp, a little larger than the O. glaberrima (Table 1). The mitogenome gene region is relatively more stable than the non-gene region, and the length of the gene region increased from 1017 to 3666 bp, whereas the non-gene region increased from 6362 to 35,359 bp (Table 1), implying the non-gene region as the main cause for the expansion of BIL mitogenomes. Prediction analysis revealed that each BIL mitogenome contained identical sets of protein and rRNA genes as the O. glaberrima did, including 35 protein genes, 3 rRNA genes, and 18 tRNA genes (Supplementary Fig. S5).
Mitogenome of BILs showed distinctive structural variation
When the BIL mitogenomes were aligned to the O. glaberrima, the pairwise genome alignment showed that there were 5, 3, 6, 7 and 6 homologous segments in P10, P88, P90, P91, and P92, respectively, being mapped to O. glaberrima (Fig. 2). These findings are in agreement with the 3, 2, 5, 6 and 5 rearrangement events in P10, P88, P90, P91, and P92, respectively as detected by GRIMM . Comprehensively, the rearrangement events included 1 reversal and 2 translocations in P10, 1 translocation and 1 fusion in P88, 4 reversals and 1 translocation in P90, 2 reversals, 3 translocations and 1 fission in P91, and 2 reversals, 2 translocations and 1 fission in P92. The varied and inconsistent rearrangement pattern in mitogenomes signifies their strong diversity among the BILs.
Although alignment showed that the mitochondrial gene content is the same between BILs and maternal line, there were 15 genes with increased copy-number in at least one BIL except for the trnP (TGG) that showed a coinciding rise and fall in copy-number in different lines (Supplementary Table S4). Additionally, there were 20 new formed gene clusters in the BIL mitogenomes (Supplementary Table S6) due to the rearrangement. For example, the gene cluster trnE (TTC)-rrnL-rrn5-rrnS in O. glaberrima was converted into trnE (TTC)-rrnL-atp6-nad5 in P10 and P92. Similar changes were also identified in P10, P91, and P92 (Supplementary Table S6), indicating that the mitogenome structure is globally transformed in the BIL population. Further analysis also confirmed that all mitogenomic variations in BILs were ascended from the maternal mitogenome itself, because any 93–11-specific mtDNA fragment (Supplementary Table S7) could not be detected in BIL mitogenomes.
Homologous recombination and non-homologous end joining (NHEJ) equally contributed to mitogenome variation
It is well documented that rearrangements in angiosperm mitogenomes are largely mediated by repeat sequences, and their homologous recombination makes mitogenome a multipartite structure . Accordingly, we analyzed the type and distribution of the repeats ≥50 bp which are regarded as the main contributor to mitogenomic rearrangement in plants. In total, 118 repeats including 15 large repeats (> 1000 bp) and 99 intermediate repeats (≥50 and ≤ 1000 bp) were identified in O. glaberrima mitogenome by BLASTALL, covering 24.14% of the mitogenome (Supplementary Table S9). Based on the repeats and flanking sequences, all the repeats were classified into 33 types, comprising of 30 two-copy repeats, 2 three-copy repeats, and 1 ten-copy repeat (Supplementary Table S10).
Then, the 28 two-copy mitochondrial repeats with detectable flanking sequences were selected for recombination analysis. Since the two-copy repeats generate only two products after recombination (Fig. 3a), it makes them convenient for detection . When the reads of the 20 Kb PacBio sequencing library to the mitochondrial repeats were mapped, huge variations in the repeat configurations (repeat + flanking sequence) of the same repeat were detected among BILs. This was illustrated by reduced MRS10-cd to 30.1% in P88 and increased MRS5-cd to 372.5% in P92 relative to the maternal line (Supplementary Fig. S6). It infers that homologous recombination (repeat-mediated) largely occurred in the BIL (Fig. 3b), as further verified by qPCR analysis that both MRS8-cb and MRS8-ad in P88 mitogenome were 100 times more than that in O. glaberrima (Supplementary Fig. S7).
Generally, the recombination rate in P91 and P92 is higher than the others. The intermediate repeat-mediated recombination became dominant in the BIL mitogenomes, and the range of variation caused by intermediate repeats was apparently wider than the large repeats (Fig. 3b), reflecting a key role of intermediate repeats in regulating mitogenome structure. Further analysis also showed that the repeat-mediated homologous recombination can lead to the change in gene order (Supplementary Fig. S8), as shown by GJ01 in P10 and P91, consistent with gene cluster change in BIL lines (Supplementary Table S6).
Interestingly, when we detect if repeats are overrepresented at rearrangement boundaries, the overlapping rate ranges from 42 to 83% in different BILs (Fig. 3c; Supplementary Table S8), stating that nearly half of the rearrangement events are caused by non-homologous end joining (NHEJ), which also plays an important role in determining the mitogenome structure. This is consistent with the increase in copy-number of 15 mitochondrial genes in BILs (Supplementary Table S4), as the NHEJ is usually related to mitochondrial gene copy number in plants as well as animals [23, 24]. Further analysis of the recombination patterns also showed that one-third of the reversal, half of the translocation, and most of the fission events were involved in non-homologous recombination (Supplementary Table S5). These results confer the non-homologous recombination as one of the major driving forces for mitogenome rearrangement.
Genomic variations alter mitochondrial gene expressional profile in BILs
Although the recombination mostly occurred in the non-gene region, about 57.1% coding genes were located within 2 Kb upstream or downstream of repeat sequences (Supplementary Table S9), suggesting the role of mtDNA recombination in promoter structure and activity as reported in the family of Silene vulgaris . Therefore, transcriptome analysis combined with qPCR was carried out to evaluate the expression profile of mitochondrial genes. Transcriptome analysis showed that the expression level of functional mitochondrial genes among the BILs ranged from 0.27 to 4.16 times of the O. glaberrima (Fig. 4a). Of which, the expression of atp9, rps19 and rps2 reduced to over 50% than O. glaberrima, whereas, the expression of 15 genes including atp9, ccmB, ccmC, ccmFc, cox1, nad1, nad4L, nad6, rpl16, rpl5, rps1, rps12, rps3, rps4, and rps7 was increased over 50% than O. glaberrima (Fig. 4a). Further correlation studies between gene expression and repeats showed that all the genes with over 50% down-regulation and 73% of the genes with over 50% up-regulation were all adjacent to the repeat within 2 Kb. The qPCR analysis also showed that expression levels were positively correlated with their gene contents, up to 2 Kb proximal repeats (Supplementary Fig. S9), inferring an important role of repeats in regulating mitochondrial gene expression. In general, expression of the genes in P91 is higher than that of P92, followed by P10, P88, and P90, respectively (Fig. 4a); the exact mechanism is yet to be ascertained. This reflects the complex nature of expression regulation in mitochondria.
The mitochondrial genes in BILs showed different RNA-editing pattern
RNA-editing plays a key role in coordinating mito-nuclear interaction . Mitogenomic structural variation serves as a dynamic means for mitochondria to deal with the cytonuclear interaction. Therefore, the RNA-editing site and rate of the protein-coding genes with reads detected ≥50 times were assayed in each mitogenome (Supplementary information 2). In general, a total of 525 C-to-U and 4 U-to-C editing sites involved in 35 mitochondrial genes were identified in the coding sequences (Supplementary Table S11). Notably, about half of the editing sites were located at the second codon, and 459 editing sites were related to amino acid change (Supplementary Table S11) in 35 genes (Supplementary information 1). With reference to the maternal line, the rise and fall of the RNA-editing rate in the same site were frequently observed among BILs (Fig. 4b). This is enormously demonstrated by the sites of mttB-381 having an editing rate of 35.5% in O. glaberrima, but 50.5% in P10 and only 21.7% in P90 (Fig. 4b), suggesting the dynamic nature of RNA-editing even for the same site in different lines.
Currently, a total of 339 mitochondrial RNA-editing sites are known in the coding sequence of rice (https://www.ncbi.nlm.nih.gov/). When these RNA editing sites were cross-referenced with the Nipponbare, about 143 mitochondrial RNA editing sites were found to be unique (Table 2). In order to test the accuracy of the bioinformatics analysis, three editing sites (cox2–144, mttB-354 and nad7–1110) with a large disturbance were selected for re-sequencing, and all of them showed similar editing pattern as anticipated by bioinformatics analysis (Supplementary Table S12), signifying their consistency.
This is the first study on the mitogenome variation of maternal inheritance in plants that uncovers a novel mechanism of the organelle genome corresponding to adapt a reshuffled heterogenic nuclear genome. The transmitted mitogenome revealed a global rearrangement (mostly occurring as reversal) through repeat-mediated recombination and non-homologous end joining (NHEJ) (Fig. 2f) and extensive alterations (including gene copy-number, expression, and substantial RNA-editing sites) that occurred in the progenies. These findings indicate that the inheritance of mitogenome across the descendent plants is dynamic and different progenies hold varied mitotypes. Surprisingly, expression levels were not always positively related to copy-number in the tested BIL lines, which reflects a complicated mechanism of mitochondria to cope with the nuclear-cytoplasmic interaction in plants.
Wide crossing alters mitotypes of progenies
It is well documented that cytoplasmic inheritance in rice and most of the other higher plants is typically due to maternal inheritance. However, increasing reports of the complicated mitotypes in rice seems compromising this opinion [19, 26]. In the present research, paternal mtDNA was not transmitted to the progenies (Supplementary Table S7), but the mitotypes in BIL lines were strongly differentiated (Fig. 2). MtDNA clustering analysis clarified that the paternal-backcrossed BIL population was apparently much more differentiated than that of the maternal-backcrossed BIL population (Fig. 1). Even in the same BIL population, the lines carrying more DNA fragments from the paternal line showed stronger mtDNA differentiation than the others. These findings suggest that the heterogenic nuclear genome could promote variation in the mitogenomes so as to concert the nuclear-mitogenomic interaction, keeping normal development of the BIL lines (Supplementary Fig. S2A).
Common wild rice is the direct ancestor of cultivars. The cultivars and wild rice with AA-genome can generate hybrids through outcrossing to some extent. Naturally, cultivars are frequently overlapped with habitant of wild rice, especially O. rufipogon in Asia, which inevitably lead to genetic drift . Meanwhile, wild rice with AA-genome is the most important gene pool for rice improvement. The vast number of favorable genes such as Gn8.1 for grain yield, Bph14 and Bph15 for brown plant hopper , Xa21 and Xa23 for white leaf blight  have been transferred to cultivars by crossing with wild rice. The natural and artificial out-crossing brings many of the genetic exchanges between wild rice and cultivars at genome level . As exemplified by the cytoplasmic male sterility, a detrimental mitochondria-nuclear interaction caused by interspecific hybridization in the plant system like rice , cotton , rape seed , sorghum  etc., and harmonious nuclear-cytoplasmic interaction play a critical role in maintaining plant normal growth, fertility, survival, and even speciation . In the wild or artificial hybridization, frequent introgression from wild rice not only reshuffles the nuclear genome but also dynamically sorts the mitogenomes through rearrangement. It helps cope with the introgressed nuclear elements to maintain stability in the species population. Thus, it will drive the variation and evolution of the mitochondrial genomes.
Recombination of mitogenome is far from repeat sequences
Plant has a complex and dynamic mitogenome configuration , which prefer to undergo active inter- and intramolecular DNA exchange . For mitogenome, recombination is the main contributor of its genetic variation, leading to stoichiometric variation in copy number, intragenomic rearrangements, and heteroplasmic molecular populations . Genomic rearrangement occurs through homologous recombination and non-homologous end-joining (NHEJ), where homologous recombination is mainly subjected to active repeat sequences . Plant mitogenomes are characterized by abundant repeat sequences of different sizes and number. Previous reports reveal that the recombination-active large repeats are responsible for the highly variable structure of plant mtDNA . In this study, the number of intermediate repeats was at least six times of the large repeats (Supplementary Table S9). Furthermore, the intermediate repeats involved in recombination frequency were apparently more than that of the large repeats in O. glaberrima as well as BILs (Fig. 3), indicating that majority of the mtDNA recombination are owing to intermediate repeats.
Unlike the significant co-occurrence of homologous region boundaries and repeats in other plants [39,40,41], there was at least 27% of non-homologous recombination being detected in different BIL mitogenomes (Supplementary Fig. S3C). Especially for P92, only 42% homologous region boundaries corresponded to the repeat sequence, while the majority of the mtDNA rearrangement was caused through NHEJ. This fact was underestimated in the previous studies , reflecting that NHEJ plays an equivalent role as that of homologous recombination. Overall, nine reversal, nine translocation, one fusion, and two fission rearrangement events were detected in BILs. Out of which, only 72% reversal sites, 50% translocation sites, and 25% fission sites were generated by homologous recombination events (Supplementary Table S5). Previous studies in plants and humans have confirmed that microhomology-mediated NHEJ involves the occurrence of copy-number variation of genes in mitogenomes [23, 24]. In this study, 15 genes were identified with increased copy-number (Supplementary Table S4), which corresponds well with the global NHEJ in BIL lines . The findings indicate that both of the repeat-mediated recombination and non-homologous recombination contribute almost equal importantly to the differentiation of mitogenomes in rice.
In the O. glaberrima/93–11//93–11 BIL population, the molecular diversity was obviously more than that of O. glaberrima/93–11//O. glaberrima (Fig. 1), reflecting high consistency of mitogenome in most lines of O. glaberrima/93–11//O. glaberrima population. Although the mitochondrial genome is inherited independently to some extent, the DNA replication, RNA transcription, protein translation, and ATP production cannot work without precise cytonuclear co-adaptation . In the five tested BIL lines, the mitogenomes were differentiated at different levels (Fig. 2). The BIL line P91 bearing more 93–11 genomic fragments showed stronger mtDNA rearrangement than the BIL P10 which bears less 93–11 nuclear fragment (Supplementary Fig. S3), even though the relationship was not completely positive correlated. The probable justification is that the mitogenome needs to reprogram to accommodate the reshuffled nuclear genome when the nuclear genome is shocked by alien genome.
Dosage effect and RNA-editing play predominant roles in plant Mito-nuclear interaction
Natural distant hybridization is the main driving force for species formation  and interspecies hybridization usually leads to cytoplasmic male sterility because of species isolation [44,45,46]. The artificial wide crossing is a process to mimic the evolution in nature to select elite descendent with favorable agronomic traits. Here, the sterile lines in the BIL population were deleted after artificial selection (Supplementary Fig. S2). The analysis of the mitogenomes provided a medium to understand the evolutionary trends and adapt the nuclear genome reprogramming. Strikingly, global mtDNA rearrangement changed 15 gene copy-numbers (Table 1), and at least 20 gene clusters in the mitogenome (Fig. 2; Supplementary Table S6). These changes extensively modify the expression profile of mitochondria in rice (Fig. 4). The results demonstrated that the dosage effect plays an important role in regulating mito-nuclear interaction, although the exact mechanism needs to be investigated further.
Increasing reports have evidenced that RNA-editing is involved in regulating plant growth and development [47, 48]. It is considered that RNA editing implicates in the rescue of the organelle dysfunction caused by genetic, physiological or environmental factors . Interestingly, 529 mitogenomic RNA editing sites were identified in this research, of which, 143 editing sites were detected for the first time in rice (Table 2). It is worth noting that about 21.0% of the editing sites showed disturbance in their editing rates among BIL lines in relation to the maternal lines (Fig. 4b). The difference in the RNA editing site and the rate in BIL lines are strongly in agreement with the nuclear heterogeneity. To some extent, it reflects the co-adaption of the mitochondria and nuclear genomes, because the PPR protein encoded by the nuclear genes controls the RNA editing. It is a well-known fact that the PPR gene family is mainly responsible for the RNA-editing and rapidly gets expanded as well as diversified in the higher plants during the course of evolution [49, 50]. Considering the fact that animals have a few PPR members, lesser DNA rearrangement, and low RNA editing , it is proposed that mitochondria are functionally specified for energy production, which is a crucial role in organism development. In order to maintain a harmonious nuclear-cytoplasmic interaction, the PPR was actively evolved in plants to modify or trigger the organelle transcripts for adapting the mutated subunit of complex proteins encoded by nuclear genome.
In the present study, it was found that BIL mitogenomes transmitted from the maternal line were globally reprogrammed through DNA recombination in multiple ways. The global rearrangement changed genome structure, organization, and gene copies that altered the expression level of genes. Meanwhile, the RNA editing site and editing rate were also extensively changed, along with the dosage effect that triggers the interaction between nuclear and mitochondria to a concerted state.
Plant materials and culture
The 3 K rice lines were provided by Dr. Binying Fu, institute of crop science, Chinese academy of agricultural sciences, Beijing 100,081, China. O. glaberrima (ICR650) was provided by international rice research institute, Philippines. 93–11 was provided by Yangzhou academy of agricultural sciences, Yangzhou 225,007, China. Two BC2F12 BIL populations (ICR650/93–11//93–11 and ICR650/93–11//ICR650), which were built by us, derived from a wide cross of O. glaberrima × O. sativa (ICR650 × 93–11) and the maternal line ICR650 were used in this research. All the plants were planted in the camp of Wuhan University with standard management.
Detection of mitochondrial genotypes
Total genomic DNA was extracted from 2-week-old rice seedlings by the CTAB method. PCR amplification was performed using mitotype-specific molecular markers (MSS3, MSS5, MSS7, MSS8, MSS12, MSS13, MSS20, MSS21, MSS25, MSS27, MSS28, MSS29, MSS32, MSS33, MSS34, MSS37, MSS38, MSS39, MSS41, MSS42, MSS43, MSS44, MSS46, MSS47, MSS50, MSS51, MSS52, MSS53, MSS54, MSS55, MSS56 and MSS57)  in a 10 μL final reaction mixture containing 3.7 μL of distilled water, 5 μL 2 × Taq Plus Master Mix II (Vazyme, China), 0.4 μL of primer (10 pmole/μL) and 0.5 μL of genomic DNA (50 ng). All the reactions were set at an initial denaturation step (95 °C for 5 min), 30 cycles of denaturation (95 °C for 30 s), annealing (50 ~ 60 °C for 30 s), and extension (72 °C for 20 s ~ 70 s, depending on the size of the amplified fragment) and final extension (72 °C for 5 min). Amplification products were evaluated by running on 1.5% agarose gel. Each band produced by PCR was treated as a unit character and scored as a binary code (1/0). The phylogenetic tree was constructed by the neighbor-joining method, implemented in DPS (Data Processing System) (v17.10)  and MEGA X . The landscaping appearance of the tree was obtained by FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/).
Bin map construction
Total genomic DNA extracted from 2-week-old rice seedlings was sent to BENAGEN for sequencing on the Illumina HiSeq/MiSeq platform with the 400 bp paired-end library. Adapter and low-quality sequencing reads were trimmed by FastQC. Nipponbare genome (IRGSP-1.0) was downloaded from EnsemblPlants as a reference and all clean data were compared by BWA. GATK (Analysis Tools for Next-generation DNA Sequencers) v3.7.0 was used to call all SNPs. The maps of BILs were aligned and compared to their genotypes for a 15 SNP interval. Adjacent 15 SNP intervals across all BILs with the same genotype were combined into a recombination bin.
Mitochondrial genome sequencing, assembly, and annotation
Rice mitochondria were isolated from callus by differential centrifugation method followed by continuous Percoll gradients as described by Heazlewood et al. . The DNA was extracted from the rice mitochondria by CTAB and sequenced by NextOmics (Wuhan, China) on an Illumina HiSeq Xten machine with the 400 bp paired-end library, and on a PacBio RSII machine with 20 kb paired-end library.
The raw reads from PacBio RSII longer than 26 Kb were used as seed reads. The reads shorter than 26 Kb were corrected by the RS_PreAssembler.1 protocol with default settings from the Pacific Biosciences SMRT analysis (v2.3.0) software package. In addition, the raw reads from Illumina HiSeq Xten were also used to correct the genome by pilon (−-changes --vcf --fix bases --threads 5 --mindepth 10). Mitochondrial genomes were annotated using MITOFY . The tRNA genes were searched using tRNAscan-SE . Finally, the genome maps were drawn by OGDRAW .
93–11 special sequence detection
Mitochondrial genome (mitogenome) of O. glaberrima was aligned against 93–11 (NC_007886) using BLASTN  with a word size of 20 nucleotides and expectation value of 1. Specific sequences in the 93–11 mitogenome were extracted.
To identify mitochondrial synteny blocks between the pairs of representative BIL lines, all mitochondrial genomes were aligned using Mauve version 2.3.1  with an LCB cutoff of 500, respectively. Pairwise rearrangement distances in terms of minimum number of rearrangements were inferred using GRIMM .
Repeats and repeat-mediated recombination
Repeats were detected using BLASTN to search O. glaberrima mitogenome against itself with a word size of 50 nucleotides and an expectation value of 1. Each repeat copy was numbered and distinguished by the flanking sequences . The repeat sequences having only two copies were shorter than 1 Kb and not at the end of the scaffold were chosen to design primers (Supplementary information 1) by Primer 5  for detection.
The frequency of repeat-mediated recombination was evaluated for all repeats in mitogenomes that had only two copies and were not at ends of scaffold using procedures described by Mower et al. . The parameters and calculation methods reported by Guo et al.  were used for our analysis.
The content of one repeat configuration is equal to the mapping read pair number divided by all mapping read pairs of repeat configurations. The relative content of one repeat configuration of BIL lines is equal to the content of repeat configuration in BIL lines divided by the content of repeat configuration in O. glaberrima.
Adjacency of rearrangement breakpoints and repeats
HRBs (homologous-region boundaries) verified by pairwise Mauve alignments were considered as agencies for rearrangement sites following the method of Cole, et al. . In brief, HRBs, which were at the extremity of scaffolds, were eliminated from our analyses. The proportion of HRBs that was close to repeats was equal to the number of HRBs within 50 bp of repeats divided by the total number of HRBs.
RNA extraction and transcriptome analysis
Total RNA was isolated using TRIzol reagent (Invitrogen) from the 2-week-old rice seedlings. Ribosomal RNA content was reduced from total RNA following the method of Guo, et al. . The Ribo-Zero™-treated organellar RNAs of all samples were sent to BENAGEN for 75 bp single-read sequencing on the Illumina HiSeq/MiSeq platform.
Adapter and low-quality bases of RNA sequencing data were trimmed according to Guo et al. . Clean reads were mapped to mitochondrial genome by bowtie2 v18.104.22.168  (−-no-mixed--no-discordant --gbar 1000 --end-to-end -k 200 -q -X 800). RSEM v1.2.25  was used to calculate FPKM of the gene.
Reverse transcription, fluorescence quantification, and RNA-editing verification
Approximately 4 μg of RNA was treated with DNase I (NEB) and cDNA was obtained by reverse transcription with SuperScriptII as described in the manufacturer’s instructions. Real-time quantitative PCR was performed using LightCycler 480 (Roche) and the SYBR Green I Master PCR kit (Roche).
RNA-editing sites were identified following the method of Guo et al. . The editing sites with a cover depth lesser than 50 were removed by SAM tools mpileup . The cDNA containing the editing site was ligated into pMD18-T and the constructs were co-transformed into E. coli. For Sanger sequencing, 50 monoclonal constructs of every editing site were taken and the RNA editing rate was calculated.
Availability of data and materials
The data of all DNA-Seq and RNA-Seq has being uploaded to the NCBI SRA database (http://www.ncbi.nlm.nih.gov/sra), with the accession number: PRJNA598996, Other data that supporting the conclusions of the article have been uploaded as additional files. Seeds are available upon request from the corresponding author.
Backcross inbred lines
Non-homologous end joining
Burton RS, Pereira RJ, Barreto FS. Cytonuclear genomic interactions and hybrid breakdown. Annu Rev Ecol Evol Syst. 2013;44(1):281–302.
Budar F, Roux F. The role of organelle genomes in plant adaptation. Plant Signal Behav. 2014;6(5):635–9.
Roux F, Mary-Huard T, Barillot E, Wenes E, Botran L, Durand S, Villoutreix R, Martin-Magniette ML, Camilleri C, Budar F. Cytonuclear interactions affect adaptive traits of the annual plant Arabidopsis thaliana in the field. Proc Natl Acad Sci U S A. 2016;113(13):3687–92.
Greiner S, Bock R. Tuning a menage a trois: co-evolution and co-adaptation of nuclear and organellar genomes in plants. BioEssays. 2013;35(4):354–65.
Greiner S, Sobanski J, Bock R. Why are most organelle genomes transmitted maternally? BioEssays. 2015;37(1):80–94.
Sato M, Sato K. Maternal inheritance of mitochondrial DNA by diverse mechanisms to eliminate paternal mitochondrial DNA. Biochim Biophys Acta. 2013;1833(8):1979–84.
Barr CM, Neiman M, Taylor DR. Inheritance and recombination of mitochondrial genomes in plants, fungi and animals. New Phytol. 2005;168(1):39–50.
McCauley DE. Paternal leakage, heteroplasmy, and the evolution of plant mitochondrial genomes. New Phytol. 2013;200(4):966–77.
Wicke S, Schneeweiss GM, De Pamphilis CW, Muller KF, Quandt D. The evolution of the plastid chromosome in land plants: gene content, gene order, gene function. Plant Mol Biol. 2011;76(3–5):273–97.
Ward B, Anderson R, Bendich A. The mitochondrial genome is large and variable in a family of plants (cucurbitaceae). Cell. 1981;25:793–803.
Elansary HO, Müller K, Olson MS, Štorchová H. Transcription profiles of mitochondrial genes correlate with mitochondrial DNA haplotypes in a natural population of Silene vulgaris. BMC Plant Biol. 2010;10(1):11.
Wang S, Li D, Yao X, Song Q, Wang Z, Zhang Q, Zhong C, Liu Y, Huang H. Evolution and diversification of kiwifruit mitogenomes through extensive whole-genome rearrangement and mosaic loss of intergenic sequences in a highly variable region. Genome Biol Evol. 2019;11(4):1192–206.
Woloszynska M. Heteroplasmy and stoichiometric complexity of plant mitochondrial genomes--though this be madness, yet there's method in't. J Exp Bot. 2010;61(3):657–71.
Gao LZ, Innan H. Nonindependent domestication of the two rice subspecies, Oryza sativa ssp. indica and ssp. japonica, demonstrated by multilocus microsatellites. Genetics. 2008;179(2):965–76.
Huang X, Kurata N, Wei X, Wang ZX, Wang A, Zhao Q, Zhao Y, Liu K, Lu H, Li W, et al. A map of rice genome variation reveals the origin of cultivated rice. Nature. 2012;490(7421):497–501.
Londo JP, Chiang YC, Hung KH, Chiang TY, Schaal BA. Phylogeography of Asian wild rice, Oryza rufipogon, reveals multiple independent domestications of cultivated rice, Oryza sativa. Proc Natl Acad Sci U S A. 2006;103(25):9578–83.
Civáň P, Craig H, Cox CJ, Brown TA. Three geographically separate domestications of Asian rice. Nat Plants. 2015;1(11):15164.
Wang W, Mauleon R, Hu Z, Chebotarov D, Tai S, Wu Z, Li M, Zheng T, Fuentes RR, Zhang F, et al. Genomic variation in 3,010 diverse accessions of Asian cultivated rice. Nature. 2018;557(7703):43–9.
Xie H, Wang J, Qian M, Li N, Zhu Y, Li S. Mitotype-specific sequences related to cytoplasmic male sterility in Oryza species. Mol Breed. 2013;33(4):803–11.
Tesler G. GRIMM: genome rearrangements web server. Bioinformatics. 2002;18(3):492–3.
Cole LW, Guo W, Mower JP, Palmer JD. High and variable rates of repeat-mediated mitochondrial genome rearrangement in a genus of plants. Mol Biol Evol. 2018;35(11):2773–85.
Arrieta-Montiel MP, Shedge V, Davila J, Christensen AC, Mackenzie SA. Diversity of the Arabidopsis mitochondrial genome occurs via nuclear-controlled recombination activity. Genetics. 2009;183(4):1261–8.
Tadi SK, Sebastian R, Dahal S, Babu RK, Choudhary B, Raghavan SC. Microhomology-mediated end joining is the principal mediator of double-strand break repair during mitochondrial DNA lesions. Mol Biol Cell. 2016;27(2):223–35.
Hastings PJ, Ira G, Lupski JR. A microhomology-mediated break-induced replication model for the origin of human copy number variation. PLoS Genet. 2009;5(1):e1000327.
Chateigner-Boutin AL, Small I. Organellar RNA editing. Wiley Interdiscip Rev RNA. 2011;2(4):493–506.
Luan J, Liu T, Luo W, Liu W, Peng M, Li W, Dai X, Liang M, Chen L. Mitochondrial DNA genetic polymorphism in thirteen rice cytoplasmic male sterile lines. Plant Cell Rep. 2013;32(4):545–54.
Song Z, Lu B, Wang B, Chen J. Fitness estimation through performance comparison of F1 hybrids with their parental species Oryza rufipogon and O. sativa. Ann Bot. 2004;93:311–6.
Fan F, Li N, Chen Y, Liu X, Sun H, Wang J, He G, Zhu Y, Li S. Development of elite BPH-resistant wide-spectrum restorer lines for three and two line hybrid rice. Front Plant Sci. 2017;8:986.
Li ZK, Sanchez A, Angeles E, Singh S, Domingo J, Huang N, Khush GS. Are the dominant and recessive plant disease resistance genes similar? A case study of rice R genes and Xanthomonas oryzae pv. oryzae races. Genetics. 2001;159(2):757–65.
Tanksley SD, McCouch SR. Seed banks and molecular maps: unlocking genetic potential from the wild. Science. 1997;277(5329):1063–6.
Li SQ, Yang DC, Zhu YG. Characterization and use of male sterility in hybrid rice breeding. J Integr Plant Biol. 2007;49(6):791–804.
Meyer VG. A study of reciprocal hybrids between upland cotton (Gossypium hirsutum L.) and experimental lines with cytoplasms from seven other species. Crop Sci. 1973;13(4):439–44.
Hu Q, Hansen N, Laursen J, Dixelius C, Andersen B. Intergeneric hybrids between Brassica napus and Orychophragmus violaceus containing traits of agronomic importance for oilseed rape breeding. Theor Appl Genet. 2002;105(6–7):834–40.
Price HJ, Hodnett GL, Burson BL, Dillon SL, Rooney WL. A Sorghum bicolor × S. macrospermum hybrid recovered by embryo rescue and culture. Aust J Botany. 2005;53(6):579.
Sloan DB, Warren JM, Williams AM, Wu Z, Abdel-Ghany SE, Chicco AJ, Havird JC. Cytonuclear integration and co-evolution. Nat Rev Genet. 2018;19(10):635–48.
Gualberto JM, Newton KJ. Plant mitochondrial genomes: dynamics and mechanisms of mutation. Annu Rev Plant Biol. 2017;68:225–52.
Arrieta-Montiel MP. Mackenzie SA. Plant mitochondrial genomes and recombination. Advances in Plant Biology. 2011:65–82.
Stern D, Palmer J. Extensive and widespread homologies between mitochondrial DNA and chloroplast DNA in plants. Proc Natl Acad Sci U S A. 1984;81(7):1946–50.
Palmer JD, Shields CR. Tripartite structure of the Brassica campestris mitochondrial genome. Nature. 1984;307:437–40.
Lonsdale DM, Hodge TP, Fauron CMR. The physical map and organisation of the mitochondrial genome from the fertile cytoplasm of maize. Nucleic Acids Res. 1984;12:9249–61.
Weng ML, Blazier JC, Govindu M, Jansen RK. Reconstruction of the ancestral plastid genome in Geraniaceae reveals a correlation between genome rearrangements, repeats, and nucleotide substitution rates. Mol Biol Evol. 2014;31(3):645–59.
Seol JH, Shim EY, Lee SE. Microhomology-mediated end joining: good, bad and ugly. Mutat Res. 2018;809:81–7.
Schumer M, Rosenthal GG, Andolfatto P. How common is homoploid hybrid speciation? Evolution. 2014;68(6):1553–60.
Budar F, Touzet P, Pelletier G. Cytoplasmic male sterility. Annual Plant Reviews Online. 2018:147–80.
Kamiński P, Podwyszyńska M, Starzycki M, Starzycka-Korbas E. Interspecific hybridisation of cytoplasmic male-sterile rapeseed with Ogura cytoplasm and Brassica rapa var. pekinensis as a method to obtain male-sterile Chinese cabbage inbred lines. Euphytica. 2015;208(3):519–34.
Chen G, Ye X, Zeng F, Wang J, Yuan L, Zhu S, Hou J, Cheng Y, Wang C. Characterization and utilization of a cytoplasmic male sterility line of Wucai (Brassica campestris L.). Hortic. Environ. Biotechnol. 2019;60(3):373–82.
Xiao H, Xu Y, Ni C, Zhang Q, Zhong F, Huang J, Liu W, Peng L, Zhu Y, Hu J. A rice dual-localized pentatricopeptide repeat protein is involved in organellar RNA editing together with OsMORFs. J Exp Bot. 2018;69(12):2923–36.
Xiao H, Zhang Q, Qin X, Xu Y, Ni C, Huang J, Zhu L, Zhong F, Liu W, Yao G, et al. Rice PPS1 encodes a DYW motif-containing pentatricopeptide repeat protein required for five consecutive RNA-editing sites of nad3 in mitochondria. New Phytol. 2018;220(3):878–92.
O'Toole N, Hattori M, Andres C, Iida K, Lurin C, Schmitz-Linneweber C, Sugita M, Small I. On the expansion of the pentatricopeptide repeat gene family in plants. Mol Biol Evol. 2008;25(6):1120–8.
Barkan A, Small I. Pentatricopeptide repeat proteins in plants. Annu Rev Plant Biol. 2014;65:415–42.
Tang QY, Zhang CX. Data processing system (DPS) software with experimental design, statistical analysis and data mining developed for use in entomological research. Insect Sci. 2013;20(2):254–60.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9.
Heazlewood JL, Howell KA, Whelan J, Millar AH. Towards an analysis of the rice mitochondrial proteome. Plant Physiol. 2003;132(1):230–42.
Alverson AJ, Wei X, Rice DW, Stern DB, Barry K, Palmer JD. Insights into the evolution of mitochondrial genome size from complete sequences of Citrullus lanatus and Cucurbita pepo (Cucurbitaceae). Mol Biol Evol. 2010;27(6):1436–48.
Schattner P, Brooks AN, Lowe TM. The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res. 2005;33(Web Server issue):W686–9.
Lohse M, Drechsel O, Kahlau S, Bock R. OrganellarGenomeDRAW--a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucleic Acids Res. 2013;41(W1):W575–81.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Darling AE, Mau B, Perna NT. ProgressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One. 2010;5(6):e11147.
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG. Primer3--new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):e115.
Mower JP, Case AL, Floro ER, Willis JH. Evidence against equimolarity of large repeat arrangements and a predominant master circle structure of the mitochondrial genome from a monkeyflower (Mimulus guttatus) lineage with cryptic CMS. Genome Biol Evol. 2012;4(5):670–86.
Guo W, Grewe F, Fan W, Young GJ, Knoop V, Palmer JD, Mower JP. Ginkgo and Welwitschia mitogenomes reveal extreme contrasts in gymnosperm mitochondrial evolution. Mol Biol Evol. 2016;33(6):1448–60.
Guo W, Grewe F, Mower JP. Variable frequency of plastid RNA editing among ferns and repeated loss of uridine-to-cytidine editing from vascular plants. PLoS One. 2015;10(1):e0117075.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. Genome project data processing S. the sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
We thank Dr. Binying Fu for his kind providing the 3 K rice lines used in this paper.
This work was partly granted from the National Key Research and Development Program (2016YFD0100903), the National Transgenic Research Program (2016ZX08001004–001-002), the National Science Foundation (31571630) of China and Huanghe Patent Program of Wuhan City. The funders had no role in the study design, data analysis and interpretation, and manuscript writing, but just provided the financial.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Polymorphism of mitochondrial genomes in landraces was detected using 32 mitochondria-specific molecular markers. The green represents Oryza sativa japonica. The blue represents Oryza sativa indica.
Gross plant morphology of the maternal parent and BILs. A. Morphologies of the maternal line and BILs. Scale bars, 10 cm. B. Basic agronomic traits of BILs.
The bin map of five BIL lines. The green represents the genome fragment which from O. glaberrima. The red represents the genome fragment which from 93-11. The blue represents the genome fragment which is heterozygotic. 15 SNPs for a window.
Validation of the mitochondrial genome of O. glaberrima.. A. Comparison between Nipponbare mitochondrial genome and O. glaberrima mitochondrial genome. B. Verification of different structural sites.
Complete mitogenomes of the maternal parent O. glaberrima and five BIL lines.
Percentage of different configurations of mitochondrial two-copy repeats in BIL lines compared to the O. glaberrima.
Relative content of different configurations of some mitochondrial repeats in BIL lines compared to O. glaberrima by fluorescence quantification.
The change of gene order based on homologous recombination.
DNA content and RNA levels of some mitochondrial genes in BIL lines compared to O. glaberrima by fluorescence quantification. A. Relative DNA and RNA level to O. glaberrima of mitochondrial genes near repeat within 2 Kb. B. Relative DNA and RNA level to O. glaberrima of mitochondrial genes far away from repeat over 2 Kb.
Summary of sequencing nuclear genomes by Illumina. Table S2. Summary of sequencing mitogenomes by Illumina. Table S3. Summary of sequencing mitogenomes by PacBio. Table S4. Multi-copy genes in mitogenomes of the BIL lines. Table S5. Rearrangement types and recombination manners in BIL mitogenomes. Table S6. Variation of some representative gene clusters in mitogenome of BIL lines. Table S7. Special mitogenome sequences in 93–11. Table S8. Homologous region boundaries of BIL lines relative to O. glaberrima. Table S9. All repeats detected in O. glaberrima mitogenome. Table S10. Category of multi-copy repeats in O. glaberrima mitogenome. Table S11. RNA editing events detected in mitochondria of BIL lines. Table S12. Confirmation of RNA editing using Sanger sequencing
About this article
Cite this article
Yang, W., Zou, J., Wang, J. et al. Wide crossing diversify mitogenomes of rice. BMC Plant Biol 20, 159 (2020). https://doi.org/10.1186/s12870-020-02380-w
- Genomic rearrangement
- Maternal inheritance
- Nucleo-cytoplasmic interaction