Over-expression of miR158 causes pollen abortion in Brassica campestris ssp. chinensis

Artificial miRNAs and target-mimics as potential tools for crop improvement

MicroRNAs (miRNAs) are endogenous, small molecules that negatively regulate gene expression to control the normal development and stress response in plants. They mediate epigenetic changes and regulate gene expression at both transcriptional and post-transcriptional levels. Synthetic biology approaches have been utilized to design efficient artificial miRNAs (amiRNAs) or target-mimics to regulate specific gene expression for understanding the biological function of genes and crop improvement. The amiRNA based gene silencing is an effective technique to "turn off" gene expression, while miRNA target-mimics or decoys are used for efficiently down regulating miRNAs and "turn on" gene expression. In this context, the development of endogenous target-mimics (eTMs) and short tandem target mimics (STTMs) represent promising biotechnological tools for enhancing crop traits like stress tolerance and disease resistance. Through this review, we present the recent developments in understanding plant miRNA biogenesis, which is utilized for the efficient design and development of amiRNAs. This is important to incorporate the artificially synthesized miRNAs as internal components and utilizing miRNA biogenesis pathways for the programming of synthetic circuits to improve crop tolerance to various abiotic and biotic stress factors. The review also examines the recent developments in the use of miRNA target-mimics or decoys for efficiently down regulating miRNAs for trait improvement. A perspective analysis and challenges on the use of amiRNAs and STTM as potent tools to engineer useful traits in plants have also been presented.

miRNA Sequencing Analysis in Maize Roots Treated with Neutral and Alkaline Salts

Soil salinization/alkalization is a complex environmental factor that includes not only neutral salt NaCl but also other components like Na2CO3. miRNAs, as small molecules that regulate gene expression post-transcriptionally, are involved in plant responses to abiotic stress. In this study, maize seedling roots were treated for 5 h with 100 mM NaCl, 50 mM Na2CO3, and H2O, respectively. Sequencing analysis of differentially expressed miRNAs under these conditions revealed that the Na2CO3 treatment group had the most differentially expressed miRNAs. Cluster analysis indicated their main involvement in the regulation of ion transport, binding, metabolism, and phenylpropanoid and flavonoid biosynthesis pathways. The unique differentially expressed miRNAs in the NaCl treatment group were related to the sulfur metabolism pathway. This indicates a significant difference in the response patterns of maize to different treatment groups. This study provides theoretical evidence and genetic resources for further analysis of the molecular mechanisms behind maize's salt-alkali tolerance.

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.

Genomic survey of high-throughput RNA-Seq data implicates involvement of long intergenic non-coding RNAs (lincRNAs) in cytoplasmic male-sterility and fertility restoration in pigeon pea

BACKGROUND

Long-intergenic non-coding RNAs (lincRNAs) originate from intergenic regions and have no coding potential. LincRNAs have emerged as key players in the regulation of various biological processes in plant development. Cytoplasmic male-sterility (CMS) in association with restorer-of-fertility (Rf) systems makes it a highly reliable tool for exploring heterosis for producing commercial hybrid seeds. To date, there have been no reports of lincRNAs during pollen development in CMS and fertility restorer lines in pigeon pea.

OBJECTIVE

Identification of lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines.

METHODS

We employed a computational approach to identify lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines using RNA-Seq data.

RESULTS

We predicted a total of 2145 potential lincRNAs of which 966 were observed to be differentially expressed between the sterile and fertile pollen. We identified, 927 cis-regulated and 383 trans-regulated target genes of the lincRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the target genes revealed that these genes were specifically enriched in pathways like pollen and pollen tube development, oxidative phosphorylation, etc. We detected 23 lincRNAs that were co-expressed with 17 pollen-related genes with known functions. Fifty-nine lincRNAs were predicted to be endogenous target mimics (eTMs) for 25 miRNAs, and found to be associated with pollen development. The, lincRNA regulatory networks revealed that different lincRNA-miRNA-mRNA networks might be associated with CMS and fertility restoration.

CONCLUSION

Thus, this study provides valuable information by highlighting the functions of lincRNAs as regulators during pollen development in pigeon pea and utilization in hybrid seed production.

Indian Himalayan natural Arabidopsis thaliana accessions with abolished miR158 levels exhibit robust miR173-initiated trans-acting cascade silencing

