Comparative expression profiling of miRNAs between the cytoplasmic male sterile line MeixiangA and its maintainer line MeixiangB during rice anther development

Identification of differentially expressed miRNAs between male sterile and fertile floral buds in watermelon (Citrullus lanatus L.) via high-throughput sequencing

This experiment used floral buds from watermelon genic male sterile dual-purpose lines as materials to explore the differentially expressed miRNAs (DEMs) between male fertile and sterile floral buds of watermelon. Paraffin sectioning technology was employed for a cytological analysis, and small RNA sequencing was used to explore miRNAs related to anther or pollen development. Cytological analysis indicated that the abnormal development of tapetal cells may cause microspore abortion. Small RNA sequencing identified a total of 314 miRNAs (29 known and 285 novel, which belonged to 12 and 61 miRNA families, respectively) in floral buds. Differential expression revealed 36 (5 known and 31 novel) DEMs between male fertile and sterile buds, 7 and 29 of which were up-regulated and down-regulated, respectively. Target genes analysis showed that the 36 DEMs were predicted to target 577 genes, and these targets might participate in various biological processes, such as response to metal ions, floral organ development, stamen development, anther development, pollen maturation, and programmed cell death. Moreover, pathway analysis indicated that these genes were mainly enriched in purine metabolism, starch and sucrose metabolism, RNA transport, and other pathways. In addition, the 55 miRNA-target modules, including 3 known and 16 novel miRNAs with 30 target genes, might be related to anther or pollen development in watermelon. Our findings provide important miRNA-target modules related to watermelon anther or pollen development and can lay the foundation for biological functional analysis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04084-6.

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.

Insights into microRNA-mediated interaction and regulation of metabolites in tomato

microRNAs direct regulation of various metabolic pathways in plants and animals. miRNAs may be useful in developing novel/elite genotypes, with enhanced metabolites and disease resistance. We examined miRNAs in tomato. In tomato, miRNAs in the carotenoid pathway have not been fully elucidated. We examined the potential role of miRNAs in biosynthesis of carotenoids, transcript profiling of miRNAs and their possible targets (genes and transcription factors) at different development stages of tomato using stem-loop PCR and RT-qPCR. We also identified miRNAs targeting key flavonoid genes, such as chalcone isomerase (CHI), and dihydroflavonol-4-reductase (DFR). Distinct expression profiles of miRNAs and their targets were found in fruits of three tomato accessions, suggesting carotenoid regulation by miRNAs at various stages of fruit development. This was also confirmed using HPLC of the carotenoids. The present study may help in understanding possible regulation of carotenoid biosynthesis. The identified miRNAs can be exploited to enhance biosynthesis of different carotenoids in plants.

Cytological observation and RNA-seq analysis reveal novel miRNAs high expression associated with the pollen fertility of neo-tetraploid rice

BACKGROUND

Neo-tetraploid rice lines exhibit high fertility and strong heterosis and harbor novel specific alleles, which are useful germplasm for polyploid rice breeding. However, the mechanism of the fertility associated with miRNAs remains unknown. In this study, a neo-tetraploid rice line, termed Huaduo21 (H21), was used. Cytological observation and RNA-sequencing were employed to identify the fertility-related miRNAs in neo-tetraploid rice.

RESULTS

H21 showed high pollen fertility (88.08%), a lower percentage of the pollen mother cell (PMC) abnormalities, and lower abnormalities during double fertilization and embryogenesis compared with autotetraploid rice. A total of 166 non-additive miRNAs and 3108 non-additive genes were detected between H21 and its parents. GO and KEGG analysis of non-additive genes revealed significant enrichments in the DNA replication, Chromosome and associated proteins, and Replication and repair pathways. Comprehensive multi-omics analysis identified 32 pairs of miRNA/target that were associated with the fertility in H21. Of these, osa-miR408-3p and osa-miR528-5p displayed high expression patterns, targeted the phytocyanin genes, and were associated with high pollen fertility. Suppression of osa-miR528-5p in Huaduo1 resulted in a low seed set and a decrease in the number of grains. Moreover, transgenic analysis implied that osa-MIR397b-p3, osa-miR5492, and osa-MIR5495-p5 might participate in the fertility of H21.

CONCLUSION

Taken together, the regulation network of fertility-related miRNAs-targets pairs might contribute to the high seed setting in neo-tetraploid rice. These findings enhance our understanding of the regulatory mechanisms of pollen fertility associated with miRNAs in neo-tetraploid rice.

High-temperature stress in crops: male sterility, yield loss and potential remedy approaches

Global food security is one of the utmost essential challenges in the 21st century in providing enough food for the growing population while coping with the already stressed environment. High temperature (HT) is one of the main factors affecting plant growth, development and reproduction and causes male sterility in plants. In male reproductive tissues, metabolic changes induced by HT involve carbohydrates, lipids, hormones, epigenetics and reactive oxygen species, leading to male sterility and ultimately reducing yield. Understanding the mechanism and genes involved in these pathways during the HT stress response will provide a new path to improve crops by using molecular breeding and biotechnological approaches. Moreover, this review provides insight into male sterility and integrates this with suggested strategies to enhance crop tolerance under HT stress conditions at the reproductive stage.

Transcriptional Regulation of zma-MIR528a by Action of Nitrate and Auxin in Maize

