Skip to main content

Genome-wide characterization of PEBP family genes in nine Rosaceae tree species and their expression analysis in P. mume

Abstract

Background

Phosphatidylethanolamine-binding proteins (PEBPs) constitute a common gene family found among animals, plants and microbes. Plant PEBP proteins play an important role in regulating flowering time, plant architecture as well as seed dormancy. Though PEBP family genes have been well studied in Arabidopsis and other model species, less is known about these genes in perennial trees.

Results

To understand the evolution of PEBP genes and their functional roles in flowering control, we identified 56 PEBP members belonging to three gene clades (MFT-like, FT-like, and TFL1-like) and five lineages (FT, BFT, CEN, TFL1, and MFT) across nine Rosaceae perennial species. Structural analysis revealed highly conserved gene structure and protein motifs among Rosaceae PEBP proteins. Codon usage analysis showed slightly biased codon usage across five gene lineages. With selection pressure analysis, we detected strong purifying selection constraining divergence within most lineages, while positive selection driving the divergence of FT-like and TFL1-like genes from the MFT-like gene clade. Spatial and temporal expression analyses revealed the essential role of FT in regulating floral bud breaking and blooming in P. mume. By employing a weighted gene co-expression network approach, we inferred a putative FT regulatory module required for dormancy release and blooming in P. mume.

Conclusions

We have characterized the PEBP family genes in nine Rosaceae species and examined their phylogeny, genomic syntenic relationship, duplication pattern, and expression profiles during flowering process. These results revealed the evolutionary history of PEBP genes and their functions in regulating floral bud development and blooming among Rosaceae tree species.

Background

Proper timing of flowering is a key adaptive strategy in plant species, especially temperate woody perennials [1,2,3]. The flowering time in annual or biennials is largely determined by the timing of the transition from vegetative growth to reproductive growth [4, 5]. However, in temperate tree species, flower buds initiate and develop during summer, undergo a short period of dormancy, exit dormancy after exposure to chilling temperatures and finally bloom in suitable environments [6]. Therefore, the blooming time of temperate woody perennials is mainly determined by intrinsic state of flower buds and external environment [7, 8]. Within the context of global climate change, warm winters and irregular occurrences of extreme weather have disrupted the timing of spring phenological events in tree species, increased the risk of frost damage, and caused abnormal fertility and poor fruit setting due to insufficient winter chill [9,10,11,12]. Therefore, it is important to study the flowering time control in perennial species and understand their adaptation mechanisms in synchronizing the timing of floral bud breaking and reproduction with local climate [10, 13, 14].

Phosphatidylethanolamine-binding proteins (PEBPs) form a superfamily of genes containing a PEBP domain, which is highly conserved across taxa, from bacteria and insects to mammals and plants [15,16,17]. Mammalian PEBPs are globular proteins composed of a functional binding site for acetate, phosphate groups and phosphorylethanolamine [18, 19]. Plant PEBP homologs share similar conserved motifs, except their C-terminal part is deleted [20, 21]. Animal PEBP proteins were reported to function as serine proteases or Raf kinase inhibitors, controlling cell growth and differentiation [22,23,24,25]. In plants, PEBP genes are central regulators in determining the flowering time, plant architecture and seed germination [26,27,28,29,30]. In angiosperms, members of the PEBP family fall into three clades of genes: FLOWERING LOCUS T (FT), TERMINAL FLOWER 1 (TFL1) and MOTHER OF FT AND TFL1 (MFT) [31, 32]. It was reported that MFT-like genes exist in both basal land plants and seed plants, while FT-like and TFL1-like genes were only found in gymnosperms and angiosperms, indicating that the MFT clade might be the evolutionary ancestor to FT-like and TFL1-like genes [32, 33]. Despite extensive sequence similarity among PEBP members, their functions have diverged from each other [34].

FT and TFL1 are two major PEBP proteins that are well studied in Arabidopsis and in many other plant species [35,36,37,38]. In Arabidopsis, FT acts as a floral signal transducer, moving from leaves to the shoot apical meristem to promote flowering, while TFL1 maintains inflorescence meristem identity in shoot apex by antagonizing FT functions [39,40,41]. The balance of FT and TFL1 modulates floral transition and inflorescence architecture by affecting determinacy of meristem identity [30, 42]. FT and TFL1 share  60% of their amino acid sequence identity, but only a few amino acid changes can convert FT from a floral promoter to a TFL1-like floral repressor [37, 43]. In addition to FT and TFL1, the Arabidopsis PEBP gene family includes MOTHER OF FT AND TFL1 (MFT), TWIN SISTER OF FT (TSF), BROTHER OF FT AND TFL1 (BFT), and CENTRORADIALIS (CEN) [27]. MFT integrates abscisic acid (ABA) and gibberellic acid (GA) signaling pathways and acts in a PIF1-dependent manner to repress seed germination under far-red light [28, 44]. TSF encodes the closest homolog of FT and resembles FT as a floral inducer under non-inductive SD conditions [45]. BFT and CEN are two floral repressors in Arabidopsis, and the overexpression of either one resulted in a late flowering phenotype similar to plants overexpressing TFL1 [46,47,48].

Although the PEBP gene family has been recognized as key floral regulators in model species, their molecular evolution and function remains less clear in woody perennials. The Rosaceae family consists of over 2500 species from approximately 90 genera, most of which are native to temperate zones around the world [49,50,51]. Prunus is a large genus belonging to the tribe Amygdaleae and contains about 430 species, many of which are important fruit crops, such as plums, cherries, apricots and peaches [52]. Additionally, Prunus includes a large number of spring-blooming trees with high ornamental and economic value. Prunus mume is one of the earliest flowering species, which blooms in late winter or early spring, followed by apricots, peaches, cherries and plums that flower during March to April. Apple and pear trees from the tribe Maleae bloom much later, around April to May in Northern China [53]. With the divergent flowering times among Rosaceae tree species, it is of great interest to investigate the evolution of PEBP family genes and their functional roles in governing flowering time among Rosaceae tree species.

Here, we provide a systematic study on the molecular evolution and function of the PEBP gene family in Rosaceae tree species. We identified 56 PEBP family genes across nine Rosaceae species and analyzed the sequence conservation, protein motifs, gene structures, and codon usage patterns of these genes. We then performed genome synteny and duplication analysis, along with nonsynonymous/synonymous substitution (dN/dS) ratio tests, to determine the evolutionary trajectory of PEBP family genes. We also analyzed the spatial and temporal expression patterns of PEBPs across tissues and in floral buds from floral initiation to bud blooming. Furthermore, we performed weighted gene co-expression network analysis (WGCNA) to determine the FT coexpressed genes in P. mume. In summary, our study provides insight into the molecular evolution of PEBP genes among Rosaceae tree species and adds information regarding their function in regulating floral bud development and blooming in woody perennials.

Results

Characterization of PEBP genes in Rosaceae species

By combining HMM and BLAST searches, we identified 56 PEBP-like proteins across nine Rosaceae tree species (Table 1). Each putative gene was validated by blasting against SMART, Pfam and NCBI CDD to ensure that they contained complete PEBP domain. We then assigned all Rosaceae PEBPs to their closest Arabidopsis homologs (Fig. 1; Table 1). In total, these Rosaceae PEBPs included 12 FT/TSF-like, 11 TFL1-like, 11 CEN-like, 10 MFT-like and 12 BFT-like genes (Table 1). TFL1 and CEN-like proteins showed the highest identities of 72.25–80.0% with their Arabidopsis orthologs, while BFT-like proteins showed the lowest identities of 62.07 to 67.82% compared with AtBFT. Five to six PEBPs were detected among Prunus species, while the average number of PEBPs almost doubled in M. domestica and Pyrus communis (Table 1). The duplicated paralogous gene pairs, such as MdTFL1 and MdTFL2, PcTFL1 and PcTFL2, were retained in the genomes of M. domestica and Pyrus communis, while only one copy of MFT was present in both species (Table 1).

Table 1 Detailed information of PEBP genes from A. thaliana and nine Rosaceae species
Fig. 1
figure1

Sequence alignment of 62 PEBP family proteins from nine Rosaceae species and A. thaliana. The sequences were aligned using Muscle. The conserved protein motif 14-3-3 interaction interface and anion-binding site are underlined in pink and purple, respectively [29]. A, B, C, and D represent four segments in exon 4 [35], which are underlined in orange, blue, green and brown

Phylogenetic analyses

Phylogenetic trees were constructed based on protein sequence alignment of Arabidopsis and Rosaceae PEBPs using three approaches: the neighbor-joining, maximum likelihood, and Bayesian inference methods (Fig. 2; Additional file 1: Fig. S1). All three phylogenetic trees shared similar topologies (Fig. 2; Additional file 1: Fig. S1). The phylogenetic trees showed that the 62 PEBP proteins can be clustered into three major clades, which are the FT-clade, TFL1-clade, and MFT-clade (Fig. 2). The FT-clade could be further split into FT/TSF-like genes and BFT-like genes, and the TFL1-clade can be split into TFL1-like and CEN-like subfamily genes (Fig. 2). Within each subfamily, the genes of Prunus species first group closely together, then group with the genes of Maleae species including M. domestica and Pyrus communis, and finally group with genes of R. occidentalis and Arabidopsis (Fig. 2). Among Prunus PEBPs, proteins within the same subgenus tend to group together, for example, P. dulcis and P. persica from the Amygdalus subgenus, P. armeniaca and P. mume from the Prunus subgenus, and proteins of P. yedoensis and P. avium from the Cerasus subgenus (Fig. 2). The duplicated paralogous gene pairs from M. domestica and Pyrus communis within the TFL1, CEN, and BFT subfamilies were grouped separately, for example, PcTFL1-MdTFL1 and PcTFL2-MdTFL2 form separate clusters, rather than genes of the same species grouping together (Fig. 2).

Fig. 2
figure2

Phylogenetic tree of PEBPs from Rosaceae species and A. thaliana constructed by the neighbor-joining method. All PEBP proteins can be clustered into three clades and five subfamilies

Structural analysis of PEBP family genes

Rosaceae PEBP family genes displayed conserved gene structures and high amino acid sequence similarity (Fig. 3; Additional file 2: Fig. S2). The length of the coding regions of PEBPs ranged from 507 to 576 bps, with FT-like genes falling between 522 to 543 bps, MFT-like genes between 507 to 576 bps, BFT-like genes between 519 to 525 bps, TFL1-like genes between 516 to 519 bps, and CEN-like genes between 519 to 522 bps. All PEBP genes have a rather loose gene structure consisting of four exons and three introns (Additional file 2: Fig. S2). For example, BFT-like genes harbor the shortest intron total lengths, ranging from 522 to 534 bp (Additional file 2: Fig. S2).

Fig. 3
figure3

Sequence conservation of the PEBP protein domains within 62 PEBP genes. a Phylogenetic analysis of sequences within the PEBP domain; b Five major motifs were predicted with MEME and visualized with TBtools

