Identification of miRNAs and their targets from Brassica napus by high-throughput sequencing and degradome analysis

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.

Clone and Function Verification of the OPR gene in Brassica napus Related to Linoleic Acid Synthesis

BACKGROUND

Fatty acid composition and content affect rapeseed oil quality. Fatty acid synthesis-related genes in rapeseed have been studied globally by researchers. Nevertheless, rapeseed oil is mainly composed of seven different fatty acids (FA), and each fatty acid was regulated by different genes. Furthermore, different FA affect each other, which needs continuous and in-depth research to obtain more clear results in Brassica napus.

RESULTS

In this paper, broad-scale miRNA expression profiles were constructed and 21 differentially expressed miRNAs were detected. GO enrichment analysis showed that most up-regulated proteins were involved in transcription factor activity and catalytic activity. KEGG pathway enrichment analysis indicated that 20 pathways involving 36 target genes were enriched, of which the bna00592 pathway may be involved in fatty acid metabolism. The results were verified using a quantitative real-time PCR (RT-qPCR) analysis, we found that the target gene of bna-miR156b > c > g was the OPR (12-oxo-phytodienoic acid reductase). Four copies of OPR gene were found, and the over-expression vectors (pCAMBIA1300-35 s-OPR and pCAMBIA1300-RNAi-OPR) were constructed to verify their functions. In T1 and T2 generation, the content of linoleic acid (LA) increased significantly in OE but deceased in OPRi.

CONCLUSIONS

This is the first study to provide four copies of the OPR gene that regulates LA metabolism, can be used for the molecular mechanism of LA and optimizing fatty acid profiles in oilseed for breeding programs.

The miR169n-NF-YA8 regulation module involved in drought resistance in Brassica napus L

As an ancient and conserved plant microRNA (miRNA) family, miR169 targets nuclear factor Y subunit alpha (NF-YA) family members. The miR169-NF-YA module is associated with plant development and various stress responses. However, the function of miR169 in response to drought stress in rapeseed (Brassica napus L.) is unclear. In the present study, we showed that miR169n acted as a negative regulator of drought resistance in rapeseed by targeting a nuclear factor Y-A gene, NF-YA8. miR169n was strongly down-regulated by drought stress. Expression of a miR169n target mimicry construct (MIM169n) which functioned as a sponge to trap miR169n resulted in enhanced resistance of transgenic plants to both osmotic stress at the post-germination stage and drought stress at the seedling stage. MIM169n plants had a higher relative water content (RWC) and proline content, lower relative electrolyte leakage (REL), and showed higher antioxidative capability compared with those of control (CK) plants under drought stress. Moreover, NF-YA8 was verified as a target of miR169n, and overexpression of NF-YA8 led to improved tolerance of rapeseed to osmotic stress at the post-germination stage. Overall, our findings implied that the miR169n-NF-YA8 regulatory module could serve as a potential target for genetic improvement of drought resistance in B. napus.

Genome-wide miRNA analysis and integrated network for flavonoid biosynthesis in Osmanthus fragrans

Background

Osmanthus fragrans is an important economical plant containing multiple secondary metabolites including flavonoids and anthocyanins. During the past years, the roles of miRNAs in regulating the biosynthesis of secondary metabolites in plants have been widely investigated. However, few studies on miRNA expression profiles and the potential roles in regulating flavonoid biosynthesis have been reported in O. fragrans.

Results

In this study, we used high-throughput sequencing technology to analyze the expression profiles of miRNAs in leaf and flower tissues of O. fragrans. As a result, 106 conserved miRNAs distributed in 47 families and 88 novel miRNAs were identified. Further analysis showed there were 133 miRNAs differentially expressed in leaves and flowers. Additionally, the potential target genes of miRNAs as well as the related metabolic pathways were predicted. In the end, flavonoid content was measured in flower and leaf tissues and potential role of miR858 in regulating flavonoid synthesis was illustrated in O. fragrans.

Conclusions

This study not only provided the genome-wide miRNA profiles in the flower and leaf tissue of O. fragrans, but also investigated the potential regulatory role of miR858a in flavonoid synthesis in O. fragrans. The results specifically indicated the connection of miRNAs to the regulation of secondary metabolite biosynthesis in non-model economical plant.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07439-y.

Turnip mosaic virus in oilseed rape activates networks of sRNA-mediated interactions between viral and host genomes

Virus-induced plant diseases in cultivated plants cause important damages in yield. Although the mechanisms of virus infection are intensely studied at the cell biology level, only little is known about the molecular dialog between the invading virus and the host genome. Here we describe a combinatorial genome-wide approach to identify networks of sRNAs-guided post-transcriptional regulation within local Turnip mosaic virus (TuMV) infection sites in Brassica napus leaves. We show that the induction of host-encoded, virus-activated small interfering RNAs (vasiRNAs) observed in virus-infected tissues is accompanied by site-specific cleavage events on both viral and host RNAs that recalls the activity of small RNA-induced silencing complexes (RISC). Cleavage events also involve virus-derived siRNA (vsiRNA)–directed cleavage of target host transcripts as well as cleavage of viral RNA by both host vasiRNAs and vsiRNAs. Furthermore, certain coding genes act as virus-activated regulatory hubs to produce vasiRNAs for the targeting of other host genes. The observations draw an advanced model of plant-virus interactions and provide insights into the complex regulatory networking at the plant-virus interface within cells undergoing early stages of infection.

Pitzalis et al. use replicative RNAseq, small RNA (sRNA)seq, and parallel analysis of RNA ends (PARE)seq analysis to identify networks of sRNAs-guided post-transcriptional regulation within local Turnip mosaic virus infection sites. This study provides insights into the complex regulatory networking at the plantvirus interface within cells undergoing early stages of infection.

Genome-Wide Differential DNA Methylation and miRNA Expression Profiling Reveals Epigenetic Regulatory Mechanisms Underlying Nitrogen-Limitation-Triggered Adaptation and Use Efficiency Enhancement in Allotetraploid Rapeseed

Improving crop nitrogen (N) limitation adaptation (NLA) is a core approach to enhance N use efficiency (NUE) and reduce N fertilizer application. Rapeseed has a high demand for N nutrients for optimal plant growth and seed production, but it exhibits low NUE. Epigenetic modification, such as DNA methylation and modification from small RNAs, is key to plant adaptive responses to various stresses. However, epigenetic regulatory mechanisms underlying NLA and NUE remain elusive in allotetraploid B. napus. In this study, we identified overaccumulated carbohydrate, and improved primary and lateral roots in rapeseed plants under N limitation, which resulted in decreased plant nitrate concentrations, enhanced root-to-shoot N translocation, and increased NUE. Transcriptomics and RT-qPCR assays revealed that N limitation induced the expression of NRT1.1, NRT1.5, NRT1.7, NRT2.1/NAR2.1, and Gln1;1, and repressed the transcriptional levels of CLCa, NRT1.8, and NIA1. High-resolution whole genome bisulfite sequencing characterized 5094 differentially methylated genes involving ubiquitin-mediated proteolysis, N recycling, and phytohormone metabolism under N limitation. Hypermethylation/hypomethylation in promoter regions or gene bodies of some key N-metabolism genes might be involved in their transcriptional regulation by N limitation. Genome-wide miRNA sequencing identified 224 N limitation-responsive differentially expressed miRNAs regulating leaf development, amino acid metabolism, and plant hormone signal transduction. Furthermore, degradome sequencing and RT-qPCR assays revealed the miR827-NLA pathway regulating limited N-induced leaf senescence as well as the miR171-SCL6 and miR160-ARF17 pathways regulating root growth under N deficiency. Our study provides a comprehensive insight into the epigenetic regulatory mechanisms underlying rapeseed NLA, and it will be helpful for genetic engineering of NUE in crop species through epigenetic modification of some N metabolism-associated genes.

Rolling Circle Amplification (RCA)-Mediated Genome-Wide ihpRNAi Mutant Library Construction in Brassica napus

With the successful completion of genomic sequencing for Brassica napus, identification of novel genes, determination of functions performed by genes, and exploring the molecular mechanisms underlying important agronomic traits were challenged. Mutagenesis-based functional genomics techniques including chemical, physical, and insertional mutagenesis have been used successfully in the functional characterization of genes. However, these techniques had their disadvantages and inherent limitations for allopolyploid Brassica napus, which contained a large number of homologous and redundant genes. Long intron-spliced hairpin RNA (ihpRNA) constructs which contained inverted repeats of the target gene separated by an intron, had been shown to be very effective in triggering RNAi in plants. In the present study, the genome-wide long ihpRNA library of B. napus was constructed with the rolling circle amplification (RCA)-mediated technology. Using the phytoene desaturase (PDS) gene as a target control, it was shown that the RCA-mediated long ihpRNA construct was significantly effective in triggering gene silence in B. napus. Subsequently, the resultant long ihpRNA library was transformed into B. napus to produce corresponding RNAi mutants. Among the obtained transgenic ihpRNA population of B. napus, five ihpRNA lines with observable mutant phenotypes were acquired including alterations in the floral model and the stamen development. The target genes could be quickly identified using specific primers. These results showed that the RCA-mediated ihpRNA construction method was effective for the genome-wide long ihpRNA library of B. napus, therefore providing a platform for study of functional genomics in allopolyploid B. napus.

Identification and expression analysis of miRNAs and elucidation of their role in salt tolerance in rice varieties susceptible and tolerant to salinity

Soil salinization is a serious problem for cultivation of rice, as among cereals rice is the most salt sensitive crop, and more than 40% of the total agricultural land amounting to approximately 80 million ha the world over is salt affected. Salinity affects a plant in a varieties of ways, including ion toxicity, osmotic stress and oxidative damage. Since miRNAs occupy the top place in biochemical events determining a trait, understanding their role in salt tolerance is highly desirable, which may allow introduction of the trait in the rice cultivars of choice through biotechnological interventions. High throughput sequencing of sRNAs in the root and shoot tissues of the seedlings of the control and NaCl treated Pokkali, a salt-tolerant rice variety, identified 75 conserved miRNAs and mapped 200 sRNAs to the rice genome as novel miRNAs. Expression of nine novel miRNAs and two conserved miRNAs were confirmed by Northern blotting. Several of both conserved and novel miRNAs that expressed differentially in root and/or shoot tissues targeted transcription factors like AP2/EREBP domain protein, ARF, NAC, MYB, NF-YA, HD-Zip III, TCP and SBP reported to be involved in salt tolerance or in abiotic stress tolerance in general. Most of the novel miRNAs expressed in the salt tolerant wild rice Oryza coarctata, suggesting conservation of miRNAs in taxonomically related species. One of the novel miRNAs, osa-miR12477, also targeted L-ascorbate oxidase (LAO), indicating build-up of oxidative stress in the plant upon salt treatment, which was confirmed by DAB staining. Thus, salt tolerance might involve miRNA-mediated regulation of 1) cellular abundance of the hormone signaling components like EREBP and ARF, 2) synthesis of abiotic stress related transcription factors, and 3) antioxidative component like LAO for mitigation of oxidative damage. The study clearly indicated importance of osa-miR12477 regulated expression of LAO in salt tolerance in the plant.

Identification of conserved and novel miRNAs responsive to heat stress in flowering Chinese cabbage using high-throughput sequencing

Plant microRNAs (miRNAs) are noncoding and endogenous key regulators that play significant functions in regulating plant responses to stress, and plant growth and development. Heat stress is a critical abiotic stress that reduces the yield and quality of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). However, limited information is available on whether miRNAs are involved in the regulation of heat stress in B. campestris. A high-throughput sequencing approach was used to identify novel and conserved heat-responsive miRNAs in four small RNA libraries of flowering Chinese cabbage using leaves collected at 0 h, 1 h, 6 h and 12 h after a 38 °C heat-stress treatment. The analysis identified 41 conserved miRNAs (belonging to 19 MIR families), of which MIR156, MIR159, MIR168, MIR171 and MIR1885 had the most abundant molecules. Prediction and evaluation of novel miRNAs using the unannotated reads resulted in 18 candidate miRNAs. Differential expression analysis showed that most of the identified miRNAs were downregulated in heat-treated groups. To better understand functional importance, bioinformatic analysis predicted 432 unique putative target miRNAs involved in cells, cell parts, catalytic activity, cellular processes and abiotic stress responses. Furthermore, the Kyoto Encyclopedia of Genes and Genomes maps of flowering Chinese cabbage identified the significant role of miRNAs in stress adaptation and stress tolerance, and in several mitogen-activated protein kinases signaling pathways including cell death. This work presents a comprehensive study of the miRNAs for understanding the regulatory mechanisms and their participation in the heat stress of flowering Chinese cabbage.

Identification of Regulatory Networks of MicroRNAs and Their Targets in Response to Colletotrichum gloeosporioides in Tea Plant (Camellia sinensis L.)

Anthracnose disease is caused by Colletotrichum gloeosporioides, and is common in leaves of the tea plant (Camellia sinensis). MicroRNAs (miRNAs) have been known as key modulators of gene expression in response to environmental stresses, disease resistance, defense responses, and plant immunity. However, the role of miRNAs in responses to C. gloeosporioides remains unexplored in tea plant. Therefore, in the present study, six miRNA sequencing data sets and two degradome data sets were generated from C. gloeosporioides-inoculated and control tea leaves. A total of 485 conserved and 761 novel miRNAs were identified. Of those, 239 known and 369 novel miRNAs exhibited significantly differential expression under C. gloeosporioides stress. One thousand one hundred thirty-four and 596 mRNAs were identified as targets of 389 conserved and 299 novel miRNAs by degradome analysis, respectively. Based on degradome analysis, most of the predicted targets are negatively correlated with their corresponding conserved and novel miRNAs. The expression levels of 12 miRNAs and their targets were validated by quantitative real-time PCR. A negative correlation between expression profiles of five miRNAs (PC-5p-80764_22, csn-miR160c, csn-miR828a, csn-miR164a, and csn-miR169e) and their targets (WRKY, ARF, MYB75, NAC, and NFY transcription factor) was observed. The predicted targets of five interesting miRNAs were further validated through 5'RLM-RACE. Furthermore, Gene Ontology and metabolism pathway analysis revealed that most of the target genes were involved in the regulation of auxin pathway, ROS scavenging pathway, salicylic acid mediated pathway, receptor kinases, and transcription factors for plant growth and development as well as stress responses in tea plant against C. gloeosporioides stress. This study enriches the resources of stress-responsive miRNAs and their targets in C. sinensis and thus provides novel insights into the miRNA-mediated regulatory mechanisms, which could contribute to the enhanced susceptibility of C. gloeosporioides in tea plant.

