Identification of miRNAs and their targets by high-throughput sequencing and degradome analysis in cytoplasmic male-sterile line NJCMS1A and its maintainer NJCMS1B of soybean

The roles of non-coding RNAs in male reproductive development and abiotic stress responses during this unique process in flowering plants

Successful male reproductive development is the guarantee for sexual reproduction of flowering plants. Male reproductive development is a complicated and multi-stage process that integrates physiological processes and adaptation and tolerance to a myriad of environmental stresses. This well-coordinated process is governed by genetic and epigenetic machineries. Non-coding RNAs (ncRNAs) play pleiotropic roles in the plant growth and development. The identification, characterization and functional analysis of ncRNAs and their target genes have opened a new avenue for comprehensively revealing the regulatory network of male reproductive development and its response to environmental stresses in plants. This review briefly addresses the types, origin, biogenesis and mechanisms of ncRNAs in plants, highlights important updates on the roles of ncRNAs in regulating male reproductive development and emphasizes the contribution of ncRNAs, especially miRNAs and lncRNAs, in responses to abiotic stresses during this unique process in flowering plants.

Heat-responsive microRNAs participate in regulating the pollen fertility stability of CMS-D2 restorer line under high-temperature stress

Anther development and pollen fertility of cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) restorer lines are susceptible to continuous high-temperature (HT) stress in summer, which seriously hinders the large-scale application of "three-line" hybrids in production. Here, integrated small RNA, transcriptome, degradome, and hormone profiling was performed to explore the roles of microRNAs (miRNAs) in regulating fertility stability in mature pollens of isonuclear alloplasmic near-isogenic restorer lines NH and SH under HT stress at two environments. A total of 211 known and 248 novel miRNAs were identified, of which 159 were differentially expressed miRNAs (DEMs). Additionally, 45 DEMs in 39 miRNA clusters (PmCs) were also identified, and most highly expressed miRNAs were significantly induced in SH under extreme HT, especially four MIR482 and six MIR6300 family miRNAs. PmC28 was located in the fine-mapped interval of the Rf1 gene and contained two DEMs, gra-miR482_L-2R + 2 and gma-miR2118a-3p_R + 1_1ss18TG. Transcriptome sequencing identified 6281 differentially expressed genes, of which heat shock protein (HSP)-related genes, such as HSP70, HSP22, HSP18.5-C, HSP18.2 and HSP17.3-B, presented significantly reduced expression levels in SH under HT stress. Through integrating multi-omics data, we constructed a comprehensive molecular network of miRNA-mRNA-gene-KEGG containing 35 pairs of miRNA/target genes involved in regulating the pollen development in response to HT, among which the mtr-miR167a_R + 1, tcc-miR167c and ghr-miR390a, tcc-miR396c_L-1 and ghr-MIR169b-p3_1ss6AG regulated the pollen fertility by influencing ARF8 responsible for the auxin signal transduction, ascorbate and aldarate metabolism, and the sugar and lipid metabolism and transport pathways, respectively. Further combination with hormone analysis revealed that HT-induced jasmonic acid signaling could activate the expression of downstream auxin synthesis-related genes and cause excessive auxin accumulation, followed by a cascade of auxin signal transduction, ultimately resulting in pollen abortion. The results provide a new understanding of how heat-responsive miRNAs regulate the stability of fertility restoration for CMS-D2 cotton under heat stress.

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.

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.

A Multi-Level Iterative Bi-Clustering Method for Discovering miRNA Co-regulation Network of Abiotic Stress Tolerance in Soybeans

Although growing evidence shows that microRNA (miRNA) regulates plant growth and development, miRNA regulatory networks in plants are not well understood. Current experimental studies cannot characterize miRNA regulatory networks on a large scale. This information gap provides an excellent opportunity to employ computational methods for global analysis and generate valuable models and hypotheses. To address this opportunity, we collected miRNA-target interactions (MTIs) and used MTIs from Arabidopsis thaliana and Medicago truncatula to predict homologous MTIs in soybeans, resulting in 80,235 soybean MTIs in total. A multi-level iterative bi-clustering method was developed to identify 483 soybean miRNA-target regulatory modules (MTRMs). Furthermore, we collected soybean miRNA expression data and corresponding gene expression data in response to abiotic stresses. By clustering these data, 37 MTRMs related to abiotic stresses were identified, including stress-specific MTRMs and shared MTRMs. These MTRMs have gene ontology (GO) enrichment in resistance response, iron transport, positive growth regulation, etc. Our study predicts soybean MTRMs and miRNA-GO networks under different stresses, and provides miRNA targeting hypotheses for experimental analyses. The method can be applied to other biological processes and other plants to elucidate miRNA co-regulation mechanisms.

Genome-wide identification of PME genes, evolution and expression analyses in soybean (Glycine max L.)

BACKGROUND

Pectin methylesterase (PME) is one of pectin-modifying enzyme that affects the pectin homeostasis in cell wall and regulates plant growth and diverse biological processes. The PME genes have been well explored and characterized in different plants. Nevertheless, systematic research on the soybean (Glycine max L.) PME genes remain lacking.

