Next Article in Journal
mTOR: A Potential New Target in Nonalcoholic Fatty Liver Disease
Previous Article in Journal
Non-Fasting Hypertriglyceridemia Burden as a Residual Risk of the Progression of Carotid Artery Stenosis
Previous Article in Special Issue
A Rapid and Efficient Method for Isolation and Transformation of Cotton Callus Protoplast
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 16
10.3390/ijms23169198
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Mitochondrial Genomes between a Cytoplasmic Male-Sterile Line and Its Restorer Line for Identifying Candidate CMS Genes in Gossypium hirsutum

by
Lisha Xuan
1,†,
Guoan Qi
1,2,†,
Xiaoran Li
1,2,
Sunyi Yan
1,
Yiwen Cao
1,2,
Chujun Huang
1,
Lu He
1,
Tianzhen Zhang
1,
Haihong Shang
3,* and
Yan Hu
1,*
1
Institute of Crop Science, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 310058, China
2
Hainan Institute of Zhejiang University, Yazhou Bay Science and Technology City, Sanya 572025, China
3
Zhengzhou Research Base, State Key Laboratory of Cotton Biology, School of Agricultural Sciences, Zhengzhou University, Zhengzhou 450001, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(16), 9198; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23169198
Submission received: 4 July 2022 / Revised: 12 August 2022 / Accepted: 12 August 2022 / Published: 16 August 2022
(This article belongs to the Special Issue Cotton Molecular Genomics and Genetics)
Browse Figures
Review Reports Versions Notes

Abstract

:
As the core of heterosis utilization, cytoplasmic male sterility (CMS) has been widely used in hybrid seed production. Previous studies have shown that CMS is always closely related to the altered programming of mitochondrial genes. To explore candidate CMS genes in cotton (Gossypium hirsutum), sequencing and de novo assembly were performed on the mitochondrial genome of the G. hirsutum CMS line SI3A, with G. harknessii CMS-D2 cytoplasm, and the corresponding G. hirsutum restorer line 0-613-2R. Remarkable variations in genome structure and gene transcripts were detected. The mitochondrial genome of SI3A has three circle molecules, including one main circle and two sub-circles, while 0-613-2R only has one. RNA-seq and RT-qPCR analysis proved that orf606a and orf109a, which have a chimeric structure and transmembrane domain, were highly expressed in abortive anthers of SI3A. In addition, comparative analysis of RNA-seq and full-length transcripts revealed the complex I gene nad4 to be expressed at a lower level in SI3A than in its restorer and that it featured an intron retention splicing pattern. These two novel chimeric ORFs and nad4 are potential candidates that confer CMS character in SI3A. This study provides new insight into the molecular basis of the nuclear–cytoplasmic interaction mechanism, and that putative CMS genes might be important sources for future precise design cross-breeding of cotton.

1. Introduction

Cytoplasmic male sterility (CMS), a maternally inherited trait that inhibits efficient pollen production through a nuclear–cytoplasmic interaction has been widely exploited in the hybrid breeding of eudicot plants such as vegetable brassica, pepper, eggplant, tomato, fennel, carrot, and cotton [1,2,3,4,5,6,7,8]. In particular, CMS lines are usually employed as the female parents in hybridization, thus avoiding self-pollination, which produces non-hybrid seeds. Thus, CMS is of great agronomic importance. Mitochondria are semi-autonomous organelles with their own independent genome, the majority of genes in which are related to the complex subunit of the electron transport chain [9]. Previous studies in higher plants have shown that the CMS phenomenon is usually caused by aberrant recombination of the mitochondrial genome that leads to a highly-expressed chimeric open reading frame (ORF). For example, orfH79 in CMS-HL rice is chimeric with cox1 and an unknown sequence [10], WA352 in CMS-WA rice is composed of three chimeric structures [11], and the chimeric gene, atp6c, in CMS-C maize contains a 481-bp leader sequence that is completely different from wild-type atp6 [12]. It is generally believed that CMS chimeric genes are co-transcribed with genes involved in mitochondrial function, and furthermore, compete with or even inhibit the expression of genuine genes. For instance, orf79 in CMS-BT is homologous with orfH79 and co-transcribed with atp6, and orf108 in Brassica juncea is co-transcribed with atpA, resulting in the CMS phenomenon [13,14,15]. Most CMS chimeric ORFs contain one or more transmembrane domains, such as orf288 in Brassica juncea CMS-hau, atp1-orf522 in CMS-PET1 sunflower, and orf79 and orfH79 in CMS-BT and CMS-HL rice [10,14,15,16,17]. In some cases, the mitochondrial retrograde regulation of nuclear genes causes pollen abortion. For example, in CMS-S maize, the nuclear-encoded gene, ZmDREB1.7, which is induced by mitochondrial retrograde signaling, promotes orf355 accumulation that leads to subsequent male sterility [18]. Currently, the main strategy for identifying CMS genes is to compare CMS lines and their restorer lines at the mitochondrial genome and transcriptome levels.
Cotton is an important cash crop worldwide that provides natural fiber and vegetable oil. As early as 1894, heterosis was reported in the progeny of an interspecific hybrid between Gossypium hirsutum and G. barbadense [19]. Subsequently, more studies demonstrated heterosis performance in cotton to relate to improvements in fiber yield, fiber quality, disease resistance, and stress resistance [20,21,22]. Along with studies examining heterosis, research on cytoplasmic male sterility has also been carried out in cotton. The study of CMS in cotton began in the United States, mainly through distant hybridization and nuclear replacement with the goal of obtaining cytoplasmic male sterility. Since 1965, many cytoplasmic male-sterile lines have been cultivated across the world from the germplasm of G. harknessii, G. arboreum, G. anomalum, G. tuilobum, and G. barbadense. CMS lines of G. hirsutum, including 104-7A and Zhongmiansuo 12-A, have also been obtained by backcross hybridization in China [23]. To uncover the molecular mechanism of CMS in cotton, Wang et al. analyzed mitochondrial DNA from multiple male-sterile lines using the RAPD technology and determined the main driver of CMS in cotton to be the abnormality of the mitochondrial DNA [24]. Using the RFLP method, Feng et al. compared the mitochondrial genomes of CMS G. harknessii and normal fertile G. hirsutum and identified significant differences [25]. Wang et al. further analyzed the mitochondrial proteins and DNA from the anthers and etiolated seedlings of a CMS line and the corresponding restorer line of G. harknessii by SDS-PAGE, RAPD, and RFLP. They identified a 31 kDa polypeptide present in the anther mitochondria of the CMS line at the abortion stage. The RFLP analysis also showed that the mitochondrial DNA of the CMS line lacked a 1.9 kb fragment that has homology with the coxII gene that was present in the restorer line [26]. Li et al. compared the sterile line, 2074A and 2074S, as well as their restorer and maintainer lines and revealed four chimeric ORFs that were differentially transcribed in 2074A [27]. Khan et al. analyzed mitochondrial transcripts in the cotton CMS line, H276A, and identified a protein-coding gene, cox3, for which expression was lower than in the corresponding maintainer line H276B [28]. You et al. integrated methylome and transcriptome analyses for G. barbadense and revealed five key genes that may be associated with CMS [29]. However, although these experimental results have proven that CMS lines and restorer lines differ at the level of mitochondrial DNA, so far, no specific ORF with a chimeric structure similar to those implicated in the sterility of other plants has been identified in the mitochondrial genomes of CMS cotton lines.
In this study, we first reported the assembly of two complete mitochondrial genome sequences for the cotton CMS line, SI3A, and its restorer line, 0-613-2R, and further sequenced their full-length transcriptomes. The restoring line 0-613-2R completely restored fertility of the CMS line in the G. harknessii cytoplasm. A comparative analysis of the mitochondrial genomes and transcriptomes revealed differences in the sequences and expression patterns between the CMS and restorer lines and yielded several CMS-associated predicted chimeric ORFs that may be responsible for male sterility in SI3A. All told, the assembled cotton mitochondrial genome, the full-length transcriptome, and the set of candidate CMS genes in SI3A will contribute to our understanding of the molecular mechanism of CMS in cotton.

2. Results

2.1. Flower Morphology of the CMS-D2 Line SI3A and Its Restorer Line 0-613-2R

Plants of the G. harknessii CMS-D2 cytoplasm line, SI3A, grew normally, and no morphological differences in vegetative growth were observed between SI3A and its restorer line, 0-613-2R. During the flowering period, the anthers of 0-613-2R were plump and the ovary enlarged, while in SI3A the stigma was elongated, the filaments were significantly shorter, and the anther sacs were completely wrinkled and did not dehisce (Figure 1A–D). The pollen grains of 0-613-2R were stained darkly by the I2-KI solution, while no pollen grains were stained for SI3A (Figure 1E,F). These findings prove that the CMS line, SI3A, is of the complete stamen degeneration type.

2.2. Mitochondrial Genome Assembly and Annotation