Small RNAs (sRNAs) such as microRNAs (miRNAs) and small interfering RNAs (siRNAs) are short 20-24-nucleotide non-coding RNAs. They are key regulators of gene expression in plants and other organisms. Several 22-nucleotide miRNAs trigger biogenesis cascades of trans-acting secondary siRNAs, which are involved in various developmental and stress responses. Here we show that Himalayan Arabidopsis thaliana accessions having natural mutations in the miR158 locus exhibit robust cascade silencing of the pentatricopeptide repeat (PPR)-like locus. Furthermore, we show that these cascade sRNAs trigger tertiary silencing of a gene involved in transpiration and stomatal opening. The natural deletions or insertions in MIR158 led to improper processing of miR158 precursors, thereby blocking synthesis of mature miR158. Reduced miR158 levels led to increased levels of its target, a pseudo-PPR gene that is targeted by tasiRNAs generated by the miR173 cascade in other accessions. Using sRNA datasets derived from Indian Himalayan accessions, as well as overexpression and knockout lines of miR158, we show that absence of miR158 led to buildup of pseudo-PPR-derived tertiary sRNAs. These tertiary sRNAs mediated robust silencing of a gene involved in stomatal closure in Himalayan accessions lacking miR158 expression. We functionally validated the tertiary phasiRNA that targets NHX2, which encodes a Na⁺ -K⁺ /H⁺ antiporter protein, thereby regulating transpiration and stomatal conductance. Overall, we report the role of the miRNA-TAS-siRNA-pseudogene-tertiary phasiRNA-NHX2 pathway in plant adaptation.

Identification of microRNA158 from Anthurium andraeanum and Its Function in Cold Stress Tolerance

Anthurium andraeanum is a tropical flower with high ornamental and economic value. Cold stress is one of the major abiotic stresses affecting the quality and value of A. andraeanum; thus, improving the cold tolerance of this species is an important breeding objective. MicroRNAs (miRNAs) have a critical role in plant abiotic stress responses, but their specific molecular regulatory mechanisms are largely unknown, including those related to the cold stress response in A. andraeanum. Here, we identified and cloned the precursor of miR158 from A. andraeanum (Aa-miR158). Both Aa-miR158 and its target gene (c48247) had higher expression levels in strong leaves than in other tissues or organs. Further study revealed that the transcript level of Aa-miR158 was increased by cold stress. Heterologous overexpression of Aa-miR158 improved cold stress tolerance in Arabidopsis, which was associated with decreases in the malondialdehyde (MDA) concentration and relative electrical conductivity (REC) as well as increases in peroxidase (POD) and catalase (CAT) activity. Moreover, overexpressing Aa-miR158 significantly increased the expression of endogenous genes related to cold stress tolerance and reactive oxygen species (ROS) levels in transgenic Arabidopsis under cold stress. Overall, our results demonstrate that Aa-miR158 is significantly involved in the cold stress response and provide a new strategy for cold tolerance breeding of A. andraeanum.

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

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

Identification of miRNAs and their target genes in genic male sterility lines in Brassica napus by small RNA sequencing

BACKGROUND

Brassica napus is the third leading source of edible oil in the world. Genic male sterility (GMS) lines provide crucial material for harnessing heterosis for rapeseed. GMS lines have been used successfully for rapeseed hybrid production in China. MicroRNAs (miRNAs) play crucial regulatory roles in various plant growth, development, and stress response processes. However, reports on miRNAs that regulate the pollen development of GMS lines in B. napus are few.

RESULTS

In this study, 12 small RNA and transcriptome libraries were constructed and sequenced for the flower buds from the fertile and sterile lines of two recessive GMS (RGMS) lines, namely, "6251AB" and "6284AB". At the same time, 12 small RNA and transcriptome libraries were also constructed and sequenced for the flower buds from the fertile and sterile lines of two dominant GMS (DGMS) lines, namely, "4001AB" and "4006AB". Based on the results, 46 known miRNAs, 27 novel miRNAs on the other arm of known pre-miRNAs, and 44 new conserved miRNAs were identified. Thirty-five pairs of novel miRNA-3p/miRNA-5p were found. Among all the identified miRNAs, fifteen differentially expressed miRNAs with over 1.5-fold change between flower buds of sterile and fertile lines were identified, including six differentially expressed miRNAs between "4001A" and "4001B", two differentially expressed miRNAs between "4006A" and "4006B", four differentially expressed miRNAs between "6251A" and "6251B", and ten differentially expressed miRNAs between "6284A" and "6284B". The correlation analysis of small RNA and transcriptome sequencing was conducted. And 257 candidate target genes were predicted for the 15 differentially expressed miRNAs. The results of 5' modified RACE indicated that BnaA09g48720D, BnaA09g11120D, and BnaCnng51960D were cleaved by bna-miR398a-3p, bna-miR158-3p and bna-miR159a, respectively. Among the differentially expressed miRNAs, miR159 was chosen to analyze its function. Overexpression of bna-miR159 in Arabidopsis resulted in decreased seed setting rate, and shortened siliques, illustrating that miR159 may regulate the fertility and silique development in rapeseed.

