Genome-wide Identification and Expression Analysis of DIR Gene Family in Tomato (Solanum lycopersicum)
SUN Hui-Ru1,2,*, ZHANG Jia-Ning1, REN Min1, WANG Yan-Feng1,2,*
1 College of Life Sciences, Yan'an University, Yan'an 716000, China; 2 Shaanxi Key Laboratory of Chinese Jujube, Yan'an University, Yan'an 716000, China
Abstract Dirigent genes (DIRs) are important regulator in the biosynthesis of plant lignin and lignans, and play pivotal roles in plant growth and development, and stress response. In order to explore the function of SlDIR family genes in tomato (Solanum lycopersicum), bioinformatics analysis of SlDIR members were carried out and expression characteristics of them in different organs, various stresses and abscisic acid (ABA) treatments were detected by released RNA-seq in TFGD website and qPCR. The results showed that 31 SlDIR genes were identified in tomato, which were unevenly distributed on 11 chromosomes (except Chr03) with 3 tandem duplication clusters. Most of SlDIR genes had no intron. The phylogenetic analysis displayed that the SlDIR family mumbers were classified into 3 subfamilies, the member numbers of which were various. All of SlDIR proteins had dirigent domain, and most of them also contained signal peptide. The distribution and arrangement of the motifs in SlDIRs from the same subfamily were similar. In addition, plant hormone and stress response elements were widely distributed in the promoter regions of SlDIRs. The expression characteristics of SlDIRs in different organs and different resistant materials under Pseudomonas syringae pv. tomato DC3000 (Pst DC3000) treatment were further analyzed. It was found that most SlDIRs were predominance expressed in roots, and the expressions of SlDIR5, SlDIR9, SlDIR13, SlDIR15 and SlDIR27 were significantly higher in the early stages than that in the later stages of tomato fruit development. The significantly different expressions of 11 SlDIRs were displayed between resistant and sensitive tomato varieties under Pst DC3000 treatment. Under ABA treatment for 12 h, the expressions of 7 SlDIRs were significantly up-regulated. Salt stress significantly induced the expression of SlDIR27. Meanwhile, SlDIR13 and SlDIR27 showed up-regulated expressions significantly under drought treatment. The results of this study provide a reference for further study on the biological functions of SlDIR family genes.
SUN Hui-Ru,ZHANG Jia-Ning,REN Min, et al. Genome-wide Identification and Expression Analysis of DIR Gene Family in Tomato (Solanum lycopersicum)[J]. 农业生物技术学报, 2023, 31(1): 36-49.
[1] 陈家璐, 张智俊, 刘笑雨, 等. 2019. 毛竹Dirigent基因家族的全基因组鉴定与分析[J]. 植物生理学报, 55(9): 1406-1417. (Chen J L, Zhang Z J, Liu X Y, et al.2019. Genome-wide identification and analysis of Dirigent gene family in moso bamboo (Phyllostachys edulis)[J]. Plant Physiology Journal, 55(9): 1406-1417.) [2] 刘辉, 王涛涛, 张俊红,等. 2015. 番茄NAC转录因子SlNAC80的克隆及表达分析[J]. 华北农学报, 30(1):77-83. (Li H, Wang T T, Zhang J H, et al.2015. Clonging and expression analysis of a NAC transcription factor SlNAC80 in tomato[J]. Acta Agriculturae Boreali-Sinica, 30(1): 77-83.) [3] 翟明明, 刘娜, 徐任园, 等. 2022. 马铃薯VQ基因家族鉴定与表达分析[J]. 农业生物技术学报, 30(1): 25-37. (Zhai M M, Liu N, Xu R Y, et al.2022. Identification and expression analysis of VQ gene family in Solanum tuberosum[J]. Journal of Agricultural Biotechnology, 30(1): 25-37.) [4] Armenteros J J A, Tsirigos K D, Sønderby C K, et al.2019. SignalP 5.0 improves signal peptide predictions using deep neural networks[J]. Nature Biotechnology, 37(4): 420-423. [5] Aoki K, Yano K, Suzuki A, et al.2010. Large-scale analysis of full-length cDNAs from the tomato (Solanum lycopersicum) cultivar Micro-Tom, a reference system for the Solanaceae genomics[J]. BMC Genomics, 11: 210-225. [6] Arasan S K T, Park J I, Ahmed N U, et al.2013. Characterization and expression analysis of dirigent family genes related to stresses in Brassica[J]. Plant Physiology and Biochemistry, 67: 144-153. [7] Artimo P, Jonnalagedda M, Arnold K, et al.2012. ExPASy: SIB bioinformatics resource portal[J]. Nucleic Acids Research, 40(W1): W597-W603. [8] Bailey T L, Mikael B, Buske F A, et al.2009. MEME SUITE: Tools for motif discovery and searching[J]. Nucleic Acids Research, 37: W202-W208. [9] Berim A, Ebel R, Schneider B, et al.2008. UDP-glucose: (6-methoxy) podophyllotoxin 7-O-glucosyltransferase from suspension cultures of Linum nodiflorum[J]. Phytochemistry, 69(2): 374-381. [10] Burlat V, Mi K, Davin L B, et al.2001. Dirigent proteins and dirigent sites in lignifying tissues[J]. Phytochemistry, 57(6): 883-897. [11] Chen C J, Chen H, Zhang Y, et al.2020. TBtools: An integrative toolkit developed for interactive analyses of big biological data[J]. Molecular Plant, 13(8): 1194-1202. [12] Cheng X, Su X Q, Muhammad A, et al.2018. Molecular characterization, evolution, and expression profiling of the dirigent (DIR) family genes in chinese white pear (Pyrus bretschneideri)[J]. Frontiers in Genetics, 9: 136. [13] Davin L B, Jourdes M, Patten A M, et al.2008. Dissection of lignin macromolecular configuration and assembly: Comparison to related biochemical processes in allyl/propenyl phenol and lignan biosynthesis[J]. Cheminform, 25(6): 1015-1090. [14] Davin L B, Lewis N G.2005. Dirigent phenoxy radical coupling: Advances and challenges[J]. Current Opinion in Biotechnology, 16(4): 398-406. [15] Davin L B, Wang H B, Crowell A L, et al.1997. Stereoselective bimolecular phenoxy radical coupling by an auxiliary (dirigent) protein without an active center[J]. Science, 275(5298): 362-366. [16] Finn R D, Penelope C, Eberhardt R Y, et al.2016. The Pfam protein families database: Towards a more sustainable future[J]. Nucleic Acids Research, 40(D1): D279-D285. [17] Gu Z l, Cavalcanti A, Chen F C, et al.2002. Extent of gene duplication in the genomes of drosophila, nematode, and yeast[J]. Molecular Biology and Evolution, 19(3): 256-262. [18] Horton P, Park K J, Obayashi T, et al.2007. WoLF PSORT: Protein localization predictor[J]. Nucleic Acids Research, 35: 585-587. [19] Hu B, Jin J, Guo A Y, et al.2015. GSDS 2.0: An upgraded gene feature visualization server[J]. Bioinformatics, 31(8): 1296. [20] Jones J, Dangl J L.2006. The plant immune system[J]. Nature, 444(7117): 323-329. [21] Khan A, Li R J, Sun J T, et al.2018. Genome-wide analysis of dirigent gene family in pepper (Capsicum annuum L.) and characterization of CaDIR7 in biotic and abiotic stresses[J]. Scientific Reports, 8(1): 5500. [22] Kunkel B N, Brooks D M.2002. Cross talk between signaling pathways in pathogen defense[J]. Current Opinion in Plant Biology, 5(4): 325-331. [23] Li L L, Sun W B, Zhou P J, et al.2021. Genome-wide characterization of dirigent proteins in populus: Gene expression variation and expression pattern in response to Marssonina brunnea and phytohormones[J]. Forests, 12(4): 507. [24] Li N H, Zhao M, Liu T F, et al.2017. A novel soybean dirigent gene GmDIR22 contributes to promotion of lignan biosynthesis and enhances resistance to Phytophthora sojae[J]. Frontiers in Plant Science, 8: 1185. [25] Liao Y R, Liu S B, Jiang Y Y, et al.2016. Genome-wide analysis and environmental response profiling of dirigent family genes in rice (Oryza sativa)[J]. Genes & Genomics, 39(1): 47-62. [26] Librado P, Rozas R.2009. DnaSP ver. 5: A software for comprehensive analyses of DNA polymorphism data[J]. Bioinformatics, 25(11): 1451-1452. [27] Liu Z W, Wang X F, Sun Z W, et al.2021. Evolution, expression and functional analysis of cultivated allotetraploid cotton DIR genes[J]. BMC Plant Biology, 21(1): 89. [28] Ma Q H, Liu Y C.2014. TaDIR13, a dirigent protein from wheat, promotes lignan biosynthesis and enhances pathogen resistance[J]. Plant Molecular Biology Reporter, 33(1): 143-152. [29] Ma X F, Xu W Y, Liu T, et al.2021. Functional characterization of soybean (Glycine max) DIRIGENT genes reveals an important role of GmDIR27 in the regulation of pod dehiscence[J]. Genomics, 113(1 Pt 2): 979-990. [30] Mitton F M, Pinedo M L, de la Canal L.2009. Phloem sap of tomato plants contains a DIR1 putative ortholog[J]. Journal of Plant Physiology, 166(5): 543-7. [31] Moura J C M S, Bonine C A V, Viana J D F, et al.2010. Abiotic and biotic stresses and changes in the lignin content and composition in plants[J]. Journal of Integrative Plant Biology, 52(4): 360-376. [32] Nobile P M, Bottcher A, Mayer J L S, et al.2017. Identification, classification and transcriptional profiles of dirigent domain-containing proteins in sugarcane[J]. Molecular Genetics and Genomics, 292(6): 1323-1340. [33] Paniagua C, Bilkova A, Jackson P, et al.2017. Dirigent proteins in plants: Modulating cell wall metabolism during abiotic and biotic stress exposure[J]. Journal of Experimental Botany, 68(13): 3287-3301. [34] Qin M Y, Luo W B, Zheng Y, et al.2019. Genome-wide identification and expression analysis of the PHD-finger gene family in Solanum tuberosum[J]. PLOS ONE, 14(12): e0226964. [35] Ralph S, Park J Y, Bohlmann J, et al.2006. Dirigent proteins in conifer defense: Gene discovery, phylogeny, and differential wound- and insect-induced expression of a family of DIR and DIR-like genes in spruce (Picea spp.)[J]. Plant Molecular Biology, 60(1): 21-40. [36] Ralph S G, Jancsik S, Bohlmann J.2007. Dirigent proteins in conifer defense Ⅱ: Extended gene discovery, phylogeny, and constitutive and stress-induced gene expression in spruce (Picea spp.)[J]. Phytochemistry, 68(14): 1975-1991. [37] Sahani P A, Ujinwal M, Singh S, et al.2018. Genome wide in silico characterization of dirigent protein family in flax (Linum usitatissmum L.)[J]. Plant Archives, 18: 61-68. [38] Shi H Y, Liu Z H, Li Z, et al.2012. Overexpression of cotton (Gossypium hirsutum) dirigent1 gene enhances lignification that blocks the spread of Verticillium dahliae[J]. Acta Biochimica et Biophysica Sinica, 44(007): 555-564. [39] Singh A, Jha S K, Bagri J, et al.2015. ABA inducible rice protein phosphatase 2C confers ABA insensitivity and abiotic stress tolerance in Arabidopsis[J]. PLOS ONE, 10(4): e0125168. [40] Song M, Peng X Y.2019. Genome-Wide identification and characterization of DIR genes in Medicago truncatula[J]. Biochemical Genetics, 57(4): 487-506. [41] Tamura K, Stecher G, Peterson D, et al.2013. MEGA6: Molecular evolutionary genetics analysis version 6.0[J]. Molecular Biology & Evolution, 30(12): 2725-2729. [42] Tiburzy R, Reisener H J.1990. Resistance of wheat to Puccinia graminis f.sp. tritici: Association of the hypersensitive reaction with the cellular accumulation of lignin-like material and callose[J]. Physiological & Molecular Plant Pathology, 36(2): 109-120. [43] Vanholme R, Morreel K, Darrah C, et al.2012. Metabolic engineering of novel lignin in biomass crops[J]. New Phytologist, 196(4): 978-1000. [44] Wang Y, Tang H, Debarry J D, et al.2012. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity[J]. Nucleic Acids Research, 40(7): e49. [45] Yadav C B, Bonthala V S, Muthamilarasan M, et al.2015. Genome-wide development of transposable elements-based markers in foxtail millet and construction of an integrated database[J]. DNA Research, 22(1): 79-90. [46] Zhou J, Lee C, Zhong R, et al.2009. MYB58 and MYB63 are transcriptional activators of the lignin biosynthetic pathway during secondary cell wall formation in Arabidopsis[J]. The Plant Cell, 21(1): 248-266.
