Comprehensive analysis of NAC transcription factors and their expression during fruit spine development in cucumber (Cucumis sativus L.)

Abstract

The cucumber (Cucumis sativus L.) is an important vegetable crop worldwide, and fruit trichomes or spines are an important trait for external fruit quality. The mechanisms underlying spine formation are not well understood, but the plant-specific NAC family of transcription factors may play important roles in fruit spine initiation and development. In this study, we conducted a genome-wide survey and identified 91 NAC gene homologs in the cucumber genome. Clustering analysis classified these genes into six subfamilies; each contained a varying number of NAC family members with a similar intron–exon structure and conserved motifs. Quantitative real-time PCR analysis revealed tissue-specific expression patterns of these genes, including 10 and 12 that exhibited preferential expression in the stem and fruit, respectively. Thirteen of the 91 NAC genes showed higher expression in the wild-type plant than in its near-isogenic trichome mutant, suggesting their important roles in fruit spine development. Exogenous application of four plant hormones promoted spine formation and increased spine density on the cucumber fruits; several NAC genes showed differential expression over time in response to phytohormone treatments on cucumber fruit, implying their essential roles in fruit-trichome development. Among the NAC genes identified, 12 were found to be targets of 13 known cucumber micro-RNAs. Collectively, these findings provide a useful resource for further analysis of the interactions between NAC genes and genes underlying trichome organogenesis and development during fruit spine development in cucumber.

References

[1] Chapman, E. J. & Estelle, M. Mechanism of auxin-regulated gene expression in plants. Annu. Rev. Genet. 43, 265-285 (2009).

[2] Hu, R. et al. Comprehensive analysis of NAC domain transcription factor gene family in Populus trichocarpa. BMC Plant Biol. 10, 145 (2010).

[3] Shah, S. et al. Molecular cloning and functional analysis of NAC family genes associated with leaf senescence and stresses in Gossypium hirsutum L.. Plant Cell. Tiss. Organ Cult. 117, 167-186 (2014).

[4] Guo, Y. & Gan, S. AtNAP, a NAC family transcription factor, has an important role in leaf senescence. Plant J. 46, 601-612 (2006).

[5] Wang, L. et al. Genome-wide survey of pseudogenes in 80 fully re-sequenced Arabidopsis thaliana accessions. PLoS One. 7(12), e51769 (2012).

[6] Yamasaki, K. et al. Structures and evolutionary origins of plant-specific transcription factor DNA-binding domains. Plant Physiol. Biochem. 46, 394-401 (2008).

[7] Wang, Z. et al. Genome-wide analysis of the basic leucine zipper (bZIP) transcription factor gene family in six legume genomes. BMC Genom. 16, 1053 (2015).

[8] Kavas, M. et al. Genome-wide characterization and expression analysis of common bean bHLH transcription factors in response to excess salt concentration. Mol. Genet. Genom. 29, 129-143 (2016).

[9] Souer, E., van Houwelingen, A., Kloos, D., Mol, J. & Koes, R. The no apical meristem gene of Petunia is required for pattern formation in embryos and flowers and is expressed at meristem and primordia boundaries. Cell 85, 159-170 (1996).

[10] Aida, M., Ishida, T., Fukaki, H., Fujisawa, H. & Tasaka, M . Genes involved in organ separation in Arabidopsis: an analysis of the cup-shaped cotyledon mutant. Plant. Cell. 9, 841-857 (1997).

[11] Gong, W. et al. Genome-wide ORFeome cloning and analysis of Arabidopsis transcription factor genes. Plant Physiol. 135, 773-782 (2004).

[12] Xiong, Y. et al. Transcription factors in rice: a genome-wide comparative analysis between monocots and eudicots. Plant Mol. Biol. 59, 191-203 (2005).

[13] Pinheiro, G. et al. Complete inventory of soybean NAC transcription factors: sequence conservation and expression analysis uncover their distinct roles in stress response. Gene 444, 10-23 (2009).

[14] Huang, S. et al. Genome-wide analysis of WRKY transcription factors in Solanum lycopersicum. Mol. Genet. Genom. 287, 495-513 (2012).

