Large-scale comparative analysis of the nuclear factor-Y transcription factors across 320 horticultural and other plants

Abstract

Nuclear factor Y (NF-Y), evolutionarily conserved heterotrimeric transcription factors (TFs), are found throughout eukaryotic organisms. Comprising the NF-YA, NF-YB, and NF-YC subfamilies, this family is established as playing critical roles in plant growth and development. While earlier research has mainly centered on the functional and evolutionary characteristics of NF-Y within individual plant species, large-scale analyses and evolutionary patterns of these genes across major plant lineages remain largely unexplored. Here, we systematically identified 15 392 nonredundant genes of NF-Y family from 320 horticultural and representative plant species. Our findings showed that this gene family originated from charophytes. In bryophytes, pteridophytes, and gymnosperms, dispersed duplication served as the predominant mode of NF-Y gene expansion, whereas in angiosperms, their expansion was driven by whole genome duplication/segmental, dispersed, and tandem duplication. Conserved motif analysis revealed that highly conserved motifs are present within each NF-Y subfamily across eight representative plant species. However, some NF-Y genes in higher plants exhibited motif loss, indicating sequence variations during their evolutionary history. Transcriptomic profiling analysis of NF-Y genes in Arabidopsis thaliana under various conditions, including hormonal treatments, abiotic/biotic stresses, as well as various developmental stages, revealed their functional versatility. Furthermore, an interaction network comprising 36 NF-Y genes along with 2473 downstream and 261 upstream genes was constructed in A. thaliana. Enrichment analysis revealed interactions between NF-Y genes and other TFs, particularly those from the Myb_DNA-binding and APETALA2 (AP2) families, which were consistently enriched among both upstream and downstream regulatory genes. This work provides the first comprehensive and large-scale investigation into the evolutionary dynamics of NF-Y genes, encompassing taxa from basal algae to advanced horticultural plants, thereby offering novel insights into their evolutionary and lineage-specific expansion.

Acknowledgements

This research was supported by the National Natural Science Foundation of China (32060392), Qian Ke He Cheng Guo ([2022] Zhong Dian 005), and the GZMARS-Forage Industry Technology System of Guizhou Province.

Authors contributions

R.D. designed the project and was responsible for its execution. Data collection was performed by K.L. X.J., Mi.L., Ma.L., Z.L., and X.Z. Bioinformatics analysis was performed by K.L., Mi.L., Ma.L., X.J., J.C., X.G., J.H., and C.C. The manuscript was written and revised by R.D. and K.L. All authors read and revised the manuscript.

Data availability

All materials and related datasets in this study are available in Table S1.

Conflicts of interest statement

The authors declare no competing interests.

Supplementary material

Supplementary material is available at Horticulture Research online.

References

[1] Kaufmann K, Airoldi CA. Master regulatory transcription fac-tors in plant development:a blooming perspective. In: Yam-aguchi N,ed. Plant Transcription Factors:Methods and Proto-cols. Springer New York: New York, NY, 2018,3-22

[2] Franco-Zorrilla JM, López-Vidriero I, Carrasco JL. et al. DNA binding specificities of plant transcription factors and their potential to define target genes. Proc Natl Acad Sci USA. 2014; 111:2367-72

[3] de Mendoza A, Sebé-Pedrós A, šestak MS. et al. Transcription factor evolution in eukaryotes and the assembly of the reg-ulatory toolkit in multicellular lineages. P Natl Acad Sci USA. 2013;110:E4858-66

[4] Sinha S, Kim IS, Sohn KY. et al. Three classes of mutations in the A subunit of the CCAAT-binding factor CBF delineate functional domains involved in the three step assembly of the CBF-DNA complex. Mol Cell Biol. 1996; 16:328-37

[5] Bucher P, Trifonov EN. CCAAT box revisited: bidirectionality, location and context. J Biomol Struct Dyn. 1988; 5:1231-6

[6] Myers ZA, Holt BF. NUCLEAR FACTOR-Y: still complex after all these years? Curr Opin Plant Biol. 2018; 45:96-102

[7] Maity SN, Crombrugghe B. Role of the CCAAT binding protein CBF/NF-Y in transcription. Trends Biochem Sci. 1998; 23:174-8

[8] Romier C, Cocchiarella F, Mantovani R. et al. The NF-YB/NF-YC structure gives insight into DNA binding and transcription reg-ulation by CCAAT factor NF-Y. J Biol Chem. 2003; 278:1336-45

[9] Nardini M, Gnesutta N, Donati G. et al. Sequence-specific tran-scription factor NF-Y displays histone-like DNA binding and H2B like ubiquitination. Cell. 2013; 152:132-43

[10] Oldfield AJ, Yang P, Conway AE. et al. Histone-fold domain protein NF-Y promotes chromatin accessibility for cell type-specific master transcription factors. Mol Cell. 2014; 55:708-22

[11] Laloum T, de Mita S, Gamas P. et al. CCAAT-box binding transcription factors in plants: Y so many? Trends Plant Sci. 2013; 18:157-66

[12] Siefers N, Dang KK, Kumimoto RW. et al. Tissue-specific expression patterns of Arabidopsis NF-Y transcription fac-tors suggest potential for extensive combinatorial complexity. Plant Physiol. 2009; 149:625-41

[13] Yang W, Lu Z, Xiong Y. et al. Genome-wide identification and co-expression network analysis of the OsNF-Y gene family in rice. Crop J. 2017; 5:21-31

[14] Zhang F, Han M, Lv Q. et al. Identification and expression pro-file analysis of NUCLEAR FACTOR-Y families in Physcomitrella patens. Front Plant Sci. 2015; 6:642

[15] Li X, Li C, Shi L. et al. Jasmonate signaling pathway confers salt tolerance through a NUCLEAR FACTOR-Y trimeric transcription factor complex in Arabidopsis. Cell Rep. 2024; 43:113825

[16] Petroni K, Kumimoto RW, Gnesutta N. et al. The promiscuous life of plant nuclear factor Y transcription factors. Plant Cell. 2012; 24:4777-92

[17] Baud S, Kelemen Z, Thévenin J. et al. Deciphering the molec-ular mechanisms underpinning the transcriptional control of gene expression by master transcriptional regulators in Ara-bidopsis seed. Plant Physiol. 2016; 171:1099-112

[18] Braybrook SA, Harada JJ. LECs go crazy in embryo develop-ment. Trends Plant Sci. 2008; 13:624-30

[19] Jo L, Pelletier JM, Harada JJ. Central role of the LEAFY COTYLEDON1 transcription factor in seed development. J Integr Plant Biol. 2019; 61:564-80

[20] Kumimoto RW, Zhang Y, Siefers N. et al. NF-YC3, NF-YC4 and NF-YC9 are required for CONSTANS-mediated, photoperiod-dependent flowering in Arabidopsis thaliana. Plant J. 2010; 63:379-91

[21] Lotan T, Ohto M, Yee KM. et al. Arabidopsis LEAFY COTYLE-DON1 is sufficient to induce embryo development in vegeta-tive cells. Cell. 1998; 93:1195-205

[22] Swain S, Myers ZA, Siriwardana CL. et al. The multifaceted roles of NUCLEAR FACTOR-Y in Arabidopsis thaliana develop-ment and stress responses. Biochim Biophys Acta Gene Regul Mech. 2017; 1860:636-44

[23] West MA, Yee KM, Danao J. et al. LEAFY COTYLEDON1 is an essential regulator of late embryogenesis and COTYLEDON identity in Arabidopsis. Plant Cell. 1994; 6:1731-45

[24] Liu G, Zhang R, Li S. et al. TaMADS29 interacts with TaNF-YB1 to synergistically regulate early grain development in bread wheat. Sci China Life Sci. 2023; 66:1647-64

[25] Wei Q, Ma C, Xu Y. et al. Control of chrysanthemum flowering through integration with an aging pathway. Nat Commun. 2017; 8:829-39

[26] Bai A, Lu X, Li D. et al. NF-YB1-regulated expression of sucrose transporters in aleurone facilitates sugar loading to rice endosperm. Cell Res. 2016; 26:384-8

[27] Lee D, Kim HI, Jang G. et al. The NF-YA transcription factor OsNF-YA7 confers drought stress tolerance of rice in an abscisic acid independent manner. Plant Sci. 2015; 241:199-210

[28] Baudin M, Laloum T, Lepage A. et al. A phylogenetically con-served group of NF-Y transcription factors interact to control nodulation in legumes. Plant Physiol. 2015; 169:2761-73

[29] Combier JP, Frugier F, de Billy F. et al. MtHAP2-1 isakey transcriptional regulator of symbiotic nodule development regulated by microRNA169 in Medicago truncatula. Genes Dev. 2006; 20:3084-8

[30] Soyano T, Kouchi H, Hirota A. et al. Nodule inception directly targets NF-Y subunit genes to regulate essential processes of root nodule development in Lotus japonicus. PLoS Genet. 2013; 9:e1003352

[31] Zanetti ME, Blanco FA, Beker MP. et al. A C subunit of the plant nuclear factor NF-Y required for rhizobial infection and nodule development affects partner selection in the common bean-rhizobium etli symbiosis. Plant Cell. 2010; 22: 4142-57

[32] Alam MM, Tanaka T, Nakamura H. et al. Overexpression of a rice heme activator protein gene (OsHAP2E) confers resistance to pathogens, salinity and drought, and increases photosyn-thesis and tiller number. Plant Biotechnol J. 2015; 13:85-96

[33] O’Conner S, Zheng W, Qi M. et al. GmNF-YC4- 2 increases pro-tein, exhibits broad disease resistance and expedites maturity in soybean. Int J Mol Sci. 2021; 22:3586

[34] Rey T, Laporte P, Bonhomme M. et al. MtNF-YA1, a central transcriptional regulator of symbiotic nodule development, is also a determinant of Medicago truncatula susceptibility toward a root pathogen. Front Plant Sci. 2016; 7:1837

[35] Lv M, Cao H, Wang X. et al. Identification and expression anal-ysis of maize NF-YA subunit genes. PeerJ. 2022; 10:e14306

[36] Zhang Z, Li X, Zhang C. et al. Isolation, structural analysis, and expression characteristics of the maize nuclear factor Y gene families. Biochem Biophys Res Commun. 2016; 478: 752-8

[37] Quach TN, Nguyen HTM, Valliyodan B. et al. Genome-wide expression analysis of soybean NF-Y genes reveals poten-tial function in development and drought response. Mol Gen Genomics. 2015; 290:1095-115

[38] Xuanyuan G, Lian Q, Jia R. et al. Genome-wide screening and identification of nuclear factor-Y family genes and exploration their function on regulating abiotic and biotic stress in potato (Solanum tuberosum L.). Gene. 2022; 812:146089

[39] Ren C, Zhang Z, Wang Y. et al. Genome-wide identification and characterization of the NF-Y gene family in grape (Vitis vinifera L.). BMC Genomics. 2016; 17:605

[40] Zheng X, Zhang H, Zhang L. et al. Identification and compre-hensive analysis of the nuclear factor-Y family genes reveal their multiple roles in response to nutrient deficiencies in Brassica napus. Int J Mol Sci. 2021; 22:10354

[41] Li M, Du Q, Li J. et al. Genome-wide identification and chilling stress analysis of the NF-Y gene family in melon. Int J Mol Sci. 2023; 24:6934

[42] Li S, Li K, Ju Z. et al. Genome-wide analysis of tomato NF-Y fac-tors and their role in fruit ripening. BMC Genomics. 2016; 17:36

[43] Li M, Li G, Liu W. et al. Genome-wide analysis of the NF-Y gene family in peach (Prunus persica L.). BMC Genomics. 2019; 20:612

[44] An Y, Suo X, Niu Q. et al. Genome-wide identification and anal-ysis of the NF-Y transcription factor family reveal its potential roles in salt stress in alfalfa (Medicago sativa L.). Int J Mol Sci. 2022; 23:6426

[45] Qu Y, Wang Y, Zhu J. et al. Genomic organization, phy-logenetic comparison, and differential expression of the nuclear factor-Y gene family in apple (Malus domestica). Plants. 2020; 10:16

[46] Liu R, Wu M, Liu HL. et al. Genome-wide identification and expression analysis of the NF-Y transcription factor family in Populus. Physiol Plant. 2021; 171:309-27

[47] Wang T, Zou H, Ren S. et al. Genome-wide identification, char-acterization, and expression analysis of NF-Y gene family in Ginkgo biloba seedlings and GbNF-YA6 involved in heat-stress response and tolerance. Int J Mol Sci. 2023; 24:12284

[48] Freeling M. Bias in plant gene content following different sorts of duplication: tandem, whole-genome, segmental, or by transposition. Annu Rev Plant Biol. 2009; 60:433-53

[49] Innan H, Kondrashov F. The evolution of gene duplica-tions: classifying and distinguishing between models. Nat Rev Genet. 2010; 11:97-108

[50] Song X, Ma X, Li C. et al. Comprehensive analyses of the BES1 gene family in Brassica napus and examination of their evo-lutionary pattern in representative species. BMC Genomics. 2018; 19:346

[51] Yu T, Bai Y, Liu Z. et al. Large-scale analyses of heat shock transcription factors and database construction based on whole-genome genes in horticultural and representative plants. Hortic Res. 2022;9:uhac035

[52] Wang Y, Tang H, Debarry JD. et al. MCScanX: a toolkit for detec-tion and evolutionary analysis of gene synteny and collinear-ity. Nucleic Acids Res. 2012; 40:e49

[53] Gu C, Guo ZH, Hao PP. et al. Multiple regulatory roles of AP2/ERF transcription factor in angiosperm. Bot Stud. 2017; 58:6

[54] 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

[55] Feng S, Li N, Chen H. et al. Large-scale analysis of the ARF and Aux/IAA gene families in 406 horticultural and other plants. Mol Hortic. 2024; 4:13

[56] Kahle J, Baake M, Doenecke D. et al. Subunits of the het-erotrimeric transcription factor NF-Y are imported into the nucleus by distinct pathways involving importin beta and importin 13. Mol Cell Biol. 2005; 25:5339-54

[57] Hou X, Zhou J, Liu C. et al. Nuclear factor Y-mediated H3K27me 3 demethylation of the SOC1 locus orchestrates flowering responses of Arabidopsis. Nat Commun. 2014; 5:4601

[58] Tang Y, Liu X, Liu X. et al. Arabidopsis NF-YCs mediate the light-controlled hypocotyl elongation via modulating histone acetylation. Mol Plant. 2017; 10:260-73

[59] Cannon SB, Mitra A, Baumgarten A. et al. The roles of segmen-tal and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004; 4:10

[60] Kong H, Landherr LL, Frohlich MW. et al. Patterns of gene duplication in the plant SKP1 gene family in angiosperms: evidence for multiple mechanisms of rapid gene birth. Plant J. 2007; 50:873-85

[61] Panchy N, Lehti-Shiu M, Shiu SH. Evolution of gene duplica-tion in plants. Plant Physiol. 2016; 171:2294-316

[62] Soltis PS, Soltis DE. Ancient WGD events as drivers of key inno-vations in angiosperms. Curr Opin Plant Biol. 2016; 30:159-65

[63] Moharana KC, Venancio TM. Polyploidization events shaped the transcription factor repertoires in legumes (Fabaceae). Plant J. 2020; 103:726-41

[64] Song XM, Wang JP, Sun PC. et al. Preferential gene retention increases the robustness of cold regulation in Brassicaceae and other plants after polyploidization. Hortic Res. 2020; 7:20

[65] Mantovani R. The molecular biology of the CCAAT binding factor NF-Y. Gene. 1999; 239:15-27

[66] Zhao H, Wu D, Kong F. et al. The Arabidopsis thaliana nuclear factor Y transcription factors. Front Plant Sci. 2017; 07:2045

[67] Cao S, Kumimoto RW, Gnesutta N. et al. A distal CCAAT/nuclear factor Y complex promotes chromatin looping at the flowering locus T promoter and regulates the timing of flowering in Arabidopsis. Plant Cell. 2014; 26:1009-17

[68] Wenkel S, Turck F, Singer K. et al. CONSTANS and the CCAAT box binding complex share a functionally important domain and interact to regulate flowering of Arabidopsis. Plant Cell. 2006; 18:2971-84

[69] Lee H, Fischer RL, Goldberg RB. et al. Arabidopsis LEAFY COTYLEDON1 represents a functionally specialized subunit of the CCAAT binding transcription factor. Proc Natl Acad Sci USA. 2003; 100:2152-6

[70] Sato H, Suzuki T, Takahashi F. et al. NF-YB2 and NF-YB3 have functionally diverged and differentially induce drought and heat stress-specific genes. Plant Physiol. 2019; 180:1677-90

[71] De Clercq I, Van de Velde J, Luo X. et al. Integrative inference of transcriptional networks in Arabidopsis yields novel ROS signalling regulators. Nat Plants. 2021; 7:500-13

[72] Hwang K, Susila H, Nasim Z. et al. Arabidopsis ABF3 and ABF4 transcription factors act with the NF-YC complex to regulate SOC1 expression and mediate drought-accelerated flowering. Mol Plant. 2019; 12:489-505

[73] Zhou Y, Zhang Y, Wang X. et al. Root-specific NF-Y family transcription factor, PdNF-YB21, positively regulates root growth and drought resistance by abscisic acid-mediated indoylacetic acid transport in Populus. New Phytol. 2020; 227:407-26

[74] Xu JJ, Zhang XF, Xue HW. Rice aleurone layer specific OsNF-YB1 regulates grain filling and endosperm development by interacting with an ERF transcription factor. J Exp Bot. 2016; 67:6399-411

[75] Ribeiro C, de Melo BP, Lourenço-Tessutti IT. et al. The regen-eration conferring transcription factor complex ERF115-PAT1 coordinates a wound-induced response in root-knot nema-tode induced galls. New Phytol. 2024; 241:878-95

[76] Chen C, Wu Y, Li J. et al. TBtools-II: a "one for all, all for one" bioinformatics platform for biological big-data mining. Mol Plant. 2023; 16:1733-42

[77] Eddy SR. Accelerated profile HMM searches. PLoS Comput Biol. 2011; 7:e1002195

[78] Paysan-Lafosse T, Blum M, Chuguransky S. et al. InterPro in 2022. Nucleic Acids Res. 2023;51:D418-27

[79] Lu S, Wang J, Chitsaz F. et al. CDD/SPARCLE:the conserved domain database in 2020. Nucleic Acids Res. 2020;48:D265-8

[80] Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004; 32:1792-7

[81] Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009; 26:1641-50

[82] Letunic I, Bork P. Interactive tree of life (iTOL) v6: recent updates to the phylogenetic tree display and annotation tool. Nucleic Acids Res. 2024;52:W78-82

[83] Stolzer M, Lai H, Xu M. et al. Inferring duplications, losses, transfers and incomplete lineage sorting with nonbinary species trees. Bioinformatics. 2012;28:i409-15

[84] Bailey TL, Williams N, Misleh C. et al. MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006;34:W369-73

[85] Winter D, Vinegar B, Nahal H. et al. An “electronic fluorescent pictograph” browser for exploring and analyzing large-scale biological data sets. PLoS One. 2007; 2:e718

[86] Jacomy M, Venturini T, Heymann S. et al. ForceAtlas2, a continuous graph layout algorithm for handy network visualization designed for the Gephi software. PLoS One. 2014; 9:e98679

[87] Cantalapiedra CP, Hernández-Plaza A, Letunic I. et al. eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol Biol Evol. 2021; 38:5825-9

[88] Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. JRStat Soc Ser B. 1995; 57:289-300