Functional characterization of the tomato HAIRPLUS gene reveals the implication of the epigenome in the control of glandular trichome formation

Abstract

Trichomes are specialised epidermal cells developed in the aerial surface of almost every terrestrial plant. These structures form physical barriers, which combined with their capability of synthesis of complex molecules, prevent plagues from spreading and confer trichomes a key role in the defence against herbivores. In this work, the tomato gene HAIRPLUS (HAP) that controls glandular trichome density in tomato plants was characterised. HAP belongs to a group of proteins involved in histone tail modifications although some also bind methylated DNA. HAP loss of function promotes epigenomic modifications in the tomato genome reflected in numerous differentially methylated cytosines and causes transcriptomic changes in hap mutant plants. Taken together, these findings demonstrate that HAP links epigenome remodelling with multicellular glandular trichome development and reveal that HAP is a valuable genomic tool for pest resistance in tomato breeding.

References

[1] Chang J, Yu T, Gao S et al. Fine mapping of the dialytic gene that controls multicellular trichome formation and stamen development in tomato. Theor Appl Genet. 2016; 129: 1531-9.

[2] Kang JH, Shi F, Jones AD et al. Distortion of trichome morphology by the hairless mutation of tomato affects leaf surface chemistry. J Exp Bot. 2010; 61: 1053-64.

[3] Martin C, Glover BJ . Functional aspects of cell patterning in aerial epidermis. Curr Opin Plant Biol. 2007; 10: 70-82.

[4] Hülskamp M, Schnittger A, Folkers U . Pattern formation and cell differentiation: trichomes in Arabidopsis as a genetic model system. Int Rev Cytol. 1999; 186: 147-78.

[5] Marks MD, Feldmann KA . Trichome development in Arabidopsis thaliana. I. T-DNA tagging of the GLABROUS1 gene. Plant Cell. 1989; 1: 1043-50.

[6] Serna L, Martin C . Trichomes: different regulatory networks lead to convergent structures. Trends Plant Sci. 2006; 11: 274-80.

[7] Zhang F, Gonzalez A, Zhao M et al. A network of redundant bHLH proteins functions in all TTG1-dependent pathways of Arabidopsis. Development. 2003; 130: 4859-69.

[8] Esch JJ, Chen M, Sanders M et al. A contradictory GLABRA3 allele helps define gene interactions controlling trichome development in Arabidopsis. Development. 2003; 130: 5885-94.

[9] Glover BJ, Bunnewell S, Martin C . Convergent evolution within the genus Solanum: the specialised anther cone develops through alternative pathways. Gene. 2004; 331: 1-7.

[10] Payne T, Clement J, Arnold D et al. Heterologous myb genes distinct from GL1 enhance trichome production when overexpressed in Nicotiana tabacum. Development. 1999; 126: 671-82.

[11] Galdon-Armero J, Arce-Rodriguez L, Downie M et al. Scanning electron micrograph-based resource for identification of loci involved in epidermal development in tomato: elucidation of a new function for the Mixta-like transcription factor in leaves. Plant Cell. 2020; 32: 1414-33.

[12] Glover BJ, Perez-Rodriguez M, Martin C . Development of several epidermal cell types can be specified by the same MYB-related plant transcription factor. Development. 1998; 125: 3497-508.

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

[14] Luckwill LC . The Genus Lycopersicon: An Historical, Biological, And Taxonomic survery Of The Wild And Cultivated Tomatoes. Aberdeen: Aberdeen University Press; 1943.

[15] Vendemiatti E, Zsogon A, Silva GFFFE et al. Loss of type-IV glandular trichomes is a heterochronic trait in tomato and can be reverted by promoting juvenility. Plant Sci. 2017; 259: 35-47.

[16] Simmons AT, Gurr GM . Trichomes of Lycopersicon species and their hybrids: effects on pests and natural enemies. Agric Forest Ent. 2005; 7: 265-76.

[17] Kang JH, Campos ML, Zemelis-Durfee S et al. Molecular cloning of the tomato hairless gene implicates actin dynamics in trichome-mediated defense and mechanical properties of stem tissue. J Exp Bot. 2016; 67: 5313-24.

[18] 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.

[19] 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.

[20] Huang PC, Paddock EF . The time and site of the semidominant lethal action of “Wo” in Lycopersicon esculentum. Am J Bot. 1962; 49: 388-93.

[21] 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.

[22] 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.

[23] 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.

[24] Lachner M, Jenuwein T . The many faces of histone lysine methylation. Curr Opin Cell Biol. 2002; 14: 286-98.

[25] Marmorstein R . Structure of SET domain proteins: a new twist on histone methylation. Trends Biochem Sci. 2003; 28: 59-62.

[26] Capel C, del Carmen AF, Alba JM et al. Wide-genome QTL mapping of fruit quality traits in a tomato RIL population derived from the wild-relative species Solanum pimpinellifolium L. Theor Appl Genet. 2015; 128: 2019-35.

[27] Jacobs SA, Harp JM, Devarakonda S et al. The active site of the SET domain is constructed on a knot. Nat Struct Biol. 2002; 9: 833-83.

[28] Law JA, Jacobsen SE . Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat Rev Genet. 2010; 11: 204-20.

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

[30] Hu Y, Xie Q, Chua NH . The Arabidopsis auxin-inducible gene ARGOS controls lateral organ size. Plant Cell. 2003; 15: 1951-61.

[31] Shuai B, Reynaga-Pena CG, Springer PS . The lateral organ boundaries gene defines a novel, plant-specific gene family. Plant Physiol. 2002; 129: 747-61.