In recent years, miR528, a monocot-specific miRNA, has been assigned multifaceted roles during development and stress response in several plant species. However, the transcription regulation and the molecular mechanisms controlling MIR528 expression in maize are still poorly explored. Here we analyzed the zma-MIR528a promoter region and found conserved transcription factor binding sites related to diverse signaling pathways, including the nitrate (TGA1/4) and auxin (AuxRE) response networks. Accumulation of both pre-miR528a and mature miR528 was up-regulated by exogenous nitrate and auxin treatments during imbibition, germination, and maize seedling establishment. Functional promoter analyses demonstrated that TGA1/4 and AuxRE sites are required for transcriptional induction by both stimuli. Overall, our findings of the nitrogen- and auxin-induced zma-MIR528a expression through cis-regulatory elements in its promoter contribute to the knowledge of miR528 regulome.

Multiple Functions of MiRNAs in Brassica napus L

The worldwide climate changes every year due to global warming, waterlogging, drought, salinity, pests, and pathogens, impeding crop productivity. Brassica napus is one of the most important oil crops in the world, and rapeseed oil is considered one of the most health-beneficial edible vegetable oils. Recently, miRNAs have been found and confirmed to control the expression of targets under disruptive environmental conditions. The mechanism is through the formation of the silencing complex that mediates post-transcriptional gene silencing, which pairs the target mRNA and target cleavage and/or translation inhibition. However, the functional role of miRNAs and targets in B. napus is still not clarified. This review focuses on the current knowledge of miRNAs concerning development regulation and biotic and abiotic stress responses in B. napus. Moreover, more strategies for miRNA manipulation in plants are discussed, along with future perspectives, and the enormous amount of transcriptome data available provides cues for miRNA functions in B. napus. Finally, the construction of the miRNA regulatory network can lead to the significant development of climate change-tolerant B. napus through miRNA manipulation.

Whole-Transcriptome Analysis Reveals Autophagy Is Involved in Early Senescence of zj-es Mutant Rice

Senescence is a necessary stage of plant growth and development, and the early senescence of rice will lead to yield reduction and quality decline. However, the mechanisms of rice senescence remain obscure. In this study, we characterized an early-senescence rice mutant, designated zj-es (ZheJing-early senescence), which was derived from the japonica rice cultivar Zhejing22. The mutant zj-es exhibited obvious early-senescence phenotype, such as collapsed chloroplast, lesions in leaves, declined fertility, plant dwarf, and decreased agronomic traits. The ZJ-ES gene was mapped in a 458 kb-interval between the molecular markers RM5992 and RM5813 on Chromosome 3, and analysis suggested that ZJ-ES is a novel gene controlling rice early senescence. Subsequently, whole-transcriptome RNA sequencing was performed on zj-es and its wild-type rice to dissect the underlying molecular mechanism for early senescence. Totally, 10,085 differentially expressed mRNAs (DEmRNAs), 1,253 differentially expressed lncRNAs (DElncRNAs), and 614 differentially expressed miRNAs (DEmiRNAs) were identified, respectively, in different comparison groups. Based on the weighted gene co-expression network analysis (WGCNA), the co-expression turquoise module was found to be the key for the occurrence of rice early senescence. Furthermore, analysis on the competing endogenous RNA (CeRNA) network revealed that 14 lncRNAs possibly regulated 16 co-expressed mRNAs through 8 miRNAs, and enrichment analysis showed that most of the DEmRNAs and the targets of DElncRNAs and DEmiRNAs were involved in reactive oxygen species (ROS)-triggered autophagy-related pathways. Further analysis showed that, in zj-es, ROS-related enzyme activities were markedly changed, ROS were largely accumulated, autophagosomes were obviously observed, cell death was significantly detected, and lesions were notably appeared in leaves. Totally, combining our results here and the remaining research, we infer that ROS-triggered autophagy induces the programmed cell death (PCD) and its coupled early senescence in zj-es mutant rice.

Comparative transcriptome analysis reveals differential gene expression in sterile and fertile rubber tree varieties during flower bud differentiation

Plant male sterility (MS) is an important agronomic trait that provides an efficient tool for hybridization and heterosis utilization of crops. Based on phenotypic and cytological observations, our study performed a multi-comparison transcriptome analysis strategy on multiple sterile and fertile rubber tree varieties using RNA-seq. Compared with the male-fertile varieties, a total of 1590 differentially expressed genes (DEGs) were detected in male-sterile varieties, including 970 up-regulated and 620 down-regulated transcripts in sterile varieties. Key DEGs were further assessed focusing on anther development, microsporogenesis and plant hormone metabolism. Twenty DEGs were selected randomly to validate transcriptome data using quantitative real-time PCR (qRT-PCR). Eleven key genes were subjected to expression pattern analysis using qRT-PCR and fluorescence in situ hybridization. Among them, nine genes, i.e., A6, GAI1, ACA7, TKPR1, CYP704B1, XTH26, MS1, MS35 and MYB33, that regulate callose metabolism, pollen wall formation, tapetum and microspores development were identified as candidate male-sterile genes. These findings provide insights into the molecular mechanism of male sterility in rubber tree.

Dynamic expression of miRNAs and functional analysis of target genes involved in the response to male sterility of the wheat line YS3038

