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

The miR3367-lncRNA67-GhCYP724B module regulates male sterility by modulating brassinosteroid biosynthesis and interacting with Aorf27 in Gossypium hirsutum

Cytoplasmic male sterile (CMS) lines play a crucial role in utilization of heterosis in crop plants. However, the mechanism underlying the manipulation of male sterility in cotton by long non-coding RNA (lncRNA) and brassinosteroids (BRs) remains elusive. Here, using an integrative approach combining lncRNA transcriptomic profiles with virus-induced gene silencing experiments, we identify a flower bud-specific lncRNA in the maintainer line 2074B, lncRNA67, negatively modulating with male sterility in upland cotton (Gossypium hirsutum). lncRNA67 positively regulates cytochrome P274B (GhCYP724B), which acted as an eTM (endogenous target mimic) for miR3367. The suppression of GhCYP724B induced symptoms of BR deficiency and male semi-sterility in upland cotton as well as in tobacco, which resulted from a reduction in the endogenous BR contents. GhCYP724B regulates BRs synthesis by interacting with GhDIM and GhCYP90B, two BRs biosynthesis proteins. Additionally, GhCYP724B suppressed a unique chimeric open reading frame (Aorf27) in 2074A mitochondrial genome. Ectopic expression of Aorf27 in yeast inhibited cellular growth, and over expression of Aorf27 in tobacco showed male sterility. Overall, the results proved that the miR3367-lncRNA67-GhCYP724B module positively regulates male sterility by modulating BRs biosynthesis. The findings uncovered the function of lncRNA67-GhCYP724B in male sterility, providing a new mechanism for understanding male sterility in upland cotton.

Orf116b Induces Pollen Abortion in a Novel Cotton (Gossypium hirsutum L.) Cytoplasmic Male Sterile Line J4A

Cytoplasmic male sterile (CMS) systems have great potential in hybrid seed production in many plants. However, the incompatibility between the cytoplasmic and nuclear genes and the availability of fewer CMS lines limit the system's application in cotton heterosis. Therefore, the present study sequenced the mitochondrial (mt) genomes of a novel cotton (Gossypium hirsutum L.) CMS line J4A and its cytoplasmic nuclear homologous maintainer line J4B to investigate the mechanism underlying CMS and improve its application. A novel CMS gene, orf116b, was identified in the CMS line. Ectopic expression of orf116b in Escherichia coli suppressed growth, while its overexpression in Arabidopsis, rice, tobacco, and cotton led to complete or partial male sterility. Further analysis of anthers revealed mt disruption in J4A with higher levels of hydrogen peroxide (H2O2) and malondialdehyde (MDA), but lower levels of ATP and ribosomal protein (RP) than in J4B. Finally, a membrane-based yeast two-hybrid (MYTH) assay and bimolecular fluorescence complementation (BiFC) assays demonstrated that orf116b probably interacts with an anther-specific protein, LAT52. These observations collectively proved that orf116b is associated with early and stable pollen abortion in cotton, providing a foundation for further research on cotton fertility restoration and heterosis breeding.

Cytoplasmic male sterility-based hybrids: mechanistic insights

MAIN CONCLUSION

A comprehensive understanding of the nucleocytoplasmic interactions that occur between genes related to the restoration of fertility and cytoplasmic male sterility (CMS) provides insight into the development of hybrids of important crop species. Modern biotechnological techniques allow this to be achieved in an efficient and quick manner. Heterosis is paramount for increasing the yield and quality of a crop. The development of hybrids for achieving heterosis has been well-studied and proven to be robust and efficient. Cytoplasmic male sterility (CMS) has been explored extensively in the production of hybrids. The underlying mechanisms of CMS include the role of cytotoxic proteins, PCD of tapetal cells, and improper RNA editing of restoration factors. On the other hand, the restoration of fertility is caused by the presence of restorer-of-fertility (Rf) genes or restorer genes, which inhibit the effects of sterility-causing genes. The interaction between mitochondria and the nuclear genome is crucial for several regulatory pathways, as observed in the CMS-Rf system and occurs at the genomic, transcriptional, post-transcriptional, translational, and post-translational levels. These CMS-Rf mechanisms have been validated in several crop systems. This review aims to summarize the nucleo-mitochondrial interaction mechanism of the CMS-Rf system. It also sheds light on biotechnological interventions, such as genetic engineering and genome editing, to achieve CMS-based hybrids.

Deciphering the complex cotton genome for improving fiber traits and abiotic stress resilience in sustainable agriculture

BACKGROUND

Understanding the complex cotton genome is of paramount importance in devising a strategy for sustainable agriculture. Cotton is probably the most economically important cash crop known for its cellulose-rich fiber content. The cotton genome has become an ideal model for deciphering polyploidization due to its polyploidy, setting it apart from other major crops. However, the main challenge in understanding the functional and regulatory functions of many genes in cotton is still the complex cotton polyploidy genome, which is not limited to a single role. Cotton production is vulnerable to the sensitive effects of climate change, which can alter or aggravate soil, pests, and diseases. Thus, conventional plant breeding coupled with advanced technologies has led to substantial progress being made in cotton production.

GENOMICS APPROACHES IN COTTON

In the frontier areas of genomics research, cotton genomics has gained momentum accomplished by robust high-throughput sequencing platforms combined with novel computational tools to make the cotton genome more tractable. Advances in long-read sequencing have allowed for the generation of the complete set of cotton gene transcripts giving incisive scientific knowledge in cotton improvement. In contrast, the integration of the latest sequencing platforms has been used to generate multiple high-quality reference genomes in diploid and tetraploid cotton. While pan-genome and 3D genomic studies are still in the early stages in cotton, it is anticipated that rapid advances in sequencing, assembly algorithms, and analysis pipelines will have a greater impact on advanced cotton research.

CONCLUSIONS

This review article briefly compiles substantial contributions in different areas of the cotton genome, which include genome sequencing, genes, and their molecular regulatory networks in fiber development and stress tolerance mechanism. This will greatly help us in understanding the robust genomic organization which in turn will help unearth candidate genes for functionally important agronomic traits.

ROS accumulation-induced tapetal PCD timing changes leads to microspore abortion in cotton CMS lines

BACKGROUND

Cytoplasmic male sterility (CMS) is the basis of heterosis exploitation. CMS has been used to hybrid production in cotton, but its molecular mechanism remains unclear. CMS is associated with advanced or delayed tapetal programmed cell death (PCD), and reactive oxygen species (ROS) may mediate this process. In this study, we obtained Jin A and Yamian A, two CMS lines with different cytoplasmic sources.

RESULTS

Compared with maintainer Jin B, Jin A anthers showed advanced tapetal PCD with DNA fragmentation, producing excessive ROS which accumulated around the cell membrane, intercellular space and mitochondrial membrane. The activities of peroxidase (POD) and catalase (CAT) enzymes which can scavenge ROS were significantly decreased. However, Yamian A tapetal PCD was delayed with lower ROS content, and the activities of superoxide dismutase (SOD) and POD were higher than its maintainer. These differences in ROS scavenging enzyme activities may be caused by isoenzyme gene expressions. In addition, we found the excess ROS generated in Jin A mitochondria and ROS overflow from complex III might be the source in parallel with the reduction of ATP content.

CONCLUSION

ROS accumulation or abrogation were mainly caused by the joint action of ROS generation and scavenging enzyme activities transformation, which led to the abnormal progression of tapetal PCD, affected the development of microspores, and eventually contributed to male sterility. In Jin A, tapetal PCD in advance might be caused by mitochondrial ROS overproduction, accompanied by energy deficiency. The above studies will provide new insights into the cotton CMS and guide the follow-up research ideas.

Identification of fertility restoration candidate genes from a restorer line R186 for Gossypium harknessii cytoplasmic male sterile cotton

BACKGROUND

The utilization of heterosis based on three-line system is an effective strategy in crop breeding. However, cloning and mechanism elucidation of restorer genes for cytoplasmic male sterility (CMS) in upland cotton have yet been realized.

RESULTS

