MicroRNAome Profile of Euphorbia kansui in Response to Methyl Jasmonate

miRNAs play vital regulatory roles in different plant developmental stages and in plant response to biotic and abiotic stresses. However, information is limited on the miRNA regulatory mechanism to methyl jasmonate (MeJA). In this study, we used the microRNAome profile to illustrate the relevant regulatory mechanisms of Euphorbia kansui in response to methyl jasmonate (MeJA) through Illumina RNA-Seq. As a result, we identified 875 miRNAs corresponding to 11,277 target mRNAs, among them, 168 known miRNA families representing 6019 target mRNAs sequences were obtained. 452 miRNA-mRNA pairs presented an anti-correlationship (Cor < −0.50 and p-value of correlation ≤ 0.05). The miRNA with a fold change ≥ 2 and a p (p-Value) < 0.05 in pairwise comparison were identified as significant differentially expressed miRNAs (DEMs). The DEMs in MeJA treatment of 0, 24, 36 and 48 h were compared by using Short Time Expression Miner (STEM) cluster and 4 significant gene profiles (p-value ≤ 0.02) were identified. Through the kyoto encyclopedia of genes and genomes (KEGG) pathway and gene ontology (GO) enrichment analysis on all miRNA targets, we identified 33 mRNAs in terpenoid biosynthesis, which were regulated by miRNAs under MeJA treatment, so the miRNA maybe involved in the response of E. kansui plant to exogenous MeJA and the results would provide very useful information on illustrating the regulatory mechanism of E. kansui and also provide an overall view of the miRNAs response to MeJA stress of a non-model plant.