Small RNA and degradome profiling involved in seed development and oil synthesis of Brassica napus

MicroRNAs (miRNAs) play a prominent role in post-transcriptional gene expression regulation and have been involved in various biological and metabolic processes to regulate gene expression. For Brassica napus, improving seed-weight and oil-content is the main breeding goal. In order to better understand the regulation mechanism of miRNAs during seed-weight formation and oil-content accumulation in B. napus, in this study, a high-throughput sequencing technology was used to profile miRNAs expression of Brassica napus immature seeds from one to six weeks after flowering. A total of 1,276 miRNAs, including 1,248 novel and 28 known miRNAs, were obtained from both the high-seed-weight with low-oil-content RNA pool (S03) and the low-seed-weight with high-oil-content RNA pool (S04). Analysis of their expression profiles disclosed that 300 novel and two known miRNAs were differentially expressed between S03 and S04. For degradome analysis, 57 genes with 64 degradation sites were predicted to be targeted for degradation by these miRNAs. Further bioinformatics analysis indicated that these differentially expressed miRNAs might participate in regulation of myriad cellular and molecular processes, during seed development and oil synthesis. Finally, 6 target genes with potential roles in regulation of seed development and 9 other targets in seed oil synthesis, were further confirmed as candidate genes from small RNA and degradome sequencing.

miRNAs play important roles in aroma weakening during the shelf life of 'Nanguo' pear after cold storage

Cold storage is commonly employed to delay senescence in 'Nanguo' pear after harvest. However, this technique also causes fruit aroma weakening. MicroRNAs are regulators of gene expression at the post-transcriptional level that play important roles in plant development and in eliciting responses to abiotic environmental stressors. In this study, the miRNA transcript profiles of the fruit on the first day (C0, LT0) move in and out of cold storage and the optimum tasting period (COTP, LTOTP) during shelf life at room temperature and after cold storage were analyzed, respectively. 314 known miRNAs were identified in 'Nanguo' pear; 176 and 135 miRNAs were significantly differentially expressed on the C0 vs. LT0 and on the COTP vs. LTOTP, respectively. After prediction the target genes of these differentially expressed miRNAs, 9 s-lipoxygenase (LOX2S), 13 s-lipoxygenase (LOX1_5), hydroperoxide lyase (HPL), and alcohol dehydrogenase (ADH1) were found differentially expressed, which were the key genes during aroma formation. The expression pattern of these target genes and the related miRNAs were identified by RT-PCR. mdm-miR172a-h, mdm-miR159a/b, mdm-miR160a-e, mdm-miR395a-i, mdm-miR399a/b/c, mdm/ppe-miR535a/b, and mdm-miR7120a/b may negatively regulate the target genes expression. These results indicate that miRNAs may play key roles in aroma weakening in cold storage 'Nanguo' pear and provide valuable information for studying the molecular mechanisms of miRNAs in the aroma weakening of fruit due to low temperature.

In silico identification and computational characterization of endogenous small interfering RNAs from diverse grapevine tissues and stages

Small interfering RNAs (siRNAs) are effectors of regulatory pathways underlying plant development, metabolism, and stress- and nutrient-signaling regulatory networks. The endogenous siRNAs are generally not conserved between plants; consequently, it is necessary and important to identify and characterize siRNAs from various plants. To address the nature and functions of siRNAs, and understand the biological roles of the huge siRNA population in grapevine (Vitis vinifera L.). The high-throughput sequencing technology was used to identify a large set of putative endogenous siRNAs from six grapevine tissues/organs. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed to classify the target genes of siRNA. In total, 520,519 candidate siRNAs were identified and their expression profiles exhibited typical temporal characters during grapevine development. In addition, we identified two grapevine trans-acting siRNA (TAS) gene homologs (VvTAS3 and VvTAS4) and the derived trans-acting siRNAs (tasiRNAs) that could target grapevine auxin response factor (ARF) and myeloblastosis (MYB) genes. Furthermore, the GO and KEGG analysis of target genes showed that most of them covered a broad range of functional categories, especially involving in disease-resistance process. The large-scale and completely genome-wide level identification and characterization of grapevine endogenous siRNAs from the diverse tissues by high throughput technology revealed the nature and functions of siRNAs in grapevine.

Integrated mRNA, sRNA, and degradome sequencing reveal oilseed rape complex responses to Sclerotinia sclerotiorum (Lib.) infection

Sclerotinia stem rot (SSR), caused by the fungal pathogen Sclerotinia sclerotiorum, is a devastating disease resulting in yield losses and decreases in seed quality in oilseed rape (Brassica napus) worldwide. However, the molecular mechanisms underlying the response of oilseed rape to S. sclerotiorum infection at the transcriptional and post-transcriptional levels are poorly understood. Here, we used an integrated omics approach (transcriptome, sRNAome, and degradome sequencing) on the Illumina platform to compare the RNA expression and post-transcriptional profiles of oilseed rape plants inoculated or not with S. sclerotiorum. In total, 7,065 differentially expressed genes (DEGs) compared with the mock-inoculated control at 48 hours post inoculation were identified. These DEGs were associated with protein kinases, signal transduction, transcription factors, hormones, pathogenesis-related proteins, secondary metabolism, and transport. In the sRNA-Seq analysis, 77 known and 176 novel miRNAs were identified; however, only 10 known and 41 novel miRNAs were differentially expressed between the samples inoculated or not with S. sclerotiorum. Degradome sequencing predicted 80 cleavage sites with 64 miRNAs. Integrated mRNA, sRNA and degradome sequencing analysis reveal oilseed rape complex responses to S. sclerotiorum infection. This study provides a global view of miRNA and mRNA expression profiles in oilseed rape following S. sclerotiorum infection.

Identification of miRNAs that regulate silique development in Brassica napus

MicroRNAs (miRNAs) are a class of non-coding small RNAs (sRNAs) that play crucial regulatory roles in various developmental processes. Silique length indirectly influences seed yield in rapeseed (Brassica napus); however, the molecular roles of miRNAs in silique length are largely unknown. Here, backcross progenies of rapeseed with long siliques (LS) and short siliques (SS) were used to elucidate these roles. Four small RNA libraries from siliques in an early stage of development were sequenced, and a total of 814 non-redundant miRNA precursors were identified, representing 65 known and 394 novel miRNAs. Expression analyses revealed that 17 miRNAs were differentially expressed in LS and SS lines. Furthermore, through degradome sequencing, we identified 522 cleavage events. Correlation analysis of the differentially expressed miRNAs and their targets suggested that miR159 and miR319 represses cell proliferation and miR160 regulates auxin signal transduction to control silique length. Additionally, the upregulation of miR2111, miR399, miR827, and miR408 reflected restricted silique development due to inorganic phosphate/copper deficiency. More significantly, high expression of miR160 in rapeseed may repress auxin response factors and result in increased silique length, illustrating that silique length might be regulated via an auxin-response pathway.

Genome-wide identification of microRNAs responsive to Ectropis oblique feeding in tea plant (Camellia sinensis L.)