This research is based on CMS line 2074A with the cytoplasm from Gossypium harknessii (D_2-2) and restorer line R186. The offspring of 2074A × R186 were used to conduct genetic analysis. The fertility mechanism of 2074A can be speculated to be governed by multiple genes, since neither the single gene model nor the double genes model could be used. The bulked segregant analysis (BSA) for (2074A × R186) F₂ determined the genetic interval of restorer genes on a region of 4.30 Mb on chromosome D05 that contains 77 annotated genes. Four genes were identified as candidates for fertility restoration using the RNA-seq data of 2074A, 2074B, and R186. There are a number of large effect variants in the four genes between 2074A and R186 that could cause amino acid changes. Evolutionary analysis and identity analysis revealed that GH_D05G3183, GH_D05G3384, and GH_D05G3490 have high identity with their homologs in D_2-2, respectively. Tissue differential expression analysis revealed that the genes GH_D05G3183, GH_D05G3384, and GH_D05G3490 were highly expressed in the buds of the line R186. The predicted results demonstrated that GH_D05G3183, GH_D05G3384 and GH_D05G3490 might interact with GH_A02G1295 to regulate orf610a in mitochondria.

CONCLUSION

Our study uncovered candidate genes for fertility restoration in the restorer line R186 and predicted the possible mechanism for restoring the male fertility in 2074A. This research provided valuable insight into the nucleoplasmic interactions.

Cotton heterosis and hybrid cultivar development

Cotton, the most important economic crop in the world, displays strong hybrid vigor, and has long been a subject for hybrid cultivar breeding. Here, advances in the theoretical and applied research in cotton heterosis along with its hybrid cultivar development by hand-emasculation and pollination (HEP), cytoplasmic (CMS) and genic male sterile lines (GMS) mainly in China during the past few decades are presented in this review. Three types of hybrids produced by HEP, CMS and GMS facilitating hybrid seed production with hand-pollination have been developed and are being planted simultaneously in cotton production. However, most hybrids commercially planted in production are produced by HEP, therefore, F₂ seeds are being extensively planted due to the high cost to produce F₁ seed. F₂ generations of these combinations exceed the check cultivars in yield usually up to 5~15%. GMS genes (ms₂ and ms₅ms₆) used in hybrid seed production and casual mitochondrial genes for G. harknessii CMS have been cloned. Challenges and opportunities in cotton heterosis and future hybrid cultivar development in cotton are discussed.

Plant non-coding RNAs function in pollen development and male sterility

Male sterility is classified as either cytoplasmic male sterility (CMS) or genic male sterility (GMS). Generally, CMS involves mitochondrial genomes interacting with the nuclear genome, while GMS is caused by nuclear genes alone. Male sterility is regulated by multilevel mechanisms in which non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and phased small interfering RNAs (phasiRNAs), which have been proven to be critical elements. The development of high-throughput sequencing technology offers new opportunities to evaluate the genetic mechanism of ncRNAs in plant male sterility. In this review, we summarize the critical ncRNAs that regulate gene expression in ways dependent on or independent of hormones, which involve the differentiation of the stamen primordia, degradation of the tapetum, formation of microspores, and the release of pollen. In addition, the key mechanisms of the miRNA-lncRNA-mRNA interaction networks mediating male sterility in plants are elaborated. We present a different perspective on exploring the ncRNA-mediated regulatory pathways that control CMS in plants and create male-sterile lines through hormones or genome editing. A refined understanding of the ncRNA regulatory mechanisms in plant male sterility for the development of new sterile lines would be conducive to improve hybridization breeding.

Polymeric structure of the Cannabis sativa L. mitochondrial genome identified with an assembly graph model

Cannabis sativa L. belongs to the family Cannabaceae in Rosales. It has been widely used as medicines, building materials, and textiles. Elucidating its genome is critical for molecular breeding and synthetic biology study. Many studies have shown that the mitochondrial genomes (mitogenomes) and even chloroplast genomes (plastomes) had complex polymeric structures. Using the Nanopore sequencing platform, we sequenced, assembled, and analyzed its mitogenome and plastome. The resulting unitig graph suggested that the mitogenome had a complex polymeric structure. However, a gap-free, circular sequence was further assembled from the unitig graph. In contrast, a circular sequence representing the plastome was obtained. The mitogenome major conformation was 415,837 bp long, and the plastome was 153,927 bp long. To test if the repeat sequences promote recombination, which corresponds to the branch points in the structure, we tested the sequences around repeats by long-read mapping. Among 208 pairs of predicted repeats, the mapping results supported the presence of cross-over around 25 pairs of repeats. Subsequent PCR amplification confirmed the presence of cross-over around 15 of the 25 repeats. By comparing the mitogenome and plastome sequences, we identified 19 mitochondria plastid DNAs, including seven complete genes (trnW-CCA, trnP-UGG, psbJ, trnN-GUU, trnD-GUC, trnH-GUG, trnM-CAU) and nine gene fragments. Furthermore, the selective pressure analysis results showed that five genes (atp1, ccmB, ccmC, cox1, nad7) had 19 positively selected sites. Lastly, we predicted 28 RNA editing sites. A total of 8 RNA editing sites located in the coding regions were successfully validated by PCR amplification and Sanger sequencing, of which four were synonymous, and four were nonsynonymous. In particular, the RNA editing events appeared to be tissue-specific in C. sativa mitogenome. In summary, we have confirmed the major confirmation of C. sativa mitogenome and characterized its structural features in detail. These results provide critical information for future variety breeding and resource development for C. sativa.

