- Research article
- Open Access
Identification and differential regulation of microRNAs in response to methyl jasmonate treatment in Lycoris aurea by deep sequencing
BMC Genomics volume 17, Article number: 789 (2016)
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.
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.
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.
Lycoris aurea (L’ Hér) Herb is a popular ornamental species of Amaryllidaceae plant widely distributed in China. It belongs to the genus Lycoris and the bulbs are very durable, tolerating the extremes of drought and waterlogging, as well as poor soil conditions. Lycoris species possess the plentiful flower colours and shapes , and are also important of medical values [2, 3]. A variety of secondary metabolites known as Amaryllidaceae alkaloids are found in Lycoris, exhibiting many important pharmacological properties, such as immunostimulatory, anti-tumoral, anti-viral and anti-malarial activities [4, 5]. For example, galanthamine, as a selective and reversible acetylcholinesterase inhibitor [6, 7], has been used worldwide in medicine for symptomatic treatment of Alzheimer’s disease ; lycorine, a pyrrolophenanthridine alkaloid, has been demonstrated to display very promising anti-tumor activity in animal and human cell lines [9, 10].
Many detailed insights in biosynthetic steps of Amaryllidaceae alkaloids production are revealed by the biochemical approaches labeling intermediates. In general, Amaryllidaceae alkaloids are regarded as derivatives of the common precursor 4′-O-methylnorbelladine via intramolecular oxidative phenol-coupling [11, 12]. There are three different groups of Amaryllidaceae alkaloids that are biosynthesized by three modes of intramolecular oxidative C–C phenol coupling (para-ortho’, para-para’ and ortho-para’). The para-ortho’ oxidative coupling leads to galanthamine, para-para’ coupling gives rise to marithidine or crinine derivatives, and ortho-para’ coupling yields lycorine . Additionally, in some plants of Amaryllidaceae family, the improved production of Amaryllidaceae alkaloids was observed when treated with methyl jasmonate (MJ), or other elicitor like ethylene [13–17].
Jasmonic acid (JA) system is a key component in the complex plant hormone signaling systems. It is biosynthesized from α-linolenic acid by the octadecanoid pathway and oxygenation of α-linolenic acid is the initial step. The oxygen has to be inserted in the C-13 position by a lipoxygenase (LOX) to form (13S)-hydroperoxy-octadecatrienoic acid (13-HPOT) . 13-HPOT is then utilized as substrates in the following reaction initiated by allene oxide synthase (AOS) and further converted to (9S, 13S)-12-oxo-phytodienoic acid [(9S, 13S)-OPDA] by allene oxide cyclase (AOC) . For further conversion, OPDA is translocated from chloroplasts to peroxisomes, where JA biosynthesis continues and a series of enzymes including OPDA reductase 3 (OPR3), acyl-CoA oxidase (ACX), multifunctional protein, and L-3-ketoacyl CoA thiolase are involved in. JA and its derivatives, collectively referred to as jasmonates, are important signal molecules detected in a wide spectrum of plant species and function on a lot of biological processes including growth inhibition, senescence, wound response, plant defense and secondary mechanism [18, 19]. For example, it was first demonstrated that an endogenous rise in JA levels upon elicitation with a rough yeast elicitor was associated with the induction of alkaloid synthesis in plant cell cultures . Later, by inducing taxoid biosynthesis with MJ in Taxus cuspidata cells, two cytochrome P450 cDNA clones that encoded hydroxylases of the taxol biosynthetic pathway were identified . Additionally, exogenous MJ elicits massive accumulation of caffeoylputrescine in tomato and this introduction is tightly controlled by the constitutive activation of the jasmonate signaling pathway . Broadly, three major classes of plant secondary metabolites defined as the terpenoids, alkaloids and phenylpropanoids were induced by jasmonates [18, 23].
In recent years, plant sRNAs are getting more and more attention for their regulatory roles in growth, development as well as defense process [24, 25]. They are classified into two major categories: small interfering RNAs (siRNAs) and microRNAs (miRNAs) . miRNAs negatively modulate the expression of a wide range of genes, at the post-transciptional levels by directing the mRNA cleavage or by repressing translation . miRNAs have been shown to play crucial roles in plant responses to a variety of abiotic and biotic stresses, such as nutritional deficiency , drought , salinity , cold , heat , heavy metal stress , and disease . Additionally, many lines of evidences support that miRNA-mediated gene regulation also plays a significant role in plant secondary metabolites biosynthesis [35, 36]. Previously, by using high-throughput sequencing, the taxoid elicitor MJ-regulated Taxus chinensis and Arabidopsis miRNA expression was observed [37, 38]. However, little is sensed about MJ-mediated mechanism on miRNAs.
For investigating roles of miRNAs in the regulation of biological process in L. aurea, we employed high-throughput sequencing technology to survey small RNA pools in L. aurea experimentally treated with MJ. Our results indicate a complex and diverse small RNA population exists in L. aurea, thereby providing a meaningful database to microRNA study in Lycoris species. Moreover, our observation suggests that there be a correlation between the MJ-elicited miRNAs and L. aurea biological process such as Amaryllidaceae alkaloids biosynthesis via a possible transcriptional or post-transcriptional regulation.
Results and discussion
Transcriptome sequencing and de novo assembly analysis
We generated transcriptome database with the Illumina Hiseq2000 system, after sequencing the two L. aurea cDNA libraries from the MJ-treated sample (MJ 100) and its untreated control (CK). In total, each library comprises about 53.45 million raw sequence reads. After removing poly (A) tails, short and low-quality tags, and adaptor contamination, 26,809,842 (CK) and 25,874,478 (MJ100) clean reads with a total of 2,412,885,780 and 2,328,703,020 nucleotides were obtained respectively, for the two pools. Due to the unavailability of the full genome sequences of L. aurea, the assembly software Trinity  was performed for de novo assembly of all the 52,684,320 clean reads. Combined with our GS FLX titanium platform of 454 pyrosequencing transcriptome data reported previously , the entire integrated transcriptome was subsequently used to analyze L. aurea sRNAs and degradome libraries.
High-throughput sequencing of CK and MJ-treated small RNA libraries
To perform a wide discovery of miRNAs in L. aurea, we sequenced six small RNA libraries constructed from control (CK1, CK2 and CK3) and MJ-treated (MJ1, MJ2 and MJ3) samples using the Illumina sequencing platform. A total of more than 29.8 million sequence reads were generated from all the samples, in the range 3.0–6.9 million for individual sample (Additional file 1: Table S1). After filtering out the adapter sequences as well as sequences with low quality, and further removing poly-A sequences and short RNA reads smaller than 18 nucleotides and larger than 40 nucleotides, the total number of clean sequences were reduced to about 21.0 million (Table 1 and Additional file 1: Table S1). The number of unique sRNAs ranged from 0.4 to 1.6 million for individual sample (Table 1). Figure 1 shows size distribution of the sRNAs from the MJ and CK libraries. The size distribution of the filtered sequence reads indicated the high-quality of the data and the majority of the small RNA reads were distributed between 18 to 30 nt. The largest fraction of total sRNAs was 21 nt long in all the samples analyzed (Fig. 1a). This result was consistent with previous reports for other plant species, such as Taxus mairei , Panax ginseng , Pinus contorta , Vitis vinifera  and Saccharina japonica . Additionally, size distribution analysis of the 1,781,274 unique small RNA sequences of all the samples showed that the 24 nt group was the biggest, which accounted for 27.64 % of total unique sequences (Fig. 1b). It suggests that 24 nt sRNAs are the most diverse in L. aurea, which is similar to the results observed in other plant species [42, 46, 47]. Overall, more than 75 % of the unique sRNAs in length were within the range of 21–24 nt. These observations were consistent with DICER-LIKE (DCL) protein cleavage products and those reported in previous studies [48–50].
As a complete L. aurea genome sequence is not yet available, we mapped the sequenced sRNAs to the L. aurea transcriptome. A total of 8,849,610 reads (with an average of 42.07 % total small RNA) representing 986,818 unique reads from all the samples were mapped to the transcriptome of L. aurea (Table 1 and Additional file 1: Table S1). However, only a small proportion of the unique sRNAs (15.31 %; Table 1) mapped to the transcriptome of L. aurea. Former study has found that a majority of miRNA genes in plants were located in intergenic regions . For example, among 26 miRNA sequences identified in Ectocarpus, 17 miRNA sequences are located in introns, eight located in 335 intergenic regions, and one in a local antisense to a transposable element . So it is possible that a high proportion of miRNA loci in L. aurea is also located in the intergenic region and/or introns rather than in the gene exons in L. aurea.
Identification of conserved and novel miRNAs in L. aurea
According to alignment of the remaining sequences from the six libraries against all conserved/known miRNAs in miRBase (Release 20.0) with no mismatches and without gaps, a total of 821,348 reads from all the samples were mapped on the miRBase plant miRNAs, which resulted in the identification of a total of 402, 408, 322, 318, 389, and 385 unique conserved miRNAs from CK1, CK2, CK3, MJ1, MJ2, and MJ3, respectively (Table 1). Many of these miRNAs were present in one or more samples analyzed (Fig. 2a). A comparative analysis of miRNAs identified from each sample led to the identification of a total of 342 non-redundant distinct conserved miRNAs of size 18–24 nt in L. aurea. Of 342 miRNAs, 153 were present in all the samples analyzed (Fig. 2a). A significant fraction of the miRNAs was also identified only from a specific sample (Fig. 2a). The size distribution analysis showed the predominant (61.40 %) representation of 21-nt-long miRNAs (Fig. 2b). About both 15.50 % of miRNAs were 20-nt- and 22-nt-long, whereas only 7.60 % was of the rest of the miRNAs (Fig. 2b). Further, most miRNAs of different lengths harboured a uridine residue at the 5′-end (Fig. 2b). Higher abundance of 21-nt-long miRNAs with uridine as a 5′ terminal nucleotide in L. aurea correlates with other plant species reported before . The 21-nt-long miRNA with 5′-uridine is a characteristic feature of DCL1 cleavage and AGO1 association, which has been found in most known miRNAs [49, 53, 54]. Further, the presence of multiple numbers of DCL and AGO proteins can produce miRNAs with different lengths, first nucleotide specificity and diverse functionality [53, 55–57].
We also analyzed the nucleotide composition of mature miRNAs in L. aurea. Most (77.78 %) of the miRNAs had a GC content within the range of 41–60 % (Additional file 2: Figure S1). The average GC content of mature miRNAs in L. aurea (51.89 %) was similar to that observed in Populus (50 %), grapevine (50 %) and maize (52 %) . It has been suggested that the GC content of miRNAs can be used as critical parameter for determining their biological roles .
Additionally, a total of 7.0 % reads matched to structural non-coding RNAs (snRNA, tRNA, and rRNA), repeat sequences and other sequences (Table 1). After removal of these reads, the 9,892,526 unannotated reads (about 47.20 % of total sRNAs) representing 5,173,454 unique reads were used for novel miRNA prediction using the miRDeep2 pipeline. The processing involved extraction of precursor sequences from the transcriptome by extending mapped reads in the flanking regions, and its potential to form a stem-loop secondary structure. After processing, a total of only 23 novel miRNA precursor candidates for CK and MJ libraries were obtained (Additional file 3: Table S2). We also found that the average expression level of the novel miRNAs was much lower than that of the conserved miRNAs (196 reads per miRNA versus 2402 reads per miRNA), implying that most of novel miRNAs are young miRNAs with imprecise processing and lack of targets . Our sequence analysis showed that the putative pre-miRNAs of each library greatly varied from 55 to 372 nucleotides in length (Additional file 3: Table S2). Some of the stem-loop secondary structures of predictive pre-miRNAs of L. aurea could be found in Additional file 2: Figure S2. In addition, miRDeep2 implemented other criteria, such as seed conservation and presence of miRNA* evidence, which ensured high-confidence prediction of miRNAs.
Identification of miRNA families and expression analysis of candidate miRNAs in L. aurea
Based on sequence similarity, we clustered all the identified L. aurea miRNAs into families. 342 conserved miRNAs were clustered into 60 known families and the number of miRNAs in each family varied (Additional file 4: Table S3). Among them, a few (32) of miRNAs were represented only by a single member (Additional file 4: Table S3). This result indicates a kind of diversity of miRNAs appeared in L. aurea. Moreover, 28 miRNA families out of 310 evolutionarily conserved miRNAs were represented by more than one member, and most of the families (20) comprised at least five members. Such large gene families have also previously been reported in plants [50, 51, 59, 60]. Heatmap by HemI software showing conserved miRNA and corresponding miRNA family expression patterns were also presented. As shown in Fig. 3a, for miRNA family expression, three samples for CK and three samples for MJ were clustered into one group respectively, while the clustering of each miRNA expression showed a different clustering pattern (CK1, CK2 and MJ2 clustering in one group; CK3, MJ1 and MJ3 clustering in one group) (Fig. 3b). Furthermore, we performed principal components analysis (PCA) to assess the variability of samples (Additional file 2: Figure S3). The results showed a similar tendency appeared in Fig. 3b.
Among the 60 known miRNA families, 20 differentially expressed miRNA families were identified according to the average expression levels between CK and MJ (Additional file 4: Table S3). Of these, ten miRNA families were up-regulated in MJ100 compared with CK, whereas ten miRNA families showed down-regulated patterns (Additional file 4: Table S3). MJ has been proven to play a key role in changing small RNA expression profile in plants [37, 38, 50]. For example, 14 miRNAs from seven different families (miR156, miR168, miR169, miR172, miR396, miR480 and miR1310) were down-regulated whereas three miRNAs from two families (miR164 and miR390) were up regulated in Chinese yew . Interestingly, among above mentioned miRNAs, we observed no miRNA was induced by MJ (Additional file 4: Table S3), which exhibit the same expression pattern under MJ treatment in Arabidopsis . In this study, the average expression of miR398 in all the samples was down-regulated by 1.33-fold in MJ treatment when compared with the CK (Additional file 4: Table S3). It has been shown that miR398 play an important role in the response to the various abiotic stresses. For example, down-regulation of Arabidopsis miR398 in response to oxidative stresses is important for two posttranscriptional CSD mRNA accumulation and oxidative stress tolerance . Heat stress rapidly induces Arabidopsis miR398, and this induction triggers a regulatory loop that is critical for thermotolerance . Under freezing stress, Arabidopsis miR398 was repressed and Chrysanthemum dichrum ICE1 [INDUCE OF C-REPEAT BINDING FACTOR (CBF) EXPRESSION 1] regulated freezing tolerance of Arabidopsis partly through the miR398-CSD pathway . On the other hand, it has been performed that JA acts with salicylic acid to confer basal thermotolerance in Arabidopsis thaliana . Recently, jasmonate functions as a critical upstream signal of the ICE-CBF/DRE BINDING FACTOR1 (DREB1) transcriptional pathway to positively regulate Arabidopsis freezing tolerance has also been revealed . Based on these results, we could speculate that miRNA398 probably regulates a complex process in the response to MJ in L. aurea. Additionally, we also observed that the average expression of miR528 in all the samples was up-regulated by 2.33-fold in MJ treatment when compared with the CK (Additional file 4: Table S3). Evidence has proven that constitutive expression of miR528, a conserved monocot-specific small RNA alters plant development and enhances tolerance to salinity stress and nitrogen starvation . It also involves in the regulation of arsenite tolerance .
On the other hand, a few miRNA families, such as miR166, miR396, and miR159, showed extraordinarily high expression levels in all the libraries, whereas others were present in much lower abundance. For example, miR159 family was the most abundant, with a total of 134,227 (CK) and 127,914 (MJ100) reads accounting respectively (Additional file 4: Table S3). We also noticed that miR156, which has been shown to be involved in the regulation of flowering time and floral development [67, 68] appeared decreased but not significant expression changes in MJ treatment when compared with the CK (Additional file 4: Table S3). miR172, another small RNA regulating flowering time, showed almost no differential expression changes under MJ treatemt in L. aurea (Additional file 4: Table S3).
By calculating log2 RPM value, the expression level of the miRNAs including novel miRNAs (Additional file 5: Table S4) was also identified (Fig. 3b). In all, 143 out of 342 conserved miRNAs and 11 out of 23 novel miRNAs were identified as differently expressed miRNAs based on the average expression levels of CK and MJ samples (Additional file 4: Table S3; Additional file 5: Table S4). Among them, most are down-regulated in response to MJ. After restricting our analysis to differentially expressed miRNAs with the following three criteria: Read counts of each sample >1, FDR <0.01 and |Fold change| >1, 15 differentially expressed miRNAs indicating 7 up-regulated and 8 down-regulated miRNAs were observed after the MJ100 treatment (Fig. 4 and Additional file 6: Table S5). These data suggest that the differential regulation of miRNAs account for functions in response to the MJ treatment in L. aurea.
qRT-PCR validation of expression profiles of selected miRNAs
To validate the miRNA expression analysis results, we performed qRT-PCR of eight randomly selected conserved miRNAs in all the samples analyzed (Fig. 5). The qRT-PCR data showed a high degree of the agreement with the expression profiles obtained by small RNA sequencing between the CK and MJ100 libraries under MJ treatment at 6 h. For example, a comparative analysis showed similar (correlation >0.70) expression patterns of half (4 of 8) of these miRNAs obtained via small RNA-seq and qRT-PCR analyses (Fig. 5a). Additionally, we also observed a very good concordance in the expression patterns of miRNAs obtained by both small RNA-seq and qRT-PCR as indicated by the overall correlation coefficient (0.86) (Fig. 5b). Similar result was also observed in other plant species .
Identification and annotation of targets for L. aurea miRNAs
Because miRNAs function by regulating their target genes , identifying the potential target genes of miRNAs is crucially important for understanding miRNA-mediated regulation mechanism during MJ treatment. High-throughput degradome sequencing has been shown to be a valuable and efficient approach to validate and characterize target genes of miRNAs in a variety of plant species . The target genes of the miRNAs under MJ treatment were further studied by degradome library sequencing. In total, ~18.71 million (CK) and ~15.93 million (MJ100) clean reads were obtained. Theses clean reads represent 4 611 090 unique reads from CK and 6 626 101 from MJ100 degradome library (Additional file 7: Table S6). A total of 2,848,018 (61.76 %) and 4,190,778 (63.25 %) unique reads from CK and MJ100 degradome libraries were mapped to the reference database (Additional file 7: Table S6). After identifying degraded targets for each of the miRNAs, the abundance of each sequence was plotted, and the degradation products were grouped into five categories (0–4) according to their relative abundance (Additional file 8: Table S7). In all, 133 targets for the 108 miRNAs were identified, involving 37 miRNA families and 13 novel miRNAs (Additional file 8: Table S7). Among the miRNA families, 11 targeted a single transcript, and the highest number of targets cleaved by a single miRNA was 11 (miR156). Some mRNA targets (such as CL5262.Contig1l), was possibly targeted by two miRNAs (miR156 and miR157).
A BlastX search of the Nr (non-redundant protein sequences) database showed that the miRNA targets shared homology with other plant proteins (Additional file 8: Table S7). A number of the identified targets for the known L. aurea miRNAs were transcription factors, such as the auxin response factor (ARF) family, ethylene-responsive transcription factor, growth-regulating factor, squamosa promoter-binding like (SPL) protein, and GRAS family transcription factor (Fig. 6 and Additional file 8: Table S7). Previous reports have also revealed transcription factors as the predominant targets of miRNAs [48–50, 71, 72]. For example, miR156 and miR157 targeted SPLs, miR172 targeted RAP2 transcirption factors in L. aurea, suggesting that these three conserved miRNAs might play significant roles in regulating flowering time and floral development . Some targets also appeared to be involved in signal transduction pathways, such as ARFs and ethylene-responsive transcription factor (Fig. 6 and Additional file 8: Table S7). During Arabidopsis adventitious root formation controlled by auxin-regulated JA homeostasis, ARF6 and ARF8, targets of the miR167, are positive regulators, whereas ARF17, a target of miR160, is a negative regulator . Our result also showed that miR167 targeted ARF6 (Additional file 8: Table S7). In addition, JA is also demonstrated to play important roles in plant signal transduction, metabolism, disease resistance, and response to environmental stresses . A few transcripts targeted by conserved miRNAs were involved in plant response to biotic and abiotic stresses, such as those encoding heat shock protein (CL11156.Contig1; miR166 and novmiR19), ubiquitin-conjugating enzyme E2 (CL273.Contig1; miR399) and copper/zinc superoxide dismutase (CL4528.Contig1, CL4528.Contig2 and CL2610.Contig4; miR398, miR528, miR6300, and novmiR9). Ubiquitin–proteasome-mediated proteolysis has also been shown to be involved in jasmonate signaling system . So it could be possible that miRNA-mediated regulation pathway affects jasmonate signaling indirectly.
Previously, the MJ-induced Amarylideceae alkaloids accumulation has been reported [13–15]. According to the transcriptome sequencing, we identified and summarized some genes such as phenylalanine ammonia-lyase, TYDC, O-methyltransferase (OMT), CYPs, and N-methyltransferase (NMT), might be implicated in Amarylideceae alkaloid galanthamine biosynthesis [40, 75]. In the present study, a few of them including cytochrome P450 (CYP) (CL8816.Contig2_All) and tyrosine decarboxylase (TYDC) (CL7028.Contig1) were indentified to be targeted by miRNAs (miR396). TYDC is responsible for the conversion of tyrosine to tyramine . CYPs are able to conduct an intramolecular C–C phenol coupling reaction in plant alkaloid biosynthesis, and at least two CYP-families (CYP80 and CYP719) have acquired this ability [76, 77]. Whether the two CYPs identified here are involved in the C–C phenol coupling reaction will need further investigation. Another central feature in the biosynthesis of Amarylideceae alkaloids is methyl transduction at multiple positions of the norbelladine core structure and norgalanthamine, and the reactions are considered to be mediated by OMT and NMT respectively [12, 76]. In our condition, neither miRNA-targeted OMT nor NMT was found. At the moment, it is difficult to conclude whether miRNAs are also involved in the post-transcriptional regulation of Amarylideceae alkaloids biosynthetic genes, as many L. aurea specific miRNA might not be found in the present analysis pipeline.
Plant miRNAs generally direct endonucleolytic cleavage of mRNAs, consistent with the suggestion that plant miRNAs enable rapid clearance of target mRNAs at specific points during plant development [78, 79]. This hypothesis predicts a negative correlation between the expression of a miRNA and its target mRNAs within a given tissue or organ. By comparing the expression levels of the differentially expressed miRNAs with the expression levels of their known and predicted targets, former study has showed that miRNA expression is generally anticorrelated with that of targeted mRNAs . In our study, to investigate the correlation between miRNAs and its target mRNA, a total of 172 data points (average expression levels of miRNAs and their corresponding targets in the control and MJ-treated samples) were shown in the scatter plot (Fig. 7). We also observed a negative but not obvious correlation (correlation = −0.14) in the expression patterns of miRNAs and their targets (Fig. 7).
A number of miRNAs with diverse expression patterns, and complex relationships between expression of miRNAs and targets were identified. This is the first report on the transcriptome/degradome-wide identification of known and novel miRNAs and their targets in response to MJ in L. aurea. The findings could advance our understanding of the functional characterization of miRNAs and their targets in regulating plant response to MJ treatment. This study also provides the first glimpse of the complex regulation of the Amarylideceae alkaloid biosynthetic pathway in L. aurea.
Plant culture and MJ treatment
Seeds of Lycoris aurea were surface sterilized and germinated on half-strength Murashige and Skoog (MS) medium (pH 5.8) in the dark at room temperature for 10 d. Afterwards, the seedlings were transferred into plastic pots containing a mixture of soil and vermiculite (3:1, v/v) and cultured in a growth chamber under 14 h light (25 °C)/10 h dark (22 °C). After 12 months growth, the seedlings were treated with 100 μmol l−1 methyl jasmonate (MJ100) for 0, 1, 3, 6, 12 and 24 h, respectively. Seedlings grown in MJ-free solution were used as controls (CK). At each indicated time, three biological replicates for the samples of seedlings were harvested and immediately frozen in liquid nitrogen and stored at −80 °C. The sample exposed to MJ at 6 h was further used for MJ-treated transcriptome, small RNA and degradome sequencing.
Transcriptome and small RNA sequencing
Total RNA from three biological replicates for the samples was isolated independently using Trizol reagent (Invitrogen) according to the manufacturer’s protocols. Equal amounts of RNA from CK and MJ100 seedlings (6 h) were used for small RNA sequencing. For transcriptome and degradome library construction, the pooled total RNA from three biological replicates of each treatment (CK and MJ) was used. The transcriptome library was prepared using an Illumina TruSeq RNA Sample PrepKit following the manufacturer’s instructions. After removing low-quality reads and reads containing only 3′-RNA adaptors, mRNA transcriptome de novo assembly was performed using the SOAP2 program .
Small RNA library construction and sequencing were carried out as described previously . Briefly, the small RNA fragments with a length of 18–30 nt were separated and purified by polyacrylamide gel electrophoresis, and ligated first to a 5′ adaptor and then to a 3′ RNA adaptor. The adaptor-ligated sRNAs were subsequently reverse-transcribed to single-stranded cDNA using SuperScript II Reverse Transcriptase (Invitrogen). Finally, deep sequencing was performed on the Illumina Genome Analyzer II (Illumina, San Diego, CA, USA) at Shanghai OE Biotech Company.
Sequencing data processing and analysis
Small RNA reads were obtained from Illumina HiSeq™ analysis. After removing the low quality reads, reads shorter than 18 nt, reads larger than 40 nt, and trimming adaptor sequences, clean sequencing reads from small RNA libraries were summarized for length distribution and common/specific sequences between samples, and mapped to the L. aurea transcriptome sequences. Reads with the same sequence were grouped and termed ‘unique reads’. The clean sequences were first used for mapping to an entire integrated transcriptome of L. aurea. The remaining unique sequences were performed for identifying the known/conserved miRNAs from miRBase 20.0 (http://www.mirbase.org/index.shtml) without mismatches and gaps. After that, the sequences matching non-coding RNAs included rRNAs, tRNAs, small nuclear RNAs (snRNAs), and small nucleolar RNAs (snoRNAs) deposited in the Rfam (http://rfam.xfam.org/) and NCBI GenBank (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/genbank/) databases were removed. miRDeep2 software  was then used to predict novel miRNAs from the left unannotated sRNAs and the stem–loop structures of pre-miRNAs were also constructed.
Comparison of miRNA expression profiles between CK and MJ-treated L. aurea seedlings
The abundance of all miRNAs was normalized to the transcript expression level per million reads (RPM). The log2-transformed RPM value for miRNA expression was used to generate heat map by Heml . For fold change of MJ to CK, average abundance of three MJ samples and three CK samples was calculated respectively. Then the formula: fold change = log2 (MJ100/CK) were performed. The significance (P-value) of the miRNA expression difference between the CK and MJ100 libraries was calculated with the formula described previously . The false discovery rate (FDR) was also estimated, to determine the threshold of P-value. Fold change (|log2 (MJ100/CK)| >1) and FDR (<0.01) were combined to identify differentially expressed sequences, which were defined as MJ responsive miRNAs. If the normalized expression was zero, it was adjusted to 0.01.
Construction and analysis of degradome libraries
Two L. aurea degradome libraries were constructed following a previously described method . In brief, poly (A+) RNA was isolated from 200 μg of total RNA using the Oligotex mRNA mini kit (Qiagen), and then a 5′ RNA oligonucleotide adaptor containing an MmeI recognition site was added to the 5′-phosphate of the poly (A+) RNA by T4 RNA ligase. After reverse transcription using oligod(T) and PCR enrichment, the PCR products were digested with MmeI and ligated to a 3′ double-stranded DNA adaptor. Finally, the ligation products were amplified with 20 PCR cycles, gel purified, and sequenced on an Illumina Genome Analyzer II.
Raw sequencing reads were pre-processed to remove adaptor sequences and low-quality sequencing reads, and only 20–21 nt sequences with high quality scores were retained for subsequent analysis. The degradome reads were mapped to the L. aurea transcriptome. Perfect matching sequences were used to validate the predicted miRNA targets using CleveLand4 with default parameters . Alignments with no more than five mismatches and no mismatches at the cleavage site (between the 10th and 11th nucleotides) were retained and scored. Finally, miRNA targets with a P-value of ≤ 0.05 were retained.
Quantitative real-time PCR (qRT-PCR) for miRNAs was used to validate the MJ responsive miRNAs. Six conserved and six novel miRNAs were selected and subjected to qPCR analysis. Addition of poly (A+) tails to sRNAs by poly (A+) polymerase and cDNA synthesis were performed using the One Step Primer Script™ miRNA qPCR Starter Kit (TakaRa). PCRs were carried out in a 20 μl reaction mixture consisting of 2 μl of diluted cDNA, 0.2 μM forward and reverse primer, and 10 μl of 2× SYBR Green PCR Master Mix. The reactions were carried out on an iCycler iQ real-time PCR detection system (BIO-RAD) at 95 °C for 30 s, and 40 cycles of 95 °C for 5 s, 56 °C for 15 s, and 72 °C for 20 s. The expression of U6 snRNA was used as an internal control to normalize for variance in the quantity of RNA and input cDNA. The specificity of each PCR reaction was determined by melting curve analysis. At least two independent biological replicates of each sample and three technical replicates of each biological replicate were analyzed by qRT-PCR. The correlation between sequencing and qRT-PCR based expression analysis results was calculated using the SPSS Version 13.0 software (SPSS Institute, Cary, NC, USA). The primers for selected miRNAs were listed in Additional file 9: Table S8.
- (9S, 13S)-OPDA:
(9S, 13S)-12-oxo-phytodienoic acid
allene oxide cyclase
allene oxide synthase
auxin response factor
copper/zinc superoxide dismutase
DRE binding factor1
false discovery rate
induce of C-repeat binding factor (CBF) expression 1
OPDA reductase 3
small nuclear RNAs
small interfering RNAs
small nucleolar RNAs
squamosa promoter-binding like
Hsu PS, Kurita S, Yu ZZ, Lin JZ. Synopsis of the genus Lycoris (Amaryllidaceae). SIDA. 1994;16:301–31.
Howes MJ, Houghton PJ. Plants used in Chinese and Indian traditional medicine for improvement of memory and cognitive function. Pharmacol Biochem Behav. 2003;75:513–27.
Song JH, Zhang L, Song Y. Alkaloids from Lycoris aurea and their cytotoxicities against the head and neck aquamous cell carcinoma. Fitoterapia. 2014;95:121–6.
Jin Z. Amaryllidaceae and Sceletium alkaloids. Nat Prod Rep. 2009;26:363–81.
Son M, Kim A, Lee J, Park CH, Heo JC, Lee HJ, et al. Ethanol extract of Lycoris radiata induces cell death in B16F10 melanoma via p38-mediated AP-1 activation. Oncol Rep. 2010;24:473–8.
Bores GM, Huger FP, Petko W, Mutlib AE, Camacho F, Rush DK, et al. Pharmacological evaluation of novel Alzheimer’s disease therapeutics: acetylcholinesterase inhibitors related to galanthamine. J Pharmacol Exp Ther. 1996;277:728–38.
Lilienfeld S. Galantamine--a novel cholinergic drug with a unique dual mode of action for the treatment of patients with Alzheimer’s disease. CNS Drug Rev. 2002;8:159–76.
Marco L, do Carmo Carreiras M. Galanthamine, a natural product for the treatment of Alzheimer’s disease. Recent Pat CNS Drug Discov. 2006;1:105–11.
Li Y, Liu J, Tang LJ, Shi YW, Ren W, Hu WX. Apoptosis induced by lycorine in KM3 cells is associated with the G0/G1 cell cycle arrest. Oncol Rep. 2007;17:377–84.
Lamoral-Theys D, Decaestecker C, Mathieu V, Dubois J, Kornienko A, Kiss R, et al. Lycorine and its derivatives for anticancer drug design. Mini Rev Med Chem. 2010;10:41–50.
Eichhorn J, Takada T, Kita Y, Zenk MH. Biosynthesis of the Amaryllidaceae alkaloid galanthamine. Phytochemistry. 1998;49:1037–47.
Bastida J, Berkov S, Torras L, Pigni NB, de Andrade JP, Martinez V, et al. Chemical and biological aspects of Amaryllidaceae alkaloids. In: Muñoz-Torrero D, editor. Recent Advances in Pharmaceutical Sciences. Kerala: Transworld Research Network; 2011. p. 65–100.
Colque R, Viladomat F, Bastida J, Codina C. Improved production of galanthamine and related alkaloids by methyl jasmonate in Narcissus confuses shoot-clumps. Planta Med. 2004;70:1180–8.
Mu HM, Wang R, Li XD, Jiang YM, Wang CY, Quan JP, et al. Effect of abiotic and biotic elicitors on growth and alkaloid accumulation of Lycoris chinensis seedlings. Z Naturforsch C. 2009;64:541–50.
Ptak A, El Tahchy A, Wyzgolik G, Henry M, Laurain-Mattar D. Effects of ethylene on somatic embryogenesis and galanthamine content in Leucojum aestivum L. cultures. Plant Cell Tissue Organ Cult. 2010;102:61–7.
Jiang Y, Xia N, Li X, Shen W, Liang L, Wang C, et al. Molecular cloning and characterization of a phenylalanine ammonia-lyase gene (LrPAL) from Lycoris radiata. Mol Biol Rep. 2011;38:1935–40.
Jiang Y, Xia B, Liang L, Li X, Xu S, Peng F, et al. Molecular and analysis of a phenylalanine ammonia-lyase gene (LrPAL2) from Lycoris radiata. Mol Biol Rep. 2013;40:2293–300.
Wasternack C, Hause B. Jasmonates: biosynthesis, perception, signal transduction and action in plant stress response, growth and development. An update to the 2007 review in Annals of Botany. Ann Bot. 2013;111:1021–58.
Wasternack C. Jasmonates: an update on biosynthesis, signal transduction and action in plant stress response, growth and development. Ann Bot. 2007;100:681–97.
Gundlach H, Müller M, Kutchan TM, Zenk MH. Jasmonic acid is a signal transducer in elicitor-induced plant cell cultures. Proc Natl Acad Sci U S A. 1992;89:2389–93.
Jennewein S, Wildung MR, Chau M, Walker K, Croteau R. Random sequencing of an induced Taxus cell cDNA library for identification of clones involved in taxol biosynthesis. Proc Natl Acad Sci U S A. 2004;101:9149–54.
Chen H, Jones AD, Howe GA. Constitutive activation of the jasmonate signaling pathway enhances the production of secondary metabolites in tomato. FEBS Lett. 2006;580:2540–6.
De Geyter N, Gholami A, Goormachtig S, Goossens A. Transcriptional machineries in jasmonate-elicited plant secondary metabolism. Trends Plant Sci. 2012;17:349–59.
Chen X. Small RNAs and their roles in plant development. Annu Rev Cell Dev Biol. 2009;25:21–44.
Weiberg A, Wang M, Bellinger M, Jin H. Small RNAs: a new paradigm in plant-micobe interactions. Annu Rev Phytopathol. 2014;52:495–516.
Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, et al. Criteria for annotation of plant MicroRNAs. Plant Cell. 2008;20:3186–90.
Li J, Reichel M, Li Y, Millar AA. The functional scope of plant microRNA-mediated silencing. Trends Plant Sci. 2014;19:750–6.
Liang G, Yang F, Yu D. MicroRNA395 mediates regulation of sulfate accumulation and allocation in Arabidopsis thaliana. Plant J. 2010;62:1046–57.
Ni Z, Hu Z, Jiang Q, Zhang H. GmNFYA3, a target gene of miR169, is a positive regulator of plant tolerance to drought stress. Plant Mol Biol. 2013;82:113–29.
Zhou M, Li D, Li Z, Hu Q, Yang C, Zhu L, et al. Constitutive expression of a miR319 gene alters plant development and enhances salt and drought tolerance in transgenic creeping bentgrass. Plant Physiol. 2013;161:1375–91.
Thiebaut F, Rojas CA, Almeida KL, Grativol C, Domiciano GC, Lamb CR, et al. Regulation of miR319 during cold stress in sugarcane. Plant Cell Environ. 2012;35:502–12.
Guan Q, Lu X, Zeng H, Zhang Y, Zhu J. Heat stress induction of miR398 triggers a regulatory loop that is critical for thermotolerance in Arabidopsis. Plant J. 2013;74:840–51.
Zhang LW, Song JB, Shu XX, Zhang Y, Yang ZM. miR395 is involved in detoxification of cadmium in Brassica napus. J Hazard Mater. 2013;250–251:204–11.
Li Y, Lu YG, Shi Y, Wu L, Xu YJ, Huang F, et al. Multiple rice microRNAs are involved in immunity against the blast fungus Magnaporthe oryzae. Plant Physiol. 2014;164:1077–92.
Gou JY, Felippes FF, Liu CJ, Weigel D, Wang JW. Negative regulation of anthocyanin biosynthesis in Arabidopsis by a miR156-targeted SPL transcription factor. Plant Cell. 2011;23:1512–22.
Ng DW, Zhang C, Miller M, Palmer G, Whiteley M, Tholl D, et al. cis- and trans-Regulation of miR163 and target genes confers natural variation of secondary metabolites in two Arabidopsis species and their allopolyploids. Plant Cell. 2011;23:1729–40.
Qiu D, Pan X, Wilson IW, Li F, Liu M, Teng W, et al. High throughput sequencing technology reveals that the taxoid elicitor methyl jasmonate regulates microRNA expression in Chinese yew (Taxus chinensis). Gene. 2009;436:37–44.
Zhang B, Xie D, Jin Z. Global analysis of non-coding small RNAs in Arabidopsis in response to jasmonate treatment by deep sequencing technology. J Integr Plant Biol. 2012;54:73–86.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.
Wang R, Xu S, Jiang Y, Jiang J, Li X, Liang L, et al. De novo sequence assembly and characterization of Lycoris aurea transcriptome using GS FLX titanium platform of 454 pyrosequencing. PLoS One. 2013;8:e60449.
Hao DC, Yang L, Xiao PG, Liu M. Identification of Taxus microRNAs and their targets with high-throughput sequencing and degradome analysis. Physiol Plant. 2012;146:388–403.
Wu B, Wang M, Ma Y, Yuan L, Lu S. High-throughput sequencing and characterization of the small RNA transcriptome reveal features of novel and conserved microRNAs in Panax ginseng. PLoS One. 2012;7:e44385.
Morin RD, Aksay G, Dolgosheina E, Ebhardt HA, Magrini V, Mardis ER, et al. Comparative analysis of the small RNA transcriptomes of Pinus contorta and Oryza sativa. Genome Res. 2008;18:571–84.
Pantaleo V, Szittya G, Moxon S, Miozzi L, Moulton V, Dalmay T, et al. Identification of grapevine microRNAs and their targets using high-throughput sequencing and degradome analysis. Plant J. 2010;62:960–76.
Liu F, Wang W, Sun X, Liang Z, Wang F. Conserved and novel heat stress responsive microRNAs were identified by deep sequencing in Saccharina japonica (Laminariales, Phaeophyta). Plant Cell Environ. 2015;38:1357–67.
Carnavale Bottino M, Rosario S, Grativol C, Thiebaut F, Farrinelli L, Hemerly AS, et al. High-throughput sequencing of small RNA transcriptome reveals salt stress regulated microRNAs in sugarcane. PLoS One. 2013;8:e59423.
Zhang R, Marshall D, Bryan GJ, Hornyik C. Identification and characterization of miRNA transcriptome in potato by high-throughput sequencing. PLoS One. 2013;8:e57233.
Jeong DH, Park S, Zhai J, Gurazada SG, De Paoli E, Meyers BC, et al. Massive analysis of rice small RNAs: mechanistic implications of regulated microRNAs and variants for differential target RNA cleavage. Plant Cell. 2011;23:4185–207.
Breakfield NW, Corcoran DL, Petricka JJ, Shen J, Sae-Seaw J, Rubio-Somoza I, et al. High-resolution experimental and computational profiling of tissue-specific known and novel miRNAs in Arabidopsis. Genome Res. 2012;22:163–76.
Jain M, Chevala VVSN, Garg R. Genome-wide discovery and differential regulation of conserved and novel microRNAs in chickpea via deep sequencing. J Exp Bot. 2014;65:5945–58.
Nozawa M, Miura S, Nei M. Origins and evolution of microRNA genes in plant species. Genome Biol Evol. 2012;4:230–9.
Cock JM, Sterck L, Rouzé P, Scornet D, Allen AE, Amoutzias G, et al. The Ectocarpus genome and the independent evolution of multicellularity in brown algae. Nature. 2010;465:617–21.
Rajagopalan R, Vaucheret H, Trejo J, Bartel DP. A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Gene Dev. 2006;20:3407–25.
Mi S, Cai T, Hu Y, Chen Y, Hodges E, Ni F, et al. Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5′ terminal nucleotide. Cell. 2008;133:116–27.
Wu L, Zhou H, Zhang Q, Zhang J, Ni F, Liu C, et al. DNA methylation mediated by a microRNA pathway. Mol Cell. 2010;38:465–75.
Xie Z, Johansen LK, Gustafson AM, Kasschau KD, Lellis AD, Zilberman D, et al. Genetic and functional diversification of small RNA pathways in plants. PLoS Biol. 2004;2:E104.
Cuperus JT, Fahlgren N, Carrington JC. Evolution and functional diversification of MIRNA genes. Plant Cell. 2011;23:431–42.
Mishra AK, Agarwal S, Jain CK, Rani V. High GC content: critical parameter for predicting stress regulated miRNAs in Arabidopsis thaliana. Bioinformation. 2009;4:151–4.
Maher C, Stein L, Ware D. Evolution of Arabidopsis microRNA families through duplication events. Genome Res. 2006;16:510–9.
Li A, Mao L. Evolution of plant microRNA gene families. Cell Res. 2007;17:212–8.
Sunkar R, Kapoor A, Zhu JK. Posttranscriptional induction of two Cu/Zn superoxide dismutase genes in Arabidopsis is mediated by downregulation of miR398 and important for oxidative stress tolerance. Plant Cell. 2006;18:2051–65.
Chen Y, Jiang J, Song A, Chen S, Shan H, Luo H, et al. Ambient temperature enhanced freezing tolerance of Chrysanthemum dichrum CdICE1 Arabidopsis via miR398. BMC Biol. 2013;11:121.
Clarke SM, Cristescu SM, Miersch O, Harren FJM, Wasternack C, Mur LA. Jasmonates act with salicylic acid to confer basal thermotolerance in Arabidopsis thaliana. New Phytol. 2009;182:175–87.
Hu Y, Jiang L, Wang F, Yu D. Jasmonate regulates the inducer of cbf expression-C-repeat binding factor/DRE binding factor1 cascade and freezing tolerance in Arabidopsis. Plant Cell. 2013;25:2907–24.
Yuan S, Li Z, Li D, Yuan N, Hu Q, Luo H. Constitutive expression of rice microRNA528 alters plant development and enhances tolerance to 17 salinity stress and nitrogen starvation in creeping bentgrass. Plant Physiol. 2015;169:576–93.
Liu Q, Hu H, Zhu L, Li R, Feng Y, Zhang L, et al. Involvement of miR528 in the regulation of arsenite tolerance in rice (Oryza sativa L.). J Agric Food Chem. 2015;63:8849–61.
Wang JW, Czech B, Weigel D. miR156-regulated SPL transcription factors define an endogenous flowering pathway in Arabidopsis thaliana. Cell. 2009;138:738–49.
Wu G, Park MY, Conway SR, Wang JW, Weigel D, Poethig RS. The sequential action of miR156 and miR172 regulates developmental timing in Arabidopsis. Cell. 2009;138:750–9.
Rogers K, Chen X. Biogenesis, turnover, and mode of action of plant microRNAs. Plant Cell. 2013;25:2383–99.
German MA, Pillay M, Jeong DH, Hetawal A, Luo S, Janardhanan P, et al. Global identification of microRNA–target RNA pairs by parallel analysis of RNA ends. Nature Biotechnol. 2008;26:941–6.
Xu L, Wang Y, Zhai L, Xu Y, Wang L, Zhu X, et al. Genome-wide identification and characterization of cadmium-responsive microRNAs and their target genes in radish (Raphanus sativus L.) roots. J Exp Bot. 2013;64:4271–87.
Xu MY, Dong Y, Zhang QX, Zhang L, Luo YZ, Sun J, et al. Identification of miRNAs and their targets from Brassica napus by high-throughput sequencing and degradome analysis. BMC Genomics. 2012;13:421.
Gutierrez L, Mongelard G, Floková K, Pacurar DI, Novák O, Staswick P, et al. Auxin controls Arabidopsis adventitious root initiation by regulating jasmonic acid homeostasis. Plant Cell. 2012;24:2515–27.
Xu L, Liu F, Lechner E, Genschik P, Crosby WL, Ma H, et al. The SCFCOI1 ubiquitin-ligase complexes are required for jasmonate response in Arabidopsis. Plant Cell. 2002;14:1919–35.
Takos AM, Rook F. Towards a molecular understanding of the biosynthesis of Amaryllidaceae alkaloids in support of their expanding medical use. Int J Mol Sci. 2013;14:11713–41.
Ikezawa N, Iwasa K, Sato F. Molecular cloning and characterization of CYP80G2, a cytochrome P450 that catalyzes an intramolecular C–C phenol coupling of (S)-reticuline in magnoflorine biosynthesis, from cultured Coptis japonica cells. J Biol Chem. 2008;283:8810–21.
Gesell A, Rolf M, Ziegler J, Diaz Chavez ML, Huang FC, Kutchan TM. CYP719B1 is salutaridine synthase, the C–C phenol-coupling enzyme of morphine biosynthesis in opium poppy. J Biol Chem. 2009;284:24432–42.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Axtell MJ, Bowman JL. Evolution of plant microRNAs and their targets. Trends Plant Sci. 2005;13:343–9.
Axtell MJ, Bartel DP. Antiquity of microRNAs and their targets in land plants. Plant Cell. 2005;17:1658–73.
Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, et al. SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009;25:1966–7.
Hafner M, Landgraf P, Ludwig J, Rice A, Ojo T, Lin C, et al. Identification of microRNAs and other small egulatory RNAs using cDNA library sequencing. Methods. 2008;44:3–12.
Friedländer MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40:37–52.
Deng W, Wang Y, Liu Z, Cheng H, Xue Y. HemI: a toolkit for illustrating heatmaps. PLoS One. 2014;9:e111988.
Eldem V, Çelikkol Akçay U, Ozhuner E, Bakır Y, Uranbey S, Unver T. Genome-wide identification of miRNAs responsive to drought in peach (Prunus persica) by high-throughput deep sequencing. PLoS One. 2012;7:e50298.
Addo-Quaye C, Eshoo TW, Bartel DP, Axtell MJ. Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome. Curr Biol. 2008;18:758–62.
Brousse C, Liu Q, Beauclair L, Deremetz A, Axtell MJ, Bouché N. A non-canonical plant microRNA target site. Nucleic Acids Res. 2014;42:5270–9.
We are grateful to Prof. WB Shen and Dr. WK Chen for their critical reviews during the manuscript preparation. We also thank L. Pan from Shanghai OE Biotech Company for the data analysis. This research was financially supported by the National Natural Science Foundation of China (nos. 31270339, 31301798, 31572151), the Industry-Academy Joint Technological Innovations Foundation of Jiangsu Province (BY2014131), the Jiangsu Provincial Public Institutions Program for Research Conditions and Building Capacity (BM2015019), and Jiangsu Provincial Key Laboratory Foundation for Plant EX SITU Conservation (QD201302).
The authors declare that they have no competing interests.
RW designed the research; SX, YMJ, NW performed most of the experiments and data analysis; SX, RW, BX wrote the draft of the paper. YMJ and YLJ carried out part of material collection, RNA extraction and data analysis. XDL, ZZZ, and YKL participated in the preparation of the manuscript. All authors have read and approved the final manuscript.
Statistics of small RNA sequences from six libraries of L.aurea. (DOCX 31 kb)
Nucleotide composition of miRNAs in L. aurea. (A) GC content distribution in mature miRNAs in L. aurea. (B) Average GC and AT content in mature miRNAs in L. aurea. Figure S2. Secondary structure prediction of novel L. aurea miRNA precursors (A) NovmiR3 and (B) NovmiR5. Figure S3. Principal components analysis (PCA) of miRNA expression of CK and MJ samples in L. aurea. (DOCX 1306 kb)
Novel miRNAs predicted by miRDeep2. (XLSX 18 kb)
Conserved miRNA and miRNA families identified by similarity. (XLSX 105 kb)
Expression levels of novel miRNAs in all samples. (XLSX 25 kb)
Validation and comparative relative expression of differentially expressed miRNAs between the CK and MJ100 libraries in L. aurea. (XLSX 19 kb)
Abundance (percentage) of different RNAs from L.aurea CK and MJ100 degradomes. (DOC 31 kb)
miRNA targets identified from the degradome of CK and MJ100. (XLSX 31 kb)
List of primer sequences used for qRT-PCR experiments. (DOCX 25 kb)
About this article
Cite this article
Xu, S., Jiang, Y., Wang, N. et al. Identification and differential regulation of microRNAs in response to methyl jasmonate treatment in Lycoris aurea by deep sequencing. BMC Genomics 17, 789 (2016) doi:10.1186/s12864-016-2645-y
- Amaryllidaceae alkaloids
- Deep sequencing
- Lycoris aurea
- Methyl jasmonate