Thermosensitive cytoplasmic male sterile (TCMS) lines play an important role in wheat breeding, heterosis utilization, and germplasm innovation. MicroRNAs (miRNAs) can regulate the expression level of target genes by inhibiting the translation of these genes. YS3038 is a wheat TCMS line. In this study, the fertility conversion mechanism of YS3038 was studied by examining the abortion characteristics of YS3038, the regulation pattern of miRNAs and the target genes of miRNAs in YS3038. MiRNA-seq was performed on three important stages of YS3038 under sterile and fertile conditions. Then, the clean reads were aligned with some databases to filter other ncRNAs and repeats. The known miRNAs and novel miRNAs were predicted by sequence comparison with known miRNAs from miRbase. Differential expression of miRNAs between different stages and between different fertile conditions was analyzed, and functional analysis of target genes with opposite expression patterns as those of the miRNAs was conducted. The Ubisch bodies and microspores of sterile anthers were covered with filamentous materials. The degradation of the tapetum cells, the chloroplast structure of endothecium cells, and the microspore structure were abnormal. Microspore development was hindered from the late uninucleate stage to the binucleate stage. Twenty, 52, and 68 differentially expressed miRNAs (DEmiRs) were identified at the early uninucleate, late uninucleate, and binucleate stages, respectively, and there were 0, 7, and 72 differentially expressed target genes (DETGs), respectively, at these three stages. At the binucleate stage, 29 DEmiRs had 41 target mRNAs in total, and the expression patterns of the 41 target mRNAs were opposite to those of the 29 miRNAs. Fifteen significantly enriched KEGG pathways were associated with the 41 target mRNAs. Leucine-rich repeat receptor-like kinases (LRR-RLKs) play important roles in plant developmental and physiological processes. Some studies have shown that the expression of LRR-RLKs is related to the differentiation of microsporocytes and tapetum cells and to male sterility. An LRR-RLK (TaeRPK) gene was silenced by the barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) method, and the seed setting rates of the TaeRPK-silenced plants (3.51%) were significantly lower than those of the negative control plants (88.78%) (P < 0.01). Thus, the TaeRPK gene is likely to be involved in the fertility conversion of YS3038.

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.

Whole-Transcriptome RNA Sequencing Reveals the Global Molecular Responses and CeRNA Regulatory Network of mRNAs, lncRNAs, miRNAs and circRNAs in Response to Salt Stress in Sugar Beet (Beta vulgaris)

Sugar beet is an important sugar-yielding crop with some tolerance to salt, but the mechanistic basis of this tolerance is not known. In the present study, we have used whole-transcriptome RNA-seq and degradome sequencing in response to salt stress to uncover differentially expressed (DE) mRNAs, microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) in both leaves and roots. A competitive endogenous RNA (ceRNA) network was constructed with the predicted DE pairs, which revealed regulatory roles under salt stress. A functional analysis suggests that ceRNAs are implicated in copper redistribution, plasma membrane permeability, glycometabolism and energy metabolism, NAC transcription factor and the phosphoinositol signaling system. Overall, we conducted for the first time a full transcriptomic analysis of sugar beet under salt stress that involves a potential ceRNA network, thus providing a basis to study the potential functions of lncRNAs/circRNAs.

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.

Profiling miRNA expression in photo-thermo-sensitive male genic sterility line (PTGMS) PA64S under high and low temperature

Peiai64S (PA64S) is a photo-thermo-sensitive genic male sterile line (PTGMS), with wide application in hybrid seed production in rice (Oryza sativa L.). Micro-RNAs are 21-24 nt, endogenously expressed small RNAs that have been characterized in various developmental stages of rice, but none have been studied with respect to the regulation of TGMS in rice. Here, we employed high-throughput sequencing to identify expression profiles of miRNAs in the anthers of PA64S at high (PA64S-H) and low temperature (PA64S-L). Two small RNA libraries from PA64S-H and PA64-L anthers were sequenced, and 263 known and 321 novel candidate miRNAs were identified. Based on the number of sequencing reads, a total of 133 known miRNAs were found to be differentially expressed between PA64S-H and PA64S-L. Target prediction showed that the target genes encode MYB and TCP transcription factors, and bHLH proteins. These target genes are related to pollen development and male sterility, suggesting that miRNA/targets may play roles in regulating TGMS in rice. Further, starch and sucrose metabolism pathways, sphingolipid metabolism, arginine and proline metabolism, and plant hormone signal transduction pathways were enriched by KEGG pathway annotation. These findings contribute to our understanding of the role of miRNAs during anther development and TGMS occurrence in rice.

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.

Small RNA profiling from meiotic and post-meiotic anthers reveals prospective miRNA-target modules for engineering male fertility in sorghum

Understanding male gametophyte development is essential to augment hybrid production in sorghum. Although small RNAs are known to critically influence anther/pollen development, their roles in sorghum reproduction have not been deciphered yet. Here, we report small RNA profiling and high-confidence annotation of microRNAs (miRNAs) from meiotic and post-meiotic anthers in sorghum. We identified 262 miRNAs (82 known and 180 novel), out of which 58 (35 known and 23 novel) exhibited differential expression between two stages. Out of 35 differentially expressed known miRNAs, 13 are known to regulate anther/pollen development in other plant species. We also demonstrated conserved spatiotemporal patterns of 21- and 24-nt phasiRNAs and their respective triggers, miR2118 and miR2275, in sorghum anthers as evidenced in other monocots. miRNA target identification yielded 5622 modules, of which 46 modules comprising 16 known and 8 novel miRNA families with 38 target genes are prospective candidates for engineering male fertility in grasses.

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.

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 Transcriptome Analysis between Fertile and CMS Flower Buds in Wucai (Brassica campestris L.)

Background

Wucai (Brassica campestris L. ssp. chinensis var. rosularis Tsen) is a variant of nonheading Chinese cabbage (Brassica campestris L.), which is one of the major vegetables in China. Cytoplasmic male sterility (CMS) has been used for Wucai breeding in recent years. However, the underlying molecular mechanism of Wucai CMS remains unclear. In this study, the phenotypic and cytological features of Wucai CMS were observed by anatomical analysis, and a comparative transcriptome analysis was carried out to identify genes related to male sterility using Illumina RNA sequencing technology (RNA-Seq).

Results