Homology-based identification of candidate genes for male sterility editing in upland cotton (Gossypium hirsutum L.)

Upland cotton (Gossypium hirsutum L.) accounts for more than 90% of the world's cotton production, providing natural material for the textile and oilseed industries worldwide. One strategy for improving upland cotton yields is through increased adoption of hybrids; however, emasculation of cotton flowers is incredibly time-consuming and genetic sources of cotton male sterility are limited. Here we review the known biochemical modes of plant nuclear male sterility (NMS), often known as plant genetic male sterility (GMS), and characterized them into four groups: transcriptional regulation, splicing, fatty acid transport and processing, and sugar transport and processing. We have explored protein sequence homology from 30 GMS genes of three monocots (maize, rice, and wheat) and three dicots (Arabidopsis, soybean, and tomato). We have analyzed evolutionary relationships between monocot and dicot GMS genes to describe the relative similarity and relatedness of these genes identified. Five were lowly conserved to their source species, four unique to monocots, five unique to dicots, 14 highly conserved among all species, and two in the other category. Using this source, we have identified 23 potential candidate genes within the upland cotton genome for the development of new male sterile germplasm to be used in hybrid cotton breeding. Combining homology-based studies with genome editing may allow for the discovery and validation of GMS genes that previously had no diversity observed in cotton and may allow for development of a desirable male sterile mutant to be used in hybrid cotton production.

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

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.

Integrated Methylome and Transcriptome Analysis Widen the Knowledge of Cytoplasmic Male Sterility in Cotton (Gossypium barbadense L.)

DNA methylation is defined as a conserved epigenetic modification mechanism that plays a key role in maintaining normal gene expression without altering the DNA sequence. Several studies have reported that altered methylation patterns were associated with male sterility in some plants such as rice and wheat, but global methylation profiles and their possible roles in cytoplasmic male sterility (CMS), especially in cotton near-isogenic lines, remain unclear. In this study, bisulfite sequencing technology and RNA-Seq were used to investigate CMS line 07-113A and its near-isogenic line 07-113B. Using integrated methylome and transcriptome analyses, we found that the number of hypermethylated genes in the differentially methylated regions, whether in the promoter region or in the gene region, was more in 07-113A than the number in 07-113B. The data indicated that 07-113A was more susceptible to methylation. In order to further analyze the regulatory network of male sterility, transcriptome sequencing and DNA methylation group data were used to compare the characteristics of near-isogenic lines 07-113A and 07-113B in cotton during the abortion stage. Combined methylation and transcriptome analysis showed that differentially expressed methylated genes were mainly concentrated in vital metabolic pathways including the starch and sucrose metabolism pathways and galactose metabolism. And there was a negative correlation between gene methylation and gene expression. In addition, five key genes that may be associated with CMS in cotton were identified. These data will support further understanding of the effect of DNA methylation on gene expression and their potential roles in cotton CMS.

Organelle genome assembly uncovers the dynamic genome reorganization and cytoplasmic male sterility associated genes in tomato