The mitochondrial genomes of the CMS-D2 line, SI3A, and its restorer line, 0-613-2R, were sequenced with average depths of 167× and 312× using Illumina Hiseq and Nanopore sequencing technology. In the Illumina sequencing, 44.34- and 40.35-million clean reads were obtained from 0-613-2R and SI3A, respectively, and the corresponding Q20 values were 97.42% and 97.39%. In the Nanopore sequencing, 2.6- and 2.8-million clean reads were obtained from 0-613-2R and SI3A with total lengths of 14,071 Mb and 14,224 Mb, respectively; the lengths obtained for scaffold N50 were 8199 bp and 7142 bp. To verify the accuracy of the sequencing results, both Illumina and Nanopore sequencing reads were mapped to the cotton G. hirsutum genome (GCA_007990345.1) and the mitochondrial genome (NC_027406.1). Supplementary Table S1 gives the average percentage of cytoplasm DNA (ctDNA) and nuclear DNA (ncDNA) identified in the Illumina and Nanopore sequence reads. Collinearity analysis revealed the SI3A and 0-613-2R mitochondrial genomes to be highly conserved with G. hirsutum (NC_027406.1) and G. harknessii (NC_027407.1); the 95.66% and 94.96% sequence identity (Figure 2) further indicates the conservation of the two mitochondrial genomes.
The features of the mitochondrial genomes of the SI3A CMS line and its restorer line 0-613-2R are listed in Table 1. The SI3A mtDNA exhibited a more complex composition than that of the restorer line. Namely, the SI3A mtDNA consisted of one main genomic circle (SI3Amt2, 396,206 bp) and two sub-circle molecules (SI3Amt1, 63,927 bp; SI3Amt3, 174,070 bp). In contrast, the 0-613-2R mitochondrial genome assembled into a single circular molecule with a size of 607,367 bp. The guanine and cytosine (GC) content of the two mitochondrial genomes was approximately 45%, which is close to the median values of other fully-sequenced mitochondrial genomes from seed plants, such as 45.08% in Brassica juncea, 44.3% in wheat, 44.8% in Arabidopsis thaliana, and 45% in soybean [30,31,32,33].
Mitochondrial genes were annotated using the Mitofy software and compared against the protein-coding sequences from closely-related species to locate encoding regions. Within the SI3A mtDNA, a total of 33 protein-coding genes, 6 rRNAs, 25 tRNAs, and 59 unidentified ORFs were annotated, while the 0-613-2R mtDNA yielded annotations for 33 protein-coding genes, 6 rRNAs, 27 tRNAs, and 48 unidentified ORFs (Figure 3A,B). Notably, the two mitochondrial genomes were identical in terms of protein-coding gene number, while SI3A differed from 0-613-2R in having the atp1 gene split into SI3Amt1 (190 bp) and SI3Amt3 (1335 bp), as well as lacking 87 bp from the end of atp9. Overall, protein-coding regions accounted for 25.8%, 19.0%, 7.4%, and 8.4% of the 0-613-2Rmt, SI3Amt1, SI3Amt2, and SI3Amt3 molecules, respectively. Detailed information describing the tRNA and rRNA gene content of the CMS and restorer line mitotypes is given in Supplementary Table S3. Interestingly, although the mitochondrial genome of SI3A was 26.8 Kb larger overall than that of the restorer line, 0-613-2R, it also featured fewer large repeats (>1 Kb). The presumed reason for this is that the repeated sequence that mediated gene rearrangements in the sterile line SI3A during assembly of the mitochondrial genome created its three circular DNA molecules, and thereby reduced the number of repeat sequences.

2.3. Comparative Analysis of the SI3A and 0-613-2R Mitochondrial Genomes

2.3.1. Repeat Sequences

CMS in plants is maternally inherited and is generally controlled by mitochondrial genes. In most cases of CMS, rearrangements in the mitochondrial genome that are mediated by repeat sequences produce a new chimeric ORF that encodes proteins containing transmembrane domains [34,35,36]. Repeated sequences of greater than 50 bp annotated in the mitochondrial genomes of the SI3A CMS line and its restorer line are presented in Table 2 and Supplementary Table S4. Almost threefold as many repeats were detected in the restorer line than in the CMS line, at 77 repeats (77,020 bp, 12.8%) and 20 repeats (26,460 bp, 4.3%), respectively. All 20 repeats in SI3A occurred in SI3Amt2, the largest circular component of its mitochondrial genome; some shorter repeats (<50 bp) were present in the sub-genomes, SI3Amt1 and SI3Amt3. Five inverted repeats (2000 bp < IR > 1000 bp) were detected in 0-613-2R, all of which were absent from SI3Amt2 (Table 2). Two large repeats (>10 kb) were also identified in 0-613-2R, comprising one invert repeat (IR) and one direct repeat (DR) (Supplementary Table S4); only one large DR (>10 kb) was found in SI3Amt2, which showed 99% identity with those in the restorer line. The three large repeats also contained a few small repeats (<100 bp). Some sequences had multiple repeats with different loci in the two mitochondrial genomes.
Notably, although the CMS line harbored fewer repeat sequences than the restorer line, and no long repeats (>2000 bp), direct repeats accounted for 50% of its repetitive sequences, twofold the proportion in the restorer line. This phenomenon may be related to gene recombination. Homologous recombination across large direct repeats separates the genome into sub-genomic molecules of different sizes, while recombination of invert repeats causes a sequence insertion [37]. The significant reduction of repeat sequences in the CMS line may be due to intragenomic homologous recombination of the largest direct repeat (10,637 bp), especially at non-coding sequences within the gene region, thereby causing the splitting of the SI3A mitochondrial genome into one main circle and two sub-circle structures.

2.3.2. Structural Variation, SNPs, and InDels

Genome structural variations, SNPs, and InDels between CMS line SI3A and its restorer line, 0-613-2R, were identified in order to explore structural and sequence variations in known mitochondrial genes and the predicted ORFs. In terms of structural variation, a total of 15 translocations, 1 inversion, and 11 translocation + inversion regions were identified (Figure 4A). Two insertions longer than 100 bp were also identified in syntenic regions, and 5 complex-InDel regions were detected.
With regard to sequence variation, a total of 714 SNPs and 138 InDels (73 insertions and 65 deletions) were detected between the 2 mitochondrial genomes. The vast majority of identified SNPs occurred in non-coding regions, of which 700 (98% of the total SNPs) were distributed in intergenic regions, and only 14 (2% of the total SNPs) were identified in known mitochondrial genes, including cox1, cox3, nad7, atp4, atp8, sdh3, matR, rps4, rpl2, rpl5, rpl10, and rpl16. Five were synonymous mutations, and not expected to cause the CMS, while nine were non-synonymous mutations that altered the encoded amino acid residues in nine mitochondrial genes (cox3, atp8, sdh3, matR, rps4, rpl2, rpl5, and rpl10) (Table 3). These non-synonymous changes are candidates for functional scrutiny in investigating the molecular mechanisms of CMS, since protein-coding genes in plant mtDNA are extraordinarily conserved and their rate of evolution is very low [38,39,40]. These SNPs primarily concerned A and C bases, and most were actually transversions rather than transitions.
The identified InDels had lengths in the range of 1–10 bp, with almost all being 1 bp (Figure 4B). Most were sited in intergenic regions; one insertion of a “T” in the coding region of SI3A atp1 resulted in a frameshift mutation affecting 44 amino acids (Figure 4C). Mutations in atp1 may impair ATP synthesis in the mitochondrial respiratory chain. No other InDels resulted in a frameshift.

2.4. Unique ORFs

To date, CMS-associated genes have been mainly identified by whole-genome sequencing of mitochondria and subsequent comparison with a standard cultivar to identify unique ORFs. In this study, 30 and 38 unique ORFs encoding >100 amino acids were identified in the 0-613-2R restorer line and the SI3A CMS line. respectively (Supplementary Tables S5 and S6). Moreover, 22 specific ORFs in SI3A were determined to be distinct from ORFs in 0-613-2R (Figure 5A). CMS-associated genes are known to be characterized by chimeric structures or to be co-transcribed with functional genes, and said functional genes usually contain transmembrane domains. Here, only two ORFs were predicted to specifically exist in the CMS line, have a chimeric structure, and be co-transcribed with genes containing transmembrane domains. In addition, at least one transmembrane domain was found in these two ORFs. The first ORF, orf119a-2, shared 73 bp of sequence with the known gene nad7 and also featured 284 bp of unknown sequence, while the second, orf606a, shared 1336 bp at its 5′ end with the known gene atp1 and contained another 485 bp of unique sequence (Figure 5B). We also chose to examine two other unique ORFs identified in SI3A, orf131a and orf109a, which featured chimeric structures. Of those, orf131a contained a 92 bp sequence corresponding to the 3′ end of nad6 plus 304 bp of unknown sequence, while orf109a fused 35 bp corresponding to the 3′ end of cox2 along with a 291 bp sequence of unknown origin (Figure 5B). Since having a transmembrane domain is another characteristic of CMS genes, we performed a structural analysis of the four selected ORFs; this revealed that orf119a-2 and orf606a each contain a transmembrane structural domain, while orf131a and orf109a do not (Figure 5C–F).

2.5. Selection of Candidate Genes Resulting in Male Sterility of SI3A