Microscopic observation demonstrated that tapetum development was abnormal in the CMS line, which failed to produce fertile pollen. Bioinformatics analysis detected 4430 differentially expressed genes (DEGs) between the fertile and sterile flower buds. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to better understand the functions of these DEGs. Among the DEGs, 35 genes (53 DEGS) were implicated in anther and pollen development, and 11 genes were involved in pollen cell wall formation and modification; most of these showed downregulated expression in sterile buds. In addition, several genes related to tapetum development (A6, AMS, MS1, MYB39, and TSM1) and a few genes annotated to flowering (CO, AP3, VIN3, FLC, FT, and AGL) were detected and confirmed by qRT-PCR as being expressed at the meiosis, tetrad, and uninucleate microspore stages, thus implying possible roles in specifying or determining the fate and development of the tapetum, male gametophyte and stamen. Moreover, the top four largest transcription factor families (MYB, bHLH, NAC and WRKY) were analyzed, and most showed reduced expression in sterile buds. These differentially expressed transcription factors might result in abortion of pollen development in Wucai.

Conclusion

The present comparative transcriptome analysis suggested that many key genes and transcription factors involved in anther development show reduced gene expression patterns in the CMS line, which might contribute to male sterility in Wucai. This study provides valuable information for a better understanding of CMS molecular mechanisms and functional genome studies in Wucai.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5331-4) contains supplementary material, which is available to authorized users.

Transcriptome and Hormone Comparison of Three Cytoplasmic Male Sterile Systems in Brassica napus

The interaction between plant mitochondria and the nucleus markedly influences stress responses and morphological features, including growth and development. An important example of this interaction is cytoplasmic male sterility (CMS), which results in plants producing non-functional pollen. In current research work, we compared the phenotypic differences in floral buds of different Brassica napus CMS (Polima, Ogura, Nsa) lines with their corresponding maintainer lines. By comparing anther developmental stages between CMS and maintainer lines, we identified that in the Nsa CMS line abnormality occurred at the tetrad stage of pollen development. Phytohormone assays demonstrated that IAA content decreased in sterile lines as compared to maintainer lines, while the total hormone content was increased two-fold in the S₂ stage compared with the S₁ stage. ABA content was higher in the S₁ stage and exhibited a two-fold decreasing trend in S₂ stage. Sterile lines however, had increased ABA content at both stages compared with the corresponding maintainer lines. Through transcriptome sequencing, we compared differentially expressed unigenes in sterile and maintainer lines at both (S₁ and S₂) developmental stages. We also explored the co-expressed genes of the three sterile lines in the two stages and classified these genes by gene function. By analyzing transcriptome data and validating by RT-PCR, it was shown that some transcription factors (TFs) and hormone-related genes were weakly or not expressed in the sterile lines. This research work provides preliminary identification of the pollen abortion stage in Nsa CMS line. Our focus on genes specifically expressed in sterile lines may be useful to understand the regulation of CMS.

Cellular and molecular characteristics of pollen abortion in chrysanthemum cv. Kingfisher

Microspore degeneration at the tetrad stage is associated with tapetum degeneration retardation. Some genes and proteins related to cell senescence and death are the key factors for pollen abortion. Chrysanthemum (Chrysanthemum morifolium) is a major floriculture crop in the world, but pollen contamination is an urgent problem to be solved in chrysanthemum production. C. morifolium 'Kingfisher' is a chrysanthemum cultivar that does not contain any pollen in mature anthers, thus it is a very important material for developing chrysanthemum without pollen contamination. However, the mechanism of its pollen abortion remains unclear. In this study, the cellular and molecular mechanisms of 'Kingfisher' pollen abortion were investigated using transmission electron microscopy, RNA sequencing, isobaric tags for relative and absolute quantitation, and bioinformatics. It was found that the meiosis of microspore mother cells was normal before the tetrad stage, the microspores began to degenerate at the tetrad stage, and no microspores were observed in the anthers after the tetrad stage. In addition, transcriptomic and proteomic analyses showed that some genes and proteins related to cell senescence and death were identified to be implicated in chrysanthemum pollen abortion. These results indicated that the tetrad stage was the main period of pollen abortion, and the genes and proteins related to cell senescence and death contributed to pollen abortion. These add to our understanding of chrysanthemum pollen abortion and will be helpful for development of flowers without pollen contamination in the future.

Genome-wide analysis of the rice PPR gene family and their expression profiles under different stress treatments

Background

Pentatricopeptide-repeat proteins (PPRs) are characterized by tandem arrays of a degenerate 35-amino-acid (PPR motifs), which can bind RNA strands and participate in post-transcription. PPR proteins family is one of the largest families in land plants and play important roles in organelle RNA metabolism and plant development. However, the functions of PPR genes involved in biotic and abiotic stresses of rice (Oryza sativa L.) remain largely unknown.

Results

In the present study, a comprehensive genome-wide analysis of PPR genes was performed. A total of 491 PPR genes were found in the rice genome, of which 246 PPR genes belong to the P subfamily, and 245 genes belong to the PLS subfamily. Gene structure analysis showed that most PPR genes lack intron. Chromosomal location analysis indicated that PPR genes were widely distributed in all 12 rice chromosomes. Phylogenetic relationship analysis revealed the distinct difference between the P and PLS subfamilies. Many PPR proteins are predicted to target chloroplasts or mitochondria, and a PPR protein (LOC_Os10g34310) was verified to localize in mitochondria. Furthermore, three PPR genes (LOC_Os03g17634,LOC_Os07g40820,LOC_Os04g51350) were verified as corresponding miRNA targets. The expression pattern analysis showed that many PPR genes could be induced under biotic and abiotic stresses. Finally, seven PPR genes were confirmed with their expression patterns under salinity or drought stress.

Conclusions

