Abstract
MicroRNAs (miRNAs) play crucial roles in plant development and secondary metabolism through different modes of sequence-specific interaction with their targets. Artemisinin biosynthesis is extensively regulated by phytohormones. However, the function of phytohormone-responsive miRNAs in artemisinin biosynthesis remains enigmatic. Thus, we combined the analysis of transcriptomics, small RNAs, and the degradome to generate a comprehensive resource for identifying key miRNA-target circuits involved in the phytohormone-induced process of artemisinin biosynthesis in Artemisia annua. In total, 151 conserved and 52 novel miRNAs and their 4132 targets were determined. Based on the differential expression analysis, miR160 was selected as a potential miRNA involved in artemisinin synthesis. Overexpressing MIR160 significantly impaired glandular trichome formation and suppressed artemisinin biosynthesis in A. annua, while repressing its expression resulted in the opposite effect, indicating that miR160 negatively regulates glandular trichome development and artemisinin biosynthesis. RNA ligase-mediated 5' RACE and transient transformation assays showed that miR160 mediates the RNA cleavage of Auxin Response Factor 1 (ARF1) in A. annua. Furthermore, ARF1 was shown to increase artemisinin synthesis by activating AaDBR2 expression. Taken together, our results reveal the intrinsic link between the miR160-ARF1 module and artemisinin biosynthesis, and may expedite the innovation of metabolic engineering approaches for high and stable production of artemisinin in the future.