CONCLUSIONS

Our findings provide an overview of miRNAs that are potentially involved in GMS and pollen development. New information on miRNAs and their related target genes are provided to exploit the GMS mechanism and reveal the miRNA networks in B. napus.

Global Analysis of the Genetic Variations in miRNA-Targeted Sites and Their Correlations With Agronomic Traits in Rapeseed

MicroRNAs (miRNAs) and their target genes play vital roles in crops. However, the genetic variations in miRNA-targeted sites that affect miRNA cleavage efficiency and their correlations with agronomic traits in crops remain unexplored. On the basis of a genome-wide DNA re-sequencing of 210 elite rapeseed (Brassica napus) accessions, we identified the single nucleotide polymorphisms (SNPs) and insertions/deletions (INDELs) in miRNA-targeted sites complementary to miRNAs. Variant calling revealed 7.14 million SNPs and 2.89 million INDELs throughout the genomes of 210 rapeseed accessions. Furthermore, we detected 330 SNPs and 79 INDELs in 357 miRNA target sites, of which 33.50% were rare variants. We also analyzed the correlation between the genetic variations in miRNA target sites and 12 rapeseed agronomic traits. Eleven SNPs in miRNA target sites were significantly correlated with phenotypes in three consecutive years. More specifically, three correlated SNPs within the miRNA-binding regions of BnSPL9-3, BnSPL13-2, and BnCUC1-2 were in the loci associated with the branch angle, seed weight, and silique number, respectively; expression profiling suggested that the variation at these 3 miRNA target sites significantly affected the expression level of the corresponding target genes. Taken together, the results of this study provide researchers and breeders with a global view of the genetic variations in miRNA-targeted sites in rapeseed and reveal the potential effects of these genetic variations on elite agronomic traits.

Circular RNAs Repertoire and Expression Profile during Brassica rapa Pollen Development

Circular RNAs (circRNAs) are covalently closed RNA molecules generated by the back-splicing of exons from linear precursor mRNAs. Though various linear RNAs have been shown to play important regulatory roles in many biological and developmental processes, little is known about the role of their circular counterparts. In this study, we performed high-throughput RNA sequencing to delineate the expression profile and potential function of circRNAs during the five stages of pollen development in Brassica rapa. A total of 1180 circRNAs were detected in pollen development, of which 367 showed stage-specific expression patterns. Functional enrichment and metabolic pathway analysis showed that the parent genes of circRNAs were mainly involved in pollen-related molecular and biological processes such as mitotic and meiotic cell division, DNA processes, protein synthesis, protein modification, and polysaccharide biosynthesis. Moreover, by predicting the circRNA–miRNA network from our differentially expressed circRNAs, we found 88 circRNAs with potential miRNA binding sites, suggesting their role in post-transcriptional regulation of the genes. Finally, we confirmed the back-splicing sites of nine selected circRNAs using divergent primers and Sanger sequencing. Our study presents the systematic analysis of circular RNAs during pollen development and forms the basis of future studies for unlocking complex gene regulatory networks underpinning reproduction in flowering plants.

Utilization of microRNAs and their regulatory functions for improving biotic stress tolerance in tea plant [Camellia sinensis (L.) O. Kuntze]

MicroRNAs play a central role in responses to biotic stressors through their interactions with their target mRNAs. Tea plant (Camellia sinensis L.), an important beverage crop, is vulnerable to tea geometrid and anthracnose disease that causes considerable crop loss and tea production worldwide. Sustainable production of tea in the current scenario to biotic factors is major challenges. To overcome the problem of biotic stresses, high-throughput sequencing (HTS) with bioinformatics analyses has been used as an effective approach for the identification of stress-responsive miRNAs and their regulatory functions in tea plant. These stress-responsive miRNAs can be utilized for miRNA-mediated gene silencing to enhance stress tolerance in tea plant. Therefore, this review summarizes the current understanding of miRNAs regulatory functions in tea plant responding to Ectropis oblique and Colletotrichum gloeosporioides attacks for future miRNA research. Also, it highlights the utilization of miRNA-mediated gene silencing strategies for developing biotic stress-tolerant tea plant.

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.

