Identification and phylogenetic analysis of PgCBF family members
Based on Pennisetum multi-omics database, a total of 23 CBF family members were identified and named PgCBF1 to PgCBF23. Using ExPASy, the physicochemical properties of these 23 PgCBF proteins were predicted (Table S2). The amino acid numbers ranged from 187 to 324, and the molecular weights ranged from 19,688.36 to 34,271.70 Da. The theoretical isoelectric points (pI) of PgCBF proteins ranged from 4.57 to 9.82, with only PgCBF13 and PgCBF23 having pI values greater than 7, indicating alkaline proteins, while the others were acidic. All 23 PgCBF proteins were predicted as unstable (instability index > 40) and showed a range of lipid solubility coefficients from 50.90 to 77.30. The average hydrophilicity of these proteins was negative, suggesting they are hydrophilic. Additionally, subcellular localization predictions indicated that most PgCBF genes are localized in the nucleus. For phylogenetic analysis of PgCBF family members, amino acid sequences of CBF family members from Arabidopsis were downloaded from the NCBI database. Using MEGA7, a phylogenetic tree was constructed to analyze the evolutionary relationship between PgCBF and AtCBF family members (Fig. 1). The results revealed that the CBF family members could be classified into 4 groups: Group A included 8 PgCBFs, Group B included 10 PgCBFs, Group C included 3 PgCBFs, and Group D included 2 PgCBFs, respectively.
Conserved motif of PgCBF family members
To further investigate the amino acid structural features of PgCBF family members, MEME online was used to analyze the conserved motifs (Fig. 2). A total of 10 conserved motifs were identified among the 23 PgCBF family members. The analysis revealed significant variation in motif numbers among different family members, ranging from 6 to 9 motifs. All PgCBF members contained motifs 1, 2, 3, 4, and 6 in the same order, indicating high conservation of these 5 motifs. Some motifs were specific to certain PgCBF members; for instance, motif 9 was present only in PgCBF14, PgCBF19, PgCBF20, and PgCBF21, while motif 10 was found in PgCBF2, PgCBF3, and PgCBF13. Further amino acid conserved domain analysis revealed that all PgCBF family members contained the AP2 domain (Fig. 2). Specifically, 18 PgCBF members had the AP2 domain, and 5 PgCBF members (PgCBF4, PgCBF14, PgCBF16, PgCBF22, and PgCBF23) had the AP2 superfamily domain.
Chromosomal localization analysis
Visualization analysis of the chromosomal positions of 23 PgCBF genes was conducted using TBtools software, revealing uneven distribution across 7 chromosomes (Fig. S1). Specifically, PgCBF1 is located on chromosome 1; PgCBF2, PgCBF3, and PgCBF4 are on chromosome 2. Chromosome 3 harbors 7 PgCBF genes, the highest among all chromosomes containing PgCBF genes. Chromosome 4 contains 3 PgCBF genes, including PgCBF13, PgCBF14, and PgCBF15. Chromosome 5 hosts 2 PgCBF genes, PgCBF16 and PgCBF17. Chromosome 6 only includes PgCBF18 gene. Chromosome 7 accommodates 5 PgCBF genes: PgCBF19, PgCBF20, PgCBF21, PgCBF22, and PgCBF23, forming a gene cluster at the upper region of the chromosome.
Cis -acting element analysis
Cis-acting elements play a crucial role in determining gene function and expression patterns. Using the Plant CARE website, we analyzed all cis-acting elements within the 2000 bp upstream promoter regions of PgCBF genes. The promoter regions of the PgCBFs contain multiple cis-acting elements primarily associated with light response, plant response to stress conditions, plant hormones, and plant growth and development. As shown in Fig. 3, all PgCBF genes possess light-responsive elements, stress-responsive elements (including anaerobic induction, drought-responsive elements, and low-temperature responsive elements), plant hormone-related elements (abscisic acid, auxin, gibberellin, and salicylic acid-responsive elements), and elements related to plant growth and development (including seed-specific regulation and endosperm expression elements). These results indicate that PgCBF genes play roles throughout the growth and development processes. Additionally, the abundance of stress-related cis-elements in the promoter regions of the PgCBF genes suggests their involvement in complex mechanisms during plant responses to stress conditions.
Tissue-specific expression patterns of PgCBF family genes
To explore the tissue-specific expression pattern of PgCBF genes, Pennisetum hybrid 'Bond 1' seedlings at different stages were used as materials, and RT-qPCR was employed to detect the expression levels of 6 PgCBF genes in leaves, stems, and roots at four-week-old seedlings (Fig. 4). The results indicate significant differences in the expression patterns of PgCBF genes across different tissues. PgCBF10, PgCBF11, PgCBF19, and PgCBF20 show markedly higher expression level in leaves compared to stems and roots. PgCBF12 and PgCBF21 exhibit highest expression in roots. Interestingly, PgCBF21 shows the lowest expression in leaves, other genes generally exhibit the lowest expression in stems. Analysis of the tissue-specific expression patterns of PgCBF genes reveals differential expression patterns during different stages of plant growth, suggesting their involvement across various developmental stages of plant growth and development.
Expression patterns of PgCBF genes under cold stress
To investigate whether PgCBF genes are induced under cold stress, Pennisetum hybrid 'Bond 1' seedlings were subjected to 4°C for 8 h, 16 h, 24 h, and 32 h to simulate low temperature stress condition. We measured the expression levels of 6 PgCBF genes in roots, stems, and leaves of seedlings under cold stress, respectively (Fig. 5–7). Specifically, after 32 h of cold treatment, the expression level of PgCBF10 showed significant upregulation in stems, with expression increasing more than 7-fold, while its expression in leaves showed a downward trend. The expression level of PgCBF11 exhibited significant upregulation in stems and roots after cold stress treatment, with expression in stems increasing by over 40-fold after 32 h of cold stress treatment and in roots increasing by around 20-fold after 8 h of cold stress treatment. PgCBF12 showed highest expression in leaves and stems after 8 h of cold stress treatment, with expression increasing by more than 3-fold and 5-fold, respectively, compared to untreated condition. After 16 h of cold stress treatment, the expression level of PgCBF12 in roots reached its peak, increasing by more than 15-fold. PgCBF19 exhibited highest expression in leaves and stems after 8 h of cold stress treatment, with expression increasing by approximately 3-fold and 2.2-fold, respectively, compared to untreated conditions. PgCBF20 showed highest expression in leaves and stems after 32 h of cold stress treatment, with expression increasing by approximately 3-fold and 2.8-fold, respectively, compared to untreated condition. Interestingly, the expression level of PgCBF21 in leaf exhibited significant changes after 8 h of cold stress treatment, showing upregulation by more than 37-fold, while its expression in roots showed a downward trend. Under cold stress condition, significant expression differences of the 6 PgCBF genes were observed in roots, stems, and leaves, suggesting their potential roles in various aspects of cold stress response.
Function validation of PgCBF21 in yeast
Based on the expression patterns of six PgCBF genes under cold Stress, PgCBF21 showed higher expression level, so PgCBF21 was used as a candidate gene to analyze its biological function. To further scrutinize the evolutionary relationships, amino acid sequences of CBF21 from fifteen different plant species were utilized to construct a phylogenetic tree. The resulting cladogram delineated two major groups, with the amino acid sequence of PgCBF21 displaying a high degree of identity to SiCBF21 from Setaria italica (Fig. S2A). Next, we analyze the conserved domains of PgCBF21 protein and found that PgCBF21 protein harbors AP2 conserved domain (Fig. S2B).
In addition, we transferred the pYES2-NTB-PgCBF21 vector into yeast and investigated the function of PgCBF21 under cold stress. We observed that the bacterial solution of the negative control (pYES2-NTB) was able to grow under conditions of SG-Ura treated with at 4℃, 0℃, -10℃, -20℃ for 24 h, 48 h, respectively. However, it failed to thrive on the SG-Ura treated with at -20℃ for 72 h and 96 h, respectively. In contrast, the growth pattern of the bacterial solution from the experimental group (pYES2-NTB-PgCBF21) differed from that of the negative control group, as it exhibited growth on the SG-Ura treated with at -20℃ for 72 h and 96 h, respectively (Fig. 8).