Sequence alignment revealed a high degree of conservation across the entire protein and within the PEBP domains present in all 65 genes (Figs. 1, 3). The phylogeny structure inferred from the alignment of PEBP domains was generally in accordance with that of the whole protein sequence alignment, suggesting the PEBP domain as the major factor driving the evolution of Rosaceae PEBPs (Fig. 3a). Five motifs covering 160 amino acids were identified by the MEME program among Rosaceae PEBP proteins (Fig. 3b). Among these, Motifs 1, 2, 4, and 5 together spread over the whole PEBP domain (Fig. 3b). Motifs DPDXP (Asp-Pro-Asp-X-Pro) and GIHR (Gly-Ile-His-Arg), which are essential for anion-binding activity, were present in the fourth exon of the PEBPs (Fig. 1). We also found that residues distinguishing FT-like from TFL1-like proteins were conserved among the two gene lineages (Additional file 3: Fig. S3). Previously reported key residues conferring the flowering-promoting role of FT including V76, Y91, E115, L134, Y140, G144, W145, Q147, and N159 were present in all Rosaceae FT-like proteins (Additional file 3: Fig. S3) [20, 43, 54]. The corresponding residues (I/T)76, H91, E115, (K/N/T)134, (F/N)140, (P/S)144, S145, D147, and D159 were found in all TFL1-like proteins (Additional file 3: Fig. S3). Residues determining the 14-3-3 receptor binding interface (R68, F107, R137) were shared by both protein types (Additional file 3: Fig. S3).

Microsynteny and duplication analysis of PEBP genes

To understand the evolution origin of PEBP family genes, we performed inter- and intra-genomic synteny analysis with MCScanX for Arabidopsis and seven Rosacea species with chromosome-level genome assemblies. We observed large interspecies collinear blocks between four Prunus species, P. avium, P. persica, P. armeniaca, and P. mume, which indicates high level of macrosynteny among Prunus species (Additional file 4: Fig. S4). The genome comparisons between R. occidentalis and M. domestica and between M. domestica and P. avium revealed large-scale chromosomal rearrangements including translocation and fusion-fission events that possibly occurred during the genome evolution of Rubus, Malus, and Prunus genera (Additional file 4: Fig. S4). Based on intra-genomic comparisons, we classified the duplication origin of orthologous gene pairs for Arabidopsis and other Rosacea species (Additional file 11: Table S1). Among all duplication types, whole-genome duplication (WGD)/segmental duplication was the major type for M. domestica, tandem duplicated genes were mostly found in A. thaliana, P. armeniaca, and P. persica, and dispersed duplication events were enriched in the genomes of R. occidentalis, P. mume, P. avium, and P. dulcis (Additional file 12: Table S1).

Furthermore, we characterized the duplication modes of PEBP family genes across species (Additional file 12: Table S2; Fig. 4; Additional file 5: Fig. S5). In Arabidopsis, R. occidentalis, and P. armeniaca, all PEBP gene members were predicted to be originated from dispersed duplications (Additional file 5: Fig. S5; Additional file 12: Table S2). In four Prunus species, FTs, MFTs, and BFTs were classified as having dispersed duplication, while TFL1-like and CEN-like genes were classified as exhibiting WGD/segmental duplication (Additional file 5: Fig. S5; Additional file 12: Table S2). The inter-genomic comparison of Prunus species confirmed that TFL1 and CEN genes were within shared syntenic blocks between species, indicating a shared duplication origin of TFL1-like and CEN-like genes in Prunus species (Fig. 4). Within the genome of M. domestica, we detected seven syntenic blocks consisting of three WGD/segmental duplication gene pairs, including MdTFL1-MdTFL2, MdCEN1-MdCEN2, and MdBFT1-MdBFT2, and two dispersed duplication events related to MdFT and MdMFT (Additional file 5: Fig. S5; Additional file 12: Table S2). The inter-genomic comparisons between M. domestica and R. occidental and between M. domestica and P. avium also confirmed that the duplicated gene pairs MdCEN1-MdCEN2, MdTFL1-MdTFL2, and MdBFT1-MdBFT2 are likely resulted from an independent WGD event unique to the Malus tribe (Fig. 4).

Fig. 4
figure4

Inter-genomic synteny blocks related to PEBP family genes in A. thaliana, R. occidentalis, M. domestica, P. avium, P. persica, P. armeniaca and P. mume. Chromosomes of Rosaceae species are labeled as Ro, Md, Pv, Pp, Pa, and Pm and are colored differently. We used purple, red, orange, green, and blue lines to connect collinear blocks containing MFTs, FTs, CENs, TFL1s, and BFTs, respectively

Codon usage bias and other gene parameters

We observed differentially preferred codons and different gene features across five Rosaceae PEBP gene lineages (Fig. 5; Additional file 6: Fig. S6; Table 2). For arginine, codons AGA and AGG were most frequently used by all lineages compared with other codons (Additional file 6: Fig. S6). Codon UCC encoding serine was mostly used in CEN-like and MFT-like proteins, while codon UCU was mostly employed by the FT and TFL1 lineages (Additional file 6: Fig. S6). We also observed significant differences in other gene parameters among different PEBP lineages (all Kruskal–Wallis tests pval < 0.01) (Fig. 5; Table 2). The codon adaptive index (CAI) of MFT and TFL1 genes is significantly larger than those of the other gene groups (Kruskal–Wallis test pval = 2.75e−7) (Fig. 5; Table 2). In contrast, the effective number of codons (ENC) estimated for MFT and TFL1 genes is much lower than those of the other groups (Kruskal–Wallis test pval = 0.005) (Fig. 5; Table 2). The average ENC values ranging from 51.24 to 54.95 indicated weak codon bias among PEBP genes. Analysis of the GC content revealed that MFT lineage genes had much higher GC content indices compared to other genes (Fig. 5; Table 2). In contrast, TFL1 and BFT lineage genes appear to have lower GC1% and GC3% but relatively higher GC2% compared to other groups (Fig. 5; Table 2). All gene parameters showed no variation among species (Additional file 7: Fig. S7; Kruskal–Wallis test pval > 0.05). Strong pairwise correlations between gene parameters were observed (Additional file 13: Table S3). For example, the CAI was positively correlated with the total GC% and GC3% (both correlation coefficient r ≥ 0.56), but was negatively correlated with the ENC (r = − 0.62) (Additional file 13: Table S3). On the other hand, the ENC displayed a negative correlation with the total GC content and GC3% (Additional file 13: Table S3).

Fig. 5
figure5

Comparison of gene parameters including the a CAI, b total GC%, c ENC, d GC1%, e GC2%, and f GC3% estimated for BFT, CEN, FT, MFT, and TFL1 genes

Table 2 Average gene parameters estimated for FT, TFL1, CEN, BFT, and MFT gene lineages

Molecular evolution of different PEBPs lineages

To investigate the evolution of PEBP genes in Rosaceae species, we performed selection scans on coding sequences of all PEBPs using the branch model, site model, and branch site model in the CODEML program of PAML (Table 3; Additional file 14: Table S4; Additional file 15: Table S5). Branch models with different ω parameters specified for foreground lineages (i.e., FT-like, TFL1-like, CEN-like, and BFT-like lineages and FT/TFL1 clades) were compared with the fixed ratio model (Additional file 14: Table S4). The likelihood ratio tests (LRT) on models specifying individual lineages of FT, TFL1, CEN, and BFT genes as the foreground branch showed no significant difference in ω between the foreground and background branch (P > 0.05) (Additional file 14: Table S4). However, the LRT test on the branch model specifying the FT and TFL1 clades as the foreground branch suggested significant divergence among FT/TFL1 and MFT clade genes (P < 0.001) (Additional file 14: Table S4). We then applied the site model LRT test and detected signs of positive selection among sites of PEBP proteins (Additional file 15: Table S5). The branch-site LRT tests further revealed strong positive selection within TFL1 lineage and slight positive selection within FT lineages at specific protein sites (Table 3). The Empirical Bayes model suggested modest selection at positions 19 and 106 when FT lineage was set as the foreground branch and at positions 11 and 18 when TFL1 lineage was set as the foreground branch (Table 3). We further validated the results by performing selective pressure analysis on five gene lineages separately with the software Selecton. Only the FT and TFL1 lineages showed the signature of positive selection, in which residues 40N, 56N, 128S, and 181L in the FT lineage (with RoFT1 as the reference gene) and 4T, 73V, 134P, 141S, 157L, and 161S in the TFL1 lineage (with PyTFL as the reference gene) were mostly selected (Fig. 6). In contrast, the genes of the other three lineages all showed signs of purifying selection across most sites (Additional file 8: Fig. S8).

Table 3 Parameter estimates and likelihood values for branch-site models among sites and lineages of PEBP
Fig. 6
figure6

Selective pressure analysis of FT and TFL1 lineage genes identified positively selected sites. Amino acids colored with yellow/purple indicate sites of positive/purifying selection, respectively

Cis-acting element analysis of the FT promoter

We extracted the 2000 bp region of FT genes and scanned for putative cis-elements by searching against the PlanPan and the PlantCARE databases (Table 4). We compared the type and copy number of cis-elements for 11 FT genes from A. thaliana, P. trichocarpa, M. domestica, Pyrus communis, R. occidentalis, P. armeniaca, P. avium, P. mume, and P. persica (Table 4). Within the promoter region of the investigated FTs, three to ten CCACA boxes (binding site for CO) were identified across nine species, while none were found for PtFT2 (Table 4). CArG boxes, the binding site for the MADS-box transcription factor, were present in all FT promoters, among which the AtFT promoter contained the most (Table 4). Light-response elements including the G-box, AE-box, GATA-motif, GT1-motif, and TCT-motif were present within all FT genes but in different types (Table 4). In addition, binding sites for MYB, MYC transcription factor, ethylene-responsive transcription factor, and abscisic-acid responsive element (ABRE) were present in all FT promoters (Table 4). Gibberellin-responsive elements of different types were detected in FT promoters, with GARE-motif in the promoters of MdFT, PcFT1, RoFT1, PvFT and P-box in the promoters of AtFT, PtFT2, PaFT, PvFT, PmFT, and PpFT. We also observed some cis-elements with species-specific distribution patterns. For example, the low-temperature responsiveness (LTR) element was only detected within the promoters of AtFT, PcFT1, RoFT1, PvFT and PmFT (Table 4). The W-box, which is the binding site for WRKY transcription factor, was detected exclusively in RoFTs, PtFT2, and Prunus FT promoter regions (Table 4).

Table 4 Summary statistics of putative cis-elements present in the 2 kb upstream promoter region of FTs across nine species

Tissue-specific expression patterns of PEBPs

To explore the functional roles of PEBP genes, we examined their expression patterns in different tissues of four Rosaceae species, P. persica, P. mume, P. yedoensis, and R. occidentalis (Fig. 7a–d). In general, we observed a differentiated expression preference of PEBP genes across different tissues (Fig. 7). Among the five PEBP subfamilies, FT-like and TFL1-like genes were expressed in both vegetative tissues such as leaf and stem, and reproductive organs such as flower bud and fruit (Fig. 7). The transcription of CENs, as the closest paralogs of TFL1, was barely detected in any organs, except in the root tissues of P. mume (Fig. 7). MFT was only detected in seed embryos of P. persica and fruit tissues of P. yedoensis and R. occidentalis (Fig. 7). BFT was detected in the fruit tissues of all species but was relatively highly expressed in leaf and stem tissues in P. yedoensis and R. occidentals, respectively (Fig. 7). We validated the tissue-specific expressions of five PEBP genes by real-time quantitative PCR (qRT-PCR) in P. mume (Additional file 9: Fig. S9). PmFT is highly expressed in floral buds compared with its expression in leaf and stem, which is consistent with result of the above tissue transcriptome sequencing in P. mume (Fig. 7; Additional file 9: Fig. S9). PmTFL and PmCEN were relatively highly expressed in root tissues (Additional file 9: Fig. S9). PmBFT and PmMFT was barely detected in the four examined tissue types (Additional file 9: Fig. S9). The somewhat inconsistent tissue-specific expression patterns of PEBP orthologs across examined species are likely a result of non-uniformity in the sampling time, plant physiological state, and tissue specificity across four independent studies. Despite the inconsistency, the divergent expression of PEBP members across different tissue types indicates significant functional differentiation of PEBP gene lineages.

Fig. 7
figure7

Tissue-specific expression of PEBP genes in four Rosaceae species including a P. persica, b P. yedoensis, c P. mume and d Rubus occidentalis

Expression analysis of PEBP genes during floral bud development in P. mume

We analyzed the expression of PEBP genes in flower buds of different developmental stages from July 10th, 2019 to January 12th, 2020 by qRT-PCR analysis. The expression of PmFT first decreased as the bud initiated the floral meristem from July to August, increased as floral organ initiated and developed (from August to October), slightly decreased during bud dormancy, and then significant increased as the floral bud exited dormancy and bloomed (Fig. 8). PmBFT maintained a low expression level throughout the whole process, with only a minor increase during floral bud development in August and September (Fig. 8). The other PEBP members retained barely detected expression levels in floral buds of all developmental stages (Fig. 8). These results imply that PmFT is possibly the primary PEBP member participating in regulating floral bud development and bud flushing in P. mume.

Fig. 8
figure8

Relative expression levels of PEBP genes during floral bud development and bud blooming in P. mume

Co-expression network analysis of FT during the blooming process in P. mume

To explore the regulatory network of FT in flowering regulation in trees, we reanalyzed the transcriptome changes of P. mume during dormancy release and the floral bud opening process [55] and performed a weighted co-expression network analysis (WGCNA). We identified 23 modules with distinct expression patterns (Additional file 10: Fig. S10a). Module-trait association analysis revealed four modules, ‘brown’, 'turquoise', ‘dark green’, and ‘salmon’, associated with the progression of bud flushing (R2 > 0.8). Among them, module ‘brown’ showed the strongest correlation with the FPKM of PmFT (Additional file 10: Fig. S10b). The ‘brown’ module genes were significantly enriched in biological processes including cell cycle (GO: 0007049), flower development (GO: 0009908), glucan metabolic process (GO: 0009251), auxin transport (GO: 0060918), and responses to abiotic stimulus (GO: 0009628). We further identified the top 50 genes most associated with PmFT and 15 known flowering-related genes such as PmLFY, PmAP1, and PmCOL (Additional file 16: Table S6) [56, 57]. Among genes in the ‘brown’ module, SVP (SHORT VEGETATIVE PHASE), SOC1 (SUPPRESSOR OF OVEREXPRESSION OF CO 1), GI (GIGANTEA), and CIB1 (CRYPTOCHROME-INTERACTING BASIC-HELIX-LOOP-HELIX 1) were previously identified as key players in the FT-dependent floral regulation in Arabidopsis [58, 59] (Fig. 9a). Four tandem-duplicated PmDAMs (PmDAM1, PmDAM4, PmDAM5, PmDAM6) from the ‘brown’ module also exhibited expression patterns negatively correlated with that of PmFT (Fig. 9a, b). The expression patterns of other known floral regulators such as COL (CONSTANS-LIKE) from the ‘turquoise’, LHY1 (LATE ELONGATED HYPOCOTYL 1) and AP1 (APETALA1) from ‘dark green’ module were not highly correlated with PmFT (R2 < 0.62) (Fig. 9b). PmFT showed a relatively weak transcription level in endodormant floral buds (Fig. 9b). As the floral bud continued accumulating chilling units and exiting dormancy, PmFT expression significantly increased and showed the highest expression in flushing buds (Fig. 9b). PmCIB1 and 37 other genes showed similar expression patterns to that of PmFT, while PmPHYB (Pm008367), PmGI, PmLHY, PmCOL, PmSVP, PmSOC1, and four PmDAMs displayed contrasting expression patterns, with their expression decreasing as the floral buds exited endodormancy (Fig. 9b). The expression patterns of FT and its coexpressed genes were further verified by qRT-PCR analysis (Fig. 9c).

Fig. 9
figure9

Co-expression network of FT during floral bud blooming in P. mume. a Cytoscape visualization of candidate genes co-expressed with PmFT during dormancy release. Candidate genes from the ‘brown’, ‘dark green, ‘green-yellow, ‘turquoise’, and ‘cyan’ modules are colored in brown, green, green-yellow, turquoise, and cyan, respectively. The circle size represents the significance of gene expression correlation with PmFT. b Expression patterns of PmFT and putative co-expressed genes during floral bud blooming. c Relative expression of PmFT and putative co-expressed genes verified by qRT-PCR analysis

Discussion

Evolution trajectory of PEBP family genes in Rosaceae genomes

PEBPs form an ancient gene family central to many plant developmental processes, including floral transition, plant architecture, and seed germination [30, 32, 60]. In Arabidopsis, the PEBP family constitutes six genes grouped into three distinct clades, FT-like (FT and TSF), TFL1-like (TFL1 and CEN), and MFT-like genes [31]. Though previous studies have characterized the functions of PEBP family genes in model plants, none have focused on a comparative analysis of the PEBP family in tree species. Our study conducted a systematic search across nine Rosaceae genomes and identified 56 PEBP family genes orthologous to six Arabidopsis genes, FT/TSF, TFL1, CEN, BFT, and MFT. The number of PEBP family members in Prunus species (chromosome 2n = 2x = 16) was approximately the same as that in Arabidopsis (five to six copies), while PEBP members were expanded in M. domestica and Pyrus communis (chromosome 2n = 2x = 34). Genome synteny and duplication analyses together supported that duplicated ortholog pairs MdTFL1-MdTFL2, MdCEN1-MdCEN2, and MdBFT1-MdBFT2 are likely originated from a recent whole-genome duplication (WGD) event that occurred in the Maleae clade after splitting from Prunus [61]. However, only one copy of MdFT, MdMFT, and PcMFT was retained in apples and pears, indicating that the duplicated copy may have been lost during species evolution after the WGD [62]. The duplication mode analysis also suggested a shared origin of TFL1 and CEN from segmental or WGD duplication in Prunus species (Additional file 12: Table S2). Previous studies reported that the angiosperm TFL1-like gene experienced duplication after splitting from basal angiosperms, followed by functional divergence, resulting in TFL1 and CEN gene lineages in eudicots [63]. Given the conserved sequence alignment of Prunus TFL1/CEN orthologs with other Rosaceae species, it is unlikely that Prunus TFL1/CEN arose from a recent segmental duplication or WGD unique to Prunus species. Therefore, the syntenic relationship may have been caused by the preservation of genomic segments containing TFL1, CEN, and their neighboring genes through rounds of chromosome rearrangements during Prunus species evolution. In Arabidopsis, the TSF gene, which is a homolog of FT, highly resembles FT in its coding sequence and flowering promoting role [64]. The absence of TSF in the Rosaceae genome suggests that the gene duplication of FT/TSF possibly occurred in Brassicaceae after splitting from their common ancestors [65].

The PEBP gene family experienced two ancient duplications, giving rise to three types: FT-like genes promoting flowering, TFL1-like genes repressing flowering and maintaining indeterminate state of meristems, and MFT-like genes controlling seed germination [17, 27, 32]. The phylogenetic analysis suggests that Rosaceae PEBPs can be clustered into three distinct clades (FT, TFL1, and MFT), which is consistent with other species [17, 27, 32]. The FT-like clade can be further divided into FT and BFT lineages, and the TFL1-like clade can be divided into TFL1 and CEN lineages. Based on maximum-likelihood test on branch models specifying different gene lineages (FT, TFL1, CEN, MFT, and BFT) as the foreground branch, we detected no evidence of positive selection acting on any of them. However, we observed significant selection acting on FT/TFL1 clade genes with the MFT clade specified as the background branch, which supports the theory that functional divergence of the FT/TFL1 clade occurred after splitting from the MFT clade [33]. Through likelihood ratio tests on branch-site models, we detected a few slightly selected codons within the FT lineage and a few strongly selected codons in the TFL1 lineage, which is consistent with results of Selecton analysis on individual lineages. In summary, these results indicate that adaptive evolution is driving the divergence of the FT and TFL1 clades from the MFT clade, as well as the diversification among FTs and TFL1s in Rosaceae species. These result are consistent with a previous study reporting that positive selection on FT-like genes especially within the fourth exon is driving their divergence from MFT and TFL1 clade [17]. We also observed strong purifying selection constraining protein evolution within the MFT, CEN, and BFT lineages in Rosaceae species. However, this does not rule out the possibility of positive selection acting on a few codons masked by strong purifying selection in preserving the other sites [17].

Additionally, we examined the codon usage patterns of PEBP genes across Rosaceae species. Codon usage bias refers to the nonrandom choice of synonymous codons in specific genes or species and can affect the translation efficiency and accuracy, protein folding, and biological functions [66, 67]. The codon usage pattern usually reflects the balanced effect of mutation pressure and selection constraints during long-term evolution [68, 69]. Several codons for amino acids were differentially preferred across five PEBP lineages. Among all codons, the most frequently used codon for arginine was AGG for FT, CEN, and MFT and AGA for the BFT and TFL1 lineages (Additional file 6: Fig. S6). Several other codons, including TCC, TCA, TCT for serine and CCT for proline, were preferred by specific PEBP gene lineages, indicating differentially selected codons by different PEBP gene lineages. To further understand the factors influencing codon usage patterns, we compared the GC content, gene length, CAI, and ENC of different PEBP lineages and species. The CAI measures the optimal codon usage for a gene and is commonly used as an index for the expression level [70]. The ENC has been widely used to determine the level of codon bias for individual genes [71]. We observed significant differences in these gene features estimated for different gene lineages but not for species. Despite the differences, all genes had a relatively high CAI (range 0.81–0.87) and moderate ENC (above 47), indicating high translational efficiency and slightly biased codon usage among PEBP genes. Furthermore, the strong pairwise correlations between ENC and GC content, ENC and CAI indicate that the nucleotide composition and gene expression level are two factors possibly contributing to the differentiated codon preference among different PEBP gene lineages [69].

Functional role of FT/TFL1 genes in Rosaceae tree species

Structural analysis of Rosaceae PEBP proteins revealed a highly conserved gene structure and amino acid sequence, especially within the PEBP functional domain (Figs. 1, 3; Additional file 2: Fig. S2). All PEBP family genes shared a common gene structure with exactly four exons of similar sizes. Among the conserved protein motifs, the anion-binding D-P-D-x-P and G-x-H-R motifs are important for the conformation of the ligand binding site in PEBP proteins [72]. Mutations close to this region may affect the binding of FT protein with phosphate ions and thus alter its interaction with FD (FLOWERING LOCUS D) [73]. Segment B on exon 4 encodes an external loop, and together with its adjacent segment C, determines the opposite functions of FT and TFL1 in Arabidopsis [35]. Another key protein motif is the 14-3-3 binding domain that is essential for FT/TFL1 interaction with 14-3-3 receptors to modulate flowering [20]. Key residues within these motifs are critical in determining FT/TFL1 functions. For example, the substitution of an amino acid (replacing His-88 in TFL1 with Tyr) can convert TFL1 into a floral promoter [37]. In another study, specific mutations at four residues—Glu-109, Trp-138, Gln-140, and Asn-152—converted FT into a TFL1-like repressor [43]. The amino acids at each of these critical positions were highly conserved and specific to FT-like and TFL1-like proteins, which suggests that the floral promoting and repressing role of FT/TFL1 genes in Rosaceae species is possibly conserved.