RESULTS

We identified 127 Glycine max PME genes (GmPME) from the soybean Wm82.a2.v1 genome, which unevenly distributed on 20 soybean chromosomes. Phylogenetic analysis classified the GmPME genes into four clades (Group I, Group II, Group III and Group IV). GmPME gene members in the same clades displayed similar gene structures and motif patterns. The gene family expansion analysis demonstrated that segmental duplication was the major driving force to acquire novel GmPME genes compared to the tandem duplication events. Further synteny and evolution analyses showed that the GmPME gene family experienced strong purifying selective pressures during evolution. The cis-element analyses together with the expression patterns of the GmPME genes in various tissues suggested that the GmPME genes broadly participate in distinct biological processes and regulate soybean developments. Importantly, based on the transcriptome data and quantitative RT-PCR validations, we examined the potential roles of the GmPME genes in regulating soybean flower bud development and seed germination.

CONCLUSION

In conclusion, we provided a comprehensive characterization of the PME genes in soybean, and our work laid a foundation for the functional study of GmPME genes in the future.

A small RNA of miR2119b from soybean CMS line acts as a negative regulator of male fertility in transgenic Arabidopsis

The miR2119 is involved in the growth, development and abiotic stress response of some legumes, including Medicago truncatula, Phaseolus vulgaris and soybean (Glycine max (L.) Merr.). Our previous small RNA sequencing analysis showed that miR2119b was up-regulated in the flower buds of soybean cytoplasmic male sterile (CMS) line compared with its maintainer line, but the role and mechanism of miR2119b in the regulation of soybean male fertility are still unclear. In this study, the gma-miR2119b and its target gene alcohol dehydrogenase 1.3b (ADH1.3b) were characterized and found to be highly expressed in the flowers of soybean CMS line and its maintainer. Transgenic Arabidopsis plants overexpressing gma-miR2119b exhibit male fertility abnormalities, including pollen fertility and germination rate decreased. Enzyme activity detection found the ADH and catalase (CAT) enzyme activities in inflorescence of gma-miR2119b overexpressed plants were lower than those of wild-type. Bioinformatics and gene expression analysis showed that gma-miR2119b/GmADH1.3b module was responsive to high temperature (HT) stress during flowering. After HT stress, the gma-miR2119b overexpressed plants showed male sterility, including shorter filament, sterile pollen, indehiscent anther and non seed. Moreover, some key genes involved in HT response and reactive oxygen species (ROS) signal regulation pathway, including heat shock protein70, galactinol synthase 1 and CAT, showed down-regulated expression in transgenic plants under HT stress, suggesting that gma-miR2119b regulates male fertility via HT-ROS signaling pathway under HT stress. It was speculated that the gma-miR2119b acted as a negative regulator of male fertility in plants by regulating ADH1, HT-induced and ROS scavenging genes expression.

Non-Coding RNAs in Legumes: Their Emerging Roles in Regulating Biotic/Abiotic Stress Responses and Plant Growth and Development

Noncoding RNAs, including microRNAs (miRNAs), small interference RNAs (siRNAs), circular RNA (circRNA), and long noncoding RNAs (lncRNAs), control gene expression at the transcription, post-transcription, and translation levels. Apart from protein-coding genes, accumulating evidence supports ncRNAs playing a critical role in shaping plant growth and development and biotic and abiotic stress responses in various species, including legume crops. Noncoding RNAs (ncRNAs) interact with DNA, RNA, and proteins, modulating their target genes. However, the regulatory mechanisms controlling these cellular processes are not well understood. Here, we discuss the features of various ncRNAs, including their emerging role in contributing to biotic/abiotic stress response and plant growth and development, in addition to the molecular mechanisms involved, focusing on legume crops. Unravelling the underlying molecular mechanisms and functional implications of ncRNAs will enhance our understanding of the coordinated regulation of plant defences against various biotic and abiotic stresses and for key growth and development processes to better design various legume crops for global food security.

MicroRNAs Involved in Regulatory Cytoplasmic Male Sterility by Analysis RNA-seq and Small RNA-seq in Soybean

Cytoplasmic male sterility (CMS) is an important plant characteristic for exploiting heterosis to enhance crop traits during breeding. However, the CMS regulatory network remains unclear in plants, even though researchers have attempted to isolate genes associated with CMS. In this study, we performed high-throughput sequencing and degradome analyses to identify microRNAs (miRNAs) and their targets in a soybean CMS line (JLCMS9A) and its maintainer line (JLCMS9B). Additionally, the differentially expressed genes during reproductive development were identified using RNA-seq data. A total of 280 miRNAs matched soybean miRNA sequences in miRBase, including mature miRNAs and pre-miRNAs. Of the 280 miRNAs, 30, 23, and 21 belonged to the miR166, miR156, and miR171 families, respectively. Moreover, 410 novel low-abundant miRNAs were identified in the JLCMS9A and JLCMS9B flower buds. Furthermore, 303 and 462 target genes unique to JLCMS9A and JLCMS9B, respectively, as well as 782 common targets were predicted based on the degradome analysis. Target genes differentially expressed between the CMS line and the maintainer line were revealed by an RNA-seq analysis. Moreover, all target genes were annotated with diverse functions related to biological processes, cellular components, and molecular functions, including transcriptional regulation, the nucleus, meristem maintenance, meristem initiation, cell differentiation, auxin-activated signaling, plant ovule development, and anther development. Finally, a network was built based on the interactions. Analyses of the miRNA, degradome, and transcriptome datasets generated in this study provided a comprehensive overview of the reproductive development of a CMS soybean line. The data presented herein represent useful information for soybean hybrid breeding. Furthermore, the study results indicate that miRNAs might contribute to the soybean CMS regulatory network by modulating the expression of CMS-related genes. These findings lay the foundation for future studies on the molecular mechanisms underlying soybean CMS.

