Next Article in Journal
Temporal Variation in Community Composition of Root Associated Endophytic Fungi and Carbon and Nitrogen Stable Isotope Abundance in Two Bletilla Species (Orchidaceae)
Previous Article in Journal
Assessment of Genetic Diversity and Symbiotic Efficiency of Selected Rhizobia Strains Nodulating Lentil (Lens culinaris Medik.)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 10
Issue 1
10.3390/plants10010016
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic Organization, Phylogenetic Comparison, and Differential Expression of the Nuclear Factor-Y Gene Family in Apple (Malus Domestica)

by
Yanjie Qu
1,2,†,
Yaping Wang
1,2,†,
Jun Zhu
1,2,
Yugang Zhang
1,2 and
Hongmin Hou
1,2,*
1
College of Horticulture, Qingdao Agricultural University, Qingdao 266109, Shandong, China
2
Qingdao Key Laboratory of Genetic Development and Breeding in Horticultural Plants, Qingdao Agricultural University, Qingdao 266109, Shandong, China
*
Author to whom correspondence should be addressed.
†
These two authors contributed equally to this work.
Plants 2021, 10(1), 16; https://doi.org/10.3390/plants10010016
Submission received: 28 November 2020 / Revised: 16 December 2020 / Accepted: 22 December 2020 / Published: 24 December 2020
(This article belongs to the Section Plant Genetics, Genomics and Biotechnology)
Browse Figures
Versions Notes

Abstract

:
The nuclear factor Y (NF-Y) as a transcription factor plays an important role in plants growth and development, and response to stress. However, few genome-wide analyzes and functional research of the NF-Y family has been undertaken in apple (Malus domestica Borkh.) so far. In this study, we comprehensively identified the 43 MdNF-Y genes in apple, which dispersedly distributed among the three subgroups based on their sequence alignment analysis, including 11 MdNF-YAs, 22 MdNF-YBs and 10 MdNF-YCs. The members in the same subgroups had similar evolution relationships, gene structures, and conserved motifs. The gene duplication analysis suggested that all the genes were dispersed followed by 27 segmental duplication. Moreover, based on synteny analysis of MdNF-Ys with eight plant species results suggested that some ortholog genes were preserved during the evolution of these species. Cis-element analysis showed potential functions of MdNF-Ys in apple growth and development and responded to abiotic stress. Furthermore, the interaction among MdNF-Ys protein were investigated in yeast two-hybrid assays. The expression patterns of MdNF-Ys in tissue-specific response reveled divergence and might play important role in apple growth and development. Subsequently, whole MdNF-Y genes family was carried out for RT-PCR in response to five abiotic stress (ABA, drought, heat, cold, and salinity) to identify their expression patterns. Taken together, our study will provide a foundation for the further study to the molecular mechanism of apple in growing development and response to abiotic stresses.

1. Introduction

Transcription factors (TFs) control the transcription or expression of downstream target genes by interacting with cis-elements through covalent binding to the DNA binding domain. Nuclear factor Y (NF-Y) TFs, known as CCAAT-binding factors (CBFs) or heme activator proteins (HAPs), play a critical regulatory role in plant vital movement by binding to the CCAAT element. NF-Y consists of three distinct subunits including NF-YA (also known as CBF-B/HAP2), NF-YB (CBF-A/HAP3), and NF-YC (CBF-C/HAP5) [1]. NF-Y is ubiquitously expressed in most eukaryotes and in mammals and yeast, each subunit is encoded by only one or two NF-Y genes. However, in plants, each NF-Y subunit has evolutionarily formed relatively large gene families expressed from multiple NF-Y genes. For example, in the model plant Arabidopsis thaliana, a total of 30 AtNF-Y genes exist including 10 AtNF-YAs, 10 AtNF-YBs, and 10 AtNF-YCs [2,3]. In addition, the three subunits each possess conserved DNA-binding domains and mutual interaction domains to form heterotrimeric complexes. It is also worth noting that NF-YA and NF-YC family members have a nuclear localization signal (NLS), which NF-YB members generally lack [2]. The NF-YB proteins are first translocated from the cytoplasm to the nucleus as part of a dimer complex with NF-YC. Subsequently, the dimer complex combines with NF-YA to form a mature heterotrimeric NF-Y complex transcription factor that binds to CCAAT boxes in the promoters of target genes [4,5]. Moreover, some of the NF-Y subunits also interact with other proteins besides the NF-YA subunit to regulate plant growth. For instance, AtNF-YB3 and AtNF-YC2 interact with post-proteolytic bZIP28 and AtNF-YA4, respectively, to assemble a complex that regulates the response to endoplasmic reticulum (ER) stress in Arabidopsis [6]. AtNF-YC9 is a positive regulator involved in abscisic acid (ABA) signaling by interacting with ABA-responsive bZIP transcription factor ABA-INSENSITIVE5(ABI5) [7,8]. The AtNF-YC4 gene interacts with the Qua-Quine Starch (QQS) orphan gene to regulate carbon and nitrogen allocation and reduce susceptibility to pathogens and pests [9,10].
Since the discovery of the NF-Y family in Arabidopsis, orthologous genes have been isolated and identified at genome-wide levels from other crops, including but not limited to walnut (Juglans regia L.) (17 JrNF-YAs, 9 JrNF-YBs, and 7 JrNF-YCs) [11], citrus (Citrus sinensis) (6 CsNF-YAs, 11 CsNF-YBs, and 5 CsNF-YCs) [12], chickpea (Cicer arietinum L.) (8 CaNF-YAs, 21 CaNF-YBs, and 11 CaNF-YCs) [13], castor bean (Ricinus communis) (6 RcNF-YAs, 12 RcNF-YBs, and 7 RcNF-YCs) [14], grape (Vitis vinifera L.) (8 VvNF-YAs, 18 VvNF-YBs, and 8 VvNF-YCs) [15], and peach (Prunus persica L.) (6 PpNF-YAs, 12 PpNF-YBs, and 6 PpNF-YCs) [16]. Many reports have demonstrated that these NF-Y genes play various vital functions in plant development and stress resistance, including seed germination [17], endosperm development [18], flowering time [19,20], photosynthesis [21], root nodule formation [10,22], hormone response [23], ER stress response, and abiotic and biotic stress response [6,24].
Abiotic stresses, such as osmotic stresses (drought or salinity) and temperature stresses (cold or heat), are increasingly becoming environmental factors that are limiting fruit productivity and quality worldwide. Previous reports have revealed that the NF-Y genes are involved in regulating abiotic stress via ABA-dependent or independent pathways. In Arabidopsis, overexpression of AtNF-YA5, AtNF-YB1, or AtNF-YB2 has been shown to enhance drought stress tolerance by regulating the expression of stress-responsive genes [25,26,27]. Similarly, overexpression of OsNF-YA7 has been shown to improve drought tolerance of transgenic rice, via an ABA-independent pathway [24]. The heterologous overexpression of CdtNF-YC1 from hybrid bermudagrass improved transgenic rice tolerance to drought and salt stress [28]. Additionally, ZmNF-YB2 was demonstrated to have a significant role in transgenic maize drought stress tolerance [27,29]. Furthermore, ZmNF-YA3 also elevated transgenic plant tolerance to heat and drought stress by binding to ZmMYC4, ZmbHLH92, and ZmFAMA in the promoter region [29]. Overexpression of PdNF-YB21 increased the drought resistance of transgenic poplars by positively regulating root growth and enlarging xylem vessels [30]. Recently, some studies have demonstrated that miRNA responded to abiotic stresses by regulating the expression of their target gene NF-Y [31]. The overexpression of soybean miR169c in transgenic Arabidopsis conferred increased drought stress sensitivity via inhibiting the expression of its target AtNF-YA gene [32]. The expression of most of Zma-miR169 genes and their target ZmNF-YA genes were diverse changes during drought, salinity and ABA treatments [33,34].
As one of the most important economic crops, apples are widely cultivated in the temperate regions of the world. Unfortunately, with global climate change, apple trees also face many abiotic stresses during their growth and development. However, the role of NF-Y genes in stress tolerance in apple trees remains elusive. Based on the complete apple (Malus domestica Borkh.) genome released in 2010, the possibility of investigating the MdNF-Y gene family was offered in the species [35]. In this work, we analyzed the members of the apple NF-Y gene family (MdNF-Y) and determined their chromosomal location and detailed genetic information. Further, we characterized the apple NF-Y protein sequences, including the construction of a phylogenetic tree, detection of conserved motifs, and gene structure and synteny analysis. Then, to assess the function of MdNF-Y genes, we analyzed the cis-regulatory elements in the promoters and their transcriptional expression in different tissues. Yeast two-hybrid (Y2H) assays were conducted to study the interaction between various MdNF-Y subunits. The transcription profiles of the MdNF-Y genes were detected under various abiotic stresses. Our results provide a foundation for further study of the functional and regulatory mechanisms controlled by the MdNF-Y gene family.

2. Results

2.1. Identification and Characterization of NF-Y Family Genes in Apple