[15] Hu, W. et al. Genome-wide identification and expression analysis of the NAC transcription factor family in cassava. PLoS One. 10, e0136993 (2015).

[16] Cenci, A., Guignon, V., Roux, N. & Rouard, M. Genomic analysis of NAC transcription factors in banana (Musa acuminata) and definition of NAC orthologous groups for monocots and dictos. Plant Mol. Biol. 85, 63-80 (2014).

[17] Wang, N., Zheng, Y., Xin, H., Fang, L. & Li, S. Comprehensive analysis of NAC domain transcription factor gene family in Vitis vinifera. Plant Cell. Rep. 32, 61-75 (2013).

[18] Ernst, H. A., Olsen, A. N., Skriver, K., Larsen, S. & Leggio, L. L. Structure of the conserved domain of ANAC, a member of the NAC family of transcription factors. EMBO Rep. 5, 297-303 (2004).

[19] Olsen, A. N., Ernst, H. A., Leggio, L. L. & Skriver, K. NAC transcription factors: structurally distinct, functionally diverse. Trends Plant. Sci. 10, 79-87 (2005).

[20] Jensen, M. et al. The Arabidopsis thaliana NAC transcription factor family: structure-function relationships and determinants of ANAC019 stress signaling. Biochem. J. 426, 183-196 (2010).

[21] Ooka, H. et al. Comprehensive analysis of NAC family genes in Oryza sativa and Arabidopsis thaliana. DNA Res. 10, 239-247 (2003).

[22] Chen, Q., Wang, Q., Xiong, L. & Lou, Z. A structural view of the conserved domain of rice stress-responsive NAC1. Protein Cell 2, 55-63 (2011).

[23] Lu, M. et al. A maize stress-responsive NAC transcription factor, ZmSNAC1, confers enhanced tolerance to dehydration in transgenic Arabidopsis. Plant Cell Rep. 31, 1701-1711 (2012).

[24] Puranik, S., Sahu, P. P., Srivastava, P. S. & Prasad, M. NAC proteins: regulation and role in stress tolerance. Trends Plant. Sci. 17, 369-381 (2012).

[25] Fang, Y., You, J., Xie, K., Xie, W. & Xiong, L. Systematic sequence analysis and identification of tissue-specific or stress-responsive genes of NAC transcription factor family in rice. Mol. Genet. Genom. 280, 547-563 (2008).

[26] Xie, Q. et al. SINAT5 promotes ubiquitin-related degradation of NAC1 to attenuate auxin signals. Nature 419, 167-170 (2002).

[27] Hao, Y. et al. Plant NAC-type transcription factor proteins contain a NARD domain for repression of Transcriptional activation. Planta 232, 1033-1043 (2010).

[28] Seo, P. J., Kim, S. G. & Park, C. M. Membrane-bound transcription factors in plants. Trends Plant. Sci. 13, 550-556 (2008).

[29] Shen, H., Yin, Y., Chen, F., Xu, Y. & Dixon, R. A. A bioinformatic analysis of NAC genes for plant cell wall development in relation to lignocellulosic bioenergy production. Bioenerg. Res. 2, 217-232 (2009).

[30] Kim, S. et al. Genome-scale screening and molecular characterization of membrane-bound transcription factors in Arabidopsis and rice. Genomics 95, 56-65 (2010).

[31] Le, D. et al. Genome-wide survey and expression analysis of the plant-specific NAC transcription factor family in soybean during development and dehydration stress. DNA Res. 18, 263-276 (2011).

[32] Nakashima, K., Takasaki, H., Mizoi, J., Shinozaki, K. & Yamaguchi-Shinozaki, K. NAC transcription factors in plant abiotic stress responses. Biochem. Biophys. Acta 1819, 97-103 (2012).

[33] Mao, H. et al. A transposable element in a NAC gene is associated with drought tolerance in maize seedlings. Nat. Commun. 6, 8326 (2015).

[34] Podzimska-Sroka, D., O’Shea, C., Gregersen, P. L. & Skriver, K. NAC transcription factors in senescence: from molecular structure to function in crops. Plants 4, 412-448 (2015).