Title: Hypermethylation of miRNA Genes During Nodule Development

DNA methylation has recently emerged as a powerful regulatory mechanism controlling the expression of key regulators of various developmental processes, including nodulation. However, the functional role of DNA methylation in regulating the expression of microRNA (miRNA) genes during the formation and development of nitrogen-fixing nodules remains largely unknown. In this study, we profiled DNA methylation patterns of miRNA genes during nodule formation, development, and early senescence stages in soybean (Glycine max) through the analysis of methylC-seq data. Absolute DNA methylation levels in the CG, CHH, and CHH sequence contexts over the promoter and primary transcript regions of miRNA genes were significantly higher in the nodules compared with the corresponding root tissues at these three distinct nodule developmental stages. We identified a total of 82 differentially methylated miRNAs in the nodules compared with roots. Differential DNA methylation of these 82 miRNAs was detected only in the promoter (69), primary transcript region (3), and both in the promoter and primary transcript regions (10). The large majority of these differentially methylated miRNAs were hypermethylated in nodules compared with the corresponding root tissues and were found mainly in the CHH context and showed stage-specific methylation patterns. Differentially methylated regions in the promoters of 25 miRNAs overlapped with transposable elements, a finding that may explain the vulnerability of miRNAs to DNA methylation changes during nodule development. Gene expression analysis of a set of promoter-differentially methylated miRNAs pointed to a negative association between DNA methylation and miRNA expression. Gene Ontology and pathways analyses indicate that changes in DNA methylation of miRNA genes are reprogrammed and contribute to nodule development through indirect regulation of genes involved in cellular processes and pathways with well-established roles in nodulation.

Heat-Responsive miRNAs Participate in the Regulation of Male Fertility Stability in Soybean CMS-Based F1 under High Temperature Stress

MicroRNAs (miRNAs), a class of noncoding small RNAs (sRNAs), are widely involved in the response to high temperature (HT) stress at both the seedling and flowering stages. To dissect the roles of miRNAs in regulating male fertility in soybean cytoplasmic male sterility (CMS)-based F1 under HT, sRNA sequencing was performed using flower buds from HT-tolerant and HT-sensitive CMS-based F1 combinations (NF1 and YF1, respectively). A total of 554 known miRNAs, 59 new members of known miRNAs, 712 novel miRNAs, and 1145 target genes of 580 differentially expressed miRNAs (DEMs) were identified under normal temperature and HT conditions. Further integrated analysis of sRNA and transcriptome sequencing found that 21 DEMs and 15 differentially expressed target genes, such as gma-miR397a/Laccase 2, gma-miR399a/Inorganic phosphate transporter 1-4, and gma-miR4413a/PPR proteins, mitochondrial-like, were negatively regulated under HT stress. Furthermore, all members of the gma-miR156 family were suppressed by HT stress in both NF1 and YF1, but were highly expressed in YF1 under HT condition. The negative correlation between gma-miR156b and its target gene squamosa promoter-binding protein-like 2b was confirmed by expression analysis, and overexpression of gma-miR156b in Arabidopsis led to male sterility under HT stress. With these results, we proposed that miRNAs play an important role in the regulation of male fertility stability in soybean CMS-based F1 under HT stress.

Identification of microRNAs and their gene targets in cytoplasmic male sterile and fertile maintainer lines of pigeonpea

Comparative analysis of genome-wide miRNAs and their gene targets between cytoplasmic male sterile (CMS) and fertile lines of pigeonpea suggests a possible role of miRNA-regulated pathways in reproductive development. Exploitation of hybrid vigor using CMS technology has delivered nearly 50% yield gain in pigeonpea. Among various sterility-inducing cytoplasms (A₁-A₉) reported so far in pigeonpea, A₂ and A₄ are the two major sources that facilitate hybrid seed production. Recent evidence suggests involvement of micro RNA in vast array of biological processes including plant reproductive development. In pigeonpea, information about the miRNAs is insufficient. In view of this, we sequenced six small RNA libraries of CMS line UPAS 120A and isogenic fertile line UPAS 120B using Illumina technology. Results revealed 316 miRNAs including 248 known and 68 novel types. A total of 637 gene targets were predicted for known miRNAs, while 324 genes were associated with novel miRNAs. Degradome analysis revealed 77 gene targets of predicted miRNAs, which included a variety of transcription factors playing key roles in plant reproduction such as F-box family proteins, apetala 2, auxin response factors, ethylene-responsive factors, homeodomain-leucine zipper proteins etc. Differential expression of both known and novel miRNAs implied roles for both conserved as well as species-specific players. We also obtained several miRNA families such as miR156, miR159, miR167 that are known to influence crucial aspects of plant fertility. Gene ontology and pathway level analyses of the target genes showed their possible implications for crucial events during male reproductive development such as tapetal degeneration, pollen wall formation, retrograde signaling etc. To the best of our knowledge, present study is first to combine deep sequencing of small RNA and degradome for elucidating the role of miRNAs in flower and male reproductive development in pigeonpea.