We found 491 PPR genes in the rice genome, and our genes structure analysis and syntenic analysis indicated that PPR genes might be derived from amplification by retro-transposition. The expression pattern present here suggested that PPR proteins have crucial roles in response to different abiotic stresses in rice. Taken together, our study provides a comprehensive analysis of the PPR gene family and will facilitate further studies on their roles in rice growth and development.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5088-9) contains supplementary material, which is available to authorized users.

The key genes and pathways related to male sterility of eggplant revealed by comparative transcriptome analysis

BACKGROUND

Male sterility (MS) is an effective tool for hybrid production. Although MS has been widely reported in other plants, such as Arabidopsis and rice, the molecular mechanism of MS in eggplant is largely unknown. To understand the mechanism, the comparative transcriptomic file of MS line and its maintainer line was analyzed with the RNA-seq technology.

RESULTS

A total of 11,7695 unigenes were assembled and 19,652 differentially expressed genes (DEGs) were obtained. The results showed that 1,716 DEGs were shared in the three stages. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that these DEGs were mainly involved in oxidation-reduction, carbohydrate and amino acid metabolism. Moreover, transcriptional regulation was also the impact effector for MS and anther development. Weighted correlation network analysis (WGCNA) showed two modules might be responsible for MS, which was similar to hierarchical cluster analysis.

CONCLUSIONS

A number of genes and pathways associated with MS were found in this study. This study threw light on the molecular mechanism of MS and identified several key genes related to MS in eggplant.

Uncovering key small RNAs associated with gametocidal action in wheat

Gametocidal (Gc) chromosomes can kill gametes that lack them by causing chromosomal breakage to ensure their preferential transmission, and they have been exploited in genetic breeding. The present study investigated the possible roles of small RNAs (sRNAs) in Gc action. By sequencing two small RNA libraries from the anthers of Triticum aestivum cv. Chinese Spring (CS) and the Chinese Spring-Gc 3C chromosome monosomic addition line (CS-3C), we identified 239 conserved and 72 putative novel miRNAs, including 135 differentially expressed miRNAs. These miRNAs were predicted to target multiple genes with various molecular functions relevant to the features of Gc action, including sterility and genome instability. The transgenic overexpression of miRNA, which was up-regulated in CS-3C, reduced rice fertility. The CS-3C line exhibited a genome-wide reduction in 24 nt siRNAs compared with that of the CS line, particularly in transposable element (TE) and repetitive DNA sequences. Corresponding to this reduction, the bisulfite sequencing analysis of four retro-TE sequences showed a decrease in CHH methylation, typical of RNA-directed DNA methylation (RdDM). These results demonstrate that both miRNA-directed regulation of gene expression and siRNA-directed DNA methylation of target TE loci could play a role in Gc action.

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.

Comprehensive analysis of small RNAs expressed in developing male strobili of Cryptomeria japonica

Deep sequencing of small RNAs (sRNAs) in developing male strobili of second-generation offspring originating from a nuclear genic male sterile tree of Cryptomeria japonica were performed to characterize sRNA populations in the male strobili at early pollen developmental stages. Comparing to sequences of microRNA (miRNA) families of plant species and sRNAs expressed in the reproductive organs of representative vascular plants, 37 conserved miRNA families were detected, of which eight were ubiquitously expressed in the reproductive organs of land plant species. In contrast, miR1083 was common in male reproductive organs of gymnosperm species but absent in angiosperm species. In addition to conserved miRNAs, 199 novel miRNAs candidates were predicted. The expression patterns of the obtained sRNAs were further investigated to detect the differentially expressed (DE) sRNAs between genic male sterile and fertile individuals. A total of 969 DE sRNAs were obtained and only three known miRNA families were included among them. These results suggest that both conserved and species-specific sRNAs contribute to the development of male strobili in C. japonica.

Effective enhancement of resistance to Phytophthora infestans by overexpression of miR172a and b in Solanum lycopersicum

Overexpression of miR172a and b in tomato ( Solanum lycopersicum ) Zaofen No. 2 increased resistance to Phytophthora infestans infection by suppressing of an AP2/ERF transcription factor. The miR172 family has been shown to participate in the growth phase transition, flowering time control, abiotic and biotic stresses by regulating the expression of a small group of AP2/ERF transcription factors. In this study, the precursors of miR172a and b were cloned from tomato, Solanum pimpinellifolium L3708. We used the degradome sequencing to determine the cleavage site of miR172 to a member of the AP2/ERF transcription factor family (Solyc11g072600.1.1). qRT-PCR results showed that the expression of AP2/ERF was negatively correlated with the expression of miR172 in S. pimpinellifolium L3708 infected with Phytophthora infestans. Overexpression of miR172a and b in S. lycopersicum Zaofen No. 2 conferred greater resistance to P. infestans infection, as evidenced by decreased disease index, lesion sizes, and P. infestans abundance. The SOD and POD play important roles in scavenging late massive ROS in plant-pathogen interaction. Malonaldehyde (MDA) is widely recognized as an indicator of lipid peroxidation. Membrane damage in plants can be estimated by measuring leakage of electrolytes, which is evaluated by determining relative electrolyte leakage (REL). Less H₂O₂ and O₂^-, higher activities of POD and SOD, less MDA content and REL, and higher chlorophyll content and photosynthetic rate were also shown in transgenic plants after inoculation with P. infestans. Our results constitute the first step towards further investigations into the biological function and molecular mechanism of miR172-mediated silencing of AP2/ERF transcription factors in S. lycopersicum-P. infestans interaction and provide a candidate gene for breeding to enhance biotic stress-resistance in S. lycopersicum.

Molecular Control of Male Fertility for Crop Hybrid Breeding

