Skip to main content

Regulation of terpenoid biosynthesis by miRNA in Persicaria minor induced by Fusarium oxysporum

  • The Correction to this article has been published in BMC Genomics 2019 20:627



Persicaria minor (kesum) is an herbaceous plant with a high level of secondary metabolite compounds, particularly terpenoids. These terpenoid compounds have well-established roles in the pharmaceutical and food industries. Although the terpenoids of P. minor have been studied thoroughly, the involvement of microRNA (miRNA) in terpenoid regulation remains poorly understood and needs to be explored. In this study, P. minor plants were inoculated with the pathogenic fungus Fusarium oxysporum for terpenoid induction.


SPME GC-MS analysis showed the highest terpenoid accumulation on the 6th day post-inoculation (dpi) compared to the other treatment time points (0 dpi, 3 dpi, and 9 dpi). Among the increased terpenoid compounds, α-cedrene, valencene and β-bisabolene were prominent. P. minor inoculated for 6 days was selected for miRNA library construction using next generation sequencing. Differential gene expression analysis showed that 58 miRNAs belonging to 30 families had significantly altered regulation. Among these 58 differentially expressed genes (DEGs), 33 miRNAs were upregulated, whereas 25 miRNAs were downregulated. Two putative novel pre-miRNAs were identified and validated through reverse transcriptase PCR. Prediction of target transcripts potentially involved in the mevalonate pathway (MVA) was carried out by psRobot software, resulting in four miRNAs: pmi-miR530, pmi-miR6173, pmi-miR6300 and a novel miRNA, pmi-Nov_13. In addition, two miRNAs, miR396a and miR398f/g, were predicted to have their target transcripts in the non-mevalonate pathway (MEP). In addition, a novel miRNA, pmi-Nov_12, was identified to have a target gene involved in green leaf volatile (GLV) biosynthesis. RT-qPCR analysis showed that pmi-miR6173, pmi-miR6300 and pmi-nov_13 were downregulated, while miR396a and miR398f/g were upregulated. Pmi-miR530 showed upregulation at 9 dpi, and dynamic expression was observed for pmi-nov_12. Pmi-6300 and pmi-miR396a cleavage sites were detected through degradome sequence analysis. Furthermore, the relationship between miRNA metabolites and mRNA metabolites was validated using correlation analysis.


Our findings suggest that six studied miRNAs post-transcriptionally regulate terpenoid biosynthesis in P. minor. This regulatory behaviour of miRNAs has potential as a genetic tool to regulate terpenoid biosynthesis in P. minor.


Humans exploit the secondary metabolites from medicinal plants in the form of flavouring agents, perfumes, insecticides, dyes and drugs, among other products. Previous studies have shown that more than 100,000 phytochemicals have been isolated from different plant sources [1]. According to the British Nutrition Foundation, secondary metabolites (SMs) are divided into four major classes: terpenoids (volatile compounds, carotenoids and glycosides), phenolic compounds (phenolic acids, tannins and flavonoids), nitrogen-containing compounds (cyanogenic glucosides and alkaloids) and sulfur-containing compounds (thionine, defensin and lectin) [2]. The terpenoid classes consist of volatile and non-volatile compounds. Though these compounds exist in complex structures, they all are made up of isoprene units (C5). Most volatile terpenoids are monoterpenes (C10) and sesquiterpenes (C15), while non-volatile terpenoids include diterpenes (C20), triterpenes (C30) and tetraterpenes (C40) [3]. Terpenoids are widely used in pharmaceuticals, cosmetic fragrances and the food industry [4]. To date, more than 40,000 known terpenoid compounds have been identified, mostly from plants [5].

Persicaria minor (locally known as kesum) is an aromatic plant that is widespread in Southeast Asia, especially in Malaysia, Thailand, Laos, Indonesia and Vietnam. This plant has been used in dishes and medical purposes for many years. The Malaysian government enlisted P. minor in the National Agro Food Policy to ensure its adequate supply and as a platform to strengthen the agricultural-based economy [6]. P. minor has also gained attention due to its abundance of secondary metabolites, especially flavonoids and terpenoids, which are used in the food, fragrance and pharmaceutical industries [7]. P. minor leaves produce the highest content of terpenoids compared to other organs [8]. These compounds play important roles in attracting pollinators, plant defence, and interactions with unfavourable environments. Two different pathways of terpenoid biosynthesis have been identified: the mevalonate (MVA) and non-mevalonate, or methylerythritol 4-phosphate (MEP), pathways [9]. These pathways operate in the cytoplasm and plastid, respectively [9].

Past achievements in unlocking the genes and enzymes involved in each of these pathways has opened up new possibilities for the assessment of terpenoid biosynthesis [10]. At the molecular level, molecules such as DNA, mRNA, transcription factors, and non-coding RNA play roles in regulating the genes involved in these pathways. In addition, these molecules may interact with each other, which can affect the target gene, thus influencing the production of the secondary metabolite [11]. However, the level of secondary metabolites in plants is relatively low because their production depends on plant species, environmental factors, and nutritional sources [12, 13]. In addition, the type of secondary metabolites produced by the plant is stimulus-dependent. For example, in the majority of cases, terpenoid compounds are produced in high abundance under biotic stress compared to abiotic [14]. Hence, the elicitation technique was introduced to enhance secondary metabolite production in plants using a biotic elicitor [13, 15]. Elicitation using fungi, especially F. oxysporum, has been reported to enhance terpenoid contents in plants [16,17,18,19].

To understand the changes at the molecular level, the relationship between miRNA and mRNA needs to be explored, because the interaction between these two RNA species may play a significant role in plant secondary metabolite production. miRNAs are a group of non-coding RNAs with small sizes (~ 20 bp) that act as a gene regulators by negatively regulating mRNAs via cleavage or translational inhibition [11, 20]. The first miRNA, lin-4 miRNA, was discovered in C. elegans through forward genetic analysis. Later, several methods were developed to discover these tiny RNAs, such as cloning and computational methods [20]. However, these methods had their own limitations. Due to the ability of several miRNAs to target single mRNAs, using forward genetic techniques for miRNA discovery is very difficult unless the miRNA of interest has only one target [11, 20]. The cloning method is also limited to highly expressed miRNAs [20]. Computational approaches could be useful for miRNA discovery but have limitations in detecting novel miRNAs in other species, since these approaches are based on homology searches. In addition, computational approaches still require experimental validation [20]. Next generation sequencing (NGS), a more recent approach, has also been introduced for miRNA study, and this method has been found to be more advanced and efficient. This approach can detect low copy number miRNAs at levels of one transcript per million [21, 22]. To date, a total of 38,589 miRNAs from animals, plants, and viruses have been discovered and deposited in public miRNA databases [23]. These miRNAs have the ability to fine-tune many biological and physiological processes such as plant development and stress response by regulating a number of target genes in a time-sensitive manner [11].

Interestingly, miRNAs were reported to be involved in plant secondary metabolite regulation, such as terpenoid, phenolic and nitrogen-containing compound biosynthesis [1]. These interactions could hold potential for future genetic manipulation. In addition, using miRNA as a tool for genetic engineering could be highly beneficial because single miRNAs can target more than single mRNAs and vice versa. Authenticated in silico miRNA discovery and wet lab approaches for miRNA target validation are required prior to their application to plant genetic systems. Thus, in this study, the regulatory roles of miRNAs in the terpenoid biosynthesis pathway were explored in P. minor induced by F. oxysporum. Metabolite profiling was carried out for the F. oxysporum-inoculated plants to determine the changes in the accumulation of terpenoid compounds. In addition, small RNA libraries were constructed, which led to the characterization and analysis of miRNAs involved in the terpenoid biosynthesis pathway.

Results and discussion

Testing of F. oxysporum inoculation and terpenoid profiling in P. minor by SPME-GCMS

To the best of our knowledge, no work on F. oxysporum inoculation to induce terpenoid contents has been reported to date in P. minor. Therefore, we carried out an inoculation test to determine the elicitation ability of this fungus to increase terpenoid content in P. minor. Comparison between mock-inoculated (C) and Fusarium-treated (F) plants showed different leaf morphologies over the time point of the treatment (Fig. 1). From the beginning of the treatment (0 dpi) to 3 dpi, there were no morphological changes observed in either of the C or F samples. At 6 dpi, no change was observed in the C samples, while the F samples showed wilting symptoms beginning at the tip of the leaf. At 9 dpi, the wilting spread to the middle of the leaves in the F samples. At 12 dpi, the wilting spread to the whole leaf, which led to the death of the plant. In comparison, no changes were observed in the C samples until the end point, 12 dpi. A previous report mentioned an accumulation of terpenoids when plants were attacked by pathogens [24]. Based on Fig. 1, the first symptom appeared at 6 dpi. Theoretically, the test suggested the highest terpenoid content at 6 dpi in P. minor inoculated with F. oxysporum. However, the SPME-GCMS result was carried out to reveal the terpenoid content at each time point. The initial SPME-GCMS data are provided in Additional file 1. The heatmap generated (Fig. 2) indicated a significant increase in three volatile compounds: α-cedrene (1H-3a,7-methanoazulene, 2,3,4,7,8,8a-hexahydro-3,6,8,8-tetramethyl-, [3R-(3.α.,3a.beta.,7.beta.,8a.α.)]-), valencene (naphthalene, 1,2,3,5,6,7,8,8a-octahydro-1,8a-dimethyl-7-(1-methylethenyl)-, [1S-(1.α.,7.α.,8a.α.)]-) and β-bisabolene on 6 dpi. Hence, samples at the 6 dpi time point were selected for small RNA library construction.

Fig. 1

Morphological changes in P. minor leaves inoculated with F. oxysporum in a time series (0 dpi, 3 dpi, 6 dpi, 9 dpi and 12 dpi)

Fig. 2

Heat map representing changes in relative metabolite contents of Fusarium-inoculated and control plants detected by SPME GC-MS experiments. Asterisks indicate statistically significant differences (P < 0.05) among all time points by ANOVA and Fisher’s LSD test. The compounds in the blue and red boxes represent the statistically significant GLV and terpenoids, respectively

High-throughput small RNA sequencing analysis

Small RNA sequencing for the C and F libraries generated approximately 12,409,685 and 42,985,084 reads, respectively. The average numbers and proportions of the different categories are presented in Table 1. Trimming of adaptor index sequences was carried out, and low-quality reads were removed to produce clean reads. Reads with lengths between 18 nt and 30 nt were used for annotation, while the rest were discarded. In total, approximately 7,724,932 and 25,726,524 unique reads were produced for the C and F libraries, respectively. The length distribution of the small RNAs in the C and F libraries is shown in Fig. 3. The most abundant small RNAs were 22 nt in length, followed by 20 nt, which belonged to the F library. Previous studies have reported that more than 60% of all plant miRNAs are 21 nt in length [25]. The clean reads from each P. minor library were annotated according to the miRBase version 21 and Rfam databases. From the total clean reads, the C library generated 111,235 (1.44%) miRNA sequences, which represented 2,610 (0.18%) unique reads. On the other hand, the F library generated 711,498 (2.77%) miRNA sequences, representing 6,131 (0.18%) unique reads. Mapping with Rfam resulted in 15.48 and 18.13% sequences in the C and F libraries, respectively, belonging to non-coding RNAs other than miRNA. The rest of the sequences were marked as unannotated sequences in both libraries.

Table 1 Average statistic for deep sequencing results in C and F sample
Fig. 3

Length distributions of small RNAs in the C and F libraries

Differential expression of miRNAs in C and F libraries

The overall expression of miRNAs in both libraries is shown in the volcano plot (Fig. 4). The plot showed significant changes in the regulation of 58 miRNAs. Among these significantly regulated miRNAs, 31 were upregulated and 27 were downregulated. Additionally, a heatmap was generated to compare the expression of 58 significantly regulated miRNAs in the C and F libraries (Fig. 5). Further details about the significantly regulated miRNAs are documented in Table 2.

Fig. 4

Volcano plot showing overall miRNA expression. The plot was constructed based on the log2 fold change on the x-axis and –log 10 P-values on the y-axis. The blue and red dots in the plot represent miRNAs. The blue dots at positive values on the x-axis show miRNAs that were not significantly upregulated, whereas red dots at positive values on the x-axis showed miRNAs that were significantly upregulated. The blue dots at negative values on the x-axis showed miRNAs that were not significantly downregulated, whereas the red dots at negative values on the x-axis showed miRNAs that were significantly downregulated

Fig. 5

Heatmap representing miRNAs significantly altered by Baggerley’s test (P < 0.05) in the C and F libraries Green colour indicates low expression of miRNA, while red colour indicates high expression of miRNA

Table 2 Statistically significant miRNA under inoculation of F. oxysporum (P < 0.05)

Common and specific miRNAs

A Venn diagram was generated to show the common and specific miRNAs in the C and F libraries (Fig. 6). Of the 58 miRNAs responsive to F. oxysporum treatment, 42 miRNAs were common to both libraries. Two miRNAs, pmi-miR168b and pmi-miR2111, were found to be specific to the C library, whereas 14 miRNAs were specific to the F library.

Fig. 6

Common and specific miRNA sequences in the C and F libraries

Discovery of novel miRNAs

To discover novel miRNA sequences in P. minor, all unannotated small RNAs were searched against the P. minor transcriptome data under accession number SRX669305 [26]. After searching for hairpin structures (Fig. 7) and performing MFEI calculations (Table 3), two unique sequences were identified as putative novel miRNAs in P. minor. These putative novel miRNAs were named pmi-nov_miR12 and pmi-nov_miR13. The secondary structures of both miRNAs were validated through reverse transcriptase PCR (Fig. 8). The sizes of the pmi-nov_miR12 precursor (pre-nov_12) and pmi-nov_miR13 precursor (pre-nov_13) are ~ 80 bp and ~ 75 bp, respectively. The sizes of the validated precursors were slightly smaller than the predicted sizes. These smaller sizes might be due to the primers, which were designed to avoid the bulge structure in the hairpin structures. The bulge region in the hairpin structure could have reduced the efficiency of primer binding to the DNA template [27].

Fig. 7

Stem-loop or pre-miRNA structures for putative novel miRNAs

Table 3 List of novel miRNAs discoveries and their information
Fig. 8

Experimental validation of miRNA stem-loop structure through RT-PCR

Prediction and functional annotation of target genes

In plants, miRNAs are involved in many biological processes by negatively regulating their target genes/transcripts [11, 28]. In this study, we searched for putative target genes against the P. minor transcriptomic library software with mismatch score of 4.0 [29]. This prediction resulted in a total of 111 potential target genes in P. minor. Among them, there were 108 target genes for 58 conserved miRNAs and 3 for two novel miRNAs (Table 4). The majority of the miRNAs have more than one target gene. The target genes were classified according to the WEGO database into cellular component, biological process and molecular function categories (Fig. 9).

Table 4 Target prediction for responsive miRNA in P. minor against F. oxysporum
Fig. 9

Functional classification of target transcripts by WEGO software

Analysis of target transcripts involved in terpenoid pathway

The terpenoid pathway consists of two distinct pathways, MVA and MEP. From the target prediction (Table 4), a total of seven miRNA targets were discovered that seemed to be directly involved in these terpenoid biosynthesis pathways. In the MVA pathway, the target involved were diphosphomevalonate decarboxylase (MVD), targeted by pmi-miR530; sesquiterpene synthase and farnesyl diphosphate synthase (FDS), targeted by pmi-miR6173; 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR), targeted by pmi-miR6300; and mevalonate kinase (MVK), targeted by pmi-nov_13. In the MEP or non-mevalonate pathway, the identified targets were 1-deoxy-d-xylulose-5-phosphate synthase (DXS) and 1-deoxy-d-xylulose-5-phosphate reductoisomerase (DXR), which were targeted by pmi-miR396a and pmi-miR398f/g, respectively. In addition, two more targets, peroxidase and alcohol dehydrogenase (ADH), were targeted by pmi-miR396a and pmi-nov_12, respectively. Both of these targets were included in this analysis, because peroxidase is involved in the early signalling of secondary metabolite biosynthesis, and ADH is an enzyme involved in green leaf volatile (GLV) biosynthesis. GLV compounds have been reported to have similar physiological and functional properties to those of terpenoids [30, 31].

Expression profiles of miRNAs and their targets by RT-qPCR

Real-time quantitative polymerase chain reaction (RT-qPCR) was performed to experimentally validate the expression of five conserved miRNAs (pmi-miR396a, pmi-miR398f/g, pmi-miR530, pmi-miR6173, and pmi-miR6300) and two novel miRNAs (pmi-nov_12 and pmi-nov_13) and their target transcripts. The expression profiles are shown in Fig. 10. Pmi-miR396a showed upregulation to 2.5-fold at 9 dpi. Its target, peroxidase57, was upregulated at 3 dpi and then downregulated at 6 dpi and 9 dpi. Peroxidase is an important enzyme that acts as a scavenger for reactive oxygen species (ROS) [32]. ROS are produced by plants as an early response to stress. Simultaneously, plants require a mechanism to prevent ROS from damaging the cell. Hence, high production of peroxidase during early inoculation (3 dpi) may help neutralize the excessive ROS inside the plant cell. In addition, peroxidase may also be involved in terpenoid biosynthesis. In cucumber inoculated with red spider mites (Tetranychus urticae Koch), increased production of (E, E)-α-farnesene resulted [33]. In A. thaliana, miR396 is encoded by two different gene loci, which lead to the biogenesis of two miR396 members, miR396a and miR396b, targeting the GRF transcription factor, which regulates the number of cells in the leaf [34, 35]. Overexpression of miR396a and miR396b resulted in reduced leaf size [34]. In addition, recent studies have revealed that the application of a target mimic approach in miR396a contributes to plant defence against fungal pathogens. The study suggested that the low activity of miR396a induces the plant defence mechanism through the accumulation of hydrogen peroxidase (H2O2) and callus formation [36]. In rice, a similar approach was used, leading to decreased accumulation of miR396a and increased GRF transcription factor activity. That kind of interaction activates auxin biosynthesis and ARF, thus resulting in greater yield [37]. In addition to peroxidase, another target of pmi-miR396a in P. minor was DXS. Negative correlations between these two were observed at 3 dpi. DXS is involved in the MEP pathway, which produces monoterpene and diterpene as the main products [9, 10]. The descending pattern of DXS expression corresponds with the increase in pmi-miR396a, especially at 3 dpi and 9 dpi. For the target transcript peroxidase57, the target was not strongly regulated with high expression of pmi-miR396a at 3 dpi. This kind of regulation was also demonstrated for miR398 in A. thaliana [38]. In A. thaliana, miR398 targeted two types of copper superoxide dismutase (CSD), CSD1 and CSD2. The upregulation of miR398 consistently leads to the downregulation of CSD1, whereas CSD2 did not show any correlation [38].

Fig. 10

Relative expression of miRNAs with respect to their target transcripts

In P. minor, the expression profile of pmi-miR398f/g showed an increasing pattern from 3 dpi to 9 dpi, and the maximum was observed at 9 dpi. The upregulation of pmi-miR398f/g led to the downregulation of the DXR target transcript. Similar to DXS, DXR is involved in the MEP pathway of terpenoid biosynthesis [10]. The miR398 family is thought to be involved in various plant stress responses by participating in the oxidative burst process [39, 40]. In A. thaliana, the expression of miR398 was downregulated under salinity, oxidative stress, Pseudomonas syringae infection and phosphate deficiency. This downregulation leads to the accumulation of a target transcript, CSD, which acts as a ROS scavenger [39]. In P. minor, pmi-miR398f/g was upregulated upon treatment with F. oxysporum. The expression pattern seems to contrast with that of miR398 in A. thaliana. However, in Medicago truncatula, miR398a/b, which targets copper-containing proteins involved in copper homeostasis, was upregulated under drought conditions [41].

Pmi-miR6300 showed a decreasing pattern from 3 dpi until 9 dpi, resulting in the accumulation of the target transcript, HMGR. The highest accumulation of HMGR was observed at 9 dpi, when RT-qPCR showed an upregulation of more than three-fold. HMGR is a vital enzyme involved in sesquiterpene biosynthesis in the MVA pathway [10]. A very limited amount is known about the role of the miR6300 family in plants. In barley, miR6300 was downregulated under drought treatment [42]. In addition, miR6300 has also been discovered in chickpea [43]. However, the target transcripts of miR6300 in both the mentioned studies are still unknown.

Pmi-miR530 targeted MVD, which is involved in the MVA pathway. No significant change was observed in pmi-miR539 expression until 6 dpi. On the other hand, the target transcript MVD showed a gradual increase from 3 dpi until 6 dpi. At 9 dpi, pmi-miR530 was drastically upregulated (exceeding 12-fold), while MVD was repressed. The spike in pmi-miR530 expression was supported by a previous study in chickpea. In chickpea inoculated with F. oxysporum, there was a 17-fold upregulation of miR530 compared to that in control plants. In addition, miR530 in chickpea targeted zinc knuckle protein, which is involved in regulating plant development [44].

Two target transcripts, FDP and sesquiterpene synthase, which are involved in the MVA pathway, were targeted by pmi-miR6173. RT-qPCR analysis showed a decreasing pattern in pmi-miR6173 expression. The expression patterns for both targets were similar. The expression of target transcripts increased at 3 dpi. However, the expression declined at 6 dpi, and no significant change was observed from 6 dpi to 9 dpi. A decreasing expression pattern of miR6173 was also observed in the herbs Sedum alfredii and Medicago sativa. In S. alfredii, miR6173 was downregulated under cadmium treatment. Target prediction showed that miR6173 in S. alfredii targeted a number of targets, such as calcium-binding EF-hand family protein, ATP synthase subunit α, and aspartic protease. However, further work is needed to discover the role of miR6173 in S. alfredii [45]. In M. sativa, miR6173 was among the miRNAs downregulated under drought treatment. miR6173 in M. sativa targeted splicing factor 3A subunit 2, which plays an important role in mRNA splicing [46].

The single gene involved in GLV biosynthesis, ADH, was targeted by the novel miRNA pmi-nov_12. Our results showed that pmi-nov_12 was downregulated from 3 dpi to 6 dpi and then upregulated at 9 dpi, where a 1.2-fold increase was observed. The expression pattern of pmi-nov_12 was downregulated and upregulated at a series of time points, showing a dynamic expression pattern. This type of expression was previously reported in Populus tremula. In P. tremula treated with salinity and abscisic acid, miR398 was downregulated in the first 48 h and then showed high expression at 72 h. In P. minor, although dynamic regulation occurs in pmi-nov_12, a negative correlation can still be observed between pmi-nov_12 and ADH. ADH was reported to be involved in GLV biosynthesis by catalysing the conversion between alcohol and aldehyde via the hydroperoxide lyase pathway [47]. Another novel miRNA in P. minor, pmi-nov_13, was downregulated in response to F. oxysporum treatment, and its target transcript (MVK) was upregulated, especially at 6 dpi. MVK is involved in the MVA pathway of terpenoid biosynthesis. Notably, pmi-nov_12 and pmi-nov_13 are unique miRNAs in P. minor that have never been reported in another plant species before. However, for pmi-miR6173, pmi-nov_12 and pmi-nov_13, the expression of the target genes was comparable to miRNA expression even with decreasing metabolite content at the end of the treatment (9 dpi). Hence, we suggest that these miRNAs may inhibit mRNA translation, which is another mode of miRNA action [48, 49].

miRNAs as post-transcriptional regulator in P. minor

In plants, miRNAs play roles in various biological processes, such as plant development, signal transduction, stress response, and secondary metabolite regulation [50]. In this study, miRNAs were discovered that regulate both the MVA and MEP pathways in terpenoid biosynthesis in P. minor (Fig. 11). Four miRNAs, pmi-miR6300, pmi-nov_12, pmi-miR530, and pmi-miR6173, were found to be involved in the MVA pathway, whereas two miRNAs, pmi-miR396a and pmi-miR398f/g, were found to be involved in the MEP pathway. Most of the miRNAs regulate upstream genes in the terpenoid biosynthesis pathway.

Fig. 11

Involvement of miRNAs in the terpenoid pathway in P. minor EC 1-deoxy-D-xylulose-5-phosphate synthase; EC 1-deoxy-D-xylulose-5-phosphate reductoisomerase; EC hydroxymethylglutaryl-CoA reductase; EC mevalonate kinase; EC diphosphomevalonate decarboxylase; EC farnesyl diphosphate synthase. The terpenoid biosynthesis backbone pathway was constructed using KEGG software. Suppression symbol (continuous line) indicate the miRNAs had displayed negative relationship against their own target, while dashed suppression symbol indicate the hypothetical effect of miRNAs to inhibit the target via translational inhibition

In the MVA pathway, among the four responsive miRNAs, three (pmi-miR6300, pmi-6173 and pmi-nov_13) were downregulated, whereas the fourth miRNA (pmi-miR530) showed a sharp increase at 9 dpi. The Pmi-miR6300-regulated gene encoded HMGR, a rate-limiting enzyme in the MVA pathway and the first rate-limiting enzyme recognized in the terpenoid biosynthesis pathway [51]. In A. thaliana, the HMGR gene has undergone a duplication process that led to HMG1 and HMG2. HMG1 is expressed throughout the plant, whereas the expression of HMG2 was observed in only the meristem and floral parts of A. thaliana plants [51]. In addition, the expression of HMGR is affected by internal factors such as plant development and phytohormone [52]. HMGR is also affected by external stimuli, such as light, wounding, elicitor treatment and pathogen invasion [52, 53]. The mutant hmg1 in A. thaliana exhibited infertility, premature senescence, and dwarfness [54]. In the MVA pathway, HMGR catalyses the rate-limiting steps that directly affect downstream products. In A. thaliana, the hmg1 mutant showed a 65% reduction in triterpene compound accumulation compared to the wild type [55]. Moreover, elicitor treatments such as fungal inoculation and wounding also increased the expression level of HMGR and subsequently led to higher accumulation of sesquiterpene compounds [53, 56]. In P. minor, terpenoids are produced in small quantities under normal conditions. RT-qPCR results showed that pmi-miR6300 was upregulated at 0 dpi, which led to the suppression of its target, HMGR. Indeed, after treatment, pmi-miR6300 was downregulated and led to the high accumulation of HMGR. In X. strumarium, HMGR was targeted by miR1134 and miR5021 [50]. However, the study of both miRNA and HMGR in X. strumarium is still at the prediction level and requires further experimental validation. In the MVA pathway, the HMGR enzyme catalyses the conversion of HMG-CoA to mevalonate, which is later converted into mevalonate-5-phosphate through the enzyme MVK. Until now, no clear interaction between miRNA and MVK has been reported in plants, while in mice, miR122 targets MVK, which is involved in cholesterol biosynthesis [1, 57]. In this study, we investigated the interaction and reported that pmi-nov_13 targets MVK in P. minor. The second miRNA, pmi-miR530, was found to regulate the MVD enzyme, which catalyses mevalonate diphosphate into IPP. Pmi-miR530 was increased drastically at 9 dpi, which led to the downregulation of its target. As for pmi-nov_13 and MVK, no miRNA in plants has previously been reported to target MVD. In mice, miR124 is involved in hypocholesterolaemia by targeting MVD [58]. In plants, information about the roles of MVK and MVD is quite limited compared to yeast [59]. However, the latest study of MVK in Gingko biloba discovered that overexpression of MVK led to the accumulation of the terpene trilactone [60]. In P. minor, the accumulation of MVK at 6 dpi may contribute to sesquiterpene biosynthesis, although the rate-limiting enzyme had its highest expression at 9 dpi.

In the MEP pathway, two miRNAs (pmi-miR396a and miR398f/g) were discovered to be involved by targeting DXS and DXR. DXS is involved in the upstream reaction in the MEP pathway by catalysing the condensation process between pyruvate and D-glyceraldehyde 3-phosphate to produce DXP [61]. The discovery of DXS was first reported in E. coli and then in A. thaliana via mutation of the cla gene [61, 62]. Transgenic A. thaliana overexpressing the DXS gene resulted in an abundance of terpenoids compared to wild plants, and DXS was categorized as a rate-limiting enzyme in the MEP pathway [61]. Further study of DXS in A. thaliana showed that DXS exists as the paralogues DSX1 and DXS2. Of these two paralogues, only DXS1 was found to be involved in terpenoid biosynthesis [63].

In A. thaliana, the majority of the DXS genes were expressed in photosynthetic organs, such as leaves, and floral parts, while homologues to DXS2 in tomato were expressed in roots and trichomes [64]. In the second step of the MEP pathway, the DXR enzyme converts DXP to MEP via rearrangement of the molecular structure followed by a reduction process by NADPH. This step is reversible. Unlike DXS, the role of DXR as a rate-limiting enzyme is still unclear, and its function may depend on plant species, organ and development stage [65]. However, similar to DXS, DXR genes are distributed in different plant parts and stimulated by light response [65]. The mutant dxr exhibited an albino phenotype, defects in gibberellin and abscisic acid biosynthesis and improper formation of trichomes and stomatal closure [66].

Interestingly, in P. minor, the majority of miRNAs involved in the MVA pathway were downregulated except for pmi-miR530, which was upregulated at 9 dpi. The decreasing pattern exhibited by pmi-miR6300, pmi-nov_13 and pmi-miR6173 resulted in the accumulation of their target transcripts, which most likely regulate sesquiterpene production. These findings were supported by the metabolite profile in Fig. 2, which shows prominent production of α-cedrene, valencene, and β-bisabolene after treatment with F. oxysporum. Meanwhile, the increasing pattern shown by pmi-miR396a and pmi-miR398f/g resulted in the suppression of DXS and DXR in the MEP pathway. The differences in miRNA expression involved in regulating MVA and MEP are probably due to the effect of F. oxysporum, which stimulates terpenoid production in the MVA rather than the MEP pathway. Downstream products in the MVA pathway (sesquiterpene), which play crucial roles in plant defence mechanisms, may influence the selection of the terpenoid pathway in the plant system. In addition, the interaction of two regulators, miRNAs and transcription factors, can also affect terpenoid biosynthesis. In A. thaliana and P. cablin, expression of miR156 reduced the production of sesquiterpene by suppression of the target, SPL [67]. In P. minor, pmi-miR156b/c was predicted to target the SPL gene. However, there is no evidence that this interaction would affect sesquiterpene production, because a previous report showed that the interaction between miR156 and SPL is conserved only in regulating floral development [11, 68].

miRNA cleavage site determination by degradome analysis

Degradome sequencing was carried out to determine the miRNA cleavage site. The sequences were deposited under accession number SRX3921398. The statistical data are shown in Table 5. There were 17,532,759 degradome sequencing reads. Filtering and removal of adaptors revealed 15,493,710 clean reads. Mapping to the P. minor transcriptome resulted in 3,079,840 sequence tags. However, to increase the hit search, another P. minor degradome library, with the accession number SRX3921610, was used together with this degradome library. The results showed that out of the 7 studied miRNAs involved in the terpenoid biosynthesis pathway, cleavage sites for pmi-miR396a and pmi-miR6300 and their targets, peroxidase57 and HMGR, were detected between positions 10 and 11 of the miRNA sequences (Fig. 12), and this result was also supported by previous findings [28, 29].

Table 5 Statistic for degradome sequencing
Fig. 12

Detection of miRNA cleavage sites. The number of sequence tags is shown on the right side. The arrows indicate the positions where the miRNAs cleave the targets. a shows the cleavage site for pmi-miR396a, while b shows the cleavage site for pmi-miR6300

Spearman correlation analysis

Overall, correlation analysis revealed that pmi-miR530 and pmi-miR398f/g had positive correlations with metabolite content, whereas the remaining miRNAs (pmi-396a, pmi-miR6300, pmi-miR6173 and pmi-nov_13) showed negative correlations (Fig. 13a). Pmi-nov_12 was excluded from this analysis because it was not involved in the terpenoid biosynthetic pathway. For the target genes, the majority of them displayed positive correlations with metabolite content, except for HMGR (Fig. 13b). The low abundance of α-cedrene, valencene and β-bisabolene at 9 dpi may have contributed to this negative correlation because HMGR expression continued to increase during the fungus treatment. The Spearman coefficients are included in Additional file 3.

Fig. 13

Spearman correlation analysis. Correlation analysis between miRNA-metabolites (a) and mRNA-metabolites (b). In A, miRNAs are represented by squares, whereas metabolites are represented by circles. In B, mRNAs are represented by squares, whereas metabolites are represented by circles. In both figures, yellow lines indicate positive correlations, and grey lines indicate negative correlations


In this study, a total of 58 conserved and two novel miRNAs responsive to F. oxysporum treatment in P. minor were identified. Five of the 58 conserved miRNAs and a novel miRNA were found to be directly involved with the MVA and MEP pathways. However, some of the miRNAs, such as pmi_miR6173, pmi-nov_12 and pmi-nov_13, might need further confirmation at the protein level against their own targets. This work provides the framework for further exploration of miRNA-mediated regulatory mechanisms in terpenoid biosynthesis in P. minor. Moreover, our study also found some miRNAs targeting mRNAs encoding transcription factors, which suggested the role of miRNAs in regulating various biological processes, including plant growth and development. The above seven studied miRNAs could be utilized to regulate secondary metabolite biosynthesis by manipulating the MVA and MEP pathways through the RNAi mechanism.


Plant and fungus cultures

P. minor plants were cultured for approximately 6 weeks on Murashige and Skoog media in a controlled chamber room at Kompleks Rumah Tumbuhan, Universiti Kebangsaan Malaysia (3° 16′ 14.63″ N, 101° 41′ 11.32″ E) [69]. The F. oxysporum fungal culture was obtained from a microbe collection room at the School of Biosciences and Biotechnology, Universiti Kebangsaan Malaysia. Small pieces of the fungal culture were cut and transferred to PDA plates. Two-week-old F. oxysporum culture was used for the inoculation.

F. oxysporum inoculation and terpenoid profiling analysis

Prior to treatment, P. minor plants were transferred from solid to liquid MS media. The plant culture was maintained in a controlled chamber overnight. F. oxysporum spores were collected by rinsing the PDA with sterile distilled water. Approximately 0.05% of Tween 20 was added to the spore solution. The concentration of F. oxysporum was adjusted to 5 × 107/mL before being added to the plant culture. The treated plants were kept in a controlled room with 16 h light and 8 h dark.

Metabolite profiling in P. minor was carried out at a series of time points, beginning with 3 days post-inoculation (dpi), 6 dpi, and 9 dpi, and using SPME-GCMS. In addition, P. minor plants treated with sterile distilled water were used as a control (0 dpi). Every time point had three biological replicates. The SPME method was adapted for the collection of volatiles. P. minor leaves were harvested (approximately 1 g) in liquid nitrogen and ground into small pieces. The ground tissues were then promptly transferred into labelled SPME vials to avoid evaporation of volatile compounds. The vials were closed with screw caps and then heated at 65 °C for 15 min in a water bath. Volatile compounds were collected by introducing an SPME needle via a septum cap. SPME fibre (100 μM polydimethylsiloxane, PDMS), which was located inside the SPME needle, was injected into the vials to absorb the volatiles. Alternatively, polyacrylate fiber also could be used for volatiles detection [70]. The equilibrium time for SPME fibre in the collection vial was set for 30 min at 55 °C–60 °C. GC-MS analysis was performed on an Agilent 7890A gas chromatograph (GC) directly coupled to the mass spectrophotometer (MS) of an Agilent 5975C inert MSD with a triple axis detector. The column used was a non-polar column, HP-5MS (30 m length × 0.25 mm) diameter and film thickness 0.25 μM. Helium was used as the carrier gas, with a flow rate of 1.3 mL/min. A splitless injection was set at 50 °C hold for 3 min, increased to 100 °C at a rate of 20 °C/ minutes, and held at 250 °C for 3 min. The peaks were identified by searching the NIST/EPA/NIH mass spectral library (version 11), and the results were combined in a GC-MS chromatogram.

Statistical analysis of metabolite profiling

The metabolite profiling was analysed using MetaboAnalyst software ( version 3.0 [71]. The data were normalized according to the sum from each sample. One-way ANOVA was carried out followed by Fisher’s LSD test. The P-value was set at P < 0.05.

Small RNA library construction

To generate high-quality small RNA libraries, the crucial step is to obtain high-quality total RNA. Low-quality starting material can reduce the quantity of small RNA [72]. Total RNA was extracted, and the quality (purity and concentration) was measured using a Nanodrop spectrophotometer (ND-1000) and Qubit, respectively. The integrity of the extracted RNA was determined by Bioanalyzer analysis (Agilent 2100) using an RNA 6000 chip. RNA samples with RNA integrity number (RIN) values over 7.0 were selected for small RNA library construction. All quality control steps were enclosed in Additional file 2. Mock-inoculated (C) and Fusarium-treated (F) small RNA libraries were constructed with two biological replicates, C1 and C2, and F1 and F2, respectively. For the F library, the time point at which P. minor emitted volatiles at the maximum level was selected for small RNA library construction. Mock-inoculated P. minor were prepared by adding sterile distilled water to the plant and used as a control. RNA from both the C and F samples was extracted using PureLink® Plant RNA Reagent (Thermo Fisher Scientific, USA) using the recommended protocol. Briefly, approximately 0.1 g of leaf sample from P. minor was ground and transferred to 1.5 mL. A volume of 0.5 mL of PureLink® Plant RNA Reagent was added and vortexed. The samples were incubated at room temperature for 5 min. Next, the samples were centrifuged at 12,000×g for 2 min at 4 °C. Approximately 450 μL of supernatant was transferred into new tubes, and 0.1 mL of 5 M NaCl and 0.3 mL of chloroform were added and mixed to the supernatant. The mixtures were centrifuged at 12,000×g for 2 min at 4 °C, and approximately 350 μL of supernatant was transferred to new tubes, mixed, and incubated for 10 min at room temperature. The tubes were then centrifuged at 12,000×g for 2 min at 4 °C, and a small pellet was observed at the bottom of each tube. The supernatant was discarded, and 75% ethanol was added to wash the pellet. Then, the tubes were centrifuged again at 12,000×g for 1 min at 4 °C. The ethanol was discarded, and the pellet was diluted with 20 μL of nuclease-free water. Then, DNase treatment was carried out for each RNA sample. The concentration of RNA was determined using a Nanodrop spectrophotometer (Thermo Scientific, USA). The integrity of the RNA was assessed by running a 1% agarose gel and Bioanalyzer analysis (Agilent 2100 Bioanalyzer, USA). A small RNA library was constructed using the NEBNext Small RNA Library preparation kit as reported by a previous study [22]. Then, the libraries were sent to Universiti Malaya, Malaysia for sequencing using an Illumina HiSeq 2500™ in Rapid Run mode.

Bioinformatic analysis

Raw data were analysed using CLC Genomics Workbench version 8 [69]. The adaptors were removed, and the low-quality reads were filtered. Only clean reads with lengths between 18 nt and 30 nt were selected for further analysis. The annotation process was carried out by mapping the clean reads to the miRNA and non-coding-RNA databases miRBase version 21 and Rfam version 12 [23, 73]. A maximum of two mismatches were allowed [74]. Then, each library was 11normalized to transcripts per million (TPM); normalized expression = (actual miRNA count/total count of clean reads) × 1,000,000 [75]. Differential gene expression was calculated using Baggerley’s test [76]. Fold change was calculated using |log2 fold change| ≥ 2, and the P-value was set at P < 0.05 [77].

Novel miRNA prediction

Novel miRNA prediction was carried out based on a previous study [78]. For secondary structure validation, forward and reverse primers were designed using the PrimerQuest Tool by Integrated DNA Technologies (, and the parameters were set to PCR 2 Primer (Table 6). The primer was synthesized at First Base (Malaysia). The PCR products were measured and visualized using 4% agarose gel.

Table 6 List of primer for pre-miRNA validation of putative novel miRNA

Prediction of target transcript and gene ontology

PsRobot software was deployed to predict miRNA targets with a score of 4.0 to obtain more targets [29, 79]. The miRNA libraries generated in this study and the transcriptomic library of P. minor ( were used in this target prediction. Then, the target genes were classified through gene ontology analysis using WEGO software ( [80]. The target genes involved in the terpenoid pathway were simply determined by the gene annotations in the transcriptomic library [81]. The genes were mapped with the KEGG database ( to identify the involvement of these genes in the terpenoid pathway [82].

Validation and expression profile using real-time quantitative PCR (RT-qPCR)

To validate the existence and expression of those miRNAs and target genes involved in the terpenoid pathway under F. oxysporum inoculation, RT-qPCR was carried out using Maxima SYBR Green qPCR Master Mix (Thermo Fisher, USA) at a series of time points: 0 dpi, 3 dpi, 6 dpi, and 9 dpi. For miRNA, the miRNA sequence was used as the forward primer, and the universal primer from Qiagen was used as the reverse primer (Table 7). For the target genes, forward and reverse primers were designed using the PrimerQuest Tool by Integrated DNA Technologies (, and the parameters were set to PCR 2 Primer intercalating dye (Table 8). 5.8 s rRNA and tubulin were used as reference genes for the miRNAs and target genes, respectively. Relative gene expression was calculated using the Livak method [83].

Table 7 List of miRNA primers for RT-qPCR
Table 8 List of target genes primer for RT-qPCR

Degradome sequencing

To determine the miRNA cleavage site, RNA samples from C and F were pooled together and sent to BGI (China) for degradome sequencing. After retrieving the raw data, low-quality reads were discarded, and adaptors were removed using Skewer software [84]. Once again, psRobot was deployed to map the degradome sequence against the reference gene (P. minor transcriptome) and to determine the cleavage site [29].

Correlation analysis

To obtain overall visualization of the relationship between miRNA metabolites and mRNA metabolites, non-parametric Spearman’s rank correlation analysis was carried out using miRTarVis software ( with the default settings [85].

Availability of data and materials

The raw data of small RNA libraries were deposited in NCBI under accession number SRX2645686, SRX2645687, SRX2645684, SRX2645685. The transcriptome data can be access in NCBI under accession number SRX669305. The degradome sequencing data also were submitted to NCBI under accession number SRX3921398 and SRX3921610.

Change history

  • 01 August 2019

    Following publication of the original article [1], the authors reported a number of errors, which are listed in this Correction article. The corrections are marked in bold.



Alcohol dehydrogenase


1-deoxy-d-xylulose-5-phosphate reductoisomerase


1-deoxy-d-xylulose-5-phosphate synthase


Farnesyl diphosphate synthase


Green leaf volatile


3-hydroxy-3-methylglutaryl-coenzyme A


Methylerythritol 4-phosphate






Diphosphomevalonate decarboxylase


Mevalonate kinase




  1. 1.

    Gupta OP, Karkute SG, Banerjee S, Meena NL, Dahuja A. Contemporary understanding of miRNA-based regulation of secondary metabolites biosynthesis in plants. Front Plant Sci. 2017;8(374).

  2. 2.

    Mazid M, Khan TA, Mohammad F. Role of secondary metabolites in defense mechanisms of plants. Biol Med. 2011;3(2):232–49.

    CAS  Google Scholar 

  3. 3.

    Dubey VS, Bhalla R, Luthra R. An overview of the non-mevalonate pathway for terpenoid biosynthesis in plants. J Biosci. 2003;28(5):637–46.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  4. 4.

    Vikram P, Chiruvella KK, Ripain IHA, Arifullah M. A recent review on phytochemical constituents and medicinal properties of kesum (Polygonum minus Huds.). Asian Pac J Trop Biomed. 2014;4(6):430–5.

    PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Abbas F, Ke Y, Yu R, Yue Y, Amanullah S, Jahangir MM, Fan Y. Volatile terpenoids: multiple functions, biosynthesis, modulation and manipulation by genetic engineering. Planta. 2017;246(5):803–16.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  6. 6.

    Ahmad R, Baharum SN, Bunawan H, Lee M, Mohd Noor N, Rohani ER, Ilias N, Zin NM. Volatile profiling of aromatic traditional medicinal plant, Polygonum minus in different tissues and its biological activities. Molecules. 2014;19(11):19220–42.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  7. 7.

    Christapher PV, Parasuraman S, Christina JM, Asmawi MZ, Vikneswaran M. Review on Polygonum minus. Huds, a commonly used food additive in Southeast Asia. Pharm Res. 2015;7(1):1–6.

    CAS  Google Scholar 

  8. 8.

    Baharum SN, Bunawan H, Ghani MA, Mustapha WA, Noor NM. Analysis of the chemical composition of the essential oil of Polygonum minus Huds. using two-dimensional gas chromatography-time-of-flight mass spectrometry (GC-TOF MS). Molecules. 2010;15(10):7006–15.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    Nagegowda DA. Plant volatile terpenoid metabolism: biosynthetic genes, transcriptional regulation and subcellular compartmentation. FEBS Lett. 2010;584(14):2965–73.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  10. 10.

    Tholl D. Biosynthesis and biological functions of terpenoids in plants. Adv Biochem Eng Biotechnol. 2015;148:63–106.

    CAS  PubMed  Google Scholar 

  11. 11.

    AFA S, Sajad M, Nazaruddin N, Fauzi IA, AMA M, Zainal Z, Ismail I. MicroRNA and transcription factor: key players in plant regulatory network. Front Plant Sci. 2017;8(565).

  12. 12.

    Ramakrishna A, Ravishankar GA. Influence of abiotic stress signals on secondary metabolites in plants. Plant Signal Behav. 2011;6(11):1720–31.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Ramirez-Estrada K, Vidal-Limon H, Hidalgo D, Moyano E, Golenioswki M, Cusido RM, Palazon J. Elicitation, an effective strategy for the biotechnological production of bioactive high-added value compounds in plant cell factories. Molecules. 2016;21(2):182.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  14. 14.

    Ee S-F, Mohamed-Hussein Z-A, Othman R, Shaharuddin NA, Ismail I, Zainal Z. Functional characterization of Sesquiterpene synthase from Polygonum minus. Sci World J. 2014;2014:840592.

    Article  CAS  Google Scholar 

  15. 15.

    Zhao J, Davis LC, Verpoorte R. Elicitor signal transduction leading to production of plant secondary metabolites. Biotechnol Adv. 2005;23(4):283–333.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  16. 16.

    Gadzovska Simic S, Tusevski O, Maury S, Hano C, Delaunay A, Chabbert B, Lamblin F, Lainé E, Joseph C, Hagège D. Fungal elicitor-mediated enhancement in phenylpropanoid and naphtodianthrone contents of Hypericum perforatum L. cell cultures. Plant Cell Tissue Organ Cult. 2015;122(1):213–26.

    CAS  Article  Google Scholar 

  17. 17.

    Pereira PS, Ticli FK, SdC F, CMdS B, Lourenço MV. Enhanced triterpene production in Tabernaemontana catharinensis cell suspension cultures in response to biotic elicitors. Quím Nova. 2007;30:1849–52.

    CAS  Article  Google Scholar 

  18. 18.

    Silvar C, Merino F, Diaz J. Resistance in pepper plants induced by Fusarium oxysporum f. sp. lycopersici involves different defence-related genes. Plant Biol (Stuttg). 2009;11(1):68–74.

    CAS  Article  Google Scholar 

  19. 19.

    Naik PM, Al–Khayri JM. Abiotic and Biotic Elicitors–Role in Secondary Metabolites Production through In Vitro Culture of Medicinal Plants. In: Shanker AK, Rijeka SC, editors. Abiotic and Biotic Stress in Plants - Recent Advances and Future Perspectives; 2016. p. Ch. 10. InTech.

    Google Scholar 

  20. 20.

    Jones-Rhoades M, Bartel D, Bartel B. MicroRNAS and their regulatory roles in plants. Annu Rev Plant Biol. 2006;57:19–53.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Motameny S, Wolters S, Nurnberg P, Schumacher B. Next generation sequencing of miRNAs - strategies, Resources and Methods. Genes. 2010;1(1):70–84.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  22. 22.


    Article  Google Scholar 

  23. 23.

    Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42:D68–73.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  24. 24.

    Piasecka A, Jedrzejczak-Rey N, Bednarek P. Secondary metabolites in plant innate immunity: conserved function of divergent chemicals. New Phytol. 2015;206(3):948–64.

    PubMed  Article  PubMed Central  Google Scholar 

  25. 25.

    Xuan P, Guo M, Huang Y, Li W, Huang Y. MaturePred: efficient identification of microRNAs within novel plant pre-miRNAs. PLoS One. 2011;6(11):e27422.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Loke K-K, Rahnamaie-Tajadod R, Yeoh C-C, Goh H-H, Mohamed-Hussein Z-A, Mohd Noor N, Zainal Z, Ismail I. RNA-seq analysis for secondary metabolite pathway gene discovery in Polygonum minus. Genom Data. 2016;7:12–3.

    PubMed  Article  PubMed Central  Google Scholar 

  27. 27.

    Stadhouders R, Pas SD, Anber J, Voermans J, Mes THM, Schutten M. The effect of primer-template mismatches on the detection and quantification of nucleic acids using the 5′ nuclease assay. J Mol Diagn. 2010;12(1):109–17.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Bartel D. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  29. 29.

    Wu HJ, Ma YK, Chen T, Wang M, Wang XJ. PsRobot: a web-based plant small RNA meta-analysis toolbox. Nucleic Acids Res. 2012;40(Web Server issue):W22–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    War AR, Paulraj MG, Ahmad T, Buhroo AA, Hussain B, Ignacimuthu S, Sharma HC. Mechanisms of plant defense against insect herbivores. Plant Signal Behav. 2012;7(10):1306–20.

    PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Niinemets Ü, Kännaste A, Copolovici L. Quantitative patterns between plant volatile emissions induced by biotic stresses and the degree of damage. Front Plant Sci. 2013;4:262.

    PubMed  PubMed Central  Article  Google Scholar 

  32. 32.

    Gor M, Ismail I, Mustapha W, Zainal Z, Noor N, Othman R, Hussein Z. Identification of cDNAs for jasmonic acid-responsive genes in Polygonum minus roots by suppression subtractive hybridization. Acta Physiol Plant. 2011;33(2):283–94.

    CAS  Article  Google Scholar 

  33. 33.

    Mercke P, Kappers IF, Verstappen FW, Vorst O, Dicke M, Bouwmeester HJ. Combined transcript and metabolite analysis reveals genes involved in spider mite induced volatile formation in cucumber plants. Plant Physiol. 2004;135(4):2012–24.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Liu D, Song Y, Chen Z, Yu D. Ectopic expression of miR396 suppresses GRF target gene expression and alters leaf growth in Arabidopsis. Physiol Plant. 2009;136(2):223–36.

    CAS  PubMed  Article  Google Scholar 

  35. 35.

    Baucher M, Moussawi J, Vandeputte OM, Monteyne D, Mol A, Perez-Morga D, El Jaziri M. A role for the miR396/GRF network in specification of organ type during flower development, as supported by ectopic expression of Populus trichocarpa miR396c in transgenic tobacco. Plant Biol (Stuttg). 2013;15(5):892–8.

    CAS  Article  Google Scholar 

  36. 36.

    Soto-Suárez M, Baldrich P, Weigel D, Rubio-Somoza I, San Segundo B. The Arabidopsis miR396 mediates pathogen-associated molecular pattern-triggered immune responses against fungal pathogens. Sci Rep. 2017;7:44898.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  37. 37.

    Gao F, Wang K, Liu Y, Chen Y, Chen P, Shi Z, Luo J, Jiang D, Fan F, Zhu Y, et al. Blocking miR396 increases rice yield by shaping inflorescence architecture. Nat Plants. 2015;2:15196.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  38. 38.

    Jagadeeswaran G, Saini A, Sunkar R. Biotic and abiotic stress down-regulate miR398 expression in Arabidopsis. Planta. 2009;229(4):1009–14.

    CAS  PubMed  Article  Google Scholar 

  39. 39.

    Sunkar R, Li Y-F, Jagadeeswaran G. Functions of microRNAs in plant stress responses. Trends Plant Sci. 2012;17:196–203.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    Zhu C, Ding Y, Liu H. MiR398 and plant stress responses. Physiol Plant. 2011;143(1):1–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  41. 41.

    Trindade I, Capitao C, Dalmay T, Fevereiro MP, Santos DM. miR398 and miR408 are up-regulated in response to water deficit in Medicago truncatula. Planta. 2010;231(3):705–16.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  42. 42.

    Hackenberg M, Gustafson P, Langridge P, Shi B-J. Differential expression of microRNAs and other small RNAs in barley between water and drought conditions. Plant Biotechnol J. 2015;13(1):2–13.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  43. 43.

    Srivastava S, Zheng Y, Kudapa H, Jagadeeswaran G, Hivrale V, Varshney RK, Sunkar R. High throughput sequencing of small RNA component of leaves and inflorescence revealed conserved and novel miRNAs as well as phasiRNA loci in chickpea. Plant Sci. 2015;235:46–57.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  44. 44.

    Kohli D, Joshi G, Deokar AA, Bhardwaj AR, Agarwal M, Katiyar-Agarwal S, Srinivasan R, Jain PK. Identification and characterization of wilt and salt stress-responsive MicroRNAs in chickpea through high-throughput sequencing. PLoS One. 2014;9(10):e108851.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  45. 45.

    Han X, Yin H, Song X, Zhang Y, Liu M, Sang J, Jiang J, Li J, Zhuo R. Integration of small RNAs, degradome and transcriptome sequencing in hyperaccumulator Sedum alfredii uncovers a complex regulatory network and provides insights into cadmium phytoremediation. Plant Biotechnol J. 2016;14(6):1470–83.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Li Y, Wan L, Bi S, Wan X, Li Z, Cao J, Tong Z, Xu H, He F, Li X. Identification of drought-responsive MicroRNAs from roots and leaves of alfalfa by high-throughput sequencing. Genes. 2017;8(4):119.

    PubMed Central  Article  CAS  Google Scholar 

  47. 47.

    Ul Hassan MN, Zainal Z, Ismail I. Green leaf volatiles: biosynthesis, biological functions and their applications in biotechnology. Plant Biotechnol J. 2015;13(6):727–39.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  48. 48.

    Zhang B, Pan X, Cobb GP, Anderson TA. Plant microRNA: a small regulatory molecule with big impact. Dev Biol. 2006;289(1):3–16.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  49. 49.

    Aukerman M, Sakai H. Regulation of flowering time and floral organ identity by a MicroRNA and its APETALA2-like target genes. Plant Cell. 2003;15:2730–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  50. 50.

    Fan R, Li Y, Li C, Zhang Y. Differential microRNA analysis of glandular trichomes and young leaves in Xanthium strumarium L. Reveals their putative roles in regulating terpenoid biosynthesis. PLoS One. 2015;10(9):e0139002.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  51. 51.

    Enjuto M, Lumbreras V, Marin C, Boronat A. Expression of the Arabidopsis HMG2 gene, encoding 3-hydroxy-3-methylglutaryl coenzyme a reductase, is restricted to meristematic and floral tissues. Plant Cell. 1995;7(5):517–27.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Tholl D, Lee S. Terpene specialized metabolism in Arabidopsis thaliana. Arabidopsis Book. 2011;9:e0143.

    PubMed  PubMed Central  Article  Google Scholar 

  53. 53.

    Kondo K, Uritani I, Oba K. Induction mechanism of 3-hydroxy-3-methylglutaryl-CoA reductase in potato tuber and sweet potato root tissues. Biosci Biotechnol Biochem. 2003;67(5):1007–17.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  54. 54.

    Suzuki M, Kamide Y, Nagata N, Seki H, Ohyama K, Kato H, Masuda K, Sato S, Kato T, Tabata S, et al. Loss of function of 3-hydroxy-3-methylglutaryl coenzyme a reductase 1 (HMG1) in Arabidopsis leads to dwarfing, early senescence and male sterility, and reduced sterol levels. Plant J. 2004;37(5):750–61.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  55. 55.

    Ohyama K, Suzuki M, Masuda K, Yoshida S, Muranaka T. Chemical phenotypes of the hmg1 and hmg2 mutants of Arabidopsis demonstrate the in-planta role of HMG-CoA reductase in triterpene biosynthesis. Chem Pharm Bull. 2007;55(10):1518–21.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

    Chappell J, VonLanken C, Vögeli U. Elicitor-inducible 3-Hydroxy-3-Methylglutaryl coenzyme a reductase activity is required for Sesquiterpene accumulation in tobacco cell suspension cultures. Plant Physiol. 1991;97(2):693–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  57. 57.

    Krutzfeldt J, Rajewsky N, Braich R, Rajeev KG, Tuschl T, Manoharan M, Stoffel M. Silencing of microRNAs in vivo with ‘antagomirs’. Nature. 2005;438(7068):685–9.

    PubMed  Article  CAS  Google Scholar 

  58. 58.

    Michihara A, Ide N, Mizutani Y, Okamoto M, Uchida M, Matsuoka H, Akasaki K. Involvement of microRNA214 and transcriptional regulation in reductions in mevalonate pyrophosphate decarboxylase mRNA levels in stroke-prone spontaneously hypertensive rat livers. Biosci Biotechnol Biochem. 2015;79(11):1759–70.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  59. 59.

    Singh B, Sharma RA. Plant terpenes: defense responses, phylogenetic analysis, regulation and clinical applications. 3 Biotech. 2015;5(2):129–51.

    PubMed  Article  PubMed Central  Google Scholar 

  60. 60.

    Chen Q, Yan J, Meng X, Xu F, Zhang W, Liao Y, Qu J. Molecular cloning, characterization, and functional analysis of acetyl-CoA C-acetyltransferase and mevalonate kinase genes involved in terpene trilactone biosynthesis from Ginkgo biloba. Molecules. 2017;22(1):74.

    PubMed Central  Article  CAS  Google Scholar 

  61. 61.

    Estevez JM, Cantero A, Reindl A, Reichler S, Leon P. 1-deoxy-D-xylulose-5-phosphate synthase, a limiting enzyme for plastidic isoprenoid biosynthesis in plants. J Biol Chem. 2001;276(25):22901–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  62. 62.

    Mandel MA, Feldmann KA, Herrera-Estrella L, Rocha-Sosa M, Leon P. CLA1, a novel gene required for chloroplast development, is highly conserved in evolution. Plant J. 1996;9(5):649–58.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  63. 63.

    Phillips MA, Leon P, Boronat A, Rodriguez-Concepcion M. The plastidial MEP pathway: unified nomenclature and resources. Trends Plant Sci. 2008;13(12):619–23.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  64. 64.

    Paetzold H, Garms S, Bartram S, Wieczorek J, Uros-Gracia EM, Rodriguez-Concepcion M, Boland W, Strack D, Hause B, Walter MH. The isogene 1-deoxy-D-xylulose 5-phosphate synthase 2 controls isoprenoid profiles, precursor pathway allocation, and density of tomato trichomes. Mol Plant. 2010;3(5):904–16.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  65. 65.

    Carretero-Paulet L, Galstyan A, Roig-Villanova I, Martinez-Garcia JF, Bilbao-Castro JR, Robertson DL. Genome-wide classification and evolutionary analysis of the bHLH family of transcription factors in arabidopsis, poplar, rice, moss, and algae. Plant Physiol. 2010;153:1398–412.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  66. 66.

    Xing S, Miao J, Li S, Qin G, Tang S, Li H, Gu H, Qu LJ. Disruption of the 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) gene results in albino, dwarf and defects in trichome initiation and stomata closure in Arabidopsis. Cell Res. 2010;20(6):688–700.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  67. 67.

    Yu ZX, Wang LJ, Zhao B, Shan CM, Zhang YH, Chen DF, Chen XY. Progressive regulation of Sesquiterpene biosynthesis in Arabidopsis and patchouli (Pogostemon cablin) by the miR156-targeted SPL transcription factors. Mol Plant. 2014;8:98–110.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  68. 68.

    Hong Y, Jackson S. Floral induction and flower formation—the role and potential applications of miRNAs. Plant Biotechnol J. 2015;13(3):282–92.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  69. 69.

    Samad AFA, Nazaruddin N, Sajad M, Jani J, Murad AMA, Zainal Z, Ismail I. Small RNA sequencing for secondary metabolite analysis in Persicaria minor. Genom Data. 2017;13:3–4.

    PubMed  PubMed Central  Article  Google Scholar 

  70. 70.

    Câmara JS, Arminda Alves M, Marques JC. Development of headspace solid-phase microextraction-gas chromatography–mass spectrometry methodology for analysis of terpenoids in Madeira wines. Analytica Chimica Acta 2006;555(2):191–200.

    Article  CAS  Google Scholar 

  71. 71.

    Xia J, Wishart DS. Using MetaboAnalyst 3.0 for comprehensive metabolomics data analysis. Curr Protoc Bioinformatics. 2016;55:14.10.11–91.

    Article  Google Scholar 

  72. 72.

    Accerbi M, Schmidt SA, De Paoli E, Park S, Jeong DH, Green PJ. Methods for isolation of total RNA to recover miRNAs and other small RNAs from diverse species. Methods Mol Biol. 2010;592:31–50.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  73. 73.

    Nawrocki EP, Burge SW, Bateman A, Daub J, Eberhardt RY, Eddy SR, Floden EW, Gardner PP, Jones TA, Tate J, et al. Rfam 12.0: updates to the RNA families database. Nucleic Acids Res. 2015;43(D1):D130–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  74. 74.

    Yang L, Jue D, Li W, Zhang R, Chen M, Yang Q. Identification of MiRNA from eggplant (Solanum melongena L.) by small RNA deep sequencing and their response to Verticillium dahliae infection. PLoS One. 2013;8(8):e72840.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  75. 75.

    Wang Y, Li L, Tang S, Liu J, Zhang H, Zhi H, Jia G, Diao X. Combined small RNA and degradome sequencing to identify miRNAs and their targets in response to drought in foxtail millet. BMC Genet. 2016;17:57.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  76. 76.

    Baggerly KA, Deng L, Morris JS, Aldaz CM. Differential expression in SAGE: accounting for normal between-library variation. Bioinformatics. 2003;19(12):1477–83.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  77. 77.

    Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7:986–95.

    CAS  PubMed  Article  Google Scholar 

  78. 78.

    AFA S, Nazaruddin N, AMA M, Jani J, Zainal Z, Ismail I. Deep sequencing and in silico analysis of small RNA library reveals novel miRNA from leaf Persicaria minor transcriptome. 3 Biotech. 2018;8(3):136.

    Google Scholar 

  79. 79.

    Yuan J, Zhang Y, Dong J, Sun Y, Lim BL, Liu D, Lu ZJ. Systematic characterization of novel lncRNAs responding to phosphate starvation in Arabidopsis thaliana. BMC Genomics. 2016;17(1):655.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  80. 80.

    Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L, et al. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34(suppl_2):W293–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  81. 81.

    Rahnamaie-Tajadod R, Loke K-K, Goh H-H, Mohd Noor N. Differential gene expression analysis in Polygonum minus leaf upon 24 hours of methyl jasmonate elicitation. Front Plant Sci. 2017;8:109.

    PubMed  PubMed Central  Article  Google Scholar 

  82. 82.

    Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–61.

    CAS  Article  Google Scholar 

  83. 83.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C (T)) method. Methods. 2001;25:402–8.

    CAS  Article  Google Scholar 

  84. 84.

    Jiang H, Lei R, Ding S-W, Zhu S. Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinformatics. 2014;15(1):182.

    PubMed  PubMed Central  Article  Google Scholar 

  85. 85.

    Jung D, Kim B, Freishtat RJ, Giri M, Hoffman E, Seo J. miRTarVis: an interactive visual analysis tool for microRNA-mRNA expression profile data. BMC Proc. 2015;9(Suppl 6 Proceedings of the 5th Symposium on Biological Data):S2.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

Download references


We would like to thank Centre for Research and Instrumentation (CRIM) Universiti Kebangsaan Malaysia for providing SPME GC-MS analysis. We also want to thank Malaysian Genome Institute (MGI) and Universiti Putra Malaysia (UPM) for their Bioanalyzer service.


The author AS is sponsored by MyBrain 15 (MyPhD) scholarship. This study on miRNA is funded by Dana Impak Perdana UKM (DIP-2015-018) from Universiti Kebangsaan Malaysia (UKM), Malaysia. The funding body provide us the space and lab facilities to carry out the experiments. Thus, allowed the data to be analyzed and interpreted by authors as well as preparing the manuscript. This also included recruiting the graduate student who perform the experimental work.

Author information




AS designed and performed the laboratory work and interpreted the results. JJ carried out the bioinformatics and degradome analysis. AM and II advised the experimental design. RT and MN worked on the transcriptomic library. The first draft of this manuscript was written by AS and revised by MS and II. The revised version were made by AS, MS and II. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ismanizan Ismail.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

SPME GC-MS data from the inoculation test. These data show the effect of F. oxysporum treatment on terpenoid content in P. minor plants at 0, 3, 6, and 9 dpi. (CSV 1 kb)

Additional file 2:

Measurement of RNA integrity for each sample. An RNA integrity check was carried out prior to small RNA library construction. (PPTX 148 kb)

Additional file 3:

Spearman coefficients for miRNA-metabolite and mRNA-metabolite. Positive values indicate positive correlations, whereas negative values indicate negative correlations. (XLSX 8 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Samad, A.F.A., Rahnamaie-Tajadod, R., Sajad, M. et al. Regulation of terpenoid biosynthesis by miRNA in Persicaria minor induced by Fusarium oxysporum. BMC Genomics 20, 586 (2019).

Download citation


  • Terpenoid biosynthesis
  • miRNA
  • P. minor
  • F. Oxysporum
  • Deep sequencing
  • Post-transcriptional regulation