To identify cytoplasmic male sterility (CMS)-associated genes in tomato, we determined the genome sequences of mitochondria and chloroplasts in three CMS tomato lines derived from independent asymmetric cell fusions, their nuclear and cytoplasmic donors, and male fertile weedy cultivated tomato and wild relatives. The structures of the CMS mitochondrial genomes were highly divergent from those of the nuclear and cytoplasmic donors, and genes of the donors were mixed up in these genomes. On the other hand, the structures of CMS chloroplast genomes were moderately conserved across the donors, but CMS chloroplast genes were unexpectedly likely derived from the nuclear donors. Comparative analysis of the structures and contents of organelle genes and transcriptome analysis identified three genes that were uniquely present in the CMS lines, but not in the donor or fertile lines. RNA-sequencing analysis indicated that these three genes transcriptionally expressed in anther, and identified different RNA editing levels in one gene, orf265, that was partially similar to ATP synthase subunit 8, between fertile and sterile lines. The orf265 was a highly potential candidate for CMS-associated gene. This study suggests that organelle reorganization mechanisms after cell fusion events differ between mitochondria and chloroplasts, and provides insight into the development of new F1 hybrid breeding programs employing the CMS system in tomato.

Unraveling the Genetic Basis of Fertility Restoration for Cytoplasmic Male Sterile Line WNJ01A Originated From Brassica juncea in Brassica napus

Crosses that lead to heterosis have been widely used in the rapeseed (Brassica napus L.) industry. Cytoplasmic male sterility (CMS)/restorer-of-fertility (Rf) systems represent one of the most useful tools for rapeseed production. Several CMS types and their restorer lines have been identified in rapeseed, but there are few studies on the mechanisms underlying fertility restoration. Here, we performed morphological observation, map-based cloning, and transcriptomic analysis of the F₂ population developed by crossing the CMS line WNJ01A with its restorer line Hui01. Paraffin-embedded sections showed that the sporogenous cell stage was the critical pollen degeneration period, with major sporogenous cells displaying loose and irregular arrangement in sterile anthers. Most mitochondrial electron transport chain (mtETC) complex genes were upregulated in fertile compared to sterile buds. Using bulked segregant analysis (BSA)-seq to analyze mixed DNA pools from sterile and fertile F₂ buds, respectively, we identified a 6.25 Mb candidate interval where Rfw is located. Using map-based cloning experiments combined with bacterial artificial chromosome (BAC) clone sequencing, the candidate interval was reduced to 99.75 kb and two pentatricopeptide repeat (PPR) genes were found among 28 predicted genes in this interval. Transcriptome sequencing showed that there were 1679 DEGs (1023 upregulated and 656 downregulated) in fertile compared to sterile F₂ buds. The upregulated differentially expressed genes (DEGs) were enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG) lysine degradation pathway and phenylalanine metabolism, and the downregulated DEGs were enriched in cutin, suberine, and wax biosynthesis. Furthermore, 44 DEGs were involved in pollen and anther development, such as tapetum, microspores, and pollen wall development. All of them were upregulated except a few such as POE1 genes (which encode Pollen Ole e I allergen and extensin family proteins). There were 261 specifically expressed DEGs (9 and 252 in sterile and fertile buds, respectively). Regarding the fertile bud-specific upregulated DEGs, the ubiquitin-proteasome pathway was enriched. The top four hub genes in the protein-protein interaction network (BnaA09g56400D, BnaA10g18210D, BnaA10g18220D, and BnaC09g41740D) encode RAD23d proteins, which deliver ubiquitinated substrates to the 26S proteasome. These findings provide evidence on the pathways regulated by Rfw and improve our understanding of fertility restoration.

Transcriptome and MiRNAomics Analyses Identify Genes Associated with Cytoplasmic Male Sterility in Cotton (Gossypium hirsutum L.)

Cytoplasmic male sterility (CMS) is important for large-scale hybrid seed production. Rearrangements in the mitochondrial DNA (mtDNA) for the cotton (Gossypium hirsutum L.) CMS line J4A were responsible for pollen abortion. However, the expression patterns of nuclear genes associated with pollen abortion and the molecular basis of CMS for J4A are unknown, and were the objectives of this study by comparing J4A with the J4B maintainer line. Cytological evaluation of J4A anthers showed that microspore abortion occurs during meiosis preventing pollen development. Changes in enzyme activity of mitochondrial respiratory chain complex IV and mitochondrial respiratory chain complex V and the content of ribosomal protein and ATP during anther abortion were observed for J4A suggesting insufficient synthesis of ATP hindered pollen production. Additionally, levels of sucrose, starch, soluble sugar, and fructose were significantly altered in J4A during the meiosis stage, suggesting reduced sugar metabolism contributed to sterility. Transcriptome and miRNAomics analyses identified 4461 differentially expressed mRNAs (DEGs) and 26 differentially expressed microRNAs (DEMIs). Pathway enrichment analysis indicated that the DEMIs were associated with starch and sugar metabolism. Six deduced target gene regulatory pairs that may participate in CMS were identified, ghi-MIR7484-10/mitogen-activated protein kinase kinase 6 (MAPKK6), ghi-undef-156/agamous-like MADS-box protein AGL19 (AGL19), ghi-MIR171-1-22/SNF1-related protein kinase regulatory subunit gamma-1 and protein trichome birefringence-like 38, and ghi-MIR156-(8/36)/WRKY transcription factor 28 (WRKY28). Overall, a putative CMS mechanism involving mitochondrial dysfunction, the ghi-MIR7484-10/MAPKK6 network, and reduced glucose metabolism was suggested, and ghi-MIR7484-10/MAPKK6 may be related to abnormal microspore meiosis and induction of excessive sucrose accumulation in anthers.