In many plant species, male-sterile female lines with cytoplasmic male sterility (CMS) or nuclear-controlled environment-sensitive genic male sterility (EGMS) have long been used to efficiently produce hybrids that harness hybrid vigor or heterosis. However, the underlying molecular mechanisms for these applications have only recently been uncovered in a few species. We provide here an update on the understanding of cytoplasmic-nuclear communication based on the discovery of mitochondrial CMS genes and their corresponding nuclear fertility determinants. Recent findings that uncover diverse mechanisms such as epigenetic, transcriptional, and post-transcriptional controls of EGMS by temperature and photoperiod signals are also reviewed. Furthermore, translational research that applies basic knowledge of plant male fertility control to hybrid seed production practice is highlighted.

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.

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.

Genetic Subtraction Profiling Identifies Candidate miRNAs Involved in Rice Female Gametophyte Abortion

The female gametophyte is an important participant in the sexual reproduction of plants. The molecular mechanism of its development has received much attention in recent years. As important regulators of gene expression, miRNAs have been certified to play a significant role in many biological processes of plants, including sexual reproduction. In this study, to investigate the potential regulatory effects of miRNAs on rice female gametophyte abortion, we used the high-throughput sequencing method to compare the miRNA transcriptome in ovules of a high frequency female-sterile line (fsv1) and a rice wild-type line (Gui 99) during ovule development. As a result, 522 known miRNAs and 295 novel miRNAs were expressed in the developing ovule of rice, while 100 known miRNAs were significantly differentially expressed between these two rice lines during ovule development. Combining with gene expression information, a total of 627 coherent target genes of these differential expressed known miRNAs between fsv1 and Gui 99 were identified. The functional analyses of these coherent target genes revealed that the coherent target genes of differential expressed known miRNAs between the two rice lines are involved in many biological pathways, such as protein degradation, auxin signal transduction, and transcription factor regulation. These results provide us with important clues to investigate the regulatory roles of miRNAs in rice female gametophyte abortion.

MicroRNAs in crop improvement: fine-tuners for complex traits

One of the most common challenges for both conventional and modern crop improvement is that the appearance of one desirable trait in a new crop variety is always balanced by the impairment of one or more other beneficial characteristics. The best way to overcome this problem is the flexible utilization of regulatory genes, especially genes that provide more efficient and precise regulation in a targeted manner. MicroRNAs (miRNAs), a type of short non-coding RNA, are promising candidates in this area due to their role as master modulators of gene expression at the post-transcriptional level, targeting messenger RNAs for cleavage or directing translational inhibition in eukaryotes. We herein highlight the current understanding of the biological role of miRNAs in orchestrating distinct agriculturally important traits by summarizing recent functional analyses of 65 miRNAs in 9 major crops worldwide. The integration of current miRNA knowledge with conventional and modern crop improvement strategies is also discussed.

Analysis of small RNAs revealed differential expressions during pollen and embryo sac development in autotetraploid rice

Background

Partial pollen and embryo sac sterilities are the two main reasons for low fertility in autotetraploid rice. Our previous study revealed that small RNAs changes may associate with pollen fertility in autotetraploid rice. However, knowledge on comparative analysis between the development of pollen and embryo sac by small RNAs in autotetraploid rice is still unknown. In the present study, WE-CLSM (whole-mount eosin B-staining confocal laser scanning microscopy) and high-throughput sequencing technology was employed to examine the cytological variations and to analyze small RNAs changes during pollen and embryo sac development in autotetraploid rice compared with its diploid counterpart.

Results

A total of 321 and 368 differentially expressed miRNAs (DEM) were detected during pollen and embryo sac development in autotetraploid rice, respectively. Gene Ontology enrichment analysis on the targets of DEM associated with embryo sac and pollen development revealed 30 prominent functional gene classes, such as cell differentiation and signal transduction during embryo sac development, while only 7 prominent functional gene classes, such as flower development and transcription factor activity, were detected during pollen development in autotetraploid rice. The expression levels of 39 DEM, which revealed interaction with meiosis-related genes, showed opposite expression patterns during pollen and embryo sac development. Of these DEM, osa-miR1436_L + 3_1ss5CT and osa-miR167h-3p were associated with the female meiosis, while osa-miR159a.1 and osa-MIR159a-p5 were related with the male meiosis. 21 nt-phasiRNAs were detected during both pollen and embryo sac development, while 24 nt-phasiRNAs were found only in pollen development, which displayed down-regulation in autotetraploid compared to diploid rice and their spatial-temporal expression patterns were similar to osa-miR2275d. 24 nt TEs-siRNAs were found to be up-regulated in embryo sac but down-regulated in pollen development.

Conclusion

The above results not only provide the small RNAs changes during four landmark stages of pollen and embryo sac development in autotetraploid rice but also have identified specifically expressed miRNAs, especially meiosis-related miRNAs, pollen-specific-24 nt-phasiRNAs and TEs-siRNAs in autotetraploid rice. Together, these findings provide a foundation for understanding the effect of polyploidy on small RNAs expression patterns during pollen and embryo sac development that may lead to different abnormalities in autotetraploid rice.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3526-8) contains supplementary material, which is available to authorized users.

Characterization and Function of MicroRNA∗s in Plants

MicroRNAs, a group of non-coding RNA molecules, play essential roles in a wide range of cellular processes in different molecules, cells, and organisms. In plants, microRNAs are a class of 20- to 24-nucleotides endogenous small RNAs that repress gene expression. The microRNA guide strand (miRNA) and its complementary strand (miRNA^∗) both originate from the miRNA/miRNA^∗ duplex. Generally, the guide strands act as post-transcriptional regulators that suppress gene expression by cleaving their target mRNA transcripts, whereas the complementary strands were thought to be degraded as 'passenger strands.' However, the complementary strand has been confirmed to possess significant biological functionality in recent reports. In this review, we summarized the binding characteristics of the miRNA^∗ strands with ARGONAUTE proteins, their tissue-specific accumulations and their biological functions, illustrating the essential roles of miRNA^∗s in biological processes and therefore providing directions for further exploration.

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.

