Genome-wide identification and expression profile analysis of nuclear factor Y family genes in Sorghum bicolor L. (Moench)

Authors: P. Maheshwari aff001;  Divya Kummari aff002;  Sudhakar Reddy Palakolanu aff002;  U. Nagasai Tejaswi aff003;  M. Nagaraju aff001;  G. Rajasheker aff001;  G. Jawahar aff001;  N. Jalaja aff003;  P. Rathnagiri aff005;  P. B. Kavi Kishor aff001
Authors place of work: Department of Genetics, Osmania University, Hyderabad, India aff001;  International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Patancheru, Hyderabad, India aff002;  Department of Biotechnology, Vignan’s Foundation for Science, Technology and Research, Vadlamudi, Guntur, Andhra Pradesh, India aff003;  Department of Biochemistry, ICMR-National Institute of Nutrition, Hyderabad, India aff004;  Genomix CARL Pvt. Ltd. Rayalapuram Road, Pulivendula, Kadapa, Andhra Pradesh, India aff005;  Genomix Molecular Diagnostics Pvt Ltd., Kukatpally, Hyderabad, India aff006;  Genomix Biotech Inc., Atlanta, GA, United States of America aff007
Published in the journal: PLoS ONE 14(9)
Category: Research Article
doi: 10.1371/journal.pone.0222203


Members of the plant Heme Activator Protein (HAP) or NUCLEAR FACTOR Y (NF-Y) are trimeric transcription factor complexes composed of the NF-YA, NF-YB and NF-YC subfamilies. They bind to the CCAAT box in the promoter regions of the target genes and regulate gene expressions. Plant NF-Ys were reported to be involved in adaptation to several abiotic stresses as well as in development. In silico analysis of Sorghum bicolor genome resulted in the identification of a total of 42 NF-Y genes, among which 8 code for the SbNF-YA, 19 for SbNF-YB and 15 for the SbNF-YC subunits. Analysis was also performed to characterize gene structures, chromosomal distribution, duplication status, protein subcellular localizations, conserved motifs, ancestral protein sequences, miRNAs and phylogenetic tree construction. Phylogenetic relationships and ortholog predictions displayed that sorghum has additional NF-YB genes with unknown functions in comparison with Arabidopsis. Analysis of promoters revealed that they harbour many stress-related cis-elements like ABRE and HSE, but surprisingly, DRE and MYB elements were not detected in any of the subfamilies. SbNF-YA1, 2, and 6 were found upregulated under 200 mM salt and 200 mM mannitol stresses. While NF-YA7 appeared associated with high temperature (40°C) stress, NF-YA8 was triggered by both cold (4°C) and high temperature stresses. Among NF-YB genes, 7, 12, 15, and 16 were induced under multiple stress conditions such as salt, mannitol, ABA, cold and high temperatures. Likewise, NF-YC 6, 11, 12, 14, and 15 were enhanced significantly in a tissue specific manner under multiple abiotic stress conditions. Majority of the mannitol (drought)-inducible genes were also induced by salt, high temperature stresses and ABA. Few of the high temperature stress-induced genes are also induced by cold stress (NF-YA2, 4, 6, 8, NF-YB2, 7, 10, 11, 12, 14, 16, 17, NF-YC4, 6, 12, and 13) thus suggesting a cross talk among them. This work paves the way for investigating the roles of diverse sorghum NF-Y proteins during abiotic stress responses and provides an insight into the evolution of diverse NF-Y members.


Biology and life sciences – Plant science – Plant physiology – Plant defenses – Plant resistance to abiotic stress – Plant pathology – Plant ecology – Plant-environment interactions – Ecology – Organisms – Eukaryota – Plants – Grasses – Sorghum – Brassica – Arabidopsis thaliana – Computational biology – Introns – Genetics – Genomics – Genome complexity – Gene expression – Gene regulation – Biochemistry – Proteins – DNA-binding proteins – Transcription factors – Regulatory proteins – Nucleic acids – RNA – Non-coding RNA – Natural antisense transcripts – MicroRNAs – Ecology and environmental sciences – Research and analysis methods – Database and informatics methods – Bioinformatics – Sequence analysis – Sequence motif analysis – Animal studies – Experimental organism systems – Model organisms – Plant and algal models – Physical sciences – Physics – Classical mechanics – Mechanical stress – Thermal stresses


Nuclear Factor Y (NF-Y), also known as heme activator protein (HAP) or CCAAT-binding factor (CBF) is a ubiquitous, complex, heterotrimeric transcription factor. It is evolutionarily conserved in all plants with three distinct subunits called NF-YA or HAP2, NF-YB or HAP3/CBF-A and NF-YC or HAP5/CBF-C [1]. The assembly of NF-Y is complex and occurs both in cytoplasm and nucleus. While NF-YA and NF-YC family members have a nuclear localization signal (NLS), NF-YB members generally lack the same and hence cannot be transported to the nucleus [2]. The NF-YA subunits are localised to the nucleus and bind with varying affinities to the CCAAT cis-elements in the promoter regions of the target genes [3, 4]. On the other hand, NF-YB and NF-YC subunits contain the conserved Histone Fold Domain (HFD) or Histone Fold Motif (HFM) and help in protein-DNA and protein-protein interactions [5, 4]. The HFD domain of NF-YB/YC is formed by 3 α-helices separated by 2 loops [6]. The α1-helices of both NF-YB and YC contain putative DNA-binding domains [4]. Via the HFDs, NF-YB and NF-YC form a heterodimer [5], which is critical for the translocation of NF-YB from the cytoplasm to the nucleus [7]. Several members of the NF-Y subfamilies play a vital role not only in a wide array of developmental processes but in tolerance to abiotic stresses as well. For example, they are involved in embryogenesis [8], ABA response and seed germination [9], abiotic stress tolerance [10, 11, 12], flowering time [13], primary root elongation [14], photosynthesis [15], endosperm development [16], and photomorphogenesis [17]. On the other hand, in leguminous plants, they are the key regulators of symbiotic root nodule development [18]. Ni et al. [19] reported that GmNF-YA3, a target gene of miR169, is a positive regulator of plant tolerance to drought stress in A. thaliana. Alam et al. [20] demonstrated that overexpression of OsHAP2E gene confers tolerance to drought and salt stresses with enhanced photosynthesis and tiller (stems produced in grass plants) number in comparison with wild-type rice plants. Transgenic rice showed resistance to Magnaporthe oryzae and Xanthomonas oryzae infections. While overexpression of NF-YA5 conferred drought stress tolerance in Arabidopsis [21], NF-YA1 in Arabidopsis resulted in post germinative growth arrest under salt stress [11]. Further, incorporation of several NF-YB and NF-YC genes improved drought stress tolerance in diverse plants like Arabidopsis, maize, poplar and rice [10, 22, 23, 24]. Thus, evidence is accumulating that NF-Y subunits act as key regulators of drought stress tolerance. In plants, NF-Y gene families comprise several paralogs. In Triticum aestivum, 37 paralogs have been described (10 NF-YA, 11 NF-YB, 14 NF-YC and 2 Dr1) [25]; in rice 28 NF-Ys (10 NF-YA, 11 NF-YB and 7 NF-YC) [26]; but later [27] reported 34 members in the same species. In Arabidopsis thaliana, NF-Y members exclude AtNF-YC11, B12, B13 (NC2 subfamily) and AtNF-YC10, C13, and B11 (Dpb3/4 subfamily), but include AtNF-YC12. Accordingly, Petroni et al. [3] pointed out that Arabidopsis contains 30 members of the NF-Y family, 10 from each family (NF-YA, NF-YB and NF-YC) though 36 were originally reported Siefers et al [13]. So, in the updated scheme, 30 members of NF-Y have been considered by Zhao et al. [28] in A. thaliana. In Brachypodium distachyon 36 (7 NF-YA, 17 NF-YB, and 12 NF-YC) [29]; in Brassica napus 33 (14 NF-YA, 14 NF-YB, 5 NF-YC) [30]; in Setaria italica 39 (10 NF-YA, 15 NF-YB and 14 NF-YC) [31]; in Glycine max 68 (21 NF-YA, 32 NF-YB, 15 NF-YC) [32]; in Prunus mume 29 [33]; in Ricinus communis 25 (6 NF-YA, 12 NF-YB and 7 NF-YC) [34]; in Citrus sinensis and C. clemantia 22 (6 NF-YA, 11 NF-YB, 5 NF-YC) [35] were characterised.

Sorghum bicolor is a semi-arid and the second most important staple food grain crop. It provides feed, fodder and fuel and shows genetic diversity [36]. Being a C4 photosynthetic plant, it is adapted to moderate drought and high temperature. But, salinity and drought coupled with high temperature limit the production and yield stability in sorghum. Enhancing the final yields and productivity of crop plants is especially challenging to the researchers due to the unpredictable nature of drought stress conditions during the growing season and complex drought stress biology [37]. Identification and expression of various transcription factors for abiotic stress tolerance using qRT-PCR and their validation by overexpression or knockouts is therefore critical for developing improved crop varieties with tolerance to water limited conditions. Members of NF-Y subfamilies impart tolerance to a very wide spectrum of both biotic and abiotic stresses as mentioned above. The number of NF-Y genes that exist in sorghum and their detailed biological roles for multiple stress tolerance and ABA-responsiveness remains unexplored. In this study, we identified 42 NF-Y genes using in silico approaches and examined their expression patterns under salt, drought, ABA, cold and high temperature stresses. Our gene expression studies reveal that majority of NF-Y genes (39) exhibited response to high temperature stress. A large number of them (24) were also expressed under multiple stresses like cold, salt (22) and drought (20). Further, 20 SbNF-Ys showed upregulation under ABA stress which indicates their role in ABA-related pathway. Keeping in view of the aforementioned reasons, we aimed to understand how the SbNF-Y members regulate abiotic stresses in an ABA-dependent or independent manner which would further delve into investigating their detailed roles during stress.

Materials and methods

Plant material and abiotic stress treatments

Sorghum bicolor variety BTx623 is an agronomically important inbred line. It is a model variety with known genome sequencing information and moderately tolerant to drought stress. The gene space of the sorghum genome sequence has also been updated by resequencing [38]. Keeping these criteria in mind, seeds of S. bicolor variety BTx623 were obtained from ICRISAT, Patancheru, Hyderabad, and sown in pots filled with 5 kg of black soil and seedlings were grown in glass house conditions at 28/20°C day/night temperatures. Sixty-day-old seedlings were subjected to 200 mM NaCl solution, 200 mM mannitol solution, and 100 μM ABA for 4 h separately. Cold stress was imposed by keeping the plants at 4°C and high temperature stress by exposing the plants to 40°C for 4 h. Corresponding controls (without any treatment) were maintained under identical conditions. After 4 h of exposure, roots, stems and leaves were collected from treated and control plants and snap frozen in liquid nitrogen and stored at -80°C for subsequent use. Three biological and three technical replicates were used for qRT-PCR analysis.

Identification and characterization of NF-Y transcription factors

NFY gene sequences of Oryza, Zea, Setaria were retrieved from plantTFDB ( (S1 Table) database and searched against Sorghum bicolor genome in Gramene database ( to find out their homologs. Genscan ( program was used to retrieve the gene and their respective protein sequences. The identified putative Sorghum nuclear factors were scanned using HMMER ( corresponding to the Pfam database and queried against the Oryza, Setaria and Zea. The identified Sorghum NF-Ys were confirmed by searching against Oryza, Setaria and Zea genomes in Gramene database ( Based on homology, the identified putative protein sequences were subjected to Motif Search ( analysis to check the reliability and to identify their conserved domains [39]. The identified NF-Y genes were mapped to their respective chromosomes based on the information provided in the Gramene Genome Database by employing MapInspect software ( Gene Structure Display Server ( software was used for obtaining the NF-Y gene structures—exons, introns, and untranslated sequence regions (UTRs) based on the alignments of their coding sequences [40]. MEME software was employed to analyze new sequence patterns and their significance [41]. The software helps to identify the nature of motifs by setting different default parameters, number of motifs from 1–10 with a motif width of 5–50, and the number of motif sites from 5–10.

NF-Y protein analysis, prediction of potential cis-regulatory elements, identification of miRNA target sites and phylogenetic analysis of NF-Ys

Molecular weight (MW), isoelectric point (pI), and GRAVY (grand average of hydropathy) of NF-Ys were identified for all NF-Y proteins by using ProtParam of Expasy tools ( [42]. Phosphorylation sites of proteins were predicted using NetPhos3.1 software of Expasy tools [43]. Subcellular localization of NF-Y proteins was carried out by WOLFPSORT program ( [44]. To predict the putative cis-acting elements of NF-Y promoter regions, 2000 bp genomic sequences upstream of start codons were analysed using PLANTCARE software [45]. The pSRNATarget software [46] was employed to identify the potential miRNA target sites in identified SbNF-Ys. Finally, the neighbor-joining (NJ) phylogenetic tree was constructed with the NF-Y protein sequences of Sorghum bicolor with the plants as shown in S1 Table using MEGA 6.2 software [47]. The NJ is a recursive algorithm, a fast method which is suited for large datasets and does not require ultra metric data and permits correction for multiple substitutions. The Poisson correction, pairwise deletion, and bootstrap value (1,000 replicates) parameters were used to draw the NJ phylogenetic tree.

Phylogenetic divergence and co-expression analysis

Gene duplication events were found [48, 49] using phylogenetic tree based on 70% similarity and 80% coverage of sequences aligned. PAL2NAL program [50] was followed for finding out synonymous and non-synonymous substitutions rates. Protein-protein interaction (PPI) map of NF-Y proteins was generated from the STRING database [51].

RNA isolation and quantitative real-time PCR analysis

From the stress exposed and control (without any stress) samples, total RNA was isolated using Macherey-Nagel NucleoSpin RNA plant kit by following the instructions given in the manual. To eliminate any genomic DNA contamination in the RNA samples, the purity of RNA was checked using Eppendorf BioPhotometer. Two micrograms of RNA sample was used as template for first strand cDNA synthesis using RevertAid First Strand cDNA Synthesis Kit (#K1622, Thermo Scientific EU, Reinach, Switzerland). To find out the relative gene expression levels of SbNF-Ys, 2X Applied Biosystems (ABI) Master Mix with gene specific primers was used (S2 Table). For qRT-PCR analysis, thermal cycling conditions of 95°C for 5 min followed by 40 cycles of 95°C for 30 s, 57°C for 30 s and 72°C f or 30 s were applied to the ABI 7500 real-time PCR system (Applied Biosystems, Foster City, CA, USA). Expression of SbNF-Y genes in control and treated samples was normalized with EIF4a (Eukaryotic Initiation Factor 4A) and PP2A (protein phosphatase2A subunit A3) reference genes [52]. qRT-PCR was carried out with three biological and three technical replicates for each sample. The PCR reaction specificity was confirmed by melting curve analysis of the amplicons. Comparative 2-DDCT method [53] was used to calculate the relative quantities of each transcript.


Identification and characterization of SbNF-Y transcription factors, motif analysis and subcellular localization

A total of 42 homologous genes comprising 8 NF-YA, 19 NF-YB and 15 NF-YC from the whole genome of sorghum were identified and confirmed. Later, they were crosschecked by using the HMM profile and searching SbNF-Ys against Oryza, Setaria and Zea for further confirm to check their reliability (S3 Table). The predicted 8 NF-YA, 19 NF-YB and 15 NF-YC genes were named as SbNF-YA1 to SbNF-YA8, SbNF-YB1 to SbNF-YB19 and SbNF-YC1 to SbNF-YC15 respectively. Based on the presence of conserved NF-YA, NF-YB and NF-YC domains, the predicted SbNF-Y family of proteins was considered for identification as a member. The exon-intron structures of all the 42 annotated NF-Y genes were analysed (Fig 1). While 17 exons and 16 introns (highest) were detected in SbNF-YA2, 3 exons (least) and 2 introns were found in SbNF-YA4 gene. Among the SbNF-YB family members, SbNF-YB11 showed a maximum of 16 exon and 15 intron regions, and 5 of the members displayed 1 intron. A maximum of 18 exons and 17 introns were noticed on SbNF-YC8. Also, six of the members were intronless and no member exhibited one intron (Fig 1). The sub-cellular localization of SbNF-Y proteins based on consensus sequence showed a majority of them to be localized to nucleolus and chloroplast although a few of them localized to cytoplasm, mitochondria and plastids (Table 1). All the NF-Ys showed nuclear localization signals (NLS); NF-YA holding LRRR sequence (motif 2 in Fig 2A), KRK motif in NF-YB (motif 1 in Fig 2B) and KRR in NF-YC (motif 1 in Fig 2C). Though they contain nuclear localization signals, their subcellular localizations were different. Majority of the SbNF-YAs showed chloroplast as the important target site. Few of them have been found localized in chloroplasts (NF-YA1, NF-YA7 and NF-YA8), mitochondria (NF-YA4) and plastid (NF-YA6). The number of phosphorylation sites in each NF-Y protein is represented in the S3 Table. All the NFYs exhibited higher number of PKC than CK1, CK2, and PKA types. The PKC number is higher in NF-YA subfamily members than in NF-YB (S4 Table). No transmembrane helices were observed except in NF-YA3 protein.

Fig. 1.

Gene structure analysis of SbNFYs (A. SbNF-YA; B. SbNF-YB and C. SbNF-YC).

Fig. 2.

Distribution of 1–10 conserved motifs in A) SbNF-YA, B) SbNF-YB and C) SbNF-YC groups. Gene clusters and p values are shown on the left side and motif sizes at the bottom of the figure.

Tab. 1.

List of identified SbNF-Ys exhibiting chromosomal location, sub group, length, DNA binding domains (DBD), molecular weight (MW), iso-electric point (pI), GRAVY, number of exons, protein localization, and instability index.

<h2>List of identified <i>SbNF-Ys</i> exhibiting chromosomal location, sub group, length, DNA binding domains (DBD), molecular weight (MW), iso-electric point (pI), GRAVY, number of exons, protein localization, and instability index.</h2>

The identified SbNF-Y genes encoded polypeptides with amino acids ranging from 130 to 1430 and pI values varied from 4.26 to 10.83. Characteristically, they showed DNA binding domains. Molecular weights of the proteins ranged from 10.21 to 83.52 kDa (Table 1). Among the SbNF-YA subfamily members, RKPYLHESRHLHAMKRARGSGGRFLNTKQ and EEPIYVNAKQYNAILRRRQARAKLEAZNK large contiguous motifs were found ubiquitous, while rest of the 8 large contiguous motifs showed variability in their distribution (see Figs 2A and S1). Similarly, SbNF-YB proteins revealed the uniform presence of one, highly conserved large contiguous motif, i.e. AKETVQECVSEFISFVTGEASDKCQREKRKTINGDDLLWALATLGFEDYY (Figs 2B and S2). On the other hand, analysis of NF-YC proteins revealed a highly conserved large contiguous motif APVVFAKACEFIQELTLRAWHEENKRRTLQKSSDIAAAIARTEVYDFL (see Figs 2C and S3). Motif analysis of complete SbNF-Y family representing conserved motifs (S4 and S5 Figs) reflect typical diagnostic features for different subunits of NF-Y family proteins in general and hence provide confirmatory identification of SbNF-Y proteins from the sorghum genome.

Phylogeny, divergence, and physical genome mapping of SbNF-Ys

The phylogenetic tree of SbNF-Y proteins was constructed using MEGA 6.2 software. It showed 2 clades, which are subdivided into 4. The phylogenetic analysis displayed a total of 11 paralogous duplication events of which 3 tandem/segmental (SbNF-YB4/SbNF-YB6; SbNF-YB18/SbNF-YB19; SbNF-YB12/SbNF-YC11) and remaining regional duplication events. Interestingly, NF-YBs exhibited paralogous events with SbNF-YCs (SbNF-YB13/SbNF-YC1; SbNF-YB12/SbNF-YC11), indicating their evolutionary relatedness (Figs 3 and S6 and Table 2). To find out the orthologous and the evolutionary relationships of SbNF-Ys, NF-YA, NF-YB and NF-YCs have been compared with other plant genomes (see Fig 4A, 4B and 4C respectively). Not surprisingly, a majority of the identified SbNF-Ys showed orthologous relationship with Zea, few with Setaria and one with Hordeum (Sorbi009G166200 (SbNF-YB15)/MLOC_36879.2). Of the 8 SbNF-YAs, 5 showed orthologous events with Zea and 1 with Setaria (Fig 4 and S5 Table). The major subfamily SbNF-YBs exhibited 11 orthologous events, of which 6 showed with Zea, 4 with Setaria and 1 with Hordeum (Fig 5 and S6 Table). On the other hand, SbNF-YCs showed 10 orthologous events, of which 8 with Zea and 2 with Setaria (Fig 6 and S7 Table).

<h2>Phylogenetic tree of <i>SbNF-Y</i> gene family.</h2>
Fig. 3.

Phylogenetic tree of SbNF-Y gene family.

<h2>Phylogenetic tree showing the relationship between <i>Sorghum</i>, <i>Oryza</i>, <i>Setaria</i>, <i>Zea</i>, <i>Medicago</i>, <i>Glycine max</i>, <i>Daucus</i>, <i>Solanum</i>, <i>Brachypodium</i>, <i>Cajanus</i>, <i>Populus</i> and <i>Arabidopsis NF-YAs</i>.</h2>
Fig. 4.

Phylogenetic tree showing the relationship between Sorghum, Oryza, Setaria, Zea, Medicago, Glycine max, Daucus, Solanum, Brachypodium, Cajanus, Populus and Arabidopsis NF-YAs.

<h2>N-j phylogenetic tree representing <i>NF-YB</i> family genes of <i>Sorghum</i>, <i>Oryza</i>, <i>Setaria</i>, <i>Zea</i>, <i>Medicago</i>, <i>Glycine max</i>, <i>Dacus</i>, <i>Solanum</i>, <i>Brachypodium</i>, <i>Cajanus</i>, <i>Populus</i> and <i>Arabidopsis</i>.</h2>
Fig. 5.

N-j phylogenetic tree representing NF-YB family genes of Sorghum, Oryza, Setaria, Zea, Medicago, Glycine max, Dacus, Solanum, Brachypodium, Cajanus, Populus and Arabidopsis.

<h2>Evolutionary relationship of the NF-YC family genes of <i>Sorghum</i>, <i>Oryza</i>, <i>Setaria</i>, <i>Zea</i>, <i>Medicago</i>, <i>Glycine max</i>, <i>Dacus</i>, <i>Solanum</i>, <i>Brachypodium</i>, <i>Cajanus</i>, <i>Populus</i> and <i>Arabidopsis</i>.</h2>
Fig. 6.

Evolutionary relationship of the NF-YC family genes of Sorghum, Oryza, Setaria, Zea, Medicago, Glycine max, Dacus, Solanum, Brachypodium, Cajanus, Populus and Arabidopsis.

Tab. 2.

Calculation of non-synonymous to synonymous substitution ratios in paralogs of SbNFYs.

<h2>Calculation of non-synonymous to synonymous substitution ratios in paralogs of <i>SbNFYs</i>.</h2>

The identified SbNF-Ys were distributed across all the 10 chromosomes. A maximum number of 7 genes each were located on chromosome 1, and 7, 5 genes each on 2, 3, and 9, 4 each on 8 and 10, 3 genes on 4, and 1 each on 5 and 6 chromosomes (Fig 7 and Table 1). Among the NF-YA subfamily, SbNF-YA1, SbNF-YA2 and SbNF-YA3 genes were located on chromosome 1. Chromosomes 2 and 8 have two genes each located on them, i.e. SbNF-YA4, SbNF-YA5 and SbNF-YA7, SbNF-YA8 respectively. Chromosome 4 is having only SbNF-YA6 localized on it. Among the SbNF-YB genes, a maximum of 4 genes each were located on chromosomes 3 and 9. While SbNF-YB4, SbNF-YB5, SbNF-YB6 and SbNF-YB7 were observed on chromosome 3, SbNF-YB13, SbNF-YB14, SbNF-YB15 and SbNF-YB16 were noticed on chromosome 9. Three genes each SbNF-YB10, SbNF-YB11, SbNF-YB12 and SbNF-YB17, SbNF-YB18, SbNF-YB19 were seen on chromosomes 7 and 10 respectively. Chromosomes 2 and 4 have 2 genes each, i.e., SbNF-YB2 and SbNF-YB3 on 2, SbNF-YB8 and SbNF-YB9 on 4. Chromosome 1 has only SbNF-YB1 located on it. Majority of the SbNF-YC genes were located on chromosome 7, and it accommodates 4 genes (SbNF-YC8, SbNF-YC9, SbNF-YC10, and SbNF-YC11). While chromosome 1 accommodates three genes (SbNF-YC1, SbNF-YC2, and SbNF-YC3), chromosome 8 contains SbNF-YC12 and SbNF-YC13. One gene each SbNF-YC4, SbNF-YC5, SbNF-YC6, SbNF-YC7, SbNF-YC14, and SbNF-YC15 was noticed on chromosomes 2, 3, 5, 6, 9, 10 respectively (Fig 7 and Table 1).

<h2>Chromosomal location of <i>NF-Y genes in Sorghum</i>.</h2>
Fig. 7.

Chromosomal location of NF-Y genes in Sorghum.

Estimation of non-synonymous and synonymous substitution rates

The non-synonymous (dN) and synonymous (dS) substitution (dN/dS) rates were calculated for genes which showed duplication events within Sorghum as paralogs and between other genomes as orthologs. The 11 paralogs (S6 Fig) exhibited the dN/dS between 0.0010 (SbNF-YB18/SbNF-YB19)-93.9760 (SbNF-YC8/SbNF-YC12). Of the 11 paralogs, only 5 showed the purifying/stabilizing selection (<1), while the remaining exhibited positive/Darwinian selection (>1) (Table 2). The SbNF-YA orthologs exhibited dN/dS rates ranging from 0.4000 (Sorbi001G340200/Seita.9G367200) and 28.6781 (Sorbi001G486000/Zm2G000686). This indicates that 2 were following purifying selection and the remaining 4 positive selection (S5 Table). The dN/dS rates of SbNF-YB orthologs varied from 0.0031 (Sorbi001G338700/Seita. 9G365700) and 10.2227 (Sorbi010G119200/ZM2G167576). While majority of them (7) evolved through purifying selection, the remaining evolved by positive selection mechanism (S6 Table). The orthologs of SbNF-YCs showed dN/dS rates between 0.0797 (Sorbi002G241500/Seita. 2G247700) and 1.9761 (Sorbi007G063200/ZM2G311316) (S7 Table).

Promoter analysis

Analysis of promoter sequences revealed cis-acting elements such as ABA-responsive (ABRE), drought-responsive (DRE, DPBF, MYB and MYC), heat shock-responsive (HSE), and low temperature-responsive (LTR) elements. Aside, methyl jasmonic acid- (MeJA-RE), salicylic acid- (SARE), and defence-responsive elements (TC-rich repeats), associated with biotic stress were detected. Majority of the genes in SbNF-YA family have G-BOX (CACGTG) cis-acting elements (S8 Table). In the SbNF-YB family, G-BOX and Sp1 were observed as major cis-acting elements except in SbNF-YB18 (S9 Table). On the other hand, Skn-1 motif is the most common cis-acting element in SbNF-YC family. The least expressed cis-acting elements in SbNF-YC are AUXRR-CORE and GCN4 motif I (S10 Table).

Protein-protein interaction (PPI) prediction analysis

The PPI mapping of SbNF-Ys showed that a cohort of proteins involved in various cellular, metabolic and molecular pathways are associated with miRNA surveillance pathway, DNA replication, base excision, nucleotide excision repair pathway, purine and pyramidine metabolism. They interacted with core histones, calcineurins, kelch motifs, serine-threonine phosphatases, histone lysine N-methyl transferase, and metal-dependant phosphatase (S7 and S8 Figs).

In silico prediction of miRNA target sites

The SbNF-YAs exhibited different miRNA target sites such as sbi-miR169, sbi-miR5389, sbi-miR6225, sbi-miR5568, sbi-miR6220, sbi-miR5567 and sbi-miR6232, SbNF-YBs showed sbi-miR5565, sbi-miR5568, sbi-miR6232, sbi-miR6220, sbi-miR821, sbi-miR437, sbi-miR528, sbi-miR395, sbi-miR169, sbi-miR171 and sbi-miR172. Likewise, SbNF-YCs showed sbi-miR6232, sbi-miR6235, sbi-miR395, sbi-miR437, sbi-miR395, sbi-miR6220, sbi-miR6225, sbi-miR6227, sbi-miR5569, sbi-miR5568, sbi-miR164, sbi-miR156, sbi-miR160, sbi-miR6218, sbi-miR6230, and sbi-miR821. Interestingly, all of them are known to be associated with translation and cleavage events (S11, S12 and S13 Tables).

Gene expression analysis of SbNF-Ys in different tissues treated with diverse abiotic stresses

Expression of all the 42 NF-Y family of genes was studied at the transcriptional level in different tissues and abiotic stress conditions besides ABA, and the heat map is presented in Fig 6A and 6B. Over all, the gene expressions were higher in leaf tissues in comparison with stem and root (Fig 8A and S8 Table). Among the 8 SbNF-YA genes, NF-YA6 appeared to be associated with salt, drought (imposed by mannitol), cold and high temperature stresses. It was also strongly triggered by ABA. On the other hand, NF-YA1 was induced by multiple stresses like salt, mannitol and high temperature stresses, but NF-YA3 and NF-YA5 were expressed only under high temperature. While salt and high temperature influenced the expression of NF-YA7, NF-YA8 was upregulated by cold and high temperature stresses. Among the 19 SbNF-YB members, SbNF-YB7, B12, B15 and B16 were strongly induced by different stresses like salt, mannitol, ABA, cold and high temperature. In contrast, seven out of fifteen SbNF-YC members, YC1, YC3, YC4, YC7, YC8, YC9, and YC13 were expressed only under high temperature stress. Further, five members (YC6, YC11, YC12, YC14 and YC15) were upregulated by multiple stresses including ABA (Fig 8B and S14 Table).

Fig. 8.

Expression analysis of SbNF-Ys A), in root, stem and leaf tissues B) across diverse tissues under different abiotic stresses. (R; Root, S; Stem, L; Leaf, MR; Mannitol Root, SR; Salt root, HR; High temperature root, CR; Cold root, ABAR; Abscisic acid root, MS; Mannitol stem, SS; Salt stem, HS; High temperature stem, CS; Cold stem, ABAS; Abscisic acid stem, ML; Mannitol leaf, SL; Salt leaf, HL; High temperature leaf, CL; Cold leaf, ABAL; Abscisic acid leaf).


Identification and structural analysis of SbNF-Y genes

It is observed from this study that the number of genes are variable in each of the three distinct subfamilies of NF-Y (NF-YA, NF-YB and NF-YC) across different taxa. NF-Ys are evolutionarily conserved in eukaryotes and each subunit is encoded by a single gene in yeast and animals [4]. But, the same is encoded by a family of genes varying from 8 to 39 in plants. A total of 33 NF-Ys were reported in Brassica napus (14 NF-YA, 14 NF-YB, 5 NF-YC) [30]; while 39 in Setaria italica (10 NF-YA, 15 NF-YB and 14 NF-YC) [31]; whereas 68 in Glycine max (21 NF-YA, 32 NF-YB, 15 NF-YC) [32]. Multiple members of NF-Y subunits in plants reflect the redundancy and differentiated functions of these genes which need to be explored by expression profiling. Pereira et al. [35] identified 22 NF-Y genes in Citrus sinensis and C. clementina (6 NF-YA, 11 NF-YB and 5 NF-YC). Xu et al. [54] pointed out that such a low number of NF-Y genes found in Citrus genome could be due to the whole genome duplication events occurred in A. thaliana when compared to Citrus. It appears that subunit YB has more number of genes in comparison with YA and YC members in most of the species including the present report. Single genes, but with multiple splicing isoforms (that encode NF-Y subunits) is generally noticed in the yeast and mammals. Contrary to this, in higher plants, multigene families are noticed which can encode each subunit. Such multiple genes are vital for plant systems for tissue specific expressions at various stages of growth and development. Further, such a subunit combination can assist the plant systems in performing diverse roles during stress/development. Malaviya et al. [55] searched the plant transcription factor database (Plant TFDB, ( version 3.0 Jin et al. [56]), for identifying NF-Y genes in sorghum. They identified a total of 33 NF-Y transcription factors comprising 8 NF-YA, 11 NF-YB, and 14 NF-YC subunits in sorghum using Plant TFDB. In contrast, in the present study, a total of 42 NF-Y genes, among which 8 code for SbNF-YA, 19 for SbNF-YB and 15 for SbNF-YC subunits were identified. This discrepancy is because, in the present study, sorghum genome sequence available in the public domain has been searched.

Koralewski and Krutovsky [57] pointed that finding out exon-intron organization is crucial since it provides an insight into evolutionary relationships among genes and organisms. Malviya et al. [55] reported no introns in 18 of the genes (out of 33), and 5 of them have only one intron. They reported 2 in NF-YA3, 5 in NF-YA5, 4 in NF-YA6, 3 in NF-YA8, 4 in NF-YB1, 5 in NF-YC4 and 3 in YC7. Interestingly, in the present study, introns were absent in 12 out of 42 NF-Y TFs. While NF-YC18 contained 17 (highest number), YA2 16, YB11 15, YA8 13, YB16 11, YC10 10, YA7, YB15, YC2, YC7 and YC15 8 introns each. Only YB4, YB8, YB14, YB17, and YB18 (in all 5) contained one intron. Like in S. bicolor, single intron genes were also noted in Medicago truncatula which lead to alternative spliced variants of NF-YA1 [58]. Similarly, one intron in the 5′-UTRs of the NF-YA members was observed in A. thaliana, O. sativa, C. sinensis [28, 27, 35]. This suggests that such a post-transcriptional regulatory mechanism is retained among NF-YA genes. Single intron NF-Ys were not observed in SbNF-YC subtype in the present analysis. Chen et al. [59] reported that most of the NF-YB contained only one exon, and the genes from the same clade displayed a similar motif pattern in Gossypium hirsuyum. Chu et al. [60] reported 5 exons and 4 introns (6 genes) or 6 exons and 5 introns (2 genes) in CaNF-YA gene family members in Cicer arietinum. Further, they noticed 1 to 6 exons in CaNF-YB family, and 7 intronless out of 11 members in the CaNF-YC family. They reported 1 intron in NF-YC1, and 3 in NF-YC9. This suggests that a post-transcriptional regulatory mechanism is retained among NF-YA genes. Thus, the presence of multiple exon/intron gene organizations have been found in all the NF-Y family members in other species like B. napus [30], and S. lycopersicum [61] also. This infers that the presence of exon/intron is an attribute and typical of NF-Ys in higher plants. Loss or gain of spliceosomal introns led to the progress in our understanding of the molecular mechanisms associated with intron evolution and variation in gene function [62]. Fusion of exons and intron loss, might play a key role in the evolution of larger families like NF-Ys. Further, several members of the NF-YB and NF-YC have been found without any introns like in S. bicolor [55], Ricinus cummunis [34] and chickpea [60]. Introns are essential parts of all eukaryotic genes. In eukaryotic systems, introns are known to execute several functions like exon shuffling [63], gene expression alterations [64] and also tune the evolutionary rate of genes [62].

Motif identification and chromosomal localization of SbNF-Ys

NF-Y proteins display both conserved and non-conserved regions in Arabidopsis and others. Such conserved sequences may be vital for DNA interactions at CCAAT sites as pointed out by Siefers et al. [13], Romier et al. [2], and Testa et al. [65]. Hahn et al. [66] demonstrated that the yeast CCAAT box factor is a heteromer that contains HAP2 and HAP3 proteins. Xing et al. [67] showed that HAP2 is a 21 residue region with 3 histidines and arginines. Both SbNF-YB and SbNF-YC proteins have histone domains, but not SbNF-YAs. Besides, they also contain centromere kinetochore components and chromatin reorganizing domains. These residues are conserved in all the 8 SbNF-YA proteins (present study) as well as in Oryza sativa, and Triticum aestivum [23, 25]. Romier et al. [2] demonstrated that NF-YC/NF-YB sub-complex interacts through histone fold motifs. The role of the alpha C-helix of NF-YC appears to be vital for trimerization as well as a target for regulatory proteins like that of MYC and p53. It looks that heterotrimeric NF-Y proteins recognize the CCAAT regulatory elements represented in promoter and enhancer regions and modulate the genes. Steidl et al. [68], Liu and Howell [69] pointed out that NF-YB and NF-YC form a dimer in the cytoplasm and then translocated to the nucleus to interact with that of NF-YA to form a heterotrimer complex. Further, it has been demonstrated that bZIP28 and NF-Y transcription factors are activated by endoplasmic reticulum stress and assemble into a transcriptional complex to regulate downstream stress response genes in A. thaliana [69]. Alpha helix transmembrane spans with average hydrophobicity were predicted in 12 of the NF-Y proteins in S. bicolor. Anchoring of protein to glycosylphosphatidyl inositol (GPI) via the C-terminal attachment was predicted in three of the NF-YB proteins namely SbNF-YB1, YB3 and YB6. This appears rational since earlier gene fusion experiments conducted by Caras et al. [70] demonstrated that the C-terminal signal sequence has GPI-anchoring residues. The distribution of NF-Y genes appears to be widespread among different chromosomes. While OsHAP genes were dispersed on 11 out of the 12 rice chromosomes [26], in S. bicolor, they are distributed on 10 chromosomes.

Phylogenetic assessment, divergence and promoter analysis

Among the NF-YA family members, A2 and A5 appeared on the same clade indicating that they are closer to each other compared to others. While Malviya et al. [55] found that SbNF-YB8 was closer to SbNF-YA and SbNF-YC proteins, we could not observe such a correlation. On the contrary, B12 was observed closer to C11 and B13 to C1 in the present study than YA family members. It appears therefore YB and YC members might have close correlations in comparison with other members. Malviya et al. [55] noticed several ortholog and paralog groups through the phylogenetic analysis of SbNF-Y proteins along with 36 Arabidopsis and 28 rice NF-Y proteins. Malviya et al. [55] reported that Sorghum NF-Y family gene expansion is due to segmental duplication events. It appears that SbNF-Y genes retained their function even after duplication. Generally, gene family expansion occurs through segmental, tandem duplications, and transposition events [71]. In the present investigation, 11 paralogs were observed due to 3 regional duplications, and 8 segmental duplications, inferring that segmental duplications are responsible for SbNF-Y gene family expansion. Six duplication events were observed in SbNF-YB family, and this is a large number when compared to other subfamilies. SbNF-YB4, B5, B13 and B14 were phylogenetically distinct from other SbNF-YBs, and might have formed by recent duplications. SbNF-Ys exhibited 20 orthologous events with Zea, 7 with Setaria and 1 with Hordeum, which indicates their monocot ancestors. The synonymous (dS) and nonsynonymous (dN) substitutions reveal the selective pressure on duplicated genes. In the present study, phylogenetic relationships and ortholog predictions displayed that sorghum has additional NF-YB genes with unknown functions in comparison with Arabidopsis. The synonymous (dS) and nonsynonymous (dN) substitutions reveal the selective pressure on duplicated genes. Nekrutenko et al. [72] pointed out that greater than 1 dN/dS value represents positive selection, less than 1 functional constraint, and equal to 1 neutral selection. In the present study, it appears that majority of the duplicated genes evolved through purifying selection. The phytohormone-responsive cis-acting elements make the plants to tolerate various environmental changes. The ABRE play an important role in ABA signalling and abiotic stress tolerance. In the present investigation, large number of ABA-responsive elements were observed in majority of NF-Ys besides Skn elements that participate in endosperm expression [73]. Further, the presence of light-responsive elements like SP1, I-Box, and G-BOX indicate their roles in the regulation of gene responses to light. Interestingly, all the elements are rich with heat shock elements (HSE), which indicates their diverse roles in various stress response mechanisms.

miRNA analysis and protein-protein interactions

It is known that stress-responsive miRNAs target the transcription factors, which regulate the plant growth and development. The rationale behind finding out miRNA target sites is to know if any miRNAs associated in the regulation of SbNF-Ys exist in the genome. Identifying the target sites would subsequently help us in elucidating the regulation of SbNF-Ys during salt, drought and high temperature stress conditions. The miRNAs may also involve in gene networks regulated by transcription factors like NF-Ys. Identifying the interactions between miRNAs and transcription factors like NF-Ys will serve to screen their roles in stress tolerance, signal transduction, different developmental stages and synthesis of secondary metabolites, which will help to develop desired phenotypes with stress tolerance. While Fang et al. [74] reported targeting of the NAC mRNA by miRNA for abiotic stress responses, Stief et al. [75] noticed down regulation of heat stress memory by another miRNA. In the present investigation, miR169 identified was known to participate in post transcriptional regulation [76, 77, 19, 78]. Furthermore, miR169 and NF-YA5 knockout plants showed hypersensitivity to drought indicating their importance in drought tolerance [21]. Overexpression of miR169c in tomato enhanced the drought tolerance by reducing stomatal opening [79]. Therefore, in silico screening for miRNAs and their validation for abiotic stress response is highly important especially in the context of non-coding RNAs playing a gamut of regulatory roles. In addition, the PPI analysis revealed that they interact with calcineurins, the calcium sensors that usually confer spatial specificity in Ca2+ signalling, and play important roles in abiotic stress tolerance [80]. NF-Ys also participate in circadian clock and flowering time regulation, serine/threonine phosphatases and metal-dependant phosphatases and control the dephosphorylation of phosphoprotein substrates [81].

Gene expression analysis in different sorghum tissues under abiotic stress conditions

Analysis of NF-Y gene expressions by qRT-PCR indicated tissue-specific and stress-inducible expression profile. NF-YA5, A6, B7, B12, B15, B16, C6, C11, C12, C14 and C15 revealed significant differential expression patterns in response to the abiotic stresses in S. bicolor. Such a tissue-specific expression pattern was earlier noticed in several plants [30, 82]. This may indicate a sub-functionalization of different members in specific tissues under different abiotic stress conditions. Pereira et al. [35] pointed out that CsNF-YA2, CsNF-YB5/11 and CsNF-YC2/3 could form potential complexes in the citrus fruit. Many NF-Y genes were reported to be associated with both biotic and abiotic stresses. Xu et al. [83] reported high expression of BnNF-YA10 and BnNF-YB3, BnNF-YB7, BnNF-YB10 and BnNF-YB14 under NaCl stress. Under polyethylene glycol treatment, expression of BnNF-YA9, 10, 11 and 12 genes increased in B. napus. Malviya et al. [55] performed in silico gene expression analysis under abiotic stress conditions using rice transcriptome data. This revealed several of the sorghum NF-Y genes are associated with salt, drought, cold and temperature stresses. Since such an analysis is based on rice transcriptome database, this cannot give accurate results. But, in the present study, detailed gene expression studies were carried out and the results indicate that SbNF-YA1, 2, and 6 are upregulated under 200 mM salt and 200 mM mannitol stresses. NF-YA7 has been found associated with high temperature (40°C) stress, but NF-YA8 is triggered by both cold (4°C) and high temperature stresses. Among NF-YB genes, 7, 12, 15, and 16 are induced under multiple stress conditions such as salt, mannitol, ABA, cold and high temperatures. Likewise, NF-YC 6, 11, 12, 14, and 15 have been found enhanced significantly in a tissue specific manner under multiple abiotic stress conditions. Thus, the present analysis revealed that several of the NF-Ys are implicated in abiotic stresses and also modulated by ABA. Such a modulation of the NF-Ys by ABA was not shown by Malviya et al. [55]. Zhang et al. [84] found that many Physcomitrella patens NF-Y genes were responsive to abiotic stresses through ABA-dependent or independent pathways. In the present study, several genes were upregulated when treated with ABA, indicating that they are ABA-dependent. It has been observed from the present study that majority of the mannitol (drought)-inducible genes were also induced by salt, high temperature stresses and ABA. Few of the high temperature stress-induced genes are also induced by cold stress (NF-YA2, 4, 6, 8, NF-YB2, 7, 10, 11, 12, 14, 16, 17, NF-YC4, 6, 12, and 13). Seki et al. [85] noticed that drought-inducible genes are also inducible by salt stress and ABA treatments in A. thaliana. Ha et al. [86] observed that diverse transcription factor families modulate plant responses to abiotic stresses independent of ABA or dependent of ABA [87, 88]. Several members of the TFs also function in both ABA-dependent and independent ways [8991]. Interestingly, such a crosstalk can be achieved via indirect interactions between TFs and cis-elements present in the same promoter regions of the target genes [92].

Quach et al. [32] reported involvement of soybean NF-Y genes in specific developmental stages and also stress responses. In Prunus mume, Yang et al. [33] observed high expression of PmNF-YA1/2/4/5/6, PmNF-YB3/4/8/10/11/13, and PmNF-YC1/2/4/5/6/8 under osmotic stress and ABA. In citrus, CsNF-YA5 and CsNF-YB1/2/4/5/11 were found upregulated by drought stress [35]. Such a finding was proved later by overexpression of AtNF-YB1 in Arabidopsis and its ortholog ZmNF-YB2 in maize which showed enhanced drought tolerance [93]. Similarly, overexpression of osmotic and ABA-inducible NF-YB genes PwNF-YB3 from Picea and PdNF-YB7 from poplar in Arabidopsis exhibited improved drought tolerance activity [94, 95]. Transgenic rice plants harbouring bermudagrass NF-YC gene showed tolerance under drought [95]. NF-Y genes participate in stress tolerance mechanism by interacting with other stress inducible genes like antioxidants. The connection between NF-Ys and antioxidants was observed in previous reports; CsNF-YA5 [35], AtNF-YA5 interacts with glutathione S-transferase, peroxidases and an oxidoreductase [94] and SiNF-YA1 enhance the activity of superoxide dismutase, peroxidase and catalase [94]. Expression profiles exhibited by paralogous SbNF-Y genes in different tissues of sorghum under stress treatments suggest a clear functional redundancy among this gene family members. It is interesting to study how these NF-Ys regulate the expression of downstream genes that perform a wide spectrum of functions. Siefers et al. [13] pointed out that some transcription factors control gene expression by binding to cis-regulatory elements as individual subunits. But, it also appears that others are deployed in a combinatorial fashion both spatially and temporally.


Genome-wide screening revealed the existence of a total of 42 NF-Y genes (8 SbNF-YA, 19 SbNF-YB and 15 SbNF-YC subunit members) in Sorghum bicolor. In silico analysis of promoters revealed that they comprise many stress-related cis-elements such as ABRE and HSE indicating their role in salt, drought and high temperature stress responsiveness. The tissue specific expression of NF-Y transcription factors under salt, drought, ABA, cold and high temperature indicated their role in multiple stress tolerance. In view of this, we firmly believe that our studies have allowed identifying the candidate genes for further validation under an array of abiotic stress conditions in a crop species.

Compliance with ethical requirement

Authors do not have any other interests that influence the results and discussion of this paper. The authors have read the Journal’s policies and the authors of this paper have the following competing interests. RP is the President & CEO of Genomix Molecular Diagnostics Pvt Ltd., Kukatpally, Hyderabad, India. RP is the CEO of Genomix CARL Pvt. Ltd., Andhra Pradesh, India, but does not receive a salary in this capacity. RP is the President and CEO of Genomix Biotech Inc., 2620 Braithwood Road, Atlanta, GA 30345, USA, but does not receive any salary in this capacity. There are no patents, or products in development or marketed products associated with this research to declare. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Supporting information

S1 Fig [ppt]
1–10 MEME identified motif sequences of .

S2 Fig [ppt]
1–10 MEME identified motif sequences of .

S3 Fig [ppt]
1–10 MEME identified motif sequences of .

S4 Fig [ppt]
Distribution of 1–10 MEME identified SbNFY-A, B, and C conserved motifs.

S5 Fig [ppt]
1–10 MEME identified motif sequences of , , and .

S6 Fig [pptx]
Gene duplications of .

S7 Fig [ppt]
Network analysis of SbNF-Ys.

S8 Fig [ppt]
Functional partners of SbNFYs and their role in interaction network.

S1 Table [doc]
List of plants searched against .

S2 Table [doc]
gene specific primers used in the gene expression analysis.

S3 Table [xlsx]
The reliability of identified .

S4 Table [doc]
Number of phosphorylation sites in NFYs.

S5 Table [doc]
Non-synonymous to synonymous substitution ratios of orthologs.

S6 Table [doc]
Non-synonymous to synonymous substitution ratios of orthologs.

S7 Table [doc]
Non-synonymous to synonymous substitution ratios of orthologs.

S8 Table [doc]
Conserved -acting elements in promoters.

S9 Table [doc]
Conserved -acting elements in promoters.

S10 Table [doc]
Conserved -acting elements in promoters.

S11 Table [doc]
analysis of miRNAs for .

S12 Table [doc]
analysis of miRNAs for .

S13 Table [doc]
analysis of miRNAs for .

S14 Table [doc]
Native and relative expression values of .


1. Kim IS, Sinha S, De Crombrugghe B, Maity SN. Determination of functional domains in the C subunit of the CCAAT-binding factor (CBF) necessary for formation of a CBF-DNA complex: CBF-B interacts simultaneously with both the CBF-A and CBF-C subunits to form a heterotrimeric CBF molecule. Molecular and Cellular Biology. 1996;16(8): 4003–13. doi: 10.1128/mcb.16.8.4003 8754798

2. Romier C, Cocchiarella F, Mantovani R, Moras D. The NF-YB/NF-YC structure gives insight into DNA binding and transcription regulation by CCAAT factor NF-Y. Journal of Biological Chemistry. 2003; 278(2): 1336–45. doi: 10.1074/jbc.M209635200 12401788

3. Petroni K, Kumimoto RW, Gnesutta N, Calvenzani V, Fornari M, Tonelli C, et al. The promiscuous life of plant NUCLEAR FACTOR Y transcription factors. The Plant Cell. 2012; 24(12): 4777–92. doi: 10.1105/tpc.112.105734 23275578

4. Laloum T, De Mita S, Gamas P, Baudin M, Niebel A. CCAAT-box binding transcription factors in plants: Y so many?. Trends in Plant Science. 2013; 18(3): 157–66. doi: 10.1016/j.tplants.2012.07.004 22939172

5. Frontini M, Imbriano C, Manni I, Mantovani R. Cell-cycle regulation of NF-YC nuclear localization. Cell Cycle. 2004; 3(2): 205–10. doi: 10.4161/cc.3.2.654

6. Mantovani R. The molecular biology of the CCAAT-binding factor NF-Y. Gene. 1999; 239(1): 15–27. doi: 10.1016/s0378-1119(99)00368-6 10571030

7. Hackenberg D, Wu Y, Voigt A, Adams R, Schramm P, Grimm B. Studies on differential nuclear translocation mechanism and assembly of the three subunits of the Arabidopsis thaliana transcription factor NF-Y. Molecular Plant. 2012; 5(4): 876–88. doi: 10.1093/mp/ssr107 22199235

8. Kwong RW, Bui AQ, Lee H, Kwong LW, Fischer RL, Goldberg RB, et al. LEAFY COTYLEDON1-LIKE defines a class of regulators essential for embryo development. The Plant Cell. 2003; 15(1): 5–18. doi: 10.1105/tpc.006973 12509518

9. Warpeha KM, Upadhyay S, Yeh J, Adamiak J, Hawkins SI, Lapik YR, et al. The GCR1, GPA1, PRN1, NF-Y signal chain mediates both blue light and abscisic acid responses in Arabidopsis. Plant Physiology. 2007; 143(4): 1590–600. doi: 10.1104/pp.106.089904 17322342

10. Nelson DE, Repetti PP, Adams TR, Creelman RA, Wu J, Warner DC, et al. Plant nuclear factor Y (NF-Y) B subunits confer drought tolerance and lead to improved corn yields on water-limited acres. Proceedings of the National Academy of Sciences. 2007; 104(42): 16450–5.

11. Li L, Yu Y, Wei J, Huang G, Zhang D, Liu Y, et al. Homologous HAP5 subunit from Picea wilsonii improved tolerance to salt and decreased sensitivity to ABA in transformed Arabidopsis. Planta. 2013; 238(2): 345–56. doi: 10.1007/s00425-013-1894-0 23703145

12. Yan DH, Xia X, Yin W. NF-YB family genes identified in a poplar genome-wide analysis and expressed in Populus euphratica are responsive to drought stress. Plant Molecular Biology Reporter. 2013; 31(2): 363–70.

13. Siefers N, Dang KK, Kumimoto RW, Bynum WE, Tayrose G, Holt BF. Tissue-specific expression patterns of Arabidopsis NF-Y transcription factors suggest potential for extensive combinatorial complexity. Plant Physiology. 2009; 149(2): 625–41 doi: 10.1104/pp.108.130591 19019982

14. Ballif J, Endo S, Kotani M, MacAdam J, Wu Y. Over-expression of HAP3b enhances primary root elongation in Arabidopsis. Plant Physiology and Biochemistry. 2011; 49(6): 579–83. doi: 10.1016/j.plaphy.2011.01.013 21316979

15. Stephenson TJ, McIntyre CL, Collet C, Xue GP. TaNF-YB3 is involved in the regulation of photosynthesis genes in Triticum aestivum. Functional & Integrative Genomics. 2011; 11(2): 327–40.

16. Sun X, Ling S, Lu Z, Ouyang YD, Liu S, Yao J. OsNF-YB1, a rice endosperm-specific gene, is essential for cell proliferation in endosperm development. Gene. 2014; 551(2): 214–21. doi: 10.1016/j.gene.2014.08.059 25178525

17. Huang M, Hu Y, Liu X, Li Y, Hou X. Arabidopsis LEAFY COTYLEDON1 mediates postembryonic development via interacting with PHYTOCHROME-INTERACTING FACTOR4. The Plant Cell. 2015; 27(11): 3099–111. doi: 10.1105/tpc.15.00750 26566918

18. Soyano T, Kouchi H, Hirota A, Hayashi M. Nodule inception directly targets NF-Y subunit genes to regulate essential processes of root nodule development in Lotus japonicus. PLoS Genetics. 2013; 9(3): e1003352. doi: 10.1371/journal.pgen.1003352 23555278

19. Ni Z, Hu Z, Jiang Q, Zhang H. GmNFYA3, a target gene of miR169, is a positive regulator of plant tolerance to drought stress. Plant Molecular Biology. 2013; 82(1–2): 113–29. doi: 10.1007/s11103-013-0040-5 23483290

20. Alam MM, Tanaka T, Nakamura H, Ichikawa H, Kobayashi K, Yaeno T, et al. Overexpression of a rice heme activator protein gene (Os HAP 2E) confers resistance to pathogens, salinity and drought, and increases photosynthesis and tiller number. Plant Biotechnology Journal. 2015; 13(1): 85–96. doi: 10.1111/pbi.12239 25168932

21. Li WX, Oono Y, Zhu J, He XJ, Wu JM, Iida K, et al. The Arabidopsis NFYA5 transcription factor is regulated transcriptionally and post transcriptionally to promote drought resistance. The Plant Cell. 2008; 20(8): 2238–51. doi: 10.1105/tpc.108.059444 18682547

22. Han X, Tang S, An Y, Zheng DC, Xia XL, Yin WL. Overexpression of the poplar NF-YB7 transcription factor confers drought tolerance and improves water-use efficiency in Arabidopsis. Journal of Experimental Botany. 2013; 64(14): 4589–601. doi: 10.1093/jxb/ert262 24006421

23. Chen M, Zhao Y, Zhuo C, Lu S, Guo Z. Overexpression of a NF‐YC transcription factor from Bermuda grass confers tolerance to drought and salinity in transgenic rice. Plant Biotechnology Journal. 2015; 13(4): 482–91. doi: 10.1111/pbi.12270 25283804

24. Zhang F, Han M, Lv Q, Bao F, He Y. Identification and expression profile analysis of NUCLEAR FACTOR-Y families in Physcomitrella patens. Frontiers in Plant Science. 2015; 6: 642. doi: 10.3389/fpls.2015.00642 26347760

25. Stephenson TJ, McIntyre CL, Collet C, Xue GP. Genome-wide identification and expression analysis of the NF-Y family of transcription factors in Triticum aestivum. Plant Molecular Biology. 2007; 65(1–2): 77–92. doi: 10.1007/s11103-007-9200-9 17598077

26. Thirumurugan T, Ito Y, Kubo T, Serizawa A, Kurata N. Identification, characterization and interaction of HAP family genes in rice. Molecular Genetics and Genomics. 2008; 279(3): 279–89. doi: 10.1007/s00438-007-0312-3 18193457

27. Yang W, Lu Z, Xiong Y, Yao J. Genome-wide identification and co-expression network analysis of the OsNF-Y gene family in rice. The Crop Journal. 2017; 5(1): 21–31.

28. Zhao H, Wu D, Kong F, Lin K, Zhang H, Li G. The Arabidopsis thaliana nuclear factor Y transcription factors. Frontiers in Plant Science. 2017; 7:2045. doi: 10.3389/fpls.2016.02045 28119722

29. Cao S, Kumimoto RW, Siriwardana CL, Risinger JR, Holt BF III. Identification and characterization of NF-Y transcription factor families in the monocot model plant Brachypodium distachyon. PloS ONE. 2011; 6(6): e21805. doi: 10.1371/journal.pone.0021805 21738795

30. Liang M, Yin X, Lin Z, Zheng Q, Liu G, Zhao G. Identification and characterization of NF-Y transcription factor families in Canola (Brassica napus L.). Planta. 2014; 239(1): 107–26. doi: 10.1007/s00425-013-1964-3 24097262

31. Feng ZJ, He GH, Zheng WJ, Lu PP, Chen M, Gong YM, et al. Foxtail millet NF-Y families: genome-wide survey and evolution analyses identified two functional genes important in abiotic stresses. Frontiers in Plant Science. 2015; 6: 1142. doi: 10.3389/fpls.2015.01142 26734043

32. Quach TN, Nguyen HT, Valliyodan B, Joshi T, Xu D, Nguyen HT. Genome-wide expression analysis of soybean NF-Y genes reveals potential function in development and drought response. Molecular Genetics and Genomics. 2015; 290(3): 1095–115. doi: 10.1007/s00438-014-0978-2 25542200

33. Yang J, Wan XL, Guo C, Zhang JW, Bao MZ. Identification and expression analysis of nuclear factor Y families in Prunus mume under different abiotic stresses. Biologia Plantarum. 2016; 60(3): 419–26.

34. Wang Y, Xu W, Chen Z, Han B, Haque ME, Liu A. Gene structure, expression pattern and interaction of Nuclear Factor-Y family in castor bean (Ricinus communis). Planta. 2018; 247(3): 559–72. doi: 10.1007/s00425-017-2809-2 29119268

35. Pereira SL, Martins CP, Sousa AO, Camillo LR, Araujo CP, Alcantara GM, et al. Genome-wide characterization and expression analysis of citrus NUCLEAR FACTOR-Y (NF-Y) transcription factors identified a novel NF-YA gene involved in drought-stress response and tolerance. PloS ONE. 2018; 13(6): e0199187. doi: 10.1371/journal.pone.0199187 29906271

36. Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, et al. The Sorghum bicolor genome and the diversification of grasses. Nature. 2009; 457(7229): 551. doi: 10.1038/nature07723 19189423

37. Tuberosa R, Salvi S, Sanguineti MC, Landi P, Maccaferri M, Conti S. Mapping QTLs regulating morpho‐physiological traits and yield: Case studies, shortcomings and perspectives in drought‐stressed maize. Annals of Botany. 2002; 89(7): 941–63.

38. Mace ES, Tai S, Gilding EK, Li Y, Prentis PJ, Bian L, et al. Whole-genome sequencing reveals untapped genetic potential in Africa’s indigenous cereal crop sorghum. Nature Communications. 2013; 4: 2320. doi: 10.1038/ncomms3320 23982223

39. Letunic I, Copley RR, Schmidt S, Ciccarelli FD, Doerks T, Schultz J, et al. SMART 4.0: towards genomic data integration. Nucleic Acids Research. 2004; 32 (suppl_1): D142–4.

40. Guo AY, Zhu QH, Chen X, Luo JC. GSDS: a gene structure display server. Yi chuan = Hereditas. 2007; 29(8): 1023–6. 17681935

41. Bailey TL, Williams N, Misleh C, Li WW. MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Research. 2006; 34(suppl_2): W369–73.

42. Gasteiger E, Hoogland C, Gattiker A, Wilkins MR, Appel RD, Bairoch A. Protein identification and analysis tools on the ExPASy server. InThe proteomics protocols handbook 2005 (pp. 571–607). Humana press.

43. Blom N, Sicheritz‐Pontén T, Gupta R, Gammeltoft S, Brunak S. Prediction of post‐translational glycosylation and phosphorylation of proteins from the amino acid sequence. Proteomics. 2004; 4(6): 1633–49. doi: 10.1002/pmic.200300771 15174133

44. Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, et al. WoLF PSORT: protein localization predictor. Nucleic Acids Research. 2007; 35(suppl_2): W585–7.

45. Lescot M, Déhais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, et al. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Research. 2002; 30(1): 325–7. doi: 10.1093/nar/30.1.325 11752327

46. Dai X, Zhao PX. psRNATarget: a plant small RNA target analysis server. Nucleic Acids Research. 2011; 39(suppl_2): W155–9.

47. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Molecular Biology and Evolution. 2013; 30(12): 2725–9. doi: 10.1093/molbev/mst197 24132122

48. Gu Z, Cavalcanti A, Chen FC, Bouman P, Li WH. Extent of gene duplication in the genomes of Drosophila, nematode, and yeast. Molecular Biology and Evolution. 2002; 19(3): 256–62. doi: 10.1093/oxfordjournals.molbev.a004079 11861885

49. Yang S, Zhang X, Yue JX, Tian D, Chen JQ. Recent duplications dominate NBS-encoding gene expansion in two woody species. Molecular Genetics and Genomics. 2008; 280(3): 187–98. doi: 10.1007/s00438-008-0355-0 18563445

50. Suyama M, Torrents D, Bork P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Research. 2006; 34(suppl_2): W609–12. doi: 10.1093/nar/gkl315 16845082

51. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Research. 2014; 43(D1): D447–52.

52. Sudhakar Reddy P, Srinivas Reddy D, Sivasakthi K, Bhatnagar-Mathur P, Vadez V, Sharma KK. Evaluation of sorghum [Sorghum bicolor (L.)] reference genes in various tissues and under abiotic stress conditions for quantitative real-time PCR data normalization. Frontiers in Plant Science. 2016; 7: 529. doi: 10.3389/fpls.2016.00529 27200008

53. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C T method. Nature Protocols. 2008; 3(6): 1101. 18546601

54. Xu Q, Chen LL, Ruan X, Chen D, Zhu A, Chen C, et al. The draft genome of sweet orange (Citrus sinensis). Nature genetics. 2013; 45(1): 59. doi: 10.1038/ng.2472 23179022

55. Malviya N, Jaiswal P, Yadhav D. Genome-wide characterization of Nuclear Factor Y (NF-Y) gene family of sorghum [Sorghum bicolor (L.) Moench]: a bioinformatics approach. Physiology and Molecular Biology of Plants. 2016; 22(1): 33–49. doi: 10.1007/s12298-016-0349-z 27186017

56. Jin J, Zhang H, Kong L, Gao G, Luo J. PlantTFDB 3.0: a portal for the functional and evolutionary study of plant transcription factors. Nucleic acids research. 2013; 42(D1): D1182–7.

57. Koralewski TE, Krutovsky KV. Evolution of exon-intron structure and alternative splicing. PLoS ONE. 2011; 6(3): e18055. doi: 10.1371/journal.pone.0018055 21464961

58. Combier JP, Frugier F, De Billy F, Boualem A, El-Yahyaoui F, Moreau S, et al. MtHAP2-1 is a key transcriptional regulator of symbiotic nodule development regulated by microRNA169 in Medicago truncatula. Genes & Development. 2006; 20(22): 3084–8.

59. Li T, Zhang H, Liu Z, Deng H, Sharma S, Wei X, et al. A group of nuclear factor Y transcription factors are sub-functionalized during endosperm development in monocots. Journal of experimental botany. 2018; 69(10): 2495–510. doi: 10.1093/jxb/ery087 29514259

60. Chu H, Nguyen K, Watanabe Y, Le D, Pham T, Mochida K, et al. Identification, Structural Characterization and Gene Expression Analysis of Members of the Nuclear Factor-Y Family in Chickpea (Cicer arietinum L.) under Dehydration and Abscisic Acid Treatments. International journal of molecular sciences. 2018; 19(11): 3290.

61. Li S, Li K, Ju Z, Cao D, Fu D, Zhu H, et al. Genome-wide analysis of tomato NF-Y factors and their role in fruit ripening. BMC Genomics. 2016; 17(1): 36.

62. Roy SW, Gilbert W. The evolution of spliceosomal introns: patterns, puzzles and progress. Nature Reviews Genetics. 2006; 7(3): 211. doi: 10.1038/nrg1807 16485020

63. Fedorova L, Fedorov A. Introns in gene evolution. InOrigin and evolution of new gene functions. Springer, Dordrecht. 2003; 123–131.

64. Le Hir H, Nott A, Moore MJ. How introns influence and enhance eukaryotic gene expression. Trends in Biochemical Sciences. 2003; 28(4): 215–20. doi: 10.1016/S0968-0004(03)00052-5 12713906

65. Testa A, Donati G, Yan P, Romani F, Huang TH, Vigano MA, et al. Chromatin immunoprecipitation (ChIP) on chip experiments uncover a widespread distribution of NF-Y binding CCAAT sites outside of core promoters. Journal of Biological Chemistry. 2005; 280(14): 13606–15. doi: 10.1074/jbc.M414039200 15647281

66. Hahn ST, Pinkham JE, Wei R, Miller RE, Guarente LE. The HAP3 regulatory locus of Saccharomyces cerevisiae encodes divergent overlapping transcripts. Molecular and Cellular Biology. 1988; 8(2): 655–63. doi: 10.1128/mcb.8.2.655 2832732

67. Xing Y, Fikes JD, Guarente L. Mutations in yeast HAP2/HAP3 define a hybrid CCAAT box binding domain. The EMBO Journal. 1993; 12(12): 4647–55. 8223474

68. Steidl S, Tuncher A, Goda H, Guder C, Papadopoulou N, Kobayashi T, et al. A single subunit of a heterotrimeric CCAAT-binding complex carries a nuclear localization signal: piggy back transport of the pre-assembled complex to the nucleus. Journal of Molecular Biology. 2004; 342(2): 515–24. doi: 10.1016/j.jmb.2004.07.011 15327951

69. Liu JX, Howell SH. bZIP28 and NF-Y transcription factors are activated by ER stress and assemble into a transcriptional complex to regulate stress response genes in Arabidopsis. The Plant Cell. 2010; 22(3): 782–96. doi: 10.1105/tpc.109.072173 20207753

70. Caras IW, Weddell GN, Davitz MA, Nussenzweig V, Martin DW. Signal for attachment of a phospholipid membrane anchor in decay accelerating factor. Science. 1987; 238(4831): 1280–3. doi: 10.1126/science.2446389 2446389

71. Kong H, Landherr LL, Frohlich MW, Leebens‐Mack J, Ma H, DePamphilis CW. Patterns of gene duplication in the plant SKP1 gene family in angiosperms: evidence for multiple mechanisms of rapid gene birth. The Plant Journal. 2007; 50(5): 873–85. doi: 10.1111/j.1365-313X.2007.03097.x 17470057

72. Nekrutenko A, Baker RJ. Sub genome-specific markers in allopolyploid cotton Gossypium hirsutum: implications for evolutionary analysis of polyploids. Gene. 2003; 306: 99–103. doi: 10.1016/s0378-1119(03)00427-x 12657471

73. Washida H, Wu CY, Suzuki A, Yamanouchi U, Akihama T, Harada K, et al. Identification of cis-regulatory elements required for endosperm expression of the rice storage protein glutelin gene GluB-1. Plant Molecular Biology. 1999; 40(1): 1–2. doi: 10.1023/a:1026459229671 10394940

74. Fang Y, Xie K, Xiong L. Conserved miR164-targeted NAC genes negatively regulate drought resistance in rice. Journal of Experimental Botany. 2014; 65(8): 2119–35. doi: 10.1093/jxb/eru072 24604734

75. Stief A, Altmann S, Hoffmann K, Pant BD, Scheible WR, Baurle I. Arabidopsis miR156 regulates tolerance to recurring environmental stress through SPL transcription factors. The Plant Cell. 2014; 26(4): 1792–807. doi: 10.1105/tpc.114.123851 24769482

76. Zhao B, Ge L, Liang R, Li W, Ruan K, Lin H, et al. Members of miR-169 family are induced by high salinity and transiently inhibit the NF-YA transcription factor. BMC Molecular Biology. 2009; 10(1): 29.

77. Han X, Tang S, An Y, Zheng DC, Xia XL, Yin WL. Overexpression of the poplar NF-YB7 transcription factor confers drought tolaerance and improves water-use efficiency in Arabidopsis. Ournal of Experimental Botany 2013; 64: 4589–4601.

78. Sorin C, Declerck M, Christ A, Blein T, Ma L, Lelandais‐Briere C, et al. A miR169 isoform regulates specific NF‐YA targets and root architecture in Arabidopsis. New Phytologist. 2014; 202(4): 1197–211. doi: 10.1111/nph.12735 24533947

79. Zhang X, Zou Z, Gong P, Zhang J, Ziaf K, Li H, et al. Over-expression of microRNA169 confers enhanced drought tolerance to tomato. Biotechnology Letters. 2011; 33(2): 403–9. doi: 10.1007/s10529-010-0436-0 20960221

80. Batistic O, Kudla J. Plant calcineurin B-like proteins and their interacting protein kinases. Biochimica et Biophysica Acta (BBA)-Molecular Cell Research. 2009; 1793(6): 985–92.

81. Shi Y. Serine/threonine phosphatases: mechanism through structure. Cell. 2009; 139(3): 468–84. doi: 10.1016/j.cell.2009.10.006 19879837

82. Feng ZJ, He GH, Zheng WJ, Lu PP, Chen M, Gong YM, et al. Foxtail millet NF-Y families: genome-wide survey and evolution analyses identified two functional genes important in abiotic stresses. Frontiers in Plant Science. 2015; 6: 1142. doi: 10.3389/fpls.2015.01142 26734043

83. Xu L, Lin Z, Tao Q, Liang M, Zhao G, Yin X, et al. Multiple NUCLEAR FACTOR Y transcription factors respond to abiotic stress in Brassica napus L. PloS ONE. 2014; 9(10): e111354. doi: 10.1371/journal.pone.0111354 25356551

84. Zhang T, Zhang D, Liu Y, Luo C, Zhou Y, Zhang L. Overexpression of a NF-YB3 transcription factor from Picea wilsonii confers tolerance to salinity and drought stress in transformed Arabidopsis thaliana. Plant Physiology and Biochemistry. 2015; 94: 153–16 doi: 10.1016/j.plaphy.2015.05.001 26093308

85. Seki M, Narusaka M, Ishida J, Nanjo T, Fujita M, Oono Y, et al. Monitoring the expression profiles of 7000 Arabidopsis genes under drought, cold and high-salinity stresses using a full-length cDNA microarray. The PlantJournal. 2002; 31(3): 279–92.

86. Van Ha C, Esfahani MN, Watanabe Y, Tran UT, Sulieman S, Mochida K, et al. Genome-wide identification and expression analysis of the CaNAC family members in chickpea during development, dehydration and ABA treatments. PLoS One. 2014; 9(12):e114107. doi: 10.1371/journal.pone.0114107 25479253

87. Barbosa EG, Leite JP, Marin SR, Marinho JP, Carvalho JD, Fuganti-Pagliarini R, et al. Overexpression of the ABA-dependent AREB1 transcription factor from Arabidopsis thaliana improves soybean tolerance to water deficit. Plant molecular biology reporter. 2013; 31(3): 719–30.

88. Yan H, Jia H, Chen X, Hao L, An H, Guo X. The cotton WRKY transcription factor GhWRKY17 functions in drought and salt stress in transgenic Nicotiana benthamiana through ABA signaling and the modulation of reactive oxygen species production. Plant and Cell Physiology. 2014; 55(12):2060–76. doi: 10.1093/pcp/pcu133 25261532

89. Yang X, Yang YN, Xue LJ, Zou MJ, Liu JY, Chen F, et al. Rice ABI5-Like1 regulates abscisic acid and auxin responses by affecting the expression of ABRE-containing genes. Plant physiology. 2011; 156(3):1397–409. doi: 10.1104/pp.111.173427 21546455

90. Kim JS, Mizoi J, Yoshida T, Fujita Y, Nakajima J, Ohori T, et al. An ABRE promoter sequence is involved in osmotic stress-responsive expression of the DREB2A gene, which encodes a transcription factor regulating drought-inducible genes in Arabidopsis. Plant and Cell Physiology. 2011; 52(12): 2136–46. doi: 10.1093/pcp/pcr143 22025559

91. Zhang H, Zhu H, Pan Y, Yu Y, Luan S, Li L. A DTX/MATE-type transporter facilitates abscisic acid efflux and modulates ABA sensitivity and drought tolerance in Arabidopsis. Molecular plant. 2014; 7(10): 1522–32. doi: 10.1093/mp/ssu063 24851876

92. Lan Thi Hoang X, Du Nhi NH, Binh Anh Thu N, Phuong Thao N, Phan Tran LS. Transcription factors and their roles in signal transduction in plants under abiotic stresses. Current genomics. 2017; 18(6):483–97. doi: 10.2174/1389202918666170227150057 29204078

93. Nelson DE, Repetti PP, Adams TR, Creelman RA, Wu J, Warner DC, et al. Plant nuclear factor Y (NF-Y) B subunits confer drought tolerance and lead to improved corn yields on water-limited acres. Proceedings of the National Academy of Science, USA. 2007; 104: 16450–16455.

94. Li WX, Oono Y, Zhu J, He XJ, Wu JM, Iida K, et al. The Arabidopsis NFYA5 transcription factor is regulated transcriptionally and posttranscriptionally to promote drought resistance. Plant Cell. 2008; 20: 2238–2251. doi: 10.1105/tpc.108.059444 18682547

95. Chen M, Zhao Y, Zhuo C, Lu S, Guo Z. Overexpression of a NF-YC transcription factor from bermudagrass confers tolerance to drought and salinity in transgenic rice. Plant Biotechnology Journal. 2015, 13: 482–491. doi: 10.1111/pbi.12270 25283804

Článok vyšiel v časopise


2019 Číslo 9