Recent molecular studies have characterized the function of Rosaceae FT/TFL-like genes in several Rosaceae perennials [33]. The overexpression of MdFT in both Arabidopsis and apple lead to precocious flowering [74]. The ectopic expression of PmFT and RoFT in tobacco leads to extremely advanced flowering [75]. Similarly, the late-flowering phenotype of Arabidopsis ft mutant can be rescued by overexpressing PpFT, indicating the conserved floral promoting role of FT in examined Rosaceae species [76]. On the other hand, prolonged vegetative growth and a late-flowering phenotype were observed for transgenic Arabidopsis/tobacco overexpressing PpTFL1, PmTFL1, RoTFL1, MdTFL1-1/2, suggesting that the Rosaceae TFL1-like genes can complement the TFL1 function in Arabidopsis [77,78,79].

Despite the conservative function of Rosaceae FT/TFL1-like genes in herbaceous plant systems, their regulatory roles in perennial trees may differ. For example, two homologs of PcFTs showed differed annual expression patterns in the apical buds of Pyrus communis [80]. The ectopic expression of PcFT2 caused early flowering in tobacco but delayed dormancy and leaf senescence in M. domestica [80]. Another study in pears reported that the expression of FTs was not induced in the reproductive meristem prior to floral initiation, while the transcripts of TFL1s rapidly decreased and maintained a very low level, indicating the essential role of TFL1 in floral induction in Pyrus pyrifolia [36]. In our study, the minimal level of TFL1 throughout all floral bud stages may indicate that the repression of TFL1 is necessary for determinate floral meristem identity and terminal flower formation during floral bud development in P. mume. The multifaceted role of FT/TFL1-like genes was also observed in other tree species [33]. In poplar, PtFT1 functions as a floral promoter activated by chilling temperatures, while vegetative growth and dormancy breaking are promoted by PtFT2 [81]. Plum trees transformed with PtFT1 displayed a shrub-like growth habit, a reduced chilling requirement, and insensitivity to short-day signals [82]. In gymnosperms, FT-like genes exhibited contrasting roles in regulating growth cycling and bud setting [83]. For example, expression of FT/TFL1-like genes in Norway spruce (PaFTL2) and Scots pine (PsFTL2) increase during bud setting in autumn and decrease during bud bursting in the next spring [84,85,86]. Thus, FT/TFL1-like genes may undertake some novel functions concerning floral transition, plant architecture, and growth-dormancy cycling during the evolution of tree species.

Regulatory role of FT in promoting bud break and blooming in perennial trees

Flowering is a major developmental process that is key to the fitness and reproduction of higher plants [87]. Plants have synchronized their seasonal timing of flowering with favorable environmental conditions to ensure sexual reproduction success and seed production [87, 88]. The regulation of flowering times requires an intricate network of signaling pathways, which has been studied in many plant species but is best characterized in Arabidopsis [57, 87, 89]. FT functions as a gene hub integrating five major floral induction pathways, including the photoperiodic pathway, vernalization pathway, autonomous pathway, gibberellin pathway, and age pathway [56, 59]. In Arabidopsis, the transcription of FT is activated by the transcription factor CONSTANS (CO), which is affected by the circadian regulatory GI [90, 91]. The GI-CO-FT module not only is used to regulate photoperiod-dependent flowering in Arabidopsis and temperate cereals [92, 93] but also showed a conserved function in regulating short-day induced bud dormancy in poplar [94]. In addition to CO, SVP, FLC (FLOWERING LOCUS C), and PIF4 (PHYTOCHROME INTERACTING FACTOR 4) from the vernalization pathway can also regulate FT transcription through directly binding the FT promoter or intronic regions [90, 95,96,97]. Upon induction by long-photoperiod signals, FT, together with other floral pathway integrators SOC1 and LFY (LEAFY), activates floral meristem identity genes such as AP1, APETALA2 (AP2), FRUITFULL (FUL), CAULIFLOWER (CAL), and LFY, which convert the vegetative meristem to floral meristem in Arabidopsis [59, 98, 99].

Though flowering regulation is well understood in model species, it is still unclear in temperature tree species. Unlike annual or biennials, many trees in temperate environments initiate floral buds in the preceding summer, cease growth in autumn, with floral buds remaining dormant during winter, and then bloom early in spring after exposure to chilling temperatures [6, 10]. Therefore, perennial flowering marks the event of the floral bud exiting dormancy and flushing instead of the time of floral meristem initiation in annual species [6]. So far, many studies on floral bud breaking regulation have been reported; however, the molecular mechanism is still far from complete. Apart from regulating floral initiation, FT has been suggested to participate in regulating bud dormancy in temperate trees [100]. Poplar exhibited constitutive expression of FT1 initiated flower-like structures directly from tissue culture and showed delayed growth cessation in short-days [81, 94], while FT2 was predominantly expressed during vegetative growth and is likely responsible for growth cessation and vegetative bud set [81]. Moreover, Rinne et al. (2011) reported that FT is hyper-induced during bud breaking in poplar, indicating that FT may also participate in regulating dormancy release in poplar [101]. In pear, chilling reduces the expression of DAM genes, which are well-known floral repressors, releasing the repression of FT and promoting floral bud breaking [102, 103]. Our expression analysis confirmed that FT is significantly induced during chilling-mediated floral bud breaking in P. mume.

To further understand the regulatory module of FT during floral bud breaking, we used WGCNA and identified a number of candidate genes whose expression patterns strongly correlated with FT in P. mume. Among these candidates, PmDAM1, PmDAM4, PmDAM5, and PmDAM6 were found to be downregulated during the progression of bud breaking. Another MADS-box gene PmSVP displaying a similar expression pattern to that of PmDAMs was reported to maintain bud dormancy in apples [104]. Thus, PmDAMs and PmSVP may function as FT repressors in the same manner as in Arabidopsis by binding to the CArG box in the promoter region of PmFT [105]. A number of genes previously identified upstream FT, including PmCOL, PmGI, and PmCIB1, were found to be induced by chilling in endodormant buds before the activation of PmFT. These genes may act directly or indirectly to activate FT expression during dormancy release in P. mume. We also observed that some known FT regulated genes, namely, AP1, SOC1, and LFY, peaked before the induction of FT, indicating their functional role during flower bud development prior to bud breaking [6, 106]. Additionally, a number of FT co-expressed genes were annotated to pathways that did not show relatedness to bud breaking or flowering in previous studies. Future functional studies are required to characterize the regulatory mechanisms of FT in floral induction and bud breaking in Rosaceae tree species.

Conclusions

In this study, we systemically characterized the PEBP gene family in nine Rosaceae species and examined their gene structure, protein features, evolutionary trajectories, and expression profiles. The 56 PEBP genes can be divided into three major clades, namely, FT-like, TFL1-like, and MFT-like genes. We observed highly conserved protein motifs and gene structure among PEBP genes. Selection scans showed that positive selection is driving the divergence of the FT and TFL1 clades, while strong purifying selection is restraining diversification within most lineages. Expression analysis of PEBP genes suggested the essential role of FT in floral bud development and blooming. Furthermore, we identified a number of FT co-expressed genes, revealing a FT-related regulatory model in Prunus species different from those in annual or biennial plants. In summary, the comprehensive analysis of the PEBP family in our study provided evidence of structural and functional conservation of PEBP genes among Rosaceae woody perennials and provided insight into the adaptive evolution of the PEBP gene family over the evolutionary history of perennial trees.

Methods

Identification of the PEBP gene family