The tea plant (Camellia sinensis L.) is vulnerable to the geometrid Ectropis oblique; although microRNAs (miRNAs) are important for plant growth, development and stress response, the function of miRNAs in the response of C. sinensis to stress from E. oblique is unclear. To identify E. oblique stress-responsive miRNAs and their target genes in tea plant, three small RNA libraries were constructed from leaves subjected to mechanical wounding (MW), geometrid attack (GA) and from healthy control (CK) leaves. Using high-throughput sequencing, 130 known miRNAs and 512 novel miRNAs were identified; of these, differential expression under GA stress was observed for 36 known and 139 novel miRNAs. Furthermore, 169 GA-responsive and 173 MW-responsive miRNAs were detected by miRNA microarray. The expression patterns of six GA-responsive miRNAs were validated by qRT-PCR. Several target genes for these miRNAs encode various transcription factors, including ethylene-responsive transcription factors and squamosa promoter-binding-like proteins, which suggests that these miRNAs may regulate stress-responsive transcriptional processes in tea plant. The present findings provide novel insights into miRNA-mediated regulatory mechanisms underlying the response to GA stress, and also offer valuable information for development of pest resistance using RNA interference-based strategies in tea plants.

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

We identified and cloned the two precursors of miR158 and its target gene in Brassica campestris ssp. chinensis, which both had high relative expression in the inflorescences. Further study revealed that over-expression of miR158 caused reduced pollen varbility, which was caused by the degradation of pollen contents from the binucleate microspore stage. These results first suggest the role of miR158 in pollen development of Brassica campestris ssp. chinensis. MicroRNAs (miRNAs) play crucial roles in many important growth and development processes both in plants and animals by regulating the expression of their target genes via mRNA cleavage or translational repression. In this study, miR158, a Brassicaceae specific miRNA, was functionally characterized with regard to its role in pollen development of non-heading Chinese cabbage (Brassica campestris ssp. chinensis). Two family members of miR158 in B. campestris, namely bra-miR158a1 and bra-miR158a2, and their target gene bra027656, which encodes a pentatricopeptide repeat (PPR) containing protein, were identified. Then, qRT-PCR analysis and GUS-reporter system revealed that both bra-miR158 and its target gene had relatively high expression levels in the inflorescences. Further study revealed that over-expression of miR158 caused reduced pollen varbility and pollen germination ratio, and the degradation of pollen contents from the binucleate microspore stage was also found in those deformed pollen grains, which led to pollen shrinking and collapse in later pollen development stage. These results first shed light on the importance of miR158 in pollen development of Brassica campestris ssp. chinensis.

Identification and differential regulation of microRNAs in response to methyl jasmonate treatment in Lycoris aurea by deep sequencing

Background

Lycoris aurea is a medicine-valuable and ornamental herb widely distributed in China. Former studied have showed that methyl jasmonate (MJ) treatment could increase the content of glanthamine-a worldwide medicine for symptomatic treatment of Alzheimer’s disease in genus Lycoris plants. To explore the possible role of miRNAs in the regulation of jasmonic acid signaling pathway and uncover their potential correlations, we investigated the expression profiles of small RNAs (sRNAs) and their targets in Lycoris aurea, with MJ treatment by using next-generation deep sequencing.

Results

A total of 365 miRNAs were identified, comprising 342 known miRNAs (representing 60 miRNA families) and 23 novel miRNAs. Among them, 143 known and 11 novel miRNAs were expressed differently under MJ treatment. Quantitative real-time PCR of eight selected miRNAs validated the expression pattern of these loci in response to MJ treatment. In addition, degradome sequencing analysis showed that 32 target genes were validated to be targeted by the 49 miRNAs, respectively. Gene function and pathway analyses showed that these targets such as auxin response factors (ARFs), squamosa promoter-binding like (SPL) proteins, basic helix-loop-helix (bHLH) proteins, and ubiquitin-conjugating enzyme E2 are involved in different plant processes, indicating miRNAs mediated regulation might play important roles in L. aurea response to MJ treatment. Furthermore, several L. aurea miRNAs associated with their target genes that might be involved in Amaryllidaceae alkloids biosynthehsis were also analyzed.

Conclusions

A number of miRNAs with diverse expression patterns, and complex relationships between expression of miRNAs and targets were identified. This study represents the first transcriptome-based analysis of miRNAs in Lycoris and will contribute to understanding the potential roles of miRNAs involved in regulation of MJ response.

Electronic supplementary material

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

Genome-Wide Identification and Characterization of MicroRNAs and Target Genes in Lonicera japonica

MiRNAs function in post-transcriptional regulation of gene expression and play very important roles in plant development. Lonicera japonica is one of the important medicinal plants in China. However, few studies on the discovery of conserved and novel miRNAs from L. japonica were reported. In this study, we employed deep sequencing technology to identify miRNAs in leaf and flower tissues of L. japonica. A total of 22.97 million clean reads from flower and leaf tissues were obtained, which generated 146 conserved miRNAs distributed in 20 families and 110 novel miRNAs. Accordingly, 72 differentially expressed miRNAs (P≤0.001) between leaves and flowers and their potential target genes were identified and validated. The qRT-PCR validation showed that majority of the differentially expressed miRNAs showed significant tissue-specific expression in L. japonica. Furthermore, the miRNA-mRNA and mRNA-mRNA regulatory networks were constructed using Cytoscape software. Taken together, this study identified a large number of miRNAs and target genes in L. japonica, which not only provides the first global miRNA expression profiles, but also sheds light on functional genomics research on L. japonica in the future.

Identification of miRNAs and their targets in wild tomato at moderately and acutely elevated temperatures by high-throughput sequencing and degradome analysis

MicroRNAs (miRNAs) are 19-24 nucleotide (nt) noncoding RNAs that play important roles in abiotic stress responses in plants. High temperatures have been the subject of considerable attention due to their negative effects on plant growth and development. Heat-responsive miRNAs have been identified in some plants. However, there have been no reports on the global identification of miRNAs and their targets in tomato at high temperatures, especially at different elevated temperatures. Here, three small-RNA libraries and three degradome libraries were constructed from the leaves of the heat-tolerant tomato at normal, moderately and acutely elevated temperatures (26/18 °C, 33/33 °C and 40/40 °C, respectively). Following high-throughput sequencing, 662 conserved and 97 novel miRNAs were identified in total with 469 conserved and 91 novel miRNAs shared in the three small-RNA libraries. Of these miRNAs, 96 and 150 miRNAs were responsive to the moderately and acutely elevated temperature, respectively. Following degradome sequencing, 349 sequences were identified as targets of 138 conserved miRNAs, and 13 sequences were identified as targets of eight novel miRNAs. The expression levels of seven miRNAs and six target genes obtained by quantitative real-time PCR (qRT-PCR) were largely consistent with the sequencing results. This study enriches the number of heat-responsive miRNAs and lays a foundation for the elucidation of the miRNA-mediated regulatory mechanism in tomatoes at elevated temperatures.

Tight regulation of the interaction between Brassica napus and Sclerotinia sclerotiorum at the microRNA level