Comparative analysis of mitochondrial genomes of soybean cytoplasmic male-sterile lines and their maintainer lines

In soybean, only one mitochondrial genome of cultispecies has been completely obtained. To explore the effect of mitochondrial genome on soybean cytoplasmic male sterility (CMS), two CMS lines and three maintainer lines were used for sequencing. Comparative analysis showed that mitochondrial genome of the CMS line was more compact than that of its maintainer line, but genes were highly conserved. Conserved and unique sequence coexisted in the genomes. Mitochondrial genomes contained different sequence lengths and copy numbers of repeats between CMS line and maintainer line. Large and short repeats mediated intramolecular and intermolecular recombination in mitochondria. Unique sequences and genes were also involved in recombination process and constituted a complex network. orf178 and orf261 were identified as CMS-associated candidate genes. They had sequence characteristics of reported CMS genes in other crops and could be transcribed in CMS lines but not in maintainer lines. This report reveals mitochondrial genome of soybean CMS lines and compares complete mitochondrial sequence between CMS lines and their maintainer lines. The information will be helpful in further understanding the characteristics of soybean mitochondrial genome and the mechanism underlying CMS.

The Mitochondrial Genome Assembly of Fennel (Foeniculum vulgare) Reveals Two Different atp6 Gene Sequences in Cytoplasmic Male Sterile Accessions

Cytoplasmic male sterility (CMS) has always aroused interest among researchers and breeders, being a valuable resource widely exploited not only to breed F1 hybrid varieties but also to investigate genes that control stamen and pollen development. With the aim of identifying candidate genes for CMS in fennel, we adopted an effective strategy relying on the comparison between mitochondrial genomes (mtDNA) of both fertile and sterile genotypes. mtDNA raw reads derived from a CMS genotype were assembled in a single molecule (296,483 bp), while a draft mtDNA assembly (166,124 nucleotides, 94 contigs) was performed using male fertile sample (MF) sequences. From their annotation and alignment, two atp6-like sequences were identified. atp6-, the putative mutant copy with a 300 bp truncation at the 5'-end, was found only in the mtDNA of CMS samples, while the wild type copy (atp6+) was detected only in the MF mtDNA. Further analyses (i.e., reads mapping and Sanger sequencing), revealed an atp6+ copy also in CMS samples, probably in the nuclear DNA. However, qPCRs performed on different tissues proved that, despite its availability, atp6+ is expressed only in MF samples, while apt6- mRNA was always detected in CMS individuals. In the light of these findings, the energy deficiency model could explain the pollen deficiency observed in male sterile flower. atp6- could represent a gene whose mRNA is translated into a not-fully functional protein leading to suboptimal ATP production that guarantees essential cellular processes but not a high energy demand process such as pollen development. Our study provides novel insights into the fennel mtDNA genome and its atp6 genes, and paves the way for further studies aimed at understanding their functional roles in the determination of male sterility.

Mitochondrial proteomic analysis reveals that proteins relate to oxidoreductase activity play a central role in pollen fertility in cotton