Noncoding RNA: An Insight into Chloroplast and Mitochondrial Gene Expressions

Regulation of gene expression in any biological system is a complex process with many checkpoints at the transcriptional, post-transcriptional and translational levels. The control mechanism is mediated by various protein factors, secondary metabolites and a newly included regulatory member, i.e., noncoding RNAs (ncRNAs). It is known that ncRNAs modulate the mRNA or protein profiles of the cell depending on the degree of complementary and context of the microenvironment. In plants, ncRNAs are essential for growth and development in normal conditions by controlling various gene expressions and have emerged as a key player to guard plants during adverse conditions. In order to have smooth functioning of the plants under any environmental pressure, two very important DNA-harboring semi-autonomous organelles, namely, chloroplasts and mitochondria, are considered as main players. These organelles conduct the most crucial metabolic pathways that are required to maintain cell homeostasis. Thus, it is imperative to explore and envisage the molecular machineries responsible for gene regulation within the organelles and their coordination with nuclear transcripts. Therefore, the present review mainly focuses on ncRNAs origination and their gene regulation in chloroplasts and plant mitochondria.

Identification of Differentially Methylated miRNA Genes During Compatible and Incompatible Interactions Between Soybean and Soybean Cyst Nematode

DNA methylation is a widespread epigenetic mark that affects gene expression and transposon mobility during plant development and stress responses. However, the role of DNA methylation in regulating the expression of microRNA (miRNA) genes remains largely unexplored. Here, we analyzed DNA methylation changes of miRNA genes using a pair of soybean (Glycine max) near-isogenic lines (NILs) differing in their response to soybean cyst nematode (SCN; Heterodera glycines). Differences in global DNA methylation levels over miRNA genes in response to SCN infection were observed between the isogenic lines. miRNA genes with significant changes in DNA methylation levels in the promoter and primary transcript-coding regions were detected in both lines. In the susceptible isogenic line (NIL-S), 82 differentially methylated miRNAs were identified in response to SCN infection whereas, in the resistant isogenic line (NIL-R), only 16 differentially methylated miRNAs were identified. Interestingly, gma-miR5032, gma-miR5043, gma-miR1520b, and gma-2107-ch16 showed opposite methylation patterns in the isogenic lines. In addition, the miRNA paralogs gma-miR5770a and gma-miR5770b showed hypermethylation and hypomethylation in NIL-S and NIL-R, respectively. Gene expression quantification of gma-miR5032, gma-miR5043, gma-miR1520b, and gma-miR5770a/b and their confirmed targets indicated a role of DNA methylation in regulating miRNA expression and, thus, their targets upon SCN infection. Furthermore, overexpression of these four miRNAs in NIL-S using transgenic hairy root system enhanced plant resistance to SCN to various degrees with a key role observed for miR5032. Together, our results provide new insights into the role of epigenetic mechanisms in controlling miRNA regulatory function during SCN-soybean interactions.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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 Transcriptomics Analysis and Functional Study Reveal Important Role of High-Temperature Stress Response Gene GmHSFA2 During Flower Bud Development of CMS-Based F1 in Soybean

High-temperature (HT) is one of the most important environmental factors that negatively impact the yield of some soybean cytoplasmic male sterility (CMS)-based hybrid (F₁) combinations. The response of soybean to HT, especially at the male organ development stage, is poorly understood. To investigate the molecular mechanisms of the response from soybean CMS-based F₁ male organ to HT, a detailed transcriptomics analysis was performed during flower bud development of soybean HT-tolerant and HT-sensitive CMS-based F₁ combinations (NF₁ and YF₁) under normal-temperature and HT conditions. Obvious HT damage was observed by subjecting YF₁ with HT, such as indehiscent anthers and decreased pollen fertility, whereas the male fertility of NF₁ was normal. In total, 8,784 differentially expressed genes (DEGs) were found to respond to HT stress, which were mainly associated with anther/pollen wall development, carbohydrate metabolism and sugar transport, and auxin signaling. The quantitative real-time PCR (qRT-PCR) analysis and substance content detection also revealed that HT caused male fertility defects in YF₁ by altering pectin metabolism, auxin, and sugar signaling pathways. Most importantly, the sugar signaling-PIF-auxin signaling pathway may underlie the instability of male fertility in YF₁ under HT. Furthermore, HT induced the expression of heat shock factor (HSF) and heat shock protein (HSP) gene families. Overexpression of GmHSFA2 in Arabidopsis can promote the expression of HT protective genes (such as HSP20) by binding to the HSE motifs in their promoters, so as to improve the HT tolerance during flowering. Our results indicated that GmHSFA2 acted as a positive regulator, conferring HT tolerance improvement in soybean CMS-based F₁. GmHSFA2 may be directly involved in the activation of male fertility protection mechanism in the soybean CMS-based F₁ under HT stress.