Small RNA and degradome profiling reveals miRNA regulation in the seed germination of ancient eudicot Nelumbo nucifera

BACKGROUND

MicroRNAs (miRNAs) play important roles in plant growth and development. MiRNAs and their targets have been widely studied in model plants, but limited knowledge is available concerning this small RNA population and their targets in sacred lotus (Nelumbo nucifera Gaertn.).

RESULTS

In this study, a total of 145 known miRNAs belonging to 47 families and 78 novel miRNAs were identified during seed germination using high-throughput small RNA sequencing. Furthermore, some miRNA families which have not yet been reported in monocot or eudicot species were detected in N. nucifera, indicating that these miRNAs was divergence from monocots and core eudicots during evolution. Using degradome sequencing, 2580 targets were detected for all the miRNAs. GO (Gene Ontology) and KEGG pathway analyses showed that many target genes enriched in "regulation of transcription" and involved in "carbohydrate", "amino acid and energy metabolism". Nine miRNAs and three corresponding targets of them were further validated by quantitative RT-PCR.

CONCLUSIONS

The results present here suggested that many miRNAs were involved in the regulation of seed germination of sacred lotus, providing a foundation for future studies of sacred lotus seed longevity. Comparative analysis of miRNAs from different plants also provided insight into the evolutionary gains and losses of miRNAs in plants.

Dual transcriptome analysis reveals insights into the response to Rice black-streaked dwarf virus in maize

Maize rough dwarf disease (MRDD) is a viral infection that results in heavy yield losses in maize worldwide, particularly in the summer maize-growing regions of China. MRDD is caused by the Rice black-streaked dwarf virus (RBSDV). In the present study, analyses of microRNAs (miRNAs), the degradome, and transcriptome sequences were used to elucidate the RBSDV-responsive pathway(s) in maize. Genomic analysis indicated that the expression of three non-conserved and 28 conserved miRNAs, representing 17 known miRNA families and 14 novel miRNAs, were significantly altered in response to RBSDV when maize was inoculated at the V3 (third leaf) stage. A total of 99 target transcripts from 48 genes of 10 known miRNAs were found to be responsive to RBSDV infection. The annotations of these target genes include a SQUAMOSA promoter binding (SPB) protein, a P450 reductase, an oxidoreductase, and a ubiquitin-related gene, among others. Characterization of the entire transcriptome suggested that a total of 28 and 1085 differentially expressed genes (DEGs) were detected at 1.5 and 3.0 d, respectively, after artificial inoculation with RBSDV. The expression patterns of cell wall- and chloroplast-related genes, and disease resistance- and stress-related genes changed significantly in response to RBSDV infection. The negatively regulated genes GRMZM2G069316 and GRMZM2G031169, which are the target genes for miR169i-p5 and miR8155, were identified as a nucleolin and a NAD(P)-binding Rossmann-fold superfamily protein in maize, respectively. The gene ontology term GO:0003824, including GRMZM2G031169 and other 51 DEGs, was designated as responsive to RBSDV.

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.

Small RNA Profiles of the Rice PTGMS Line Wuxiang S Reveal miRNAs Involved in Fertility Transition

MicroRNAs (miRNAs) play key roles in the regulation of plant growth and developmental processes. In this study, RNA-seq was used to examine the expression profiles of miRNAs in a novel, photo-thermosensitive genic male sterile (PTGMS) rice line, Wuxiang S (WXS), during fertility transition. A total of 497 known miRNAs and 273 novel miRNAs were identified. In a differential expression analysis, 26 miRNAs exhibited significant differential expression between WXS (Sterile, S) and WXS (Fertile, F). Some of these miRNAs were validated by quantitative real-time PCR. Among these miRNAs, 11 showed decreased expression levels, and 15 showed increased expression levels in WXS (S) compared to WXS (F). Some of these miRNAs, such as osa-miR156a-j, osa-miR164d, and osa-miR528, were shown to be negatively correlated with their targets. These targets have previously been reported to be related to pollen development and male sterility, suggesting that these miRNAs may be involved in the regulation of pollen development in the rice PTGMS line WXS. Furthermore, miRNA-mediated editing events were also observed. In this study, a possible model for the control of signaling pathways during the process of fertility transition in the rice PTGMS line WXS by miRNAs was developed. These findings contribute to our understanding of the roles of miRNAs during anther development in PTGMS lines in rice.

Comparative Small RNA Analysis of Pollen Development in Autotetraploid and Diploid Rice

MicroRNAs (miRNAs) play key roles in plant reproduction. However, knowledge on microRNAome analysis in autotetraploid rice is rather limited. Here, high-throughput sequencing technology was employed to analyze miRNAomes during pollen development in diploid and polyploid rice. A total of 172 differentially expressed miRNAs (DEM) were detected in autotetraploid rice compared to its diploid counterpart, and 57 miRNAs were specifically expressed in autotetraploid rice. Of the 172 DEM, 115 and 61 miRNAs exhibited up- and down-regulation, respectively. Gene Ontology analysis on the targets of up-regulated DEM showed that they were enriched in transport and membrane in pre-meiotic interphase, reproduction in meiosis, and nucleotide binding in single microspore stage. osa-miR5788 and osa-miR1432-5p_R+1 were up-regulated in meiosis and their targets revealed interaction with the meiosis-related genes, suggesting that they may involve in the genes regulation associated with the chromosome behavior. Abundant 24 nt siRNAs associated with transposable elements were found in autotetraploid rice during pollen development; however, they significantly declined in diploid rice, suggesting that 24 nt siRNAs may play a role in pollen development. These findings provide a foundation for understanding the effect of polyploidy on small RNA expression patterns during pollen development that cause pollen sterility in autotetraploid rice.