Cotton (Gossypium hirsutum L.) is an important economic crop. Cytoplasm male sterility (CMS) has been used to develop hybrid system and to produce hybrid seeds in cotton, but the molecular mechanism of CMS remains unclear. Mitochondria are semi-autonomous organelles, which play an important role in the reproduction of flowering plants. Male sterility has been proved associated with mitochondrial dysfunction in plants. In present study, a new strategy of proteomic sequencing data-independent acquisition (DIA) was used to analysis protein abundance across CMS lines 2074A (cytoplasm of Gossypium harknessii, D_2-2) and 2074S (cytoplasm of G. hirsutum, AD₁), and their maintainer 2074B. Comparing with transcriptome results showed that there is little consistence between proteome and transcriptome. A total of 2095 protein species were identified in three materials, and 186 and 161 differentially proteins were detected in the comparisons of 2074A vs 2074B, and 2074S vs 2074B, respectively. Among them, 49 and 50 proteins were specific existed in anther, and mainly participated in oxidoreductase activity, carbohydrate metabolism, fatty acid metabolism, cell aging, wax or cutin deposition and signal transduction. Gh_A07G0770 and Gh_D05G1908 were specific up-regulated in sterility lines, and the other genes Gh_D08G1196, Gh_D12G1971, Gh_A11G1250, Gh_D08G0388 were down-regulated, which presented similar expression tendency verified by qRT-PCR, transcriptome and proteome results. These six genes related to lipid synthesis, response to oxidative stress and cell aging, suggested them being involved in CMS occurrence. Using virus-induced gene silencing (VIGS) system, sterility obtained demonstrated the silencing Gh_A11G1250 in maintainer 2074B led to partial anthers abortion. Gh_A11G1250 encoded a mitochondrial localization of peroxisomal-like protein, participated in response to reactive oxygen species (ROS). Twenty-two proteins interacting with Gh_A11G1250 mainly related to chlorophyll biosynthetic process, photoperiodism and flowering, which showed different expression pattern between the male sterile line 2074A and maintainer 2074B. This novel research based on mitochondrial proteomics comparison confirmed that DAPs related to oxidative stress are critical to pollen abortion. BIOLOGICAL SIGNIFICANCE: Cytoplasm male sterility (CMS) system is utilized widely for hybrid production in cotton. However, the genetic and molecular mechanisms of CMS still need to be further elucidated. Up till now, fewer comprehensive comparisons of the mitochondrial proteomes from cotton CMS line and maintainer line have been reported. In this study, we performed a novel comparison of mitochondrial protein profiles in two CMS lines and their common maintainer line. Based on our results, we found a potential protein related to oxidative stress led to the anthers abortion. These results accumulate data to interpret the molecular mechanisms of CMS in cotton.

Integrated Methylome and Transcriptome Analysis between the CMS-D2 Line ZBA and Its Maintainer Line ZB in Upland Cotton

DNA methylation is an important epigenetic modification involved in multiple biological processes. Altered methylation patterns have been reported to be associated with male sterility in some plants, but their role in cotton cytoplasmic male sterility (CMS) remains unclear. Here, integrated methylome and transcriptome analyses were conducted between the CMS-D2 line ZBA and its near-isogenic maintainer line ZB in upland cotton. More methylated cytosine sites (mCs) and higher methylation levels (MLs) were found among the three sequence contexts in ZB compared to ZBA. A total of 4568 differentially methylated regions (DMRs) and 2096 differentially methylated genes (DMGs) were identified. Among the differentially expressed genes (DEGs) associated with DMRs (DMEGs), 396 genes were upregulated and 281 genes were downregulated. A bioinformatics analysis of these DMEGs showed that hyper-DEGs were significantly enriched in the “oxidative phosphorylation” pathway. Further qRT-PCR validation indicated that these hypermethylated genes (encoding the subunits of mitochondrial electron transport chain (ETC) complexes I and V) were all significantly upregulated in ZB. Our biochemical data revealed a higher extent of H₂O₂ production but a lower level of adenosine triphosphate (ATP) synthesis in CMS-D2 line ZBA. On the basis of the above results, we propose that disrupted DNA methylation in ZBA may disrupt the homeostasis of reactive oxygen species (ROS) production and ATP synthesis in mitochondria, triggering a burst of ROS that is transferred to the nucleus to initiate programmed cell death (PCD) prematurely, ultimately leading to microspore abortion. This study illustrates the important role of DNA methylation in cotton CMS.

Intergenomic gene transfer in diploid and allopolyploid Gossypium

BACKGROUND