CircRNA Expression Pattern and ceRNA and miRNA-mRNA Networks Involved in Anther Development in the CMS Line of Brassica campestris

Male-sterile plants provide an important breeding tool for the heterosis of hybrid crops, such as Brassicaceae. In the last decade, circular RNAs (circRNAs), as a novel class of covalently closed and single-stranded endogenous non-coding RNAs (ncRNAs), have received much attention because of their functions as “microRNA (miRNA) sponges” and “competing endogenous RNAs” (ceRNAs). However, the information about circRNAs in the regulation of male-sterility and anther development is limited. In this study, we established the Polima cytoplasm male sterility (CMS) line “Bcpol97-05A”, and the fertile line, “Bcajh97-01B”, in Brassicacampestris L. ssp. chinensis Makino, syn. B. rapa ssp. chinensis, and performed RNA expression profiling comparisons between the flower buds of the sterile line and fertile line by whole-transcriptome sequencing. A total of 31 differentially expressed (DE) circRNAs, 47 DE miRNAs, and 4779 DE mRNAs were identified. By using Cytoscape, the miRNA-mediated regulatory network and ceRNA network were constructed, and the circRNA A02:23507399|23531438 was hypothesized to be an important circRNA regulating anther development at the post-transcriptional level. The gene ontology (GO) analysis demonstrated that miRNAs and circRNAs could regulate the orderly secretion and deposition of cellulose, sporopollenin, pectin, and tryphine; the timely degradation of lipids; and the programmed cell death (PCD) of tapetum cells, which play key roles in anther development. Our study revealed a new circRNA–miRNA–mRNA network, which is involved in the anther development of B. campestris, which enriched the understanding of CMS in flowering plants, and laid a foundation for further study on the functions of circRNAs and miRNAs during anther development.

Identification of Differentially Expressed miRNAs between a Wheat K-type Cytoplasmic Male Sterility Line and Its Near-Isogenic Restorer Line

K-type cytoplasmic male sterility (KCMS) lines were ideal material for three-line hybrid wheat system due to the major role in hybrid wheat production. In this study, the morphology of developing microspore and mature pollen was compared between a KCMS line and its near-isogenic restorer line (KCMS-NIL). The most striking difference is that the microspore was unable to develop into tricellular pollen in the KCMS line. MicroRNA plays vital roles in flowering and gametophyte development. Small RNA sequencing identified a total of 274 known and 401 novel miRNAs differentially expressed between two lines or two developmental stages. Most of miRNAs with high abundance were differentially expressed at the uninucleate stage, and their expression level recovered or remained at the binucleate stage. Further degradome sequencing identified target genes which were mainly enriched in transcription regulation, phytohormone signaling and RNA degradation pathways. Combining with the transcriptome data, a correlation was found between the abnormal anther development, such as postmeiotic mitosis cessation, deformative pollen wall and the chromosome condensation of the vegetative cell, and the alterations in the related miRNA and their targets expression profiles. According to the correlation and pathway analysis, we propose a hypothetic miRNA-mediated network for the control of KCMS restoration.

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

Ogura-type cytoplasmic male sterility (Ogura-CMS) has been widely used in the hybrid breeding industry for cruciferous vegetables. Turnip (Brassica rapa ssp. rapifera) is one of the most important local cruciferous vegetables in China, cultivated for its fleshy root as a flat disc. Here, morphological characteristics of an Ogura-CMS line ‘BY10-2A’ and its maintainer fertile (MF) line ‘BY10-2B’ of turnip were investigated. Ogura-CMS turnip showed a reduction in the size of the fleshy root, and had distinct defects in microspore development and tapetum degeneration during the transition from microspore mother cells to tetrads. Defective microspore production and premature tapetum degeneration during microgametogenesis resulted in short filaments and withered white anthers, leading to complete male sterility of the Ogura-CMS line. Additionally, the mechanism regulating Ogura-CMS in turnip was investigated using inflorescence transcriptome analyses of the Ogura-CMS and MF lines. The de novo assembly resulted in a total of 84,132 unigenes. Among them, 5,117 differentially expressed genes (DEGs) were identified, including 1,339 up- and 3,778 down-regulated genes in the Ogura-CMS line compared to the MF line. A number of functionally known members involved in anther development and microspore formation were addressed in our DEG pool, particularly genes regulating tapetum programmed cell death (PCD), and associated with pollen wall formation. Additionally, 185 novel genes were proposed to function in male organ development based on GO analyses, of which 26 DEGs were genotype-specifically expressed. Our research provides a comprehensive foundation for understanding anther development and the CMS mechanism in turnip.

Mitochondrial genome and transcriptome analysis of five alloplasmic male-sterile lines in Brassica juncea

Background