To further explore the expression of CMS-associated genes and their intron splicing patterns in SI3A, Illumina and full-length PacBio transcriptome sequencing were performed. Specifically, RNA from abortion stage tissue (4–5 mm floral bud) of CMS line SI3A and its restorer line, 0-613-2R, was sequenced using the Illumina NovaSeq 6000 platform and PacBio Sequel II System. PacBio sequencing yielded a total of 42.2 M and 28.0 M subreads from 0-613-2R and SI3A, and the length of scaffold N50 was 2111 bp and 2122 bp, respectively. For Illumina RNA-Seq, 150-bp paired-reads were generated and the Q20 percentages exceeded 97% in each sample. In total, 82.2 Mb and 77 Mb of raw reads were obtained for 0-613-2R and SI3A, respectively, corresponding to 78.8 Mb and 73.5 Mb of clean reads after filtering. Clean reads were mapped to the two mitochondrial genomes using DEG-seq2 for a transcript-level examination of mitochondrial protein-coding genes (Table 4). Novel transcript prediction subsequently identified 99 and 91 novel transcripts in SI3A and 0-613-2R, respectively, including 34 and 33 protein-coding transcripts. Of protein-coding genes in SI3A, we found rps10 to be the most highly expressed and nad4 the least expressed (Figure 6A). In addition to quantifying expression levels from the RNA-seq data and visual inspection in IGV, RT-PCR analysis was performed for the four ORFs in abortion floral buds from SI3A and 0-613-2R. The ORFs did not show obvious differences between the two lines except for orf606a, which was significantly more highly expressed in SI3A (Figure 6C). We also analyzed relative gene expression using RT-qPCR and found that orf109a and especially orf606a are highly expressed in SI3A (Figure 6B), which is consistent with the RT-PCR results. Sequence comparison revealed orf606a to be partially identical to orf610a, which was previously reported as the CMS candidate gene in cotton; in particular, a transgenic Arabidopsis thaliana line harboring orf610a with a mitochondrial signal peptide exhibits partial pollen abortion, shortened siliques, shortened filaments, and pollen tubes [41]. Accordingly, orf606a is a candidate sterility gene for validation in future work using transgenic cotton.
The present work identified a total of seven intron-containing protein-coding genes in the two mitochondrial genomes. Full-length PacBio sequencing revealed SI3A to exhibit an intron retention splicing pattern for nad4, in which intron 1 was almost completely preserved and intron 3 was partially preserved (Figure 6D); corresponding transcripts were detected in RNA-seq data from SI3A but not in data from 0-613-2R. Therefore, the differential transcript caused by intron retention may be related to the low expression of nad4 in the CMS.

3. Discussion

3.1. Characteristics of the SI3A and 0-613-2R Mitochondrial Genomes

The complete mitochondrial genomes of CMS line SI3A and the corresponding restorer line, 0-613-2R, were constructed through sequencing and de novo assembly. The mtDNA sequences of SI3A and 0-613-2R were nearly identical to those of the 2074A and 2074S G. hirsutum CMS lines for which mitochondrial genomes were previously reported [27]. The differences in mitochondrial genome size between CMS and restorer lines may be caused by unique sequences, as has been observed in tomato and wheat [32,42]. Here, unique sequences in SI3Amt2 (21.3 kb) and SI3Amt3 (1.1 kb) were the main factors leading to the size difference between SI3A and 0-613-2R. Notably, high-frequency homologous recombination mediated by large repeat sequences results in sub-circle structures, such as those observed in SI3A. Higher plant mitochondrial genomes are frequently composed of a main circle structure and multiple sub-genomic molecules [43]; for example, the mitochondrial genome of monkeyflower (Mimulus guttatus) contains three large repeats for which four possible recombination relationships have been deduced, ultimately producing eight sub-ring structures [44]. In soybean, the CMS-K line, NJCMS1B, features a main circle, a linear sequence, and a sub-genomic circle [31]. Here, one and two large (>10 kb) repeat sequences were identified in SI3A and 0-613-2R, respectively (Table 2 and Supplementary Table S4). Intra-genomic homologous recombination of the large repeat sequence in SI3A resulted in its main circle and two sub-circle structures (Figure 3B). Notably, SI3A exhibited a larger genome but significantly fewer repeats than 0-613-2R, which may be due to the increased frequency of gene rearrangements caused by the generation of the sub-circle structures, ultimately disrupting the original repeat sequences.
Functional mitochondrial genes are, in general, quite conserved. Between-species differences in the number of mitochondrial coding genes are mainly caused by multi-copy genes, ribosomal RNA, and tRNA [45]. Here, both examined genomes harbored 33 protein-coding genes (Table 1), the finding of which is consistent with 2074A and 2074S—except for the loss of nad1 and nad2 [27]. This discrepancy may be related to the sequencing depth and annotation method. Additionally, a comparative analysis revealed nine non-synonymous SNP mutations between SI3A and 0-613-2R, which occurred in eight protein-coding genes (cox3, atp8, sdh3, matR, rps4, rpl2, rpl5, and rpl10), along with just one exonic InDel site, which caused a frameshift mutation in atp1 (Table 3 and Figure 4A,B). Mutations in atp1 may disrupt the assembly of complex V, and hence ATP synthesis in the mitochondrial electron transport chain.

3.2. orf606a and nad4 May Be CMS Candidate Genes for the CMS-D2 Line SI3A

CMS genes are always associated with novel chimeric ORFs created by mitogenome recombination and are often co-transcribed with upstream or downstream protein-coding genes. For example, in Raphanus sativus DCGMS, orf463 represents a partial fusion of cox1 with an unidentified sequence [46], while in CMS-S maize, orf77 contains three segments derived from atp9 [47]. This study identified four ORFs specific to SI3A, orf109a, orf119a-2, orf606a and orf131a (Figure 5B), which were all chimeric with protein-coding genes and contained one or more transmembrane domains (Figure 5C–F). Of those, orf606a, which was chimeric with atp1 and an unknown sequence, was found to be especially highly expressed in abortion floral buds of SI3A compared with 0-613-2R (Figure 6B,C). In addition, the sequence of orf606a is homologous with the CMS gene, orf610a, previously reported in the CMS-D2 line ZBA, which causes ROS accumulation and ATP reduction [41], and its expression pattern is consistent with other known CMS genes, such as orf256 in wheat, cox3 in H276A cotton, and orf326 in Nau radish [28,48,49].
Most reported restorer-of-fertility (Rf) genes that counteract male sterility encode members of the pentatricopeptide repeat (PPR) family; these include petunia Rf [50], Brassica napus Rfk1 [51], Rahanus stivus Rfo [52], rice Rf1 [53], sorghum Rf1 [54], and others. Most PPR proteins localize to the mitochondria and chloroplasts and play important roles in post-transcriptional regulation within those organelles, mainly in RNA splicing, RNA cleavage, RNA editing, and RNA stabilization [55,56]. For example, the PPR protein, Rf6, is known to function with OsHXK6 to promote processing of the atp6–orfH79 transcript and restore fertility in rice [57], while OTP43 is essential for trans-splicing of nad1 intron 1 in Arabidopsis thaliana [58]. In this study, RNA-seq and full-length transcript data revealed nad4 to have low expression in SI3A, and furthermore that retention of intron 1 resulted in the production of abnormal transcripts (Figure 6A,D). In the current literature, most known PPR proteins are involved in cis–trans-splicing of nad1, nad2, and nad4. For example, the PPR protein, MTSF3, is involved in 3’ end processing and intron trans-splicing of mitochondrial nad2 mRNA [59]. Similarly, in CMS-S maize, PPRK2 (Rf3) targets mitochondria, binds to orf355, and inhibits its editing and degradation to restore fertility [60]. The PPR protein, Ghlm, has likewise been demonstrated to be necessary for splicing of nad7 mRNA during cotton fiber development [61]. In light of these reports, we speculate that the abnormal transcript generated from nad4 might be the target of the cotton restorer protein, Rf1. Thus, orf606a and nad4 may be potential candidate CMS-associated genes in cotton. We hope these results will contribute to increasing our understanding of the cause of CMS in cotton and provide candidate CMS genes for hybrid breeding.

4. Materials and Methods

4.1. Plant Materials and Phenotypic Analysis

The plant materials used in this study were the Gossypium hirsutum nuclear restorer line, 0-613-2R, with viable pollen, and SI3A, an upland CMS-D2 line with cytoplasm from G. harknessii. Plants were grown under standard growing conditions in the experimental field of Zhejiang University (30.1° N, 120.1° E), Hangzhou, in which the average daily temperature ranged from 30 °C to 35 °C. Floral buds of about 4–5 mm in length (the abortion stage) from 100 plants, each of CMS-D2 line SI3A and its restorer line 0-613-2R, were collected and combined for RNA extraction. All harvested samples were frozen in liquid nitrogen and stored at −80 °C before use. Flower morphology was observed by stereoscopic microscopy (Stemi 508, Suzhou, China). Pollen viability was evaluated on the flowering day by staining with 1% I2-KI solution after crushing anthers, with the staining observed under a microscope (XSP-37XC, Shanghai, China).

4.2. Library Construction, Sequencing, and Assembly of Mitochondrial Genomes

High-quality DNA was extracted from leaves using a modified CTAB method [62]. The purity, concentration, and integrity of the DNA samples were confirmed by Nanodrop (OD260/280 >1.8) and agarose gel electrophoresis. Afterwards, the purified DNA was fragmented with a Covaris instrument and used to construct a paired-end library with an average fragment length of 350 bp using the NexteraXT DNA Library Preparation Kit (Illumina, San Diego, CA, USA). Sequencing was then performed on the Illumina Novaseq 6000 platform (Huitong Biotechnology Co. Ltd., Shenzhen, China). In parallel, fragmentation was performed using a BluePippin system (Sage Science, Beverly, MA, USA), then a single-end library (average length 8 kb), was constructed using the SQK-LSK109 kit (Oxford Nanopore, Oxford, UK) following the manufacturer’s instructions, and was sequenced on the Nanopore PromethION sequencing platform (Huitong biotechnology Co. Ltd., Shenzhen, China).
Before assembly, the original Illumina reads were filtered to remove reads with adaptors, reads showing a quality score below 20 (Q < 20), short fragments of less than 50 bp after adapter removal, and reads containing greater than 10% uncalled bases (“N” characters). The mitochondrial genome was reconstructed using the Oxford Nanopore data and Illumina NovaSeq data in combination as follows. First, the clean Illumina sequencing data were de novo assembled by SPAdes v3.11.0 (Bankevich A, Petersburg, Russia) using default parameters without a cutoff. Second, the scaffolds were filtered using BlastN and Exnerate with an alignment threshold of e-value 1e-10 and protein similarity threshold of 70%. Selected scaffolds having homology to the Gossypium arboreum (NCBI: KR736342) mitochondrial genome and protein-coding gene sequences were removed as fragments that were low coverage and obviously not the target genome. The collected fragmented target sequences were extended and merged using PRICE and MITObim [63] to minimize the number of scaffolds. Then, the Nanopore reads were aligned back to our filtered SI3A and 0-613-2R scaffolds with minimap2 [64], segregated by aligned reads, and reassembled de novo with Flye [65]. The final genome sequences of the CMS-D2 line, SI3A, and its restorer line, 0-613-2R, were obtained by polishing with pilon [66] using Illumina sequencing reads and made publicly available through the National Center for Biotechnology Information (NCBI) under the GenBank accession numbers: SI3Amt1, ON931349; SI3Amt2, ON931350; SI3Amt3, ON931351; and 0-613-2R, ON931352.

4.3. Mitochondrial Genome Annotation and Identification of ORFs

The mitochondrial genes were annotated by Motify [67]. Protein sequences were aligned with the cotton mitochondrial reference genomes from Gossypium arboreum (NCBI: KR736342), G. hirsutum (NC_027406.1), and G. harknessii (NC_027407.1), and the resulting gene sets were then integrated and manually corrected to obtain the final list of protein-coding genes. Transfer RNA (tRNA) and ribosomal RNA (rRNA) genes were predicted using the homology prediction methods of tRNAscan-E and rRNAmmer [68]. Gene function annotation was carried out using BLAST with e-value 1 × 10−10 and the following databases: Kyoto Encylopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), Swiss-Prot, Gene Ontology (GO), and the NCBI non-redundant protein sequences (NR). ORFs were predicted with the NCBI ORF Finder (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/orffinder/) (accessed on 9 April 2022) and retained if the hypothetical protein encoded by the ORF was longer than 100 amino acids. ORFs were named according to the number of amino acids encoded. Hypothetical protein transmembrane domains were predicted by the TMHMM online server v2.0 http://www.cbs.dtu.dk/services/TMHMM/ (accessed on 6 May 2022).

4.4. Differential Analysis of the SI3A and 0-613-2R Mitochondrial Genomes

Using restorer line 0-613-2R as the reference genome, syntenic analysis was performed to compare the sequence information of these two mitochondrial genomes in detail. First, large-scale synteny between the genomes was determined by MUMmer 4.0 (Alekley Zimin, Baltimore, MD, USA) [69]; then, LASTZ https://lastz.github.io/lastz/ (accessed on 31 September 2020) was used to compare blocks and confirm the local position alignment. MUMmer 4.0 and BLAST were employed for global alignment to identify differential single nucleotide polymorphism (SNP) loci. InDel sequences were obtained from the LASTZ alignment and verified by BWA [70,71] and SAMtools [72] filtering to arrive at a reliable set of InDels. Repeat sequences were predicted using RepeatMasker 4.1.1 https://www.repeatmasker.org/ (accessed on 15 October 2020).

4.5. RNA Extraction and Transcriptome Sequencing

