Environmental metal pollution is a common problem threatening sustainable and safe crop production. Heavy metals (HMs) cause toxicity by targeting key molecules and life processes in plant cells. Plants counteract excess metals in the environment by enhancing defense responses, such as metal chelation, isolation to vacuoles, regulating metal intake through transporters, and strengthening antioxidant mechanisms. In recent years, microRNAs (miRNAs), as a small non-coding RNA, have become the central regulator of a variety of abiotic stresses, including HMs. With the introduction of the latest technologies such as next-generation sequencing (NGS), more and more miRNAs have been widely recognized in several plants due to their diverse roles. Metal-regulated miRNAs and their target genes are part of a complex regulatory network. Known miRNAs coordinate plant responses to metal stress through antioxidant functions, root growth, hormone signals, transcription factors (TF), and metal transporters. This article reviews the research progress of miRNAs in the stress response of plants to the accumulation of HMs, such as Cu, Cd, Hg, Cr, and Al, and the toxicity of heavy metal ions.

Ilex verticillata is not only an excellent ornamental tree species for courtyards, but it is also a popular bonsai tree. ‘Oosterwijk’ and ‘Red sprite’ are two varieties of Ilex verticillata. The former has a long stem with few branches, while the latter has a short stem. In order to explain the stem growth differences between the two cultivars ‘Oosterwijk’ and ‘Red sprite’, determination of the microstructure, transcriptome sequence and IAA content was carried out. The results showed that the xylem thickness, vessel area and vessel number of ‘Oosterwijk’ were larger than in ‘Red sprite’. In addition, our analysis revealed that the differentially expressed genes which were enriched in phenylpropanoid biosynthesis; phenylalanine metabolism and phenylalanine, tyrosine and tryptophan biosynthesis in the black and tan modules of the two varieties. We found that AST, HCT and bHLH 94 may be key genes in the formation of shoot difference. Moreover, we found that the IAA content and auxin-related DEGs GH3.6, GH3, ATRP5, IAA27, SAUR36-like, GH3.6-like and AIP 10A5-like may play important roles in the formation of shoot differences. In summary, these results indicated that stem growth variations of ‘Oosterwijk’ and ‘Red sprite’ were associated with DEGs related to phenylpropanoid biosynthesis, phenylalanine metabolism and phenylalanine, tyrosine and tryptophan biosynthesis, as well as auxin content and DEGs related to the auxin signaling pathway.