MicroRNAs (miRNAs) are multifunctional non-coding short nucleotide molecules. Nevertheless, the role of miRNAs in the interactions between plants and necrotrophic pathogens is largely unknown. Here, we report the identification of the miRNA repertoire of the economically important oil crop oilseed rape (Brassica napus) and those involved in interacting with its most devastating necrotrophic pathogen Sclerotinia sclerotiorum. We identified 280 B. napus miRNA candidates, including 53 novel candidates and 227 canonical members or variants of known miRNA families, by high-throughput deep sequencing of small RNAs from both normal and S. sclerotiorum-inoculated leaves. Target genes of 15 novel candidates and 222 known miRNAs were further identified by sequencing of degradomes from the two types of samples. MiRNA microarray analysis revealed that 68 miRNAs were differentially expressed between S. sclerotiorum-inoculated and uninoculated leaves. A set of these miRNAs target genes involved in plant defense to S. sclerotiorum and/or other pathogens such as nucleotide binding site-leucine-rich repeat (NBS-LRR) R genes and nitric oxygen and reactive oxygen species related genes. Additionally, three miRNAs target AGO1 and AGO2, key components of post-transcriptional gene silencing (PTGS). Expression of several viral PTGS suppressors reduced resistance to S. sclerotiorum. Arabidopsis mutants of AGO1 and AGO2 exhibited reduced resistance while transgenic lines over-expressing AGO1 displayed increased resistance to S. sclerotiorum in an AGO1 expression level-dependent manner. Moreover, transient over-expression of miRNAs targeting AGO1 and AGO2 decreased resistance to S. sclerotiorum in oilseed rape. Our results demonstrate that the interactions between B. napus and S. sclerotiorum are tightly regulated at miRNA level and probably involve PTGS.

miR156a Mimic Represses the Epithelial-Mesenchymal Transition of Human Nasopharyngeal Cancer Cells by Targeting Junctional Adhesion Molecule A

MicroRNAs (miRNAs) have been documented as having an important role in the development of cancer. Broccoli is very popular in large groups of the population and has anticancer properties. Junctional adhesion molecule A (JAMA) is preferentially concentrated at tight junctions and influences cell morphology and migration. Epithelial-mesenchymal transition (EMT) is a developmental program associated with cancer progression and metastasis. In this study we aimed to investigate the role of miRNAs from broccoli in human nasopharyngeal cancer (NPC). We demonstrated that a total of 84 conserved miRNAs and 184 putative novel miRNAs were found in broccoli by sequencing technology. Among these, miR156a was expressed the most. In addition, synthetic miR156a mimic inhibited the EMT of NPC cells in vitro. Furthermore, it was confirmed that JAMA was the target of miR156a mimic as validated by 3' UTR luciferase reporter assays and western blotting. Knockdown of JAMA was consistent with the effects of miR156a mimic on the EMT of NPC, and the up-regulation of JAMA could partially restore EMT repressed by miR156a mimic. In conclusion, these results indicate that the miR156a mimic inhibits the EMT of NPC cells by targeting the 3' UTR of JAMA. These miRNA profiles of broccoli provide a fundamental basis for further research. Moreover, the discovery of miR156a may have clinical implications for the treatment of patients with NPC.

High-throughput sequencing and degradome analysis reveal neutral evolution of Cercis gigantea microRNAs and their targets

High-throughput sequencing and degradome analysis for Cercis gigantea identified 194 known miRNAs and 23 novel miRNAs with 61 targets. The comparison results of highly conserved miRNAs and non-conserved miRNAs implied that C. gigantea miRNAs were subjected to the neutral evolution. MicroRNAs play a key role in post-transcriptionally regulating gene expression during plant growth, development and other various biological processes. Although numerous miRNAs have been identified and documented, to date, there are no reports on Cercis gigantea (C. gigantea) miRNAs. In this study, C. gigantea miRNAs and their target genes were investigated by extracting RNA from young roots, tender stems, young leaves, and flower buds of C. gigantea to establish a small RNA and a degradome library to further sequence. This study identified 194 known miRNAs belonging to 52 miRNA families and 23 novel miRNAs. Among these, 158 miRNAs from 27 miRNA families were highly conserved and existed in a plurality of plants. In addition, 60 different targets for 30 known families and one target for novel miRNA were identified by high-throughput sequencing and degradome analysis in C. gigantea. The comparison results revealed that highly conserved miRNAs have higher expression levels, more family members and more targets than non-conserved miRNAs, indicating that C. gigantea miRNAs were subjected to the neutral evolution. Meanwhile, these conserved miRNAs were also found to be involved in auxin signal transduction, regulation of transcription, and other developmental processes, which will help further understanding regulatory mechanisms of C. gigantea miRNAs.

Identification of Rapeseed MicroRNAs Involved in Early Stage Seed Germination under Salt and Drought Stresses

Drought and salinity are severe and wide-ranging abiotic stresses that substantially affect crop germination, development and productivity, and seed germination is the first critical step in plant growth and development. To comprehensively investigate small-RNA targets and improve our understanding of miRNA-mediated post-transcriptional regulation networks during Brassica napus seed imbibition under drought and salt stresses, we constructed three small-RNA libraries from B. napus variety ZS11 embryos exposed to salt (200 mM NaCl, denoted "S"), drought (200 g L(-1) PEG-6000, denoted "D"), and distilled water (denoted "CK") during imbibition and sequenced them using an Illumina Genome Analyzer. A total of 11,528,557, 12,080,081, and 12,315,608 raw reads were obtained from the CK, D, and S libraries, respectively. Further analysis identified 85 known miRNAs belonging to 31 miRNA families and 882 novel miRNAs among the three libraries. Comparison of the D and CK libraries revealed significant down-regulation of six miRNA families, miR156, miR169, miR860, miR399, miR171, and miR395, whereas only miR172 was significantly up-regulated. In contrast, comparison of the S library with the CK library showed significant down-regulation of only two miRNA families: miRNA393 and miRNA399. Putative targets for 336, 376, and 340 novel miRNAs were successfully predicted in the CK, D, and S libraries, respectively, and 271 miRNA families and 20 target gene families [including disease resistance protein (DIRP), drought-responsive family protein (DRRP), early responsive to dehydration stress protein (ERD), stress-responsive alpha-beta barrel domain protein (SRAP), and salt tolerance homolog2 (STH2)] were confirmed as being core miRNAs and genes involved in the seed imbibition response to salt and drought stresses. The sequencing results were partially validated by quantitative RT-PCR for both conserved and novel miRNAs as well as the predicted target genes. Our data suggest that diverse and complex miRNAs are involved in seed imbibition, indicating that miRNAs are involved in plant hormone regulation, and may play important roles during seed germination under salt- or drought-stress conditions.

Novel and conserved miRNAs in the halophyte Suaeda maritima identified by deep sequencing and computational predictions using the ESTs of two mangrove plants

BACKGROUND

Although miRNAs are reportedly involved in the salt stress tolerance of plants, miRNA profiling in plants has largely remained restricted to glycophytes, including certain crop species that do not exhibit any tolerance to salinity. Hence, this manuscript describes the results from the miRNA profiling of the halophyte Suaeda maritima, which is used worldwide to study salt tolerance in plants.

RESULTS

A total of 134 conserved miRNAs were identified from unique sRNA reads, with 126 identified using miRBase 21.0 and an additional eight identified using the Plant Non-coding RNA Database. The presence of the precursors of seven conserved miRNAs was validated in S. maritima. In addition, 13 novel miRNAs were predicted using the ESTs of two mangrove plants, Rhizophora mangle and Heritiera littoralis, and the precursors of seven miRNAs were found in S. maritima. Most of the miRNAs considered for characterization were responsive to NaCl application, indicating their importance in the regulation of metabolic activities in plants exposed to salinity. An expression study of the novel miRNAs in plants of diverse ecological and taxonomic groups revealed that two of the miRNAs, sma-miR6 and sma-miR7, were also expressed in Oryza sativa, whereas another two, sma-miR2 and sma-miR5, were only expressed in plants growing under the influence of seawater, similar to S. maritima.