Total RNA was extracted from floral buds (4–5 mm size, abortion stage) of SI3A and 0-613-2R using the Spectrum Plant Total RNA Kit (Sigma-Aldrich, St. Louis, MO, USA) with three replicates for each sample. RNA concentration and integrity were assessed with the Fragment Analyzer 5400 (Agilent Technologies, CA, USA). Sequencing libraries were generated using the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA), and index codes were added in order to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads and a fragmentation reaction was carried out under elevated temperature to break the mRNA into short template pieces. Second-strand cDNA synthesis was subsequently performed using DNA polymerase I and RNA H, after which the library fragments were purified using the AMPure XP system (Beckman Coulter, Beverly, MA, USA). After adaptor-ligase and end repair reactions, suitable fragments were selected for PCR amplification, the resulting PCR products were purified (AMPure XP system) (Beckman Coulter, Beverly, MA, USA), and the library quality was assessed on the Agilent Bioanalyzer 2100 system. Clustering was then performed on a cBot Cluster Generation System using the TruSeq PE Cluster Kit v3-cBot-HS (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. After cluster generation, the prepared libraries were sequenced on an Illumina NovaSeq 6000 Platform, generating 150 bp paired-end reads. Raw reads were filtered by discarding any pair in which one read contained adapter contamination, more than 10% uncertain bases, or more than 50% low quality (Phred quality <5) bases. The clean sequence files of these six samples (FASTQ files) are accessible from the NCBI through the Sequence Read Archive (SRA) under accession number SRR20067583-SRR20067588.

4.6. Full-Length Transcriptome Sequencing

The total RNA from 4–5 mm floral buds was first assessed for quality by measuring RNA Integrity Number (RIN) and concentration using an Agilent 2100 instrument; then, 2 µg of total RNA was used as the input material for RNA sample preparation. The SMRTbellTM Template Prep Kit (Pacific Biosciences, Menlo Park, CA, USA) was used to generate SMRTbellTM libraries by following the manufacturer’s instructions for end-repair, A-tailing, and adaptor ligation, after which the libraries were sequenced on the PacBio Sequel II System (Pacific Biosciences, CA, USA). Sequence data were processed using the SMRT Analysis software (ISO-seq version 3.0) (Pacific Biosciences, CA, USA) with parameters: min length, 300; min accuracy, 0.75; and min passes, 0. Finally, high-quality isoforms and polished low-quality isoforms were mixed for subsequent alignment to the reference genome using GMAP [73,74] with parameters -min-trimmed-coverage 0.85 and -min-identity 0.9 against reference genome. Any redundancy in polished consensus reads was removed by TOFU (-min-coverage 0.85 and -min-identity 0.9 against reference genome) to obtain final transcripts for subsequent analysis.

4.7. RT-PCR and RT-qPCR Analysis of Candidate ORFs

Floral buds (4–5 mm, abortion stage) of the CMS line, SI3A, and its restorer line, 0-613-2R, were used in measuring the expression of candidate ORFs. First, the total RNA was extracted with the Spectrum Plant Total RNA Kit (Sigma-Aldrich, St. Louis, Mo, USA); then, reverse transcription was performed using the HiScript QRT SuperMix for RT-qPCR kit (Vazyme, Nanjing, China) following the manufacturer’s protocol. Semiquantitative RT-PCR was then performed with the following temperature program: 95 °C for 5 min, then 32 cycles of 95 °C for 15 s, 58 °C for 15 s, 72 °C for 30 s, and 5 min at 95 °C. RT-qPCR was performed with the ABI StepOnePlus system using AceQ RT-qPCR SYBR Green Master Mix (Vazyme, Nanjing, China). Values were standardized according to the smallest sample threshold cycle number (Ct value) and the highest fluorescence value, and relative gene expression levels were calculated using the 2−ΔΔCT method [75]. GhHis3 was used as the internal reference gene for RT-PCR and RT-qPCR. Supplementary Table S2 lists the primers used in RT-PCR and RT-qPCR.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms23169198/s1.

Author Contributions

L.X. performed the experiments, analysed the data and drafted the manuscript under the supervision of Y.H.; G.Q. made some figures and modified the manuscript; X.L., S.Y., Y.C., C.H. and L.H. participated in collecting floral buds, DNA extraction, etc. T.Z., H.S. and Y.H. revised the manuscript. All authors discussed results and provide constructive advices. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China grant number 321702005 to Yan Hu. Leading Innovative and Entrepreneur Team Introduction Program of Zhejiang, grant number 2019R01002 to Tianzhen Zhang.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The mt genome of 0-613-2R and SI3A are available at NCBI GenBank, under the accession number: SI3Amt1, ON931349; SI3Amt2, ON931350; SI3Amt3, ON931351; 0-613-2R, ON931352. The Illumina RNA-seq data are available at the Sequence read Archive under accession number SRR20067583-SRR20067588.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Singh, S.; Dey, S.S.; Bhatia, R.; Kumar, R.; Behera, T.K. Current understanding of male sterility systems in vegetable Brassicas and their exploitation in hybrid breeding. Plant Reprod. 2019, 32, 231–256. [Google Scholar] [CrossRef] [PubMed]
  2. Singh, S.; Dey, S.S.; Bhatia, R.; Kumar, R.; Sharma, K.; Behera, T.K. Heterosis and combining ability in cytoplasmic male sterile and doubled haploid based Brassica oleracea progenies and prediction of heterosis using microsatellites. PLoS ONE 2019, 14, e0210772. [Google Scholar] [CrossRef]
  3. Wang, P.; Lu, Q.; Ai, Y.; Wang, Y.; Li, T.; Wu, L.; Liu, J.; Cheng, Q.; Sun, L.; Shen, H. Candidate gene selection for cytoplasmic male sterility in pepper (Capsicum annuum L.) through whole mitochondrial genome sequencing. Int. J. Mol. Sci. 2019, 20, 578. [Google Scholar] [CrossRef] [PubMed]
  4. Khan, M.M.R.; Isshiki, S. Cytoplasmic male sterility in eggplant. Hortic. J. 2016, 85, 1–7. [Google Scholar] [CrossRef]
  5. Kuwabara, K.; Arimura, S.-I.; Shirasawa, K.; Ariizumi, T. orf137 triggers cytoplasmic male sterility in tomato. Plant Physiol. 2022, 189, 465–468. [Google Scholar] [CrossRef]
  6. Palumbo, F.; Vitulo, N.; Vannozzi, A.; Magon, G.; Barcaccia, G. The mitochondrial genome assembly of fennel (Foeniculum vulgare) reveals two different atp6 gene sequences in cytoplasmic male sterile. Int. J. Mol. Sci. 2020, 21, 4664. [Google Scholar] [CrossRef]
  7. Linke, B.; Alessandro, M.S.; Galmarini, C.R.; Nothnagel, T. Carrot floral development and reproductive biology. In The Carrot Genome; Simon, P., Lorizzo, M., Grzebelus, D., Baranski, R., Eds.; Springer: New York, NY, USA, 2019; Volume 126, pp. 415–423. [Google Scholar]
  8. Wu, J.Y.; Zhang, M.; Zhang, B.B.; Zhang, X.X.; Guo, L.P.; Qi, T.X.; Wang, H.L.; Zhang, J.F.; Xing, C.Z. Genome-wide comparative transcriptome analysis of CMS-D2 and its maintainer and restorer lines in upland cotton. BMC Genom. 2017, 18, 454. [Google Scholar] [CrossRef]
  9. Singh, M.; Brown, G.G. Suppression of cytoplasmic male sterility by nuclear genes alters expression of a novel mitochondrial gene region. Plant Cell 1991, 3, 1349–1362. [Google Scholar] [CrossRef]
  10. Peng, X.J.; Wang, K.; Hu, C.F.; Zhu, Y.L.; Wang, T.; Yang, J.; Tong, J.P.; Li, S.Q.; Zhu, Y.G. The mitochondrial gene orfH79 plays a critical role in impairing both male gametophyte development and root growth in CMS-Honglian rice. BMC Plant Biol. 2010, 10, 125. [Google Scholar] [CrossRef]
  11. Luo, D.P.; Xu, H.; Liu, Z.L.; Guo, J.X.; Li, H.Y.; Chen, L.T.; Fang, C.; Zhang, Q.Y.; Bai, M.; Yao, N.; et al. A detrimental mitochondrial-nuclear interaction causes cytoplasmic male sterility in rice. Nat. Genet. 2013, 45, 573–577. [Google Scholar] [CrossRef]
  12. Yang, H.L.; Xue, Y.D.; Li, B.; Lin, Y.N.; Li, H.C.; Guo, Z.Y.; Li, W.H.; Fu, Z.Y.; Ding, D.; Tang, J. The chimeric gene atp6c confers cytoplasmic male sterility in maize by impairing the assembly of the mitochondrial ATP synthase complex. Mol. Plant 2022, 15, 872–886. [Google Scholar] [CrossRef]
  13. Ashutosh; Kumar, P.; Dinesh Kumar, V.; Sharma, P.C.; Prakash, S.; Bhat, S.R. A novel orf108 co-transcribed with the atpA gene is associated with cytoplasmic male sterility in Brassica juncea carrying Moricandia arvensis cytoplasm. Plant Cell Physiol. 2008, 49, 284–289. [Google Scholar] [CrossRef]
  14. Kojima, H.; Kazama, T.; Fujii, S.; Toriyama, K. Cytoplasmic male sterility-associated ORF79 is toxic to plant regeneration when expressed with mitochondrial targeting sequence of ATPase γ subunit. Plant Biotechnol. 2010, 27, 111–114. [Google Scholar] [CrossRef]
  15. Kazama, T.; Itabashi, E.; Fujii, S.; Nakamura, T.; Toriyama, K. Mitochondrial ORF79 levels determine pollen abortion in cytoplasmic male sterile rice. Plant J. 2016, 85, 707–716. [Google Scholar] [CrossRef]
  16. Jing, B.; Heng, S.P.; Tong, D.; Wan, Z.J.; Fu, T.D.; Tu, J.X.; Ma, C.Z.; Yi, B.; Wen, J.; Shen, J. A male sterility-associated cytotoxic protein ORF288 in Brassica juncea causes aborted pollen development. J. Exp. Bot. 2012, 63, 1285–1295. [Google Scholar] [CrossRef]
  17. Moneger, F.; Smart, C.J.; Leaver, C.J. 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. [Google Scholar] [CrossRef]
  18. Xiao, S.L.; Zang, J.; Pei, Y.R.; Liu, J.; Liu, J.; Song, W.; Shi, Z.; Su, A.G.; Zhao, J.R.; Chen, H. Activation of mitochondrial orf355 gene expression by a nuclear-encoded DREB transcription factor causes cytoplasmic male sterility in maize. Mol. Plant 2020, 13, 1270–1283. [Google Scholar] [CrossRef]
  19. Mell, P.H. Experiments in Crossing for the Purpose of Improving the Cotton Fiber; Alabama Agricultural Experiment Station (AAES) Report; Agricultural Experiment Station of the Agricultural and Mechanical College: Auburn, AL, USA, 1894; Volume 6. [Google Scholar]
  20. Li, C.; Yu, H.R.; Li, C.; Zhao, T.L.; Dong, Y.T.; Deng, X.L.; Hu, J.H.; Zhang, Y.; Zhang, F.; Li, J. QTL mapping and heterosis analysis for fiber quality traits across multiple genetic populations and environments in upland cotton. Front. Plant Sci. 2018, 9, 1364. [Google Scholar] [CrossRef]
  21. Munir, S.; Hussain, S.B.; Manzoor, H.; Quereshi, M.K.; Zubair, M.; Nouman, W.; Shehzad, A.N.; Rasul, S.; Manzoor, S.A. Heterosis and correlation in interspecific and intraspecific hybrids of cotton. Genet. Mol. Res. 2016, 15, gmr8083. [Google Scholar] [CrossRef]
  22. Zaidi, S.S.A.; Naqvi, R.Z.; Asif, M.; Strickler, S.; Shakir, S.; Shafiq, M.; Khan, A.M.; Amin, I.; Mishra, B.; Mukhtar, M.S.; et al. Molecular insight into cotton leaf curl geminivirus disease resistance in cultivated cotton (Gossypium hirsutum). Plant Biotechnol. J. 2020, 18, 691–706. [Google Scholar] [CrossRef]
  23. Zhang, T.Z.; Jing, S.R. Theory and Practice of Selection and Breeding of Male Sterile Hybrids in Cotton; China Agriculture Press: Beijing, China, 1998. [Google Scholar]
  24. Wang, X.D.; Zhang, T.Z.; Pan, J.J. Cytological observation of microsporogenesis and RAPD analysis of mitochondrial DNAs for cytoplasmic male sterile cotton lines. Sci. Agric. Sin. 1998, 31, 70–75. [Google Scholar]
  25. Feng, C.D.; Guo, J.H.; Nie, Y.C.; Wu, Z.B.; Zhang, X.L.; Zhang, J.F.; Stewart, J.M. Cytoplasmic-nuclear male sterility in cotton: Comparative RFLP analysis of mitochondrial DNA. In Proceedings of the 2000 Proceedings Beltwide Cotton Conferences, San Antonio, TX, USA, 4–8 January 2000. [Google Scholar]
  26. Wang, X.D. Analyses of mitochondrial protein and DNA from cytoplasmic male sterile cotton. Acta Agron. Sin. 2000, 26, 35–39. [Google Scholar]
  27. Li, S.S.; Chen, Z.W.; Zhao, N.; Wang, Y.M.; Nie, H.S.; Hua, J.P. The comparison of four mitochondrial genomes reveals cytoplasmic male sterility candidate genes in cotton. BMC Genom. 2018, 19, 775. [Google Scholar] [CrossRef]
  28. Khan, A.; Kong, X.; Liao, X.; Zheng, J.; You, J.; Li, M.; Hussain, R.M.; Raza, H.; Zhou, R. Mitochondrial gene expression analysis reveals aberrant transcription of cox3 in Gossypium barbadense CMS line H276A. Dev. Genes Evol. 2022, 232, 15–23. [Google Scholar] [CrossRef]
  29. You, J.Y.; Li, M.; Li, H.W.; Bai, Y.L.; Zhu, X.; Kong, X.J.; Chen, X.Y.; Zhou, R.Y. Integrated methylome and transcriptome analysis widen the knowledge of cytoplasmic male sterility in cotton (Gossypium barbadense L.). Front. Plant Sci. 2022, 13, 77098. [Google Scholar] [CrossRef]
  30. Heng, S.P.; Wei, C.; Jing, B.; Wan, Z.J.; Wen, J.; Yi, B.; Ma, C.Z.; Tu, J.X.; Fu, T.D.; Shen, J. Comparative analysis of mitochondrial genomes between the hau cytoplasmic male sterility (CMS) line and its iso-nuclear maintainer line in Brassica juncea to reveal the origin of the CMS-associated gene orf288. BMC Genom. 2014, 15, 322. [Google Scholar] [CrossRef]
  31. He, T.T.; Ding, X.L.; Zhang, H.; Li, Y.W.; Chen, L.F.; Wang, T.L.; Yang, L.S.; Nie, Z.X.; Song, Q.J.; Gai, J.; et al. Comparative analysis of mitochondrial genomes of soybean cytoplasmic male-sterile lines and their maintainer lines. Funct. Integr. Genom. 2021, 21, 43–57. [Google Scholar] [CrossRef] [PubMed]
  32. Liu, H.; Cui, P.; Zhan, K.; Lin, Q.; Zhuo, G.; Guo, X.; Ding, F.; Yang, W.; Liu, D.; Hu, S.; et al. Comparative analysis of mitochondrial genomes between a wheat K-type cytoplasmic male sterility (CMS) line and its maintainer line. BMC Genom. 2011, 12, 163. [Google Scholar] [CrossRef] [PubMed]
  33. Unseld, M.; Marienfeld, J.R.; Brandt, P.; Brennicke, A. The mitochondrial genome of Arabidopsis thaliana contains 57 genes in 366,924 nucleotides. Nat. Genet 1997, 15, 57–61. [Google Scholar] [CrossRef] [PubMed]
  34. Hu, J.; Huang, W.C.; Huang, Q.; Qin, X.J.; Yu, C.C.; Wang, L.L.; Li, S.Q.; Zhu, R.S.; Zhu, Y.G. Mitochondria and cytoplasmic male sterility in plants. Mitochondrion 2014, 19, 282–288. [Google Scholar] [CrossRef]
  35. Toriyama, K. Molecular basis of cytoplasmic male sterility and fertility restoration in rice. Plant Biotechnol. 2021, 38, 285–295. [Google Scholar] [CrossRef]
  36. Ivanov, M.K.; Dymshits, G.M. Cytoplasmic male sterility and restoration of pollen fertility in higher plants. Genetika 2007, 43, 451–468. [Google Scholar] [CrossRef]
  37. Sloan, D.B.; Alverson, A.J.; Štorchová, H.; Palmer, J.D.; Taylor, D.R. Extensive loss of translational genes in the structurally dynamic mitochondrial genome of the angiosperm Silene latifolia. BMC Evol. Biol. 2010, 10, 274. [Google Scholar] [CrossRef]
  38. Arrieta-Montiel, M.P.; Mackenzie, S.A. Plant mitochondrial genomes and recombination. In Plant Mitochondrial; Kempken, F., Ed.; Springer: New York, NY, USA, 2011; Volume 1, pp. 65–82. [Google Scholar]
  39. Tang, H.W.; Zheng, X.M.; Li, C.L.; Xie, X.R.; Chen, Y.L.; Chen, L.T.; Zhao, X.C.; Zheng, H.Q.; Zhou, J.J.; Ye, S.; et al. Multi-step formation, evolution, and functionalization of new cytoplasmic male sterility genes in the plant mitochondrial genomes. Cell Res. 2016, 27, 130–146. [Google Scholar] [CrossRef]
  40. Bentolila, S.; Stefanov, S. A reevaluation of rice mitochondrial evolution based on the complete sequence of male-fertile and male-sterile mitochondrial genomes. Plant Physiol. 2012, 158, 996–1017. [Google Scholar] [CrossRef]
  41. Zhang, Y.J.; Han, Y.; Zhang, M.; Zhang, X.X.; Guo, L.P.; Qi, T.X.; Li, Y.Q.; Feng, J.J.; Wang, H.L.; Tang, H.; et al. The cotton mitochondrial chimeric gene orf610a causes male sterility by disturbing the dynamic balance of ATP synthesis and ROS burst. Crop. J. 2022. [Google Scholar] [CrossRef]
  42. Kuwabara, K.; Harada, I.; Matsuzawa, Y.; Ariizumi, T.; Shirasawa, K. Organelle genome assembly uncovers the dynamic genome reorganization and cytoplasmic male sterility associated genes in tomato. Hortic Res. 2021, 8, 250–261. [Google Scholar] [CrossRef]
  43. Alverson, A.J.; Rice, D.W.; Dickinson, S.; Barry, K.; Palmer, J.D. Origins and recombination of the bacterial-sized multichromosomal mitochondrial genome of cucumber. Plant Cell 2011, 23, 2499–2513. [Google Scholar] [CrossRef]
  44. Mower, J.P.; Case, A.L.; Floro, E.R.; Willis, J.H. 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–686. [Google Scholar] [CrossRef]
  45. Tian, X.J.; Zheng, J.; Hu, S.N.; Yu, J. The rice mitochondrial genomes and their variations. Plant Physiol. 2006, 140, 401–410. [Google Scholar] [CrossRef]
  46. Park, J.Y.; Lee, Y.P.; Lee, J.; Choi, B.S.; Kim, S.; Yang, T.J. Complete mitochondrial genome sequence and identification of a candidate gene responsible for cytoplasmic male Sterility in Radish (Raphanus sativus L.) containing DCGMS cytoplasm. Theor. Appl. Genet. 2013, 126, 1763–1774. [Google Scholar] [CrossRef]
  47. Zabala, G.; Gabay-Laughnan, S.; Laughnan, J.R. The nuclear gene Rf3 affects the expression of the mitochondrial chimeric sequence R implicated in S-type male sterility in maize. Genetics 1997, 147, 847–860. [Google Scholar] [CrossRef]
  48. Sang, S.F.; Cheng, H.T.; Hao, M.Y.; Ding, B.L.; Mei, D.S.; Wang, H.; Wang, W.X.; Liu, J.; Fu, L.; Liu, K.; et al. Mitochondrial localization of ORF346 causes pollen abortion in alloplasmic male sterility. Crop. J. 2021, 9, 1320–1329. [Google Scholar] [CrossRef]
  49. Song, J.; Hedgcoth, C. A chimeric gene (orf256) is expressed as protein only in cytoplasmic male-sterile lines of wheat. Plant Mol. Biol. 1994, 26, 535–539. [Google Scholar] [CrossRef] [PubMed]
  50. Dahan, J.; Mireau, H. The Rf and Rf-like PPR in higher plants, a fast-evolving subclass of PPR genes. RNA Biol. 2013, 10, 1469–1476. [Google Scholar] [CrossRef] [PubMed]
  51. Bentolila, S.; Alfonso, A.A.; Hanson, M.R. A pentatricopeptide repeat-containing gene restores fertility to cytoplasmic male-sterile plants. Proc. Natl. Acad. Sci. USA 2002, 99, 10887–10892. [Google Scholar] [CrossRef]
  52. Brown, G.G.; Formanová, N.; Jin, H.; Wargachuk, R.; Dendy, C.; Patil, P.; Laforest, M.; Zhang, J.; Cheung, W.Y.; Landry, B.S. The radish Rfo restorer gene of ogura cytoplasmic male sterility encodes a protein with multiple pentatricopeptide repeats. Plant J. 2003, 35, 262–272. [Google Scholar] [CrossRef] [PubMed]
  53. Komori, T.; Ohta, S.; Murai, N.; Takakura, Y.; Kuraya, Y.; Suzuki, S.; Hiei, Y.; Imaseki, H.; Nitta, N. Map-based cloning of a fertility restorer gene, Rf-1, in rice (Oryza sativa L.). Plant J. 2004, 37, 315–325. [Google Scholar] [CrossRef]
  54. Klein, R.R.; Klein, P.E.; Mullet, J.E.; Minx, P.; Rooney, W.L.; Schertz, K.F. Fertility restorer locus Rf1 of sorghum (Sorghum bicolor L.) encodes a pentatricopeptide repeat protein not present in the colinear region of rice chromosome 12. Appl Genet. 2005, 111, 994–1012. [Google Scholar] [CrossRef]
  55. Barkan, A.; Small, I. Pentatricopeptide repeat proteins in plants. Annu Rev. Plant. Biol 2014, 65, 415–442. [Google Scholar] [CrossRef]
  56. Li, X.L.; Sun, M.D.; Liu, S.J.; Teng, Q.; Li, S.H.; Jiang, Y.S. Functions of PPR proteins in plant growth and development. Int. J. Mol. Sci. 2021, 22, 11274. [Google Scholar] [CrossRef]
  57. Huang, W.C.; Yu, C.C.; Hu, J.; Wang, L.L.; Dan, Z.W.; Zhou, W.; He, C.L.; Zeng, Y.F.; Yao, G.X.; Qi, J.Z.; et al. Pentatricopeptide-repeat family protein RF6 functions with hexokinase 6 to rescue rice cytoplasmic male sterility. Proc. Natl. Acad. Sci. USA 2015, 112, 14984–14989. [Google Scholar] [CrossRef]
  58. Falcon de Longevialle, A.F.; Meyer, E.H.; Andrés, C.; Taylor, N.L.; Lurin, C.; Millar, A.H.; Smalla, I.D. The pentatricopeptide repeat gene OTP43 is required for trans-splicing of the mitochondrial nad1 intron 1 in Arabidopsis thaliana. Plant Cell 2007, 19, 3256–3265. [Google Scholar] [CrossRef]
  59. Wang, C.D.; Blondel, L.; Quadrado, M.; Dargel-Graffin, C.; Mireau, H. Pentatricopeptide repeat protein MTSF3 ensures mitochondrial RNA stability and embryogenesis. Plant Physiol. 2022, 1–13. [Google Scholar] [CrossRef]
  60. Qin, X.E.; Tian, S.K.; Zhang, W.L.; Zheng, Q.; Wang, H.; Feng, Y.; Lin, Y.N.; Tang, J.H.; Wang, Y.; Yan, J.B.; 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, 1961–1964. [Google Scholar] [CrossRef]
  61. Zhang, D.Y.; Chen, C.; Wang, H.T.; Niu, E.; Zhao, P.Y.; Fang, S.; Zhu, G.Z.; Shang, X.G.; Guo, W.Z. Cotton fiber development requires the pentatricopeptide repeat protein GhIm for splicing of mitochondrial Nad7 mRNA. Genetics 2021, 217, 1–17. [Google Scholar] [CrossRef]
  62. Paterson, A.H.; Brubaker, C.L.; Wendel, J.F. A rapid method for extraction of cotton (Gossypium spp.) genomic DNA suitable for RFLP or PCR analysis. Plant Mol. Biol. Rep. 1993, 11, 122–127. [Google Scholar] [CrossRef]
  63. Hahn, C.; Bachmann, L.; Chevreux, B. Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads—A baiting and iterative mapping approach. Nucleic Acids Res. 2013, 41, e129. [Google Scholar] [CrossRef]
  64. Li, H. Minimap2: Pairwise alignment for nucleotide sequences. Bioinformatics 2018, 34, 3094–3100. [Google Scholar] [CrossRef]
  65. Kolmogorov, M.; Yuan, J.; Lin, Y.; Pevzner, P.A. Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol. 2019, 37, 540–546. [Google Scholar] [CrossRef] [PubMed]
  66. Walker, B.J.; Abeel, T.; Shea, T.; Priest, M.; Abouelliel, A.; Sakthikumar, S.; Cuomo, C.A.; Zeng, Q.; Wortman, J.; Young, S.K.; et al. Pilon: An integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 2014, 9, e112963. [Google Scholar] [CrossRef] [PubMed]
  67. Alverson, A.J.; Wei, X.; Rice, D.W.; Stern, D.B.; Barry, K.; Palmer, J.D. Insights into the evolution of mitochondrial genome size from complete sequences of Citrullus lanatus and Cucurbita pepo (Cucurbitaceae). Mol. Biol. Evol. 2010, 27, 1436–1448. [Google Scholar] [CrossRef]
  68. Lagesen, K.; Hallin, P.; Rødland, E.A.; Stærfeldt, H.H.; Rognes, T.; Ussery, D.W. RNAmmer: Consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007, 35, 3100–3108. [Google Scholar] [CrossRef]
  69. Kurtz, S.; Phillippy, A.; Delcher, A.L.; Smoot, M.; Shumway, M.; Antonescu, C.; Salzberg, S.L. Versatile and open software for comparing large genomes. Genome Biol. 2004, 5, 12. [Google Scholar] [CrossRef]
  70. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef]
  71. Li, H. Exploring single-sample SNP and InDel calling with whole-genome de novo assembly. Bioinformatics 2012, 28, 1838–1844. [Google Scholar] [CrossRef]
  72. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef]
  73. Wu, T.D.; Nacu, S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics 2010, 26, 873–881. [Google Scholar] [CrossRef]
  74. Wu, T.D.; Watanabe, C.K. GMAP: A genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics 2005, 21, 1859–1875. [Google Scholar] [CrossRef]
  75. Rao, X.Y.; Huang, X.L.; Zhou, Z.C.; Lin, X. An improvement of the 2ˆ(–delta delta CT) method for quantitative real-time polymerase chain reaction data analysis. Biostat. Bioinform. Biomath. 2013, 3, 71. [Google Scholar]
Ijms 23 09198 g001 550
Figure 1. Floral organs morphology and pollen shed of the CMS line, SI3A, and its restorer line, 0-613-2R. (A,C), The stamens and pistils of restorer line 0-613-2R (fertile) and CMS line SI3A (sterile). (B,D) The enlargement of the anther. (E,F) Pollen grains stained with I2-KI. Scale bars correspond to 3 cm (A,C), 0.5 cm (B,D), and 100 μm (E,F).
Figure 1. Floral organs morphology and pollen shed of the CMS line, SI3A, and its restorer line, 0-613-2R. (A,C), The stamens and pistils of restorer line 0-613-2R (fertile) and CMS line SI3A (sterile). (B,D) The enlargement of the anther. (E,F) Pollen grains stained with I2-KI. Scale bars correspond to 3 cm (A,C), 0.5 cm (B,D), and 100 μm (E,F).
Ijms 23 09198 g001
Ijms 23 09198 g002 550
Figure 2. Sequence collinearity analysis of 0-613-2R and SI3A mitochondrial genome with G. hirsutum and G. harknessii mitochondrial genome. Different colors represent different blocks in each mitochondrial genome and the two blocks that are connected by corresponding line indicates high homology between them. Reverse and direct transcript orientation are indicated above and below the central line.
Figure 2. Sequence collinearity analysis of 0-613-2R and SI3A mitochondrial genome with G. hirsutum and G. harknessii mitochondrial genome. Different colors represent different blocks in each mitochondrial genome and the two blocks that are connected by corresponding line indicates high homology between them. Reverse and direct transcript orientation are indicated above and below the central line.
Ijms 23 09198 g002
Ijms 23 09198 g003 550
Figure 3. Mitochondrial genome map of the restorer line, 0-613-2R, and the CMS line, SI3A. (A) Mitochondrial genome map of restorer line, 0-613-2R. (B) Mitochondrial genome map of CMS line, SI3A. Genes inside the loop represent transcription in a clockwise direction, while genes outside the loop are in the opposite direction. Different functional genes are marked with different colors. The built-in gray histogram shows the genomic GC content, and the middle gray line is the 50% threshold line.
Figure 3. Mitochondrial genome map of the restorer line, 0-613-2R, and the CMS line, SI3A. (A) Mitochondrial genome map of restorer line, 0-613-2R. (B) Mitochondrial genome map of CMS line, SI3A. Genes inside the loop represent transcription in a clockwise direction, while genes outside the loop are in the opposite direction. Different functional genes are marked with different colors. The built-in gray histogram shows the genomic GC content, and the middle gray line is the 50% threshold line.
Ijms 23 09198 g003
Ijms 23 09198 g004 550
Figure 4. Structural variation and InDels number in mitochondrial genomes of SI3A and 0-613-2R. (A) The structural variation in SI3A and 0-613-2R. Forward chain: forward chain of the genome sequence; Reverse chain: reverse chain of the genome sequence; Collinearity: the same linear region; Tran + Inver: the area of translocation and inversion; Inversion: the area of inversion; Translocation: the area of translocation. (B) The length of InDels and its corresponding number in SI3A mitochondrial genome compared with 0-613-2R. Yellow, insertion; Blue, deletion. (C) Amino acid sequences alignment of atp1 between 0-613-2R and SI3A.
Figure 4. Structural variation and InDels number in mitochondrial genomes of SI3A and 0-613-2R. (A) The structural variation in SI3A and 0-613-2R. Forward chain: forward chain of the genome sequence; Reverse chain: reverse chain of the genome sequence; Collinearity: the same linear region; Tran + Inver: the area of translocation and inversion; Inversion: the area of inversion; Translocation: the area of translocation. (B) The length of InDels and its corresponding number in SI3A mitochondrial genome compared with 0-613-2R. Yellow, insertion; Blue, deletion. (C) Amino acid sequences alignment of atp1 between 0-613-2R and SI3A.
Ijms 23 09198 g004
Ijms 23 09198 g005 550
Figure 5. Distribution of four types of ORFs and chimeric structures, transmembrane domain of the specific ORF in SI3A. The Venn diagram (A) shows the distribution of four types of ORFs. All ORFs: numbers of all the detected ORFs in SI3A; CM-ORFs: numbers of chimeric structure ORFs in SI3A; UR-ORFs: numbers of unique region ORFs in SI3A; TM-ORFs: numbers of transmembrane domain ORFs in SI3A. (B) Chimeric structure of orf606a, orf131a, orf119a-2, and orf109a. Rectangle boxes represent coding sequences, different colors represent different coding genes, linear represents the different size of the noncoding sequences, and U represents the unknown sequences. (C–F) The location and probabilities of transmembrane domains of the four specific ORFs orf606a (C), orf131a (D), orf119a-2 (E), and orf109a (F).
Figure 5. Distribution of four types of ORFs and chimeric structures, transmembrane domain of the specific ORF in SI3A. The Venn diagram (A) shows the distribution of four types of ORFs. All ORFs: numbers of all the detected ORFs in SI3A; CM-ORFs: numbers of chimeric structure ORFs in SI3A; UR-ORFs: numbers of unique region ORFs in SI3A; TM-ORFs: numbers of transmembrane domain ORFs in SI3A. (B) Chimeric structure of orf606a, orf131a, orf119a-2, and orf109a. Rectangle boxes represent coding sequences, different colors represent different coding genes, linear represents the different size of the noncoding sequences, and U represents the unknown sequences. (C–F) The location and probabilities of transmembrane domains of the four specific ORFs orf606a (C), orf131a (D), orf119a-2 (E), and orf109a (F).
Ijms 23 09198 g005
Ijms 23 09198 g006 550
Figure 6. Different expression levels of mitochondrial genes and schematic diagram of nad4 intron splicing patterns. (A) Relative expression levels of mitochondrial genes with significant differences in RNA−seq data. Blue represents significantly downregulated genes and red represents significantly upregulated genes in SI3A. (B) RT−qPCR analysis of four specific ORFs expression in CMS line SI3A and restorer line 0-613-2R, GhHis3 was used as the internal control. (C) RT−PCR analysis of the four specific ORFs in 0-613-2R and SI3A. (D) Schematic diagram of nad4 intron splicing patterns, two transcripts formed after splicing retention of intron 1 and intron 3 are shown on the right.
Figure 6. Different expression levels of mitochondrial genes and schematic diagram of nad4 intron splicing patterns. (A) Relative expression levels of mitochondrial genes with significant differences in RNA−seq data. Blue represents significantly downregulated genes and red represents significantly upregulated genes in SI3A. (B) RT−qPCR analysis of four specific ORFs expression in CMS line SI3A and restorer line 0-613-2R, GhHis3 was used as the internal control. (C) RT−PCR analysis of the four specific ORFs in 0-613-2R and SI3A. (D) Schematic diagram of nad4 intron splicing patterns, two transcripts formed after splicing retention of intron 1 and intron 3 are shown on the right.
Ijms 23 09198 g006
Table
Table 1. Features of the assembly mitogenomes in the CMS line, SI3A, and its restorer line, 0-613-2R.
Table 1. Features of the assembly mitogenomes in the CMS line, SI3A, and its restorer line, 0-613-2R.
Genome Characteristics0613-2R mtSI3A
SI3A mt1SI3A mt2SI3A mt3
Genomic size(bp)607,36763,927396,206174,074
G + C content (%)44464445
coding sequence (%)25.3%19.0%7.4%8.4%
ORF48123314
Protein coding genes336207
tRNA genes272185
rRNA genes6 13 11 12 1
Repeat content percent coverage of total genome16.2%0.8%8.8%9.2%
Large repeats: >1 kb(number)13050
Small repeats: <1 kb(number)4271116934
1 Multiple copies of rRNAs.
Table
Table 2. Repeat sequences (>50 bp) in the mitogenomes of CMS line SI3Amt2.
Table 2. Repeat sequences (>50 bp) in the mitogenomes of CMS line SI3Amt2.
Copy1 Copy2
NoSize (bp)Start1Stop1TypeSize (bp)Start2Stop2Identity (%)Difference between Copies
AR110,637120,585131,221DR 110,636359,345369,98099copy1 9 bp indel, copy2 11 bp indel.
AR2435237,794238,228DR435334,087334,521100
AR3333211,333211,665IR 2333360,536360,868100
AR4229154,125154,353DR229254,009254,237100
AR5203237,590237,792DR203333,883334,085100
AR617583257IR175281,969282,143100
AR715456,02756,180DR154359,760359,913100
AR812156,02756,147DR121121,000121,120100
AR911340,19240,304DR113102,378102,490100
AR1011357,14157,253DR113212,119212,231100
AR11107201,247201,353DR107389,303389,409100
AR1293111,326111,418DR93331,353331,445100
AR137875,29375,370IR78358,615358,692100
AR1471264,143264,213IR71360,450360,520100
AR156675,14675,211IR66116,361116,426100
AR1666211,392211,457IR66297,562297,627100
AR1766297,562297,627DR66360,744360,809100
AR186040,30640,365DR60102,492102,551100
AR195772527308DR57324,694324,750100
AR205347,23347,285IR53170,558170,610100
1 DR: direct repeats, 2 IR: invert repeats.
Table
Table 3. The annotation of SNP in the mitogenomes of SI3A compared with 0-613-2R.
Table 3. The annotation of SNP in the mitogenomes of SI3A compared with 0-613-2R.
GeneLocation0-613-2R<->SI3A (Nucleic Acid)0-613-2R<->SI3A (Amino Acid)Mutation Type 1SNP Type
cox1495,890C<->AIle<->IleStransversion
cox1496,358C<->AIle<->IleStransversion
cox3494,152C<->ALeu<->IleNtransversion
nad7583,680A<->CIle<->IleStransversion
atp4307,826C<->TPhe<->PheStransition
atp8205,705C<->ASer<->ArgNtransversion
sdh3317,242A<->CLeu<->PheNtransversion
matR16,075C<->ACys<->LysNtransversion
rps4426,548A<->CLys<->CysNtransversion
rpl2458,973T<->GPhe<->LeuNtransversion
rpl2459,220A<->CIIe<->LeuNtransversion
rpl5457,988A<->CLys<->CysNtransversion
rpl10171,493G<->AGln<->LysNtransition
rpl16399,033A<->CVal<->ValStransversion
1 S: synonymous, N: non-synonymous.
Table
Table 4. Mapping results for Illumina RNA-seq data to mitochondrial genome.
Table 4. Mapping results for Illumina RNA-seq data to mitochondrial genome.
Sample0-613-2R-10-613-2R-20-613-2R-3SI3A-1SI3A-2SI3A-3
Total raw reads27,570,39827,645,30027,163,44924,396,97130,170,97822,451,704
Total clean reads26,157,63426,560,99626,094,84723,354,55128,705,90821,508,933
Clean reads Q20(%)97.3897.497.2297.3397.1397.38
Clean reads Q30(%)92.5992.6692.2592.4492.0892.59
Total alignment with mt genome25,765 (0.10%)25,125 (0.09%)26,713 (0.10%)16,247 (0.07%)21,734 (0.08%)18,759 (0.09%)
Exact alignment with mt genome14,708 (0.06%)14,376 (0.05%)15,446 (0.06%)12,592 (0.05%)16,988 (0.06%)14,509 (0.07%)
Average (exact)0.06%0.06%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xuan, L.; Qi, G.; Li, X.; Yan, S.; Cao, Y.; Huang, C.; He, L.; Zhang, T.; Shang, H.; Hu, Y. Comparison of Mitochondrial Genomes between a Cytoplasmic Male-Sterile Line and Its Restorer Line for Identifying Candidate CMS Genes in Gossypium hirsutum. Int. J. Mol. Sci. 2022, 23, 9198. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23169198

AMA Style

Xuan L, Qi G, Li X, Yan S, Cao Y, Huang C, He L, Zhang T, Shang H, Hu Y. Comparison of Mitochondrial Genomes between a Cytoplasmic Male-Sterile Line and Its Restorer Line for Identifying Candidate CMS Genes in Gossypium hirsutum. International Journal of Molecular Sciences. 2022; 23(16):9198. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23169198

Chicago/Turabian Style

Xuan, Lisha, Guoan Qi, Xiaoran Li, Sunyi Yan, Yiwen Cao, Chujun Huang, Lu He, Tianzhen Zhang, Haihong Shang, and Yan Hu. 2022. "Comparison of Mitochondrial Genomes between a Cytoplasmic Male-Sterile Line and Its Restorer Line for Identifying Candidate CMS Genes in Gossypium hirsutum" International Journal of Molecular Sciences 23, no. 16: 9198. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23169198

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop