Comprehensive analysis of Ogura cytoplasmic male sterility-related genes in turnip (Brassica rapa ssp. rapifera) using RNA sequencing analysis and bioinformatics

Transcription Factors and Their Regulatory Roles in the Male Gametophyte Development of Flowering Plants

Male gametophyte development in plants relies on the functions of numerous genes, whose expression is regulated by transcription factors (TFs), non-coding RNAs, hormones, and diverse environmental stresses. Several excellent reviews are available that address the genes and enzymes associated with male gametophyte development, especially pollen wall formation. Growing evidence from genetic studies, transcriptome analysis, and gene-by-gene studies suggests that TFs coordinate with epigenetic machinery to regulate the expression of these genes and enzymes for the sequential male gametophyte development. However, very little summarization has been performed to comprehensively review their intricate regulatory roles and discuss their downstream targets and upstream regulators in this unique process. In the present review, we highlight the research progress on the regulatory roles of TF families in the male gametophyte development of flowering plants. The transcriptional regulation, epigenetic control, and other regulators of TFs involved in male gametophyte development are also addressed.

Global transcriptome analysis reveals potential genes associated with genic male sterility of rapeseed (Brassica napus L.)

Rapeseed is the third leading source of edible oil in the world. Genic male sterility (GMS) lines provide crucial material for harnessing heterosis for rapeseed. GMS lines have been widely used successfully for rapeseed hybrid production. The physiological and molecular mechanism of pollen development in GMS lines of rapeseed (Brassica napus L.) need to be determined for the creation of hybrids and cultivation of new varieties. However, limited studies have focused on systematically mining genes that regulate the pollen development of GMS lines in B. napus. In the present study, to determine the stage at which pollen development begins to show abnormality in the GMS lines, we performed semi-thin section analysis of the anthers with five pollen development stages. The results indicated that the abnormal pollen development in DGMS lines might start at the meiotic stage, and abnormal pollen development in RGMS lines probably occurred before the tetrad stage. To investigate the critical genes and pathways involved in pollen development in GMS lines, we constructed and sequenced 24 transcriptome libraries for the flower buds from the fertile and sterile lines of two recessive GMS (RGMS) lines (6251AB and 6284AB) and two dominant GMS (DGMS) lines (4001AB and 4006AB). A total of 23,554 redundant DEGs with over two-fold change between sterile and fertile lines were obtained. A total of 346 DEGs were specifically related to DGMS, while 1,553 DEGs were specifically related to RGMS. A total of 1,545 DEGs were shared between DGMS and RGMS. And 253 transcription factors were found to be differentially expressed between the sterile and fertile lines of GMS. In addition, 6,099 DEGs possibly related to anther, pollen, and microspore development processes were identified. Many of these genes have been reported to be involved in anther and microspore developmental processes. Several DEGs were speculated to be key genes involved in the regulation of fertility. Three differentially expressed genes were randomly selected and their expression levels were verified by quantitative PCR (qRT-PCR). The results of qRT-PCR largely agreed with the transcriptome sequencing results. Our findings provide a global view of genes that are potentially involved in GMS occurrence. The expression profiles and function analysis of these DEGs were provided to expand our understanding of the complex molecular mechanism in pollen and sterility development in B. napus.

Mechanism and Utilization of Ogura Cytoplasmic Male Sterility in Cruciferae Crops

Hybrid production using lines with cytoplasmic male sterility (CMS) has become an important way to utilize heterosis in vegetables. Ogura CMS, with the advantages of complete pollen abortion, ease of transfer and a progeny sterility rate reaching 100%, is widely used in cruciferous crop breeding. The mapping, cloning, mechanism and application of Ogura CMS and fertility restorer genes in Brassica napus, Brassica rapa, Brassica oleracea and other cruciferous crops are reviewed herein, and the existing problems and future research directions in the application of Ogura CMS are discussed.

Transcriptomics Integrated with Metabolomics Unveil Carotenoids Accumulation and Correlated Gene Regulation in White and Yellow-Fleshed Turnip (Brassica rapa ssp. rapa)