CONCLUSION

The distribution of conserved miRNAs among only 25 families indicated the possibility of identifying a greater number of miRNAs with increase in knowledge of the genomes of more halophytes. The expression of two novel miRNAs, sma-miR2 and sma-miR5, only in plants growing under the influence of seawater suggested their metabolic regulatory roles specific to saline environments, and such behavior might be mediated by alterations in the expression of certain genes, modifications of proteins leading to changes in their activity and production of secondary metabolites as revealed by the miRNA target predictions. Moreover, the auxin responsive factor targeted by sma-miR7 could also be involved in salt tolerance because the target is conserved between species. This study also indicated that the transcriptome of one species can be successfully used to computationally predict the miRNAs in other species, especially those that have similar metabolism, even if they are taxonomically separated.

Identification, evolution, and expression partitioning of miRNAs in allopolyploid Brassica napus

The recently published genome of Brassica napus offers for the first time the opportunity to gain insights into the genomic organization and the evolution of miRNAs in oilseed rape. In this study, 12 small RNA libraries from two B. napus cultivars (Tapidor and Ningyou7) and their four double-haploid lines were sequenced, employing the newly sequenced B. napus genome, together with genomes of its progenitors Brassica rapa and Brassica oleracea. A total of 645 miRNAs including 280 conserved and 365 novel miRNAs were identified. Comparative analysis revealed a high level of genomic conservation of MIRNAs (75.9%) between the subgenomes of B. napus and its two progenitors' genomes, and MIRNA lost/gain events (133) occurred in B. napus after its speciation. Furthermore, significant partitioning of miRNA expressions between the two subgenomes in B. napus was detected. The data of degradome sequencing, miRNA-mediated cleavage, and expression analyses support specific interactions between miRNAs and their targets in the modulation of diverse physiological processes in roots and leaves, as well as in biosynthesis of, for example, glucosinolates and lipids in oilseed rape. These data provide a first genome-wide view on the origin, evolution, and genomic organization of B. napus MIRNAs.

Altered Fruit and Seed Development of Transgenic Rapeseed (Brassica napus) Over-Expressing MicroRNA394

Fruit and seed development in plants is a complex biological process mainly involved in input and biosynthesis of many storage compounds such as proteins and oils. Although the basic biochemical pathways for production of the storage metabolites in plants are well characterized, their regulatory mechanisms are not fully understood. In this study, we functionally identified rapeseed (Brassica napus) miR394 with its target gene Brassica napus LEAF CURLING RESPONSIVENESS (BnLCR) to dissect a role of miR394 during the fruit and seed development. Transgenic rapeseed plants over-expressing miR394 under the control of the cauliflower mosaic virus 35S promoter were generated. miR394 over-expression plants exhibited a delayed flowering time and enlarged size of plants, leaf blade, pods and seed body, but developed seeds with higher contents of protein and glucosinolates (GLS) and lower levels of oil accumulation as compared to wild-type. Over-expression of miR394 altered the fatty acid (FA) composition by increasing several FA species such as C16:0 and C18:0 and unsaturated species of C20:1 and C22:1 but lowering C18:3. This change was accompanied by induction of genes coding for transcription factors of FA synthesis including LEAFY COTYLEDON1 (BnLEC1), BnLEC2, and FUSCA3 (FUS3). Because the phytohormone auxin plays a crucial role in fruit development and seed patterning, the DR5-GUS reporter was used for monitoring the auxin response in Arabidopsis siliques and demonstrated that the DR5 gene was strongly expressed. These results suggest that BnmiR394 is involved in rapeseed fruit and seed development.

High-Throughput Sequencing Reveals Diverse Sets of Conserved, Nonconserved, and Species-Specific miRNAs in Jute

MicroRNAs play a pivotal role in regulating a broad range of biological processes, acting by cleaving mRNAs or by translational repression. A group of plant microRNAs are evolutionarily conserved; however, others are expressed in a species-specific manner. Jute is an agroeconomically important fibre crop; nonetheless, no practical information is available for microRNAs in jute to date. In this study, Illumina sequencing revealed a total of 227 known microRNAs and 17 potential novel microRNA candidates in jute, of which 164 belong to 23 conserved families and the remaining 63 belong to 58 nonconserved families. Among a total of 81 identified microRNA families, 116 potential target genes were predicted for 39 families and 11 targets were predicted for 4 among the 17 identified novel microRNAs. For understanding better the functions of microRNAs, target genes were analyzed by Gene Ontology and their pathways illustrated by KEGG pathway analyses. The presence of microRNAs identified in jute was validated by stem-loop RT-PCR followed by end point PCR and qPCR for randomly selected 20 known and novel microRNAs. This study exhaustively identifies microRNAs and their target genes in jute which will ultimately pave the way for understanding their role in this crop and other crops.

Identification of miRNAs and their targets in transgenic Brassica napus and its acceptor (Westar) by high-throughput sequencing and degradome analysis

MicroRNAs (miRNAs) are a class of noncoding small RNAs (sRNAs) that play many roles in plant growth, development, and the stress response.

Identification of chilling stress-responsive tomato microRNAs and their target genes by high-throughput sequencing and degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are a class of noncoding small RNAs (sRNAs) that are 20-24 nucleotides (nt) in length. Extensive studies have indicated that miRNAs play versatile roles in plants, functioning in processes such as growth, development and stress responses. Chilling is a common abiotic stress that seriously affects plants growth and development. Recently, chilling-responsive miRNAs have been detected in several plant species. However, little is known about the miRNAs in the model plant tomato. 'LA1777' (Solanum habrochaites) has been shown to survive chilling stress due to its various characteristics.

RESULTS

Here, two small RNA libraries and two degradome libraries were produced from chilling-treated (CT) and non-chilling-treated (NT) leaves of S. habrochaites seedlings. Following high-throughput sequencing and filtering, 161 conserved and 236 novel miRNAs were identified in the two libraries. Of these miRNAs, 192 increased in the response to chilling stress while 205 decreased. Furthermore, the target genes of the miRNAs were predicted using a degradome sequencing approach. It was found that 62 target genes were cleaved by 42 conserved miRNAs, while nine target genes were cleaved by nine novel miRNAs. Additionally, nine miRNAs and six target genes were validated by quantitative real-time PCR (qRT-PCR). Target gene functional analysis showed that most target genes played positive roles in the chilling response, primarily by regulating the expression of anti-stress proteins, antioxidant enzyme and genes involved in cell wall formation.

CONCLUSIONS

Tomato is an important model plant for basic biological research. In this study, numerous conserved and novel miRNAs involved in the chilling response were identified using high-throughput sequencing, and the target genes were analyzed by degradome sequencing. The work helps identify chilling-responsive miRNAs in tomato and increases the number of identified miRNAs involved in chilling stress. Furthermore, the work provides a foundation for further study of the regulation of miRNAs in the plant response to chilling stress.

Identification and characterization of microRNAs in oilseed rape (Brassica napus) responsive to infection with the pathogenic fungus Verticillium longisporum using Brassica AA (Brassica rapa) and CC (Brassica oleracea) as reference genomes

Verticillium longisporum, a soil-borne pathogenic fungus, causes vascular disease in oilseed rape (Brassica napus). We proposed that plant microRNAs (miRNAs) are involved in the plant-V. longisporum interaction. To identify oilseed rape miRNAs, we deep-sequenced two small RNA libraries made from V. longisporum infected/noninfected roots and employed Brassica rapa and Brassica oleracea genomes as references for miRNA prediction and characterization. We identified 893 B. napus miRNAs representing 360 conserved and 533 novel miRNAs, and mapped 429 and 464 miRNAs to the AA and CC genomes, respectively. Microsynteny analysis with the conserved miRNAs and their flanking protein coding sequences revealed 137 AA-CC genome syntenic miRNA pairs and 61 AA and 42 CC genome-unique miRNAs. Sixty-two miRNAs were responsive to the V. longisporum infection. We present data for specific interactions and simultaneously reciprocal changes in the expression levels of the miRNAs and their targets in the infected roots. We demonstrate that miRNAs are involved in the plant-fungus interaction and that miRNA168-Argonaute 1 (AGO1) expression modulation might act as a key regulatory module in a compatible plant-V. longisporum interaction. Our results suggest that V. longisporum may have evolved a virulence mechanism by interference with plant miRNAs to reprogram plant gene expression and achieve infection.

Novel miRNA and phasiRNA biogenesis networks in soybean roots from two sister lines that are resistant and susceptible to SCN race 4

The soybean cyst nematode (SCN), Heterodera glycines, is the most devastating pathogen of soybean worldwide. SiRNAs (small interfere RNAs) have been proven to induce the silencing of cyst nematode genes. However, whether small RNAs from soybean root have evolved a similar mechanism against SCN is unknown. Two genetically related soybean sister lines (ZP03-5373 and ZP03-5413), which are resistant and susceptible, respectively, to SCN race 4 infection were selected for small RNA deep sequencing to identify small RNAs targeted to SCN. We identified 71 less-conserved miRNAs-miRNAs* counterparts belonging to 32 families derived from 91 loci, and 88 novel soybean-specific miRNAs with distinct expression patterns. The identified miRNAs targeted 42 genes representing a wide range of enzymatic and regulatory activities. Roots of soybean conserved one TAS (Trans-acting siRNA) gene family with a similar but unique trans-acting small interfering RNA (tasiRNA) biogenesis profile. In addition, we found that six miRNAs (gma-miR393, 1507, 1510, 1515, 171, 2118) guide targets to produce secondary phasiRNAs (phased, secondary, small interfering RNAs) in soybean root. Multiple targets of these phasiRNAs were predicted and detected. Importantly, we also found that the expression of 34 miRNAs differed significantly between the two lines. Seven ZP03-5373-specific miRNAs were differentially expressed after SCN infection. Forty-four transcripts from SCN were predicted to be potential targets of ZP03-5373-specific differential miRNAs. These findings suggest that miRNAs play an important role in the soybean response to SCN.

Identification and characterization of microRNAs related to salt stress in broccoli, using high-throughput sequencing and bioinformatics analysis

BACKGROUND

MicroRNAs (miRNAs) are a new class of endogenous regulators of a broad range of physiological processes, which act by regulating gene expression post-transcriptionally. The brassica vegetable, broccoli (Brassica oleracea var. italica), is very popular with a wide range of consumers, but environmental stresses such as salinity are a problem worldwide in restricting its growth and yield. Little is known about the role of miRNAs in the response of broccoli to salt stress. In this study, broccoli subjected to salt stress and broccoli grown under control conditions were analyzed by high-throughput sequencing. Differential miRNA expression was confirmed by real-time reverse transcription polymerase chain reaction (RT-PCR). The prediction of miRNA targets was undertaken using the Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology (KO) database and Gene Ontology (GO)-enrichment analyses.

RESULTS

Two libraries of small (or short) RNAs (sRNAs) were constructed and sequenced by high-throughput Solexa sequencing. A total of 24,511,963 and 21,034,728 clean reads, representing 9,861,236 (40.23%) and 8,574,665 (40.76%) unique reads, were obtained for control and salt-stressed broccoli, respectively. Furthermore, 42 putative known and 39 putative candidate miRNAs that were differentially expressed between control and salt-stressed broccoli were revealed by their read counts and confirmed by the use of stem-loop real-time RT-PCR. Amongst these, the putative conserved miRNAs, miR393 and miR855, and two putative candidate miRNAs, miR3 and miR34, were the most strongly down-regulated when broccoli was salt-stressed, whereas the putative conserved miRNA, miR396a, and the putative candidate miRNA, miR37, were the most up-regulated. Finally, analysis of the predicted gene targets of miRNAs using the GO and KO databases indicated that a range of metabolic and other cellular functions known to be associated with salt stress were up-regulated in broccoli treated with salt.

CONCLUSION

A comprehensive study of broccoli miRNA in relation to salt stress has been performed. We report significant data on the miRNA profile of broccoli that will underpin further studies on stress responses in broccoli and related species. The differential regulation of miRNAs between control and salt-stressed broccoli indicates that miRNAs play an integral role in the regulation of responses to salt stress.

A genome-wide perspective of miRNAome in response to high temperature, salinity and drought stresses in Brassica juncea (Czern) L

Micro RNAs (miRNAs) are involved in diverse biological processes including adaptive response towards abiotic stresses. To unravel small RNAs and more specifically miRNAs that can potentially regulate determinants of abiotic stress tolerance, next generation sequencing of B. juncea seedlings subjected to high temperature, high salt and drought conditions was carried out. With the help of UEA sRNA workbench software package, 51 conserved miRNAs belonging to 30 miRNA families were identified. As there was limited genomic information available for B. juncea, we generated and assembled its genome sequence at a low coverage. Using the generated sequence and other publically available Brassica genomic/transcriptomic resources as mapping reference, 126 novel (not reported in any plant species) were discovered for the first time in B. juncea. Further analysis also revealed existence of 32 and 37 star sequences for conserved and novel miRNAs, respectively. The expression of selected conserved and novel miRNAs under conditions of different abiotic stresses was revalidated through universal TaqMan based real time PCR. Putative targets of identified conserved and novel miRNAs were predicted in B. rapa to gain insights into functional roles manifested by B. juncea miRNAs. Furthermore, SPL2-like, ARF17-like and a NAC domain containing protein were experimentally validated as targets of miR156, miR160 and miR164 respectively. Investigation of gene ontologies linked with targets of known and novel miRNAs forecasted their involvement in various biological functions.

A contig-based strategy for the genome-wide discovery of microRNAs without complete genome resources

MicroRNAs (miRNAs) are important regulators of many cellular processes and exist in a wide range of eukaryotes. High-throughput sequencing is a mainstream method of miRNA identification through which it is possible to obtain the complete small RNA profile of an organism. Currently, most approaches to miRNA identification rely on a reference genome for the prediction of hairpin structures. However, many species of economic and phylogenetic importance are non-model organisms without complete genome sequences, and this limits miRNA discovery. Here, to overcome this limitation, we have developed a contig-based miRNA identification strategy. We applied this method to a triploid species of edible banana (GCTCV-119, Musa spp. AAA group) and identified 180 pre-miRNAs and 314 mature miRNAs, which is three times more than those were predicted by the available dataset-based methods (represented by EST+GSS). Based on the recently published miRNA data set of Musa acuminate, the recall rate and precision of our strategy are estimated to be 70.6% and 92.2%, respectively, significantly better than those of EST+GSS-based strategy (10.2% and 50.0%, respectively). Our novel, efficient and cost-effective strategy facilitates the study of the functional and evolutionary role of miRNAs, as well as miRNA-based molecular breeding, in non-model species of economic or evolutionary interest.