Following the removal of redundant sequences, we initially identified 11 MdNF-YAs, 26 MdNF-YBs, and 12 MdNF-YCs through Hidden Markov Model (HMM) analysis of the M. domestica genome (https://www.rosaceae.org/). Four putative MdNF-YBs (MD03G1280100, MD06G1209300, MD11G1164600, and MD14G1219800) and two putative MdNF-YCs (MD02G1273400 and MD07G1042300) did not contain the core structure of the NF-Y domain, so we removed them from further analysis. In order to facilitate further research, these genes were renamed according to their locations in chromosomes and the Arabidopsis nomenclature [2,36] (Table 1). Finally, the MdNF-YA, MdNF-YB, and MdNF-YC subunit families contained 11 (MdNF-YA1- A11), 22 (MdNF-YB1-YB22), and 10 (MdNF-YC1 -YC10) members, respectively.
The gene length, coding sequence (CDS), and protein sequence of MdNF-Ys ranged from 365 bp (MdNF-YC3) to 9824 bp (MdNF-YB17), from 207 bp (MdNF-YA4) to 1077 bp (MdNF-YA3), and from 68 (MdNF-YA4) to 457 (MdNF-YB19) amino acids (aa), respectively. The isoelectric points (pI), and the molecular weights (MW) of MdNF-Y proteins varied from 5.15 (MdNF-YC4) to 9.63 (MdNF-YA10), and from 7449.24 Da (MdNF-YA4) to 38,855.31 Da (MdNF-YA3), respectively. The detailed information of the MdNF-Y genes including protein and CDS lengths is shown in Table S1.
In addition, the three NF-Y subunits contained highly conserved domains among the previously reported species such as Arabidopsis (A. thaliana), grape (V. vinifera), and orange (C. sinensis) [2,12,15]. Similarly, multiple alignments indicated that apple NF-Y subunits also contained the same highly conserved regions as Arabidopsis, grape, and orange NF-Y subunits (Figure 1). The conserved central domain of the MdNF-YA, MdNF-YB, and MdNF-YC subunits, respectively, consisted of about 54, 90, and 81 aas, and three MdNF-Y subunits all contained a DNA binding domain. Moreover, the MdNF-YA contained an NF-YB/C interaction domain, and the MdNF-YB and MdNF-YC both contained an NF-YA interaction domain (Figure 1). These domains were necessary for NF-YA, NF-YB, and NF-YC to form a heterotrimeric complex to bind to CCAAT boxes [37]. Furthermore, the MdNF-YB and MdNF-YC both contained a histone-fold motif (HFM) resembling the core of H2B and H2A, respectively.

2.2. Phylogenetic Tree, Conserved Motifs, and Gene Structure Analysis of MdNF-Y Family

To deduce the potential evolutionary diversity and relationship of MdNF-Y subunit genes, we constructed a phylogenetic tree and analyzed the conserved motifs and gene structure for each MdNF-Y subunit (Figure 2). For each MdNF-Y subunit the corresponding genes had close phylogenetic relationships and consistently displayed conserved motifs and gene structures. For example, MdNF-YA9 and MdNF-YA10 were closely related evolutionarily and they presented similar motifs type (motif 3, 4, 7, 12, and 19) (Figure 2A). On the other hand, differences between various MdNF-Y subunits were evident in terms of major motifs. For example, motifs 3, 4, 7, 8, and 12 were contained in most MdNF-YA genes whereas MdNF-YA9 and MdNF-YA10 lacked motif 8. Most of the MdNF-YB genes contained motif 1, 2, and 3 except MdNF-YB4, MdNF-YB14, and MdNF-YB17. Most of the MdNF-YC genes contained seven types of motifs (motif 1, 4, 5, 6, 9, 10, and 14). However, MdNF-YC3 and MdNF-YC9 presented different motif types compared with the other MdNF-YC genes. In addition, the gene structures of the various MdNF-Y subunits also displayed differences in the distribution of exons and introns, whereas the gene structures within each MdNF-Y subunit were relatively similar (Figure 2B). For example, in the MdNF-YA subfamily, the CDS were separated by four introns except MdNF-YA4, and all MdNF-YCs contained a long CDS region which is not separated by introns.

2.3. Chromosome Distribution and Synteny Analysis of MdNF-Y Family Genes

To understand the chromosomal distribution of the different MdNF-Ys, chromosomal location map was created (Figure 3). All MdNF-Y genes were unevenly distributed on 16 of the 17 apple chromosomes except chr08. The number of genes on the chromosomes that contained MdNF-Y genes varied from one to four. Six chromosomes with four MdNF-Y genes, including chr02, chr03, chr05, chr11, chr12, and chr15, while only MdNF-YB1 was located on the chr01. In addition, two MdNF-Ys were located on chr06, chr07, chr09, chr13, and chr14, and three MdNF-Ys were located on the terminus of chr04 and chr10, respectively. Distribution of MdNF-Ys was concentrated on the terminus of several chromosomes including chr03, chr04, chr05, chr10, and chr11.
To investigate the gene family expansion mechanism of the MdNF-Y genes, we analyzed their synteny relationships in apple genomes (Figure 4). The results showed that no tandem duplication event had occurred but that 27 pairs of gene segmental duplication events could be identified (Figure 4, Table S2). Overall, 11 MdNF-YA, 15 MdNF-YB, and 8 MdNF-YC genes were mapped to the 16 chromosomes except for chr 08, while the pairs of paralogous genes were 7, 14, and 6, respectively. These results suggested that segmental duplications were the cause of MdNF-Y genes amplification. Interestingly, we found that there was respectively one triangular relationship in the MdNF-YA (MdNF-YA1, MdNF-YA6, and MdNF-YA11) and MdNF-YC subunits (MdNF-YC1, MdNF-YC6, and MdNF-YC10). There were four triangular relationships in MdNF-YB subunit including MdNF-YB1, MdNF-YB11, and MdNF-YB17; MdNF-YB6, MdNF-YB16, and MdNF-YB19; MdNF-YB7, MdNF-YB11, and MdNF-YB17; and MdNF-YB8, MdNF-YB18, and MdNF-YB19. In addition, more gene pairs are listed in Table S2.
To further study the gene expansion relationship and evolution of NF-Y subunits (Figure 5), we chose eight representative plant models with widely ranging homologies to analyze the synteny of NF-Y genes. The eight species contained six dicots, including three Rosacea species (Pyrus betulifolia, Prunus persica, and Fragaria vesca), Vitis vinifera, A. thaliana, and Brassica rapa and two monocots (Oryza sativa L. and Zea Mays L). The results suggested that many NF-Y genes in apple have homology to reference plants. It is well known that apples also belong to the Rosaceae, like F. vesca, P. betulifolia, and P. persica. Furtherly, many ortholog pairs of MdNF-Y genes were found among P. betulifolia (59 orthologous gene pairs distributed on all apple chr except chr8), P. persica (37 orthologous gene pairs distributed on all apple chr except chr8 and chr17), and F. vesca (37 orthologous gene pairs distributed on all apple chr except chr8), respectively. Further, with the exception of MdNF-YA1 and MdNF-YB10, the majority of MdNF-Y subunits had orthologous pairs in pears, which indicated that the NF-Y transcription factor families are highly homologous in apples and pears. However, fewer homologous gene pairs were observed between apple and O. sativa L. (only seven) and between apple and Z. Mays. (only one) (Table S3). In addition, we found some highly homologous genes were preserved during the species evolution. For example, MdNF-YA6 had homologous pairs in seven species excepted maize, and MdNF-YB19 had homologous pairs in six species excepted Arabidopsis and B. rapa. Likewise, MdNF-YC1, MdNF-YC6, and MdNF-YC8 had homologous pairs in five species, respectively.

2.4. The Cis-Acting Regulatory Members in the Promoter of MdNF-Y Family Genes

To predict the potential function of MdNF-Y transcription factors, we chose the 1500 bp upstream sequences using plantCARE to analyze the type of cis-elements in the promoter (Figure 6 and Table S5). We found a preponderant number of TATA-boxes and CAAT-boxes, which have been analyzed in terms of their roles in transcription [38]. The results also revealed the presence of a very large number of light cis-elements. Therefore, we classified the results into four broad categories: Those participated in plant development, phytohormones, abiotic, and biotic stress-responsive, and light-responsive. The cis-elements participated in light responsiveness included a 3-AF1 binding site, G-Box, GA-motif, and a GATA-motif. The cis-elements involved in plant development included flavonoid biosynthetic genes regulation (MBSI), endosperm expression (GCN4), meristem expression (CAT-box), circadian control (Circadian), Seed (RY-element), and root-specific (motif I). The cis-elements partake in abiotic and biotic stress included anaerobic induction (ARE), heat stress responsiveness (HSE), low-temperature responsiveness (LTR), anoxic specific inducibility (GC), drought-inducibility (MBS), defense and stress responsiveness (TC-rich repeats), auxin-responsive element (TGA-element), pathogen (W-box), and wound-responsive element (WUN-motif). The cis-elements involved in phytohormone responsive included abscisic acid (ABA) responsiveness (ABRE) MeJA-responsiveness (CGTCA-motif), ethylene (ERE), and salicylic acid (SA) responsiveness (TCA-element).
In addition, most of the MdNF-Ys had ABRE cis-acting regulatory elements with a total of 169 elements, indicating that MdNF-Y genes play an important role in ABA response. All MdNF-Y genes contained more or less different types of light-responsive elements, although all MdNF-YBs possessed a G-Box element. A number of cis-acting regulatory elements involved in plant development were found in some of the MdNF-Y genes. For example, motif I, which was involved in plant root development [39], was only found in MdNF-YC6. A circadian cis-acting element, involved in circadian control [40], was only discovered in MdNF-YC4. An RY-element, participated in plant seed-specific regulation [41], was only detected in MdNF-YB13. Therefore, these MdNF-Y genes need to be further explored since they may play a critical role in apple growth. The specific cis-element analysis of the MdNF-Y genes is shown in Table S5. Actually, some reports in other plants have indicated that NF-Ys participated in response to abiotic stresses (cold [42], heat [25], drought [43], and salinity [28], root development [44], seed-specific regulation [44] and photoperiod-dependent flowering [29]. So, we further detected the expression levels of MdNF-Y genes in different tissues and under various abiotic stresses in subsequent experiments.

2.5. Protein Interaction Analysis of MdNF-Y Genes

NF-YB, which does not contain a nuclear localization signal, needs to form a tight dimer with NF-YC in order to translocate from the cytoplasm to the nucleus, and subsequently binds to NF-YA to form a heterotrimeric complex. Then, the complex can interact with other regulatory factors to activate or repress the expression of downstream genes [45,46]. To determine the potential interactions between MdNF-YB and MdNF-YC, or between MdNF-YB and MdNF-YA members, we detected a physical interaction of some MdNF-Y gene products including five MdNF-YBs, three MdNF-YCs, and one MdNF-YA in the Y2H system (Figure 7). The full CDS of five MdNF-YBs were fused with the activation domain (AD) of the pGADT7 vector (expressing the “prey”), while the CDS of three MdNF-YCs and one MdNF-YA were fused to the DNA-binding domain (BD) of the pGBKT7 vector (expressing the “bait”). After pairwise co-transforming AD- and BD-expressing vectors into the Y2H Gold yeast strain, almost all yeast cells bearing both the MdNF-YB and MdNF-YC-expressing plasmids, or both MdNF-YB and MdNF-YA-expressing plasmids (except yeast cells bearing both MdNF-YB11 and MdNF-YA5) were capable of growth on SD-Leu/-Trp/-His medium indicating an interaction between prey and bait. Then, the surviving yeast cells were transferred to SD-Leu/-Trp/-His/-Ade medium. MdNF-YC5 showed strong interaction signals with MdNF-YB7, 11, and 17. MdNF-YC8 showed strong interaction signals with MdNF-YB1 and MdNF-YB8. MdNF-YB1 showed weak interaction signals with MdNF-YA5 and MdNF-YC4, and 5. MdNF-YC4 showed no interaction signals with MdNF-YB8, 11, or 17. MdNF-YA5 showed no interaction signals with MdNF-YB11 and 17. Taken together, these results suggested that many MdNF-YBs could interact with MdNF-YAs or three MdNF-YC genes on SD-Leu/-Trp/-His/-Ade medium.

2.6. Transcript Profiles Analysis of MdNF-Y Family Genes in Different Apple Tissues

To initially understand the function of the apple MdNF-Y genes, we analyzed the transcript profiles of the MdNF-Y gene family in roots, stems, leaves, receptacles, peel, sarcocarps, young fruits, and seeds using qRT-PCR. In general, the 43 MdNF-Y genes show diverse tissue-specific expression patterns and spatiotemporal expression characteristics (Figure 8). For the MdNF-YA subfamily, many genes had higher expression levels in vegetative and reproductive organs. For example, MdNF-YA3, 4, 7, 8 were strongly expressed in both young fruit and leaves. Moreover, MdNF-YA6 was highly expressed in roots, stem, and leaves, and MdNF-YA9 were highly expressed in roots. For the MdNF-YB subfamily, MdNF-YB1, 9,11,12,17,18 was highly expressed in sarcocarps, and MdNF-YB2, 4, 6, 13, 14, 16 were highly expressed in peel. MdNF-YB8, 18, 19 were highly expressed in root tissue. In addition, all MdNF-YC genes had higher expression levels in at least one reproductive organ, such as receptacles, peel, sarcocarps, young fruits, and seeds. However, it is worth noting that MdNF-Y genes sharing very high sequence and exon–intron structure similarity in duplicated genomic regions (Figure 2), exhibited similar expression patterns (Figure 8). For example, MdNF-YB1, MdNF-YB11, and MdNF-YB17 located in the duplicated genomic regions, were all highly expressed in sarcocarp. On the other hand, MdNF-YB8, MdNF-YB18, and MdNF-YB19 were highly expressed in root tissue. Overall, the overlapping but distinct expression patterns of MdNF-Y genes indicated that the MdNF-Y family plays a critical role in different growth and development stages of apples.

2.7. Expression Levels of MdNF-Y Genes under Different Abiotic Stresses

Previous studies have indicated that NF-Y not only regulates plant growing development but also responds to abiotic stresses [47]. Similarly, in our study, a very high number of abiotic and biotic stress-responsive elements were detected in the upstream promoter of MdNF-Y genes (Figure 6). Therefore, we used qRT-PCR to investigate the response of the MdNF-Y gene family to abiotic stress. The transcriptional profiles of MdNF-Y genes under abiotic stresses from 0 to 24 h were monitored in this study (Figure 9). Under low-temperature treatment, the expression of all MdNF-YAs and MdNF-YCs was upregulated and the expression levels of MdNF-YA genes were higher than those of MdNF-YC genes as a whole (Figure 9). In MdNF-YB gene family, MdNF-YB3, MdNF-YB5, and MdNF-YB9 gene expression levels reached their peak at 24 h whereas the expression of the rest of the MdNF-YB genes reached their peak at 12 h. The expression of MdNF-YB3 decreased at first but had recovered at 24 h. Under high-temperature treatment, the expression levels of most MdNF-Y genes were upregulated except MdNF-YC7 and MdNF-YC10 that were downregulated (Figure 9C). Interestingly, of the MdNF-YA subfamily, MdNF-YA3, MdNF-YA5, and MdNF-YA8 were highly expressed under low- or high-temperature treatment. Under drought treatment, the MdNF-Y gene family members were all upregulated showing a strong response to water deficit (Figure 9). In MdNF-YB gene family, MdNF-YB10 and MdNF-YB13 genes were first downregulated after PEG treatment and recovered at 24 h post-treatment. Under salinity treatment, the expression levels of MdNF-YB2, MdNF-YB3, MdNF-YB4, MdNF-YB13, and MdNF-YB14 were extremely low (Figure 9B). Under ABA treatment, the expression levels of all MdNF-YAs were upregulated (Figure 9A). Different MdNF-YA genes reached their peak at different post-treatment times, including MdNF-YA3 at 6 h, MdNF-YA1, MdNF-YA4, and MdNF-YA7 at 24 h, and the rest of MdNF-YA genes at 12 h. The expression levels of MdNF-YB9, MdNF-YB10, MdNF-YB15, MdNF-YB17, and MdNF-YB18 were decreased at 6 h post-treatment but had recovered at 12 h post-treatment, while the rest of MdNF-YB genes was upregulated. The expression of 59% of the MdNF-YB genes reached their peak at 24 h. However, the expression of all MdNF-YC genes reached its peak at 24 h. With the exception of the MdNF-YC4 and MdNF-YC6 genes, the expression of other genes was first downregulated but then upregulated.

3. Discussion

3.1. Conservation, Evolutionary and Divergence of the MdNF-Y Gene Family in Apple

Studies on the NF-Y genes in plant species have been accumulating since the function and regulatory mechanism of the first plant NF-Y gene was identified [47,48,49]. To date, NF-Y genes have been identified from simple model plants to more complex plants, and from monocotyledons to dicotyledons [2,15,16,50,51], such as Arabidopsis, maize, rice, grapes, rubber trees, and so on. However, there is few reports concerning the NF-Y gene family in apple. The number of NF-Y genes identified from plants varied from 13 in tomato [52] to 50 in maize [53]. In this study, we first systematically identified and analyzed the NF-Y gene family in apple at the genomic level and discovered 43 MdNF-Y genes in the apple genome, which was a greater number than for most of the other plants investigated so far.
Duplications at gene, chromosomal, or entire genomic level have been considered a major source of evolution, contributing to the origin of new gene functions and expression patterns [54]. Therefore, we further identified 27 paralogous pairs from segmental duplication events in apple, including 7 MdNF-YAs, 14 MdNF-YBs, and 6 MdNF-YCs (Figure 4 and Table S2). However, only 11 or 12 paralogous segmental duplication events had happened in S. bicolor L. or maize, while a total of 42 and 50 NF-Y proteins were respectively identified in those species [53,55]. Moreover, the genome size of plant NF-Y members studied greatly varied, from 265 Mb in peach to 742 Mb in apple, and only five pairs of paralogous events were found in peach [16]. These results all suggested that a potential correlation between MdNF-Y gene duplications and genome expansion existed during the species evolution.
In addition, the majority NF-Y genes from M. domestica were located in syntenic regions of other eight species genomes (P. betulifolia, P. persica, F. vesca, V. vinifera, A. thaliana, B. rapa, O. sativa and Z. Mays) (Figure 5 and Table S3). From the evolutionary data on NF-Y genes, we found that the number of ortholog pairs of apple and other species were related to their evolutionary relationship. It is widely known that M. domestica, F. vesca, P. betulifolia, and P. persica all belong to the Rosaceae and have a closer relationship than other species selected. As expected, they have significantly more ortholog pairs with apples than other species. Therefore, we concluded that the conservation of gene duplication during species evolution also supports the great differentiation of genome evolution. Meanwhile, the ortholog pairs from other plants can provide references to determine the function and mechanisms of apple NF-Y transcription factors.
Gene structural analysis of MdNF-Y subfamilies showed that there were many similarities in each subfamily with the corresponding subfamily in other species. For example, the genes in the MdNF-YA subunit, were interrupted by at least four introns. While, many MdNF-YB genes lacked introns. Intriguingly, MdNF-YB17 and CsNF-YB16, as the longest apple or tea NF-YB gene, both have five introns [56]. The result was consistent with the gene structure of PpNF-Y and CsNF-Y families [14,16].

3.2. Differentially Expression Pattern of MdNF-Y Genes in Apple Tissues

To date, NF-Y genes have been found to play critical roles in regulating flower and fruit development, as well as various other physiological processes; however, their roles in apple have remained unclear. Therefore, in this study, we predicted the functions of the apple NF-Y genes based on their other species ortholog pairs in syntenic regions of the two genomes (Figure 5).
In general, the NF-Y genes exhibited distinct spatiotemporal expression patterns in apple tissues and organs. However, NF-Y genes with very high sequence and exon-intron structure similarity have similar expression patterns and gene function in different apple tissues. For example, MdNF-YB1 and MdNF-YB11, located in the duplicated genomic regions, were all highly expressed in the sarcocarp and receptacle (Figure 4 and Figure 8). Apple fruits are considered ’false fruits’, since the sarcocarp has developed from the receptacle. Therefore, we hypothesized that these three homologous genes may play a synergistic role in the development of the receptacle into the fruit. In addition, the interactions between MdNF-YB17 and MdNF-YC5 or MdNF-YC8 proteins were identified by the Y2H experiments (Figure 7). Interestingly, MdNF-YC5 and MdNF-YC8 were also highly expressed in the sarcocarp or peel (Figure 8), indicating that they may play an important role in apple reproductive growth. Actually, many studies of NF-Ys in other plants have supported this hypothesis. In Arabidopsis, NF-YB and NF-YC subunits could interact with CONSTANS (CO) to form complexes, and further affected FLOWERING LOCUS T (FT) expression to induce the floral transition [2,23,57]. Analogously, ZmNF-YA3, in complex with CO and flowering promoting factor1 (FPF1) could bind to the FT-like12 promoter to promote early flowering in maize [29]. In addition, overexpression of the TaNF-YB4 gene significantly improved transgenic wheat grain yield [58].
LEAFY COTYLEDON1 (LEC1), also known as NF-YB9 is a key regulator that controls the complex process of seed development in Arabidopsis [44]. AtNF-YB9 (LEC1) and AtNF-YB6 (LEC1-Like) were both expressed in seeds [59]. Previous research has suggested that the gene expression patterns are related to the complex process of seed development which is highly coordinated both temporally and spatially in cellular processes [60]. In this study, we found that MdNF-YB13 gene promoter contained an RY-element, which is involved in plant seed-specific regulation (Figure 6). Interestingly, MdNF-YB13 and MdNF-YB10, located in the duplicated genomic regions, are ortholog genes of AtNF-YB6 (LEC1-Like) and AtNF-YB9 (LEC1) (Figureas 4 and 5 and Table S3), and also show high expression levels in seeds compared with other tissues (Figure 8). Recent studies have demonstrated that LEC1 combinates with other TFs, such as ABA-RESPONSIVE ELEMENT BINDING PROTEIN3 (AREB3), bZIP67, and ABI3 to regulate the diverse stages of seed development [61]. Moreover, previous studies have shown that LEC1 acts as a regulon to regulate hormone synthesis genes and as an integrator for light- and hormone signals to play an important role in the development of plant embryo [62]. The result of cis-acting element analysis also showed that the promoter region of MdNF-YB13 contained those important plant endogenous hormone and light elements during seed development, such as ABRE (ABA), TGA-element (auxin), and G-box (light). Taken together, our study provides a direction for future research of the analysis of the specific regulatory mechanism of MdNF-YB10 and MdNF-YB13 during apple seed development.

3.3. Function of MdNF-Y Genes in Abiotic Stress

It has been reported that TFs, including MYB, WRKY, and NAC participated in the fight against abiotic stress to help plants resist or optimize the changes in the environment [63,64,65,66]. To date, numerous studies have demonstrated that NF-Y genes also played an important role in abiotic stress response. For instance, overexpression of AtNF-YB2 and AtNF-YB3 in Arabidopsis specifically conferred tolerance to drought and heat stress, respectively [25]. Allogenic overexpression of the CdtNF-YC1 transcription factor from bermudagrass enhanced the transgenic rice tolerance under the drought and salt treatment, through ABA-dependent and ABA-independent pathways [28]. Another recent study has shown that PdNF-YB21 positively regulates the tolerance to drought stress by ABA-mediated IAA transport in Populus [30].
In addition, it is widely known that root is the major organ to respond to osmotic stress caused by drought and salt [43]. In this study, MdNF-YA6, MdNF-YA9, MdNF-YB8, MdNF-YB9, MdNF-YB18, MdNF-YB19, MdNF-YC7, and MdNF-YC8 are mainly or specifically expressed in roots (Figure 8). Further, we found the promoter region of MdNF-YB8, MdNF-YB9, MdNF-YB18, and MdNF-YC7 contains drought cis-acting element MBS (CAACTG), and MdNF-YA6 contains defence cis-element TC-rich repeats (ATTCTCTAAC/GTTTTCTTAC) (Figure 6). Meanwhile, the RT-PCR results also showed that these genes were highly expressed under drought and salinity stress (Figure 9). Therefore, we hypothesized that they may response to abiotic stress by regulating root growth. Existing research has also verified this hypothesis. For example, overexpression of the transcription factor AtNF-YB3 increased the length of primary root [67], and improved drought and heat tolerance in A. thaliana [25]. Overexpression of PdNF-YB21 in poplar promoted root growth with highly lignified and enlarged xylem vessels, resulting in increased drought resistance [30]. TaNF-YB4 plays an important role in root development and in nitrogen and phosphorus uptake in wheat [58]. Higher levels CdtNF-YC1 was detected in roots of bermudagrass and overexpression of CdtNF-YC1 elevated tolerance to drought and salt stress in transgenic rice [28].
In summary, the MdNF-Y genes displayed different degrees of responses to abiotic stress. The specific mechanisms were studied in model plants. Although the detailed molecular mechanisms of the responses of MdNF-Y genes to abiotic stress remain unclear, this study points to a number of genes that deserve further exploration in future studies.

4. Materials and Methods

4.1. Plant Materials and Treatments

The apple Xinjiang No. 1 tissue culture seedlings used for stress treatment was cultured in a specific-medium with Murashigeand Skoog (MS) medium, 0.8% agar, 0.5 mmol L−1 indole-3-butytric acid (IBA), and 0.7 mmol L−1 6-benzylaminopurine (6-BA). Apple tissue culture seedlings were first grown in specific-medium formulation under 16 h light/8 h dark at 25 °C environmental conditions for thirty days. Thirty days later, 10% polyethylene glycol (PEG) 6000, 100 mmol L−1 NaCl, and 100 mmol L−1 Abscisic acid (ABA) were added to the specific medium, respectively, to induce the stress response of apple tissue culture seedlings and cultivated under the above environmental conditions. Besides, apple tissue culture seedlings were cultured in a specific medium, and the environmental temperature was adjusted to 4 °C and 40 °C, respectively. After the above five stress treatments, samples of apple tissue culture seedlings were collected at 0, 6, 12, and 24 h respectively, frozen rapidly with liquid nitrogen, and stored in −80 °C refrigerator for succedent experiments [68]. For all treatments, three biological replicates were collected.
The different apple tissues, including root, stem, leaves, peel, receptacle, young fruits (20 days after flowering), sarcocarp, and seed was obtained from seven-year-old ‘Xinjiang No.1’ apple trees at the experimental station (longitude 120°39′ E, latitude 36°27′ N) of Qing Dao Agricultural University in 2018. These tissues materials were frozen rapidly with liquid nitrogen and stored in −80 °C refrigerator for subsequent experiments. For all sample, three biological replicates were collected.

4.2. Identification of Apple MdNF-Y Genes

The amino acid sequences of Arabidopsis 30 AtNF-Y family were obtained from the TAIR database (https://www.arabidopsis.org/). After that, The NF-Y domains were searched from BlastP (https://www.rosaceae.org/blast/protein/protein) in the NCBI database in the apple genome in the NCBI database. The Hidden Markov Models (HMM) profile of the conserved domains of NF-YA (PF02045) and NF-YB (PF00808) were downloaded from the Pfam database (http://pfam.xfam.org/). Then, they were used to search for protein sequences in the apple genome. Apple candidate NF-Y gene family members obtained by the above two methods. To ensure accuracy, SMART (http://smart.embl-heidelberg.de/) and CCD programs were used to verify the existence of conserved Pfam and complete NF-Y domains. The amino acid number, PI, and molecular weight of the identified NF-Y protein sequence was obtained from the ExPASy website (http://web.expasy.org/). The location of NF-Y gene on chromosome was determined by Gene Structure Display Server (GSDS, http://gsds.cbi.pku.edu.cn) [69].

4.3. Alignments, Synteny Analysis of MdNF-Ys

Multiple sequence alignments were performed on 43 MdNF-Y protein sequences by using DNAMAN 9 with its default. MCScanx (http://chibba.pgml.uga.edu/mcscan2/) software was used to search the homologous genes of MdNF-Y and their collinearity was obtained [70]. Circos (http://circos.ca/) software was used to analyze the collinearity of the MdNF-Y gene family [71].

4.4. Phylogenetic, Conversed Motifs, and Gene Structure Analysis of MdNF-Ys

The CDS and DNA sequences of MdNF-Y were obtained from the Apple Genome Browser. Using MEGA 7 software, the Neighbor-Joining (NJ) method was used to construct the phylogenetic tree of apple and Arabidopsis NF-Y protein sequences. Bootstrap was set to 1000 replicates. The protein sequence of the candidate MdNF-Y gene was analyzed using MEME (http://meme-suite.org/) software [72]. Set to look for 20 motifs. Phylogenetic tree, conserved motifs, and gene structure of apple MdNF-Y gene family was visualized on the TBtools toolkit [73].

4.5. Prediction of Cis-Acting Elements in Promoters of MdNF-Ys

The Plant CARE database was employed to predict the potential cis-acting elements in the 1500 bp promoters upstream region of the apple MdNF-Ys gene family [74], and visualized on the TBtools toolkit [73]. Details information for the promoters used were listed in Table S4.

4.6. Yeast Two Hybrid Assays

The CDS of MdNF-YB, MdNF-YA, or MdNF-YC genes were cloned into prey vector pGADT7 with the activating domain (AD). The CDS of MdNF-YA or MdNF-YC genes were cloned into bait vector pGBKT7 with the DNA binding domain (BD). The primers were listed in Table S6. Two constructs were co-transformed into the Y2H Gold yeast strain according to the protocol (www.weidibio.com). Then, these yeast strains were cultured on SD(-Leu/-Trp) growth medium (Clontech) and selected on SD(-Leu/-Trp/-His) and SD(-Leu/-Trp/-His/-Ade) screening medium (Clontech).

4.7. Quantitative Real-Time RT-PCR Analysis

RNA from various stress treatments (ABA, drought, salinity, heat, and cold) or different apple tissues (root, stem, leaves, peel, receptacle, young fruits (20 days after flowering), sarcocarp, and seed) was extracted with the EASYspin Plant RNA Rapid Extraction Kit (YPHBIO, Beijing, China), and the extracted RNA concentration was determined by the instrument Nano Drop 2000 (Gene Company Limited, Hong Kong, China). The cDNA was obtained by using the reverse transcription kit (Takara, Dalian, China) according to the manufacturer’s methods. Quantitative real-time PCR (qRT-PCR) was performed with ChamQ SYBR Color qPCR Master Mix (Without ROX) (Vazyme, Nanjing, China) in Roche machine (Roche, Shanghai, China). The PCR program was as followed: 95 °C for 2 min, 40 cycles of 95 °C for 30 s, 56 °C for 30 s and 72 °C for 30 s with a final dissociation stage. The internal reference gene, MdActin (nucleotide, CN938023), was used to normalize the expression levels of the tested genes. The primers for qRT-PCR experiment were designed on NCBI-BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi) and listed in Table S6. MdActin was used as an internal control. The relative expression of MdNF-Y genes was calculated with 2−ΔΔCt methods [75] and figured by GraphPad.Prism.5.0 software. Three biological replicates with three technical replicates each were measured. Transcripts profiles from different apple tissues were visualized heat-mapped using the TBtools toolbox [73].

5. Conclusions

In this work, 43 MdNF-Y genes (11 MdNF-YAs, 22 MdNF-YBs, and 10 MdNF-YCs) were identified and their evolutionary, structure, biological function, and expression pattern were analyzed. Based on prediction and experimental data, MdNF-Ys might play an important role in apple development and response to five abiotic stress (ABA, drought, heat, cold, and salinity). Our findings will contribute a foundation for further study of the functional and regulatory mechanisms controlled by the NF-Y gene family in apple.

Supplementary Materials

Supplementary Materials can be found at https://www.mdpi.com/2223-7747/10/1/16/s1. Table S1: MdNF-Y protein sequences and CDS length. Table S2: The paralogous gene pairs of 43 MdNF-Ys. Table S3: The ortholog pairs of apples with eight species. Table S4: The 1500 bp promoter sequences of the 43 MdNF-Ys. Table S5: cis-element analysis of the 43 MdNF-Ys in the 1500 bp promoter sequences. Table S6: The primers used in this study.

Author Contributions

Y.Z. and H.H., designed the study; Y.Q., Y.W., and H.H., performed data analysis; Y.Z. and J.Z., provided guidance on the whole study. Y.Q., Y.W., H.H., and Y.Z., wrote and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (31801841); China Agriculture Research System Foundation (Grant No. CARS-27); Natural Science Foundation of Shandong (ZR2019MC003); Qingdao Scientific Research Foundation (Grant No. 19-6-1-60-nsh); Breeding Plan of Shandong Provincial Qingchuang Research Team (2019).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in supplementary material S1 to S6.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mantovani, R. The molecular biology of the CCAAT-binding factor NF-Y. Gene 1999, 239, 15–27. [Google Scholar] [CrossRef]
  2. Siefers, N.; Dang, K.K.; Kumimoto, R.W.; Bynum, W.E.t.; Tayrose, G.; Holt, B.F., 3rd. Tissue-specific expression patterns of Arabidopsis NF-Y transcription factors suggest potential for extensive combinatorial complexity. Plant Physiol. 2009, 149, 625–641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Petroni, K.; Kumimoto, R.W.; Gnesutta, N.; Calvenzani, V.; Fornari, M.; Tonelli, C.; Holt, B.F., 3rd; Mantovani, R. The promiscuous life of plant NUCLEAR FACTOR Y transcription factors. Plant Cell 2012, 24, 4777–4792. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Frontini, M.; Imbriano, C.; Manni, I.; Mantovani, R. Cell cycle regulation of NF-YC nuclear localization. Cell Cycle 2004, 3, 217–222. [Google Scholar] [CrossRef] [Green Version]
  5. Kahle, J.; Baake, M.; Doenecke, D.; Albig, W. Subunits of the heterotrimeric transcription factor NF-Y are imported into the nucleus by distinct pathways involving importin beta and importin 13. Mol. Cell Biol. 2005, 25, 5339–5354. [Google Scholar] [CrossRef] [Green Version]
  6. Liu, J.X.; Howell, S.H. bZIP28 and NF-Y transcription factors are activated by ER stress and assemble into a transcriptional complex to regulate stress response genes in Arabidopsis. Plant Cell 2010, 22, 782–796. [Google Scholar] [CrossRef] [Green Version]
  7. Kumimoto, R.W.; Siriwardana, C.L.; Gayler, K.K.; Risinger, J.R.; Siefers, N.; Holt, B.F., 3rd. Nuclear Factor y transcription factors have both opposing and additive roles in ABA-mediated seed germination. PLoS ONE 2013, 8, e59481. [Google Scholar] [CrossRef] [Green Version]
  8. Bi, C.; Ma, Y.; Wang, X.-F.; Zhang, D.-P. Overexpression of the transcription factor NF-YC9 confers abscisic acid hypersensitivity in Arabidopsis. Plant Mol. Biol. 2017, 95, 425–439. [Google Scholar] [CrossRef] [Green Version]
  9. Qi, M.; Zheng, W.; Zhao, X.; Hohenstein, J.D.; Kandel, Y.; O’Conner, S.; Wang, Y.; Du, C.; Nettleton, D.; MacIntosh, G.C.; et al. QQS orphan gene and its interactor NF-YC4 reduce susceptibility to pathogens and pests. Plant Biotechnol. J. 2019, 17, 252–263. [Google Scholar] [CrossRef] [Green Version]
  10. Li, L.; Zheng, W.; Zhu, Y.; Ye, H.; Tang, B.; Arendsee, Z.W.; Jones, D.; Li, R.; Ortiz, D.; Zhao, X.; et al. QQS orphan gene regulates carbon and nitrogen partitioning across species via NF-YC interactions. Proc. Natl. Acad. Sci. USA 2015, 112, 14734–14739. [Google Scholar] [CrossRef] [Green Version]
  11. Quan, S.; Niu, J.; Zhou, L.; Xu, H.; Ma, L.; Qin, Y. Identification and characterization of NF-Y gene family in walnut (Juglans regia L.). BMC Plant Biol. 2018, 18, 255. [Google Scholar] [CrossRef] [Green Version]
  12. Pereira, S.L.S.; Martins, C.P.S.; Sousa, A.O.; Camillo, L.R.; Araujo, C.P.; Alcantara, G.M.; Camargo, D.S.; Cidade, L.C.; de Almeida, A.F.; Costa, M.G.C. 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, e0199187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Chu, H.D.; Nguyen, K.H.; Watanabe, Y.; Le, D.T.; Pham, T.L.T.; Mochida, K.; Tran, L.P. 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. Int. J. Mol. Sci. 2018, 19, 3290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Wang, Y.; Xu, W.; Chen, Z.; Han, B.; Haque, M.E.; Liu, A. Gene structure, expression pattern and interaction of Nuclear Factor-Y family in castor bean (Ricinus communis). Planta 2018, 247, 559–572. [Google Scholar] [CrossRef] [PubMed]
  15. Ren, C.; Zhang, Z.; Wang, Y.; Li, S.; Liang, Z. Genome-wide identification and characterization of the NF-Y gene family in grape (vitis vinifera L.). BMC Genom. 2016, 17, 605. [Google Scholar] [CrossRef] [Green Version]
  16. Li, M.; Li, G.; Liu, W.; Dong, X.; Zhang, A. Genome-wide analysis of the NF-Y gene family in peach (Prunus persica L.). BMC Genom. 2019, 20, 612. [Google Scholar] [CrossRef] [Green Version]
  17. Liu, X.; Hu, P.; Huang, M.; Tang, Y.; Li, Y.; Li, L.; Hou, X. The NF-YC-RGL2 module integrates GA and ABA signalling to regulate seed germination in Arabidopsis. Nat. Commun. 2016, 7, 12768. [Google Scholar] [CrossRef] [Green Version]
  18. Li, T.; Zhang, H.; Liu, Z.; Deng, H.; Sharma, S.; Wei, X.; Wang, L.; Niu, B.; Chen, C. A group of nuclear factor Y transcription factors are sub-functionalized during endosperm development in monocots. J. Exp. Bot. 2018, 69, 2495–2510. [Google Scholar] [CrossRef]
  19. Gnesutta, N.; Kumimoto, R.W.; Swain, S.; Chiara, M.; Siriwardana, C.; Horner, D.S.; Holt, B.F., 3rd; Mantovani, R. Constans Imparts DNA Sequence Specificity to the Histone Fold NF-YB/NF-YC Dimer. Plant Cell 2017, 29, 1516–1532. [Google Scholar] [CrossRef] [Green Version]
  20. Tokutsu, R.; Fujimura-Kamada, K.; Matsuo, T.; Yamasaki, T.; Minagawa, J. The constans flowering complex controls the protective response of photosynthesis in the green alga Chlamydomonas. Nat. Commun. 2019, 10, 4099. [Google Scholar] [CrossRef]
  21. Xuanyuan, G.; Lu, C.; Zhang, R.; Jiang, J. Overexpression of StNF-YB3.1 reduces photosynthetic capacity and tuber production, and promotes ABA-mediated stomatal closure in potato (Solanum tuberosum L.). Plant Sci. 2017, 261, 50–59. [Google Scholar] [CrossRef] [PubMed]
  22. Zhao, M.; Ding, H.; Zhu, J.K.; Zhang, F.; Li, W.X. Involvement of miR169 in the nitrogen-starvation responses in Arabidopsis. New Phytol. 2011, 190, 906–915. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Hou, X.; Zhou, J.; Liu, C.; Liu, L.; Shen, L.; Yu, H. Nuclear factor Y-mediated H3K27me3 demethylation of the SOC1 locus orchestrates flowering responses of Arabidopsis. Nat. Commun. 2014, 5, 4601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lee, D.-K.; Kim, H.I.; Jang, G.; Chung, P.J.; Jeong, J.S.; Kim, Y.S.; Bang, S.W.; Jung, H.; Choi, Y.D.; Kim, J.-K. The NF-YA transcription factor OsNF-YA7 confers drought stress tolerance of rice in an abscisic acid independent manner. Plant Sci. Int. J. Exp. Plant Biol. 2015, 241, 199–210. [Google Scholar] [CrossRef]
  25. Sato, H.; Suzuki, T.; Takahashi, F.; Shinozaki, K.; Yamaguchi-Shinozaki, K. NF-YB2 and NF-YB3 Have Functionally Diverged and Differentially Induce Drought and Heat Stress-Specific Genes. Plant Physiol. 2019, 180, 1677–1690. [Google Scholar] [CrossRef] [Green Version]
  26. Li, W.-X.; Oono, Y.; Zhu, J.; He, X.-J.; Wu, J.-M.; Iida, K.; Lu, X.-Y.; Cui, X.; Jin, H.; Zhu, J.-K. The Arabidopsis NFYA5 transcription factor is regulated transcriptionally and posttranscriptionally to promote drought resistance. Plant Cell 2008, 20, 2238–2251. [Google Scholar] [CrossRef] [Green Version]
  27. Nelson, D.E.; Repetti, P.P.; Adams, T.R.; Creelman, R.A.; Wu, J.; Warner, D.C.; Anstrom, D.C.; Bensen, R.J.; Castiglioni, P.P.; Donnarummo, M.G.; et al. Plant nuclear factor Y (NF-Y) B subunits confer drought tolerance and lead to improved corn yields on water-limited acres. Proc. Natl. Acad. Sci. USA 2007, 104, 16450–16455. [Google Scholar] [CrossRef] [Green Version]
  28. 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 Biotechnol. J. 2015, 13, 482–491. [Google Scholar] [CrossRef]
  29. Su, H.; Cao, Y.; Ku, L.; Yao, W.; Cao, Y.; Ren, Z.; Dou, D.; Wang, H.; Ren, Z.; Liu, H.; et al. Dual functions of ZmNF-YA3 in photoperiod-dependent flowering and abiotic stress responses in maize. J. Exp. Bot. 2018, 69, 5177–5189. [Google Scholar] [CrossRef] [Green Version]
  30. Zhou, Y.; Zhang, Y.; Wang, X.; Han, X.; An, Y.; Lin, S.; Shen, C.; Wen, J.; Liu, C.; Yin, W.; 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. [Google Scholar] [CrossRef]
  31. Pawełkowicz, M.E.; Skarzyńska, A.; Sroka, M.; Szwacka, M.; Pniewski, T.; Pląder, W. Effect of Transgenesis on mRNA and miRNA Profiles in Cucumber Fruits Expressing Thaumatin II. Genes 2020, 11, 334. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Yu, Y.; Ni, Z.; Wang, Y.; Wan, H.; Hu, Z.; Jiang, Q.; Sun, X.; Zhang, H. Overexpression of soybean miR169c confers increased drought stress sensitivity in transgenic Arabidopsis thaliana. Plant Sci. Int. J. Exp. Plant Biol. 2019, 285, 68–78. [Google Scholar] [CrossRef] [PubMed]
  33. Luan, M.; Xu, M.; Lu, Y.; Zhang, L.; Fan, Y.; Wang, L. Expression of zma-miR169 miRNAs and their target ZmNF-YA genes in response to abiotic stress in maize leaves. Gene 2015, 555, 178–185. [Google Scholar] [CrossRef]
  34. Luan, M.; Xu, M.; Lu, Y.; Zhang, Q.; Zhang, L.; Zhang, C.; Fan, Y.; Lang, Z.; Wang, L. Family-wide survey of miR169s and NF-YAs and their expression profiles response to abiotic stress in maize roots. PLoS ONE 2014, 9, e91369. [Google Scholar] [CrossRef] [Green Version]
  35. Velasco, R.; Zharkikh, A.; Affourtit, J.; Dhingra, A.; Cestaro, A.; Kalyanaraman, A.; Fontana, P.; Bhatnagar, S.K.; Troggio, M.; Pruss, D.; et al. The genome of the domesticated apple (Malus × domestica Borkh.). Nat. Genet. 2010, 42, 833–839. [Google Scholar] [CrossRef] [PubMed]
  36. Laloum, T.; De Mita, S.; Gamas, P.; Baudin, M.; Niebel, A. CCAAT-box binding transcription factors in plants: Y so many? Trends Plant Sci. 2013, 18, 157–166. [Google Scholar] [CrossRef]
  37. Dolfini, D.; Gatta, R.; Mantovani, R. NF-Y and the transcriptional activation of CCAAT promoters. Crit. Rev. Biochem. Mol. Biol. 2012, 47, 29–49. [Google Scholar] [CrossRef]
  38. Matuoka, K.; Chen, K.Y. Transcriptional regulation of cellular ageing by the CCAAT box-binding factor CBF/NF-Y. Ageing Res. Rev. 2002, 1, 639–651. [Google Scholar] [CrossRef]
  39. Yamaguchi-Shinozaki, K.; Mundy, J.; Chua, N.H. Four tightly linked rab genes are differentially expressed in rice. Plant Mol. Biol. 1990, 14, 29–39. [Google Scholar] [CrossRef]
  40. Anderson, S.L.; Teakle, G.R.; Martino-Catt, S.J.; Kay, S.A. Circadian clock- and phytochrome-regulated transcription is conferred by a 78 bp cis-acting domain of the Arabidopsis CAB2 promoter. Plant J. 1994, 6, 457–470. [Google Scholar] [CrossRef]
  41. Fujiwara, T.; Beachy, R.N. Tissue-specific and temporal regulation of a beta-conglycinin gene: Roles of the RY repeat and other cis-acting elements. Plant Mol. Biol. 1994, 24, 261–272. [Google Scholar] [CrossRef] [PubMed]
  42. 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 Mol. Biol. 2013, 82, 113–129. [Google Scholar] [CrossRef] [PubMed]
  43. Kreszies, T.; Eggels, S.; Kreszies, V.; Osthoff, A.; Shellakkutti, N.; Baldauf, J.A.; Zeisler-Diehl, V.V.; Hochholdinger, F.; Ranathunge, K.; Schreiber, L. Seminal roots of wild and cultivated barley differentially respond to osmotic stress in gene expression, suberization, and hydraulic conductivity. Plant Cell Environ. 2020, 43, 344–357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Boulard, C.; Thevenin, J.; Tranquet, O.; Laporte, V.; Lepiniec, L.; Dubreucq, B. LEC1 (NF-YB9) directly interacts with LEC2 to control gene expression in seed. Biochim. Biophys. Acta Gene Regul. Mech. 2018, 1861, 443–450. [Google Scholar] [CrossRef] [PubMed]
  45. Wright, K.L.; Moore, T.L.; Vilen, B.J.; Brown, A.M.; Ting, J.P. Major histocompatibility complex class II-associated invariant chain gene expression is up-regulated by cooperative interactions of Sp1 and NF-Y. J. Biol. Chem. 1995, 270, 20978–20986. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Benatti, P.; Basile, V.; Merico, D.; Fantoni, L.I.; Tagliafico, E.; Imbriano, C. A balance between NF-Y and p53 governs the pro- and anti-apoptotic transcriptional response. Nucleic Acids Res. 2008, 36, 1415–1428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Swain, S.; Myers, Z.A.; Siriwardana, C.L.; Holt, B.F., 3rd. The multifaceted roles of NUCLEAR FACTOR-Y in Arabidopsis thaliana development and stress responses. Biochim. Biophys. Acta Gene Regul. Mech. 2017, 1860, 636–644. [Google Scholar] [CrossRef]
  48. Albani, D.; Robert, L.S. Cloning and characterization of a Brassica napus gene encoding a homologue of the B subunit of a heteromeric CCAAT-binding factor. Gene 1995, 167, 209–213. [Google Scholar] [CrossRef]
  49. Sorin, C.; Declerck, M.; Christ, A.; Blein, T.; Ma, L.; Lelandais-Briere, C.; Njo, M.F.; Beeckman, T.; Crespi, M.; Hartmann, C. A miR169 isoform regulates specific NF-YA targets and root architecture in Arabidopsis. New Phytol. 2014, 202, 1197–1211. [Google Scholar] [CrossRef]
  50. Stephenson, T.J.; McIntyre, C.L.; Collet, C.; Xue, G.P. Genome-wide identification and expression analysis of the NF-Y family of transcription factors in Triticum aestivum. Plant Mol. Biol. 2007, 65, 77–92. [Google Scholar] [CrossRef] [Green Version]
  51. Thirumurugan, T.; Ito, Y.; Kubo, T.; Serizawa, A.; Kurata, N. Identification, characterization and interaction of HAP family genes in rice. Mol. Genet. Genom. 2008, 279, 279–289. [Google Scholar] [CrossRef]
  52. Li, S.; Li, K.; Ju, Z.; Cao, D.; Fu, D.; Zhu, H.; Zhu, B.; Luo, Y. Genome-wide analysis of tomato NF-Y factors and their role in fruit ripening. BMC Genom. 2016, 17, 36. [Google Scholar] [CrossRef] [Green Version]
  53. Zhang, Z.; Li, X.; Zhang, C.; Zou, H.; Wu, Z. Isolation, structural analysis, and expression characteristics of the maize nuclear factor Y gene families. Biochem. Biophys. Res. Commun. 2016, 478, 752–758. [Google Scholar] [CrossRef] [PubMed]
  54. Lynch, M.; Conery, J.S. The evolutionary fate and consequences of duplicate genes. Science 2000, 290, 1151–1155. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Maheshwari, P.; Kummari, D.; Palakolanu, S.R.; Nagasai Tejaswi, U.; Nagaraju, M.; Rajasheker, G.; Jawahar, G.; Jalaja, N.; Rathnagiri, P.; Kavi Kishor, P.B. Genome-wide identification and expression profile analysis of nuclear factor Y family genes in Sorghum bicolor L. (Moench). PLoS ONE 2019, 14, e0222203. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, P.; Zheng, Y.; Guo, Y.; Chen, X.; Sun, Y.; Yang, J.; Ye, N. Identification, expression, and putative target gene analysis of nuclear factor-Y (NF-Y) transcription factors in tea plant (Camellia sinensis). Planta 2019, 250, 1671–1686. [Google Scholar] [CrossRef]
  57. Ben-Naim, O.; Eshed, R.; Parnis, A.; Teper-Bamnolker, P.; Shalit, A.; Coupland, G.; Samach, A.; Lifschitz, E. The CCAAT binding factor can mediate interactions between CONSTANS-like proteins and DNA. Plant J. Cell Mol. Biol. 2006, 46, 462–476. [Google Scholar] [CrossRef] [Green Version]
  58. Yadav, D.; Shavrukov, Y.; Bazanova, N.; Chirkova, L.; Borisjuk, N.; Kovalchuk, N.; Ismagul, A.; Parent, B.; Langridge, P.; Hrmova, M.; et al. Constitutive overexpression of the TaNF-YB4 gene in transgenic wheat significantly improves grain yield. J Exp. Bot. 2015, 66, 6635–6650. [Google Scholar] [CrossRef] [Green Version]
  59. Warpeha, K.M.; Upadhyay, S.; Yeh, J.; Adamiak, J.; Hawkins, S.I.; Lapik, Y.R.; Anderson, M.B.; Kaufman, L.S. The GCR1, GPA1, PRN1, NF-Y signal chain mediates both blue light and abscisic acid responses in Arabidopsis. Plant Physiol. 2007, 143, 1590–1600. [Google Scholar] [CrossRef] [Green Version]
  60. Jo, L.; Pelletier, J.M.; Harada, J.J. Central role of the leafy cotyledon1 transcription factor in seed development. J. Integr. Plant Biol. 2019, 61, 564–580. [Google Scholar] [CrossRef] [Green Version]
  61. Jo, L.; Pelletier, J.M.; Hsu, S.W.; Baden, R.; Goldberg, R.B.; Harada, J.J. Combinatorial interactions of the LEC1 transcription factor specify diverse developmental programs during soybean seed development. Proc. Natl. Acad. Sci. USA 2020, 117, 1223–1232. [Google Scholar] [CrossRef] [PubMed]
  62. Junker, A.; Mönke, G.; Rutten, T.; Keilwagen, J.; Seifert, M.; Thi, T.M.; Renou, J.P.; Balzergue, S.; Viehöver, P.; Hähnel, U.; et al. Elongation-related functions of LEAFY COTYLEDON1 during the development of Arabidopsis thaliana. Plant J. 2012, 71, 427–442. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Ambawat, S.; Sharma, P.; Yadav, N.R.; Yadav, R.C. MYB transcription factor genes as regulators for plant responses: An overview. Physiol. Mol. Biol. Plants Int. J. Funct. Plant Biol. 2013, 19, 307–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Nuruzzaman, M.; Sharoni, A.M.; Kikuchi, S. Roles of NAC transcription factors in the regulation of biotic and abiotic stress responses in plants. Front. Microbiol. 2013, 4, 248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Jiang, J.; Ma, S.; Ye, N.; Jiang, M.; Cao, J.; Zhang, J. WRKY transcription factors in plant responses to stresses. J. Integr. Plant Biol. 2017, 59, 86–101. [Google Scholar] [CrossRef]
  66. Myers, Z.A.; Holt, B.F.I. NUCLEAR FACTOR-Y: Still complex after all these years? Curr. Opin. Plant Biol. 2018, 45, 96–102. [Google Scholar] [CrossRef]
  67. Ballif, J.; Endo, S.; Kotani, M.; MacAdam, J.; Wu, Y. Over-expression of HAP3b enhances primary root elongation in Arabidopsis. Plant Physiol. Biochem. 2011, 49, 579–583. [Google Scholar] [CrossRef]
  68. Hou, H.; Lv, L.; Huo, H.; Dai, H.; Zhang, Y. Genome-Wide Identification of the ABA Receptors Genes and Their Response to Abiotic Stress in Apple. Plants 2020, 9, 1028. [Google Scholar] [CrossRef]
  69. Hu, B.; Jin, J.; Guo, A.Y.; Zhang, H.; Luo, J.; Gao, G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 2015, 31, 1296–1297. [Google Scholar] [CrossRef] [Green Version]
  70. Wang, Y.; Tang, H.; Debarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [Green Version]
  71. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. Meme suite: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef] [PubMed]
  73. Chen, C.; Xia, R.; Chen, H.; He, Y. TBtools, a Toolkit for Biologists integrating various HTS-data handling tools with a user-friendly interface. bioRxiv 2018, 289660. [Google Scholar] [CrossRef]
  74. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van de Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef] [PubMed]
  75. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Plants 10 00016 g001 550
Figure 1. Multiple sequences alignment of the conserved domains between MdNF-Y (M. domestica), NF-Y Arabidopsis (A. thaliana) and NF-Y Grape (V. vinifera). The level of these three species amino acid homology = 100%, ≥75% and ≥50% are colored by blue, cyan, and pink boxes, respectively. (A) MdNF-YA subfamily; (B) MdNF-YB subfamily; (C) MdNF-YC subfamily. Multiple sequences alignment was constructed by DNAMAN software.
Figure 1. Multiple sequences alignment of the conserved domains between MdNF-Y (M. domestica), NF-Y Arabidopsis (A. thaliana) and NF-Y Grape (V. vinifera). The level of these three species amino acid homology = 100%, ≥75% and ≥50% are colored by blue, cyan, and pink boxes, respectively. (A) MdNF-YA subfamily; (B) MdNF-YB subfamily; (C) MdNF-YC subfamily. Multiple sequences alignment was constructed by DNAMAN software.
Plants 10 00016 g001
Plants 10 00016 g002 550
Figure 2. Conserved motifs (A) and gene structure (B) analysis in each MdNF-Y subfamily. Different colors of figure A represent the type of motifs MdNF-Y genes have. The green and yellow color box of figure B respectively indicates the coding sequences (CDS) and untranslated regions (UTR). Intron are shown with black lines. Scale bars located on the bottom side of figure representing the relative position based on the kilobase scale.
Figure 2. Conserved motifs (A) and gene structure (B) analysis in each MdNF-Y subfamily. Different colors of figure A represent the type of motifs MdNF-Y genes have. The green and yellow color box of figure B respectively indicates the coding sequences (CDS) and untranslated regions (UTR). Intron are shown with black lines. Scale bars located on the bottom side of figure representing the relative position based on the kilobase scale.
Plants 10 00016 g002
Plants 10 00016 g003 550
Figure 3. Physical distribution and chromosomal location of 43 MdNF-Ys genes in apple genome. 43 MdNF-Ys genes were mapped onto 16 chromosomes. Hollow bars indicated the chromosomes. Chromosomes and MdNF-Ys genes names were shown at the top and right of the bar, respectively. Scale bars located on left side of figure representing the length of each chromosome are in megabases (Mb).
Figure 3. Physical distribution and chromosomal location of 43 MdNF-Ys genes in apple genome. 43 MdNF-Ys genes were mapped onto 16 chromosomes. Hollow bars indicated the chromosomes. Chromosomes and MdNF-Ys genes names were shown at the top and right of the bar, respectively. Scale bars located on left side of figure representing the length of each chromosome are in megabases (Mb).
Plants 10 00016 g003
Plants 10 00016 g004 550
Figure 4. Synteny analysis of MdNF-YAs (A), MdNF-YBs (B) and MdNF-YCs (C) in the apple. The chromosomal localizations were shown for apple was random colors (chr01-17). The red lines indicated the segmental duplication. Gray lines represented all synteny blocks in the apple genome.
Figure 4. Synteny analysis of MdNF-YAs (A), MdNF-YBs (B) and MdNF-YCs (C) in the apple. The chromosomal localizations were shown for apple was random colors (chr01-17). The red lines indicated the segmental duplication. Gray lines represented all synteny blocks in the apple genome.
Plants 10 00016 g004
Plants 10 00016 g005 550
Figure 5. Synteny analysis of MdNF-Ys with other eight species, including P. betulifolia (pear), P. persica (peach), F. vesca (strawberry), V. vinifera(grape), A. thaliana (Arabidopsis), O. sativa (rice), Z. Mays (maize), and B. rapa (Chinese cabbage). The red solid bars represent the chromosomes of apple, while the black solid bars represent the chromosomes of other species. The gray lines represent the collinear blocks within apple and other plant genomes, while the yellow, purple, and blue lines highlight syntenic MdNF-YA, MdNF-YB, and MdNF-YC gene pairs, respectively.
Figure 5. Synteny analysis of MdNF-Ys with other eight species, including P. betulifolia (pear), P. persica (peach), F. vesca (strawberry), V. vinifera(grape), A. thaliana (Arabidopsis), O. sativa (rice), Z. Mays (maize), and B. rapa (Chinese cabbage). The red solid bars represent the chromosomes of apple, while the black solid bars represent the chromosomes of other species. The gray lines represent the collinear blocks within apple and other plant genomes, while the yellow, purple, and blue lines highlight syntenic MdNF-YA, MdNF-YB, and MdNF-YC gene pairs, respectively.
Plants 10 00016 g005
Plants 10 00016 g006 550
Figure 6. The cis-acting regulatory elements of promoters in apple MdNF-Ys genes. The blue, green, and orange blocks represent phytohormone and abiotic stress, plant development, and light responsive cis-elements, respectively. The number of cis-acting elements was indicated by different colors and circle sizes. The size of green to yellow circle represented the number of cis-acting elements.
Figure 6. The cis-acting regulatory elements of promoters in apple MdNF-Ys genes. The blue, green, and orange blocks represent phytohormone and abiotic stress, plant development, and light responsive cis-elements, respectively. The number of cis-acting elements was indicated by different colors and circle sizes. The size of green to yellow circle represented the number of cis-acting elements.
Plants 10 00016 g006
Plants 10 00016 g007 550
Figure 7. Protein–protein interactions between MdNF-Y subunits by yeast two hybrid. MdNF-YB genes were constructed to vectors with GAL4 activation domain (AD), MdNF-YA and MdNF-YC genes with GAL4 DNA-binding domain (BD).
Figure 7. Protein–protein interactions between MdNF-Y subunits by yeast two hybrid. MdNF-YB genes were constructed to vectors with GAL4 activation domain (AD), MdNF-YA and MdNF-YC genes with GAL4 DNA-binding domain (BD).
Plants 10 00016 g007
Plants 10 00016 g008 550
Figure 8. Expression profiles of MdNF-Ys in various apple tissues including roots, stems, leaves, receptacles, sarcocarps, young fruits, and seeds. In the heat map, values were transformed to log2 (value). Green, low expression; black, medium expression; red, high expression. MdActin was used as an internal control. The results were based on three biological replicates and three technical replicates.
Figure 8. Expression profiles of MdNF-Ys in various apple tissues including roots, stems, leaves, receptacles, sarcocarps, young fruits, and seeds. In the heat map, values were transformed to log2 (value). Green, low expression; black, medium expression; red, high expression. MdActin was used as an internal control. The results were based on three biological replicates and three technical replicates.
Plants 10 00016 g008
Plants 10 00016 g009a 550Plants 10 00016 g009b 550Plants 10 00016 g009c 550
Figure 9. Expression profiles of MdNF-YAs (A), MdNF-YBs (B) and MdNF-YCs (C) in response to abiotic stress treatment, including ABA, drought, salinity, heat, and cold. The black, pink, purple, and gray bars represented the abiotic stress treatment time from 0, 6, 12, to 24 h, respectively. The black lines indicated error bars. The expression of MdNF-Ys at 0 h were set to 1 and MdActin was used as an internal control. Asterisks denote significance determined by t-test: * p < 0.05 and ** p < 0.01. The results were based on three biological replicates and three technical replicates.
Figure 9. Expression profiles of MdNF-YAs (A), MdNF-YBs (B) and MdNF-YCs (C) in response to abiotic stress treatment, including ABA, drought, salinity, heat, and cold. The black, pink, purple, and gray bars represented the abiotic stress treatment time from 0, 6, 12, to 24 h, respectively. The black lines indicated error bars. The expression of MdNF-Ys at 0 h were set to 1 and MdActin was used as an internal control. Asterisks denote significance determined by t-test: * p < 0.05 and ** p < 0.01. The results were based on three biological replicates and three technical replicates.
Plants 10 00016 g009aPlants 10 00016 g009bPlants 10 00016 g009c
Table
Table 1. Information of the MdNF-Ys genes in apple.
Table 1. Information of the MdNF-Ys genes in apple.
NameGeneChromosome Location (bp)DNA
(bp)		CDS LengthProtein (aa)pIMW(Da)
MdNF-YA1MD02G10862006760401–676392335229273089.0233,933.04
MdNF-YA2MD02G130980036538789–36544223543410023338.5236,822.80
MdNF-YA3MD03G117460023867258–23872067480910773587.0438,855.31
MdNF-YA4MD05G127330040831929–408368054876207686.427449.24
MdNF-YA5MD07G10113001023634–102910154679813268.4135,657.67
MdNF-YA6MD09G118620016163341–1616848051398972989.4332,857.83
MdNF-YA7MD10G125330034584154–3458871445606272088.1622,824.18
MdNF-YA8MD11G119240027332752–27339954720210473488.4337,923.42
MdNF-YA9MD12G10421004650199–465516949709693229.0735,109.45
MdNF-YA10MD14G10413003873166–387812949639303099.6333,939.24
MdNF-YA11MD15G121340017122823–1712580129789033009.2533,209.33
MdNF-YB1MD01G111240022658861–2266239935386452148.7022,929.78
MdNF-YB2MD02G119190017939911–179406457347352445.2527,017.25
MdNF-YB3MD02G119210018023631–18025409177810443474.8038,127.12
MdNF-YB4MD03G117960024691097–246916635665671884.8020,686.4
MdNF-YB5MD03G118340025025010–250255945845851946.7620,837.01
MdNF-YB6MD03G128370036413956–364145165605611867.0020,947.32
MdNF-YB7MD04G110470019216567–1922009035235251747.0518,780.78
MdNF-YB8MD04G120330028953974–289544264524531506.6016,863.98
MdNF-YB9MD05G136100047685520–476858463263271085.1111,910.38
MdNF-YB10MD05G136160047696633–4769897723448762916.9631,267.54
MdNF-YB11MD07G118020025950191–2595349533046482158.7123,159.02
MdNF-YB12MD09G11269009795758–97964116536542177.4422,184.36
MdNF-YB13MD10G133900041496215–414969997847442477.3527,349.57
MdNF-YB14MD11G119930028769434–287699855515521835.8920,252.83
MdNF-YB15MD11G120040028876837–288774095725731907.4320,088.35
MdNF-YB16MD11G130250041788083–417885804974981656.6918,596.63
MdNF-YB17MD12G112480019970949–1998077398245191725.8118,595.6
MdNF-YB18MD12G121710029460289–294607414524531507.8416,794.94
MdNF-YB19MD13G117020013815353–1381582647313744576.6017,505.45
MdNF-YB20MD15G11346009778391–97791017107112366.8526,721.4
MdNF-YB21MD15G133420037152321–371529476266272087.4623,010.4
MdNF-YB22MD17G111720010167264–101682409767202395.5824,810.33
MdNF-YC1MD04G120750029269117–2927140622898072686.4629,749.6
MdNF-YC2MD05G122960036298931–3630086219318072685.3629,718.64
MdNF-YC3MD06G107850019362575–193629403653661217.7813,368.54
MdNF-YC4MD06G114120028495891–2849887629857202395.1525,832.07
MdNF-YC5MD10G120880030806582–3080916625848342776.1131,174.24
MdNF-YC6MD12G122180029861745–2986413723927892626.0529,217.11
MdNF-YC7MD13G118520015766106–1576877126657562516.1928,006.56
MdNF-YC8MD14G115640025090589–2509313425457172386.052577.04
MdNF-YC9MD15G115940011916470–119171717017022335.4026,518.35
MdNF-YC10MD16G118590016034439–1603711126727382456.1927,400.74
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qu, Y.; Wang, Y.; Zhu, J.; Zhang, Y.; Hou, H. Genomic Organization, Phylogenetic Comparison, and Differential Expression of the Nuclear Factor-Y Gene Family in Apple (Malus Domestica). Plants 2021, 10, 16. https://doi.org/10.3390/plants10010016

AMA Style

Qu Y, Wang Y, Zhu J, Zhang Y, Hou H. Genomic Organization, Phylogenetic Comparison, and Differential Expression of the Nuclear Factor-Y Gene Family in Apple (Malus Domestica). Plants. 2021; 10(1):16. https://doi.org/10.3390/plants10010016

Chicago/Turabian Style

Qu, Yanjie, Yaping Wang, Jun Zhu, Yugang Zhang, and Hongmin Hou. 2021. "Genomic Organization, Phylogenetic Comparison, and Differential Expression of the Nuclear Factor-Y Gene Family in Apple (Malus Domestica)" Plants 10, no. 1: 16. https://doi.org/10.3390/plants10010016

APA Style

Qu, Y., Wang, Y., Zhu, J., Zhang, Y., & Hou, H. (2021). Genomic Organization, Phylogenetic Comparison, and Differential Expression of the Nuclear Factor-Y Gene Family in Apple (Malus Domestica). Plants, 10(1), 16. https://doi.org/10.3390/plants10010016

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop