Advances in understanding epigenetic regulation of plant trichome development: a comprehensive review

Abstract

Plant growth and development are controlled by a complex gene regulatory network, which is currently a focal point of research. It has been established that epigenetic factors play a crucial role in plant growth. Trichomes, specialized appendages that arise from epidermal cells, are of great significance in plant growth and development. As a model system for studying plant development, trichomes possess both commercial and research value. Epigenetic regulation has only recently been implicated in the development of trichomes in a limited number of studies, and microRNA-mediated post-transcriptional regulation appears to dominate in this context. In light of this, we have conducted a review that explores the interplay between epigenetic regulations and the formation of plant trichomes, building upon existing knowledge of hormones and transcription factors in trichome development. Through this review, we aim to deepen our understanding of the regulatory mechanisms underlying trichome formation and shed light on future avenues of research in the field of epigenetics as it pertains to epidermal hair growth.

Acknowledgements

This research was supported by the National Natural Science Foundation of China (31830080, 32020103014) and the Construction of Beijing Science and Technology Innovation and Service Capacity in Top Subjects (CEFF-PXM2019 014207 000032). The authors are grateful to the members of Ren’s laboratory for their technical assistance and discussions.

Data availability

All data supporting the findings of this review are available within the article.

Conflict of interest

The authors declare that they have no conflict of interest.

Supplementary data

Supplementary data is available at Horticulture Research online.

References

[1] Ma X, Li K, Wang Z et al. Research progress in regulation model in different types of plant trichome. Sheng Wu Gong Cheng Xue Bao. 2020; 36: 2051-65 In Chinese

[2] Feng Z, Bartholomew ES, Liu Z et al. Glandular trichomes: new focus on horticultural crops. Hortic Res. 2021; 8: 158

[3] Andrade MC, Da Silva AA, Neiva IP et al. Inheritance of type IV glandular trichome density and its association with whitefly resistance from Solanum galapagense accession LA1401 . Euphytica. 2017; 213: 52

[4] Atalay Z, Celep F, Bara F et al. Systematic significance of anatomy and trichome morphology in Lamium (Lamioideae; Lamiaceae). Flora Morphol Distrib Funct Ecol Plants. 2016; 225: 60-75

[5] Szymanski DB, Jilk RA, Pollock SM et al. Control of GL2 expression in Arabidopsis leaves and trichomes . Development. 1998; 125: 1161-71

[6] Basra AS, Malik CP . Development of the cotton fiber. Int Rev Cytol. 1984; 89: 65-113

[7] Zhao JL, Pan JS, Guan Y et al. Transcriptome analysis in Cucumis sativus identifies genes involved in multicellular trichome development . Genomics. 2015; 105: 296-303

[8] Zhao JL, Wang YL, Yao DQ et al. Transcriptome profiling of trichome-less reveals genes associated with multicellular trichome development in Cucumis sativus. Mol Gen Genomics. 2015; 290: 2007-18

[9] Xue S, Dong M, Liu X et al. Classification of fruit trichomes in cucumber and effects of plant hormones on type II fruit trichome development. Planta. 2019; 249: 407-16

[10] Chen C, Liu M, Jiang L et al. Transcriptome profiling reveals roles of meristem regulators and polarity genes during fruit trichome development in cucumber (Cucumis sativus L.). J Exp Bot. 2014; 65: 4943-58

[11] Chang J, Xu Z, Li M et al. Spatiotemporal cytoskeleton organizations determine morphogenesis of multicellular trichomes in tomato. PLoS Genet. 2019; 15: e1008438

[12] Yang C, Ye Z . Trichomes as models for studying plant cell differentiation. Cell Mol Life Sci. 2013; 70: 1937-48

[13] Wang Z, Yang Z, Li F . Updates on molecular mechanisms in the development of branched trichome in Arabidopsis and nonbranched in cotton . Plant Biotechnol J. 2019; 17: 1706-22

[14] Schellmann S, Hülskamp M . Epidermal differentiation: trichomes in Arabidopsis as a model system . Int J Dev Biol. 2005; 49: 579-84

[15] Dong M, Xue S, Bartholomew ES et al. Transcriptomic and functional analysis provides molecular insights into multicellular trichome development. Plant Physiol. 2022; 189: 301-14

[16] Huchelmann A, Boutry M, Hachez C . Plant glandular trichomes: natural cell factories of high biotechnological interest. Plant Physiol. 2017; 175: 6-22

[17] Wang Y, Zeng J, Xia X et al. Comparative analysis of leaf trichomes, epidermal wax and defense enzymes activities in response to Puccinia horiana in Chrysanthemum and Ajania species . Hortic Plant J. 2020; 6: 191-8

[18] Chalvin C, Drevensek S, Dron M et al. Genetic control of glandular trichome development. Trends Plant Sci. 2020; 25: 477-87

[19] Traw MB, Bergelson J . Interactive effects of jasmonic acid, salicylic acid, and gibberellin on induction of trichomes in Arabidopsis. Plant Physiol. 2003; 133: 1367-75

[20] Inthima P, Nakano M, Otani M et al. Overexpression of the gibberellin 20-oxidase gene from Torenia fournieri resulted in modified trichome formation and terpenoid metabolites of Artemisia annua L. Plant Cell Tissue Organ Cult. 2017; 129: 223-36

[21] Liu X, Bartholomew E, Cai Y et al. Trichome-related mutants provide a new perspective on multicellular trichome initiation and development in cucumber (Cucumis sativus L). Front Plant Sci. 2016; 7: 1187

[22] Pikaard CS, Mittelsten SO . Epigenetic regulation in plants. Cold Spring Harb Perspect Biol. 2014; 6: a019315

[23] Tang D, Gallusci P, Lang Z . Fruit development and epigenetic modifications. New Phytol. 2020; 228: 839-44

[24] Liang Z, Riaz A, Chachar S et al. Epigenetic modifications of mRNA and DNA in plants. Mol Plant. 2020; 13: 14-30

[25] Yue J, Wei Y, Zhao M . The reversible methylation of m6A is involved in plant virus infection. Biology (Basel). 2022; 11: 271

[26] Růžička K, Zhang M, Campilho A et al. Identification of factors required for m ⁶A mRNA methylation in Arabidopsis reveals a role for the conserved E3 ubiquitin ligase HAKAI . New Phytol. 2017; 215: 157-72

[27] Shen L, Liang Z, Gu X et al. N(6)-Methyladenosine RNA modification regulates shoot stem cell fate in Arabidopsis. Dev Cell. 2016; 38: 186-200

[28] Vespa L, Vachon G, Berger F et al. The immunophilin-interacting protein AtFIP37 from Arabidopsis is essential for plant development and is involved in trichome endoreduplication . Plant Physiol. 2004; 134: 1283-92

[29] Zhou L, Tian S, Qin G . RNA methylomes reveal the m ⁶A-mediated regulation of DNA demethylase gene SlDML2 in tomato fruit ripening . Genome Biol. 2019; 20: 156

[30] Martínez-Pérez M, Aparicio F, López-Gresa MP et al. Arabidopsis m ⁶A demethylase activity modulates viral infection of a plant virus and the m ⁶A abundance in its genomic RNAs . Proc Natl Acad Sci USA. 2017; 114: 10755-60

[31] Miao Z, Zhang T, Qi Y et al. Evolution of the RNA N ⁶ -methyladenosine methylome mediated by genomic duplication . Plant Physiol. 2020; 182: 345-60

[32] Cui X, Liang Z, Shen L et al. 5-Methylcytosine RNA methylation in Arabidopsis thaliana. Mol Plant. 2017; 10: 1387-99

[33] Tang Y, Gao CC, Gao Y et al. OsNSUN2-mediated 5-methylcytosine mRNA modification enhances rice adaptation to high temperature. Dev Cell. 2020; 53: 272-286.e7

[34] Zhang P, Wu W, Chen Q et al. Non-coding RNAs and their integrated networks. J Integr Bioinform. 2019; 16: 20190027

[35] Carrington JC, Ambros V . Role of microRNAs in plant and animal development. Science. 2003; 301: 336-8

[36] Bartel DP . MicroRNAs: target recognition and regulatory functions. Cell. 2009; 136: 215-33

[37] Laufs P, Peaucelle A, Morin H et al. MicroRNA regulation of the CUC genes is required for boundary size control in Arabidopsis meristems . Development. 2004; 131: 4311-22

[38] Mallory AC, Dugas DV, Bartel DP et al. MicroRNA regulation of NAC-domain targets is required for proper formation and separation of adjacent embryonic, vegetative, and floral organs. Curr Biol. 2004; 14: 1035-46

[39] Guo HS, Xie Q, Fei JF et al. MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for Arabidopsis lateral root development . Plant Cell. 2005; 17: 1376-86

[40] Zhang X, Zou Z, Zhang J et al. Over-expression of sly-miR156a in tomato results in multiple vegetative and reproductive trait alterations and partial phenocopy of the sft mutant . FEBS Lett. 2011; 585: 435-9

[41] Sunkar R, Zhu JK . Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell. 2004; 16: 2001-19

[42] Zhang X, Bao Y, Shan D et al. Magnaporthe oryzae induces the expression of a microRNA to suppress the immune response in rice . Plant Physiol. 2018; 177: 352-68

[43] Hülskamp M, Mis´ra S, Jürgens G . Genetic dissection of trichome cell development in Arabidopsis. Cell. 1994; 76: 555-66

[44] Oppenheimer DG, Herman PL, Sivakumaran S et al. A myb gene required for leaf trichome differentiation in Arabidopsis is expressed in stipules . Cell. 1991; 67: 483-93

[45] Schellmann S, Schnittger A, Kirik V et al. TRIPTYCHON and CAPRICE mediate lateral inhibition during trichome and root hair patterning in Arabidopsis. EMBO J. 2002; 21: 5036-46

[46] Han G, Li Y, Qiao Z et al. Advances in the regulation of epidermal cell development by C2H2 zinc finger proteins in plants. Front Plant Sci. 2021; 12: 754512

[47] Sun L, Zhang A, Zhou Z et al. GLABROUS INFLORESCENCE STEMS3 (GIS3) regulates trichome initiation and development in Arabidopsis. New Phytol. 2015; 206: 220-30

[48] Gan Y, Liu C, Yu H et al. Integration of cytokinin and gibberellin signaling by Arabidopsis transcription factors GIS, ZFP8 and GIS2 in the regulation of epidermal cell fate . Development. 2007; 134: 2073-81

[49] Matías-Hernández L, Aguilar-Jaramillo AE, Osnato M et al. TEMPRANILLO reveals the mesophyll as crucial for epidermal trichome formation. Plant Physiol. 2016; 170: 1624-39

[50] Qi T, Song S, Ren Q et al. The jasmonate-ZIM-domain proteins interact with the WD-repeat/bHLH/MYB complexes to regulate jasmonate-mediated anthocyanin accumulation and trichome initiation in Arabidopsis thaliana. Plant Cell. 2011; 23: 1795-814

[51] Vadde BVL, Challa KR, Nath U . The TCP4 transcription factor regulates trichome cell differentiation by directly activating GLABROUS INFLORESCENCE STEMS in Arabidopsis thaliana. Plant J. 2018; 93: 259-69

[52] Yu N, Cai WJ, Wang S et al. Temporal control of trichome distribution by microRNA156-targeted SPL genes in Arabidopsis thaliana. Plant Cell. 2010; 22: 2322-35

[53] Wu G, Park MY, Conway SR et al. The sequential action of miR156 and miR172 regulates developmental timing in Arabidopsis. Cell. 2009; 138: 750-9

[54] Wang L, Zhou CM, Mai YX et al. A spatiotemporally regulated transcriptional complex underlies heteroblastic development of leaf hairs in Arabidopsis thaliana. EMBO J. 2019; 38: e100063

[55] Sorin C, Declerck M, Christ A et al. A miR169 isoform regulates specific NF-YA targets and root architecture in Arabidopsis. New Phytol. 2014; 202: 1197-211

[56] Wei LH, Song P, Wang Y et al. The m ⁶A reader ECT2 controls trichome morphology by affecting mRNA stability in Arabidopsis. Plant Cell. 2018; 30: 968-85

[57] Kotak J, Saisana M, Gegas V et al. The histone acetyltransferase GCN5 and the transcriptional coactivator ADA2b affect leaf development and trichome morphogenesis in Arabidopsis. Planta. 2018; 248: 613-28

[58] Wang T, Jia Q, Wang W et al. GCN5 modulates trichome initiation in Arabidopsis by manipulating histone acetylation of core trichome initiation regulator genes . Plant Cell Rep. 2019; 38: 755-65

[59] Patra B, Pattanaik S, Yuan L . Ubiquitin protein ligase 3 mediates the proteasomal degradation of GLABROUS 3 and ENHANCER OF GLABROUS 3, regulators of trichome development and flavonoid biosynthesis in Arabidopsis. Plant J. 2013; 74: 435-47

[60] Exner V, Gruissem W, Hennig L . Control of trichome branching by chromatin assembly factor-1. BMC Plant Biol. 2008; 8: 54

[61] Liu B, Zhu Y, Zhang T . The R3-MYB gene GhCPC negatively regulates cotton fiber elongation . PLoS One. 2015; 10: e0116272

[62] Guan XY, Li QJ, Shan CM et al. The HD-zip IV gene GaHOX1 from cotton is a functional homologue of the Arabidopsis GLABRA2. Physiol Plant. 2008; 134: 174-82

[63] Shan CM, Shangguan XX, Zhao B et al. Control of cotton fibre elongation by a homeodomain transcription factor GhHOX3. Nat Commun. 2014; 5: 5519

[64] Zhang J, Huang GQ, Zou D et al. The cotton (Gossypium hirsutum) NAC transcription factor (FSN1) as a positive regulator participates in controlling secondary cell wall biosynthesis and modification of fibers . New Phytol. 2018; 217: 625-40

[65] Huang J, Guo Y, Sun Q et al. Genome-wide identification of R2R3-MYB transcription factors regulating secondary cell wall thickening in cotton fiber development. Plant Cell Physiol. 2019; 60: 687-701

[66] Gong SY, Huang GQ, Sun X et al. Cotton KNL1, encoding a class II KNOX transcription factor, is involved in regulation of fiber development . J Exp Bot. 2014; 65: 4133-47

[67] Shi YH, Zhu SW, Mao XZ et al. Transcriptome profiling, molecular biological, and physiological studies reveal a major role for ethylene in cotton fiber cell elongation. Plant Cell. 2006; 18: 651-64

[68] Hu H, He X, Tu L et al. GhJAZ2 negatively regulates cotton fiber initiation by interacting with the R2R3-MYB transcription factor GhMYB25-like. Plant J. 2016; 88: 921-35

[69] Zhou Y, Zhang ZT, Li M et al. Cotton (Gossypium hirsutum) 14-3-3 proteins participate in regulation of fiber initiation and elongation by modulating brassinosteroid signaling . Plant Biotechnol J. 2015; 13: 269-80

[70] Sun Y, Veerabomma S, Abdel-Mageed HA et al. Brassinosteroid regulates fiber development on cultured cotton ovules. Plant Cell Physiol. 2005; 46: 1384-91

[71] Sun Y, Fokar M, Asami T et al. Characterization of the Brassinosteroid insensitive 1 genes of cotton . Plant Mol Biol. 2004; 54: 221-32

[72] Guan X, Pang M, Nah G et al. miR828 and miR858 regulate homoeologous MYB2 gene functions in Arabidopsis trichome and cotton fiber development . Nat Commun. 2014; 5: 3050

[73] Xie F, Wang Q, Sun R et al. Deep sequencing reveals important roles of microRNAs in response to drought and salinity stress in cotton. J Exp Bot. 2015; 66: 789-804

[74] Zhao T, Xu X, Wang M et al. Identification and profiling of upland cotton microRNAs at fiber initiation stage under exogenous IAA application. BMC Genomics. 2019; 20: 421

[75] Wang ZM, Xue W, Dong CJ et al. A comparative miRNAome analysis reveals seven fiber initiation-related and 36 novel miRNAs in developing cotton ovules. Mol Plant. 2012; 5: 889-900

[76] Sun R, Li C, Zhang J et al. Differential expression of microRNAs during fiber development between fuzzless-lintless mutant and its wild-type allotetraploid cotton. Sci Rep. 2017; 7: 3

[77] Xie F, Jones DC, Wang Q et al. Small RNA sequencing identifies miRNA roles in ovule and fiber development. Plant Biotechnol J. 2015; 13: 355-69

[78] Xue W, Wang Z, Du M et al. Genome-wide analysis of small RNAs reveals eight fiber elongation-related and 257 novel microRNAs in elongating cotton fiber cells. BMC Genomics. 2013; 14: 629

[79] Zhang B, Wang Q, Wang K et al. Identification of cotton microRNAs and their targets. Gene. 2007; 397: 26-37

[80] Kumar V, Singh B, Singh SK et al. Role of GhHDA5 in H3K9 deacetylation and fiber initiation in Gossypium hirsutum. Plant J. 2018; 95: 1069-83

[81] Wan Q, Guan X, Yang N et al. Small interfering RNAs from bidirectional transcripts of GhMML3 A12 regulate cotton fiber development . New Phytol. 2016; 210: 1298-310

[82] Wang M, Wang P, Tu L et al. Multi-omics maps of cotton fibre reveal epigenetic basis for staged single-cell differentiation. Nucleic Acids Res. 2016; 44: 4067-79

[83] Feng H, Li X, Chen H et al. GhHUB2, a ubiquitin ligase, is involved in cotton fiber development via the ubiquitin-26S proteasome pathway . J Exp Bot. 2018; 69: 5059-75

[84] Wang Y, Li Y, Gong SY et al. GhKNL1 controls fiber elongation and secondary cell wall synthesis by repressing its downstream genes in cotton (Gossypium hirsutum). J Integr Plant Biol. 2022; 64: 39-55

[85] Pan Y, Bo K, Cheng Z et al. The loss-of-function GLABROUS 3 mutation in cucumber is due to LTR-retrotransposon insertion in a class IV HD-ZIP transcription factor gene CsGL3 that is epistatic over CsGL1. BMC Plant Biol. 2015; 15: 302

[86] Cui JY, Miao H, Ding LH et al. A new Glabrous gene (csgl3) identified in trichome development in cucumber (Cucumis sativus L.). PLoS One. 2016; 11: e0148422

[87] Wang YL, Nie JT, Chen HM et al. Identification and mapping of Tril, a homeodomain-leucine zipper gene involved in multicellular trichome initiation in Cucumis sativus. Theor Appl Genet. 2016; 129: 305-16

[88] Li Q, Cao C, Zhang C et al. The identification of Cucumis sativus Glabrous 1 (CsGL1) required for the formation of trichomes uncovers a novel function for the homeodomain-leucine zipper I gene . J Exp Bot. 2015; 66: 2515-26

[89] Zhao JL, Pan JS, Guan Y et al. Micro-trichome as a class I homeodomain-leucine zipper gene regulates multicellular trichome development in Cucumis sativus. J Integr Plant Biol. 2015; 57: 925-35

[90] Yang S, Cai Y, Liu X et al. A CsMYB6-CsTRY module regulates fruit trichome initiation in cucumbers . J Exp Bot. 2018; 69: 1887-902

[91] Chen C, Yin S, Liu X et al. The WD-repeat protein CsTTG1 regulates fruit wart formation through interaction with the homeodomain-leucine zipper I protein Mict. Plant Physiol. 2016; 171: 1156-68

[92] Zhang W, He H, Guan Y et al. Identification and mapping of molecular markers linked to the tuberculate fruit gene in the cucumber (Cucumis sativus L.). Theor Appl Genet. 2010; 120: 645-54

[93] Yang X, Zhang W, He H et al. Tuberculate fruit gene Tu encodes a C2H2 zinc finger protein that is required for the warty fruit phenotype in cucumber (Cucumis sativus L.). Plant J. 2014; 78: 1034-46

[94] Wang Z, Wang L, Han L et al. HECATE2 acts with GLABROUS3 and Tu to boost cytokinin biosynthesis and regulate cucumber fruit wart formation. Plant Physiol. 2021; 187: 1619-35

[95] Sun H, Pang B, Yan J et al. Comprehensive analysis of cucumber gibberellin oxidase family genes and functional characterization of CsGA20ox1 in root development in Arabidopsis. Int J Mol Sci. 2018; 19: 3135

[96] Xie Q, Liu P, Shi L et al. Combined fine mapping, genetic diversity, and transcriptome profiling reveals that the auxin transporter gene ns plays an important role in cucumber fruit spine development . Theor Appl Genet. 2018; 131: 1239-52

[97] Yang S, Wen C, Liu B et al. A CsTu-TS1 regulatory module promotes fruit tubercule formation in cucumber . Plant Biotechnol J. 2019; 17: 289-301

[98] Liu X, Wang T, Bartholomew E et al. Comprehensive analysis of NAC transcription factors and their expression during fruit spine development in cucumber (Cucumis sativus L.). Hortic Res. 2018; 5: 31

[99] Zhang X, Lai Y, Zhang W et al. MicroRNAs and their targets in cucumber shoot apices in response to temperature and photoperiod. BMC Genomics. 2018; 19: 819

[100] Xu J, van Herwijnen ZO, Dräger DB et al. SlMYC1 regulates type VI glandular trichome formation and terpene biosynthesis in tomato glandular cells . Plant Cell. 2018; 30: 2988-3005

[101] Yang C, Li H, Zhang J et al. A regulatory gene induces trichome formation and embryo lethality in tomato. Proc Natl Acad Sci USA. 2011; 108: 11836-41

[102] Yang C, Li H, Zhang J et al. Fine-mapping of the woolly gene controlling multicellular trichome formation and embryonic development in tomato . Theor Appl Genet. 2011; 123: 625-33

[103] Hua B, Chang J, Wu M et al. Mediation of JA signalling in glandular trichomes by the woolly/SlMYC1 regulatory module improves pest resistance in tomato . Plant Biotechnol J. 2021; 19: 375-93

[104] Zheng F, Cui L, Li C et al. Hair interacts with SlZFP8-like to regulate the initiation and elongation of trichomes by modulating SlZFP6 expression in tomato . J Exp Bot. 2022; 73: 228-44

[105] Gao S, Gao Y, Xiong C et al. The tomato B-type cyclin gene, SlCycB2, plays key roles in reproductive organ development, trichome initiation, terpenoids biosynthesis and Prodenia litura defense . Plant Sci. 2017; 262: 103-14

[106] Ying S, Su M, Wu Y et al. Trichome regulator SlMIXTA-like directly manipulates primary metabolism in tomato fruit . Plant Biotechnol J. 2020; 18: 354-63

[107] Xie Q, Gao Y, Li J et al. The HD-zip IV transcription factor SlHDZIV8 controls multicellular trichome morphology by regulating the expression of Hairless-2. J Exp Bot. 2020; 71: 7132-45

[108] Yu X, Chen G, Tang B et al. The jasmonate ZIM-domain protein gene SlJAZ2 regulates plant morphology and accelerates flower initiation in Solanum lycopersicum plants . Plant Sci. 2018; 267: 65-73

[109] Hua B, Chang J, Han X et al. H and HL synergistically regulate jasmonate-triggered trichome formation in tomato. Hortic Res. 2022; 9: uhab080

[110] Chang J, Yu T, Yang Q et al. Hair, encoding a single C2H2 zinc-finger protein, regulates multicellular trichome formation in tomato . Plant J. 2018; 96: 90-102

[111] Thines B, Katsir L, Melotto M et al. JAZ repressor proteins are targets of the SCF(COI1) complex during jasmonate signaling. Nature. 2007; 448: 661-5

[112] Hua B, Chang J, Xu Z et al. HOMEODOMAIN PROTEIN8 mediates jasmonate-triggered trichome elongation in tomato . New Phytol. 2021; 230: 1063-77

[113] Chen Y, Su D, Li J et al. Overexpression of bHLH95, a basic helix-loop-helix transcription factor family member, impacts trichome formation via regulating gibberellin biosynthesis in tomato . J Exp Bot. 2020; 71: 3450-62

[114] Deng W, Yang Y, Ren Z et al. The tomato SlIAA15 is involved in trichome formation and axillary shoot development . New Phytol. 2012; 194: 379-90

[115] Zhang X, Yan F, Tang Y et al. Auxin response gene SlARF3 plays multiple roles in tomato development and is involved in the formation of epidermal cells and trichomes . Plant Cell Physiol. 2015; 56: 2110-24

[116] Yuan Y, Xu X, Luo Y et al. R2R3 MYB-dependent auxin signalling regulates trichome formation, and increased trichome density confers spider mite tolerance on tomato. Plant Biotechnol J. 2021; 19: 138-52

[117] Gong Z, Luo Y, Zhang W et al. A SlMYB75-centred transcriptional cascade regulates trichome formation and sesquiterpene accumulation in tomato . J Exp Bot. 2021; 72: 3806-20

[118] Liao X, Wang J, Zhu S et al. Transcriptomic and functional analyses uncover the regulatory role of lncRNA000170 in tomato multicellular trichome formation. Plant J. 2020; 104: 18-29

[119] Zilong L, Jingwei F, Rui Y et al. Effects of MiR319a on the growth and development of tomato. J Beijing Agric Univ. 2015; 30: 49-53

[120] Bao N, Lye K, Barton MK . MicroRNA binding sites in Arabidopsis class III HD-ZIP mRNAs are required for methylation of the template chromosome . Dev Cell. 2004; 7: 653-62

[121] Wu L, Zhou H, Zhang J et al. DNA methylation mediated by a microRNA pathway. Mol Cell. 2010; 38: 465-75