Joint RNA-Seq and miRNA Profiling Analyses to Reveal Molecular Mechanisms in Regulating Thickness of Pod Canopy in Brassica napus

Oilseed rape (Brassica napus) is the second largest oilseed crop worldwide. As an architecture component of B. napus, thickness of pod canopy (TPC) plays an important role in yield formation, especially under high-density cultivation conditions. However, the mechanisms underlying the regulation of TPC remain unclear. RNA and microRNA (miRNA) profiling of two groups of B. napus lines with significantly different TPC at the bolting with a tiny bud stage revealed differential expressions of numerous genes involved in nitrogen-related pathways. Expression of several nitrogen-related response genes, including ASP5, ASP2, ASN3, ATCYSC1, PAL2, APT2, CRTISO, and COX15, was dramatically changed in the thick TPC lines compared to those in the thin TPC lines. Differentially expressed miRNAs also included many involved in nitrogen-related pathways. Expression of most target genes was negatively associated with corresponding miRNAs, such as miR159, miR6029, and miR827. In addition, 12 (including miR319, miR845, and miR158) differentially expressed miRNAs between two plant tissues sampled (stem apex and flower bud) were identified, implying that they might have roles in determining overall plant architecture. These results suggest that nitrogen signaling may play a pivotal role in regulating TPC in B. napus.

MicroRNAs and new biotechnological tools for its modulation and improving stress tolerance in plants

MicroRNAs (miRNAs) modulate the abundance and spatial-temporal accumulation of target mRNAs and indirectly regulate several plant processes. Transcriptional regulation of the genes encoding miRNAs (MIR genes) can be activated by numerous transcription factors, which themselves are regulated by other miRNAs. Fine-tuning of MIR genes or miRNAs is a powerful biotechnological strategy to improve tolerance to abiotic or biotic stresses in crops of economic importance. Current approaches for miRNA fine-tuning are based on the down- or up-regulation of MIR gene transcription and the use of genetic engineering tools to manipulate the final concentration of these miRNAs in the cytoplasm. Transgenesis, cisgenesis, intragenesis, artificial MIR genes, endogenous and artificial target mimicry, MIR genes editing using Meganucleases, ZNF proteins, TALENs and CRISPR/Cas9 or CRISPR/Cpf1, CRISPR/dCas9 or dCpf1, CRISPR13a, topical delivery of miRNAs and epigenetic memory have been successfully explored to MIR gene or miRNA modulation and improve agronomic traits in several model or crop plants. However, advantages and drawbacks of each of these new biotechnological tools (NBTs) are still not well understood. In this review, we provide a brief overview of the biogenesis and role of miRNAs in response to abiotic or biotic stresses, we present critically the main NBTs used for the manipulation of MIR genes and miRNAs, we show current efforts and findings with the MIR genes and miRNAs modulation in plants, and we summarize the advantages and drawbacks of these NBTs and provide some alternatives to overcome. Finally, challenges and future perspectives to miRNA modulating in important crops are also discussed.

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.

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.

Small RNA Sequencing Reveals Differential miRNA Expression in the Early Development of Broccoli (Brassica oleracea var. italica) Pollen

Pollen development is an important and complex biological process in the sexual reproduction of flowering plants. Although the cytological characteristics of pollen development are well defined, the regulation of its early stages remains largely unknown. In the present study, miRNAs were explored in the early development of broccoli (Brassica oleracea var. italica) pollen. A total of 333 known miRNAs that originated from 235 miRNA families were detected. Fifty-five novel miRNA candidates were identified. Sixty of the 333 known miRNAs and 49 of the 55 predicted novel miRNAs exhibited significantly differential expression profiling in the three distinct developmental stages of broccoli pollen. Among these differentially expressed miRNAs, miRNAs that would be involved in the developmental phase transition from uninucleate microspores to binucleate pollen grains or from binucleate to trinucleate pollen grains were identified. miRNAs that showed significantly enriched expression in a specific early stage of broccoli pollen development were also observed. In addition, 552 targets for 127 known miRNAs and 69 targets for 40 predicted novel miRNAs were bioinformatically identified. Functional annotation and GO (Gene Ontology) analysis indicated that the putative miRNA targets showed significant enrichment in GO terms that were related to plant organ formation and morphogenesis. Some of enriched GO terms were detected for the targets directly involved in plant male reproduction development. These findings provided new insights into the functions of miRNA-mediated regulatory networks in broccoli pollen development.