Alloplasmic lines, in which the nuclear genome is combined with wild cytoplasm, are often characterized by cytoplasmic male sterility (CMS), regardless of whether it was derived from sexual or somatic hybridization with wild relatives. In this study, we sequenced and analyzed the mitochondrial genomes of five such alloplasmic lines in Brassica juncea.

Results

The assembled and annotated mitochondrial genomes of the five alloplasmic lines were found to have virtually identical gene contents. They preserved most of the ancestral mitochondrial segments, and the same candidate male sterility gene (orf108) was found harbored in mitotype-specific sequences. We also detected promiscuous sequences of chloroplast origin that were conserved among plants of the Brassicaceae, and found the RNA editing profiles to vary across the five mitochondrial genomes.

Conclusions

On the basis of our characterization of the genetic nature of five alloplasmic mitochondrial genomes, we speculated that the putative candidate male sterility gene orf108 may not be responsible for the CMS observed in Brassica oxyrrhina and Diplotaxis catholica. Furthermore, we propose the potential coincidence of CMS in alloplasmic lines. Our findings lay the foundation for further elucidation of male sterility gene.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5721-2) contains supplementary material, which is available to authorized users.

Exploration of miRNA-mediated fertility regulation network of cytoplasmic male sterility during flower bud development in soybean

Cytoplasmic male sterility (CMS) plays an important role in the production of soybean hybrid seeds. MicroRNAs (miRNAs) are a class of non-coding endogenous ~ 21 nt small RNAs that play crucial roles in flower and pollen development by targeting genes in plants. To dissect the function of miRNAs in soybean CMS, a total of 558 known miRNAs, 10 novel miRNAs, and 466 target genes were identified in flower buds of the soybean CMS line NJCMS1A and its restorer line NJCMS1C through small RNA sequencing and degradome analysis. In addition, miRNA-mediated editing events were also observed, and the two most frequently observed editing types (A → G and C → U) were validated by cloning and sequencing. And as the base editing occurred, some targets were filtered, such as gma-miR2118b-P6GT with Glyma.08G122000.2. Further integrated analysis of transcriptome and small RNA found some miRNAs and their targets' expression patterns showing a negative correlation, such as gma-miR156b/GmSPL9a and gma-miR4413b/GmPPR. Furthermore, opposite expression pattern was observed between gma-miR156b and GmSPL9 during early stage of flower bud development. Taken together, the regulatory network of gma-miR156b/GmSPL9 and gma-miR4413b/GmPPR with flower bud development in soybean CMS was developed. Most importantly, previous reports showed that these targets might be related to pollen development and male sterility, suggesting that both conserved and species-specific miRNAs might act as regulators of flower bud development in soybean CMS. These findings may provide a better understanding of the miRNA-mediated regulatory networks of CMS mechanisms in soybean.

Comparative analysis of circular RNAs between soybean cytoplasmic male-sterile line NJCMS1A and its maintainer NJCMS1B by high-throughput sequencing

BACKGROUND

Cytoplasmic male sterility (CMS) is a natural phenomenon of pollen abortion caused by the interaction between cytoplasmic genes and nuclear genes. CMS is a simple and effective pollination control system, and plays an important role in crop heterosis utilization. Circular RNAs (circRNAs) are a vital type of non-coding RNAs, which play crucial roles in microRNAs (miRNAs) function and post-transcription control. To explore the expression profile and possible functions of circRNAs in the soybean CMS line NJCMS1A and its maintainer NJCMS1B, high-throughput deep sequencing coupled with RNase R enrichment strategy was conducted.

RESULTS

CircRNA libraries were constructed from flower buds of NJCMS1A and its maintainer NJCMS1B with three biological replicates. A total of 2867 circRNAs were identified, with 1009 circRNAs differentially expressed between NJCMS1A and NJCMS1B based on analysis of high-throughput sequencing. Of the 12 randomly selected circRNAs with different expression levels, 10 showed consistent expression patterns based on high-throughput sequencing and quantitative real-time PCR analyses. Tissue specific expression patterns were also verified with two circRNAs by quantitative real-time PCR. Most parental genes of differentially expressed circRNAs were mainly involved in biological processes such as metabolic process, biological regulation, and reproductive process. Moreover, 83 miRNAs were predicted from the differentially expressed circRNAs, some of which were strongly related to pollen development and male fertility; The functions of miRNA targets were analyzed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), and the target mRNAs were significantly enriched in signal transduction and programmed cell death. Furthermore, a total of 165 soybean circRNAs were predicted to contain at least one internal ribosome entry site (IRES) element and an open reading frame, indicating their potential to encode polypeptides or proteins.

CONCLUSIONS

Our study indicated that the circRNAs might participate in the regulation of flower and pollen development, which could provide a new insight into the molecular mechanisms of CMS in soybean.

Exploration of miRNAs and target genes of cytoplasmic male sterility line in cotton during flower bud development

Cytoplasmic male sterility (CMS) lines provide crucial material to harness heterosis for crop plants, which serves as an important strategy for hybrid seed production. However, the molecular mechanism remains obscure. Although microRNAs (miRNAs) play important roles in vegetative growth and reproductive growth, there are few reports on miRNAs regulating the development of male sterility in Upland cotton. In present study, 12 small RNA libraries were constructed and sequenced for two development stages of flower buds from a CMS line and its maintainer line. Based on the results, 256 novel miRNAs were allocated to 141 new miRNA families, and 77 known miRNAs belonging to 54 conserved miRNA families were identified as well. Comparative analysis revealed that 61 novel and 10 conserved miRNAs were differentially expressed. Further transcriptome analysis identified 232 target genes for these miRNAs, which participated in cellular developmental process, cell death, pollen germination, and sexual reproduction. In addition, expression patterns of typical miRNA and the negatively regulated target genes, such as PPR, ARF, AP2, and AFB, were verified by qRT-PCR in cotton flower buds. These targets were previously reported to be related to reproduction development and male sterility, suggesting that miRNAs might act as regulators of CMS occurrence. Some miRNAs displayed specific expression profiles in special developmental stages of CMS line and its fertile hybrid (F₁). Present study offers new information on miRNAs and their related target genes in exploiting CMS mechanism, and revealing the miRNA regulatory networks in Upland cotton.

Non-coding RNAs and plant male sterility: current knowledge and future prospects

Latest outcomes assign functional role to non-coding (nc) RNA molecules in regulatory networks that confer male sterility to plants. Male sterility in plants offers great opportunity for improving crop performance through application of hybrid technology. In this respect, cytoplasmic male sterility (CMS) and sterility induced by photoperiod (PGMS)/temperature (TGMS) have greatly facilitated development of high-yielding hybrids in crops. Participation of non-coding (nc) RNA molecules in plant reproductive development is increasingly becoming evident. Recent breakthroughs in rice definitively associate ncRNAs with PGMS and TGMS. In case of CMS, the exact mechanism through which the mitochondrial ORFs exert influence on the development of male gametophyte remains obscure in several crops. High-throughput sequencing has enabled genome-wide discovery and validation of these regulatory molecules and their target genes, describing their potential roles performed in relation to CMS. Discovery of ncRNA localized in plant mtDNA with its possible implication in CMS induction is intriguing in this respect. Still, conclusive evidences linking ncRNA with CMS phenotypes are currently unavailable, demanding complementing genetic approaches like transgenics to substantiate the preliminary findings. Here, we review the recent literature on the contribution of ncRNAs in conferring male sterility to plants, with an emphasis on microRNAs. Also, we present a perspective on improved understanding about ncRNA-mediated regulatory pathways that control male sterility in plants. A refined understanding of plant male sterility would strengthen crop hybrid industry to deliver hybrids with improved performance.

The Role of Non-Coding RNAs in Cytoplasmic Male Sterility in Flowering Plants

The interactions between mitochondria and nucleus substantially influence plant development, stress response and morphological features. The prominent example of a mitochondrial-nuclear interaction is cytoplasmic male sterility (CMS), when plants produce aborted anthers or inviable pollen. The genes responsible for CMS are located in mitochondrial genome, but their expression is controlled by nuclear genes, called fertility restorers. Recent explosion of high-throughput sequencing methods enabled to study transcriptomic alterations in the level of non-coding RNAs under CMS biogenesis. We summarize current knowledge of the role of nucleus encoded regulatory non-coding RNAs (long non-coding RNA, microRNA as well as small interfering RNA) in CMS. We also focus on the emerging data of non-coding RNAs encoded by mitochondrial genome and their possible involvement in mitochondrial-nuclear interactions and CMS development.

Genome-Wide Analysis of DNA Methylation During Ovule Development of Female-Sterile Rice fsv1

The regulation of female fertility is an important field of rice sexual reproduction research. DNA methylation is an essential epigenetic modification that dynamically regulates gene expression during development processes. However, few reports have described the methylation profiles of female-sterile rice during ovule development. In this study, ovules were continuously acquired from the beginning of megaspore mother cell meiosis until the mature female gametophyte formation period, and global DNA methylation patterns were compared in the ovules of a high-frequency female-sterile line (fsv1) and a wild-type rice line (Gui99) using whole-genome bisulfite sequencing (WGBS). Profiling of the global DNA methylation revealed hypo-methylation, and 3471 significantly differentially methylated regions (DMRs) were observed in fsv1 ovules compared with Gui99. Based on functional annotation and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis of differentially methylated genes (DMGs), we observed more DMGs enriched in cellular component, reproduction regulation, metabolic pathway, and other pathways. In particular, many ovule development genes and plant hormone-related genes showed significantly different methylation patterns in the two rice lines, and these differences may provide important clues for revealing the mechanism of female gametophyte abortion.

Genome-wide comparative analysis of DNA methylation between soybean cytoplasmic male-sterile line NJCMS5A and its maintainer NJCMS5B

BACKGROUND

DNA methylation is an important epigenetic modification. It can regulate the expression of many key genes without changing the primary structure of the genomic DNA, and plays a vital role in the growth and development of the organism. The genome-wide DNA methylation profile of the cytoplasmic male sterile (CMS) line in soybean has not been reported so far.

RESULTS