We obtained the most recent versions of genomes for P. persica [107], P. mume [108], P. yedoensis [109], P. avium [110], P. dulcis, P. armeniaca [111], M. domestica [112], Pyrus communis [113], and R. occidentalis [114] from GDR (Genome Database For Rosaceae) [115]. To identify the PEBP genes of each species, we retrieved the HMM model PF01161 of PBP domain from the Pfam database (https://pfam.xfam.org) and searched the genome protein databases with an e-value cutoff of 1.0 × e−5 using HMMER 3.1 software [116]. In addition, we used protein sequences AtFT (At1g65480.1), AtTSF (At4g20370.1), AtTFL1 (At5g03840.1), AtBFT (At5g62040.1), and AtMFT (At1g18100.1) downloaded from TAIR (The Arabidopsis Information Resource) (www.arabidopsis.org) as query sequences to blast against the local protein databases of nine species, and we only retained putative PEBP proteins with identities > 40% and e-values ≤ 1.0 × e−10. The genes identified by both methods were considered as candidate PEBP family genes and were then verified with SMART [117], Pfam [118], and the CDD database [119] to ensure the completeness of the PBP domain. Redundant sequences or sequences with incomplete PEBP domain were excluded from the following analyses.

Phylogenetic analysis

Multiple sequence alignment was performed using the protein sequences with software MUSCLE v3.8 [120] and was visualized with GeneDoc v2.6 [121]. Phylogenetic trees were constructed using neighbor-joining (NJ) method with MEGA7 [122], maximum likelihood (ML) analysis with RAxML v8.1 [123], and Bayesian inference (BI) with MrBayes 3.1 [124]. Bayesian inference was performed with 100,000 generations of Markov-chain Monte Carlo (MCMC) simulations, discarding the first 2500 trees as ‘burn-in’. With consistent tree topologies inferred by these three approaches, the neighbor-joining tree was chosen to display the phylogeny of Rosaceae full PEBP protein sequences. Furthermore, amino acids within the regions of predicted PEBP domains were extracted and used to construct a PEBP domain tree by the NJ method.

Gene structure and protein motif detection

The exon and intron locations of PEBP genes were analyzed by comparing the coding sequences with their genome sequences. The MEME (Multiple Expectation Maximization for Motif) online tool (http://meme-suite.org/tools/meme) was used to predict protein motifs [125]. The protein motifs were further annotated with the Pfam [118], SMART [117] and CDD [119] online tools. The chromosome distributions of PEBP genes were obtained based on genome GFF3 files. Finally, the gene structures, protein motifs, and chromosome locations were visualized with the software TBtools [126].

Microsynteny analysis and codon usage evaluation

To identify the synteny of PEBP family genes among species, we performed all-to-all BLASTP between the genomes of A. thaliana, R. occidentalis, M. domestica, P. avium, P. persica, P. armeniaca, and P. mume. We also performed self-blast by comparing protein-coding genes against their own genome using BLASTP. All BLASTP hits with e-values < 1e−10 were used as input for software MCScanX (Multiple Collinearity Scan toolkit) [127] to identify possible collinear blocks within and between genomes of different species. Based on the self-blast results, we classified the duplication origin of orthologous genes pairs including PEBP family genes with the ‘duplicate_gene_classifier’ toolkit built in MCScanX for each species. All intra/inter-genomic synteny relationships were visualized with TBtools [126].

Gene parameters including the GC content (total GC%, GC1%, GC2%, and GC3%), CAI, and ENC were computed using CAICal (http://genomes.urv.cat/CAIcal/) [128, 129]. CAI provides an estimate of directional translational selection in optimizing the codon usage patterns of genes and is used to predict highly expressed genes [70]. ENC is a number between 20 to 61 that measures the degree of codon usage bias (where ENC = 20 refers to the preference of only one codon per amino acid, while ENC = 61 refers to complete unbiased codon usage) [130]. We compared these gene parameters for FT, TFL1, CEN, BFT, and MFT gene lineages and across species using the Kruskal Wallis Test with the ‘kruskal.test’ function in R. The relative synonymous codon usage (RSCU) is defined as the ratio of the observed codon frequency to the expected frequency of all synonymous codons per amino acid and is calculated using software MEGA7 [131].

Molecular evolution of PEBP genes

To investigate the signatures of positive selection on Rosaceae PEBP genes, we extracted the coding sequences of PEBP genes and aligned them with MUSCLE v3.8 [120]. The sequence alignment was then trimmed with Gblocks [132] in ‘codon’ mode, and the resulting alignments were used to infer phylogenetic relationships with RAxML [123]. The ratios (ω) of nonsynonymous substitution sites (dN) and synonymous substitution sites (dS) were computed for each PEBP lineage gene using the branch model, site model and branch-site model with the codeml package in PAML 4.0 [133]. To test the hypothesis of adaptive evolution in specific PEBP lineages and across sites, we performed likelihood ratio tests to evaluate the fit of branch models (FT, TFL1, CEN, BFT and (FT, TFL1, CEN, BFT) set as foreground branch), site models, and branch site models. The positively selected sites were detected by Bayes Empirical Bayes analysis in PAML 4.0 [133]. To better visualize the site-specific selection on amino acids within each PEBP lineage, we performed a selection pressure test with site model M8 and visualized the results with Selecton Server [134].

Cis-element analysis of the FT promoter region

To investigate the conservation of the cis-regulatory model of FT genes across different species, we extracted the 2 kb upstream region of the start codon (ATG) and submitted the sequences to the PlantCARE [135] and PlantPan 2.0 databases [136]. The cis-acting elements predicted by both methods were integrated and considered as putative cis-acting elements.

Tissue-specific expression profiles of PEBP genes

The tissue transcriptome sequencing data of P. mume, P. yedoensis, P. persica, and R. occidentalis was retrieved from four independent studies: GSE4760162 from the GEO database [108] and SRP136962, SRA053230, and SRP149938 from the NCBI SRA database [109, 137]. The raw SRA files were first dumped to FASTQ format using SRA toolkit and preprocessed with Trimmomatic v0.38 [138] to trim off poor-quality reads. Clean paired reads were aligned with the reference genomes of P. mume, P. yedoensis, P. persica, and R. occidentalis, respectively, with software HISAT2 [139]. The genic count was computed with HTSeq [140] and normalized to RPKM with R package ‘edgeR’ [141]. The RPKM value of each PEBP gene across different tissues of P. mume, P. persica, P. yedoensis, and R. occidentalis was extracted and visualized using the ‘pheatmap’ package in R. The relative expression of PEBP genes in leaf, stem, root, and floral bud tissues was tested in P. mume using real-time PCR analysis with detailed procedure described below.

Expression analysis of PEBP genes during the flower bud development process

To further understand the functional role of PEBP genes in floral bud initiation and the bud flushing process, we performed real-time quantitative PCR analysis to examine the temporal expression patterns of PEBP genes. Lateral floral bud samples were collected from P. mume ‘Fei Lve’ tree grown in the Jiufeng sunlight greenhouse approximately every four weeks from July 10th, 2019 to January 12th, 2020. The total RNA was extracted from mixed bud samples using the E.Z.N.A.® Plant RNA Kit following the manufacturer’s instructions (Omega Bio-tek, Norcross) and was reverse-transcribed into cDNA using the PrimeScript RT reagent kit with gDNA Eraser (Takara, Japan). We performed real-time PCR experiments with at least three technical replicates on the PikoReal real-time PCR platform (Thermo Fisher Scientific, Germany). The temperature was set as follows: 95 °C for 30 s; 40 cycles of 95 °C for 5 s, 60 °C for 30 s; 60 °C for 30 s; and ending 20 °C. We used protein phosphatase 2A (PP2A) as an internal reference and calculated the relative transcription levels of target genes using the 2 − ΔΔCt method [142]. The primers used for qRT-PCR experiments are listed (Additional file 17: Table S7).

Co-expression network of FT during the blooming process in P. mume

To investigate the functional role of FT during floral bud flushing, we obtained the transcriptome data of four successive stages during dormancy release and blooming in P. mume from a previous study reported by Zhang et al. (2018) [55]. The procedure of sample collection, RNA extraction, sequencing library construction, quality control, and gene expression quantification was described in detail [55]. We normalized the gene expression and performed weighted gene co-expression network analysis with WGCNA v1.67 package in R [143]. The Dynamic Tree Cut algorithm was applied to detect gene modules (power β of 4; height cutoff of 0.3; minimal module size of 30). To identify the key modules coexpressed with FT, we calculated the module-trait association and ranked genes by their correlation with the FPKM value of PmFT. Finally, the top 50 candidate genes (R2 > 0.6) coexpressed with PmFT and 15 FT interacting factors identified in Arabidopsis flowering pathways [56, 57] were selected to construct the coexpression network of FT. The FT regulatory network was visualized with Cytoscape 3.1 [144]. The expression levels of FT and putative co-expressed genes were further validated by qRT-PCR analysis. The primers are described in the supplementary data (Additional file 18: Table S8).

Availability of data and materials

The datasets analyzed during the current study are available from the corresponding author on reasonable request.

References

  1. 1.

    Ågren J, Oakley CG, Lundemo S, Schemske DW. Adaptive divergence in flowering time among natural populations of Arabidopsis thaliana: estimates of selection and QTL mapping. Evolution. 2017;71(3):550–64.

    PubMed  Article  Google Scholar 

  2. 2.

    Gaudinier A, Blackman BK. Evolutionary processes from the perspective of flowering time diversity. New Phytol. 2019;225(5):1883–98.

    PubMed  Article  Google Scholar 

  3. 3.

    Townsend T, Albani M, Wilkinson M, Coupland G, Battey N: The diversity and significance of flowering in perennials. In: Annual plant reviews volume 20: Flowering and its manipulation. 2006. p. 181–97.

  4. 4.

    Khan MRG, Ai XY, Zhang JZ. Genetic regulation of flowering time in annual and perennial plants. Wiley Interdiscip Rev. 2014;5(3):347–59.

    CAS  Article  Google Scholar 

  5. 5.

    Huijser P, Schmid M. The control of developmental phase transitions in plants. Development. 2011;138(19):4117–29.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Wilkie JD, Sedgley M, Olesen T. Regulation of floral initiation in horticultural trees. J Exp Bot. 2008;59(12):3215–28.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Basler D, Korner C. Photoperiod and temperature responses of bud swelling and bud burst in four temperate forest tree species. Tree Physiol. 2014;34(4):377–88.

    PubMed  Article  Google Scholar 

  8. 8.

    Rohde A, Bastien C, Boerjan W. Temperature signals contribute to the timing of photoperiodic growth cessation and bud set in poplar. Tree Physiol. 2011;31(5):472–82.

    PubMed  Article  Google Scholar 

  9. 9.

    Rafferty NE, Scaven VL. Physiological effects of climate warming on flowering plants and insect pollinators and potential consequences for their interactions. Curr Zool. 2013;59(3):418–26.

    PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Abbott AG, Zhebentyayeva T, Barakat A, Liu Z. Chapter six - the genetic control of bud-break in trees. In: Plomion C, Adam-Blondon A-F, editors. Advances in botanical research. Boca Raton: Academic Press; 2015. p. 201–28.

    Google Scholar 

  11. 11.

    Dirlewanger E, Quero-García J, Le Dantec L, Lambert P, Ruiz D, Dondini L, Illa E, Quilot-Turion B, Audergon JM, Tartarini S, et al. Comparison of the genetic determinism of two key phenological traits, flowering and maturity dates, in three Prunus species: peach, apricot and sweet cherry. Heredity. 2012;109(5):280–92.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Vitasse Y, Lenz A, Körner C. The interaction between freezing tolerance and phenology in temperate deciduous trees. Front Plant Sci. 2014;5:1–8.

    Article  Google Scholar 

  13. 13.

    Ladwig LM, Chandler JL, Guiden PW, Henn JJ. Extreme winter warm event causes exceptionally early bud break for many woody species. Ecosphere. 2019;10(1):1–9.

    Article  Google Scholar 

  14. 14.

    Cattani AM, Sartor T, da Silveira Falavigna V, Porto DD, Silveira CP, Diasde Oliveira PR, Revers LF. The control of bud break and flowering time in plants. Adv Bot Res. 2018;88:277–325.

    CAS  Article  Google Scholar 

  15. 15.

    Chautard H, Jacquet M, Schoentgen F, Bureaud N, Bénédetti H. Tfs1p, a member of the PEBP family, inhibits the Ira2p but not the Ira1p Ras GTPase-activating protein in Saccharomyces cerevisiae. Eukaryot Cell. 2004;3(2):459–70.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. 16.

    Yeung K, Seitz T, Li S, Janosch P, McFerran B, Kaiser C, Fee F, Katsanakis KD, Rose DW, Mischak H, et al. Suppression of Raf-1 kinase activity and MAP kinase signalling by RKIP. Nature. 1999;401(6749):173–7.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Zheng X-M, Wu F-Q, Zhang X, Lin Q-B, Wang J, Guo X-P, Lei C-L, Cheng Z-J, Zou C, Wan J-M. Evolution of the PEBP gene family and selective signature on FT-like clade. J Syst Evol. 2016;54(5):502–10.

    Article  Google Scholar 

  18. 18.

    Serre L, Vallée B, Bureaud N, Schoentgen F, Zelwer C. Crystal structure of the phosphatidylethanolamine-binding protein from bovine brain: a novel structural class of phospholipid-binding proteins. Structure. 1998;6(10):1255–65.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Banfield MJ, Barker JJ, Perry ACF, Brady RL. Function from structure? The crystal structure of human phosphatidylethanolamine-binding protein suggests a role in membrane signal transduction. Structure. 1998;6(10):1245–54.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Tsuji H, Taoka K-I. Chapter five - florigen signaling. In: Machida Y, Lin C, Tamanoi F, editors. The enzymes. London: Academic Press; 2014. p. 113–44.

    Google Scholar 

  21. 21.

    Vallee BS, Coadou G, Labbe H, Sy D, Vovelle F, Schoentgen F. Peptides corresponding to the N- and C-terminal parts of PEBP are well-structured in solution: new insights into their possible interaction with partners in vivo. J Peptide Res. 2008;61(2):47–57.

    Article  Google Scholar 

  22. 22.

    Hengst U, Albrecht H, Hess D, Monard D. The phosphatidylethanolamine-binding protein is the prototype of a novel family of serine protease inhibitors. J Biol Chem. 2001;276(1):535–40.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Keller ET, Fu Z, Brennan M. The role of Raf kinase inhibitor protein (RKIP) in health and disease. Biochem Pharmacol. 2004;68(6):1049–53.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Lorenz K, Lohse MJ, Quitterer U. Protein kinase C switches the Raf kinase inhibitor from Raf-1 to GRK-2. Nature. 2003;426(6966):574–9.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Odabaei G, Chatterjee D, Jazirehi AR, Goodglick L, Yeung K, Bonavida B. Raf-1 kinase inhibitor protein: structure, function, regulation of cell signaling, and pivotal role in apoptosis. Adv Cancer Res. 2004;91:169–200.

    CAS  PubMed  Article  Google Scholar 

  26. 26.

    Bradley D, Carpenter R, Copsey L, Vincent C, Rothstein S, Coen E. Control of inflorescence architecture in Antirrhinum. Nature. 1996;379(6568):791–7.

    CAS  PubMed  Article  Google Scholar 

  27. 27.

    Karlgren A, Gyllenstrand N, Kallman T, Sundstrom JF, Moore D, Lascoux M, Lagercrantz U. Evolution of the PEBP gene family in plants: functional diversification in seed plant evolution. Plant Physiol. 2011;156(4):1967–77.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Vaistij FE, Barros-Galvão T, Cole AF, Gilday AD, He Z, Li Y, Harvey D, Larson TR, Graham IA. MOTHER-OF-FT-AND-TFL1 represses seed germination under far-red light by modulating phytohormone responses in Arabidopsis thaliana. Proc Natl Acad Sci USA. 2018;115(33):8442–7.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Mackenzie KK, Coelho LL, Lütken H, Müller R. Phylogenomic analysis of the PEBP gene family from Kalanchoë. Agronomy. 2019;9(4):171.

    CAS  Article  Google Scholar 

  30. 30.

    Jin S, Nasim Z, Susila H, Ahn JH: Evolution and functional diversification of FLOWERING LOCUS T/TERMINAL FLOWER 1 family genes in plants. Seminars in Cell & Developmental Biology. 2021;109:20–30.

  31. 31.

    Hedman H, Källman T, Lagercrantz U. Early evolution of the MFT-like gene family in plants. Plant Mol Biol. 2009;70(4):359–69.

    CAS  PubMed  Article  Google Scholar 

  32. 32.

    Liu Y-Y, Yang K-Z, Wei X-X, Wang X-Q. Revisiting the phosphatidylethanolamine-binding protein (PEBP) gene family reveals cryptic FLOWERING LOCUS T gene homologs in gymnosperms and sheds new light on functional evolution. New Phytol. 2016;212(3):730–44.

    CAS  PubMed  Article  Google Scholar 

  33. 33.

    Wickland Daniel P, Hanzawa Y. The FLOWERING LOCUS T/TERMINAL FLOWER 1 gene family: functional evolution and molecular mechanisms. Mol Plant. 2015;8(7):983–97.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Wang Z, Zhou Z, Liu Y, Liu T, Li Q, Ji Y, Li C, Fang C, Wang M, Wu M, et al. Functional evolution of phosphatidylethanolamine binding proteins in soybean and Arabidopsis. Plant Cell Online. 2015;27(2):323–36.

    CAS  Article  Google Scholar 

  35. 35.

    Ahn JH, Miller D, Winter VJ, Banfield MJ, Lee JH, Yoo SY, Henz SR, Brady RL, Weigel D. A divergent external loop confers antagonistic activity on floral regulators FT and TFL1. EMBO J. 2006;25(3):605–14.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Bai S, Tuan PA, Saito T, Ito A, Ubi BE, Ban Y, Moriguchi T. Repression of TERMINAL FLOWER1 primarily mediates floral induction in pear (Pyrus pyrifolia Nakai) concomitant with change in gene expression of plant hormone-related genes and transcription factors. J Exp Bot. 2017;68(17):4899–914.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Hanzawa Y, Money T, Bradley D. A single amino acid converts a repressor to an activator of flowering. Proc Natl Acad Sci USA. 2005;102(21):7748–53.

    CAS  PubMed  Article  Google Scholar 

  38. 38.

    Koskela EA, Sønsteby A, Flachowsky H, Heide OM, Hanke M-V, Elomaa P, Hytönen T. TERMINAL FLOWER1 is a breeding target for a novel everbearing trait and tailored flowering responses in cultivated strawberry (Fragaria × ananassa Duch.). Plant Biotechnol J. 2016;14(9):1852–61.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  39. 39.

    Baumann K, Venail J, Berbel A, Domenech MJ, Money T, Conti L, Hanzawa Y, Madueno F, Bradley D. Changing the spatial pattern of TFL1 expression reveals its key role in the shoot meristem in controlling Arabidopsis flowering architecture. J Exp Bot. 2015;66(15):4769–80.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Chen Q, Payyavula RS, Chen L, Zhang J, Zhang C, Turgeon R. FLOWERING LOCUS T mRNA is synthesized in specialized companion cells in Arabidopsis and Maryland Mammoth tobacco leaf veins. Proc Natl Acad Sci USA. 2018;115(11):2830–5.

    CAS  PubMed  Article  Google Scholar 

  41. 41.

    Luccioni L, Krzymuski M, Sánchez-Lamas M, Karayekov E, Cerdán PD, Casal JJ. CONSTANS delays Arabidopsis flowering under short days. Plant J. 2019;97(5):923–32.

    CAS  PubMed  Article  Google Scholar 

  42. 42.

    Lee C, Kim S-J, Jin S, Susila H, Youn G, Nasim Z, Alavilli H, Chung K-S, Yoo SJ, Ahn JH. Genetic interactions reveal the antagonistic roles of FT/TSF and TFL1 in the determination of inflorescence meristem identity in Arabidopsis. Plant J. 2019;99(3):452–64.

    CAS  PubMed  Article  Google Scholar 

  43. 43.

    Ho WWH, Weigel D. Structural features determining flower-promoting activity of Arabidopsis FLOWERING LOCUS T. Plant Cell. 2014;26(2):552–64.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Footitt S, Ölçer-Footitt H, Hambidge AJ, Finch-Savage WE. A laboratory simulation of Arabidopsis seed dormancy cycling provides new insight into its regulation by clock genes and the dormancy-related genes DOG1, MFT, CIPK23 and PHYA. Plant Cell Environ. 2017;40(8):1474–86.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  45. 45.

    Jin S, Jung HS, Chung KS, Lee JH, Ahn JH. FLOWERING LOCUS T has higher protein mobility than TWIN SISTER OF FT. J Exp Bot. 2015;66(20):6109–17.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Huang N-C, Jane W-N, Chen J, Yu T-S. Arabidopsis thaliana CENTRORADIALIS homologue (ATC) acts systemically to inhibit floral initiation in Arabidopsis. Plant J. 2012;72(2):175–84.

    CAS  PubMed  Article  Google Scholar 

  47. 47.

    Liu X, Zhang J, Abuahmad A, Franks RG, Xie D-Y, Xiang Q-Y. Analysis of two TFL1 homologs of dogwood species (Cornus L.) indicates functional conservation in control of transition to flowering. Planta. 2016;243(5):1129–41.

    CAS  PubMed  Article  Google Scholar 

  48. 48.

    Ryu JY, Lee H-J, Seo PJ, Jung J-H, Ahn JH, Park C-M. The Arabidopsis floral repressor BFT delays flowering by competing with FT for FD binding under high salinity. Mol Plant. 2014;7(2):377–87.

    CAS  PubMed  Article  Google Scholar 

  49. 49.

    Shulaev V, Korban SS, Sosinski B, Abbott AG, Aldwinckle HS, Folta KM, Iezzoni A, Main D, Arus P, Dandekar AM, et al. Multiple models for Rosaceae genomics. Plant Physiol. 2008;147(3):985–1003.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  50. 50.

    Soundararajan P, Won SY, Kim JS. Insight on Rosaceae family with genome sequencing and functional genomics perspective. Biomed Res Int. 2019;2019:1–12.

    Article  CAS  Google Scholar 

  51. 51.

    Yamamoto T, Terakami S. Genomics of pear and other Rosaceae fruit trees. Breed Sci. 2016;66(1):148–59.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  52. 52.

    Xiang Y, Huang CH, Hu Y, Wen J, Li S, Yi T, Chen H, Xiang J, Ma H. Evolution of rosaceae fruit types based on nuclear phylogeny in the context of geological times and genome duplication. Mol Biol Evol. 2017;34(2):262–81.

    CAS  PubMed  Google Scholar 

  53. 53.

    Gottschalk C, van Nocker S. Diversity in seasonal bloom time and floral development among apple species and hybrids. J Am Soc Hortic Sci. 2013;138(5):367–74.

    Article  Google Scholar 

  54. 54.

    Pin PA, Benlloch R, Bonnet D, Wremerth-Weich E, Kraft T, Gielen JJL, Nilsson O. An antagonistic pair of FT homologs mediates the control of flowering time in sugar beet. Science. 2010;330(6009):1397–400.

    CAS  PubMed  Article  Google Scholar 

  55. 55.

    Zhang Z, Zhuo X, Zhao K, Zheng T, Han Y, Yuan C, Zhang Q. Transcriptome profiles reveal the crucial roles of hormone and sugar in the bud dormancy of Prunus mume. Sci Rep. 2018;8(1):1–9.

    Article  CAS  Google Scholar 

  56. 56.

    Amasino R. Seasonal and developmental timing of flowering. Plant J. 2010;61(6):1001–13.

    CAS  PubMed  Article  Google Scholar 

  57. 57.

    Song YH, Ito S, Imaizumi T. Flowering time regulation: photoperiod- and temperature-sensing in leaves. Trends Plant Sci. 2013;18(10):575–83.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  58. 58.

    Komeda Y. Genetic regulation of time to Flower in Arabidopsis thaliana. Annu Rev Plant Biol. 2004;55(1):521–35.

    CAS  PubMed  Article  Google Scholar 

  59. 59.

    Teotia S, Tang G. To bloom or not to bloom: role of MicroRNAs in plant flowering. Mol Plant. 2015;8(3):359–77.

    CAS  PubMed  Article  Google Scholar 

  60. 60.

    Ksikiewicz M, Rychel S, Nelson MN, Wyrwa K, Naganowska B, Wolko B. Expansion of the phosphatidylethanolamine binding protein family in legumes: a case study of Lupinus angustifolius L. FLOWERING LOCUS T homologs, LanFTc1 and LanFTc2. BMC Genomics. 2016;17(1):186.

    Article  CAS  Google Scholar 

  61. 61.

    Wu J, Wang Z, Shi Z, Zhang S, Ming R, Zhu S, Khan MA, Tao S, Korban SS, Wang H, et al. The genome of the pear (Pyrus bretschneideri Rehd). Genome Res. 2012;23(2):396–408.

    PubMed  Article  CAS  Google Scholar 

  62. 62.

    Xu G, Ma H, Nei M, Kong H. Evolution of F-box genes in plants: different modes of sequence divergence and their relationships with functional diversification. Proc Natl Acad Sci USA. 2009;106(3):835–40.

    CAS  PubMed  Article  Google Scholar 

  63. 63.

    Gao J, Huang B-H, Wan Y-T, Chang J, Li J-Q, Liao P-C. Functional divergence and intron variability during evolution of angiosperm TERMINAL FLOWER1 (TFL1) genes. Sci Rep. 2017;7(1):1–13.

    Article  CAS  Google Scholar 

  64. 64.

    Yamaguchi A, Kobayashi Y, Goto K, Abe M, Araki T. TWIN SISTER OF FT (TSF) acts as a floral pathway integrator redundantly with FT. Plant Cell Physiol. 2005;46(8):1175–89.

    CAS  PubMed  Article  Google Scholar 

  65. 65.

    Hu Y, Liu B, Sun C, Liu J, Wang X, Cheng F, Liang J, Wang X, Wu J. Evolution of TWIN SISTER of FT (TSF) Genes in Brassicaceae. Hortic Plant J. 2016;2(1):16–25.

    Article  Google Scholar 

  66. 66.

    Komar AA. The Yin and Yang of codon usage. Hum Mol Genet. 2016;25(R2):R77–85.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  67. 67.

    Yu C-H, Dang Y, Zhou Z, Wu C, Zhao F, Sachs Matthew S, Liu Y. Codon usage influences the local rate of translation elongation to regulate co-translational protein folding. Mol Cell. 2015;59(5):744–54.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  68. 68.

    Eyre-Walker AC. An analysis of codon usage in mammals: selection or mutation bias? J Mol Evol. 1991;33(5):442–9.

    CAS  PubMed  Article  Google Scholar 

  69. 69.

    Guan D-L, Ma L-B, Khan MS, Zhang X-X, Xu S-Q, Xie J-Y. Analysis of codon usage patterns in Hirudinaria manillensis reveals a preference for GC-ending codons caused by dominant selection constraints. BMC Genomics. 2018;19(1):1–7.

    Article  CAS  Google Scholar 

  70. 70.

    Sharp PM, Li W-H. The codon adaptation index-a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 1987;15(3):1281–95.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  71. 71.

    Wright F. The ‘effective number of codons’ used in a gene. Gene. 1990;87(1):23–9.

    CAS  PubMed  Article  Google Scholar 

  72. 72.

    Banfield MJ, Brady RL. The structure of Antirrhinum centroradialis protein (CEN) suggests a role as a kinase regulator. J Mol Biol. 2000;297(5):1159–70.

    CAS  PubMed  Article  Google Scholar 

  73. 73.

    Si Z, Liu H, Zhu J, Chen J, Wang Q, Fang L, Gao F, Tian Y, Chen Y, Chang L, et al. Mutation of SELF-PRUNING homologs in cotton promotes short-branching plant architecture. J Exp Bot. 2018;69(10):2543–53.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  74. 74.

    Kotoda N, Hayashi H, Suzuki M, Igarashi M, Hatsuyama Y, Kidou S-I, Igasaki T, Nishiguchi M, Yano K, Shimizu T, et al. Molecular characterization of FLOWERING LOCUS T-Like genes of Apple (Malus × domestica Borkh.). Plant Cell Physiol. 2010;51(4):561–75.

    CAS  PubMed  Article  Google Scholar 

  75. 75.

    Wang Z, Yang R, Devisetty UK, Maloof JN, Zuo Y, Li J, Shen Y, Zhao J, Bao M, Ning G. The divergence of flowering time modulated by FT/TFL1 is independent to their interaction and binding activities. Front Plant Sci. 2017;8:1–9.

    Google Scholar 

  76. 76.

    Zhang X, An L, Nguyen TH, Liang H, Wang R, Liu X, Li T, Qi Y, Yu F. The cloning and functional characterization of peach CONSTANS and FLOWERING LOCUS T homologous genes PpCO and PpFT. PLoS ONE. 2015;10(4):e0124108.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  77. 77.

    Mimida N, Kotoda N, Ueda T, Igarashi M, Hatsuyama Y, Iwanami H, Moriya S, Abe K. Four TFL1/CEN-like genes on distinct linkage groups show different expression patterns to regulate vegetative and reproductive development in apple (Malus×domestica Borkh.). Plant Cell Physiol. 2009;50(2):394–412.

    CAS  PubMed  Article  Google Scholar 

  78. 78.

    Esumi T, Kitamura Y, Hagihara C, Yamane H, Tao R. Identification of a TFL1 ortholog in Japanese apricot (Prunus mume Sieb. et Zucc.). Sci Hortic. 2010;125(4):608–16.

    CAS  Article  Google Scholar 

  79. 79.

    Hanano S, Goto K. Arabidopsis TERMINAL FLOWER1 is involved in the regulation of flowering time and inflorescence development through transcriptional repression. Plant Cell. 2011;23(9):3172–84.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  80. 80.

    Freiman A, Golobovitch S, Yablovitz Z, Belausov E, Dahan Y, Peer R, Avraham L, Freiman Z, Evenor D, Reuveni M, et al. Expression of flowering locus T2 transgene from Pyrus communis L. delays dormancy and leaf senescence in Malus × domestica Borkh, and causes early flowering in tobacco. Plant Sci. 2015;241:164–76.

    CAS  PubMed  Article  Google Scholar 

  81. 81.

    Hsu CY, Adams JP, Kim H, No K, Ma C, Strauss SH, Drnevich J, Vandervelde L, Ellis JD, Rice BM, et al. FLOWERING LOCUS T duplication coordinates reproductive and vegetative growth in perennial poplar. Proc Natl Acad Sci USA. 2011;108(26):10756–61.

    CAS  PubMed  Article  Google Scholar 

  82. 82.

    Park S, Srinivasan C, Dardick C, Callahan A, Scorza R. Plum (Prunus domestica) trees transformed with poplar FT1 result in altered architecture, dormancy requirement, and continuous flowering. PLoS ONE. 2012;7(7):e40715.

    Article  CAS  Google Scholar 

  83. 83.

    Klintenäs M, Pin PA, Benlloch R, Ingvarsson PK, Nilsson O. Analysis of conifer FLOWERING LOCUS T/TERMINAL FLOWER1-like genes provides evidence for dramatic biochemical evolution in the angiosperm FT lineage. New Phytol. 2012;196(4):1260–73.

    PubMed  Article  CAS  Google Scholar 

  84. 84.

    Avia K, Kärkkäinen K, Lagercrantz U, Savolainen O. Association of FLOWERING LOCUS T/TERMINAL FLOWER 1-like gene FTL2 expression with growth rhythm in Scots pine (Pinus sylvestris). New Phytol. 2014;204(1):159–70.

    CAS  PubMed  Article  Google Scholar 

  85. 85.

    Gyllenstrand N, Clapham D, Kallman T, Lagercrantz U. A Norway spruce FLOWERING LOCUS T homolog is implicated in control of growth rhythm in conifers. Plant Physiol. 2007;144(1):248–57.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  86. 86.

    Karlgren A, Gyllenstrand N, Clapham D, Lagercrantz U. FLOWERING LOCUS T/TERMINAL FLOWER1-like genes affect growth rhythm and bud set in Norway spruce. Plant Physiol. 2013;163(2):792–803.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  87. 87.

    Mouradov A, Cremer F, Coupland G. Control of flowering time. Plant Cell. 2002;14(suppl 1):S111–30.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  88. 88.

    Putterill J, Laurie R, Macknight R. It’s time to flower: the genetic control of flowering time. BioEssays. 2004;26(4):363–73.

    CAS  PubMed  Article  Google Scholar 

  89. 89.

    Song YH, Shim JS, Kinmonth-Schultz HA, Imaizumi T. Photoperiodic flowering: time measurement mechanisms in leaves. Annu Rev Plant Biol. 2015;66(1):441–64.

    CAS  PubMed  Article  Google Scholar 

  90. 90.

    Samach A, Gover A. Photoperiodism: the consistent use of CONSTANS. Curr Biol. 2001;11(16):R651–4.

    CAS  PubMed  Article  Google Scholar 

  91. 91.

    Yang S, Weers BD, Morishige DT, Mullet JE. CONSTANS is a photoperiod regulated activator of flowering in sorghum. BMC Plant Biol. 2014;14(1):148.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  92. 92.

    Hayama R, Coupland G. The molecular basis of diversity in the photoperiodic flowering responses of Arabidopsis and rice. Plant Physiol. 2004;135(2):677–84.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  93. 93.

    Distelfeld A, Li C, Dubcovsky J. Regulation of flowering in temperate cereals. Curr Opin Plant Biol. 2009;12(2):178–84.

    CAS  PubMed  Article  Google Scholar 

  94. 94.

    Böhlenius H, Huang T, Charbonnel-Campaa L, Brunner AM, Jansson S, Strauss SH, Nilsson O. CO/FT regulatory module controls timing of flowering and seasonal growth cessation in trees. Science. 2006;312(5776):1040–3.

    PubMed  Article  CAS  Google Scholar 

  95. 95.

    Deng W, Ying H, Helliwell CA, Taylor JM, Peacock WJ, Dennis ES. FLOWERING LOCUS C (FLC) regulates development pathways throughout the life cycle of Arabidopsis. Proc Natl Acad Sci USA. 2011;108(16):6680–5.

    CAS  PubMed  Article  Google Scholar 

  96. 96.

    Hartmann U, Hohmann S, Nettesheim K, Wisman E, Saedler H, Huijser P. Molecular cloning of SVP: a negative regulator of the floral transition in Arabidopsis. Plant J. 2000;21(4):351–60.

    CAS  PubMed  Article  Google Scholar 

  97. 97.

    Kumar SV, Lucyshyn D, Jaeger KE, Alós E, Alvey E, Harberd NP, Wigge PA. Transcription factor PIF4 controls the thermosensory activation of flowering. Nature. 2012;484(7393):242–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  98. 98.

    Lee J, Lee I. Regulation and function of SOC1, a flowering pathway integrator. J Exp Bot. 2010;61(9):2247–54.

    CAS  PubMed  Article  Google Scholar 

  99. 99.

    Teper-Bamnolker P. The flowering integrator FT regulates SEPALLATA3 and FRUITFULL accumulation in Arabidopsis leaves. Plant Cell Online. 2005;17(10):2661–75.

    CAS  Article  Google Scholar 

  100. 100.

    Pin PA, Nilsson O. The multifaceted roles of FLOWERING LOCUS T in plant development. Plant Cell Environ. 2012;35(10):1742–55.

    CAS  PubMed  Article  Google Scholar 

  101. 101.

    Rinne PLH, Welling A, Vahala J, Ripel L, Ruonala R, Kangasjärvi J, van der Schoot C. Chilling of dormant buds hyperinduces FLOWERING LOCUS T and recruits GA-inducible 1,3-β-Glucanases to reopen signal conduits and release dormancy in Populus. Plant Cell. 2011;23(1):130–46.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  102. 102.

    Niu Q, Li J, Cai D, Qian M, Jia H, Bai S, Hussain S, Liu G, Teng Y, Zheng X. Dormancy-associated MADS-box genes and microRNAs jointly control dormancy transition in pear (Pyrus pyrifolia white pear group) flower bud. J Exp Bot. 2016;67(1):239–57.

    CAS  PubMed  Article  Google Scholar 

  103. 103.

    Ito A, Saito T, Sakamoto D, Sugiura T, Bai S, Moriguchi T, Näsholm T. Physiological differences between bud breaking and flowering after dormancy completion revealed by DAM and FT/TFL1 expression in Japanese pear (Pyrus pyrifolia). Tree Physiol. 2016;36(1):109–20.

    PubMed  Article  Google Scholar 

  104. 104.

    Rongmei W, Sumathi T, Sakuntala K, Tustin SD, Hellens RP, Allan AC, Macknight RC, Erika VG. SVP-like MADS box genes control dormancy and budbreak in apple. Front Plant Sci. 2017;08(55):1–9.

    Google Scholar 

  105. 105.

    Mateos JL, Madrigal P, Tsuda K, Rawat V, Richter R, Romera-Branchat M, Fornara F, Schneeberger K, Krajewski P, Coupland G. Combinatorial activities of SHORT VEGETATIVE PHASE and FLOWERING LOCUS C define distinct modes of flowering regulation in Arabidopsis. Genome Biol. 2015;16(1):31.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  106. 106.

    Kurokura T, Mimida N, Battey NH, Hytönen T. The regulation of seasonal flowering in the Rosaceae. J Exp Bot. 2013;64(14):4131–41.

    CAS  PubMed  Article  Google Scholar 

  107. 107.

    Verde I, Jenkins J, Dondini L, Micali S, Pagliarani G, Vendramin E, Paris R, Aramini V, Gazza L, Rossini L, et al. The Peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity. BMC Genomics. 2017;18(1):87.

    Article  CAS  Google Scholar 

  108. 108.

    Zhang Q, Chen W, Sun L, Zhao F, Huang B, Yang W, Tao Y, Wang J, Yuan Z, Fan G, et al. The genome of Prunus mume. Nat Commun. 2012;3(1):1–13.

    Google Scholar 

  109. 109.

    Baek S, Choi K, Kim G-B, Yu H-J, Cho A, Jang H, Kim C, Kim H-J, Chang KS, Kim J-H, et al. Draft genome sequence of wild Prunus yedoensis reveals massive inter-specific hybridization between sympatric flowering cherries. Genome Biol. 2018;19(1):1–9.

    Article  CAS  Google Scholar 

  110. 110.

    Shirasawa K, Isuzugawa K, Ikenaga M, Saito Y, Yamamoto T, Hirakawa H, Isobe S. The genome sequence of sweet cherry (Prunus avium) for use in genomics-assisted breeding. DNA Res. 2017;24(5):499–508.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  111. 111.

    Jiang F, Zhang J, Wang S, Yang L, Luo Y, Gao S, Zhang M, Wu S, Hu S, Sun H, et al. The apricot (Prunus armeniaca L.) genome elucidates Rosaceae evolution and beta-carotenoid synthesis. Hortic Res. 2019;6(1):18.

    Article  CAS  Google Scholar 

  112. 112.

    Daccord N, Celton J-M, Linsmith G, Becker C, Choisne N, Schijlen E, van de Geest H, Bianco L, Micheletti D, Velasco R, et al. High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development. Nat Genet. 2017;49(7):1099–106.

    CAS  PubMed  Article  Google Scholar 

  113. 113.

    Tinker NA, Chagné D, Crowhurst RN, Pindo M, Thrimawithana A, Deng C, Ireland H, Fiers M, Dzierzon H, Cestaro A, et al. The draft genome sequence of European Pear (Pyrus communis L. ‘Bartlett’). PLoS ONE. 2014;9(4):92644.

    Article  CAS  Google Scholar 

  114. 114.

    VanBuren R, Wai CM, Colle M, Wang J, Sullivan S, Bushakra JM, Liachko I, Vining KJ, Dossett M, Finn CE, et al. A near complete, chromosome-scale assembly of the black raspberry (Rubus occidentalis) genome. GigaScience. 2018;7(8):1–9.

    PubMed  Article  CAS  Google Scholar 

  115. 115.

    Jung S, Jesudurai C, Staton M, Du Z, Ficklin S, Cho I, Abbott A, Tomkins J, Main D. GDR (Genome Database for Rosaceae): integrated web resources for Rosaceae genomics and genetics research. BMC Bioinf. 2004;5(1):130.

    Article  CAS  Google Scholar 

  116. 116.

    Potter SC, Luciani A, Eddy SR, Park Y, Lopez R, Finn RD. HMMER web server: 2018 update. Nucleic Acids Res. 2018;46(W1):W200–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  117. 117.

    Letunic I, Doerks T, Bork P. SMART: recent updates, new developments and status in 2015. Nucleic Acids Res. 2015;43(D1):D257–60.

    CAS  PubMed  Article  Google Scholar 

  118. 118.

    Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42(D1):D222–30.

    CAS  PubMed  Article  Google Scholar 

  119. 119.

    Marchler-Bauer A, Bo Y, Han L, He J, Lanczycki CJ, Lu S, Chitsaz F, Derbyshire MK, Geer RC, Gonzales NR, et al. CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2017;45(D1):D200–3.

    CAS  PubMed  Article  Google Scholar 

  120. 120.

    Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  121. 121.

    Nicholas KB: GeneDoc: Analysis and visualization of genetic variation, EMBNEW. Embnew News 1997, 4.

  122. 122.

    Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.

    CAS  Article  Google Scholar 

  123. 123.

    Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  124. 124.

    Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19(12):1572–4.

    CAS  Article  Google Scholar 

  125. 125.

    Bailey TL, Johnson J, Grant CE, Noble WS. The MEME suite. Nucleic Acids Res. 2015;43(W1):W39–49.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  126. 126.

    Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202.

    CAS  PubMed  Article  Google Scholar 

  127. 127.

    Wang Y, Tang H, DeBarry JD, Tan X, Li J, Wang X, Lee Th, 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(7):e49–e49.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  128. 128.

    García-Gil MR, Lecoy J: Evolutionary analysis of LP3 gene family in conifers: an ASR homolog. bioRxiv 2020.

  129. 129.

    Puigbò P, Bravo IG, Garcia-Vallve S. CAIcal: a combined set of tools to assess codon usage adaptation. Biol Direct. 2008;3(1):38.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  130. 130.

    Fuglsang A. The ‘effective number of codons’ revisited. Biochem Biophys Res Commun. 2004;317(3):957–64.

    CAS  PubMed  Article  Google Scholar 

  131. 131.

    Xu X-z, Liu Q-p, Fan L-j, Cui X-f, Zhou X-p. Analysis of synonymous codon usage and evolution of begomoviruses. J Zhejiang Univ Sci B. 2008;9(9):667–74.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  132. 132.

    Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17(4):540–52.

    CAS  PubMed  Article  Google Scholar 

  133. 133.

    Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24(8):1586–91.

    CAS  PubMed  Article  Google Scholar 

  134. 134.

    Stern A, Doron-Faigenboim A, Erez E, Martz E, Bacharach E, Pupko T. Selecton 2007: advanced models for detecting positive and purifying selection using a Bayesian inference approach. Nucleic Acids Res. 2007;35:W506–11.

    PubMed  PubMed Central  Article  Google Scholar 

  135. 135.

    Lescot M. 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(1):325–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  136. 136.

    Chow C-N, Zheng H-Q, Wu N-Y, Chien C-H, Huang H-D, Lee T-Y, Chiang-Hsieh Y-F, Hou P-F, Yang T-Y, Chang W-C. PlantPAN 20: an update of plant promoter analysis navigator for reconstructing transcriptional regulatory networks in plants. Nucleic Acids Res. 2016;44(D1):D1154–60.

    CAS  PubMed  Article  Google Scholar 

  137. 137.

    Lazzari B, Caprera A, Vecchietti A, Stella A, Milanesi L, Pozzi C. ESTree db: a tool for peach functional genomics. BMC Bioinformatics. 2005;6(Suppl 4):S16.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  138. 138.

    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  139. 139.

    Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  140. 140.

    Anders S, Pyl PT, Huber W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics. 2014;31(2):166–9.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  141. 141.

    Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009;26(1):139–40.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  142. 142.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25(4):402–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  143. 143.

    Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9(1):1–9.

    Article  CAS  Google Scholar 

  144. 144.

    Saito R, Smoot ME, Ono K, Ruscheinski J, Wang P-L, Lotia S, Pico AR, Bader GD, Ideker T. A travel guide to Cytoscape plugins. Nat Methods. 2012;9(11):1069–76.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgements

We thank Dr. Jason Holliday (Virginia Tech) for his help in revising the manuscript.

Funding

This work was financially supported by the Fundamental Research Funds for the Central Universities, China (Grant No. BLX201828) and Special Fund for Beijing Common Construction Project.

Author information

Affiliations

Authors

Contributions

MZ designed and conducted the study; PL performed the blast analysis and sequence curation; XY assisted in sample collection; JW and TC provided help with transcriptome analysis; QZ supervised the project and revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Qixiang Zhang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that no competing interests exist.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file1: Fig. S1.

Phylogenetic trees were constructed using three different algorithms.

Additional file2: Fig. S2.

Exon-intron distributions of PEBP family genes from nine Rosaceae species and A. thaliana.

Additional file 3: Fig. S3.

Protein sequence alignment of FT and TFL1-like proteins in Arabidopsis and nine Rosaceae species.

Additional file4: Fig. S4.

Inter-genomic synteny blocks between species.

Additional file 5: Fig. S5.

Distribution and collinearity of PEBP family genes within genomes of nine Rosaceae species.

Additional file 6: Fig. S6.

Heatmap of average RSCU scores estimated for different codons in five Rosaceae PEBP gene lineages.

Additional file 7: Fig. S7.

Comparison of gene parameters including the (a) CAI, (b) total GC%, (c) ENC, (d) GC1%, (e) GC2%, (f) GC3% for all PEBP genes across nine Rosaceae species.

Additional file 8: Fig. S8.

Selective pressure analysis of the BFT, CEN, and MFT lineages and PEBP genes as a whole identified strong purifying selection across protein sites.

Additional file 9: Fig. S9.

Relative expression of PEBP family genes in floral bud, leaf, stem, and root tissues of P. mume.

Additional file 10: Fig. S10.

Weighted gene co-expression network analysis (WGCNA) during four stages of floral bud blooming in P. mume.

Additional file 11: TableS1.

Numbers of genes originated from different types of duplication events in the genomes of Arabidopsis and seven Rosaceae species.

Additional file 12: TableS2.

The duplication modes of PEBP genes in the genomes of Arabidopsis and seven Rosaceae species.

Additional file 13: TableS3.

Pearson correlation between gene parameters estimated for 56 Rosaceae PEBP genes.

Additional file 14: TableS4.

Parameter estimates and likelihood values for branch models among PEBP lineages.

Additional file  15: TableS5.

Parameter estimates and likelihood values for site models across PEBP protein sites.

Additional file 16: TableS6.

Candidate genes co-expressed with PmFTduring dormancy release in P. mume.

Additional file 17: TableS7.

Primers used in qRT-PCR analysis for PEBPfamily genes in P. mume.

Additional file 18: TableS8.

Primers used in qRT-PCR analysis for FTco-expressed genes in P. mume.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, M., Li, P., Yan, X. et al. Genome-wide characterization of PEBP family genes in nine Rosaceae tree species and their expression analysis in P. mume. BMC Ecol Evo 21, 32 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12862-021-01762-4

Download citation

Keywords

  • PEBP gene family
  • Rosaceae species
  • Flowering time regulation in perennials
  • Gene family evolution
  • Floral bud break