Applications and challenges of next-generation sequencing in Brassica species

Next-generation sequencing (NGS) produces numerous (often millions) short DNA sequence reads, typically varying between 25 and 400 bp in length, at a relatively low cost and in a short time. This revolutionary technology is being increasingly applied in whole-genome, transcriptome, epigenome and small RNA sequencing, molecular marker and gene discovery, comparative and evolutionary genomics, and association studies. The Brassica genus comprises some of the most agro-economically important crops, providing abundant vegetables, condiments, fodder, oil and medicinal products. Many Brassica species have undergone the process of polyploidization, which makes their genomes exceptionally complex and can create difficulties in genomics research. NGS injects new vigor into Brassica research, yet also faces specific challenges in the analysis of complex crop genomes and traits. In this article, we review the advantages and limitations of different NGS technologies and their applications and challenges, using Brassica as an advanced model system for agronomically important, polyploid crops. Specifically, we focus on the use of NGS for genome resequencing, transcriptome sequencing, development of single-nucleotide polymorphism markers, and identification of novel microRNAs and their targets. We present trends and advances in NGS technology in relation to Brassica crop improvement, with wide application for sophisticated genomics research into agronomically important polyploid crops.

Genome-wide identification and characterization of cadmium-responsive microRNAs and their target genes in radish (Raphanus sativus L.) roots

MicroRNAs (miRNAs) are endogenous non-coding small RNAs that play vital regulatory roles in plant growth, development, and environmental stress responses. Cadmium (Cd) is a non-essential heavy metal that is highly toxic to living organisms. To date, a number of conserved and non-conserved miRNAs have been identified to be involved in response to Cd stress in some plant species. However, the miRNA-mediated gene regulatory networks responsive to Cd stress in radish (Raphanus sativus L.) remain largely unexplored. To dissect Cd-responsive miRNAs and their targets systematically at the global level, two small RNA libraries were constructed from Cd-treated and Cd-free roots of radish seedlings. Using Solexa sequencing technology, 93 conserved and 16 non-conserved miRNAs (representing 26 miRNA families) and 28 novel miRNAs (representing 22 miRNA families) were identified. In all, 15 known and eight novel miRNA families were significantly differently regulated under Cd stress. The expression patterns of a set of Cd-responsive miRNAs were validated by quantitative real-time PCR. Based on the radish mRNA transcriptome, 18 and 71 targets for novel and known miRNA families, respectively, were identified by the degradome sequencing approach. Furthermore, a few target transcripts including phytochelatin synthase 1 (PCS1), iron transporter protein, and ABC transporter protein were involved in plant response to Cd stress. This study represents the first transcriptome-based analysis of miRNAs and their targets responsive to Cd stress in radish roots. These findings could provide valuable information for functional characterization of miRNAs and their targets in regulatory networks responsive to Cd stress in radish.

Construction of a genomewide RNAi mutant library in rice

Long hairpin RNA (hpRNA) transgenes are a powerful tool for gene function studies in plants, but a genomewide RNAi mutant library using hpRNA transgenes has not been reported for plants. Here, we report the construction of a hpRNA library for the genomewide identification of gene function in rice using an improved rolling circle amplification-mediated hpRNA (RMHR) method. Transformation of rice with the library resulted in thousands of transgenic lines containing hpRNAs targeting genes of various function. The target mRNA was down-regulated in the hpRNA lines, and this was correlated with the accumulation of siRNAs corresponding to the double-stranded arms of the hpRNA. Multiple members of a gene family were simultaneously silenced by hpRNAs derived from a single member, but the degree of such cross-silencing depended on the level of sequence homology between the members as well as the abundance of matching siRNAs. The silencing of key genes tended to cause a severe phenotype, but these transgenic lines usually survived in the field long enough for phenotypic and molecular analyses to be conducted. Deep sequencing analysis of small RNAs showed that the hpRNA-derived siRNAs were characteristic of Argonaute-binding small RNAs. Our results indicate that RNAi mutant library is a high-efficient approach for genomewide gene identification in plants.

Identification and characterization of novel and conserved microRNAs in radish (Raphanus sativus L.) using high-throughput sequencing

MicroRNAs (miRNAs) are endogenous, non-coding, small RNAs that play significant regulatory roles in plant growth, development, and biotic and abiotic stress responses. To date, a great number of conserved and species-specific miRNAs have been identified in many important plant species such as Arabidopsis, rice and poplar. However, little is known about identification of miRNAs and their target genes in radish (Raphanus sativus L.). In the present study, a small RNA library from radish root was constructed and sequenced using the high-throughput Solexa sequencing. Through sequence alignment and secondary structure prediction, a total of 545 conserved miRNA families as well as 15 novel (with their miRNA* strand) and 64 potentially novel miRNAs were identified. Quantitative real-time PCR (qRT-PCR) analysis confirmed that both conserved and novel miRNAs were expressed in radish, and some of them were preferentially expressed in certain tissues. A total of 196 potential target genes were predicted for 42 novel radish miRNAs. Gene ontology (GO) analysis showed that most of the targets were involved in plant growth, development, metabolism and stress responses. This study represents a first large-scale identification and characterization of radish miRNAs and their potential target genes. These results could lead to the further identification of radish miRNAs and enhance our understanding of radish miRNA regulatory mechanisms in diverse biological and metabolic processes.

MicroRNAs and their putative targets in Brassica napus seed maturation

Background

MicroRNAs (miRNAs) are 20–21 nucleotide RNA molecules that suppress the transcription of target genes and may also inhibit translation. Despite the thousands of miRNAs identified and validated in numerous plant species, only small numbers have been identified from the oilseed crop plant Brassica napus (canola) – especially in seeds.

Results

Using next-generation sequencing technologies, we performed a comprehensive analysis of miRNAs during seed maturation at 9 time points from 10 days after flowering (DAF) to 50 DAF using whole seeds and included separate analyses of radicle, hypocotyl, cotyledon, embryo, endosperm and seed coat tissues at 4 selected time points. We identified more than 500 conserved miRNA or variant unique sequences with >300 sequence reads and also found 10 novel miRNAs. Only 27 of the conserved miRNA sequences had been previously identified in B. napus (miRBase Release 18). More than 180 MIRNA loci were identified/annotated using the B. rapa genome as a surrogate for the B.napus A genome. Numerous miRNAs were expressed in a stage- or tissue-specific manner suggesting that they have specific functions related to the fine tuning of transcript abundance during seed development. miRNA targets in B. napus were predicted and their expression patterns profiled using microarray analyses. Global correlation analysis of the expression patterns of miRNAs and their targets revealed complex miRNA-target gene regulatory networks during seed development. The miR156 family was the most abundant and the majority of the family members were primarily expressed in the embryo.

Conclusions

Large numbers of miRNAs with diverse expression patterns, multiple-targeting and co-targeting of many miRNAs, and complex relationships between expression of miRNAs and targets were identified in this study. Several key miRNA-target expression patterns were identified and new roles of miRNAs in regulating seed development are suggested. miR156, miR159, miR172, miR167, miR158 and miR166 are the major contributors to the network controlling seed development and maturation through their pivotal roles in plant development. miR156 may regulate the developmental transition to germination.