- Open Access
Characterization of the mitochondrial genome of the MAX1 type of cytoplasmic male-sterile sunflower
BMC Plant Biologyvolume 19, Article number: 51 (2019)
More than 70 cytoplasmic male sterility (CMS) types have been identified in Helianthus, but only for less than half of them, research of mitochondrial organization has been conducted. Moreover, complete mitochondrion sequences have only been published for two CMS sources – PET1 and PET2. It has been demonstrated that other sunflower CMS sources like MAX1, significantly differ from the PET1 and PET2 types. However, possible molecular causes for the CMS induction by MAX1 have not yet been proposed. In the present study, we have investigated structural changes in the mitochondrial genome of HA89 (MAX1) CMS sunflower line in comparison to the fertile mitochondrial genome.
Eight significant major reorganization events have been determined in HA89 (MAX1) mtDNA: one 110 kb inverted region, four deletions of 439 bp, 978 bp, 3183 bp and 14,296 bp, respectively, and three insertions of 1999 bp, 5272 bp and 6583 bp. The rearrangements have led to functional changes in the mitochondrial genome of HA89 (MAX1) resulting in the complete elimination of orf777 and the appearance of new ORFs - orf306, orf480, orf645 and orf1287. Aligning the mtDNA of the CMS sources PET1 and PET2 with MAX1 we found some common reorganization features in their mitochondrial genome sequences.
The new open reading frame orf1287, representing a chimeric atp6 gene, may play a key role in MAX1 CMS phenotype formation in sunflower, while the contribution of other mitochondrial reorganizations seems to appear negligible for the CMS development.
Widespread use of hybrid seed production for cross-pollinated crop plants makes cytoplasmic male sterility (CMS) one of the most desired agricultural traits . More than 70 cytoplasmic male sterility types have been identified in Helianthus , but almost all sunflower hybrids are currently based on a single source of CMS – PET1, discovered 50 years ago by Leclercq  in an interspecific hybrid between Helianthus petiolaris Nutt and the cultivated sunflower - H. annuus L. Such genetic homogeneity can result in decrease of yield and disease resistance [4, 5]. However, potential new sources of CMS (CMS cytotypes) have not been sufficiently investigated at the molecular level, presenting a serious obstacle for their introduction into commercial hybrid breeding .
The CMS phenomenon provides a convenient model to study nuclear-cytoplasmic interaction and its role in pollen development [7,8,9]. Nuclear male sterility is controlled only by nuclear genes, whereas cytoplasmic sterility occurs as a result of the interaction of a male sterile cytoplasm with nuclear sterility factors [10, 11]. Nuclear control of CMS is accomplished by Rf genes performing fertility restoration .
Genetic basis of cytoplasmic male sterility is formed by structural rearrangements of mitochondrial DNA (mtDNA) [12,13,14]. New chimeric ORFs in the mtDNA of CMS lines are formed from fusion of known functional mitochondrial genes with uncharacteristic DNA sequences [1, 12, 15]. For instance, a new open reading frame (orfH522) appeared in the 3′-flanking region of the atp1 gene encoding the alpha subunit of mitochondrial F1 ATPase in the CMS PET1 of sunflower as a result of an 11-kb-inversion and a 5-kb-insertion [16, 17]. This orfH522 is co-transcribed with the atp1 gene as a polycistronic mRNA  and translated into the CMS-specific 16-kDa-protein . Molecular basis of another CMS type of sunflower - PET2 - has been studied recently [20, 21]. Both CMS types (PET1, PET2) were obtained from interspecific crosses between H. petiolaris and the cultivated sunflower H. annuus [2, 22]. However, structural rearrangements of the PET2 mtDNA leading to the emergence of new open reading frames differ from those in CMS-PET1. Multiple rearrangements have been detected in the mitochondrial sequence of PET2: two translocations (27.5 kb and 106.5 kb), two deletions (711 bp and 3780 bp), as well as insertions of 5050 bp and 15,885 bp . Some of these mitochondrial genome reorganizations resulted in formation of novel ORFs (orf645, orf2565, orf228, and orf285) with transcriptional activity specific to the CMS PET2. It was suggested that orf228 (orf231 including the stop codon) and orf285 (orf288 including the stop codon), being chimeric atp9 ORFs, may be the main driver for the development of PET2 CMS phenotype [20, 21]. Reduction of the co-transcript of the two ORFs could be specifically detected in the anthers of fertility-restored hybrids, supporting their role in the CMS phenotype .
Except for PET1 and PET2, the other CMS types of sunflower are poorly investigated. To the best of our knowledge, only one comparative study of mitochondrial genomes organizations including 28 CMS types of sunflower has been conducted . There were also several independent attempts to analyze the mtDNA of MAX1 [24, 25] and PEF1 [26, 27] CMS types. The investigations demonstrated that MAX1 CMS source significantly differs from the PET1 and PET2 mitotypes [23,24,25]. Thus MAX1 CMS source, which was initially obtained by the interspecific hybridization of an annual sunflower (Helianthus annuus L.) with the perennial species Helianthus maximiliani Schrad , is insufficiently explored and is of particular interest.
In the present study, we have investigated structural changes in the mitochondrial genome of HA89 (MAX1) CMS sunflower line and identified possible causes for the development of male sterility in the MAX1 mitotype.
Materials and methods
The CMS line HA 89 (MAX1) as well as the fertile line HA 89 with the normal cytoplasm were obtained from the genetic collection of the N. I. Vavilov Institute of Plant Genetic Resources (VIR, Russia). Sunflower plants were grown in regularly irrigated pots in the growth chamber KBWF 720 (Binder, Germany) at the following growth conditions: 26 °C (78.8 F) temperature, 70% humidity, and 10/14 h dark/light cycle. For the DNA and RNA manipulations, we used leaf tissue from five plants, which were on the same growth stage and had maximal phenotypic similarity.
Mitochondrial DNA extraction, genome library construction, and NGS sequencing
We extracted the organelle fraction with a reduced amount of nuclear DNA from leaves of 14-days-old sunflower seedlings as described by Makarenko et al. . DNA isolation was performed with the PhytoSorb kit (Syntol, Russia), according to the manufacturer’s protocol. The NGS libraries were prepared with Nextera XT DNA Library Prep Kit (Illumina, USA), following the sample preparation protocol by Illumina. The concentration of the prepared library was measured with the Qubit fluorometer (Invitrogen, USA) and qPCR. Fragment length distribution was determined with Bioanalyzer 2100 (Agilent, USA). Libraries for NGS sequencing were diluted up to the concentration of 10 pM and sequenced on NextSeq 500 sequencer using High Output v2 kit (Illumina, USA). A total number of 10,174,497,150-bp paired reads were generated.
Mitogenome assembly and annotation
Quality of reads was determined with FastQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Trimming of adapter-derived and low quality (Q-score below 25) reads was performed with the Trimmomatic software . The whole genome assembly was done using the SPAdes genome assembler v 3.10.1 . High coverage (depths > 100) contigs were then manually assembled based on the scaffolds and presence of bridge contigs. Validation was performed by mapping the initial reads on the assembled genome sequence using bowtie 2 . The assembled genome was aligned to the reference sunflower mitochondrial genome (NC_023337.1) using BLAST. Verification of the discovered rearrangements was made by PCR. Primer sequences are presented in Additional file 1. The mitochondrial genome was annotated with glimmer 3.02b  and BLAST tools. Potential ORFs were identified using ORFfinder (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/orffinder). Graphical genome map was generated using the OGDRAW tool . Transmembrane domains were predicted using the TMHMM server v.2.0 (available online: http://www.cbs.dtu.dk/services/TMHMM-2.0/). We have conducted analysis of distribution of transcription factor binding sites (TFBS) using Match  and cisExpress [36, 37]. The complete mitochondrion sequence of CMS line HA89(MAX1) has been deposited to GenBank under the accession number MH704580.1.
Comparison of HA89 (MAX1) was with HA 89 fertile line genome (available at https://doi.org/10.6084/m9.figshare.7265648.v1); this sequence will later be deposited to NCBI GenBank. In our analysis we used positions of a publicly available mitochondrial genome of a fertile line that differs from HA89 by two nucleotides - HA412 (NC_023337.1).
RNA extraction and qRT-PCR
Total RNA from leaves was extracted with guanidinium thiocyanate-phenol-chloroform reagent kit – ExtractRNA (Evrogen, Russia). RNA quality and concentration were measured using the NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, USA) and the Qubit fluorometer (Invitrogen, USA). 0.5 μg of total RNA was treated with DNAse I (Thermo Fisher Scientific, USA) per manufacturer’s instruction. The cDNA was synthesized using random primers and MMLV RT kit (Evrogen, Russia). The qPCR was performed with PCR kit with EvaGreen dye (Syntol, Russia) on Rotor-Gene 6000 (Corbett Research, Australia). Primer sequences are presented in Additional file 1.
We assembled the complete mitochondrial genome of HA89 sunflower line with MAX1 cytoplasmic sterility type. The master chromosome of HA89 (MAX1) consists of 295,586 bp (Fig. 1). HA89 (MAX1) mitochondrial genome is 5361 bp shorter than the mtDNA genome of HA89 fertile analogue. The HA89 (MAX1) mitochondrial genome showed several rearrangements as compared to the HA89 fertile genome. We have identified eight significant reorganization events of HA89 (MAX1) mtDNA: one large inverted region, four deletions, and three insertions.
Comparative analysis of HA89 (MAX1) and its fertile analogue identified one inverted region, which is approximately 110 kbp long. This rearrangement includes a region flanked by orf873 and atp8 and is highlighted in Fig. 1. We like to emphasize that this rearrangement covers more than 30% of the entire mitochondrial genome. Presence of this huge rearrangement is also associated with most of the identified deletions and insertions. We have found four long deletions (more than 100 bp) in HA89 (MAX1) as compared with the HA89 fertile line: 439 bp, 978 bp, 3183 bp and 14,296 bp. All deletions except for the smallest one (439 bp) are associated with this 110 kb region. The 978 bp deletion is flanking the 3’region of the MAX1 specific 110 kb rearrangement. Matching this deletion to the sunflower fertile line mitochondrial genome map, we determined that it was located between the genes nad4L and atp8 (35,559–36,537 positions in the HA89 fertile mitochondrion). This area of the mitochondrial genome normally codes for orf777. Thus the 978 bp deletion results in the absence of orf777 in the HA89 (MAX1) mitochondrial genome. The other two deletions (3183 and 14,296 bp) also occurred within this 110 kb region of the MAX1 mitochondrial genome. Compared to the HA89 fertile line, both deletions are in intergenic regions and do not directly impact coding sequences. The 3183 bp deletion is localized between atp6-coxII (272,051–275,233 positions in the HA89 fertile mtDNA). The 14,296 bp deletion predominantly affects the large repeat region of the sunflower mtDNA - coxII-nad2 (286,492–300,787 positions in the HA89 fertile mitochondrion). The smallest deletion (439 bp) also belongs to a large repeat intergenic region of mtDNA – nad4-ccmB (56,166–56,604 positions in the HA89 fertile mtDNA).
We have also detected several insertions in the HA89 (MAX1) CMS line. Three long (> 100 bp) insertions were found with sizes of 1999, 5272 and 6583 bp, respectively. The smallest one (1999 bp) is in the intergenic region nad4-ccmB, co-located with the 439 bp deletion. Thus, the occurrence of the 1999 bp insertion and the deletion of the mtDNA large repeat region (14,296 bp deletion) could be connected. On one hand, insertion in one copy of mtDNA repeat region could impair recombination of mtDNA molecules and eventually lead to the elimination of the other repeat copy in the mitochondrial genome. On the other hand, initial deletion of one copy of the repeat region can make the remaining copy more prone to the structural changes. The 1999 bp insertion has 100% similarity to the inverted repeat region of the sunflower chloroplast DNA. This insertion does not lead to the formation of new open reading frames; however, it includes a part of the rrn23 gene. The 5272 bp insertion is co-localized with the other deletion site (3183 bp) - atp6-coxII and makes up a part of the rearrangement encompassing the 110 kb inverted region. The 5272 bp insertion results in a single large (more than 300 bp) new open reading frame – orf645. We detected transcriptional activity for orf645 in MAX1 CMS line by RT-qPCR. The largest insertion (6583 bp) identified in the HA89 (MAX1) genome is localized in the 5′ flanking region of the 110 kb rearrangement between orf873 and nad4L (Fig. 1). The 6583 bp insertion includes 1233 bp 99% identical to another part of the mitochondrial genome (205,493–206,725 positions of MAX1 mtDNA), creating a duplication of this mitochondrial region. The duplicated region includes a partial sequence of the atp6 gene. The 6583 bp insertion results in the chimeric open reading frame orf1287. The orf1287 is a fusion of 530 bp (in the 5′-end) of an unknown sequence with 760 bp (in the 3′-end) identical to atp6. This orf1287 would encode an elongated protein of 429 amino acids with a partial homology to ATP6. The chimeric atp6 transcript of the orf1287 was detected in the male sterile HA89(MAX1) line by qPCR, but not in the fertile line. In addition to the orf1287, we predicted in the 6583 bp insertion two more potential ORFs - orf306 and orf480. The search for transmembrane domains in these ORFs showed no transmembrane helix in orf306, but one possible membrane related domain in orf480 (prediction chance was about 50%) and seven transmembrane helices in orf1287 (Fig. 2).
The analysis of transcription factor binding sites (TFBS) revealed significant differences in distribution of TBFS in the promoter region of atp8. HA89 (MAX1) CMS promoter is 166 bp long, whereas HA 89 fertile line promoter is considerably longer, 384 bp. Earlier studies in radish had shown that in the male sterile mitochondrial genome the promoter of the atp8 gene was lost . The same study demonstrated that in CMS radish, atp8 was co-transcribed with orf138 and trnfM, but transcribed from its own promoter in a fertile variety. In the CMS MAX1 mitochondrial genome the genes atp1, atp8 and coxIII are located next to each other (Fig. 1), hinting at a possibility of co-transcription using the promoter of atp1. The promoter sequences of atp1 are not conserved between HA 89 (MAX1) and HA 89 fertile lines.
Although more than 70 types of cytoplasmic male sterility have been identified in sunflower , but research on their mitochondrial organization has been conducted only for less than half of them [23,24,25,26, 39, 40]. Assembly of plant mitochondrial genomes is always a difficult task . It is especially hard to assemble a single mitochondrial DNA molecule, the so-called “master chromosome”. Most problems are associated with dynamic rates of rearrangements in plant mitochondrial genomes, resulting in the formation of a variety of subgenomic forms [42, 43]. Complete mitochondrial sequences have been investigated only for two CMS sources in sunflower . Comparing the current research data with our previous study , where we analyzed mtDNA of PET1 (MG735191.1) and PET2 (MG770607.2) CMS types of sunflower, we can formulate a hypothesis about the molecular basis of CMS.
Matching the mtDNA of PET1 and MAX1 CMS sources, we found common mitochondrial rearrangement features. Both lines have a similar “rearrangement hot spot” near the orf873 gene. A part of the 6583 bp insertion (MAX1) is complimentary (with 99% identity) to the 4732 bp insertion (PET1). Such sequence homology leads to the presence of orf306 in both lines and orf480 with partial similarity to orfH522 in MAX1. Comparison of PET2 and MAX1 mtDNA also showed similarities in reorganizations of their mitochondrial genomes. Both CMS types have deletions in the nad4L-atp8 region and are characterized by the total elimination of orf777. About 80% of the nucleotide sequence of 5272 bp insertion (MAX1) is reverse-complement (with 99% identity) to the 5050 bp insertion (PET2). Moreover, orf645 is found in both CMS lines. The 6583 bp insertion (MAX1) has similarity to the 15,885 bp insertion in PET2. Both insertions have a common region counting about 500 nucleotides and the 1233 bp duplication site (MAX1) is also identical in the PET2 mitochondrial genome. The provided data in the current study confirm the similarity of MAX1 and PET2 mtDNAs organization, which had been observed in former hybridization results . However, the fertility restoration patterns for three CMS types (MAX1, PET1, PET2) differed considerably. For CMS MAX1 all investigated maintainer and restorer lines of PET1, apart from HA 89, restored fertility. But only two maintainers and one restorer line of CMS PET1 were able to restore the fertility of plants with PET2 CMS . It indicates a different origin of male sterility in these three CMS types, which is confirmed by the molecular comparison of the CMS types presented in the current research.
CMS phenotype formation predominantly involves the appearance of new open reading frames with partial complementary to mitochondrial genes: ATPase complex subunits, respiratory-related proteins, etc. [12, 45]. Mitochondrial DNA deletions usually lead to male sterility, when they result in the creation of new ORFs or the impairing coding sequence of vital mitochondrial proteins . Elimination of the entire orf777, encoding a protein with an unknown function, is unlikely to result in male sterility, despite the absence of the orf777 in both CMS types of sunflower (PET2 and MAX1). Most CMS associated genes encode polypeptides with transmembrane domains [11, 46]. Among the discovered new ORFs are orf306 and orf645, encoding proteins lacking transmembrane regions; therefore, it is unlikely that they might cause CMS by disturbing membrane integrity. The appearance of orf645 in two CMS types (PET2 and MAX1) of sunflower is probably not accidental, and orf645 could also be a conserved ORF in mitochondrial genomes of other Helianthus genus species. However, orf480 has a higher likelihood to result in the appearance of the CMS phenotype. The orf480 translated protein has one membrane related region and shows similarity to the ORFH522 protein, which plays a major role in formation of the PET1 CMS type . However, in the case of the ORFH522 protein, there is a definite transmembrane helix that is formed by the N-terminal 22 amino acids (20–42 aa). In turn, the orf480 translated polypeptide has only a membrane related region (predominantly inside the membrane’s layer) in the N-terminus (5–27 aa) and is unlikely to form a transmembrane helix. This orf480, also called orfM160 in a prior MAX1 study , was not transcribed and was also present in the male fertile H. maximiliani population MAX30. That is why we also do not consider orf480 to be the leading cause of the MAX1 CMS type in sunflower.
The chimeric orf1287 identified in the HA89 (MAX1) mitochondrial genome, is a potential candidate gene for the CMS development. Neither the HA89 fertile line nor the CMS lines HA89 (PET1), HA89 (PET2) has mitochondrial genes homologous to orf1287. The orf1287 encodes the atp6 chimeric protein, which has 250 aa common to ATP synthase Fo subunit 6 C-terminus but has an additional extended and dissimilar N-terminal part (Fig. 2). All seven transmembrane helices present in the C-terminus of a normal ATP6 protein also exist in the protein encoded by orf1287. We searched for homologous proteins in the NCBI database and found a significant identity (more than 70%) to ATP synthase subunit 6 of sunflower line with another CMS type - ANT1 (CAA57790.1). Compared with a normal ATP synthase subunit 6 (351 aa long) in both types of CMS (MAX1, ANT1) the elongated proteins would be translated into proteins of 429 aa (MAX1) and 437 aa (ANT1), respectively. It should be noted that they have identical C-termini and highly similar N-termini. The atp6 gene is one of the most frequently rearranged mitochondrial genes in plants [12, 47, 48]. Correlation between the chimeric atp6 gene and CMS phenotypes have been found in several species, including pepper , maize , radish , monkeyflower , etc. It is interesting to note, that in sugar beet the chimeric ORF presented by the atp6 5′ elongated transcript has been translated in vivo and caused CMS phenotype . From these observations, we suggest that the predicted translation product of orf1287 may play an important role in MAX1 CMS type formation in sunflower.
Eight significant major reorganization events of HA89 (MAX1) mtDNA have been determined: one 110 kb encompassing inverted region, four deletions – 439 bp, 978 bp, 3183 bp and 14,296 bp, and three insertions - 1999 bp, 5272 bp and 6583 bp. The rearrangements have led to functional changes in mitochondrial genome of HA89 (MAX1): complete elimination of orf777 and appearance of new ORFs – orf306, orf480, orf645 and orf1287. The predicted translation product of orf1287 may play a key role in the development of the male sterile phenotype in CMS MAX1 in sunflower, while the contributions of other mtDNA reorganizations to the development of CMS phenotype appear to be unlikely.
Bohra A, Jha UC, Adhimoolam P, Bisht D, Singh NP. Cytoplasmic male sterility (CMS) in hybrid breeding in field crops. Plant Cell Rep. 2016;35:967–93.
Serieys HA. Identification, study and utlization in breeding programs of new CMS. Helia. 1996;19:144–60.
Leclercq P. Une sterilite male chez le tournesol. Ann Amelior Plant. 1969;19:99–106.
Levings CS 3rd. The Texas cytoplasm of maize: cytoplasmic male sterility and disease susceptibility. Science. 1990;250:942–7.
Usatov AV, Klimenko AI, Azarin KV, Gorbachenko OF, Markin NV, Tikhobaeva VE, Kolosov YA, Usatova OA, Bakoev S, Makarenko M, Getmantseva L. The relationship between heterosis and genetic distances based on SSR markers in Helianthus annuus. Am J Agric Biol Sci. 2014;9:270–6.
Dimitrijevic A, Horn R. Sunflower hybrid breeding: from markers to genomic selection. Front Plant Sci. 2018;8:2238.
Hanson MR, Bentolila S. Interactions of mitochondrial and nuclear genes that affect male gametophyte development. Plant Cell. 2004;16:154–69.
Zhao HF, Huang W, Ahmed SS, Gong ZH, Zhao LM. The pollen wall and tapetum are altered in the cytoplasmic male sterile line RC7 of Chinese cabbage (Brassica campestris ssp pekinensis). Genet Mol Res. 2012;11:4145–56.
Chen Z, Zhao N, Li S, Grover CE, Nie H, Wendel JF, Hua J. Plant mitochondrial genome evolution and cytoplasmic male sterility. Crit Rev Plant Sci. 2017;36:55–69.
Eckardt NA. Cytoplasmic male sterility and fertility restoration. Plant Cell. 2006;18:515–7.
Chen L, Liu YG. Male sterility and fertility restoration in crops. Annu Rev Plant Biol. 2014;65:579–606.
Horn R, Gupta KJ, Colombo N. Mitochondrion role in molecular basis of cytoplasmic male sterility. Mitochondrion. 2014;19:198–205.
Ivanov MK, Dymshits GM. Cytoplasmic male sterility and restoration of pollen fertility in higher plants. Russ J Genet. 2007;43:354–68.
Garayalde AF, Presotto A, Carrera A, Poverene M, Cantamutto M. Characterization of a new male sterility source identified in an invasive biotype of Helianthus annuus (L.). Euphytica. 2015;206:579–95.
Touzet P, Meyer EH. Cytoplasmic male sterility and mitochondrial metabolism in plants. Mitochondrion. 2014;19:166–71.
Siculella L, Palmer JD. Physical and gene organization of mitochondrial DNA in fertile and male sterile sunflower. CMS-associated alterations in structure and transcription of the atpA gene. Nucleic Acids Res. 1988;16:3787–99.
Köhler RH, Horn R, Lössl A, Zetsche K. Cytoplasmic male sterility in sunflower is correlated with the co-transcription of a new open reading frame with the atpA gene. Mol Gen Genet. 1991;227:369–76.
Monéger F, Smart CJ, Leaver CJ. Nuclear restoration of cytoplasmic male sterility in sunflower is associated with the tissue-specific regulation of a novel mitochondrial gene. EMBO J. 1994;13:8–17.
Horn R, Köhler RH. Zetsche, K. A mitochondrial 16 kDa protein is associated with cytoplasmic male sterility in sunflower. Plant Mol Biol. 1991;17:29–36.
Makarenko M, Kornienko I, Azarin K, Usatov A, Logacheva M, Markin N, Gavrilova V. Mitochondrial genomes organization in alloplasmic lines of sunflower (Helianthus annuus) with various types of cytoplasmic male sterility. PeerJ. 2018;6:e5266. https://doi.org/10.7717/peerj.5266.
Reddemann A, Horn R. Recombination events involving the atp9 gene are associated with male sterility of CMS PET2 in sunflower. Int J Mol Sci. 2018;19:806.
Whelan EDP, Dedio W. Registration of sunflower germplasm composite crosses CMG-1, CMG-2, and CMG-3. Crop Sci. 1980;20:832.
Horn R. Molecular diversity of male sterility inducing and male-fertile cytoplasms in the genus Helianthus. Theor Appl Genet. 2002;104:562–70.
Hahn V, Friedt W. Molecular analysis of the cms-inducing MAX1 cytoplasm in sunflower. Theor Appl Genet. 1994;89:379–85.
Ardila F, Echeverria MM, Rios R, Rodriguez RH. Structural features of a cytoplasmic male sterility source from Helianthus resinosus, CMS RES1. J Plant Breeding Crop Sci. 2010;2:168–72.
de la Canal L, Crouzillat D, Quetier F. Ledoigt G. A transcriptional alteration on the atp9 gene is associated with a sunflower male-sterile cytoplasm. Theor Appl Genet. 2001;102:1185–9.
Schnabel U, Engelmann U, Horn R. Development of markers for the use of the PEF1 cytoplasm in sunflower hybrid breeding. Plant Breed. 2008;127:587–91.
Whelan EDP. Hybridization between annual and perennial diploid species of Helianthus. Can J Gen Cytol. 1978;20:523–30.
Makarenko MS, Usatov AV, Markin NV, Azarin KV, Gorbachenko OF, Usatov NA. Comparative genomics of domesticated and wild sunflower: complete chloroplast and mitochondrial genomes. OnLine J Biol Sci. 2015;16:71–5.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for illumine sequence data. Bioinformatics. 2014;30:2114–20.
Nurk S, Bankevich A, Antipov D, Gurevich A, Korobeynikov A, Lapidus A, Prjibelsky A, Pyshkin A, Sirotkin A, Sirotkin Y, Stepanauskas R, McLean J, Lasken R, Clingenpeel SR, Woyke T, Tesler G, Alekseyev MA, Pevzner PA. Assembling genomes and mini-metagenomes from highly chimeric reads. Lect Notes Comput Sci. 2013;7821:158–70.
Langmead B, Salzberg S. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9.
Delcher AL, Bratke KA, Powers EC, Salzberg SL. Identifying bacterial genes and endosymbiont DNA with glimmer. Bioinformatics. 2007;23:673–9.
Lohse M, Drechsel O, Bock R. OrganellarGenomeDRAW (OGDRAW): a tool for the easy generation of high-quality custom graphical maps of plastid and mitochondrial genomes. Curr Genetics. 2007;52:267–74.
Kel AE, Gössling E, Reuter I, Cheremushkin E, Kel-Margoulis OV, Wingender E. MATCH: A tool for searching transcription factor binding sites in DNA sequences. Nucleic Acids Res. 2003;31:3576–9.
Triska M, Grocutt D, Southern J, Murphy DJ, Tatarinova T. cisExpress: motif detection in DNA sequences. Bioinformatics. 2013;29:2203–5.
Triska M, Ivliev A, Nikolsky Y, Tatarinova TV. Analysis of cis-regulatory elements in gene co-expression networks in cancer. Methods Mol Biol. 2017;1613:291–310.
Tanaka Y, Tsuda M, Yasumoto K, Yamagishi H, Terachi T. A complete mitochondrial genome sequence of Ogura-type male-sterile cytoplasm and its comparative analysis with that of normal cytoplasm in radish (Raphanus sativus L.). BMC Genomics. 2012;13:352.
Horn R, Friedt W. CMS sources in sunflower: different origin but same mechanism? Theor Appl Genet. 1999;98:195–201.
Horn R, Kusterer B, Lazarescu E, Prüfe M, Özdemir N, Friedt W. Molecular diversity of CMS sources and fertility restoration in the genus Helianthus. Helia. 2002;25:29–40.
Michael CS, Witkowski J, McCombie WR. Current challenges in de novo plant genome sequencing and assembly. Genome Biol. 2012;13:243.
Liu H, Cui P, Zhan K, Lin Q, Zhuo G, Guo X, Ding F, Yang W, Liu D, Hu S, Yu J, Zhang A. Comparative analysis of mitochondrial genomes between a wheat K-type cytoplasmic male sterility (CMS) line and its maintainer line. BMC Genomics. 2011;12:163.
Sloan DB, Alverson AJ, Chuckalovcak JP, Wu M, McCauley DE, Palmer JD, Taylor DR. Rapid evolution of enormous, multichromosomal genomes in flowering plant mitochondria with exceptionally high mutation rates. PLoS Biol. 2012;10:e1001241. https://doi.org/10.1371/journal.pbio.1001241.
Horn R, Friedt W. Fertility restoration of new CMS sources in sunflower (Helianthus annuus L.). Plant Breed. 1997;116:317–22.
Yang JH, Huai Y, Zhang MF. Mitochondrial atpA gene is altered in a new orf220-type cytoplasmic male-sterile line of stem mustard (Brassica juncea). Mol Biol Rep. 2009;36:273–80.
Tuteja R, Saxena RK, Davila J, Shah T, Chen W, Xiao YL, Fan G, Saxena KB, Alverson AJ, Spillane C, Town C, Varshney RK. Cytoplasmic males-associated chimeric open reading frames identified by mitochondrial genome sequencing of four Cajanus genotypes. DNA Res. 2013;20:485–95.
Yamamoto MP, Kubo T, Mikami T. The 5′-leader sequence of sugar beet mitochondrial atp6 encodes a novel polypeptide that is characteristic of Owen cytoplasmic male sterility. Mol Gen Genomics. 2005;273:342–9.
Kim DH, Kim BD. The organization of mitochondrial atp6 gene region in male fertile and CMS lines of pepper (Capsicum annuum L.). Curr Genet. 2006;49:59–67.
Allen JO, Fauron CM, Minx P, Roark L, Oddiraju S, Lin GN, Meyer L, Sun H, Kim K, Wang C, Du F, Xu D, Gibson M, Cifrese J, Clifton SW, Newton KJ. Comparisons among two fertile and three male-sterile mitochondrial genomes of maize. Genetics. 2007;177:1173–92.
Lee YP, Kim S, Lim H, Ahn YS, Sung SK. Identification of mitochondrial genome rearrangements unique to novel cytoplasmic male sterility in radish (Raphanus sativus L.). Theor Appl Genet. 2009;118:719–28.
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:670–86.
We thank reviewers and the editor for their insightful suggestions and comments on the paper, which has helped us to improve the manuscript. We are grateful to the cooperators of Northern Crop Science Laboratory (Fargo, USA) for providing seeds of HA89 alloplasmic lines of sunflower to the genetic collection of the N. I. Vavilov Institute of Plant Genetic Resources.
The study was supported by the Ministry of Education and Science of Russia project no. 6.929.2017/4.6. Analytical work was carried out on the equipment of centers for collective use of Southern Federal University “High Technology.” The publication costs are funded by the Ministry of Education and Science of Russia project no. 6.929.2017/4.6.
Availability of data and materials
The HA89 fertile line genome is available at https://doi.org/10.6084/m9.figshare.7265648.v1; this sequence will later be deposited to NCBI GenBank. The complete mitochondrion sequence of CMS line HA89(MAX1) has been deposited to GenBank under the accession number MH704580.1.
About this supplement
This article has been published as part of BMC Plant Biology Volume 19 Supplement 1, 2018: Selected articles from BGRS\SB-2018: plant biology. The full contents of the supplement are available online at https://0-bmcplantbiol-biomedcentral-com.brum.beds.ac.uk/articles/supplements/volume-19-supplement-1.
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.
The primers sets used for HA89 (MAX1) genome reorganizations validation and gene expression analysis. *- has been also used in previous study (Makarenko et al., 2018). (DOCX 13 kb)