Comparative analysis of microRNA profiles of rice anthers between cool-sensitive and cool-tolerant cultivars under cool-temperature stress

Plants subjected to abiotic stress can regulate gene expression post-transcriptionally by means of small RNAs such as microRNAs. Cool-temperature stress causes abnormal tapetum hypertrophy in rice anthers, leading to pollen sterility. As a first step toward understanding the molecular mechanisms of cool tolerance in developing anthers of rice, we report here a comprehensive comparative analysis of microRNAs between cool-sensitive Sasanishiki and cool-tolerant Hitomebore cultivars. High-throughput Illumina sequencing revealed 241 known and 46 novel microRNAs. Interestingly, 15 of these microRNAs accumulated differentially in the two cultivars at the uninucleate microspore stage under cool conditions. Inverse correlations between expression patterns of microRNAs and their target genes were confirmed by quantitative RT-PCR analysis, and cleavage sites of some of the target genes were determined by 5' RNA ligase-mediated RACE experiments. Thus, our data are useful resources to elucidate microRNA-mediated mechanism(s) of cool tolerance in rice anthers at the booting stage.

Identification of microRNAs and Their Target Genes Explores miRNA-Mediated Regulatory Network of Cytoplasmic Male Sterility Occurrence during Anther Development in Radish (Raphanus sativus L.)

MicroRNAs (miRNAs) are a type of endogenous non-coding small RNAs that play critical roles in plant growth and developmental processes. Cytoplasmic male sterility (CMS) is typically a maternally inherited trait and widely used in plant heterosis utilization. However, the miRNA-mediated regulatory network of CMS occurrence during anther development remains largely unknown in radish. In this study, a comparative small RNAome sequencing was conducted in floral buds of CMS line 'WA' and its maintainer line 'WB' by high-throughput sequencing. A total of 162 known miRNAs belonging to 25 conserved and 24 non-conserved miRNA families were isolated and 27 potential novel miRNA families were identified for the first time in floral buds of radish. Of these miRNAs, 28 known and 14 potential novel miRNAs were differentially expressed during anther development. Several target genes for CMS occurrence-related miRNAs encode important transcription factors and functional proteins, which might be involved in multiple biological processes including auxin signaling pathways, signal transduction, miRNA target silencing, floral organ development, and organellar gene expression. Moreover, the expression patterns of several CMS occurrence-related miRNAs and their targets during three stages of anther development were validated by qRT-PCR. In addition, a potential miRNA-mediated regulatory network of CMS occurrence during anther development was firstly proposed in radish. These findings could contribute new insights into complex miRNA-mediated genetic regulatory network of CMS occurrence and advance our understanding of the roles of miRNAs during CMS occurrence and microspore formation in radish and other crops.

The expanding world of small RNAs in plants

Plant genomes encode various small RNAs that function in distinct, yet overlapping, genetic and epigenetic silencing pathways. However, the abundance and diversity of small-RNA classes varies among plant species, suggesting coevolution between environmental adaptations and gene-silencing mechanisms. Biogenesis of small RNAs in plants is well understood, but we are just beginning to uncover their intricate regulation and activity. Here, we discuss the biogenesis of plant small RNAs, such as microRNAs, secondary siRNAs and heterochromatic siRNAs, and their diverse cellular and developmental functions, including in reproductive transitions, genomic imprinting and paramutation. We also discuss the diversification of small-RNA-directed silencing pathways through the expansion of RNA-dependent RNA polymerases, DICER proteins and ARGONAUTE proteins.

High-throughput sequencing reveals differential expression of miRNAs in tomato inoculated with Phytophthora infestans

The characterization and compare expression profiling of the miRNA transcriptome lay a solid foundation for unraveling the complex miRNA-mediated regulatory network in tomato resistance mechanisms against LB. MicroRNAs (miRNAs) are a class of small endogenous non-coding RNAs with 20-24 nt. They have been identified in many plants with their diverse regulatory roles in biotic stresses. The knowledge, that miRNAs regulate late blight (LB), caused by Phytophthora infestans, is rather limited. In this study, we used miRNA-Seq to investigate the miRNA expression difference between the tomatoes treated with and without P. infestans. A total of 42,714,516 raw reads were generated from two small RNA libraries by high-throughput sequencing. Finally, 207 known miRNAs and 67 new miRNAs were obtained. The differential expression profile of miRNAs in tomato was further analyzed with twofold change (P value ≤0.01). A total of 70 miRNAs were manifested to change significantly in samples treated with P. infestans, including 50 down-regulated miRNAs and 20 up-regulated miRNAs. Moreover, a total of 73 target genes were acquired for 28 differentially expressed miRNAs by psRNATarget analysis. By enrichment pathway analysis of target genes, plant-pathogen interaction was the most highly relevant pathway which played an important role in disease defense. In addition, 30 miRNAs were selected for qRT-PCR to validate their expression patterns. The expression patterns for targets of miR6027, miR5300, miR476b, miR159a, miR164a and miRn13 were selectively examined, and the results showed that there was a negative correlation on the expression patterns between miRNAs and their targets. The targets have previously been reported to be related with plant immune and involved in plant-pathogen interaction pathway in this study, suggesting these miRNAs might act as regulators in process of tomato resistance against P. infestans. These discoveries will provide us useful information to explain tomato resistance mechanisms against LB.