Intergenomic gene transfer (IGT) between nuclear and organellar genomes is a common phenomenon during plant evolution. Gossypium is a useful model to evaluate the genomic consequences of IGT for both diploid and polyploid species. Here, we explore IGT among nuclear, mitochondrial, and plastid genomes of four cotton species, including two allopolyploids and their model diploid progenitors (genome donors, G. arboreum: A₂ and G. raimondii: D₅).

RESULTS

Extensive IGT events exist for both diploid and allotetraploid cotton (Gossypium) species, with the nuclear genome being the predominant recipient of transferred DNA followed by the mitochondrial genome. The nuclear genome has integrated 100 times more foreign sequences than the mitochondrial genome has in total length. In the nucleus, the integrated length of chloroplast DNA (cpDNA) was between 1.87 times (in diploids) to nearly four times (in allopolyploids) greater than that of mitochondrial DNA (mtDNA). In the mitochondrion, the length of nuclear DNA (nuDNA) was typically three times than that of cpDNA. Gossypium mitochondrial genomes integrated three nuclear retrotransposons and eight chloroplast tRNA genes, and incorporated chloroplast DNA prior to divergence between the diploids and allopolyploid formation. For mitochondrial chloroplast-tRNA genes, there were 2-6 bp conserved microhomologies flanking their insertion sites across distantly related genera, which increased to 10 bp microhomologies for the four cotton species studied. For organellar DNA sequences, there are source hotspots, e.g., the atp6-trnW intergenic region in the mitochondrion and the inverted repeat region in the chloroplast. Organellar DNAs in the nucleus were rarely expressed, and at low levels. Surprisingly, there was asymmetry in the survivorship of ancestral insertions following allopolyploidy, with most numts (nuclear mitochondrial insertions) decaying or being lost whereas most nupts (nuclear plastidial insertions) were retained.

CONCLUSIONS

This study characterized and compared intracellular transfer among nuclear and organellar genomes within two cultivated allopolyploids and their ancestral diploid cotton species. A striking asymmetry in the fate of IGTs in allopolyploid cotton was discovered, with numts being preferentially lost relative to nupts. Our results connect intergenomic gene transfer with allotetraploidy and provide new insight into intracellular genome evolution.

Using Transcriptome Analysis to Screen for Key Genes and Pathways Related to Cytoplasmic Male Sterility in Cotton (Gossypium hirsutum L.)

Cotton (Gossypium hirsutum L.) is one of the most important cash crops worldwide. Cytoplasmic male sterility (CMS) is an excellent breeding system for exploitation of heterosis, which has great potential to increase crop yields. To understand the molecular mechanism of CMS in cotton, we compared transcriptome, cytomorphological, physiological and bioinformatics data between the CMS line C2P5A and its maintainer line C2P5B. By using high-throughput sequencing technology, 178,166 transcripts were assembled and 2013 differentially expression genes (DEGs) were identified at three different stages of C2P5A anther development. In this study, we identified DEGs associated with reactive oxygen species (ROS), peroxisomes, aldehyde dehydrogenases (ALDH), cytochrome oxidase subunit VI, and cytochrome P450, and DEGs associated with tapetum development, Jojoba acyl-CoA reductase-related male sterility protein, basic helix-loop-helix (bHLH) and MYB transcription factors. The abnormal expression of one of these genes may be responsible for the CMS C2P5A line. In gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, DEGs were mainly related to carbohydrate metabolism, amino acid metabolism, transport and catabolism, and signal transduction. Carbohydrate metabolism provides energy for anther development, starch and sucrose metabolism, fatty acid biosynthesis and metabolism and ascorbate and aldarate metabolism. These results showed that numerous genes and multiple complex metabolic pathways regulate cotton anther development. Weighted correlation network analysis (WGCNA) indicated that three modules, ‘turquoise,’ ‘blue,’ and ‘green,’ were specific for the CMS C2P5A line. The ‘turquoise’ and ‘blue’ modules were mainly related to carbohydrate metabolism, amino acid metabolism, energy metabolism, peroxisomes, pyruvate metabolism as well as fatty acid degradation. The ‘green’ module was mainly related to energy metabolism, carbon metabolism, translation, and lipid metabolism. RNA-sequencing and WGCNA polymerization modules were screened for key genes and pathways related to CMS in cotton. This study presents a new perspective for further research into the metabolic pathways of pollen abortion in the CMS C2P5A line and also provides a theoretical basis for its breeding and production.