In this study, genome-wide comparative analysis of DNA methylation between soybean CMS line NJCMS5A and its maintainer NJCMS5B was conducted by whole-genome bisulfite sequencing. The results showed 3527 differentially methylated regions (DMRs) and 485 differentially methylated genes (DMGs), including 353 high-credible methylated genes, 56 methylated genes coding unknown protein and 76 novel methylated genes with no known function were identified. Among them, 25 DMRs were further validated that the genome-wide DNA methylation data were reliable through bisulfite treatment, and 9 DMRs were confirmed the relationship between DNA methylation and gene expression by qRT-PCR. Finally, 8 key DMGs possibly associated with soybean CMS were identified.

CONCLUSIONS

Genome-wide DNA methylation profile of the soybean CMS line NJCMS5A and its maintainer NJCMS5B was obtained for the first time. Several specific DMGs which participated in pollen and flower development were further identified to be probably associated with soybean CMS. This study will contribute to further understanding of the molecular mechanism behind soybean CMS.

Uncovering Male Fertility Transition Responsive miRNA in a Wheat Photo-Thermosensitive Genic Male Sterile Line by Deep Sequencing and Degradome Analysis

MicroRNAs (miRNAs) are endogenous small RNAs which play important negative regulatory roles at both the transcriptional and post-transcriptional levels in plants. Wheat is the most commonly cultivated plant species worldwide. In this study, RNA-seq analysis was used to examine the expression profiles of miRNA in the spikelets of photo-thermosenisitive genic male sterile (PTGMS) wheat line BS366 during male fertility transition. Through mapping on their corresponding precursors, 917-7,762 novel miRNAs were found in six libraries. Six novel miRNAs were selected for examination of their secondary structures and confirmation by stem-loop RT-PCR. In a differential expression analysis, 20, 22, and 58 known miRNAs exhibited significant differential expression between developmental stages 1 (secondary sporogenous cells had formed), 2 (all cells layers were present and mitosis had ceased), and 3 (meiotic division stage), respectively, of fertile and sterile plants. Some of these differential expressed miRNAs, such as tae-miR156, tae-miR164, tae-miR171, and tae-miR172, were shown to be associated with their targets. These targets were previously reported to be related to pollen development and/or male sterility, indicating that these miRNAs and their targets may be involved in the regulation of male fertility transition in the PTGMS wheat line BS366. Furthermore, target genes of miRNA cleavage sites were validated by degradome sequencing. In this study, a possible signal model for the miRNA-mediated signaling pathway during the process of male fertility transition in the PTGMS wheat line BS366 was developed. This study provides a new perspective for understanding the roles of miRNAs in male fertility in PTGMS lines of wheat.

Identification and Functional Analysis of microRNAs Involved in the Anther Development in Cotton Genic Male Sterile Line Yu98-8A

Hybrid vigor contributes in a large way to the yield and quality of cotton (Gossypium hirsutum) fiber. Although microRNAs play essential regulatory roles in flower induction and development, it is still unclear if microRNAs are involved in male sterility, as the regulatory molecular mechanisms of male sterility in cotton need to be better defined. In this study, two independent small RNA libraries were constructed and sequenced from the young buds collected from the sporogenous cell formation to the meiosis stage of the male sterile line Yu98-8A and the near-isogenic line. Sequencing revealed 1588 and 1536 known microRNAs and 347 and 351 novel miRNAs from male sterile and male fertile libraries, respectively. MicroRNA expression profiles revealed that 49 conserved and 51 novel miRNAs were differentially expressed. Bioinformatic and degradome analysis indicated the regulatory complexity of microRNAs during flower induction and development. Further RT-qPCR and physiological analysis indicated that, among the different Kyoto Encyclopedia Gene and Genomes pathways, indole-3-acetic acid and gibberellic acid signaling transduction pathways may play pivotal regulatory functions in male sterility.

Cytoplasmic male sterility (CMS) in hybrid breeding in field crops

A comprehensive understanding of CMS/Rf system enabled by modern omics tools and technologies considerably improves our ability to harness hybrid technology for enhancing the productivity of field crops. Harnessing hybrid vigor or heterosis is a promising approach to tackle the current challenge of sustaining enhanced yield gains of field crops. In the context, cytoplasmic male sterility (CMS) owing to its heritable nature to manifest non-functional male gametophyte remains a cost-effective system to promote efficient hybrid seed production. The phenomenon of CMS stems from a complex interplay between maternally-inherited (mitochondrion) and bi-parental (nucleus) genomic elements. In recent years, attempts aimed to comprehend the sterility-inducing factors (orfs) and corresponding fertility determinants (Rf) in plants have greatly increased our access to candidate genomic segments and the cloned genes. To this end, novel insights obtained by applying state-of-the-art omics platforms have substantially enriched our understanding of cytoplasmic-nuclear communication. Concomitantly, molecular tools including DNA markers have been implicated in crop hybrid breeding in order to greatly expedite the progress. Here, we review the status of diverse sterility-inducing cytoplasms and associated Rf factors reported across different field crops along with exploring opportunities for integrating modern omics tools with CMS-based hybrid breeding.