[32] Hur YS, Um J-H, Kim S et al. Arabidopsis thaliana homeobox 12 (ATHB 12), a homeodomain-leucine zipper protein, regulates leaf growth by promoting cell expansion and endoreduplication. New Phytol. 2015; 205: 316-28.

[33] Song L, Carol Huang S-S, Wise A et al. A transcription factor hierarchy defines an environmental stress response network. Science. 2016; 354: aag1550.

[34] Powell AL, Nguyen C, Hill T et al. Uniform ripening encodes a Golden 2-like transcription factor regulating tomato fruit chloroplast development. Science. 2012; 336: 1711-5.

[35] Vogt T, Jones P . Glycosyltransferases in plant natural product synthesis: characterization of a supergene family. Trends Plant Sci. 2000; 5: 380-6.

[36] Harris CJ, Scheibe M, Wongpalee SP et al. A DNA methylation reader complex that enhances gene transcription. Science. 2018; 362: 1182-6.

[37] Akalin A, Kormaksson M, Li S et al. MethylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 2012; 13: R87.

[38] Xiao X, Zhang J, Li T et al. A group of SUVH methyl-DNA binding proteins regulate expression of the DNA demethylase ROS1 in Arabidopsis. J Integr Plant Biol. 2019; 61: 110-9.

[39] Baumbusch LO, Thorstensen T, Krauss V et al. The Arabidopsis thaliana genome contains at least 29 active genes encoding SET domain proteins that can be assigned to four revolutionarily conserved classes. Nucleic Acids Res. 2001; 29: 4319-33.

[40] Satish M, Nivya MA, Abhishek S et al. Computational characterization of substrate and product specificities, and functionality of S-adenosylmethionine binding pocket in histone lysine methyltransferases from Arabidopsis, rice and maize. Proteins. 2018; 86: 21-34.

[41] Jackson JP, Lindroth AM, Cao X et al. Control of CpNpG DNA methylation by the KRYPTONITE histone H3 methyltransferase. Nature. 2002; 416: 556-60.

[42] Zhao QQ, Lin RN, Li L et al. A methylated-DNA-binding complex required for plant development mediates transcriptional activation of promoter methylated genes. J Integr Plant Biol. 2019; 61: 120-39.

[43] Lioznova AV, Khamis AM, Artemov AV et al. CpG traffic lights are markers of regulatory regions in human genome. BMC Genomics. 2019; 20: 102.

[44] Simmons AT, Gurr GM, McGrath D et al. Entrapment of Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) on glandular trichomes of Lycopersicon species. Aust J Entomol. 2004; 43: 196-200.

[45] Carter CD, Snyder JC . Mite responses in relation to trichomes of Lycopersicon esculentum x L. hirsutum F₂ hybrids . Euphytica. 1985; 34: 177-85.

[46] Salinas M, Capel C, Alba JM et al. Genetic mapping of two QTL from the wild tomato Solanum pimpinellifolium L. controlling resistance against two-spotted spider mite (Tetranychus urticae Koch). Theor Appl Genet. 2013; 126: 83-92.

[47] Lucatti AF, van Heusden AW, de Vos RC et al. Differences in insect resistance between tomato species endemic to the Galapagos Islands. BMC Evol Biol. 2013; 13: 175.

[48] Lozano R, Angosto T, Gomez P et al. Tomato flower abnormalities induced by low temperatures are associated with changes of expression of MADS-box genes. Plant Physiol. 1998; 117: 91-100.

[49] Langmead B, Salzberg SL . Fast gapped-read alignment with bowtie 2. Nat Methods. 2012; 9: 357-9.

[50] DePristo MA, Banks E, Poplin R et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011; 43: 491-8.

[51] Danecek P, Auton A, Abecasis G et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8.

[52] Li H, Bob H, Wysoker A et al. The sequence alignment/map (SAM) format and SAMtools. Bioinformatics. 2009; 25: 2078-9.

[53] R Development Core Team. R Development Core Team : R: A Language And Environment For Statistical Computing. Vienna: R Foundation for Statistical Computing; 2011.

[54] Pérez-Martín F, Yuste-Lisbona FJ, Pineda B et al. A collection of enhancer trap insertional mutants for functional genomics in tomato. Plant Biotechnol J. 2017; 15: 1439-52.

[55] Vazquez-Vilar M, Bernabe-Orts JM, Del-Carmen AF et al. A modular toolbox for gRNA-Cas9 genome engineering in plants based on the GoldenBraid standard. Plant Methods. 2016; 12: 10.

[56] Oliveros JC, Franch M, Tabas-Madrid D et al. Breaking-Cas-interactive design of guide RNAs for CRISPR-Cas experiments for ENSEMBL genomes. Nucleic Acids Res. 2016; 44: W267-71.

[57] Yuste-Lisbona FJ, Fernandez-Lozano A, Pineda B et al. ENO regulates tomato fruit size through the floral meristem development network. Proc Natl Acad Sci. 2020; 117: 8187-95.

[58] Capel C, Albaladejo I, Egea I et al. The res (restored cell structure by salinity) tomato mutant reveals the role of the DEAD-box RNA helicase SlDEAD39 in plant development and salt response. Plant Cell Environ. 2020; 43: 1722-39.

[59] Xu J, Zhou S, Gong X et al. Single-base methylome analysis reveals dynamic epigenomic differences associated with water deficit in apple. Plant Biotechnol J. 2018; 16: 672-87.