Turnip (Brassica rapa ssp. rapa) is considered to be a highly nutritious and health-promoting vegetable crop, whose flesh color can be divided into yellow and white. It is widely accepted that yellow-fleshed turnips have higher nutritional value. However, reports about flesh color formation is lacking. Here, the white-fleshed inbred line, W21, and yellow-fleshed inbred line, W25, were profiled from the swollen root of the turnip at three developmental periods to elucidate the yellow color formation. Transcriptomics integrated with metabolomics analysis showed that the PSY gene was the key gene affecting the carotenoids formation in W25. The coding sequence of BrrPSY-W25 was 1278 bp and that of BrrPSY-W21 was 1275 bp, and BrrPSY was more highly expressed in swollen roots in W25 than in W21. Transient transgenic tobacco leaf over-expressing BrrPSY-W and BrrPSY-Y showed higher transcript levels and carotenoids contents. Results revealed that yellow turnip formation is due to high expression of the PSY gene rather than mutations in the PSY gene, indicating that a post-transcriptional regulatory mechanism may affect carotenoids formation. Results obtained in this study will be helpful for explaining the carotenoids accumulation of turnips.

Transcript Profiling Analysis and ncRNAs' Identification of Male-Sterile Systems of Brassica campestris Reveal New Insights Into the Mechanism Underlying Anther and Pollen Development

Male-sterile mutants are useful materials to study the anther and pollen development. Here, whole transcriptome sequencing was performed for inflorescences in three sterile lines of Chinese cabbage (Brassica campestris L. ssp. chinensis Makino, syn. B. rapa ssp. chinensis), the genic male-sterile line (A line), the Polima cytoplasmic male-sterile (CMS) line (P line), and the Ogura CMS line (O line) along with their maintainer line (B line). In total, 7,136 differentially expressed genes (DEGs), 361 differentially expressed long non-coding RNAs (lncRNAs) (DELs), 56 differentially expressed microRNAs (miRNAs) (DEMs) were selected out. Specific regulatory networks related to anther cell differentiation, meiosis cytokinesis, pollen wall formation, and tapetum development were constructed based on the abortion characteristics of male-sterile lines. Candidate genes and lncRNAs related to cell differentiation were identified in sporocyteless P line, sixteen of which were common to the DEGs in Arabidopsis spl/nzz mutant. Genes and lncRNAs concerning cell plate formation were selected in A line that is defected in meiosis cytokinesis. Also, the orthologs of pollen wall formation and tapetum development genes in Arabidopsis showed distinct expression patterns in the three different sterile lines. Among 361 DELs, 35 were predicted to interact with miRNAs, including 28 targets, 47 endogenous target mimics, and five precursors for miRNAs. Two lncRNAs were further proved to be functional precursors for bra-miR156 and bra-miR5718, respectively. Overexpression of bra-miR5718HG in B. campestris slowed down the growth of pollen tubes, caused shorter pollen tubes, and ultimately affected the seed set. Our study provides new insights into molecular regulation especially the ncRNA interaction during pollen development in Brassica crops.

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.

Genome-wide comparative transcriptome analysis of the A4-CMS line ICPA 2043 and its maintainer ICPB 2043 during the floral bud development of pigeonpea

Cytoplasmic male sterility (CMS) offers a unique system to understand cytoplasmic nuclear crosstalk, and is also employed for exploitation of hybrid vigor in various crops. Pigeonpea A4-CMS, a predominant source of male sterility, is being used for efficient hybrid seed production. The molecular mechanisms of CMS trait remain poorly studied in pigeonpea. We performed genome-wide transcriptome profiling of A4-CMS line ICPA 2043 and its isogenic maintainer ICPB 2043 at two different stages of floral bud development (stage S1 and stage S2). Consistent with the evidences from some other crops, we also observed significant difference in the expression levels of genes in the later stage, i.e., stage S2. Differential expression was observed for 143 and 55 genes within the two stages of ICPA 2043 and ICPB 2043, respectively. We obtained only 10 differentially expressed genes (DEGs) between the stage S1 of the two genotypes, whereas expression change was significant for 582 genes in the case of stage S2. The qRT-PCR assay of randomly selected six genes supported the differential expression of genes between ICPA 2043 and ICPB 2043. Further, GO and KEGG pathway mapping suggested a possible compromise in key bioprocesses during flower and pollen development. Besides providing novel insights into the functional genomics of CMS trait, our results were in strong agreement with the gene expression atlas of pigeonpea that implicated various candidate genes like sucrose-proton symporter 2 and an uncharacterized protein along with pectate lyase, pectinesterase inhibitors, L-ascorbate oxidase homolog, ATPase, β-galactosidase, polygalacturonase, and aldose 1-epimerase for pollen development of pigeonpea. The dataset presented here provides a rich genomic resource to improve understanding of CMS trait and its deployment in heterosis breeding in pigeonpea.

Transcriptomic and proteomic analyses of a new cytoplasmic male sterile line with a wild Gossypium bickii genetic background

Background

Cotton is an important fiber crop but has serious heterosis effects, and cytoplasmic male sterility (CMS) is the major cause of heterosis in plants. However, to the best of our knowledge,

no studies have investigated CMS Yamian A in cotton with the genetic background of Australian wild Gossypium bickii. Conjoint transcriptomic and proteomic analysis was first performed between Yamian A and its maintainer Yamian B.

Results

We detected 550 differentially expressed transcript-derived fragments (TDFs) and at least 1013 proteins in anthers at various developmental stages. Forty-two TDFs and 11 differentially expressed proteins (DEPs) were annotated by analysis in the genomic databases of G. austral, G. arboreum and G. hirsutum. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed to better understand the functions of these TDFs and DEPs. Transcriptomic and proteomic results showed that UDP-glucuronosyl/UDP-glucosyltransferase, 60S ribosomal protein L13a-4-like, and glutathione S-transferase were upregulated; while heat shock protein Hsp20, ATPase, F0 complex, and subunit D were downregulated at the microspore abortion stage of Yamian A. In addition, several TDFs from the transcriptome and several DEPs from the proteome were detected and confirmed by quantitative real-time PCR as being expressed in the buds of seven different periods of development. We established the databases of differentially expressed genes and proteins between Yamian A and its maintainer Yamian B in the anthers at various developmental stages and constructed an interaction network based on the databases for a comprehensive understanding of the mechanism underlying CMS with a wild cotton genetic background.

Conclusion

We first analyzed the molecular mechanism of CMS Yamian A from the perspective of omics, thereby providing an experimental basis and theoretical foundation for future research attempting to analyze the abortion mechanism of new CMS with a wild Gossypium bickii background and to realize three-line matching.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07261-y.

Identification of microRNAs and their targets in inflorescences of an Ogura-type cytoplasmic male-sterile line and its maintainer fertile line of turnip (Brassica rapa ssp. rapifera) via high-throughput sequencing and degradome analysis

Cytoplasmic male sterility (CMS) is a widely used trait in angiosperms caused by perturbations in nucleus-mitochondrion interactions that suppress the production of functional pollen. MicroRNAs (miRNAs) are small non-coding RNAs that act as regulatory molecules of transcriptional or post-transcriptional gene silencing in plants. The discovery of miRNAs and their possible implications in CMS induction provides clues for the intricacies and complexity of this phenomenon. Previously, we characterized an Ogura-CMS line of turnip (Brassica rapa ssp. rapifera) that displays distinct impaired anther development with defective microspore production and premature tapetum degeneration. In the present study, high-throughput sequencing was employed for a genome-wide investigation of miRNAs. Six small RNA libraries of inflorescences collected from the Ogura-CMS line and its maintainer fertile (MF) line of turnip were constructed. A total of 120 pre-miRNAs corresponding to 89 mature miRNAs were identified, including 87 conversed miRNAs and 33 novel miRNAs. Among these miRNAs, the expression of 10 differentially expressed mature miRNAs originating from 12 pre-miRNAs was shown to have changed by more than two-fold between inflorescences of the Ogura-CMS line and inflorescences of the MF line, including 8 down- and 2 up-regulated miRNAs. The expression profiles of the differentially expressed miRNAs were confirmed by stem-loop quantitative real-time PCR. In addition, to identify the targets of the identified miRNAs, a degradome analysis was performed. A total of 22 targets of 25 miRNAs and 17 targets of 28 miRNAs were identified as being involved in the reproductive development for Ogura-CMS and MF lines of turnip, respectively. Negative correlations of expression patterns between partial miRNAs and their targets were detected. Some of these identified targets, such as squamosa promoter-binding-like transcription factor family proteins, auxin response factors and pentatricopeptide repeat-containing proteins, were previously reported to be involved in reproductive development in plants. Taken together, our results can help improve the understanding of miRNA-mediated regulatory pathways that might be involved in CMS occurrence in turnip.

Comparative proteomic analysis of mitochondrial proteins from maize CMS-C sterile, maintainer and restorer anthers

The maize C system of cytoplasmic male sterility (CMS) and its fertility restoration gene Rf4 have been widely used for maize hybrid production; however, the underlying mechanism is still uncertain. The sterility factor functions in mitochondria, where it interacts directly or indirectly with the restorer. Mitoproteomics can capture all participants involved in CMS and restoration at the organelle level. In the present study, we identified and quantified anther mitochondrial proteins from CMS, maintainer and restorer lines. We obtained 14,528 unique peptides belonging to 3,369 proteins. Comparative analysis of 1840 high-confidence proteins revealed 68 were differentially accumulated proteins likely involved in CMS or its restoration within mitochondria. These proteins were mainly associated with fatty acid metabolism, amino acid metabolism and protein-processing pathways. These results suggest that an energy deficiency caused by the sterility factor hinders other proteins or protein complexes required for pollen development through nuclear-mitochondrial interaction. The restorer factor may boost the energy generation by activating alternative metabolic pathways and by improving the post-translation processing efficiency of proteins in energy-producing complexes to restore pollen fertility. Our findings may aid detailed molecular analysis and contribute to a better understanding of maize CMS-C restoration and sterility.

Genome-wide analysis of proline-rich extension-like receptor protein kinase (PERK) in Brassica rapa and its association with the pollen development

BACKGROUND

Proline-rich extension-like receptor protein kinases (PERKs) are an important class of receptor kinases located in the plasma membrane, most of which play a vital role in pollen development.

RESULTS

Our study identified 25 putative PERK genes from the whole Brassica rapa genome (AA). Phylogenetic analysis of PERK protein sequences from 16 Brassicaceae species divided them into four subfamilies. The biophysical properties of the BrPERKs were investigated. Gene duplication and synteny analyses and the calculation of Ka/Ks values suggested that all 80 orthologous/paralogous gene pairs between B. rapa and A. thaliana, B. nigra and B. oleracea have experienced strong purifying selection. RNA-Seq data and qRT-PCR analyses showed that several BrPERK genes were expressed in different tissues, while some BrPERKs exhibited high expression levels only in buds. Furthermore, comparative transcriptome analyses from six male-sterile lines of B. rapa indicated that 7 BrPERK genes were downregulated in all six male-sterile lines. Meanwhile, the interaction networks of the BrPERK genes were constructed and 13 PERK coexpressed genes were identified, most of which were downregulated in the male sterile buds.

CONCLUSION

Combined with interaction networks, coexpression and qRT-PCR analyses, these results demonstrated that two BrPERK genes, Bra001723.1 and Bra037558.1 (the orthologs of AtPERK6 (AT3G18810)), were downregulated beginning in the meiosis II period of male sterile lines and involved in anther development. Overall, this comprehensive analysis of some BrPERK genes elucidated their roles in male sterility.

Comparative transcriptome analysis in Chinese cabbage (Brassica rapa ssp. pekinesis) for DEGs of Ogura-, Polima-CMS and their shared maintainer

Cytoplasmic male sterility (CMS) is maternally inherited trait, which hinders the ability to produce viable pollen in plants. It serves as a useful tool for hybrid seed production via exploiting heterosis in crops. The molecular mechanism of CMS and fertility restoration has been investigated in different crops. However, limited number of reports is available on comparison of Ogura- and Polima-CMS with their shared maintainer in Chinese cabbage. We performed transcript profiling of sterile Ogura CMS (Tyms), Polima CMS (22m2) and their shared maintainer line (231-330) with an aim to identify genes associated with male sterility. In this work, we identified 912, 7199 and 6381 DEGs (Differentially Expressed Genes) in 22m2 Vs Tyms, 231-330 VS 22m2 and 231-330 Vs Tyms, respectively. The GO (Gene Ontology) annotation and KEGG pathway analysis suggested that most of the DEGs were involved in pollen development, carbon metabolism, lipase activity, lipid binding, penta-tricopeptide repeat (PPR), citrate cycle and oxidative phosphorylation, which were down-regulated in both CMS lines. This result will provide an important resource for further understanding of functional pollen development, the CMS mechanism and to improve molecular breeding in Chinese cabbage.