[35] Shiriga, K. et al. Genome-wide identification and expression pattern of drought-responsive members of the NAC family in maize. Meta Gene 2, 407-417 (2014).

[36] Kikuchi, K. et al. Molecular analysis of the NAC gene family in rice. Mol. Gen. Genet. 262, 1047-1051 (2000).

[37] Xie, Q., Frugis, G., Colgan, D. & Chua, N. H. Arabidopsis NAC1 transduces auxin signal downstream of TIR1 to promote lateral root development. Genes Dev. 14, 3024-3036 (2000).

[38] Mitsuda, N., Seki, M., Shinozaki, K. & Ohme-Takagi, M. The NAC transcription factors NST1 and NST2 of Arabidopsis regulate secondary wall thickenings and are required for anther dehiscence. Plant Cell 17, 2993-3006 (2005).

[39] Mitsuda, N. et al. NAC transcription factors, NST1 and NST3, are key regulators of the formation of secondary walls in woody tissues of Arabidopsis. Plant Cell 19, 270-280 (2007).

[40] Larsson, E., Sitbon, F., Sundström, J. & von Arnold, S. Expression of PaNAC01, a Picea abies CUP-SHAPED COTYLEDON orthologue, is regulated by polar auxin transport and associated with differentiation of the shoot apical meristem and formation of separated cotyledons. Ann. Bot. 110, 923-934 (2012).

[41] Lan, Y. et al. Identification of novel MiRNAs and MiRNA expression profiling during grain development in indica rice. BMC Genom. 13, 264 (2012).

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

[43] Zhao, J. et al. Transcriptome profiling of trichome-less reveals genes associated with multicellular trichome development in Cucumis sativus. Mol. Genet. Genom. 290, 2007-2018 (2015).

[44] Liu, X. W., Bartholonew, E., Cai, Y. L. & Ren, H. Z. Trichome-related mutants provide a new perspective on multicellular trichome initiation and development in cucumber (Cucumis sativus. L). Front. Plant Sci. 7, 1187 (2016).

[45] Li, Y. H., Wen, C. L. & Weng, Y. Q. Fine mapping of the pleiotropic locus B for black spine and orange mature fruit color in cucumber identifies a 50kb region containing a R2R3-MYB transcription factor. Theor. Appl. Genet. 126, 2187-2196 (2013).

[46] Li, Q. 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. 66, 2515-2516 (2015).

[47] Li, L., Stoeckert, C. J. & Roos, D. S. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178-2189 (2003).

[48] Krzywinski, M. et al. Circos: an information aesthetic for comparative genomics. Genome Res. 19, 1639-1645 (2009).

[49] Wang, W. et al. Cucumber CECRIFRUMI1 (CsCER1), which influences the cuticle properties and drought tolerance of cucumber, plays a key role in VLC alkanes biosynthesis. Plant Mol. Biol. 87, 219-233 (2015).

[50] Huang, S. et al. The genome of the cucumber, Cucumis sativus L.. Nat. Genet. 41, 1275-1281 (2009).

[51] An, L., Zhou, Z., Yan, A. & Gan, Y. Progress on trichome development regulated by phytohormones signaling. Plant Signal. Behav. 6, 1959-1962 (2011).

[52] Liu, J., Chen, K., Zhang, Z., Chen, X. & Wang, A. Effects of exogenous GA, Me-JA, IAA, SA and KT on trichome formation in tomato. Acta Horti Sin. (Chin.) 43, 2151-2160 (2016).

[53] Nuruzzaman, M. et al. Genome-wide analysis of NAC transcription factor family in rice. Gene 465, 30-44 (2010).

[54] Zhang, X. et al. Response gene SlARF3 plays multiple roles in tomato development and is involved in the formation of epidermal cells and trichomes. Plant Cell Physiol. 56, 2110-2124 (2015).

[55] Greenboim-Wainberg, Y. et al. Cross talk between gibberellin and cytokinin: the Arabidopsis GA response inhibitor SPINDLY plays a positive role in cytokinin signaling. Plant Cell 17, 92-102 (2005).